Corrigendum: A review of healthy and fibrotic myocardium microstructure modeling and corresponding intracardiac electrograms
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
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.
1 Introduction
Patients with cardiac arrhythmias are often treated with ablation therapy. Substrate-based ablation therapy is guided by intracardiac measurements acquired from catheters inserted into the cardiac chamber that record the extracellular potential.
The signal recorded by an electrode with respect to a distant reference is called unipolar electrogram (uEGM). EGMs of several electrodes on a catheter and/or multiple catheter locations are used to understand the dynamics of the cardiac arrhythmia. However, the recorded uEGMs are affected by different artifacts such as contraction of the heart, breathing of the patient, far-field signals from distant parts of the heart and noise from different hardware components. To alleviate these issues, bipolar electrograms (biEGM) are most frequently used, which subtract the uEGMs of two close-by electrodes. In this way, artifacts that affect both electrodes in the same way are cancelled. The difference between two potentials is called voltage and we should keep in mind that we can only measure voltages. Therefore, uEGMs always have to be considered with respect to their (distant) reference electrode. In clinical literature, also the peak-to-peak amplitude of an electrogram signal (i.e., a voltage time course) is often called “voltage”.
The mathematical model of an excitable cell proposed by Hodgkin and Huxley, (1952), the tissue homogenization approach proposed by Schmitt, (1969), and the set of bidomain equations first applied by Tung, (1978) in 1978 is the most complete and accurate model that describes the spread of electrical depolarization across the myocardium and its cells.
Computational simulations based on this mathematical model have been used to understand the phenomena of the depolarization spread in cardiac tissue and their effects on electrogram genesis and morphology (Bishop and Plank, 2011b; Oesterlein et al., 2016; Roney et al., 2016; Pollnow et al., 2017; Beheshti et al., 2018; Hwang et al., 2019). While EGMs can be extracted from the extracellular medium in a bidomain simulation, this approach is computationally expensive. Thus, different methods based on excitation propagation simulations in the monodomain model have been proposed. Another modeling approach to accelerate the computation is the so-called reaction-eikonal model Neic et al. (2017), which can simulate physiological propagation using a coarser mesh (element average length 400 μm). In the monodomain model and the reaction-eikonal model, the extracellular potential is not calculated directly. However, it can be approximated with the pseudo-bidomain approach or the infinite homogeneous volume conductor method to obtain EGMs as detailed below. The infinite homogeneous volume conductor method approximates the extracellular potential caused by a group of cells spatially distributed in space and acting as sources of the electric field (Malmivuo and Plonsey, 1995).
In this review, we give an overview of the biophysical phenomena governing wave propagation in cardiac tissue and the corresponding extracellular potentials measured as electrograms. We will particularly focus on different approaches used to model fibrotic remodelling and simulate the corresponding electrograms to reproduce and understand the clinically observed changes in electrogram amplitude and morphology.
2 Intracardiac Electrograms
The electrical activity in the myocardium originates from the coordinated opening and closing of the ion channels in the cell membrane. The time course of the difference between the potential in the intracellular and in the extracellular medium is known as the action potential. In cardiac tissue, the cells are interconnected through gap junctions that will start a cascade effect of cellular activation along the major axis in which myocytes are aligned locally (also known as fiber direction), resulting in excitation propagation across the myocardium.
The extracellular field is a consequence of the spatial distribution of the transmembrane voltage of the cells in the myocardium (Figure 1A). An advancing depolarization wave in the cardiac tissue changes the spatial distribution of the extracellular potential. The extracellular potential can be measured as the uEGM at one electrode (technically the voltage between the extracellular potential at the measuring electrode with reference to for example, Wilson central terminal). The unipolar electrogram morphology is characterized by a biphasic symmetric shape (Figure 1B) where the positive phase (R-peak) indicates the approaching of the wavefront to the measuring electrode and the fast downslope indicates the moment that the wavefront is underneath the electrode. The opposing negative phase (S-peak) indicates the movement of the wavefront away from the measuring electrode. The peak-to-peak amplitude of the signal is also called “voltage” in the clinical literature. Peak-to-peak voltage is used as a marker to distinguish healthy from pathological tissue both for biEGMs (Jadidi et al., 2020) and uEGMs (Nairn et al., 2020b). However, biEGM amplitude can be affected by to several factors (Hwang et al., 2019) such as the orientation of the catheter (Schuler et al., 2013; Gaeta et al., 2020), the electrode spacing and size (Beheshti et al., 2018; Abdi et al., 2020; Nairn et al., 2020a; Takigawa et al., 2022), depolarization patterns (Jacquemet et al., 2003), substrate remodeling (Jacquemet et al., 2003; McDowell et al., 2012; Campos et al., 2013; Mendonca Costa et al., 2014; Roney et al., 2016; Sánchez et al., 2021b) and signal filter settings (Starreveld et al., 2020).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/908069/fphys-13-908069-HTML-r1/image_m/fphys-13-908069-g001.jpg)
Figure 1. Electrical propagation in healthy cardiac tissue. (A) Extracellular field caused by the depolarization of cardiomyocytes when an excitation propagates from left to right (green arrows). Spatial transmembrane voltage distribution is color-coded on the membranes and shown in the bottom row. The leftmost cells are already depolarized an in the action potential plateau while those on the right are still at resting membrane voltage. (B) Symmetric unipolar electrogram measured at the surface of the cardiac tissue (pink). The initial positive wave (R-peak) is caused by the wavefront approaching the electrode (dark gray), the polarity changes when the wavefront passes underneath the electrode, and the S-peak is caused by the wavefront traveling away from the measuring electrode.
3 Modelling Intracardiac Signals
Computational cardiac modeling has advanced rapidly in the last years and different numerical methods to simulate the propagation of the cardiac depolarization have been proposed over the years. Finite difference approaches have been widely used (Potse et al., 2006) and can be generalized for grids with distinct spacing (Trew et al., 2005; Sánchez et al., 2019a). Also the finite element method has been used to discretize complex geometries such as the cardiac chambers to simulate cardiac electrophysiology (Vigmond et al., 2003; Cooper et al., 2015; Neic et al., 2017; Plank et al., 2021).
The bidomain model represents cardiac tissue as a homogenized medium composed of the intracellular and the extracellular domains. The two computational domains coexist in the bidomain model and occupy the same geometrical space:
where ϕ represents the electrical potential, the indices i and e refer to the intracellular and extracellular spaces, respectively. σ is the conductivity tensor, β is the surface to volume ratio of the myocytes and Iion the total transmembrane ionic current density defined by the cellular model. The latter is dependent on Vm and a vector ν of further state variables. Iintra (a transmembrane current density) and Iextra (an extracellular current density) describe external stimuli. If a bath surrounds the tissue, it is treated as an extension of the extracellular space.
Adding Eqs 1, 2 and incorporating it into Eq. 4 yields:
As mentioned before, the reference potential during an electro-anatomical mapping procedure is usually a potential in a remote site or an average of potential values such as Wilson’s central terminal. For a bidomain model, when calculating uEGMs, the reference potential can, for example, be considered as an average of the extracellular potential of the furthest surface with respect to the tissue (Colli Franzone et al., 2007; Keller et al., 2014), which is not a perfect approximation of a remote reference electrode (e.g., a surface patch on the back of the patient) but markedly reduces drift of the reference potential. The further away the reference is from the myocardial tissue in the model, the better the representation of the reference potential but also the higher the computational cost due to the extended computational domain. Considering the average potential in a remote surface or volume is numerically advantageous compared to defining a fixed reference potential as a Dirichlet boundary condition.
The monodomain model is an approximation that assumes that the anisotropy of the extracellular and intracellular conductivity are aligned. Therefore, under the assumption of equal anisotropy ratios, one needs to solve only the parabolic partial differential equation above with the monodomain conductivity set appropriately:
where the bidomain equivalent monodomain conductivity σm is given as
Potse et al. (2006) performed a thorough comparison between the results of the bidomain model and monodomain model. The authors conclude that the monodomain model, although being a simplification of the bidomain model, is sufficient to study and understand the electrical propagation in the cardiac tissue under physiological conditions as well as for electrically remodeled tissue (ionic current abnormalities). The acceleration of the wavefront at the tissue-to-blood interface due to the bath-loading effect can be represented with the augmented bidomain approach (Bishop and Plank, 2011b). One of the biggest disadvantages of the bidomain model is the long computation time that it requires. Therefore, a common modeling approach is to combine the monodomain model with independent forward calculation of extracellular potentials. The most simplistic approach is the infinite volume conductor assumption, which assumes that the cardiac tissue is immersed in a homogeneous extracellular medium with infinite extent. This approach was for example, used to study the relation of the spread of depolarization in the cardiac tissue to the genesis and morphology of the unipolar electrogram (Gima and Rudy, 2002; Ganesan et al., 2013; Ugarte et al., 2014; Cabrera-Lozoya et al., 2017; Hwang et al., 2019) but neglects the influence of the heterogeneous surrounding tissue like other cardiac chambers, the lungs or the liver.
Briefly, the source and the measuring point (electrode) for a dipole are assumed to be immersed in an unbounded (infinite) volume conductor with homogeneous properties. The time course of the potential of the dipole corresponds to the uEGM electrogram measured at a location x in a certain distance to the source located in the cardiac tissue (xsrc) with respect to a reference electrode in infinite distance using the integral solution to Poisson’s equation:
where ϕe is the extracellular potential, σ is the conductivity of the volume conductor, Isrc is the source current density and ‖x − xsrc‖ is the Euclidean distance from the source point to the measuring point.
Bishop and Plank (2011a) proposed a combined bidomain and monodomain model (pseudo-bidomain) to calculate the extracellular potential. The proposed pseudo-bidomain approach computes the elliptic bidomain equation for a given transmembrane voltage distribution only at the time instants for which the extracellular potential is sampled. This approach is suitable to reproduce extracellular signals [EGM (Keller et al., 2012) and ECG (Nagel et al., 2022)] for a finite surrounding conductive medium (bath, potentially inhomogeneous) and is computationally efficient.
3.1 Factors Affecting the Intracardiac Signals
Using the bidomain model and realistic geometries of commercially available catheters can help to better understand EGM morphology (Schuler et al., 2013; Pollnow et al., 2017; Sánchez et al., 2021a). Schuler et al. (2013) modeled a realistic 7F catheter with two electrodes such that the tip was at the center of the tissue and in direct contact with the tissue patch surface. The catheter angle was changed with respect to the surface of the tissue (elevation) and to the wavefront propagation direction (rotation). Additionally, the authors explored the impact of the tissue thickness and conduction velocity on biEGM amplitude and duration. One of their main findings was that catheter orientation greatly affects the height and ratio of the positive and negative bipolar signal amplitude, which can be traced back to changes in the proximal signal. Moreover, the authors pointed out that the substrate characteristics (thickness and conduction velocity) mainly affect the biEGM peak-to-peak amplitude.
In new highly detailed bidomain simulations for this review, we show the biophysical phenomena of the spread of depolarization in the left atrium and the EGMs from a 7F LASSO™ (Biosense Webster) catheter in a healthy left atrium. Figure 2 shows that local activation time is the main factor that impacts the biEGM amplitude and that it is less sensitive to the wavefront direction. Additionally, bidomain simulations showed that biEGMs from electrodes that are not in direct contact with the tissue have the same activation time resulting in a small biEGM amplitude, which confirms the results previously shown (Gaeta et al., 2020). In brief, the atrial anatomical model (Roney et al., 2021) has a realistic wall thickness and an average edge length of 100 μm. Tissue conductivity was tune to achieved a conduction velocity of 40 cm/s (McDowell et al., 2013). The value of conductivity of the blood were as reported by Clerc (1976), the electrode conductivity was set to 1 × 1012 S/m to represent a good conductor that yields an isopotential volume, the conductivity of the catheter insulator was set close to zero (1 × 106 S/m).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/908069/fphys-13-908069-HTML-r1/image_m/fphys-13-908069-g002.jpg)
Figure 2. Bidomain simulation of a realistically deformed LASSO™ (Biosense Webster) catheter in a left atrium to study the genesis of different EGM morphologies in healthy myocardium. (A) The wavefront approaches the electrode pair 7–8 and activates both electrodes at the same time, the resulting bipolar electrogram with a reduced peak-to-peak amplitude (0.42 mV). (B) Several wavefronts approaching electrode pair 13–14, both unipolar electrograms are asymmetrical, lacking an S-peak; the resulting bipolar electrogram has a high peak-to-peak amplitude and a positive polarity. (C) The wavefront travels almost perpendicular to electrode pair 3–4; the electrodes are activated at different times, the resulting bipolar electrogram has negative polarity and a high peak-to-peak amplitude (7.45 mV).
The amplitude of uEGMs is affected by the geometrical properties of the electrode, such as the size of the electrode. Nairn et al. (2020a) performed a series of in silico experiments to understand the effect of the electrode size on the amplitude of the measured EGM. uEGM amplitude was shown to be inversely related to the size of the electrode. biEGM amplitude is additionally affected by the electrode pair spacing. Beheshti et al. (2018) showed that biEGM amplitude was increased when the electrode spacing increased. Assuming a plane wave and a perfectly symmetric uEGM in a simple thought experiment, the biEGM amplitude is zero for electrodes that are activated at exactly the same time. When increasing the distance between the electrodes, the peak-to-peak biEGM amplitude increases up to two times the uEGM amplitude. When further increasing the interelectrode distance, the biEGM amplitude decreases again until there is no more temporal overlap between the two uEGMs and the biEGM amplitude plateaus at the uEGM amplitude.
An additional factor that impacts the EGM amplitude and morphology are the filter settings (Schneider et al., 2004; Lin et al., 2007; Starreveld et al., 2020). In clinical practice, a bandpass filter is commonly used. However, the cut-off values of the bandpass filter differs for different mapping systems, catheters or due to the noise environment present in the specific electrophysiology laboratory. During an electroanatomical mapping procedure, uEGMs are typically filtered with a highpass of 0.5–2 Hz and a lowpass filter of 300–600 Hz biEGMs are typically bandpass filtered with a highpass of 1–30 Hz and a lowpass of 300–500 Hz. Both EGM types are also filtered at the frequency of the powerline with a notch filter (50 or 60 Hz). Figure 3 depicts the effect of the filter settings on both uEGMs (panel A) and biEGMs (panel C). In particular for biEGMs, the highpass filter cut-off value affects the measure voltage (Figure 3D). The higher amplitude of these simulated EGMs compared to clinical EGMs is likely due to the chosen extracellular conductivity, perfect contact of the electrode with the tissue and absence of losses along the signal chain.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/908069/fphys-13-908069-HTML-r1/image_m/fphys-13-908069-g003.jpg)
Figure 3. Filter effects on electrograms measured in healthy myocardium. (A) Unipolar electrograms with different highpass filter values. (B) Effect of different highpass cut-off values on the unipolar electrogram amplitude. (C) Bipolar electrogram with different highpass cut-off values. (D) Effect of different highpass cut-off values on the bipolar electrogram amplitude.
Considering the numerous factors that affect the uEGM and biEGM amplitude and morphology, standardized mapping modality (uEGM or biEGM), electrode size, electrode spacing and filter settings could increase comparability between studies. For modeling the healthy myocardium and electrograms, bidomain models provide the most accurate representation of the biophysical phenomena of depolarization and the influence of the catheter inside the cardiac cavity. Monodomain models and reaction-eikonal models in combination with forward calculation approaches to obtain the EGMs provide sufficient information about the propagation in the cardiac tissue in most scenarios. After reviewing the factors that influence the EGMs in the healthy myocardium, the next section covers factors that increase the complexity of the signals due to heterogeneities of the tissue and different patterns of propagation.
4 Myocardial Structural Remodeling and Intracardiac Signals
Structural remodeling alters the cardiac substrate, and the depolarization wavefront often has to follow a zig-zag pattern (Figure 4 white arrows). The zig-zag pattern of the propagation is reflected in uEGM and biEGM as fractionation in the signal due to constantly changing orientation of the wavefront. Fractionation is defined as an increase of deflections, thus an increase in complexity of the signal as well as a prolongation of the EGM (Jacquemet and Henriquez, 2009; Verheule and Schotten, 2021). As previously mentioned, the highpass filter cut-off value affects the signal amplitude. In the presence of fibrotic tissue, uEGMs and biEGMs have a different frequency spectrum and are affected in a different manner. Figure 5D shows that there is no optimal cut-off frequency as previously reported by Starreveld et al. (2020). The filtered biEGM amplitude (orange line) drops due to the highpass cut-off but does not intersect the unfiltered amplitude (blue dashed line) as is the case for healthy myocardium (Figure 3D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/908069/fphys-13-908069-HTML-r1/image_m/fphys-13-908069-g004.jpg)
Figure 4. Electrical propagation in fibrotic cardiac tissue, the composition is heterogeneous and includes cardiac myocytes (orange), myofibroblasts/fibroblasts (green) and collagen fibers (purple). Depolarization of cardiomyocytes when an excitation propagates from left to right (brown arrows). Dotted arrows represent a conduction block, while dashed arrows represent slowed conduction. As a result of the zig-zag propagation of the wavefront (white arrows), the unipolar electrogram morphology is not symmetric, is prolonged and shows multiple deflections.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/908069/fphys-13-908069-HTML-r1/image_m/fphys-13-908069-g005.jpg)
Figure 5. Filter effect on electrograms measured in the proximity of a fibrotic area. (A) Unipolar electrograms with different highpass filter values. (B) Effect of different highpass cut-off values on the unipolar electrogram amplitude. (C) Bipolar electrogram with different highpass cut-off values. (D) Effect of different highpass cut-off values on the bipolar electrogram amplitude.
Many approaches have been proposed to model fibrotic cardiac tissue (Table 1) to understand the effect on the wavefront propagation and the corresponding electrograms (Ashihara et al., 2012; McDowell et al., 2013; Roney et al., 2016).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/908069/fphys-13-908069-HTML-r1/image_m/fphys-13-908069-t001.jpg)
Table 1. Different modeling approaches to represent fibrotic tissue in computational models and their effect on simulated EGMs.
Creating a model of cardiac fibrotic tissue is not an easy task as fibrosis formation has been associated with different diseases (myocardial infarction (Liu et al., 2017), diabetes (Russo and Frangogiannis, 2016), autoimmune diseases (Tschöpe et al., 2021) and others), which produce different patterns of structural remodeling (interstitial, compact, diffuse, and patchy) (Nguyen et al., 2014). For example, it has been described that during an ischemic episode in the ventricle, the myocardium undergoes electrical remodeling (Mendonca Costa et al., 2018). From a macroscopic view, conduction velocity is reduced in the scar area, which can be modeled by decreasing the conductivity or by including isolating barriers (Balaban et al., 2018). Additionally, at a cellular scale the cardiac myocytes undergo electrical remodeling (Mendonca Costa et al., 2018). At the border zone of the ischemic area, cardiomyocytes lack oxygen which impacts their metabolism and increase acidity. This triggers a series of effects in the cell’s ion channels. To model these effects, the maximum conductance (Rodriguez et al., 2006; Loewe et al., 2018) of certain ionic channels are modified including an ATP-sensitive potassium channel (IKATP), which has a major contribution during ischemic episodes (Dutta et al., 2017).
Moreover, computational models of pathological tissue need to include fibrosis at the tissue scale. Fibrosis patterns (Figure 6) can be modeled using different approaches by assigning different properties to the mesh using for example, a random distribution (e.g., uniform or Gaussian) (Sánchez et al., 2019b; ten Tusscher and Panfilov, 2007; Alonso and Bär, 2013; Vigmond et al., 2016), by extracting the scar area from MRI (McDowell et al., 2012; Krueger et al., 2014; Morgan et al., 2016; Beach et al., 2021) or by using algorithms that synthetically generate similar patterns as observed in histological cuts of fibrotic tissue (Jakes et al., 2019; Pezzuto et al., 2019; Sutanto et al., 2020; Sánchez et al., 2021b).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/908069/fphys-13-908069-HTML-r1/image_m/fphys-13-908069-g006.jpg)
Figure 6. Sketches of a tissue cut for healthy and fibrotic tissue (top row). Fibrotic sketches represent different fibrotic patterns (diffuse, interstitial, patchy or compact). The bottom row depicts the homogenization assumption where a hexahedral mesh element of 100 μm × 100 μm × 100 μm represents several cardiac myocytes and has to assume one average set of properties that describes the electrophysiology of this group of cells. For fibrotic tissue, homogenization implies that one element contains different types of cells (cardiac myocytes and fibroblasts/myofibroblasts) and collagen. Also the electrophysiological characteristics of this piece of tissue has to be represented by one set of effective parameters.
Fibrosis can be modeled differently and many studies reduce the tissue conductivity such as for example, informed by microstructural modeling in Gokhale et al. (2017). The conductivity of the fibrotic areas can also be reduced in the transversal direction (McDowell et al., 2012) to represent lateralization of gap junctions, close to zero in all directions (Clayton, 2018) or affected by a no flux boundary condition (ten Tusscher and Panfilov, 2007; Alonso and Bär, 2013) to represent replacement fibrosis. The specific spatial distribution of conductivities or conduction velocity can be informed by fibrosis imaging such as the pixel intensity in late gadolinium enhanced magnetic resonance images (Krueger et al., 2014; Morgan et al., 2016; Caixal et al., 2020; Beach et al., 2021) or using a mathematical function determined from EGM amplitude (Lim et al., 2020). Within the regions, either uniform conductivities can be assumed or a gradient from the center of the fibrotic area to the healthy surrounding tissue is assumed.
Furthermore, the edge splitting method has been proposed to separate the computational mesh along its edges with the aim to reproduce the effect of collagen deposition in fibrotic tissue that separates the cardiac myocytes (Mendonca Costa et al., 2014). Edge splitting consists of splitting the nodes along and edge to disconnect adjacent elements creating an alternative path for the wavefront propagation in the cardiac tissue. However, reducing the conductivity or splitting the edges of the mesh does not capture the effect of increased cellular heterogeneity in the cardiac tissue (fibroblast-myocyte coupling) and the inflammatory response. To model cellular heterogeneity, myofibroblast or fibroblast models have been introduced (MacCannell et al., 2007; Ashihara et al., 2012; Morgan et al., 2016; Roney et al., 2016; Sánchez et al., 2021b). Myofibroblasts or fibroblasts were electrically connected to the myocytes by gap junctions. There are equivocal data about the exact conductance of these gap junctions and the number of fibroblasts that a myocyte couples to. In computational models, the value of conductance ranges between 0.5 nS to 2 nS and up to 9 fibroblasts are considered Morgan et al. (2016), MacCannell et al. (2007), Maleckar et al. (2009), Rook et al. (1992), Sánchez et al. (2019a), Seemann et al. (2017). The inflammatory response (myocyte-fibroblast paracrine interactions) has been modeled by altering the maximum conductance of the sodium ion channel (reduced by 50%), the maximum conductance of the L-type calcium ion channel (reduced by 50%), and the maximum conductance of the inward potassium rectifier ion channel (reduced by 40%) (Zahid et al., 2016), as reported by in vitro experiments (Avila et al., 2007; Ramos-Mondragón et al., 2011).
Lately, Vigmond et al. (2016) proposed to represent fibrotic tissue in a monodomain model by removing the elements of the mesh to capture the effect of the low conductive extracellular medium and the absence of an intracellular current path. One advantage of the proposed modeling approach, is that there is no flux of current towards the fibrotic tissue; therefore, there are no source elements that will contribute to the calculated extracellular potential. Using this modeling approach, the authors observed that at the percolation threshold (Alonso and Bär, 2013) the fibrotic tissue was be able to trigger and maintain an arrhythmia. The EGMs calculated over the fibrotic tissue exhibit fractionation due to the zig-zag patterns of depolarization in the cardiac tissue in this modeling approach. Moreover, the study also looked at the impact of the mesh resolution when modeling fibrotic tissue and showed that in meshes with a resolution of 300 μm conduction block was reached at lower degrees of fibrosis than in meshes with finer resolution (
Using a realistic geometry Jacquemet et al. (2003) studied the morphology of uEGMs during different atrial fibrillation propagation patterns. The authors showed that different propagation patterns (plane waves, spiral waves, and wavefront collision) lead to different uEGM morphology (symmetry and amplitude) and that asymmetric signals (Figure 2B) occurred in less than 2% of the cases in homogeneous substrate. However, the increase of heterogeneities in the cardiac tissue also increases the asymmetry and reduces the amplitude of the EGM(van der Does and de Groot, 2017). Frontera et al. (2018) showed how different depolarization patterns affected the biEGM morphology. High peak-to-peak amplitude and short duration of biEGMs are wavefront collisions or pivotal points, low peak-to-peak amplitude and EGM prolongation are associated with slow conduction areas. The authors remarked how understanding the genesis of the electrograms is a key factor to improving the arrhythmia treatments.
Including heterogeneous tissue composition in the model changes the wavefront propagation in the cardiac tissue (McDowell et al., 2012; Campos et al., 2013; Mendonca Costa et al., 2014; Roney et al., 2016) (Figure 7). Roney et al. (2016) showed how different modeling approaches of cardiac fibrosis can change the propagation in the cardiac tissue and affect the morphology of EGMs. In that study, Roney et al. (2016) modeled fibrosis as conduction disturbances (lower conductivity, edge splitting, or removing elements). They included electrical remodeling of the cardiac myocyte due to inflammatory processes mediated by transforming growth factor-β1, myocyte-fibroblast coupling and combinations of the preceding. EGM morphology was mostly affected when fibrosis was modeled by edge splitting or removing the elements (Figure 7) as also shown previously. In addition, including fibroblast coupling has an organizing effect on rotor dynamics, also shown by other studies (McDowell et al., 2012; Sánchez et al., 2019b).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/908069/fphys-13-908069-HTML-r1/image_m/fphys-13-908069-g007.jpg)
Figure 7. Different fibrosis modeling methodologies and their corresponding unipolar and bipolar electrograms. Central fibrotic area (dashed line), the choice of the modeling approach affects the resulting uEGMs and biEGMs. (A) fibroblasts coupled to myocytes in the fibrotic region; (B) removing a share of elements in the fibrotic region from the computational domain; (C) severely reducing conductivity in a share of elements in the fibrotic region but the spatial transmembrane voltage gradients of fibrotic elements still contribute to the EGMs; (D) conductivity gradient from center of the fibrotic region to the surrounding healthy tissue. Absolute EGM amplitudes are smaller than in vivo due to the small size of the tissue patch.
The amplitude of the EGMs can also be affected by conduction impairment along certain axes (McDowell et al., 2012; Gharaviri et al., 2016; Irakoze and Jacquemet, 2020). Gharaviri et al. (2016) created a model of the cardiac tissue that enables the study of dissociation between transmural layers, for example, dissociation between the subendocardial and the subepicardial myocardium as can be caused by endomysial fibrosis. Moreover, Saba et al. (2009) described how the epicardial EGM amplitude varies in the ventricle with the thickness of the epicardial fat layer. The authors showed that biEGM amplitude was inversely related to epicardial fat thickness. Thus, using a voltage cut-off of 0.5 mV to define scar tissue would lead to identifying also healthy areas with overlaying fat and more information needs to be used to define epicardial tissue characteristics.
5 Other Factors Impacting Intracardiac Signals
EGM morphology and amplitude are also affected by electrode polarization, excessive contact pressure, catheter motion (Oesterlein et al., 2016), electromagnetic interference (Unger et al., 2019), near field and far field effects (Schicketanz et al., 2021), and poor grounding. However, most in silico experiments do not consider these factors, which might alter the EGM characteristics. Simulation studies have created a model of clinical noise which covers the electromagnetic interference (Sánchez et al., 2021a; Nothstein et al., 2021). However, further aspects likely need to be considered explicitly if their influence is relevant for the intended use of the model.
6 Research Gaps and Potential Future Developments
Modeling of the cardiac tissue has significantly advanced understanding of the electrical propagation and the measured intracardiac EGMs. There is consensus on how to assign the properties of the computational model to represent healthy myocardium and the advantages and limitations of the different approaches to compute the extracellular potentials are mostly characterized. However, the question how to model fibrosis is far from being ultimately answered and will most likely continue to depend on the question of interest to be answered with a specific model. Additionally, the mesh resolution used in most of the studies of ≈300 μm determines the degree of homogenization (Figure 6). Spatial discretization of the mesh at the cellular level should be considered to study the influence of microstructural heterogeneity in the tissue (e.g., fibrosis) on EGMs (Figure 4). In addition, such models with subcellular resolution would enable to investigate to which degree discontinuous propagation within a cell vs. between cells leads to fractionation in healthy tissue. Here, we presented an overview of the commonly used methods and their corresponding EGMs.
Over the last years, the human cardiac digital twin has been under development to suggest personalized treatments for cardiac arrhythmias. Gillette et al. (2021) proposed an automated framework to generate a patient’s digital twin from clinical data and Nagel et al. (2021) proposed a statistical approach to generate a population of anatomical models. While the anatomical model can be accurately generated from magnetic resonance images or statistical shape models, functional twinning can be achieved by tuning a phenomenological model or using generalized global properties for the cardiac tissue. Functional information will impact the morphology and amplitude of the EGM. However, over the years, different studies proposed distinct methodologies to extract structural and functional information from the EGM signals. One open question is still the possibility of obtaining repolarization times from EGMs as repolarization of the cardiac tissue plays a pivotal role for the initiation of arrhythmias (Rivaud et al., 2021). From simulations of atrial electrophysiology, Celotto et al. (2021) proposed a method to detect areas of parasympathetic innervation from the amplitude of the repolarization EGM. Verrier et al. (2016) showed that repolarization times can be recovered from EGMs for both the atrium and the ventricle in a controlled clinical environment. However, initial experience in other groups including our own suggest that reliably obtaining atrial repolarization information from EGMs remains a challenge.
Different studies demonstrated discrepancies when using the same voltage threshold (for example, 0.5 mV) to distinguish healthy from pathological tissue when mapping during different rhythms (sinus rhythm and AF) (Rodríguez-Mañero et al., 2018; Nairn et al., 2020b; Nairn et al., 2022). Nairn et al. (2022) looked at how the amplitude of EGMs changed when electroanatomical mapping was performed under three different rhythms (sinus rhythm, native AF, and induced AF). The authors proposed not only one single cut-off voltage value for the entire atrium but regional voltage thresholds to minimize the discrepancies between different mapping rhythms. Computer models could help to further characterize the voltage relations during different rhythms and to overcome the use of a voltage threshold to distinguish the cardiac substrate (healthy and fibrotic) by combining in vivo data and in silico data to fully exploit the information contained in EGMs (Sánchez et al., 2021a). Additionally, computer models of cardiac electrophysiology could aid the design of medical devices helping in understanding the factors that affect EGMs to raise awareness for them (Oesterlein et al., 2016; Pollnow et al., 2017; Beheshti et al., 2018; Hwang et al., 2019) as well as to inform the choice of parameters to improve the technologies as proposed for cardiac resynchronization therapy (Jolley et al., 2010).
Understanding the functional relationship between the discrete structure and continuum behaviour of cardiac tissue at microscopic and macroscopic levels is a significant challenge (Gokhale et al., 2017). At the microscopic level, Tveito et al. (2017) and Bécue et al. (2017) proposed a cell-by-cell approach that explicitly models the extracellular, membrane and intracellular domain. However, cell-by-cell models are computationally expensive and will require an increase of computational resources such that finer meshes up to cellular resolution can be handled efficiently (Potse et al., 2020). At the macroscopic level, reduced order models (Fresca et al., 2020) could help to reproduce in detail the electrophysiology of the cardiac tissue without losing important details that will determine the vulnerability of the tissue to arrhythmia. Recently, (Herrero Martin et al., 2022) explored the use of Physics Informed Neural Networks (PINN) to model the electrical propagation in the cardiac tissue. The authors introduced electrophysiology models to the neural network and were able to reconstruct the spatial-temporal dynamics of the action potential and its propagation. One of the big drawbacks of these approaches is the amount of data needed to train the network in order to predict different possible propagations patterns.
Software plays a fundamental role in cardiac modeling. Recent work demonstrated significant speedup of simulations of cardiac electrophysiology (Sundnes et al., 2006; Seemann et al., 2010; Cooper et al., 2015; Quarteroni et al., 2017; Sánchez et al., 2020; Plank et al., 2021). However, it remains to be seen how effectively GPUs can be integrated into large-scale cardiac simulations. Regardless, several numerical libraries are currently available, opening the door to accelerate cardiac electrophysiology simulations (Anzt et al., 2020; Mills et al., 2021).
7 Conclusion
Models of cardiac tissue electrophysiology have played an essential role in advancing our understanding of action potential propagation in the heart and the genesis of EGMs. Despite the significant progress of different modeling approaches and efficient numerical software, there are substantial challenges, such as modeling of the microstructure at a close-to-cellular scale, modeling the different aspects of fibrosis, electrophysiological heterogeneity as well as realistic electrode configurations. Dedicated simulation studies with refined models will help to further elucidate the different factors that contribute to EGM genesis and impact their morphology.
Author Contributions
JS and AL contributed to conception and design of the study. JS drafted the manuscript. All authors contributed to manuscript revision, read and approved the submitted version.
Funding
This work was supported by the European High-Performance Computing Joint Undertaking EuroHPC under grant agreement No 955495 (MICROCARD) co-funded by the Horizon 2020 programme of the European Union (EU), the French National Research Agency ANR, the German Federal Ministry of Education and Research, the Italian ministry of economic development, the Swiss State Secretariat for Education, Research and Innovation, the Austrian Research Promotion Agency FFG, and the Research Council of Norway.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We gratefully acknowledge the state of Baden-Württemberg through bwHPC and the KIT-Publication Fund of the Karlsruhe Institute of Technology. In addition, the authors thank Deborah Nairn, Joshua Frederic Steyer, and Mark Potse for their valuable suggestions.
References
Abdi B., Hendriks R. C., van der Veen A.-J., de Groot N. M. S. (2020). Improved Local Activation Time Annotation of Fractionated Atrial Electrograms for Atrial Mapping. Comput. Biol. Med. 117, 103590. doi:10.1016/j.compbiomed.2019.103590
Alonso S., Bär M. (2013). Reentry Near the Percolation Threshold in a Heterogeneous Discrete Model for Cardiac Tissue. Phys. Rev. Lett. 110, 158101. doi:10.1103/PhysRevLett.110.158101
Anzt H., Cojean T., Flegar G., Göbel F., Grützmacher T., Nayak P., et al. (2020). Ginkgo: A Modern Linear Operator Algebra Framework for High Performance Computing. Available at: https://arxiv.org/abs/2006.16852
Ashihara T., Haraguchi R., Nakazawa K., Namba T., Ikeda T., Nakazawa Y., et al. (2012). The Role of Fibroblasts in Complex Fractionated Electrograms During Persistent/Permanent Atrial Fibrillation. Circ. Res. 110, 275–284. doi:10.1161/CIRCRESAHA.111.255026
Avila G., Medina I. M., Jiménez E., Elizondo G., Aguilar C. I. (2007). Transforming Growth Factor-Β1 Decreases Cardiac Muscle L-type Ca2+ Current and Charge Movement by Acting on the Cav1.2 mRNA. Am. J. Physiology-Heart Circulatory Physiology 292, H622–H631. doi:10.1152/ajpheart.00781.2006
Bécue P.-E., Potse M., Coudière Y. (2017). A Three-Dimensional Computational Model of Action Potential Propagation through a Network of Individual Cells. Computing in Cardiology (CinC), 1–4. doi:10.22489/CinC.2017.022-271
Balaban G., Halliday B. P., Mendonca Costa C., Bai W., Porter B., Rinaldi C. A., et al. (2018). Fibrosis Microstructure Modulates Reentry in Non-Ischemic Dilated Cardiomyopathy: Insights From Imaged Guided 2D Computational Modeling. Front. Physiol. 9, 1–13. doi:10.3389/fphys.2018.01832
Beach M., Sim I., Mehta A., Kotadia I., O’Hare D., Whitaker J., et al. (2021). Using the Universal Atrial Coordinate System for MRI and Electroanatomic Data Registration in Patient-Specific Left Atrial Model Construction and Simulation. 11th International Conference on Functional Imaging and Modeling of the Heart. 629–638. doi:10.1007/978-3-030-78710-3_60
Beheshti M., Magtibay K., Massé S., Porta-Sanchez A., Haldar S., Bhaskaran A., et al. (2018). Determinants of Atrial Bipolar Voltage: Inter Electrode Distance and Wavefront Angle. Comput. Biol. Med. 102, 449–457. doi:10.1016/j.compbiomed.2018.07.011
Bishop M. J., Plank G. (2011a). Bidomain ECG Simulations Using an Augmented Monodomain Model for the Cardiac Source. IEEE Trans. Biomed. Eng. 58, 2297–2307. doi:10.1109/TBME.2011.2148718
Bishop M. J., Plank G. (2011b). Representing Cardiac Bidomain Bath-Loading Effects by an Augmented Monodomain Approach: Application to Complex Ventricular Models. IEEE Trans. Biomed. Eng. 58, 1066–1075. doi:10.1109/TBME.2010.2096425
Cabrera-Lozoya R., Berte B., Cochet H., Jais P., Ayache N., Sermesant M. (2017). Image-Based Biophysical Simulation of Intracardiac Abnormal Ventricular Electrograms. IEEE Trans. Biomed. Eng. 64, 1446–1454. doi:10.1109/TBME.2016.2562918
Caixal G., Alarcón F., Althoff T. F., Nuñez-Garcia M., Benito E. M., Borràs R., et al. (2020). Accuracy of Left Atrial Fibrosis Detection with Cardiac Magnetic Resonance: Correlation of Late Gadolinium Enhancement with Endocardial Voltage and Conduction Velocity. EP Eur. 23, 380–388. doi:10.1093/europace/euaa313
Campos F. O., Wiener T., Prassl A. J., Weber dos Santos R., Sanchez-Quintana D., Ahammer H., et al. (2013). Electroanatomical Characterization of Atrial Microfibrosis in a Histologically Detailed Computer Model. IEEE Trans. Biomed. Eng. 60, 2339–2349. doi:10.1109/TBME.2013.2256359
Celotto C., Sánchez C., Mountris K. A., Laguna P., Pueyo E. (2021). Location of Parasympathetic Innervation Regions From Electrograms to Guide Atrial Fibrillation Ablation Therapy: An In Silico Modeling Study. Front. Physiol. 12. Available at: https://www.frontiersin.org/articles/10.3389/fphys.2021.674197/full. doi:10.3389/fphys.2021.674197
Clayton R. H. (2018). Dispersion of Recovery and Vulnerability to Re-Entry in a Model of Human Atrial Tissue with Simulated Diffuse and Focal Patterns of Fibrosis. Front. Physiol. 9, 1052. doi:10.3389/fphys.2018.01052
Clerc L. (1976). Directional Differences of Impulse Spread in Trabecular Muscle from Mammalian Heart. J. Physiology 255, 335–346. doi:10.1113/jphysiol.1976.sp011283
Colli Franzone P., Pavarino L. F., Scacchi S., Taccardi B. (2007). Monophasic Action Potentials Generated by Bidomain Modeling as a Tool for Detecting Cardiac Repolarization Times. Am. J. Physiology-Heart Circulatory Physiology 293, H2771–H2785. doi:10.1152/ajpheart.00651.2007
Cooper J., Spiteri R. J., Mirams G. R. (2015). Cellular Cardiac Electrophysiology Modeling with Chaste and CellML. Front. Physiol. 5, 511. doi:10.3389/fphys.2014.00511
Dutta S., Mincholé A., Quinn T. A., Rodriguez B. (2017). Electrophysiological Properties of Computational Human Ventricular Cell Action Potential Models under Acute Ischemic Conditions. Prog. Biophysics Mol. Biol. 129, 40–52. doi:10.1016/j.pbiomolbio.2017.02.007
Fresca S., Manzoni A., Dedè L., Quarteroni A. (2020). Deep Learning-Based Reduced Order Models in Cardiac Electrophysiology. Plos One 15, e0239416. doi:10.1371/journal.pone.0239416
Frontera A., Takigawa M., Martin R., Thompson N., Cheniti G., Massoullié G., et al. (2018). Electrogram Signature of Specific Activation Patterns: Analysis of Atrial Tachycardias at High-Density Endocardial Mapping. Heart rhythm. 15, 28–37. doi:10.1016/j.hrthm.2017.08.001
Gaeta S., Bahnson T. D., Henriquez C. (2020). Mechanism and Magnitude of Bipolar Electrogram Directional Sensitivity: Characterizing Underlying Determinants of Bipolar Amplitude. Heart rhythm. 17, 777–785. doi:10.1016/j.hrthm.2019.12.010
Ganesan A. N., Kuklik P., Lau D. H., Brooks A. G., Baumert M., Lim W. W., et al. (2013). Bipolar Electrogram Shannon Entropy at Sites of Rotational Activation. Circ Arrhythmia Electrophysiol. 6, 48–57. doi:10.1161/CIRCEP.112.976654
Gharaviri A., Verheule S., Eckstein J., Potse M., Kuklik P., Kuijpers N. H., et al. (2016). How Disruption of Endo-Epicardial Electrical Connections Enhances Endo-Epicardial Conduction during Atrial Fibrillation. Europace 19 (2), 308–318. doi:10.1093/europace/euv445
Gillette K., Gsell M. A. F., Prassl A. J., Karabelas E., Reiter U., Reiter G., et al. (2021). A Framework for the Generation of Digital Twins of Cardiac Electrophysiology from Clinical 12-leads ECGs. Med. Image Anal. 71, 102080. doi:10.1016/j.media.2021.102080
Gima K., Rudy Y. (2002). Ionic Current Basis of Electrocardiographic Waveforms. Circulation Res. 90, 889–896. doi:10.1161/01.RES.0000016960.61087.86
Gokhale T. A., Medvescek E., Henriquez C. S. (2017). Modeling Dynamics in Diseased Cardiac Tissue: Impact of Model Choice. Chaos 27, 093909. doi:10.1063/1.4999605
Herrero Martin C., Oved A., Chowdhury R. A., Ullmann E., Peters N. S., Bharath A. A., et al. (2022). EP-PINNs: Cardiac Electrophysiology Characterisation Using Physics-Informed Neural Networks. Front. Cardiovasc. Med. 8. doi:10.3389/fcvm.2021.768419
Hodgkin A. L., Huxley A. F. (1952). A Quantitative Description of Membrane Current and its Application to Conduction and Excitation in Nerve. J. Physiology 117, 500–544. doi:10.1113/jphysiol.1952.sp004764
Hwang M., Kim J., Lim B., Song J.-S., Joung B., Shim E. B., et al. (2019). Multiple Factors Influence the Morphology of the Bipolar Electrogram: An In Silico Modeling Study. PLoS Comput. Biol. 15, e1006765. doi:10.1371/journal.pcbi.1006765
Irakoze É., Jacquemet V. (2020). Breakthrough Wave Detection in a 3D Computer Model of Atrial Endo-Epicardial Dissociation. 2020 Computing in Cardiology Conference. doi:10.22489/CinC.2020.425
Jacquemet V., Henriquez C. S. (2009). Genesis of Complex Fractionated Atrial Electrograms in Zones of Slow Conduction: A Computer Model of Microfibrosis. Heart rhythm. 6, 803–810. doi:10.1016/j.hrthm.2009.02.026
Jacquemet V., Virag N., Ihara Z., Dang L., Blanc O., Zozor S., et al. (2003). Study of Unipolar Electrogram Morphology in a Computer Model of Atrial Fibrillation. J. Cardiovasc. Electrophysiol. 14, S172–S179. doi:10.1046/j.1540.8167.90308.x
Jadidi A., Nothstein M., Chen J., Lehrmann H., Dössel O., Allgeier J., et al. (2020). Specific Electrogram Characteristics Identify the Extra-pulmonary Vein Arrhythmogenic Sources of Persistent Atrial Fibrillation - Characterization of the Arrhythmogenic Electrogram Patterns During Atrial Fibrillation and Sinus Rhythm. Sci. Rep. 10, 9147. doi:10.1038/s41598-020-65564-2
Jakes D., Burrage K., Drovandi C. C., Burrage P., Bueno-Orovio A., Santos R. W. d., et al. (2019). Perlin noise generation of physiologically realistic patterns of fibrosis. bioRxiv. doi:10.1101/668848
Jolley M., Stinstra J., Tate J., Pieper S., MacLeod R., Chu L., et al. (2010). Finite Element Modeling of Subcutaneous Implantable Defibrillator Electrodes in an Adult Torso. Heart rhythm. 7, 692–698. doi:10.1016/j.hrthm.2010.01.030
Keller M. W., Schuler S., Seemann G., Dössel O. (2012). Differences in Intracardiac Signals on a Realistic Catheter Geometry Using Mono- and Bidomain Models. Comput. Cardiol., 305–308. Available at: http://cinc.mit.edu/archives/2012/pdf/0305.pdf
Keller M. W., Schuler S., Wilhelms M., Lenis G., Seemann G., Schmitt C., et al. (2014). Characterization of Radiofrequency Ablation Lesion Development Based on Simulated and Measured Intracardiac Electrograms. IEEE Trans. Biomed. Eng. 61, 2467–2478. doi:10.1109/TBME.2014.2322515
Krueger M. W., Rhode K. S., O‘Neill M. D., Rinaldi C. A., Gill J., Razavi R., et al. (2014). Patient-specific Modeling of Atrial Fibrosis Increases the Accuracy of Sinus Rhythm Simulations and May Explain Maintenance of Atrial Fibrillation. J. Electrocardiol. 47, 324–328. doi:10.1016/j.jelectrocard.2013.11.003
Lim B., Kim J., Hwang M., Song J.-S., Lee J. K., Yu H.-T., et al. (2020). In Situ procedure for High-Efficiency Computational Modeling of Atrial Fibrillation Reflecting Personal Anatomy, Fiber Orientation, Fibrosis, and Electrophysiology. Sci. Rep. 10, 2417. doi:10.1038/s41598-020-59372-x
Lin Y.-J., Tai C.-T., Lo L.-W., Udyavar A. R., Chang S.-L., Wongcharoen W., et al. (2007). Optimal Electrogram Voltage Recording Technique for Detecting the Acute Ablative Tissue Injury in the Human Right Atrium. J. Cardiovasc Electrophysiol. 18, 617–622. doi:10.1111/j.1540-8167.2007.00803.x
Liu T., Song D., Dong J., Zhu P., Liu J., Liu W., et al. (2017). Current Understanding of the Pathophysiology of Myocardial Fibrosis and its Quantitative Assessment in Heart Failure. Front. Physiol. 8, 238. doi:10.3389/fphys.2017.00238
Loewe A., Wülfers E. M., Seemann G. (2018). Cardiac Ischemia-Insights from Computational Models. Herzschr Elektrophys 29, 48–56. doi:10.1007/s00399-017-0539-6
MacCannell K., Bazzazi H., Chilton L., Shibukawa Y., Clark R. B., Giles W. R. (2007). A Mathematical Model of Electrotonic Interactions Between Ventricular Myocytes and Fibroblasts. Biophysical J. 92, 4121–4132. doi:10.1529/biophysj.106.101410
Maleckar M. M., Greenstein J. L., Giles W. R., Trayanova N. A. (2009). Electrotonic Coupling between Human Atrial Myocytes and Fibroblasts Alters Myocyte Excitability and Repolarization. Biophysical J. 97, 2179–2190. doi:10.1016/j.bpj.2009.07.054
Malmivuo J., Plonsey R. (1995). BioelectromagnetismPrinciples and Applications of Bioelectric and Biomagnetic Fields. Oxford: Oxford University Press.
McDowell K. S., Vadakkumpadan F., Blake R., Blauer J., Plank G., Macleod R. S., et al. (2013). Mechanistic Inquiry into the Role of Tissue Remodeling in Fibrotic Lesions in Human Atrial Fibrillation. Biophysical J. 104, 2764–2773. doi:10.1016/j.bpj.2013.05.025
McDowell K. S., Vadakkumpadan F., Blake R., Blauer J., Plank G., MacLeod R. S., et al. (2012). Methodology for Patient-Specific Modeling of Atrial Fibrosis as a Substrate for Atrial Fibrillation. J. Electrocardiol. 45, 640–645. doi:10.1016/j.jelectrocard.2012.08.005
Mendonca Costa C., Campos F. O., Prassl A. J., Weber dos Santos R., Sanchez-Quintana D., Ahammer H., et al. (2014). An Efficient Finite Element Approach for Modeling Fibrotic Clefts in the Heart. IEEE Trans. Biomed. Eng. 61, 900–910. doi:10.1109/TBME.2013.2292320
Mendonca Costa C., Plank G., Rinaldi C. A., Niederer S. A., Bishop M. J. (2018). Modeling the Electrophysiological Properties of the Infarct Border Zone. Front. Physiol. 9, 356. doi:10.3389/fphys.2018.00356
Mills R. T., Adams M. F., Balay S., Brown J., Dener A., Knepley M., et al. (2021). Toward Performance-Portable PETSc for GPU-Based Exascale Systems. Parallel Comput. 108, 102831. doi:10.1016/j.parco.2021.102831
Morgan R., Colman M. A., Chubb H., Seemann G., Aslanidi O. V. (2016). Slow Conduction in the Border Zones of Patchy Fibrosis Stabilizes the Drivers for Atrial Fibrillation: Insights from Multi-Scale Human Atrial Modeling. Front. Physiol. 7, 1–15. doi:10.3389/fphys.2016.00474
Nagel C., Espinosa C. B., Gillette K., Gsell M. A. F., Sánchez J., Plank G., et al. (2022). Comparison of Propagation Models and Forward Calculation Methods on Cellular, Tissue and Organ Scale Atrial Electrophysiology. Available at: https://arxiv.org/abs/2203.07776
Nagel C., Schuler S., Dössel O., Loewe A. (2021). A Bi-atrial Statistical Shape Model for Large-Scale In Silico Studies of Human Atria: Model Development and Application to ECG Simulations. Med. Image Anal. 74, 102210. doi:10.1016/j.media.2021.102210
Nairn D., Eichenlaub M., Lehrmann H., Müller-Edenborn B., Chen J., Huang T., et al. (2022). Spatial Correlation of Left Atrial Low Voltage Substrate in Sinus Rhythm versus Atrial Fibrillation: Identifying the Pathological Substrate Irrespective of the Rhythm. medRxiv. doi:10.1101/2022.02.18.22271172
Nairn D., Hunyar D., Sánchez J., Doessel O., Loewe A. (2020a). Impact of Electrode Size on Electrogram Voltage in Healthy and Diseased Tissue. 2020 Computing in Cardiology Conference. doi:10.22489/CinC.2020.146
Nairn D., Lehrmann H., Müller-Edenborn B., Schuler S., Arentz T., Dössel O., et al. (2020b). Comparison of Unipolar and Bipolar Voltage Mapping for Localization of Left Atrial Arrhythmogenic Substrate in Patients with Atrial Fibrillation. Front. Physiol. 11, 575846. doi:10.3389/fphys.2020.575846
Neic A., Campos F. O., Prassl A. J., Niederer S. A., Bishop M. J., Vigmond E. J., et al. (2017). Efficient Computation of Electrograms and ECGs in Human Whole Heart Simulations Using a Reaction-Eikonal Model. J. Comput. Phys. 346, 191–211. doi:10.1016/j.jcp.2017.06.020
Nguyen T. P., Qu Z., Weiss J. N. (2014). Cardiac Fibrosis and Arrhythmogenesis: The Road to Repair Is Paved with Perils. J. Mol. Cell. Cardiol. 70, 83–91. doi:10.1016/j.yjmcc.2013.10.018
Nothstein M., Luik A., Jadidi A., Sánchez J., Unger L. A., Wülfers E. M., et al. (2021). CVAR-Seg: An Automated Signal Segmentation Pipeline for Conduction Velocity and Amplitude Restitution. Front. Physiol. 12, 673047. doi:10.3389/fphys.2021.673047
Oesterlein T., Frisch D., Loewe A., Seemann G., Schmitt C., Dössel O., et al. (2016). Basket-Type Catheters: Diagnostic Pitfalls Caused by Deformation and Limited Coverage. BioMed Res. Int. 2016, 1–13. doi:10.1155/2016/5340574
Pezzuto S., Quaglino A., Potse M. (2019). “On Sampling Spatially-Correlated Random Fields for Complex Geometries,” in Functional Imaging and Modeling of the Heart (Cham: Springer International Publishing). doi:10.1007/978-3-030-21949-9_12
Plank G., Loewe A., Neic A., Augustin C., Huang Y. L., Gsell M. A. F., et al. (2021). The openCARP Simulation Environment for Cardiac Electrophysiology. Comput. Methods Programs Biomed. 208, 106223. doi:10.1016/j.cmpb.2021.106223
Pollnow S., Greiner J., Oesterlein T., Wülfers E. M., Loewe A., Dössel O. (2017). Mini Electrodes on Ablation Catheters: Valuable Addition or Redundant Information?-Insights from a Computational Study. Comput. Math. Methods Med. 2017, 1–13. doi:10.1155/2017/1686290
Potse M., Dube B., Vinet A., Cardinal R. (2006). “A Comparison of Monodomain and Bidomain Propagation Models for the Human Heart,” in 2006 International Conference of the IEEE Engineering in Medicine and Biology Society (IEEE), 3895–3898. doi:10.1109/IEMBS.2006.259484
Potse M., Saillard E., Barthou D., Coudière Y. (2020). Feasibility of Whole-Heart Electrophysiological Models with Near-Cellular Resolution. 2020 Computing in Cardiology. doi:10.22489/CinC.2020.126
Quarteroni A., Manzoni A., Vergara C. (2017). The Cardiovascular System: Mathematical Modelling, Numerical Algorithms and Clinical Applications. Acta Numer. 26, 365–590. doi:10.1017/S0962492917000046
Ramos-Mondragón R., Vega A. V., Avila G. (2011). Long-term Modulation of Na+ and K+ Channels by TGF-Β1 in Neonatal Rat Cardiac Myocytes. Pflugers Arch. - Eur. J. Physiol. 461, 235–247. doi:10.1007/s00424-010-0912-3
Rivaud M. R., Bayer J. D., Cluitmans M., van der Waal J., Bear L. R., Boukens B. J., et al. (2021). Critical repolarization gradients determine the induction of reentry-based torsades de pointes arrhythmia in models of long QT syndrome. Heart rhythm 18, 278–287. doi:10.1016/j.hrthm.2020.09.020
Rodriguez B., Trayanova N., Noble D. (2006). Modeling Cardiac Ischemia. Ann. N. Y. Acad. Sci. 1080, 395–414. doi:10.1196/annals.1380.029
Rodríguez-Mañero M., Valderrábano M., Baluja A., Kreidieh O., Martínez-Sande J. L., García-Seara J., et al. (2018). Validating Left Atrial Low Voltage Areas During Atrial Fibrillation and Atrial Flutter Using Multielectrode Automated Electroanatomic Mapping. JACC Clin. Electrophysiol. 4, 1541–1552. doi:10.1016/j.jacep.2018.08.015
Roney C. H., Bayer J. D., Zahid S., Meo M., Boyle P. M. J., Trayanova N. A., et al. (2016). Modelling Methodology of Atrial Fibrosis Affects Rotor Dynamics and Electrograms. EP Eur. 18, iv146–iv155. doi:10.1093/europace/euw365
Roney C. H., Bendikas R., Pashakhanloo F., Corrado C., Vigmond E. J., McVeigh E. R., et al. (2021). Constructing a Human Atrial Fibre Atlas. Ann. Biomed. Eng. 49, 233–250. doi:10.1007/s10439-020-02525-w
Rook M. B., Van Ginneken A. C., de Jonge B., el Aoumari A., Gros D., Jongsma H. J., et al. (1992). Differences in Gap Junction Channels between Cardiac Myocytes, Fibroblasts, and Heterologous Pairs. Am. J. Physiology-Cell Physiology 263, C959–C977. doi:10.1152/ajpcell.1992.263.5.C959
Russo I., Frangogiannis N. G. (2016). Diabetes-associated Cardiac Fibrosis: Cellular Effectors, Molecular Mechanisms and Therapeutic Opportunities. J. Mol. Cell. Cardiol. 90, 84–93. doi:10.1016/j.yjmcc.2015.12.011
Saba M. M., Akella J., Gammie J., Poston R., Johnson A., Hood R. E., et al. (2009). The Influence of Fat Thickness on the Human Epicardial Bipolar Electrogram Characteristics: Measurements on Patients Undergoing Open-Heart Surgery. Europace 11, 949–953. doi:10.1093/europace/eup156
Sánchez J., Gomez J. F., Martinez-Mateu L., Romero L., Saiz J., Trenor B. (2019a). Heterogeneous Effects of Fibroblast-Myocyte Coupling in Different Regions of the Human Atria under Conditions of Atrial Fibrillation. Front. Physiol. 10, 847. doi:10.3389/fphys.2019.00847
Sánchez J., Nothstein M., Unger L., Saiz J., Trénor B., Dössel O., et al. (2019b). “Influence of Fibrotic Tissue Arrangement on Intracardiac Electrograms During Persistent Atrial Fibrillation,” in 2019 Computing in Cardiology (CinC), 1. doi:10.22489/cinc.2019.342
Sánchez J., Luongo G., Nothstein M., Unger L. A., Saiz J., Trenor B., et al. (2021a). Using Machine Learning to Characterize Atrial Fibrotic Substrate From Intracardiac Signals with a Hybrid In Silico and In Vivo Dataset. Front. Physiol. 12, 699291. doi:10.3389/fphys.2021.699291
Sánchez J., Trenor B., Saiz J., Dössel O., Loewe A. (2021b). Fibrotic Remodeling during Persistent Atrial Fibrillation: In Silico Investigation of the Role of Calcium for Human Atrial Myofibroblast Electrophysiology. Cells 10, 2852. doi:10.3390/cells10112852
Sánchez J., Nothstein M., Neic A., Huang Y.-L., J. Prassl A., Ulrich R., et al. (2020). “openCARP: An Open Sustainable Framework for In-Silico Cardiac Electrophysiology Research,” in 2020 Computing in Cardiology Conference (CinC). doi:10.22489/CinC.2020.111
Schicketanz L., Unger L. A., Sánchez J., Dössel O., Loewe A. (2021). Separating Atrial Near Fields and Atrial Far Fields in Simulated Intra-atrial Electrograms. Curr. Dir. Biomed. Eng. 7, 175–178. doi:10.1515/cdbme-2021-2045
Schmitt O. H. (1969). “Biological Information Processing Using the Concept of Interpenetrating Domains,” in Information Processing in The Nervous System: Proceedings of a Symposium held at the State University of New York at Buffalo. Editor K. N. Leibovic (Berlin, Heidelberg: Springer Berlin Heidelberg), 325–331. doi:10.1007/978-3-642-87086-6_18
Schneider M. A. E., Ndrepepa G., Weber S., Deisenhofer I., Schomig A., Schmitt C. (2004). Influence of High-Pass Filtering on Noncontact Mapping and Ablation of Atrial Tachycardias. Pacing Clin. Electrophysiol. 27, 38–46. doi:10.1111/j.1540-8159.2004.00383.x
Schuler S., Keller M. W., Oesterlein T., Seemann G., Dössel O. (2013). Influence of Catheter Orientation, Tissue Thickness and Conduction Velocity on the Intracardiac Electrogram. Biomed. Eng./Biomed. Tech. doi:10.1515/bmt-2013-4334
Seemann G., Bustamante P., Ponto S., Wilhelms M., Scholz E., Dossel O. (2010). “Atrial Fibrillation-Based Electrical Remodeling in a Computer Model of the Human Atrium,” in 2010 Computers in Cardiology, 417–420.
Seemann G., Loewe A., Wülfers E. M. (2017). “Effects of Fibroblasts Coupling on the Electrophysiology of Cardiomyocytes from Different Regions of the Human Atrium: a Simulation Study,” in 2017 Computing in Cardiology (CinC). doi:10.22489/CinC.2017.380-451
Starreveld R., Knops P., Roos-Serote M., Kik C., Bogers A. J. J. C., Brundel B. J. J. M., et al. (2020). The Impact of Filter Settings on Morphology of Unipolar Fibrillation Potentials. J. Cardiovasc. Trans. Res. 13, 953–964. doi:10.1007/s12265-020-10011-w
Sundnes J., Nielsen B. F., Mardal K. A., Cai X., Lines G. T., Tveito A. (2006). On the Computational Complexity of the Bidomain and the Monodomain Models of Electrophysiology. Ann. Biomed. Eng. 34, 1088–1097. doi:10.1007/s10439-006-9082-z
Sutanto H., Cluitmans M. J. M., Dobrev D., Volders P. G. A., Bébarová M., Heijman J. (2020). Acute Effects of Alcohol on Cardiac Electrophysiology and Arrhythmogenesis: Insights from Multiscale In Silico Analyses. J. Mol. Cell. Cardiol. 146, 69–83. doi:10.1016/j.yjmcc.2020.07.007
Takigawa M., Kitamura T., Basu S., Bartal M., Martin C. A., Martin R., et al. (2022). Effect of Electrode Size and Spacing on Electrograms: Optimized Electrode Configuration for Near-Field Electrogram Characterization. Heart rhythm. 19, 102–112. doi:10.1016/j.hrthm.2021.09.011
ten Tusscher K. H. W. J., Panfilov A. V. (2007). Influence of Diffuse Fibrosis on Wave Propagation in Human Ventricular Tissue. EP Eur. 9, vi38–vi45. doi:10.1093/europace/eum206
Trew M. L., Smaill B. H., Bullivant D. P., Hunter P. J., Pullan A. J. (2005). A Generalized Finite Difference Method for Modeling Cardiac Electrical Activation on Arbitrary, Irregular Computational Meshes. Math. Biosci. 198, 169–189. doi:10.1016/j.mbs.2005.07.007
Tschöpe C., Ammirati E., Bozkurt B., Caforio A. L. P., Cooper L. T., Felix S. B., et al. (2021). Myocarditis and Inflammatory Cardiomyopathy: Current Evidence and Future Directions. Nat. Rev. Cardiol. 18, 169–193. doi:10.1038/s41569-020-00435-x
Tung L. (1978). A Bidomain Model for Describing Ischemic Myocardial D-C Potentials. Massachusetts: Ph.D. thesis.
Tveito A., Jæger K. H., Kuchta M., Mardal K.-A., Rognes M. E. (2017). A Cell-Based Framework for Numerical Modeling of Electrical Conduction in Cardiac Tissue. Front. Phys. 5, 48. doi:10.3389/fphy.2017.00048
Ugarte J. P., Orozco-Duque A., Tobón V., Novak D., Saiz J., Oesterlein T., et al. (2014). Dynamic Approximate Entropy Electroanatomic Maps Detect Rotors in a Simulated Atrial Fibrillation Model. Plos One 9, e114577–19. doi:10.1371/journal.pone.0114577
Unger L. A., Oesterlein T. G., Loewe A., Dössel O. (2019). “Noise Quantification and Noise Reduction for Unipolar and Bipolar Electrograms,” in 2019 Computing in Cardiology (CinC). doi:10.22489/cinc.2019.288
van der Does L. J. M. E., de Groot N. M. S. (2017). Inhomogeneity and Complexity in Defining Fractionated Electrograms. Heart rhythm. 14, 616–624. doi:10.1016/j.hrthm.2017.01.021
Verheule S., Schotten U. (2021). Electrophysiological Consequences of Cardiac Fibrosis. Cells 10, 3220. doi:10.3390/cells10113220
Verrier R. L., Fuller H., Justo F., Nearing B. D., Rajamani S., Belardinelli L. (2016). Unmasking Atrial Repolarization to Assess Alternans, Spatiotemporal Heterogeneity, and Susceptibility to Atrial Fibrillation. Heart rhythm. 13, 953–961. doi:10.1016/j.hrthm.2015.11.019
Vigmond E. J., Hughes M., Plank G., Leon L. J. (2003). Computational Tools for Modeling Electrical Activity in Cardiac Tissue. J. Electrocardiol. 36, 69–74. doi:10.1016/j.jelectrocard.2003.09.017
Vigmond E., Pashaei A., Amraoui S., Cochet H., Hassaguerre M. (2016). Percolation as a Mechanism to Explain Atrial Fractionated Electrograms and Reentry in a Fibrosis Model Based on Imaging Data. Heart rhythm. 13, 1536–1543. doi:10.1016/j.hrthm.2016.03.019
Keywords: cardiac modeling, fibrosis, electrogram, multiscale, microstructure
Citation: Sánchez J and Loewe A (2022) A Review of Healthy and Fibrotic Myocardium Microstructure Modeling and Corresponding Intracardiac Electrograms. Front. Physiol. 13:908069. doi: 10.3389/fphys.2022.908069
Received: 30 March 2022; Accepted: 21 April 2022;
Published: 10 May 2022.
Edited by:
Matthijs Cluitmans, Maastricht University, NetherlandsReviewed by:
Jason D. Bayer, Université de Bordeaux, FranceMarianna Meo, EPD Solutions, A Philips Company, Netherlands
Copyright © 2022 Sánchez and Loewe. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jorge Sánchez, publications@ibt.kit.edu