Reentry and Ectopic Pacemakers Emerge in a Three-Dimensional Model for a Slab of Cardiac Tissue with Diffuse Microfibrosis near the Percolation Threshold

A Matlab Toolbox for cardiac electrophysiology simulations on patient-specific geometries

In this paper, we present CardioMat, a Matlab toolbox for cardiac electrophysiology simulation based on patient-specific anatomies. The strength of CardioMat is the easy and fast construction of electrophysiology cardiac digital twins from segmented anatomical images in a general-purpose software such as Matlab. CardioMat implements a quasi-automatic pipeline that guides the user toward the construction of anatomically detailed cardiac electrophysiology models. Importantly, the CardioMat framework includes the generation of physiologically plausible fiber orientation and Purkinje networks. The main novelty of our framework is its ability to handle voxel-based geometries as produced by segmentation procedures directly, without the need for an unstructured mesh. Indeed, the CardioMat monodomain solver uses a smoothed boundary approach and runs completely on GPU for fast simulations. We employed CardioMat in different application scenarios to show its potentialities and provide preliminary assessment of the feasibility, diagnostic performance, and accuracy of the toolbox. In particular, we showed that CardioMat simulations derived from post-infarction patients hold high sensitivity, specificity, predictive value, and accuracy for localization of deceleration zones in sinus rhythm.

Perlin noise generation of physiologically realistic cardiac fibrosis

Fibrosis, a pathological increase in extracellular matrix proteins, is a significant health issue that hinders the function of many organs in the body, in some cases fatally. In the heart, fibrosis impacts on electrical propagation in a complex and poorly predictable fashion, potentially serving as a substrate for dangerous arrhythmias. Individual risk depends on the spatial manifestation of fibrotic tissue, and learning the spatial arrangement on the fine scale in order to predict these impacts still relies upon invasive ex vivo procedures. As a result, the effects of spatial variability on the symptomatic impact of cardiac fibrosis remain poorly understood. In this work, we address the issue of availability of such imaging data via a computational methodology for generating new realisations of cardiac fibrosis microstructure. Using the Perlin noise technique from computer graphics, together with an automated calibration process that requires only a single training image, we demonstrate successful capture of collagen texturing in four types of fibrosis microstructure observed in histological sections. We then use this generator to quantitatively analyse the conductive properties of these different types of cardiac fibrosis, as well as produce three-dimensional realisations of histologically-observed patterning. Owing to the generator's flexibility and automated calibration process, we also anticipate that it might be useful in producing additional realisations of other physiological structures.

'Trapped re-entry' as source of acute focal atrial arrhythmias

AIMS

Diseased atria are characterized by functional and structural heterogeneities, adding to abnormal impulse generation and propagation. These heterogeneities are thought to lie at the origin of fractionated electrograms recorded during sinus rhythm (SR) in atrial fibrillation (AF) patients and are assumed to be involved in the onset and perpetuation (e.g. by re-entry) of this disorder. The underlying mechanisms, however, remain incompletely understood. Here, we tested whether regions of dense fibrosis could create an electrically isolated conduction pathway (EICP) in which re-entry could be established via ectopy and local block to become 'trapped'. We also investigated whether this could generate local fractionated electrograms and whether the re-entrant wave could 'escape' and cause a global tachyarrhythmia due to dynamic changes at a connecting isthmus.

METHODS AND RESULTS

To precisely control and explore the geometrical properties of EICPs, we used light-gated depolarizing ion channels and patterned illumination for creating specific non-conducting regions in silico and in vitro. Insight from these studies was used for complementary investigations in virtual human atria with localized fibrosis. We demonstrated that a re-entrant tachyarrhythmia can exist locally within an EICP with SR prevailing in the surrounding tissue and identified conditions under which re-entry could escape from the EICP, thereby converting a local latent arrhythmic source into an active driver with global impact on the heart. In a realistic three-dimensional model of human atria, unipolar epicardial pseudo-electrograms showed fractionation at the site of 'trapped re-entry' in coexistence with regular SR electrograms elsewhere in the atria. Upon escape of the re-entrant wave, acute arrhythmia onset was observed.

CONCLUSIONS

Trapped re-entry as a latent source of arrhythmogenesis can explain the sudden onset of focal arrhythmias, which are able to transgress into AF. Our study might help to improve the effectiveness of ablation of aberrant cardiac electrical signals 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.

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.

Identifying locations susceptible to micro-anatomical reentry using a spatial network representation of atrial fibre maps

Micro-anatomical reentry has been identified as a potential driver of atrial fibrillation (AF). In this paper, we introduce a novel computational method which aims to identify which atrial regions are most susceptible to micro-reentry. The approach, which considers the structural basis for micro-reentry only, is based on the premise that the accumulation of electrically insulating interstitial fibrosis can be modelled by simulating percolation-like phenomena on spatial networks. Our results suggest that at high coupling, where micro-reentry is rare, the micro-reentrant substrate is highly clustered in areas where the atrial walls are thin and have convex wall morphology, likely facilitating localised treatment via ablation. However, as transverse connections between fibres are removed, mimicking the accumulation of interstitial fibrosis, the substrate becomes less spatially clustered, and the bias to forming in thin, convex regions of the atria is reduced, possibly restricting the efficacy of localised ablation. Comparing our algorithm on image-based models with and without atrial fibre structure, we find that strong longitudinal fibre coupling can suppress the micro-reentrant substrate, whereas regions with disordered fibre orientations have an enhanced risk of micro-reentry. With further development, these methods may be useful for modelling the temporal development of the fibrotic substrate on an individualised basis.

Critical Requirements for the Initiation of a Cardiac Arrhythmia in Rat Ventricle: How Many Myocytes?

Cardiovascular disease is the leading cause of death worldwide due in a large part to arrhythmia. In order to understand how calcium dynamics play a role in arrhythmogenesis, normal and dysfunctional Ca²⁺ signaling in a subcellular, cellular, and tissued level is examined using cardiac ventricular myocytes at a high temporal and spatial resolution using multiscale computational modeling. Ca²⁺ sparks underlie normal excitation-contraction coupling. However, under pathological conditions, Ca²⁺ sparks can combine to form Ca²⁺ waves. These propagating elevations of (Ca²⁺)_i can activate an inward Na⁺-Ca²⁺ exchanger current (I_NCX) that contributes to early after-depolarization (EADs) and delayed after-depolarizations (DADs). However, how cellular currents lead to full depolarization of the myocardium and how they initiate extra systoles is still not fully understood. This study explores how many myocytes must be entrained to initiate arrhythmogenic depolarizations in biophysically detailed computational models. The model presented here suggests that only a small number of myocytes must activate in order to trigger an arrhythmogenic propagating action potential. These conditions were examined in 1-D, 2-D, and 3-D considering heart geometry. The depolarization of only a few hundred ventricular myocytes is required to trigger an ectopic depolarization. The number decreases under disease conditions such as heart failure. Furthermore, in geometrically restricted parts of the heart such as the thin muscle strands found in the trabeculae and papillary muscle, the number of cells needed to trigger a propagating depolarization falls even further to less than ten myocytes.

Multiparametric analysis of geometric features of fibrotic textures leading to cardiac arrhythmias

One of the important questions in cardiac electrophysiology is to characterise the arrhythmogenic substrate; for example, from the texture of the cardiac fibrosis, which is considered one of the major arrhythmogenic conditions. In this paper, we perform an extensive in silico study of the relationships between various local geometric characteristics of fibrosis on the onset of cardiac arrhythmias. In order to define which texture characteristics have better predictive value, we induce arrhythmias by external stimulation, selecting 4363 textures in which arrhythmia can be induced and also selecting 4363 non-arrhythmogenic textures. For each texture, we determine such characteristics as cluster area, solidity, mean distance, local density and zig-zag propagation path, and compare them in arrhythmogenic and non-arrhythmogenic cases. Our study shows that geometrical characteristics, such as cluster area or solidity, turn out to be the most important for prediction of the arrhythmogenic textures. Overall, we were able to achieve an accuracy of 67% for the arrhythmogenic texture-classification problem. However, the accuracy of predictions depends on the size of the region chosen for the analysis. The optimal size for the local areas of the tissue was of the order of 0.28 of the wavelength of the arrhythmia. We discuss further developments and possible applications of this method for characterising the substrate of arrhythmias in fibrotic textures.

Minimal Functional Clusters Predict the Probability of Reentry in Cardiac Fibrotic Tissue

Cardiac fibrosis is a well-known arrhythmogenic condition which can lead to sudden cardiac death. Physically, fibrosis can be viewed as a large number of small obstacles in an excitable medium, which may create nonlinear wave turbulence or reentry. The relation between the specific texture of fibrosis and the onset of reentry is of great theoretical and practical importance. Here, we present a conceptual framework which combines functional aspects of propagation manifested as conduction blocks, with reentry wavelength and geometrical clusters of fibrosis. We formulate them into the single concept of minimal functional cluster and through extensive simulations show that it characterizes the path of reexcitation accurately, and importantly its size distribution quantitatively predicts the reentry probability for different fibrosis densities and tissue excitabilities.

Machine Learning Identification of Pro-arrhythmic Structures in Cardiac Fibrosis

Cardiac fibrosis and other scarring of the heart, arising from conditions ranging from myocardial infarction to ageing, promotes dangerous arrhythmias by blocking the healthy propagation of cardiac excitation. Owing to the complexity of the dynamics of electrical signalling in the heart, however, the connection between different arrangements of blockage and various arrhythmic consequences remains poorly understood. Where a mechanism defies traditional understanding, machine learning can be invaluable for enabling accurate prediction of quantities of interest (measures of arrhythmic risk) in terms of predictor variables (such as the arrangement or pattern of obstructive scarring). In this study, we simulate the propagation of the action potential (AP) in tissue affected by fibrotic changes and hence detect sites that initiate re-entrant activation patterns. By separately considering multiple different stimulus regimes, we directly observe and quantify the sensitivity of re-entry formation to activation sequence in the fibrotic region. Then, by extracting the fibrotic structures around locations that both do and do not initiate re-entries, we use neural networks to determine to what extent re-entry initiation is predictable, and over what spatial scale conduction heterogeneities appear to act to produce this effect. We find that structural information within about 0.5 mm of a given point is sufficient to predict structures that initiate re-entry with more than 90% accuracy.

Simulation of diffuse and stringy fibrosis in a bilayer interconnected cable model of the left atrium

AIMS

The aim of this study is to design a computer model of the left atrium for investigating fibre-orientation-dependent microstructure such as stringy fibrosis.

METHODS AND RESULTS

We developed an approach for automatic construction of bilayer interconnected cable models from left atrial geometry and epi- and endocardial fibre orientation. The model consisted of two layers (epi- and endocardium) of longitudinal and transverse cables intertwined-like fabric threads, with a spatial discretization of 100 µm. Model validation was performed by comparison with cubic volumetric models in normal conditions. Then, diffuse (n = 2904), stringy (n = 3600), and mixed fibrosis patterns (n = 6840) were randomly generated by uncoupling longitudinal and transverse connections in the interconnected cable model. Fibrosis density was varied from 0% to 40% and mean stringy obstacle length from 0.1 to 2 mm. Total activation time, apparent anisotropy ratio, and local activation time jitter were computed during normal rhythm in each pattern. Non-linear regression formulas were identified for expressing measured propagation parameters as a function of fibrosis density and obstacle length (stringy and mixed patterns). Longer obstacles (even below tissue space constant) were independently associated with prolonged activation times, increased anisotropy, and local fluctuations in activation times. This effect was increased by endo-epicardial dissociation and mitigated when fibrosis was limited to the epicardium.

CONCLUSION

Interconnected cable models enable the study of microstructure in organ-size models despite limitations in the description of transmural structures.

Computational modeling identifies embolic stroke of undetermined source patients with potential arrhythmic substrate

Cardiac magnetic resonance imaging (MRI) has revealed fibrosis in embolic stroke of undetermined source (ESUS) patients comparable to levels seen in atrial fibrillation (AFib). We used computational modeling to understand the absence of arrhythmia in ESUS despite the presence of putatively pro-arrhythmic fibrosis. MRI-based atrial models were reconstructed for 45 ESUS and 45 AFib patients. The fibrotic substrate's arrhythmogenic capacity in each patient was assessed computationally. Reentrant drivers were induced in 24/45 (53%) ESUS and 22/45 (49%) AFib models. Inducible models had more fibrosis (16.7 ± 5.45%) than non-inducible models (11.07 ± 3.61%; p<0.0001); however, inducible subsets of ESUS and AFib models had similar fibrosis levels (p=0.90), meaning that the intrinsic pro-arrhythmic substrate properties of fibrosis in ESUS and AFib are indistinguishable. This suggests that some ESUS patients have latent pre-clinical fibrotic substrate that could be a future source of arrhythmogenicity. Thus, our work prompts the hypothesis that ESUS patients with fibrotic atria are spared from AFib due to an absence of arrhythmia triggers.

3D Electrophysiological Modeling of Interstitial Fibrosis Networks and Their Role in Ventricular Arrhythmias in Non-Ischemic Cardiomyopathy

OBJECTIVE

Interstitial fibrosis is a pathological expansion of the heart's inter-cellular collagen matrix. It is a potential complication of nonischemic cardiomyopathy (NICM), a class of diseases involving electrical and or mechanical dysfunction of cardiac tissue not caused by atherosclerosis. Patients with NICM and interstitial fibrosis often suffer from life threatening arrhythmias, which we aim to simulate in this study.

METHODS

Our methodology builds on an efficient discrete finite element (DFE) method which allows for the representation of fibrosis as infinitesimal splits in a mesh. We update the DFE method with a local connectivity analysis which creates a consistent topology in the fibrosis network. This is particularly important in nonischemic disease due to the potential presence of large and contiguous fibrotic regions and therefore potentially complex fibrosis networks.

RESULTS

In experiments with an image-based model, we demonstrate that our methodology is able to simulate reentrant electrical events associated with cardiac arrhythmias. These reentries depended crucially upon sufficient fibrosis density, which was marked by conduction slowing at high pacing rates. We also created a 2D test-case which demonstrated that fibrosis topologies can modulate transient conduction block, and thereby reentrant activations.

CONCLUSION

Ventricular arrhythmias due to interstitial fibrosis in NICM can be efficiently simulated using our methods in medical image based geometries. Furthermore, fibrosis topology modulates transient conduction block, and should be accounted for in electrophysiological simulations with interstitial fibrosis.

SIGNIFICANCE

Our study provides methodology which has the potential to predict arrhythmias and to optimize treatments non-invasively for nonischemic cardiomyopathies.

Myocardial Fibrosis in a 3D Model: Effect of Texture on Wave Propagation

Non-linear electrical waves propagate through the heart and control cardiac contraction. Abnormal wave propagation causes various forms of the heart disease and can be lethal. One of the main causes of abnormality is a condition of cardiac fibrosis, which, from mathematical point of view, is the presence of multiple non-conducting obstacles for wave propagation. The fibrosis can have different texture which varies from diffuse (e.g., small randomly distributed obstacles), patchy (e.g., elongated interstitional stria), and focal (e.g., post-infarct scars) forms. Recently, Nezlobinsky et al. (2020) used 2D biophysical models to quantify the effects of elongation of obstacles (fibrosis texture) and showed that longitudinal and transversal propagation differently depends on the obstacle length resulting in anisotropy for wave propagation. In this paper, we extend these studies to 3D tissue models. We show that 3D consideration brings essential new effects; for the same obstacle length in 3D systems, anisotropy is about two times smaller compared to 2D, however, wave propagation is more stable with percolation threshold of about 60% (compared to 35% in 2D). The percolation threshold increases with the obstacle length for the longitudinal propagation, while it decreases for the transversal propagation. Further, in 3D, the dependency of velocity on the obstacle length for the transversal propagation disappears.

Variability in electrophysiological properties and conducting obstacles controls re-entry risk in heterogeneous ischaemic tissue

Ischaemia, in which inadequate blood supply compromises and eventually kills regions of cardiac tissue, can cause many types of arrhythmia, some life-threatening. A significant component of this is the effects of the resulting hypoxia, and concomitant hyperklaemia and acidosis, on the electrophysiological properties of myocytes. Clinical and experimental data have also shown that regions of structural heterogeneity (fibrosis, necrosis, fibro-fatty infiltration) can act as triggers for arrhythmias under acute ischaemic conditions. Mechanistic models have successfully captured these effects in silico. However, the relative significance of these separate facets of the condition, and how sensitive arrhythmic risk is to the extents of each, is far less explored. In this work, we use partitioned Gaussian process emulation and new metrics for source-sink mismatch that rely on simulations of bifurcating cardiac fibres to interrogate a model of heterogeneous ischaemic tissue. Re-entries were most sensitive to the level of hypoxia and the fraction of non-excitable tissue. In addition, our results reveal both protective and pro-arrhythmic effects of hyperklaemia, and present the levels of hyperklaemia, hypoxia and percentage of non-excitable tissue that pose the highest arrhythmic risks. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

An audit of uncertainty in multi-scale cardiac electrophysiology models

Models of electrical activation and recovery in cardiac cells and tissue have become valuable research tools, and are beginning to be used in safety-critical applications including guidance for clinical procedures and for drug safety assessment. As a consequence, there is an urgent need for a more detailed and quantitative understanding of the ways that uncertainty and variability influence model predictions. In this paper, we review the sources of uncertainty in these models at different spatial scales, discuss how uncertainties are communicated across scales, and begin to assess their relative importance. We conclude by highlighting important challenges that continue to face the cardiac modelling community, identifying open questions, and making recommendations for future studies. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Understanding the transition from paroxysmal to persistent atrial fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhytmia, characterized by the chaotic motion of electrical wavefronts in the atria. In clinical practice, AF is classified under two primary categories: paroxysmal AF, short intermittent episodes separated by periods of normal electrical activity; and persistent AF, longer uninterrupted episodes of chaotic electrical activity. However, the precise reasons why AF in a given patient is paroxysmal or persistent is poorly understood. Recently, we have introduced the percolation-based Christensen-Manani-Peters (CMP) model of AF which naturally exhibits both paroxysmal and persistent AF, but precisely how these differences emerge in the model is unclear. In this paper, we dissect the CMP model to identify the cause of these different AF classifications. Starting from a mean-field model where we describe AF as a simple birth-death process, we add layers of complexity to the model and show that persistent AF arises from reentrant circuits which exhibit an asymmetry in their probability of activation relative to deactivation. As a result, different simulations generated at identical model parameters can exhibit fibrillatory episodes spanning several orders of magnitude from a few seconds to months. These findings demonstrate that diverse, complex fibrillatory dynamics can emerge from very simple dynamics in models of AF.

Anisotropic conduction in the myocardium due to fibrosis: the effect of texture on wave propagation

Cardiac fibrosis occurs in many forms of heart disease. It is well established that the spatial pattern of fibrosis, its texture, substantially affects the onset of arrhythmia. However, in most modelling studies fibrosis is represented by multiple randomly distributed short obstacles that mimic only one possible texture, diffuse fibrosis. An important characteristic feature of other fibrosis textures, such as interstitial and patchy textures, is that fibrotic inclusions have substantial length, which is suggested to have a pronounced effect on wave propagation. In this paper, we study the effect of the elongation of inexcitable inclusions (obstacles) on wave propagation in a 2D model of cardiac tissue described by the TP06 model for human ventricular cells. We study in detail how the elongation of obstacles affects various characteristics of the waves. We quantify the anisotropy induced by the textures, its dependency on the obstacle length and the effects of the texture on the shape of the propagating wave. Because such anisotropy is a result of zig-zag propagation we show, for the first time, quantification of the effects of geometry and source-sink relationship, on the zig-zag nature of the pathway of electrical conduction. We also study the effect of fibrosis in the case of pre-existing anisotropy and introduce a procedure for scaling of the fibrosis texture. We show that fibrosis can decrease or increase the preexisting anisotropy depending on its scaled texture.

Unified mechanism of local drivers in a percolation model of atrial fibrillation

The mechanisms of atrial fibrillation (AF) are poorly understood, resulting in disappointing success rates of ablative treatment. Different mechanisms defined largely by different atrial activation patterns have been proposed and, arguably, this dispute has slowed the progress of AF research. Recent clinical evidence suggests a unifying mechanism of local drivers based on sustained reentrant circuits in the complex atrial architecture. Here, we present a percolation inspired computational model showing spontaneous emergence of AF that strongly supports, and gives a theoretical explanation for, the clinically observed diversity of activation. We show that the difference in surface activation patterns is a direct consequence of the thickness of the discrete network of heart muscle cells through which electrical signals percolate to reach the imaged surface. The model naturally follows the clinical spectrum of AF spanning sinus rhythm, paroxysmal AF, and persistent AF as the decoupling of myocardial cells results in the lattice approaching the percolation threshold. This allows the model to make the prediction that, for paroxysmal AF, reentrant circuits emerge near the endocardium, but in persistent AF they emerge deeper in the bulk of the atrial wall. If experimentally verified, this may go towards explaining the lowering ablation success rate as AF becomes more persistent.

Scar shape analysis and simulated electrical instabilities in a non-ischemic dilated cardiomyopathy patient cohort

This paper presents a morphological analysis of fibrotic scarring in non-ischemic dilated cardiomyopathy, and its relationship to electrical instabilities which underlie reentrant arrhythmias. Two dimensional electrophysiological simulation models were constructed from a set of 699 late gadolinium enhanced cardiac magnetic resonance images originating from 157 patients. Areas of late gadolinium enhancement (LGE) in each image were assigned one of 10 possible microstructures, which modelled the details of fibrotic scarring an order of magnitude below the MRI scan resolution. A simulated programmed electrical stimulation protocol tested each model for the possibility of generating either a transmural block or a transmural reentry. The outcomes of the simulations were compared against morphological LGE features extracted from the images. Models which blocked or reentered, grouped by microstructure, were significantly different from one another in myocardial-LGE interface length, number of components and entropy, but not in relative area and transmurality. With an unknown microstructure, transmurality alone was the best predictor of block, whereas a combination of interface length, transmurality and number of components was the best predictor of reentry in linear discriminant analysis.

Addressing challenges of quantitative methodologies and event interpretation in the study of atrial fibrillation

Atrial fibrillation (AF) is the commonest arrhythmia, yet the mechanisms of its onset and persistence are incompletely known. Although techniques for quantitative assessment have been investigated, there have been few attempts to integrate this information to advance disease treatment protocols. In this review, key quantitative methods for AF analysis are described, and suggestions are provided for the coordination of the available information, and to develop foci and directions for future research efforts. Quantitative biologists may have an interest in this topic in order to develop machine learning and tools for arrhythmia characterization, but they may perhaps have a minimal background in the clinical methodology and in the types of observed events and mechanistic hypotheses that have thus far been developed. We attempt to address these issues via exploration of the published literature. Although no new data is presented in this review, examples are shown of current lines of investigation, and in particular, how electrogram analysis and whole-chamber quantitative modeling of the left atrium may be useful to characterize fibrillatory patterns of activity, so as to propose avenues for more efficacious acquisition and interpretation of AF data.

Self-organization of conducting pathways explains electrical wave propagation in cardiac tissues with high fraction of non-conducting cells

Cardiac fibrosis occurs in many forms of heart disease and is considered to be one of the main arrhythmogenic factors. Regions with a high density of fibroblasts are likely to cause blocks of wave propagation that give rise to dangerous cardiac arrhythmias. Therefore, studies of the wave propagation through these regions are very important, yet the precise mechanisms leading to arrhythmia formation in fibrotic cardiac tissue remain poorly understood. Particularly, it is not clear how wave propagation is organized at the cellular level, as experiments show that the regions with a high percentage of fibroblasts (65-75%) are still conducting electrical signals, whereas geometric analysis of randomly distributed conducting and non-conducting cells predicts connectivity loss at 40% at the most (percolation threshold). To address this question, we used a joint in vitro-in silico approach, which combined experiments in neonatal rat cardiac monolayers with morphological and electrophysiological computer simulations. We have shown that the main reason for sustainable wave propagation in highly fibrotic samples is the formation of a branching network of cardiomyocytes. We have successfully reproduced the morphology of conductive pathways in computer modelling, assuming that cardiomyocytes align their cytoskeletons to fuse into cardiac syncytium. The electrophysiological properties of the monolayers, such as conduction velocity, conduction blocks and wave fractionation, were reproduced as well. In a virtual cardiac tissue, we have also examined the wave propagation at the subcellular level, detected wavebreaks formation and its relation to the structure of fibrosis and, thus, analysed the processes leading to the onset of arrhythmias.

Cardiac arrhythmias are one of the major causes of death in the industrialized world. The most dangerous ones are often caused by the blocks of propagation of electrical signals. One of the common factors that contribute to the likelihood of these blocks, is a condition called cardiac fibrosis. In fibrosis, excitable cardiac tissue is partially replaced with the inexcitable and non-conducting connective tissue. The precise mechanisms leading to arrhythmia formation in fibrotic cardiac tissue remain poorly understood. Therefore, it is important to study wave propagation in fibrosis from cellular to tissue level. In this paper, we study tissues with high densities of non-conducting cells in experiments and computer simulations. We have observed a paradoxical ability of the tissue with extremely high portion of non-conducting cells (up to 75%) to conduct electrical signals and contract synchronously, whereas geometric analysis of randomly distributed cells predicted connectivity loss at 40% at the most. To explain this phenomenon, we have studied the patterns that cardiac cells form in the tissue and reproduced their self-organisation in a computer model. Our virtual model also took into account the polygonal shapes of the spreading cells and explained high arrhythmogenicity of fibrotic tissue.

The effect of complex intramural microstructure caused by structural remodeling on the stability of atrial fibrillation: Insights from a three-dimensional multi-layer modeling study

BACKGROUND

Recent researches have suggested that the complex three-dimensional structures caused by structural remodeling play a key role in atrial fibrillation (AF) substrates. Here we aimed to investigate this hypothesis using a multi-layer model representing intramural microstructural features.

METHODS

The proposed multi-layer model was composed of the endocardium, connection wall, and epicardium. In the connection wall, intramural fibrosis was simulated using fibrotic patches randomly scattered in the myocardial tissue of fibrotic layers, while endo-epicardial dissociation was simulated using myocardial patches randomly scattered in the fibrotic tissue of isolation layers. Multiple simulation groups were generated to quantitatively analyze the effects of endo-epicardial dissociation and intramural fibrosis on AF stability, including a stochastic group, interrelated groups, fibrosis-degree-controlled groups, and dissociation-degree-controlled groups.

RESULTS

1. Stable intramural re-entries were observed to move along complete re-entrant circuits inside the transmural wall in four of 65 simulations in the stochastic group. 2. About 21 of 23 stable simulations in the stochastic group were distributed in the areas with high endo-epicardial dissociation and intramural fibrosis. 3. The difference between fibrosis-degree-controlled groups and dissociation-degree-controlled groups suggested that some distributions of connection areas may affect AF episodes despite low intramural fibrosis and endo-epicardial dissociation. 4. The overview of tracking phase singularities revealed that endo-epicardial dissociation played a visible role in AF substrates.

CONCLUSION

The complex intramural microstructure is positively correlated with critical components of AF maintenance mechanisms. The occurrence of intramural re-entry further indicates the complexity of AF wave-dynamics.

Ectopic beats arise from micro-reentries near infarct regions in simulations of a patient-specific heart model

Ectopic beats are known to be involved in the initiation of a variety of cardiac arrhythmias. Although their location may vary, ectopic excitations have been found to originate from infarct areas, regions of micro-fibrosis and other heterogeneous tissues. However, the underlying mechanisms that link ectopic foci to heterogeneous tissues have yet to be fully understood. In this work, we investigate the mechanism of micro-reentry that leads to the generation of ectopic beats near infarct areas using a patient-specific heart model. The patient-specific geometrical model of the heart, including scar and peri-infarct zones, is obtained through magnetic resonance imaging (MRI). The infarct region is composed of ischemic myocytes and non-conducting cells (fibrosis, for instance). Electrophysiology is captured using an established cardiac myocyte model of the human ventricle modified to describe ischemia. The simulation results clearly reveal that ectopic beats emerge from micro-reentries that are sustained by the heterogeneous structure of the infarct regions. Because microscopic information about the heterogeneous structure of the infarct regions is not available, Monte-Carlo simulations are used to identify the probabilities of an infarct region to behave as an ectopic focus for different levels of ischemia and different percentages of non-conducting cells. From the proposed model, it is observed that ectopic beats are generated when a percentage of non-conducting cells is near a topological metric known as the percolation threshold. Although the mechanism for micro-reentries was proposed half a century ago to be a source of ectopic beats or premature ventricular contractions during myocardial infarction, the present study is the first to reproduce this mechanism in-silico using patient-specific data.

Atrial Fibrillation Mechanisms and Implications for Catheter Ablation

AF is a heterogeneous rhythm disorder that is related to a wide spectrum of etiologies and has broad clinical presentations. Mechanisms underlying AF are complex and remain incompletely understood despite extensive research. They associate interactions between triggers, substrate and modulators including ionic and anatomic remodeling, genetic predisposition and neuro-humoral contributors. The pulmonary veins play a key role in the pathogenesis of AF and their isolation is associated to high rates of AF freedom in patients with paroxysmal AF. However, ablation of persistent AF remains less effective, mainly limited by the difficulty to identify the sources sustaining AF. Many theories were advanced to explain the perpetuation of this form of AF, ranging from a single localized focal and reentrant source to diffuse bi-atrial multiple wavelets. Translating these mechanisms to the clinical practice remains challenging and limited by the spatio-temporal resolution of the mapping techniques. AF is driven by focal or reentrant activities that are initially clustered in a relatively limited atrial surface then disseminate everywhere in both atria. Evidence for structural remodeling, mainly represented by atrial fibrosis suggests that reentrant activities using anatomical substrate are the key mechanism sustaining AF. These reentries can be endocardial, epicardial, and intramural which makes them less accessible for mapping and for ablation. Subsequently, early interventions before irreversible remodeling are of major importance. Circumferential pulmonary vein isolation remains the cornerstone of the treatment of AF, regardless of the AF form and of the AF duration. No ablation strategy consistently demonstrated superiority to pulmonary vein isolation in preventing long term recurrences of atrial arrhythmias. Further research that allows accurate identification of the mechanisms underlying AF and efficient ablation should improve the results of PsAF ablation.

Killing Many Birds With Two Stones: Hypoxia and Fibrosis Can Generate Ectopic Beats in a Human Ventricular Model

During cardiac diseases many types of anatomical and functional remodeling of cardiac tissue can occur. In this work, we focus on two conditions: hypoxia and fibrosis, which are part of complex pathological modifications that take place in many cardiac diseases (hypertrophic cardiomyopathy, hypertensive heart disease, and recurrent myocardial infarction) and respiratory diseases (obstructive pulmonary disease, obstructive sleep apnea, and cystic fibrosis). Using computational models of cardiac electrophysiology, we evaluate if the interplay between hypoxia and fibrosis is sufficient to trigger cardiac arrhythmia. We study the mechanisms behind the generation of ectopic beats, an arrhythmic trigger also known as premature ventricular contractions (PVCs), in regions with high hypoxia and fibrosis. First, we modify an electrophysiological model of myocytes of the human left ventricle to include the effects of hypoxia. Second, diffuse fibrosis is modeled by randomly replacing cardiac myocytes by non-excitable and non-conducting cells. The Monte Carlo method is used to evaluate the probability of a region to generate ectopic beats with respect to different levels of hypoxia and fibrosis. In addition, we evaluate the minimum size of three-dimensional slabs needed to sustain reentries for different stimulation protocols. The observed mechanism behind the initiation of ectopic beats is unidirectional block, giving rise to sustained micro-reentries inside the region with diffuse fibrosis and hypoxia. In summary, our results suggest that hypoxia and fibrosis are sufficient for the creation of a focal region in the heart that generates PVCs.

Discretization-dependent model for weakly connected excitable media

Pattern formation has been widely observed in extended chemical and biological processes. Although the biochemical systems are highly heterogeneous, homogenized continuum approaches formed by partial differential equations have been employed frequently. Such approaches are usually justified by the difference of scales between the heterogeneities and the characteristic spatial size of the patterns. Under different conditions, for example, under weak coupling, discrete models are more adequate. However, discrete models may be less manageable, for instance, in terms of numerical implementation and mesh generation, than the associated continuum models. Here we study a model to approach discreteness which permits the computer implementation on general unstructured meshes. The model is cast as a partial differential equation but with a parameter that depends not only on heterogeneities sizes, as in the case of quasicontinuum models, but also on the discretization mesh. Therefore, we refer to it as a discretization-dependent model. We validate the approach in a generic excitable media that simulates three different phenomena: the propagation of action membrane potential in cardiac tissue, in myelinated axons of neurons, and concentration waves in chemical microemulsions.

Reentry via high-frequency pacing in a mathematical model for human-ventricular cardiac tissue with a localized fibrotic region

Localized heterogeneities, caused by the regional proliferation of fibroblasts, occur in mammalian hearts because of diseases like myocardial infarction. Such fibroblast clumps can become sources of pathological reentrant activities, e.g., spiral or scroll waves of electrical activation in cardiac tissue. The occurrence of reentry in cardiac tissue with heterogeneities, such as fibroblast clumps, can depend on the frequency at which the medium is paced. Therefore, it is important to study the reentry-initiating potential of such fibroblast clumps at different frequencies of pacing. We investigate the arrhythmogenic effects of fibroblast clumps at high- and low-frequency pacing. We find that reentrant waves are induced in the medium more prominently at high-frequency pacing than with low-frequency pacing. We also study the other factors that affect the potential of fibroblast clumps to induce reentry in cardiac tissue. In particular, we show that the ability of a fibroblast clump to induce reentry depends on the size of the clump, the distribution and percentage of fibroblasts in the clump, and the excitability of the medium. We study the process of reentry in two-dimensional and a three-dimensional mathematical models for cardiac tissue.

Ablation as targeted perturbation to rewire communication network of persistent atrial fibrillation

Persistent atrial fibrillation (AF) can be viewed as disintegrated patterns of information transmission by action potential across the communication network consisting of nodes linked by functional connectivity. To test the hypothesis that ablation of persistent AF is associated with improvement in both local and global connectivity within the communication networks, we analyzed multi-electrode basket catheter electrograms of 22 consecutive patients (63.5 ± 9.7 years, 78% male) during persistent AF before and after the focal impulse and rotor modulation-guided ablation. Eight patients (36%) developed recurrence within 6 months after ablation. We defined communication networks of AF by nodes (cardiac tissue adjacent to each electrode) and edges (mutual information between pairs of nodes). To evaluate patient-specific parameters of communication, thresholds of mutual information were applied to preserve 10% to 30% of the strongest edges. There was no significant difference in network parameters between both atria at baseline. Ablation effectively rewired the communication network of persistent AF to improve the overall connectivity. In addition, successful ablation improved local connectivity by increasing the average clustering coefficient, and also improved global connectivity by decreasing the characteristic path length. As a result, successful ablation improved the efficiency and robustness of the communication network by increasing the small-world index. These changes were not observed in patients with AF recurrence. Furthermore, a significant increase in the small-world index after ablation was associated with synchronization of the rhythm by acute AF termination. In conclusion, successful ablation rewires communication networks during persistent AF, making it more robust, efficient, and easier to synchronize. Quantitative analysis of communication networks provides not only a mechanistic insight that AF may be sustained by spatially localized sources and global connectivity, but also patient-specific metrics that could serve as a valid endpoint for therapeutic interventions.