Patient-derived models link re-entrant driver localization in atrial fibrillation to fibrosis spatial pattern

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Atrial fibrillation in cancer, anticancer therapies, and underlying mechanisms

Atrial fibrillation (AF) is a common arrhythmic complication in cancer patients and can be exacerbated by traditional cytotoxic and targeted anticancer therapies. Increased incidence of AF in cancer patients is independent of confounding factors, including preexisting myocardial arrhythmogenic substrates, type of cancer, or cancer stage. Mechanistically, AF is characterized by fast unsynchronized atrial contractions with rapid ventricular response, which impairs ventricular filling and results in various symptoms such as fatigue, chest pain, and shortness of breath. Due to increased blood stasis, a consequence of both cancer and AF, concern for stroke increases in this patient population. To compound matters, cardiotoxic anticancer therapies themselves promote AF; thereby exacerbating AF morbidity and mortality in cancer patients. In this review, we examine the relationship between AF, cancer, and anticancer therapies with a focus on the shared molecular and electrophysiological mechanisms linking these disease processes. We also explore the potential role of sodium-glucose co-transporter 2 inhibitors (SGLT2i) in the management of anticancer-therapy induced AF.

Patient-specific multi-physics simulations of fibrotic changes in left atrial tissue mechanics impact on hemodynamics

Stroke is a leading cause of death and disability worldwide. Atrial myopathy, including fibrosis, is associated with an increased risk of ischemic stroke, but the mechanisms underlying this association are poorly understood. Fibrosis modifies myocardial structure, impairing electrical propagation and tissue biomechanics, and creating stagnant flow regions where clots could form. Fibrosis can be mapped non-invasively using late gadolinium enhancement magnetic resonance imaging (LGE-MRI). However, fibrosis maps are not currently incorporated into stroke risk calculations or computational electro-mechano-fluidic models. We present multi-physics simulations of left atrial (LA) myocardial motion and hemodynamics using patient-specific anatomies and fibrotic maps from LGE-MRI. We modify tissue stiffness and active tension generation in fibrotic regions and investigate how these changes affect LA flow for different fibrotic burdens. We find that fibrotic regions and, to a lesser extent, non-fibrotic regions experience reduced myocardial strain, resulting in decreased LA emptying fraction consistent with clinical observations. Both fibrotic tissue stiffening and hypocontractility independently reduce LA function, but together, these two alterations cause more pronounced effects than either one alone. Fibrosis significantly alters flow patterns throughout the atrial chamber, and particularly, the filling and emptying jets of the left atrial appendage (LAA). The effects of fibrosis in LA flow are largely captured by the concomitant changes in LA emptying fraction except inside the LAA, where a multi-factorial behavior is observed. This work illustrates how high-fidelity, multi-physics models can be used to study thrombogenesis mechanisms in a patient-specific manner, shedding light onto the link between atrial fibrosis and ischemic stroke.

Key points

Left atrial (LA) fibrosis is associated with arrhythmogenesis and increased risk of ischemic stroke; its extent and pattern can be quantified on a patient-specific basis using late gadolinium enhancement magnetic resonance imaging.Current stroke risk prediction tools have limited personalization, and their accuracy could be improvedfib by incorporating patient-specific information like fibrotic maps and hemodynamic patterns.We present the first electro-mechano-fluidic multi-physics computational simulations of LA flow, including fibrosis and anatomies from medical imaging.Mechanical changes in fibrotic tissue impair global LA motion, decreasing LA and left atrial appendage (LAA) emptying fractions, especially in subjects with higher fibrosis burdens.Fibrotic-mediated LA motion impairment alters LA and LAA flow near the endocardium and the whole cavity, ultimately leading to more stagnant blood regions in the LAA.

Mechanisms of stretch-induced electro-anatomical remodeling and atrial arrhythmogenesis

Atrial fibrillation (AF) is the most common cardiac rhythm disorder, often occurring in the setting of atrial distension and elevated myocardialstretch. While various mechano-electrochemical signal transduction pathways have been linked to AF development and progression, the underlying molecular mechanisms remain poorly understood, hampering AF therapies. In this review, we describe different aspects of stretch-induced electro-anatomical remodeling as seen in animal models and in patients with AF. Specifically, we focus on cellular and molecular mechanisms that are responsible for mechano-electrochemical signal transduction and the development of ectopic beats triggering AF from pulmonary veins, the most common source of paroxysmal AF. Furthermore, we describe structural changes caused by stretch occurring before and shortly after the onset of AF as well as during AF progression, contributing to longstanding forms of AF. We also propose mechanical stretch as a new dimension to the concept "AF begets AF", in addition to underlying diseases. Finally, we discuss the mechanisms of these electro-anatomical alterations in a search for potential therapeutic strategies and the development of novel antiarrhythmic drugs targeted at the components of mechano-electrochemical signal transduction not only in cardiac myocytes, but also in cardiac non-myocyte cells.

Role of Machine Learning and Artificial Intelligence in Arrhythmias and Electrophysiology

Machine learning (ML), a subset of artificial intelligence (AI) centered on machines learning from extensive datasets, stands at the forefront of a technological revolution shaping various facets of society. Cardiovascular medicine has emerged as a key domain for ML applications, with considerable efforts to integrate these innovations into routine clinical practice. Within cardiac electrophysiology, ML applications, especially in the automated interpretation of electrocardiograms, have garnered substantial attention in existing literature. However, less recognized are the diverse applications of ML in cardiac electrophysiology and arrhythmias, spanning basic science research on arrhythmia mechanisms, both experimental and computational, as well as contributions to enhanced techniques for mapping cardiac electrical function and translational research related to arrhythmia management. This comprehensive review delves into various ML applications within the scope of this journal, organized into 3 parts. The first section provides a fundamental understanding of general ML principles and methodologies, serving as a foundational resource for readers interested in exploring ML applications in arrhythmia research. The second part offers an in-depth review of studies in arrhythmia and electrophysiology that leverage ML methodologies, showcasing the broad potential of ML approaches. Each subject is thoroughly outlined, accompanied by a review of notable ML research advancements. Finally, the review delves into the primary challenges and future perspectives surrounding ML-driven cardiac electrophysiology and arrhythmias research.

Computer based method for identification of fibrotic scars from electrograms and local activation times on the epi- and endocardial surfaces of the ventricles

Cardiac fibrosis stands as one of the most critical conditions leading to lethal cardiac arrhythmias. Identifying the precise location of cardiac fibrosis is crucial for planning clinical interventions in patients with various forms of ventricular and atrial arrhythmias. As fibrosis impedes and alters the path of electrical waves, detecting fibrosis in the heart can be achieved through analyzing electrical signals recorded from its surface. In current clinical practices, it has become feasible to record electrical activity from both the endocardial and epicardial surfaces of the heart. This paper presents a computational method for reconstructing 3D fibrosis using unipolar electrograms obtained from both surfaces of the ventricles. The proposed method calculates the percentage of fibrosis in various ventricular segments by analyzing the local activation times and peak-to-peak amplitudes of the electrograms. Initially, the method was tested using simulated data representing idealized fibrosis in a heart segment; subsequently, it was validated in the left ventricle with fibrosis obtained from a patient with nonischemic cardiomyopathy. The method successfully determined the location and extent of fibrosis in 204 segments of the left ventricle model with an average error of 0.0±4.3% (N = 204). Moreover, the method effectively detected fibrotic scars in the mid-myocardial region, a region known to present challenges in accurate detection using electrogram amplitude as the primary criterion.

Anti- and pro-fibrillatory effects of pulmonary vein isolation gaps in human atrial fibrillation digital twins

Although pulmonary vein isolation (PVI) gaps and extrapulmonary vein triggers contribute to recurrence after atrial fibrillation (AF) ablation, their precise mechanisms remain unproven. Our study assessed the impact of PVI gaps on rhythm outcomes using a human AF digital twin. We included 50 patients (76.0% with persistent AF) who underwent catheter ablation with a realistic AF digital twin by integrating computed tomography and electroanatomical mapping. We evaluated the final rhythm status, including AF and atrial tachycardia (AT), across 600 AF episodes, considering factors including PVI level, PVI gap number, and pacing locations. Our findings revealed that antral PVI had a significantly lower ratio of AF at the final rhythm (28% vs. 56%, p = 0.002) than ostial PVI. Increasing PVI gap numbers correlated with an increased ratio of AF at the final rhythm (p < 0.001). Extra-PV induction yielded a higher ratio of AF at the final rhythm than internal PV induction (77.5% vs. 59.0%, p < 0.001). In conclusion, our human AF digital twin model helped assess AF maintenance mechanisms. Clinical trial registration: https://www.clinicaltrials.gov ; Unique identifier: NCT02138695.

Directed Graph Mapping exceeds Phase Mapping for the detection of simulated 2D meandering rotors in fibrotic tissue with added noise

Cardiac arrhythmias such as atrial fibrillation (AF) are recognised to be associated with re-entry or rotors. A rotor is a wave of excitation in the cardiac tissue that wraps around its refractory tail, causing faster-than-normal periodic excitation. The detection of rotor centres is of crucial importance in guiding ablation strategies for the treatment of arrhythmia. The most popular technique for detecting rotor centres is Phase Mapping (PM), which detects phase singularities derived from the phase of a signal. This method has been proven to be prone to errors, especially in regimes of fibrotic tissue and temporal noise. Recently, a novel technique called Directed Graph Mapping (DGM) was developed to detect rotational activity such as rotors by creating a network of excitation. This research aims to compare the performance of advanced PM techniques versus DGM for the detection of rotors using 64 simulated 2D meandering rotors in the presence of various levels of fibrotic tissue and temporal noise. Four strategies were employed to compare the performances of PM and DGM. These included a visual analysis, a comparison of F2-scores and distance distributions, and calculating p-values using the mid-p McNemar test. Results indicate that in the case of low meandering, fibrosis and noise, PM and DGM yield excellent results and are comparable. However, in the case of high meandering, fibrosis and noise, PM is undeniably prone to errors, mainly in the form of an excess of false positives, resulting in low precision. In contrast, DGM is more robust against these factors as F2-scores remain high, yielding F2≥0.931 as opposed to the best PM F2≥0.635 across all 64 simulations.

Digital twins in medicine

Medical digital twins, which are potentially vital for personalized medicine, have become a recent focus in medical research. Here we present an overview of the state of the art in medical digital twin development, especially in oncology and cardiology, where it is most advanced. We discuss major challenges, such as data integration and privacy, and provide an outlook on future advancements. Emphasizing the importance of this technology in healthcare, we highlight the potential for substantial improvements in patient-specific treatments and diagnostics.

Distributions and number of drivers on real-time phase mapping associated with successful atrial fibrillation termination during catheter ablation for non-paroxysmal atrial fibrillation

BACKGROUND

Real-time phase mapping (ExTRa™) is useful in determining the strategy of catheter ablation for non-paroxysmal atrial fibrillation (AF). This study aimed to investigate the features of drivers of AF associated with its termination during ablation.

METHODS

Thirty-six patients who underwent catheter ablation for non-paroxysmal AF using online real-time phase mapping (ExTRa™) were enrolled. A significant AF driver was defined as an area with a non-passively activated ratio of ≥ 50% on mapping analysis in the left atrium (LA). All drivers were simultaneously evaluated using a low-voltage area, complex fractionated atrial electrogram (CFAE), and rotational activity by unipolar electrogram analysis. The electrical characteristics of drivers were compared between patients with and without AF termination during the procedure.

RESULTS

Twelve patients achieved AF termination during the procedure. The total number of drivers detected on the mapping was significantly lower (4.4 ± 1.6 vs. 7.4 ± 3.8, p = 0.007), and the drivers were more concentrated in limited LA regions (2.8 ± 0.9 vs. 3.9 ± 1.4, p = 0.009) in the termination group than in the non-termination group. The presence of drivers 2-6 with limited (≤ 3) LA regions showed a tenfold increase in the likelihood of AF termination, with 83% specificity and 67% sensitivity. Among 231 AF drivers, the drivers related to termination exhibited a greater overlap of CFAE (56.8 ± 34.1% vs. 39.5 ± 30.4%, p = 0.004) than the non-related drivers. The termination group showed a trend toward a lower recurrence rate after ablation (p = 0.163).

CONCLUSIONS

Rotors responsible for AF maintenance may be characterized in cases with concentrated regions and fewer drivers on mapping.

Fibroblast mediated dynamics in diffusively uncoupled myocytes: a simulation study using 2-cell motifs

In healthy hearts myocytes are typically coupled to nearest neighbours through gap junctions. Under pathological conditions such as fibrosis, or in scar tissue, or across ablation lines myocytes can uncouple from their neighbours. Electrical conduction may still occur via fibroblasts that not only couple proximal myocytes but can also couple otherwise unconnected regions. We hypothesise that such coupling can alter conduction between myocytes via introduction of delays or by initiation of premature stimuli that can potentially result in reentry or conduction blocks. To test this hypothesis we have developed several 2-cell motifs and investigated the effect of fibroblast mediated electrical coupling between uncoupled myocytes. We have identified various regimes of myocyte behaviour that depend on the strength of gap-junctional conductance, connection topology, and parameters of the myocyte and fibroblast models. These motifs are useful in developing a mechanistic understanding of long-distance coupling on myocyte dynamics and enable the characterisation of interaction between different features such as myocyte and fibroblast properties, coupling strengths and pacing period. They are computationally inexpensive and allow for incorporation of spatial effects such as conduction velocity. They provide a framework for constructing scar tissue boundaries and enable linking of cellular level interactions with scar induced arrhythmia.

The role of computational methods in cardiovascular medicine: a narrative review

Background and Objective

Computational models of the cardiovascular system allow for a detailed and quantitative investigation of both physiological and pathological conditions, thanks to their ability to combine clinical-possibly patient-specific-data with physical knowledge of the processes underlying the heart function. These models have been increasingly employed in clinical practice to understand pathological mechanisms and their progression, design medical devices, support clinicians in improving therapies. Hinging upon a long-year experience in cardiovascular modeling, we have recently constructed a computational multi-physics and multi-scale integrated model of the heart for the investigation of its physiological function, the analysis of pathological conditions, and to support clinicians in both diagnosis and treatment planning. This narrative review aims to systematically discuss the role that such model had in addressing specific clinical questions, and how further impact of computational models on clinical practice are envisaged.

Methods

We developed computational models of the physical processes encompassed by the heart function (electrophysiology, electrical activation, force generation, mechanics, blood flow dynamics, valve dynamics, myocardial perfusion) and of their inherently strong coupling. To solve the equations of such models, we devised advanced numerical methods, implemented in a flexible and highly efficient software library. We also developed computational procedures for clinical data post-processing-like the reconstruction of the heart geometry and motion from diagnostic images-and for their integration into computational models.

Key Content and Findings

Our integrated computational model of the heart function provides non-invasive measures of indicators characterizing the heart function and dysfunctions, and sheds light on its underlying processes and their coupling. Moreover, thanks to the close collaboration with several clinical partners, we addressed specific clinical questions on pathological conditions, such as arrhythmias, ventricular dyssynchrony, hypertrophic cardiomyopathy, degeneration of prosthetic valves, and the way coronavirus disease 2019 (COVID-19) infection may affect the cardiac function. In multiple cases, we were also able to provide quantitative indications for treatment.

Conclusions

Computational models provide a quantitative and detailed tool to support clinicians in patient care, which can enhance the assessment of cardiac diseases, the prediction of the development of pathological conditions, and the planning of treatments and follow-up tests.

Computational Modeling and Simulation of the Fibrotic Human Atria

Patient-specific modeling of atrial electrical activity enables the execution of simulations that can provide mechanistic insights and provide novel solutions to vexing clinical problems. The geometry and fibrotic remodeling of the heart can be reconstructed from clinical-grade medical scans and used to inform personalized models with detail incorporated at the cell- and tissue-scale to represent changes in image-identified diseased regions. Here, we provide a rubric for the reconstruction of realistic atrial models from pre-segmented 3D renderings of the left atrium with fibrotic tissue regions delineated, which are the output from clinical-grade systems for quantifying fibrosis. We then provide a roadmap for using those models to carry out patient-specific characterization of the fibrotic substrate in terms of its potential to harbor reentrant drivers via cardiac electrophysiology simulations.

Up digital and personal: How heart digital twins can transform heart patient care

Precision medicine is the vision of health care where therapy is tailored to each patient. As part of this vision, digital twinning technology promises to deliver a digital representation of organs or even patients by using tools capable of simulating personal health conditions and predicting patient or disease trajectories on the basis of relationships learned both from data and from biophysics knowledge. Such virtual replicas would update themselves with data from monitoring devices and medical tests and assessments, reflecting dynamically the changes in our health conditions and the responses to treatment. In precision cardiology, the concepts and initial applications of heart digital twins have slowly been gaining popularity and the trust of the clinical community. In this article, we review the advancement in heart digital twinning and its initial translation to the management of heart rhythm disorders.

Impact of noise on the instability of spiral waves in stochastic 2D mathematical models of human atrial fibrillation

Sustained spiral waves, also known as rotors, are pivotal mechanisms in persistent atrial fibrillation (AF). Stochasticity is inevitable in nonlinear biological systems such as the heart; however, it is unclear how noise affects the instability of spiral waves in human AF. This study presents a stochastic two-dimensional mathematical model of human AF and explores how Gaussian white noise affects the instability of spiral waves. In homogeneous tissue models, Gaussian white noise may lead to spiral-wave meandering and wavefront break-up. As the noise intensity increases, the spatial dispersion of phase singularity (PS) points increases. This finding indicates the potential AF-protective effects of cardiac system stochasticity by destabilizing the rotors. By contrast, Gaussian white noise is unlikely to affect the spiral-wave instability in the presence of localized scar or fibrosis regions. The PS points are located at the boundary or inside the scar/fibrosis regions. Localized scar or fibrosis may play a pivotal role in stabilizing spiral waves regardless of the presence of noise. This study suggests that fibrosis and scars are essential for stabilizing the rotors in stochastic mathematical models of AF. Further patient-derived realistic modeling studies are required to confirm the role of scar/fibrosis in AF pathophysiology.

lifex-ep: a robust and efficient software for cardiac electrophysiology simulations

BACKGROUND

Simulating the cardiac function requires the numerical solution of multi-physics and multi-scale mathematical models. This underscores the need for streamlined, accurate, and high-performance computational tools. Despite the dedicated endeavors of various research teams, comprehensive and user-friendly software programs for cardiac simulations, capable of accurately replicating both normal and pathological conditions, are still in the process of achieving full maturity within the scientific community.

RESULTS

This work introduces [Formula: see text]-ep, a publicly available software for numerical simulations of the electrophysiology activity of the cardiac muscle, under both normal and pathological conditions. [Formula: see text]-ep employs the monodomain equation to model the heart's electrical activity. It incorporates both phenomenological and second-generation ionic models. These models are discretized using the Finite Element method on tetrahedral or hexahedral meshes. Additionally, [Formula: see text]-ep integrates the generation of myocardial fibers based on Laplace-Dirichlet Rule-Based Methods, previously released in Africa et al., 2023, within [Formula: see text]-fiber. As an alternative, users can also choose to import myofibers from a file. This paper provides a concise overview of the mathematical models and numerical methods underlying [Formula: see text]-ep, along with comprehensive implementation details and instructions for users. [Formula: see text]-ep features exceptional parallel speedup, scaling efficiently when using up to thousands of cores, and its implementation has been verified against an established benchmark problem for computational electrophysiology. We showcase the key features of [Formula: see text]-ep through various idealized and realistic simulations conducted in both normal and pathological scenarios. Furthermore, the software offers a user-friendly and flexible interface, simplifying the setup of simulations using self-documenting parameter files.

CONCLUSIONS

[Formula: see text]-ep provides easy access to cardiac electrophysiology simulations for a wide user community. It offers a computational tool that integrates models and accurate methods for simulating cardiac electrophysiology within a high-performance framework, while maintaining a user-friendly interface. [Formula: see text]-ep represents a valuable tool for conducting in silico patient-specific simulations.

Comparing Inducibility of Re-Entrant Arrhythmia in Patient-Specific Computational Models to Clinical Atrial Fibrillation Phenotypes

BACKGROUND

Computational models of fibrosis-mediated, re-entrant left atrial (LA) arrhythmia can identify possible substrate for persistent atrial fibrillation (AF) ablation. Contemporary models use a 1-size-fits-all approach to represent electrophysiological properties, limiting agreement between simulations and patient outcomes.

OBJECTIVES

The goal of this study was to test the hypothesis that conduction velocity (ϴ) modulation in persistent AF models can improve simulation agreement with clinical arrhythmias.

METHODS

Patients with persistent AF (n = 37) underwent ablation and were followed up for ≥2 years to determine post-ablation outcomes: AF, atrial flutter (AFL), or no recurrence. Patient-specific LA models (n = 74) were constructed using pre-ablation and ≥90 days' post-ablation magnetic resonance imaging data. Simulated pacing gauged in silico arrhythmia inducibility due to AF-like rotors or AFL-like macro re-entrant tachycardias. A physiologically plausible range of ϴ values (±10 or 20% vs. baseline) was tested, and model/clinical agreement was assessed.

RESULTS

Fifteen (41%) patients had a recurrence with AF and 6 (16%) with AFL. Arrhythmia was induced in 1,078 of 5,550 simulations. Using baseline ϴ, model/clinical agreement was 46% (34 of 74 models), improving to 65% (48 of 74) when any possible ϴ value was used (McNemar's test, P = 0.014). ϴ modulation improved model/clinical agreement in both pre-ablation and post-ablation models. Pre-ablation model/clinical agreement was significantly greater for patients with extensive LA fibrosis (>17.2%) and an elevated body mass index (>32.0 kg/m2).

CONCLUSIONS

Simulations in persistent AF models show a 41% relative improvement in model/clinical agreement when ϴ is modulated. Patient-specific calibration of ϴ values could improve model/clinical agreement and model usefulness, especially in patients with higher body mass index or LA fibrosis burden. This could ultimately facilitate better personalized modeling, with immediate clinical implications.

AugmentA: Patient-specific augmented atrial model generation tool

Digital twins of patients' hearts are a promising tool to assess arrhythmia vulnerability and to personalize therapy. However, the process of building personalized computational models can be challenging and requires a high level of human interaction. We propose a patient-specific Augmented Atria generation pipeline (AugmentA) as a highly automated framework which, starting from clinical geometrical data, provides ready-to-use atrial personalized computational models. AugmentA identifies and labels atrial orifices using only one reference point per atrium. If the user chooses to fit a statistical shape model to the input geometry, it is first rigidly aligned with the given mean shape before a non-rigid fitting procedure is applied. AugmentA automatically generates the fiber orientation and finds local conduction velocities by minimizing the error between the simulated and clinical local activation time (LAT) map. The pipeline was tested on a cohort of 29 patients on both segmented magnetic resonance images (MRI) and electroanatomical maps of the left atrium. Moreover, the pipeline was applied to a bi-atrial volumetric mesh derived from MRI. The pipeline robustly integrated fiber orientation and anatomical region annotations in 38.4 ± 5.7 s. In conclusion, AugmentA offers an automated and comprehensive pipeline delivering atrial digital twins from clinical data in procedural time.

What determines the optimal pharmacological treatment of atrial fibrillation? Insights from in silico trials in 800 virtual atria

The best pharmacological treatment for each atrial fibrillation (AF) patient is unclear. We aim to exploit AF simulations in 800 virtual atria to identify key patient characteristics that guide the optimal selection of anti-arrhythmic drugs. The virtual cohort considered variability in electrophysiology and low voltage areas (LVA) and was developed and validated against experimental and clinical data from ionic currents to ECG. AF sustained in 494 (62%) atria, with large inward rectifier K+ current (IK1 ) and Na+ /K+ pump (INaK ) densities (IK1 0.11 ± 0.03 vs. 0.07 ± 0.03 S mF-1 ; INaK 0.68 ± 0.15 vs. 0.38 ± 26 S mF-1 ; sustained vs. un-sustained AF). In severely remodelled left atrium, with LVA extensions of more than 40% in the posterior wall, higher IK1 (median density 0.12 ± 0.02 S mF-1 ) was required for AF maintenance, and rotors localized in healthy right atrium. For lower LVA extensions, rotors could also anchor to LVA, in atria presenting short refractoriness (median L-type Ca2+ current, ICaL , density 0.08 ± 0.03 S mF-1 ). This atrial refractoriness, modulated by ICaL and fast Na+ current (INa ), determined pharmacological treatment success for both small and large LVA. Vernakalant was effective in atria presenting long refractoriness (median ICaL density 0.13 ± 0.05 S mF-1 ). For short refractoriness, atria with high INa (median density 8.92 ± 2.59 S mF-1 ) responded more favourably to amiodarone than flecainide, and the opposite was found in atria with low INa (median density 5.33 ± 1.41 S mF-1 ). In silico drug trials in 800 human atria identify inward currents as critical for optimal stratification of AF patient to pharmacological treatment and, together with the left atrial LVA extension, for accurately phenotyping AF dynamics. KEY POINTS: Atrial fibrillation (AF) maintenance is facilitated by small L-type Ca2+ current (ICaL ) and large inward rectifier K+ current (IK1 ) and Na+ /K+ pump. In severely remodelled left atrium, with low voltage areas (LVA) covering more than 40% of the posterior wall, sustained AF requires higher IK1 and rotors localize in healthy right atrium. For lower LVA extensions, rotors can also anchor to LVA, if the atria present short refractoriness (low ICaL ) Vernakalant is effective in atria presenting long refractoriness (high ICaL ). For short refractoriness, atria with fast Na+ current (INa ) up-regulation respond more favourably to amiodarone than flecainide, and the opposite is found in atria with low INa . The inward currents (ICaL and INa ) are critical for optimal stratification of AF patient to pharmacological treatment and, together with the left atrial LVA extension, for accurately phenotyping AF dynamics.

Explainable Machine Learning to Predict Anchored Reentry Substrate Created by Persistent Atrial Fibrillation Ablation in Computational Models

Background Postablation arrhythmia recurrence occurs in ~40% of patients with persistent atrial fibrillation. Fibrotic remodeling exacerbates arrhythmic activity in persistent atrial fibrillation and can play a key role in reentrant arrhythmia, but emergent interaction between nonconductive ablation-induced scar and native fibrosis (ie, residual fibrosis) is poorly understood. Methods and Results We conducted computational simulations in pre- and postablation left atrial models reconstructed from late gadolinium enhanced magnetic resonance imaging scans to test the hypothesis that ablation in patients with persistent atrial fibrillation creates new substrate conducive to recurrent arrhythmia mediated by anchored reentry. We trained a random forest machine learning classifier to accurately pinpoint specific nonconductive tissue regions (ie, areas of ablation-delivered scar or vein/valve boundaries) with the capacity to serve as substrate for anchored reentry-driven recurrent arrhythmia (area under the curve: 0.91±0.03). Our analysis suggests there is a distinctive nonconductive tissue pattern prone to serving as arrhythmogenic substrate in postablation models, defined by a specific size and proximity to residual fibrosis. Conclusions Overall, this suggests persistent atrial fibrillation ablation transforms substrate that favors functional reentry (ie, rotors meandering in excitable tissue) into an arrhythmogenic milieu more conducive to anchored reentry. Our work also indicates that explainable machine learning and computational simulations can be combined to effectively probe mechanisms of recurrent arrhythmia.

Recurrence quantification analysis for fine-scale characterisation of arrhythmic patterns in cardiac tissue

This paper uses recurrence quantification analysis (RQA) combined with entropy measures and organization indices to characterize arrhythmic patterns and dynamics in computer simulations of cardiac tissue. We performed different simulations of cardiac tissues of sizes comparable to the human heart atrium. In these simulations, we observed four classic arrhythmic patterns: a spiral wave anchored to a highly fibrotic region resulting in sustained re-entry, a meandering spiral wave, fibrillation, and a spiral wave anchored to a scar region that breaks up into wavelets away from the main rotor. A detailed analysis revealed that, within the same simulation, maps of RQA metrics could differentiate regions with regular AP propagation from ones with chaotic activity. In particular, the combination of two RQA metrics, the length of the longest diagonal string of recurrence points and the mean length of diagonal lines, was able to identify the location of rotor tips, which are the active elements that maintain spiral waves and fibrillation. By proposing low-dimensional models based on the mean value and spatial correlation of metrics calculated from membrane potential time series, we identify RQA-based metrics that successfully separate the four different types of cardiac arrhythmia into distinct regions of the feature space, and thus might be used for automatic classification, in particular distinguishing between fibrillation driven by self-sustaining chaos and that created by a persistent rotor and wavebreak. We also discuss the practical applicability of such an approach.

Personalized biomechanical insights in atrial fibrillation: opportunities & challenges

INTRODUCTION

Atrial fibrillation (AF) is an increasingly prevalent and significant worldwide health problem. Manifested as an irregular atrial electrophysiological activation, it is associated with many serious health complications. AF affects the biomechanical function of the heart as contraction follows the electrical activation, subsequently leading to reduced blood flow. The underlying mechanisms behind AF are not fully understood, but it is known that AF is highly correlated with the presence of atrial fibrosis, and with a manifold increase in risk of stroke.

AREAS COVERED

In this review, we focus on biomechanical aspects in atrial fibrillation, current and emerging use of clinical images, and personalized computational models. We also discuss how these can be used to provide patient-specific care.

EXPERT OPINION

Understanding the connection betweenatrial fibrillation and atrial remodeling might lead to valuable understanding of stroke and heart failure pathophysiology. Established and emerging imaging modalities can bring us closer to this understanding, especially with continued advancements in processing accuracy, reproducibility, and clinical relevance of the associated technologies. Computational models of cardiac electromechanics can be used to glean additional insights on the roles of AF and remodeling in heart function.

Machine Learning and the Conundrum of Stroke Risk Prediction

Stroke is a leading cause of death worldwide. With escalating healthcare costs, early non-invasive stroke risk stratification is vital. The current paradigm of stroke risk assessment and mitigation is focused on clinical risk factors and comorbidities. Standard algorithms predict risk using regression-based statistical associations, which, while useful and easy to use, have moderate predictive accuracy. This review summarises recent efforts to deploy machine learning (ML) to predict stroke risk and enrich the understanding of the mechanisms underlying stroke. The surveyed body of literature includes studies comparing ML algorithms with conventional statistical models for predicting cardiovascular disease and, in particular, different stroke subtypes. Another avenue of research explored is ML as a means of enriching multiscale computational modelling, which holds great promise for revealing thrombogenesis mechanisms. Overall, ML offers a new approach to stroke risk stratification that accounts for subtle physiologic variants between patients, potentially leading to more reliable and personalised predictions than standard regression-based statistical associations.

SGLT2i alleviates epicardial adipose tissue inflammation by modulating ketone body-glyceraldehyde-3-phosphate dehydrogenase malonylation pathway

AIMS

Inflammation in the epicardial adipose tissue (EAT) is a contributor to atrial fibrillation. Studies have reported that sodium glucose co-transporter 2 inhibitor (SGLT2i) can alleviate EAT inflammation. However, the mechanism remains elusive. This study aims to investigate the molecular mechanism of SGLT2i in reducing EAT inflammation and to explore the effects of SGLT2i on atrial fibrosis in atrial fibrillation.

METHODS

Sprague-Dawley rats were injected with angiotensin II to induce atrial fibrillation and randomly assigned to receive SGLT2i ( n  = 6) or vehicle ( n  = 6). Macrophages (RAW264.7) were treated with ketone bodies; ACC1 knockdown/overexpression and malonyl-CoA overexpression were performed in vitro . The levels of inflammatory cytokines, ACC1, and malonyl-CoA were examined by ELISA. GAPDH malonylation was measured by co-immunoprecipitation.

RESULTS

In atrial fibrillation rats, SGLT2i increased the ketone body levels and decreased the expression of ACC1 and alleviated EAT inflammation and atrial fibrosis. In RAW264.7 cells, ketone bodies decreased the levels of ACC1, malonyl-CoA, and GAPDH malonylation, accompanied by reduced inflammatory cytokines. ACC1 knockdown decreased the expression of malonyl-CoA and GAPDH malonylation and alleviated lipopolysaccharide (LPS)-induced macrophage inflammation; these effects were inhibited by malonyl-CoA overexpression. Furthermore, the protective effects of ketone bodies on macrophage inflammation were abrogated by ACC1 overexpression.

CONCLUSION

SGLT2i alleviates EAT inflammation by reducing GAPDH malonylation via downregulating the expression of ACC1 through increasing ketone bodies, thus attenuating atrial fibrosis.

microRNAs as Biomarkers of Endothelial Dysfunction and Therapeutic Target in the Pathogenesis of Atrial Fibrillation

The pathophysiology of atrial fibrillation (AF) may involve atrial fibrosis/remodeling and dysfunctional endothelial activities. Despite the currently available treatment approaches, the progression of AF, its recurrence rate, and the high mortality risk of related complications underlay the need for more advanced prognostic and therapeutic strategies. There is increasing attention on the molecular mechanisms controlling AF onset and progression points to the complex cell to cell interplay that triggers fibroblasts, immune cells and myofibroblasts, enhancing atrial fibrosis. In this scenario, endothelial cell dysfunction (ED) might play an unexpected but significant role. microRNAs (miRNAs) regulate gene expression at the post-transcriptional level. In the cardiovascular compartment, both free circulating and exosomal miRNAs entail the control of plaque formation, lipid metabolism, inflammation and angiogenesis, cardiomyocyte growth and contractility, and even the maintenance of cardiac rhythm. Abnormal miRNAs levels may indicate the activation state of circulating cells, and thus represent a specific read-out of cardiac tissue changes. Although several unresolved questions still limit their clinical use, the ease of accessibility in biofluids and their prognostic and diagnostic properties make them novel and attractive biomarker candidates in AF. This article summarizes the most recent features of AF associated with miRNAs and relates them to potentially underlying mechanisms.

Development and Validation of a Postoperative Prognostic Nomogram to Predict Recurrence in Patients with Persistent Atrial Fibrillation: A Retrospective Cohort Study

Background: Patients with persistent atrial fibrillation (PsAF) have a high risk of recurrence after catheter radiofrequency ablation. Nevertheless, no effective prognostic tools have been developed to identify these high-risk patients to date. This study sought to develop and validate a simple linear predictive model for predicting postoperative recurrence in patients with PsAF.

Methods: From June 2013 to June 2021, patients with PsAF admitted to our hospital were enrolled in this single-center, retrospective, observational study. The characteristics substantially associated with recurrence in patients with PsAF were screened through univariate and multivariate logistic regression analysis. The receiver operating characteristic curve was used to assess the predictive significance of the nomogram model after nomogram development. Furthermore, to assess the clinical value of the nomogram, we performed calibration curve and decision curve analyses.

Results: A total of 209 patients were included in the study, 42 (20.10%) of whom were monitored up to 1 year for recurrent AF. The duration of AF episodes, left atrial diameter, BMI, CKMB, and alcohol consumption were found to be independent risk factors (P<0.05) and were integrated into the nomogram model development. The area under the curve was 0.895, the sensitivity was 93.3%, and the specificity was 71.4%, thus indicating the model’s excellent predictive ability. The C-index of the predictive nomogram model was 0.906. Calibration curve and decision curve analyses further revealed that the model had robust prediction and strong discrimination ability.

Conclusion: This simple, practical, and innovative nomogram can help clinicians in evaluation of the risk of PsAF recurrence after catheter ablation, thus facilitating preoperative evaluation, postoperative monitoring and ultimately the construction of more personalized therapeutic protocols.

Treatment Planning for Atrial Fibrillation Using Patient-Specific Models Showing the Importance of Fibrillatory-Areas

Computational models have made it possible to study the effect of fibrosis and scar on atrial fibrillation (AF) and plan future personalized treatments. Here, we study the effect of area available for fibrillatory waves to sustain AF. Then we use it to plan for AF ablation to improve procedural outcomes. CARPentry was used to create patient-specific models to determine the association between the size of residual contiguous areas available for AF wavefronts to propagate and sustain AF [fibrillatory area (FA)] after ablation with procedural outcomes. The FA was quantified in a novel manner accounting for gaps in ablation lines. We selected 30 persistent AF patients with known ablation outcomes. We divided the atrial surface into five areas based on ablation scar pattern and anatomical landmarks and calculated the FAs. We validated the models based on clinical outcomes and suggested future ablation lines that minimize the FAs and terminate rotor activities in simulations. We also simulated the effects of three common antiarrhythmic drugs. In the patient-specific models, the predicted arrhythmias matched the clinical outcomes in 25 of 30 patients (accuracy 83.33%). The average largest FA (FAmax) in the recurrence group was 8517 ± 1444 vs. 6772 ± 1531 mm2 in the no recurrence group (p < 0.004). The final FAs after adding the suggested ablation lines in the AF recurrence group reduced the average FAmax from 8517 ± 1444 to 6168 ± 1358 mm2 (p < 0.001) and stopped the sustained rotor activity. Simulations also correctly anticipated the effect of antiarrhythmic drugs in 5 out of 6 patients who used drug therapy post unsuccessful ablation (accuracy 83.33%). Sizes of FAs available for AF wavefronts to propagate are important determinants for ablation outcomes. FA size in combination with computational simulations can be used to direct ablation in persistent AF to minimize the critical mass required to sustain recurrent AF.

Quantifying the impact of shape uncertainty on predicted arrhythmias

BACKGROUND

Personalised computer models are increasingly used to diagnose cardiac arrhythmias and tailor treatment. Patient-specific models of the left atrium are often derived from pre-procedural imaging of anatomy and fibrosis. These images contain noise that can affect simulation predictions. There are few computationally tractable methods for propagating uncertainties from images to clinical predictions.

METHOD

We describe the left atrium anatomy using our Bayesian shape model that captures anatomical uncertainty in medical images and has been validated on 63 independent clinical images. This algorithm describes the left atrium anatomy using Nmodes=15 principal components, capturing 95% of the shape variance and calculated from 70 clinical cardiac magnetic resonance (CMR) images. Latent variables encode shape uncertainty: we evaluate their posterior distribution for each new anatomy. We assume a normally distributed prior. We use the unscented transform to sample from the posterior shape distribution. For each sample, we assign the local material properties of the tissue using the projection of late gadolinium enhancement CMR (LGE-CMR) onto the anatomy to estimate local fibrosis. To test which activation patterns an atrium can sustain, we perform an arrhythmia simulation for each sample. We consider 34 possible outcomes (31 macro-re-entries, functional re-entry, atrial fibrillation, and non-sustained arrhythmia). For each sample, we determine the outcome by comparing pre- and post-ablation activation patterns following a cross-field stimulus.

RESULTS

We create patient-specific atrial electrophysiology models of ten patients. We validate the mean and standard deviation maps from the unscented transform with the same statistics obtained with 12,000 Monte Carlo (ground truth) samples. We found discrepancies <3% and <2% for the mean and standard deviation for fibrosis burden and activation time, respectively. For each patient case, we then compare the predicted outcome from a model built on the clinical data (deterministic approach) with the probability distribution obtained from the simulated samples. We found that the deterministic approach did not predict the most likely outcome in 80% of the cases. Finally, we estimate the influence of each source of uncertainty independently. Fixing the anatomy to the posterior mean and maintaining uncertainty in fibrosis reduced the prediction of self-terminating arrhythmias from ≃14% to ≃7%. Keeping the fibrosis fixed to the sample mean while retaining uncertainty in shape decreased the prediction of substrate-driven arrhythmias from ≃33% to ≃18% and increased the prediction of macro-re-entries from ≃54% to ≃68%.

CONCLUSIONS

We presented a novel method for propagating shape uncertainty in atrial models through to uncertainty in numerical simulations. The algorithm takes advantage of the unscented transform to compute the output distribution of the outcomes. We validated the unscented transform as a viable sampling strategy to deal with anatomy uncertainty. We then showed that the prediction computed with a deterministic model does not always coincide with the most likely outcome. Finally, we found that shape uncertainty affects the predictions of macro-re-entries, while fibrosis uncertainty affects the predictions of functional re-entries.

Contribution of atrial myofiber architecture to atrial fibrillation

BACKGROUND

The role of fiber orientation on a global chamber level in sustaining atrial fibrillation (AF) is unknown. The goal of this study was to correlate the fiber direction derived from Diffusion Tensor Imaging (DTI) with AF inducibility.

METHODS

Transgenic goats with cardiac-specific overexpression of constitutively active TGF-β1 (n = 14) underwent AF inducibility testing by rapid pacing in the left atrium. We chose a minimum of 10 minutes of sustained AF as a cut-off for AF inducibility. Explanted hearts underwent DTI to determine the fiber direction. Using tractography data, we clustered, visualized, and quantified the fiber helix angles in 8 different regions of the left atrial wall using two reference vectors defined based on anatomical landmarks.

RESULTS

Sustained AF was induced in 7 out of 14 goats. The mean helix fiber angles in 7 out of 8 selected regions were statistically different (P-Value < 0.05) in the AF inducible group. The average fractional anisotropy (FA) and the mean diffusivity (MD) were similar in the two groups with FA of 0.32±0.08 and MD of 8.54±1.72 mm2/s in the non-inducible group and FA of 0.31±0.05 (P-value = 0.90) and MD of 8.68±1.60 mm2/s (P-value = 0.88) in the inducible group.

CONCLUSIONS

DTI based fiber direction shows significant variability across subjects with a significant difference between animals that are AF inducible versus animals that are not inducible. Fiber direction might be contributing to the initiation and sustaining of AF, and its role needs to be investigated further.

Homogenisation for the monodomain model in the presence of microscopic fibrotic structures

Computational models in cardiac electrophysiology are notorious for long runtimes, restricting the numbers of nodes and mesh elements in the numerical discretisations used for their solution. This makes it particularly challenging to incorporate structural heterogeneities on small spatial scales, preventing a full understanding of the critical arrhythmogenic effects of conditions such as cardiac fibrosis. In this work, we explore the technique of homogenisation by volume averaging for the inclusion of non-conductive micro-structures into larger-scale cardiac meshes with minor computational overhead. Importantly, our approach is not restricted to periodic patterns, enabling homogenised models to represent, for example, the intricate patterns of collagen deposition present in different types of fibrosis. We first highlight the importance of appropriate boundary condition choice for the closure problems that define the parameters of homogenised models. Then, we demonstrate the technique's ability to correctly upscale the effects of fibrotic patterns with a spatial resolution of 10 µm into much larger numerical mesh sizes of 100- 250 µm . The homogenised models using these coarser meshes correctly predict critical pro-arrhythmic effects of fibrosis, including slowed conduction, source/sink mismatch, and stabilisation of re-entrant activation patterns. As such, this approach to homogenisation represents a significant step towards whole organ simulations that unravel the effects of microscopic cardiac tissue heterogeneities.

Atrial fibrillation: Insights from animal models, computational modeling, and clinical studies

Atrial fibrillation (AF) is the most common human arrhythmia, affecting millions of patients worldwide. A combination of risk factors and comorbidities results in complex atrial remodeling, which increases AF vulnerability and persistence. Insights from animal models, clinical studies, and computational modeling have advanced the understanding of the mechanisms and pathophysiology of AF. Areas of heterogeneous pathological remodeling, as well as altered electrophysiological properties, serve as a substrate for AF drivers and spontaneous activations. The complex and individualized presentation of this arrhythmia suggests that mechanisms-based personalized approaches will likely be needed to overcome current challenges in AF management. In this paper, we review the insights on the mechanisms of AF obtained from animal models and clinical studies and how computational models integrate this knowledge to advance AF clinical management. We also assess the challenges that need to be overcome to implement these mechanistic models in clinical practice.

The physics of heart rhythm disorders

The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.

Interpretable machine learning of action potential duration restitution kinetics in single-cell models of atrial cardiomyocytes

Action potential duration (APD) restitution curve and its maximal slope (Smax) reflect single cell-level dynamic instability for inducing chaotic heart rhythms. However, conventional parameter sensitivity analysis often fails to describe nonlinear relationships between ion channel parameters and electrophysiological phenotypes, such as Smax. We explored the parameter-phenotype mapping in a population of 5000 single-cell atrial cell models through interpretable machine learning (ML) approaches. Parameter sensitivity analyses could explain the linear relationships between parameters and electrophysiological phenotypes, including APD₉₀, resting membrane potential, Vmax, refractory period, and APD/calcium alternans threshold, but not for Smax. However, neural network models had better prediction performance for Smax. To interpret the ML model, we evaluated the parameter importance at the global and local levels by computing the permutation feature importance and the local interpretable model-agnostic explanations (LIME) values, respectively. Increases in I_CaL, I_NCX, and I_Kr, and decreases in I_K1, I_b,Cl, I_Kur, I_SERCA, and I_to are correlated with higher Smax values. The LIME algorithm determined that I_NaK plays a significant role in determining Smax as well as I_to and I_Kur. The atrial cardiomyocyte population was hierarchically clustered into three distinct groups based on the LIME values and the single-cell simulation confirmed that perturbations in I_NaK resulted in different behaviors of APD restitution curves in three clusters. Our combined top-down interpretable ML and bottom-up mechanistic simulation approaches uncovered the role of I_NaK in heterogeneous behaviors of Smax in the atrial cardiomyocyte population.

Rotors Drift Toward and Stabilize in Low Power Regions in Heterogeneous Models of Atrial Fibrillation

Atrial fibrillation (AF) afflicts more than 33 million people worldwide. Success of therapy strategies remains poor and better understanding of the arrhythmia and how to device more effective therapies are needed. The aim of this work is to study the role of electric power distributions in rotors and AF dynamics. For this purpose, single cell and tissue simulations were performed to study the effect of ionic currents gradients and fibrosis in rotor's drifting. The root mean square of the ionic (Pion), capacitance (Pc) and electrotonic (Pele) power was computed over action potentials. Single cell simulations were performed for different values of IK1 and ICaL and number of coupled myofibroblasts. Tissue simulations were performed in presence of IK1 and ICaL gradients and diffused fibrosis. Single cell simulations showed that Pion and Pc increased with IK1, while decreased by increasing ICaL. Increasing the number of coupled myofibroblasts reduced Pion and Pc, whereas Pele increased. Finally, in tissue simulations rotors drifted to regions with low power and anchored in regions with higher density of blunted ionic induced power gradients.

How synergy between mechanistic and statistical models is impacting research in atrial fibrillation

Atrial fibrillation (AF) with multiple complications, high morbidity and mortality, and low cure rates, has become a global public health problem. Although significant progress has been made in the treatment methods represented by anti-AF drugs and radiofrequency ablation, the therapeutic effect is not as good as expected. The reason is mainly because of our lack of understanding of AF mechanisms. This field has benefited from mechanistic and (or) statistical methodologies. Recent renewed interest in digital twin techniques by synergizing between mechanistic and statistical models has opened new frontiers in AF analysis. In the review, we briefly present findings that gave rise to the AF pathophysiology and current therapeutic modalities. We then summarize the achievements of digital twin technologies in three aspects: understanding AF mechanisms, screening anti-AF drugs and optimizing ablation strategies. Finally, we discuss the challenges that hinder the clinical application of the digital twin heart. With the rapid progress in data reuse and sharing, we expect their application to realize the transition from AF description to response prediction.

R-on-T and the initiation of reentry revisited: Integrating old and new concepts

Initiation of reentry requires 2 factors: (1) a triggering event, most commonly focal excitations such as premature ventricular complexes (PVCs); and (2) a vulnerable substrate with regional dispersion of refractoriness and/or excitability, such as occurs during the T wave of the electrocardiogram when some areas of the ventricle have repolarized and recovered excitability but others have not. When the R wave of a PVC coincides in time with the T wave of the previous beat, this timing can lead to unidirectional block and initiation of reentry, known as the R-on-T phenomenon. Classically, the PVC triggering reentry has been viewed as arising focally from 1 region and propagating into another region whose recovery is delayed, resulting in unidirectional conduction block and reentry initiation. However, more recent evidence indicates that PVCs also can arise from the T wave itself. In the latter case, the PVC initiating reentry is not a separate event from the T wave but rather is causally generated from the repolarization gradient that manifests as the T wave. We call the former an "R-to-T" mechanism and the latter an "R-from-T" mechanism, which are initiation mechanisms distinct from each other. Both are important components of the R-on-T phenomenon and need to be taken into account when designing antiarrhythmic strategies. Strategies targeting suppression of triggers alone or vulnerable substrate alone may be appropriate in some instances but not in others. Preventing R-from-T arrhythmias requires suppressing the underlying dynamic tissue instabilities responsible for producing both triggers and substrate vulnerability simultaneously. The same principles are likely to apply to supraventricular arrhythmias.

Ablation and antiarrhythmic drug effects on PITX2+/− deficient atrial fibrillation: A computational modeling study

Introduction

Atrial fibrillation (AF) is a heritable disease, and the paired-like homeodomain transcription factor 2 (PITX2) gene is highly associated with AF. We explored the differences in the circumferential pulmonary vein isolation (CPVI), which is the cornerstone procedure for AF catheter ablation, additional high dominant frequency (DF) site ablation, and antiarrhythmic drug (AAD) effects according to the patient genotype (wild-type and PITX2+/− deficient) using computational modeling.

Methods

We included 25 patients with AF (68% men, 59.8 ± 9.8 years of age, 32% paroxysmal AF) who underwent AF catheter ablation to develop a realistic computational AF model. The ion currents for baseline AF and the amiodarone, dronedarone, and flecainide AADs according to the patient genotype (wild type and PITX2+/− deficient) were defined by relevant publications. We tested the virtual CPVI (V-CPVI) with and without DF ablation (±DFA) and three virtual AADs (V-AADs, amiodarone, dronedarone, and flecainide) and evaluated the AF defragmentation rates (AF termination or changes to regular atrial tachycardia (AT), DF, and maximal slope of the action potential duration restitution curves (Smax), which indicates the vulnerability of wave-breaks.

Results

At the baseline AF, mean DF (p = 0.003), and Smax (p &lt; 0.001) were significantly lower in PITX2+/− deficient patients than wild-type patients. In the overall AF episodes, V-CPVI (±DFA) resulted in a higher AF defragmentation relative to V-AADs (65 vs. 42%, p &lt; 0.001) without changing the DF or Smax. Although a PITX2+/− deficiency did not affect the AF defragmentation rate after the V-CPVI (±DFA), V-AADs had a higher AF defragmentation rate (p = 0.014), lower DF (p &lt; 0.001), and lower Smax (p = 0.001) in PITX2+/− deficient AF than in wild-type patients. In the clinical setting, the PITX2+/− genetic risk score did not affect the AF ablation rhythm outcome (Log-rank p = 0.273).

Conclusion

Consistent with previous clinical studies, the V-CPVI had effective anti-AF effects regardless of the PITX2 genotype, whereas V-AADs exhibited more significant defragmentation or wave-dynamic change in the PITX2+/− deficient patients.

The Atrium in Atrial Fibrillation - A Clinical Review on How to Manage Atrial Fibrotic Substrates

Atrial fibrillation (AF) is the most common sustained arrhythmia in the population and is associated with a significant clinical and economic burden. Rigorous assessment of the presence and degree of an atrial arrhythmic substrate is essential for determining treatment options, predicting long-term success after catheter ablation, and as a substrate critical in the pathophysiology of atrial thrombogenesis. Catheter ablation of AF has developed into an essential rhythm-control strategy. Nowadays is one of the most common cardiac ablation procedures performed worldwide, with its success inversely related to the extent of atrial structural disease. Although atrial substrate evaluation remains complex, several diagnostic resources allow for a more comprehensive assessment and quantification of the extent of left atrial structural remodeling and the presence of atrial fibrosis. In this review, we summarize the current knowledge on the pathophysiology, etiology, and electrophysiological aspects of atrial substrates promoting the development of AF. We also describe the risk factors for its development and how to diagnose its presence using imaging, electrocardiograms, and electroanatomic voltage mapping. Finally, we discuss recent data regarding fibrosis biomarkers that could help diagnose atrial fibrotic substrates.

Serum periostin as a predictor of early recurrence of atrial fibrillation after catheter ablation

Catheter ablation is an effective method of rhythm therapy for atrial fibrillation (AF). AF recurrence is a common problem after catheter ablation. The aim of this study was to investigate influence factors of early recurrence after catheter ablation for AF. One hundred and three consecutive patients with AF were enrolled and underwent catheter ablation. Venous blood (Marked as A) was collected before ablation and left atrial blood (Marked as B) was collected after successful atrial septal puncture to detect serum periostin. After 3 months of follow-up, statistical analysis was made based on the recurrence of AF. 27 (26.2%) patients had a recurrence of atrial arrhythmia after catheter ablation. Patients with recurrent atrial arrhythmia had a larger left atrial volume (162.31 ± 47.76 vs. 141.98 ± 41.64,p = 0.039), and higher serum periostin levels (periostin A. 99.71 ± 16.475 vs. 90.36 ± 13.63, p = 0.005; periostin B. 103.95 ± 13.09 vs. 94.46 ± 15.85, p = 0.006) compared with the non-recurrent group. The numbers of patients with left atrial low-voltage areas (LVAs) were more in the recurrence group (p < 0.001). Left atrial volume, serum periostin and left atrial LVAs were included in univariate and multivariate COX regression analysis. It showed that left atrial LVAs (HR3.81; 95% CI 1.54 to 9.44; p = 0.004) and serum periostin A (HR1.07; 95% CI 1.02 to1.13; p = 0.008) were the independent predictors of AF recurrence. The cut-off value of serum periostin A was 87.95 ng/ ml (AUC, 0.681; sensitivity 88.9% and specificity 53.9%). Kaplan-Meier survival curve showed that the recurrence rate of AF was higher in patients with left atrial LVAs and higher serum periostin. The venous serum periostin level and left atrial LVAs were independent predictors of early recurrence of AF after catheter ablation.

Graph-based homogenisation for modelling cardiac fibrosis

Fibrosis, the excess of extracellular matrix, can affect, and even block, propagation of action potential in cardiac tissue. This can result in deleterious effects on heart function, but the nature and severity of these effects depend strongly on the localisation of fibrosis and its by-products in cardiac tissue, such as collagen scar formation. Computer simulation is an important means of understanding the complex effects of fibrosis on activation patterns in the heart, but concerns of computational cost place restrictions on the spatial resolution of these simulations. In this work, we present a novel numerical homogenisation technique that uses both Eikonal and graph approaches to allow fine-scale heterogeneities in conductivity to be incorporated into a coarser mesh. Homogenisation achieves this by deriving effective conductivity tensors so that a coarser mesh can then be used for numerical simulation. By taking a graph-based approach, our homogenisation technique functions naturally on irregular grids and does not rely upon any assumptions of periodicity, even implicitly. We present results of action potential propagation through fibrotic tissue in two dimensions that show the graph-based homogenisation technique is an accurate and effective way to capture fine-scale domain information on coarser meshes in the context of sharp-fronted travelling waves of activation. As test problems, we consider excitation propagation in tissue with diffuse fibrosis and through a tunnel-like structure designed to test homogenisation, interaction of an excitation wave with a scar region, and functional re-entry.

A Review of Healthy and Fibrotic Myocardium Microstructure Modeling and Corresponding Intracardiac Electrograms

Computational simulations of cardiac electrophysiology provide detailed information on the depolarization phenomena at different spatial and temporal scales. With the development of new hardware and software, in silico experiments have gained more importance in cardiac electrophysiology research. For plane waves in healthy tissue, in vivo and in silico electrograms at the surface of the tissue demonstrate symmetric morphology and high peak-to-peak amplitude. Simulations provided insight into the factors that alter the morphology and amplitude of the electrograms. The situation is more complex in remodeled tissue with fibrotic infiltrations. Clinically, different changes including fractionation of the signal, extended duration and reduced amplitude have been described. In silico, numerous approaches have been proposed to represent the pathological changes on different spatial and functional scales. Different modeling approaches can reproduce distinct subsets of the clinically observed electrogram phenomena. This review provides an overview of how different modeling approaches to incorporate fibrotic and structural remodeling affect the electrogram and highlights open challenges to be addressed in future research.

Bi‐atrial characterization of the electrical substrate in patients with atrial fibrillation

Background

Methods

Results

Conclusion

Structural Remodeling and Rotational Activity in Persistent/Long-Lasting Atrial Fibrillation: Gender-Effect Differences and Impact on Post-ablation Outcome

Background

Structural and post-ablation gender differences are reported in atrial fibrillation (AF). We analyzed the gender differences in structural remodeling and AF mechanisms in patients with persistent/long-lasting AF who underwent wide area circumferential pulmonary vein isolation (WACPVI).

Materials and Methods

Ultra-high-density mapping was used to study atrial remodeling and AF drivers in 85 consecutive patients. Focal and rotational activity (RAc) were identified with the CartoFinder system and activation sequence analysis. The impact of RAc location on post-ablation outcomes was analyzed.

Results

This study included 64 men and 21 women. RAc was detected in 73.4% of men and 38.1% of women (p = 0.003). RAc patients had higher left atrium (LA) voltage (0.64 ± 0.3 vs. 0.50 ± 0.2 mV; p = 0.01), RAc sites had higher voltage than non-RAc sites 0.77 ± 0.46 vs. 0.53 ± 0.37 mV (p < 0.001). Women had lower LA voltage than men (0.42 vs. 0.64 mV; p < 0.001), including pulmonary vein (PV) antra (0.16 vs. 0.30 mV; p < 0.001) and posterior wall (0.34 vs. 0.51 mV; p < 0.001). RAc in the posterior atrium was recorded in few women (23.8 vs. 54.7% in men; p = 0.014). AF recurrence rate was higher in patients with RAc outside WACPVI than those with all RAc inside WACPVI or no RAc (63.4 vs. 11.1 and 31.0%; p = 0.008 and p = 0.01). Comparison of selected patients using propensity score matching confirmed lower atrial voltage (0.4 ± 0.2 vs. 0.7 ± 0.3 mV; p = 0.007) and less RAc (38 vs. 75%; p = 0.02) in women.

Conclusion

Women have shown more advanced structural remodeling at ablation, which is associated with a lower incidence of RAc (usually located outside the WACPVI). These findings could explain post-ablation gender differences.

Simulation of atrial fibrillation in a non-ohmic propagation model with dynamic gap junctions

Gap junctions exhibit nonlinear electrical properties that have been hypothesized to be relevant to arrhythmogenicity in a structurally remodeled tissue. Large-scale implementation of gap junction dynamics in 3D propagation models remains challenging. We aim to quantify the impact of nonlinear diffusion during episodes of arrhythmias simulated in a left atrial model. Homogenization of conduction properties in the presence of nonlinear gap junctions was performed by generalizing a previously developed mathematical framework. A monodomain model was solved in which conductivities were time-varying and depended on transjunctional potentials. Gap junction conductances were derived from a simplified Vogel-Weingart model with first-order gating and adjustable time constant. A bilayer interconnected cable model of the left atrium with 100  μm resolution was used. The diffusion matrix was recomputed at each time step according to the state of the gap junctions. Sinus rhythm and atrial fibrillation episodes were simulated in remodeled tissue substrates. Slow conduction was induced by reduced coupling and by diffuse or stringy fibrosis. Simulations starting from the same initial conditions were repeated with linear and nonlinear gap junctions. The discrepancy in activation times between the linear and nonlinear diffusion models was quantified. The results largely validated the linear approximation for conduction velocities >20 cm/s. In very slow conduction substrates, the discrepancy accumulated over time during atrial fibrillation, eventually leading to qualitative differences in propagation patterns, while keeping the descriptive statistics, such as cycle lengths, unchanged. The discrepancy growth rate was increased by impaired conduction, fibrosis, conduction heterogeneity, lateral uncoupling, fast gap junction time constant, and steeper action potential duration restitution.

Anti-atrial Fibrillation Effects of Pulmonary Vein Isolation With or Without Ablation Gaps: A Computational Modeling Study

Background

Although pulmonary vein isolation (PVI) gaps contribute to recurrence after atrial fibrillation (AF) catheter ablation, the mechanism is unclear. We used realistic computational human AF modeling to explore the AF wave-dynamic changes of PVI with gaps (PVI-gaps).

Methods

We included 40 patients (80% male, 61.0 ± 9.8 years old, 92.5% persistent AF) who underwent AF catheter ablation to develop our realistic computational AF model. We compared the effects of a complete PVI (CPVI) and PVI-gap (2-mm × 4) on the AF wave-dynamics by evaluating the dominant frequency (DF), spatial change of DF, maximal slope of the action potential duration restitution curve (Smax), and AF defragmentation rate (termination or change to atrial tachycardia), and tested the effects of additional virtual interventions and flecainide on ongoing AF with PVI-gaps.

Results

Compared with the baseline AF, CPVIs significantly reduced extra-PV DFs (p < 0.001), but PVI-gaps did not. COV-DFs were greater after CPVIs than PVI-gaps (p < 0.001). Neither CPVIs nor PVI-gaps changed the mean Smax. CPVIs resulted in higher AF defragmentation rates (80%) than PVI-gaps (12.5%, p < 0.001). In ongoing AF after PVI-gaps, the AF defragmentation rates after a wave-breaking gap ablation, extra-PV DF ablation, or flecainide were 60.0, 34.3, and 25.7%, respectively (p = 0.010).

Conclusion

CPVIs effectively reduced the DF, increased its spatial heterogeneity in extra-PV areas, and offered better anti-AF effects than extra-PV DF ablation or additional flecainide in PVI-gap conditions.

Prednisone Ameliorates Atrial Inflammation and Fibrosis in Atrial Tachypacing Dogs

Atrial inflammation and fibrosis have long been considered culprits in the development of atrial fibrillation (AF). Prior clinical studies showed that corticosteroid therapy is beneficial in patients with AF. Here we sought to determine whether prednisone treatment prevents atrial tachypacing (ATP) induced atrial fibrosis.Dogs were randomized into the sham, ATP, ATP + low-dose prednisone (ALP), and ATP + high-dose prednisone (AHP) groups. After 6 days of recovery from surgery, dogs were subjected to ATP at 400 beats per minute for 4 weeks while being treated with prednisone (15 or 40 mg/day) or a placebo. Pacemakers were not activated in the sham group.Compared with the ATP group, the expression of collagen I, collagen III, α-smooth muscle actin, transforming growth factor-β1 and connective tissue growth factor were significantly reduced in the ALP and AHP groups. Fluorescence assays showed that reactive oxygen species formation in the right atrium was suppressed in the ALP and AHP groups compared with the ATP group. The protein level of NADPH oxidase 2 was reduced in the ALP and AHP groups' versus ATP group, while NOX4 and NOX5 were unchanged. ATP-induced downregulation of BH4 and eNOS uncoupling in the atria was partially restored in the prednisone-treated groups.Our study demonstrated that atrial fibrosis induced by ATP were suppressed by prednisone. Low-dose prednisone was also effective in suppressing the development of atrial fibrosis.

Restitution Slope Affects the Outcome of Dominant Frequency Ablation in Persistent Atrial Fibrillation: CUVIA-AF2 Post-Hoc Analysis Based on Computational Modeling Study

Introduction

Although the dominant frequency (DF) localizes the reentrant drivers and the maximal slope of the action potential duration (APD) restitution curve (Smax) reflects the tendency of the wave-break, their interaction has never been studied. We hypothesized that DF ablation has different effects on atrial fibrillation (AF) depending on Smax.

Methods

We studied the DF and Smax in 25 realistic human persistent AF model samples (68% male, 60 ± 10 years old). Virtual AF was induced by ramp pacing measuring Smax, followed by spatiotemporal DF evaluation for 34 s. We assessed the DF ablation effect depending on Smax in both computational modeling and a previous clinical trial, CUVIA-AF (170 patients with persistent AF, 70.6% male, 60 ± 11 years old).

Results

Mean DF had an inverse relationship with Smax regardless of AF acquisition timing (p &lt; 0.001). Virtual DF ablations increased the defragmentation rate compared to pulmonary vein isolation (PVI) alone (p = 0.015), especially at Smax &lt;1 (61.5 vs. 7.7%, p = 0.011). In post-DF ablation defragmentation episodes, DF was significantly higher (p = 0.002), and Smax was lower (p = 0.003) than in episodes without defragmentation. In the post-hoc analysis of CUVIA-AF2, we replicated the inverse relationship between Smax and DF (r = −0.47, p &lt; 0.001), and we observed better rhythm outcomes of clinical DF ablations in addition to a PVI than of empirical PVI at Smax &lt;1 [hazard ratio 0.45, 95% CI (0.22–0.89), p = 0.022; log-rank p = 0.021] but not at ≥ 1 (log-rank p = 0.177).

Conclusion

We found an inverse relationship between DF and Smax and the outcome of DF ablation after PVI was superior at the condition with Smax &lt;1 in both in-silico and clinical trials.