Abstract
Background: Atrial fibrillation (AF), the most common cardiac arrhythmia, is characterized by alteration of the action potential (AP) propagation. Under persistent AF, myocytes undergo electrophysiological and structural remodeling, which involves fibroblast proliferation and differentiation, modifying the substrate for AP propagation. The aim of this study was to analyze the effects on the AP of fibroblast-myocyte coupling during AF and its propagation in different regions of the atria.
Methods: Isolated myocytes were coupled to different numbers of fibroblasts using the established AP models and tissue simulations were performed by randomly distributing fibroblasts. Fibroblast formulations were updated to match recent experimental data. Major ion current conductances of the myocyte model were modified to simulate AP heterogeneity in four different atrial regions (right atrium posterior wall, crista terminalis, left atrium posterior wall, and pulmonary vein) according to experimental and computational studies.
Results: The results of the coupled myocyte-fibroblast simulations suggest that a more depolarized membrane potential and higher fibroblast membrane capacitance have a greater impact on AP duration and myocyte maximum depolarization velocity. The number of coupled fibroblasts and the stimulation frequency are determining factors in altering myocyte AP. Strand simulations show that conduction velocity tends to homogenize in all regions, while the left atrium is more likely to be affected by fibroblast and AP propagation block is more likely to occur. The pulmonary vein is the most affected region, even at low fibroblast densities. In 2D sheets with randomly placed fibroblasts, wavebreaks are observed in the low density (10%) central fibrotic zone and when fibroblast density increases (40%) propagation in the fibrotic region is practically blocked. At densities of 10 and 20% the width of the vulnerable window increases with respect to control but is decreased at 40%.
Conclusion: Myocyte-fibroblast coupling characteristics heterogeneously affect AP propagation and features in the different atrial zones, and myocytes from the left atria are more sensitive to fibroblast coupling.
Introduction
Atrial fibrillation (AF) is the most common cardiac arrhythmia and is becoming more prevalent with population aging (Iwasaki et al., 2011). In persistent atrial fibrillation (PeAF), lasting longer than 7 days (Kirchhof et al., 2016), spontaneous, pharmacological, or ablative resumption of sinus rhythm is infrequent with prompt recurrences or commonly failed cardioversions (Jalife and Kaur, 2015).
In PeAF, cardiac tissue experiences structural and electrical remodeling driven by changes of the extracellular matrix, fibroblast differentiation, fatty and inflammatory infiltration, and ion channel remodeling, among other alterations (Kirchhof et al., 2016). The heterogeneities in the atrial substrate cause fiber discontinuities and local conduction disturbances favoring re-entries (Nattel and Dobrev, 2017).
The atria have a complex anatomical structure, differences between the right and left atrium, and complex structures like the pulmonary vein or left and right appendages (Sanchez-Quintana et al., 2012; Ferrer et al., 2015). They all have their own electrophysiological characteristics, which determine specific action potential (AP) waveforms and conduction velocities (CV), which are crucial to AF pathophysiology.
Approximately 75% of the adult myocardium tissue volume is occupied by myocytes, accounting for only 30–40% of the total number of cells (Biernacka and Frangogiannis, 2011). Fibroblasts are considered as the predominant non-myocyte cell. They become active under inflammatory response and then differentiate into myofibroblasts (Chacar et al., 2017) by different pathways (e.g., oxidative stress, atrial dilatation, calcium overload, and inflammation), although the precise mechanisms have not yet been established (Jalife and Kaur, 2015). Myofibroblasts differ from fibroblasts in that they develop contractile proteins and exhibit a more depolarized resting membrane potential (Salvarani et al., 2017) and greater membrane capacitance (Sridhar et al., 2017). The differences in resting membrane potential and membrane capacitance between fibroblasts and myofibroblasts might have an important repercussion on the electrophysiology of coupled myocytes.
Recent experimental studies have shown that fibroblasts express voltage-dependent sodium channels allowing an inward current (INa) (Chatelier et al., 2012; Koivumäki et al., 2014a; Poulet et al., 2016). Although they are not electrically excitable they may affect the myocytes’ electrophysiological properties (Gaudesius et al., 2003; Kohl et al., 2005; Miragoli et al., 2006). Several animal experimental models have shown that fibroblasts electrically couple to myocytes and alter their AP (Rook et al., 1992; Gaudesius et al., 2003; Burstein et al., 2007; Rohr, 2011; Jalife and Kaur, 2015; Jousset et al., 2016; Quinn et al., 2016). Kohl et al. (1994) found that fibroblasts reacted electrically with myocytes through gap junctions and this was later corroborated in vitro in several animal species (Gaudesius et al., 2003; Camelliti et al., 2004; Kohl et al., 2005). However, the extent of (in vivo) fibroblast-myocyte electrotonic coupling in native myocardium remains controversial (Kohl and Gourdie, 2014). Recently, a study by Quinn et al. (2016) showed the existence of tunneling nanotube connections between excitable and non-excitable cells (e.g., myofibroblasts) in cardiac scar border tissue at the border zone of the scar in mice hearts. Similarly, Burstein et al. (2007) observed fibroblasts proliferation and differentiation into myofibroblasts during AF in vivo canine hearts. Fibroblast proliferation has also been observed in human hearts (Krul et al., 2015; Fukumoto et al., 2016), in regions with altered conduction velocity, although electrical coupling has not been demonstrated.
Computational simulations of the cellular electrophysiology have become a powerful tool for studying the role that fibroblasts play in the atrial substrate (MacCannell et al., 2007; Maleckar et al., 2009; Koivumäki et al., 2014a). Their flexibility permits modifying the cellular electrophysiology and varying different characteristics to analyze the impact of electrotonic coupling with electrically remodeled myocytes in PeAF (McDowell et al., 2013; Morgan et al., 2014; Sridhar et al., 2017). MacCannell et al. (2007) formulated an active and a passive model for fibroblasts electrophysiology, considering a capacitance and an ohmic resistance in the first case and a capacitance and four ionic currents in the second case. The active model is more realistic and has a major impact on myocytes, this is why in the present study the active model is implemented. Simulations also offer the possibility of studying how fibroblasts can change arrhythmia dynamics by reducing CV and introducing heterogeneities in the atrial substrate (Tanaka et al., 2007; Saha et al., 2018).
The present simulation study analyzes the effects of fibroblast-myocyte coupling in PeAF in specific structurally and electrophysiologically remodeled atrial tissues and compares the effect of fibroblasts and myofibroblasts on myocyte AP features, tissue conduction velocity and vulnerability to re-entries.
Materials and Methods
Myocyte Electrophysiological Models
The human atrial AP model proposed by Koivumäki et al. (2011) was used to simulate atrial APs. To reproduce the tissue electrophysiology of different anatomical atrial regions, we modified the conductance of five ionic currents: transient outward K+ current (Ito), potassium rapid current (IKr), potassium slow current (IKs), time independent K+ current (IK1), and L-type Ca2+ current (ICaL) as proposed in previous simulation studies (Krueger et al., 2012; Ferrer et al., 2015). The conductance values, shown in Table 1, were modified according to the experimental values to reproduce the AP waveform in all four regions: right atrium posterior wall (RA), crista terminalis (CT), left atrium posterior wall (LA), and the pulmonary vein (PV) areas.
Table 1
| Channel conductance | RA | LA | PV | CT |
|---|---|---|---|---|
| gto | 1.0 | 1.0 | 0.25 | 0.55 |
| gKr | 1.0 | 1.6 | 4.8 | 0.55 |
| gCaL | 1.0 | 0.67 | 0.3 | 1.0 |
| gK1 | 1.0 | 1.0 | 0.7425 | 1.0 |
| gKs | 1.0 | 1.0 | 5.12 | 1.0 |
Ionic conductance modifications.
Relative values for each atrial region (right atrium -RA-, left atrium -LA-, pulmonary veins -PV- and crista terminalis -CT-) with respect to the maximum conductance of the potassium transient outward current (gto), potassium rapid current (gKr), calcium L-type current (gCaL), and potassium slow current (gKs) in the Koivumäki et al. (2011) model.
Ionic Remodeling
Table 2 shows the values used to introduce AF electrical remodeling in all atrial regions by modifying ion channel conductances for ICaL, Ito, IK1, sustained outward K+ current (Isus), Na+/Ca2+ exchanger (NCX), sarcoplasmic reticulum Ca2+ ATPase (SERCA) pump, and ryanodine receptors (RyR), and specific calcium handling parameters, such as phospholamban (PLB), sarcolipin (SLN), and the baseline phosphorylation (phos). Dilation was also modeled by increasing the length of the cell by a factor of 1.1 (see Table 2) as described in Koivumäki et al. (2014b).
Table 2
| Lcell | gCaL | gto | gsus | gK1 | KNaCa | cpumps | PLB | SLN | phos | RyR |
|---|---|---|---|---|---|---|---|---|---|---|
| 1.1 | 0.41 | 0.38 | 0.62 | 1.62 | 1.5 | 0.84 | 1.18 | 0.6 | 2 | 2 |
Atrial fibrillation electrical remodeling.
Based on experimental observations/prior computational works (Koivumäki et al., 2014b), to account for persistent atrial fibrillation remodeling the ionic and cellular parameters modified were: cell length (Lcell), conductance of calcium L-type current (gCaL), conductance of potassium transient outward current (gto), conductance of potassium sustained current (gsus), conductance of potassium inward rectifier (gK1), maximum conductance of Na+/Ca2+ exchanger (KNaCa), conductance of sarcoplasmic reticulum Ca2+ ATPase (SERCA) pump (cpumps), phospholamban (PLB), sarcolipin (SLN), baseline phosphorylation (phos), and ryanodine receptors (RyR). The values indicate the factor by which the original values are multiplied.
Fibroblast Electrophysiological Model
Fibroblast electrophysiology was represented according to the Koivumäki et al. (2014a) model. Membrane sodium current was updated to match recent experimental results (Poulet et al., 2016) (see Supplementary Table S1). The uncoupled fibroblast resting membrane potential (RMPf) was set to -45 or -26 mV, as suggested by experimental data (Poulet et al., 2016; Salvarani et al., 2017). These values were obtained by shifting the gating variable for the time dependent potassium current as in previous simulation studies (Maleckar et al., 2009; Morgan et al., 2016), and shown in Supplementary Table S2. The fibroblast membrane capacitance (Cmf) values were varied within experimental ranges (6.3 and 50.4 pF) (MacCannell et al., 2007; Xie et al., 2009; McArthur et al., 2015; Poulet et al., 2016; Sridhar et al., 2017). Following previous experimental and simulation studies, we set RMPf to -45 mV and Cmf to 6.3 pF for atrial fibroblasts (MacCannell et al., 2007; Maleckar et al., 2009), and a more depolarized RMPf of -26 mV and a Cmf of 50.4 pF for atrial myofibroblasts (Morgan et al., 2016; Poulet et al., 2016; Salvarani et al., 2017; Sridhar et al., 2017). The code is available and can be found at http://hdl.handle.net/10251/120395.
Cellular Simulations (0D)
To compute the membrane potential in a cell, Equations (1, 2, and 3) were used:
where Cmyo is the membrane capacitance of the myocyte, Iion-myo is the total ionic current flowing through ionic channels of the myocyte, Vmyo is the membrane potential in the myocyte, Cfib is the membrane capacitance of the fibroblast, Iion-fib is the total ionic current flowing through ionic channels of the fibroblast, Vfib the membrane potential in the fibroblast, Igap is the total ionic current flowing through the gap junction, n is the number of coupled fibroblasts to one myocyte and Ggap is the coupling conductance between myocyte and fibroblast. Ggap was set to 0.5 nS (Morgan et al., 2016), a value within the range of 0.3–8 nS measured in cultured myocyte-fibroblast experiments as used in previous simulation studies (MacCannell et al., 2007; Rook et al., 1992). The number of fibroblasts coupled to a single myocyte varied from 0 to 9 (Koivumäki et al., 2014a), as shown in Supplementary Figure S1.
Stabilization of the cellular model of each region was achieved after 1000 stimuli of 2 ms duration and an amplitude of twice the threshold value. After this number of stimuli, APD90 values and ion concentration ranges reached steady state. The basic cycle length (BCL) was 1 s, 250 ms, or 2 s.
The following features of the myocyte AP were measured: AP duration at 90% repolarization (APD90), depolarization upstroke velocity measured as the maximum variation of the depolarization phase (dV/dtmax), and the RMP.
Strand and Tissue Simulations
The diffusion-reaction Eq. (4) of the monodomain formalism was used to simulate AP electrical propagation:
where D is the diffusion tensor, Cm is the membrane capacitance, and Iion stands for ionic currents through the membrane. Diffusion coefficients were adjusted for the different atrial regions to achieve physiological CVs along a strand of myocytes in this particular tissue.
Diffusion coefficients for myocytes (Dm) were calculated and adjusted to achieve physiological CVs, yielding a Dm of 3.84 cm2/s in RA and LA (CV = 70 cm/s), a Dm of 6.384 cm2/s in CT (CV = 100 cm/s), and a Dm of 4.081 cm2/s in PV (CV = 80 cm/s). In PeAF, the diffusion coefficient was reduced by 50% to reproduce gap junction remodeling (Harrild and Hen, 2000; Ten Tusscher et al., 2004; Krogh-madsen et al., 2012; McDowell et al., 2012).
For one-dimensional simulations (1D), a strand with a space discretization of 100 μm was used for myocytes (Satoh et al., 1996; Courtemanche et al., 1998; Nygren et al., 1998; Sachse et al., 2008) and 10 μm for fibroblasts-myofibroblasts (Sachse et al., 2008; Xie et al., 2009; Brown et al., 2015). Fibroblasts-myofibroblasts were randomly distributed along the strand at different densities (10, 20, and 40%) using a uniform random probabilistic function. One hundred random distributions were generated for each density value (see Supplementary Figure S3). The strand was paced with a BCL of 1 s for 50 s. The first element of the strand was stimulated with an amplitude of twice its threshold. CV was measured in the 50th pulse. The diffusion coefficient for elements with fibroblasts was halved with respect to the myocyte elements (Vasquez et al., 2004; Gomez et al., 2014b; Morgan et al., 2016).
Two-dimensional (2D) meshes representing cardiac tissues for RA and LA were built with a central fibrotic region of 2 cm diameter with different randomly distributed myofibroblast densities. Ten random distributions were implemented for each density (see Supplementary Figures S1, S3). The tissue grid had 524 × 524 elements with a spatial resolution of 100 μm. An anisotropic ratio of 2.86:1 was considered in RA and LA (Ferrer et al., 2015).
To study the effect of myofibroblasts on vulnerability to reentry, reentrant activity was generated using cross-field stimulation protocol S1-S2. S1 stimuli consisted of 10 pulses applied to the left border of the computational mesh to stabilize the tissue and S2 was applied to the bottom left corner of the tissue.
In order to obtain the instantaneous location of the rotors, we used phase singularity (PS) detections (Sabir et al., 1974; Gray et al., 1998; Gomez et al., 2014a; Martinez-Mateu et al., 2018) first building phase maps based on the Hilbert transform (HT) of the APs (Warren et al., 2003; Yamazaki et al., 2012) (5) by computing the instantaneous phase θ (6), whose values ranged from -π to π radians:
where V is the membrane potential. Then PSs, where all phases converge (7), were computed to track the rotor trajectory:
Phase singularity detections were also used to assess heterogeneity degree in the tissue due to the inclusion of fibrotic regions.
Implementation
The differential equations of the cellular models were solved using the Rush-Larsen method for the gating variables and the forward Euler method for all other ordinary differential equations with a time step of 10 μs. The monodomain Eq. (4) was solved using the finite elements method with “no flux” boundary condition.
Simulations were implemented in C++ and CUDA languages. Computations were performed on an Intel(R) Xeon(R) CPU E5-2603v3 processor with an NVIDIA Tesla K40c graphic card (see Supplementary Figure S2).
Results
Cellular Simulations
Figure 1 shows the APs obtained in cellular simulations with a BCL of 1000 ms under normal sinus rhythm (Panel A) and PeAF (Panel B) conditions, highlighting the differences in APD90 for each atrial region. Under physiological conditions, the APD90 was 236.4 ms for RA, 214 ms for LA, 176.1 ms for PV, and 292.8 ms for CT. However, under PeAF conditions APD90 was 165.5 ms, 139.8 ms, 120.8 ms, and 180.2 ms for the RA, LA, PV, and CT, respectively. APD90 decreased by 40% for RA, LA, and CT and by approximately 30% for PV with respect to normal conditions (see Panel C). In PaAF, in addition to APD shortening, RMP dropped from -75 to -79 mV, and dV/dtmax increased from 163 V/s to 168 V/s for RA, CT, and LA. RMP fell from -68 to -78 mV, and dV/dtmax rose from 157 V/s to 165 V/s for PV (see Panels E and D).
Figure 1
When fibroblasts were coupled to a single myocyte, we analyzed the effect of their electrophysiological characteristics (RMPf and Cmf) on myocyte electrical behavior (see Figure 2). When the number of coupled fibroblasts increased, the myocyte RMP became more depolarized and moved closer to the value of RMPf. A higher Cmf value seemed to further reduce myocyte APD when RMPf was -45 mV.
Figure 2
Figure 3 shows the effect of myocyte fibroblast coupling at different BCLs. As can be seen in Panel A, at a BCL of 300 ms myocyte AP was barely affected by coupling 1 fibroblast. APD90 was very slightly reduced at both RMPf values (purple and green traces in panel A), for all BCLs.
Figure 3
In panel B at higher Cmf the effect on APD90s was similar to Panel A, except for the fact that an RMPf of -26 mV, APD90 was slightly higher than for an RMPf of -45 mV, for all BCLs. As the number of fibroblasts increased (Panels C–F), the effects changed with the electrical characteristics of the fibroblasts and the BCL. Panels C and D show that for a different RMPf, APD90 changed with BCL. Longer APD90s were obtained for BCLs of 500 ms and 1000 ms for a RMPf of -26 mV. In Panel D, for a BCL of 300 ms, APD alternans arose, and APD90 increased with higher BCL at an RMPf of -26 mV. When the number of coupled fibroblasts increased (Panels E and F), APD90 was reduced to 50 ms for -45 mV and all BCLs and to less than 40 ms for an RMPf of -26 mV.
We also analyzed myocyte RMP and dV/dtmax (see Figure 4). Myocyte RMP (Panel B) changed similarly when Cmf increased for both RMPf values (-26 mV and -45 mV). Increasing the number of fibroblasts moved the myocyte RMP closer to RMPf. The most significant change was in the dV/dtmax (Panel A), which was much lower at an RMPf of -26 mV than the reduction obtained at an RMPf of -45 mV. When fibroblasts couple to myocytes their RMPf move closer to each other and more Na+ channels are available, although myocyte Na+ channel availability can be reduced due to a more depolarized RMP according to the number of coupled fibroblasts.
Figure 4
The experimental data agree with previous simulation studies in which fibroblasts presented an RMP of around -45 mV and a Cmf of 6.3 pF, while myofibroblasts have a more depolarized RMPf of -26 mV and a Cmf of 50.4 pF (Sridhar et al., 2017). In all the different atrial zones we analyzed the effects of coupling three fibroblasts or three myofibroblasts to a single myocyte (see Figure 5).
Figure 5
Panel A in Figure 5 shows how APD90 was reduced when fibroblasts were coupled for all BCLs. When myofibroblasts were coupled, for a BCL of 300 ms APD alternans arose and for a BCL of 1000 ms APD90 increased. Panel B shows APD90s in LA for different BCLs and similar fibroblast and myofibroblasts effects. The reduced APD90 in the PV case (Panel C) seems to be more pronounced in myofibroblast coupling. CT alternans also arose when myofibroblasts were coupled.
1D Strand Simulations
Conduction velocity in the atria changes locally according to the region’s characteristics, as shown in Figure 6. During nSR and without fibroblasts, the RA has a CV of 70 cm/s, CT has a CV of 100 cm/s, LA has a CV of 70 cm/s, and PV has a CV of 80 cm/s (blue discontinuous lines) (Ferrer et al., 2015). AP propagation along the strand was affected by PeAF electrical remodeling (red discontinuous lines) and CVs dropped significantly. At higher fibroblast (yellow) or myofibroblasts (green) densities in the strand, CV dropped. Boxplot measurements of the CV are represented for the 100 random distributions of fibroblasts for each density (10, 20, and 40%). The region with the greatest differential effect on CV was the PV, depending on whether the distribution was with fibroblasts or myofibroblasts. PV also experienced conduction blocks in some of the random distributions (indicated by yellow asterisks). Conduction block was also seen in LA at a density of 40% in some distributions.
Figure 6
2D Tissue Simulations
The 2D atrial tissue electrical activity in PeAF with different myofibroblast densities was analyzed to assess vulnerability to reentry. Figure 7 shows snapshots of phase maps (taken at the same time). Re-entrant circuits can be seen in the RA (top panels) and in the LA (bottom panels) in PeAF remodeling and increasing levels of myofibroblast density from left to right (membrane potential snapshots can be seen in Supplementary Figure S4). In the absence of myofibroblasts a functional reentry was obtained in RA and LA. The tip of these rotors (Figure 7, first column), corresponding to the PS superimposed in white on the phase maps, describes a regular circular path, which agrees with the results obtained in previous simulation studies (Wilhelms et al., 2013). In LA, the rotor tip describes a smaller path due to the shorter wavelength caused by shorter LA ERP (Fernandez-Lozano et al., 2006; Krueger et al., 2012; Ferrer et al., 2015).
Figure 7
When myofibroblasts were present in the center of the tissue (see section “Materials and Methods”), the obstacle altered the reentrant activity. Small percentages of myofibroblasts (10–20%) allowed the wavefront to propagate through the fibrotic region, but the electrophysiological heterogeneities of myofibroblasts and myocytes caused wavebreaks, which were detected as PSs (quantified in Supplementary Table S3). However, propagation in the fibrotic region was practically blocked when myofibroblast density was raised to 40%, which produced an anatomical reentry surrounding the fibrotic obstacle. Since the wavefront did not propagate through the fibrotic region the number of wavebreaks was significantly reduced, as were the number of PSs detected (see Supplementary Table S3).
The myofibroblasts in the tissue increased reentry vulnerability, measured as the vulnerable window (VW), a time interval for which premature S2 stimulation generates a reentry (Figure 8). In the RA the vulnerable window in the absence of myofibroblasts was 37 ms. When myofibroblasts density was raised to 10%, the VW increased to 38 ± 0.0 ms. VW also increased (39 ± 0.63 ms) when density was raised to 20%, but at 40% VW dropped below the control value (35 ± 0.82 ms). Interestingly, LA was more sensitive to myofibroblasts with a larger VW than the RA. The LA VW in the absence of myofibroblasts was 40 ms. When myofibroblast density was raised to 10%, VW rose to 40 ± 0.0 ms, at 20% it increased to 40.5 ± 0.53 ms and at 40% it dropped to 38 ± 0.88 ms.
Figure 8
Discussion
Computational modeling was used to investigate the effect of coupling fibroblasts and myofibroblasts to myocytes in four different regions of the atria during persistent AF. The study’s major findings can be summarized as follows: (i) myocyte-fibroblast coupling heterogeneously shortens myocyte APD and depolarizes the myocyte resting membrane potential in the 4 different atrial regions These effects strongly depend on the fibroblast electrophysiology. (ii) Fibroblasts, and specially myofibroblasts, introduce heterogeneities in the atrial substrate which alter the propagation of the AP, slowing conduction velocity. (iii) Fibroblasts-myofibroblasts change the atrial substrate during PeAF, which alters the vulnerability to re-entries and the vulnerable window presents a biphasic behavior related to myofibroblast density. These results suggest that the heterogeneity of the atrial tissue in the presence of fibroblasts/myofibroblats promotes reentrant events and alters the dynamics of arrhythmogenic propagation.
Heterogeneous Effects of PeAF Remodeling and Electrical Coupling of Fibroblasts in Atrial Tissues
Atrial substrate is differently affected by PeAF remodeling and by the presence of fibroblasts, due to the electrophysiological heterogeneity of the different atrial regions. Our results from isolated single cells show differences in the RMP, dV/dtmax, and APD for the four different atrial regions (RA, LA, CT, and PV) in nSR and in PeAF. These differences are in agreement with the simulations carried out by Krueger et al. (2013), who reported the differential effects of AF remodeling in the different atrial regions. It has to be noted that in contrast to the Krueger study, our model presents a long-term stability in all regions in single-cell and tissue simulations and also considers the effect of fibroblast coupling. To our knowledge, this is the first simulation study including the three components (atrial heterogeneity, AF remodeling, and fibroblasts) using a detailed electrophysiological AP model for fibroblasts and focusing on the analysis of the different effects exerted by fibrosis in the different atrial regions. A recent study by Roney et al. (2018) showed that high PS density in the PVs favored the effectiveness of PV isolation in ablation procedure. Their model also considered electrophysiological remodeling in AF, electrophysiological heterogeneities in different atrial regions, and fibrosis was simulated by changes in tissue conductivity. In a previous study (Roney et al., 2016), the same group modeled fibrosis by different methods and did not consider either electrophysiological heterogeneities in the atrial regions or AF remodeling to determine how different fibrosis models could affect rotor dynamics.
Different experimental studies show that atrial fibroblasts have a different electrophysiology from ventricular fibroblasts (MacCannell et al., 2007; Burstein and Nattel, 2008; Poulet et al., 2016; Salvarani et al., 2017). Morgan et al. (2016) found that fibroblast electrophysiology changes the dynamics of an arrhythmic process and provides relevant information on the effect of myocyte-fibroblast coupling in the atria. Our results indicate that RMPf, Cmf, and the number of coupled fibroblasts altered the behavior of myocytes AP, as was found in previous simulation studies (Maleckar et al., 2009; Koivumäki et al., 2014a; Sridhar et al., 2017). Furthermore, in the present study we found that introducing INa current into the fibroblast model had an interesting effect; due to the high RMP of isolated fibroblasts INa current was blocked but when fibroblasts were coupled with myocytes INa channels became available. Additionally, myocyte-fibroblast coupling led to a partial inactivation of the myocyte INa due to the higher RMP in the myocyte.
Fibroblast electrophysiology (RMPf and Cmf) changes myocyte AP characteristics (Jacquemet and Henriquez, 2008, 2009; Maleckar et al., 2009). Our simulation results also show that electrical coupling with myocytes increases atrial electrophysiological heterogeneity. Changes in the BCL altered the behavior of the coupled cells, with different responses in different regions. Interestingly, myofibroblast -myocyte coupling in regions with higher IK1 and ICaL (RA and CT) exhibited more sensitivity to changes in frequency, while regions with smaller IK1 and ICaL (PV) developed no AP for any of the BCLs. In contrast to McDowell et al. (2013), we defined different electric characteristics for the atrial myofibroblasts, which have a different effect on myocyte AP. Myofibroblasts act as the current source, raising the myocyte RMP (Jacquemet and Henriquez, 2007) according to the number of coupled myofibroblasts (Maleckar et al., 2009; Koivumäki et al., 2014a), thus leading to a partial inactivation of the myocyte INa current (see Supplementary Figure S5).
Heterogenous Effects of Fibrosis During PeAF in Atrial Tissues
Structural remodeling of the cardiac tissue contributes to reducing conduction velocity, delaying regional functional activations, and increasing structural heterogeneities, which are important factors for establishing a re-entrant driver or conduction block (Camelliti et al., 2005). Our results show that fibroblasts and myofibroblasts can alter the activation time in a 1D tissue strand, in agreement with different studies (Sachse et al., 2008; Xie et al., 2009). We implemented one hundred random configurations for different fibroblast/myofibroblast densities in the four atrial regions. Zhan et al. (2014) showed that fibroblasts can alter the CV and can lead to blocks in conduction with fibrosis densities of 40 and 45%, our results showed that a high density (40%) led to conduction blocks in the LA and that the PV was the most sensitive region to the presence of fibroblasts-myofibroblasts. Similarly, in an experimental study Miragoli et al. (2006) showed that myofibroblast proliferation changed the tissue conduction velocity and myocyte dV/dtmax. Our results showed a reduction in CV, in agreement with several other experimental and simulation studies that found that fibroblasts-myofibroblasts can establish an electric coupling with myocytes, reducing their dV/dtmax and activation time, reflected in reduced CV (Miragoli et al., 2006; Rohr, 2009; Xie et al., 2009; Yue et al., 2011). We also found a monotonic reduction in all four atrial regions.
Tanaka et al. (2007) have shown that local fibrosis distribution reduces CV in the different atrial regions, in agreement with our results, which showed a reduced CV with a tendency to homogenize in all four atrial regions. CV heterogeneity is responsible for giving the atria the characteristic activation times (Nguyen et al., 2012); if all the regions were to have a homogenous CV, this might induce the appropriate conditions for reentrant rhythms and conductions blocks (Gaspo et al., 1997). As the fast conduction systems’ (CT) conduction velocity was significantly reduced, this may be an interesting mechanism for AF in the right atrium.
Vulnerability to Reentry During PeAF in Fibrotic Tissue
Structural remodeling and endo-epicardial dissociation alter the atrial substrate and could produce macroreentries and focal activity (Everett and Olgin, 2007; Verheule et al., 2014). When, we analyzed the propagation in different regions of the atria and at different myofibroblast densities, we found that a low myofibroblasts percentage increased the number of PSs due to the wavebreaks. However, at higher percentages, propagation through the fibrotic zone was blocked, the number of wavebreaks and PSs decreased, and an anatomical reentry was anchored around the fibrotic zone, in agreement with previous studies (McDowell et al., 2013; Roney et al., 2016). Several simulation studies have also shown that reentry dynamics is altered by heterogeneities of the AP in the cardiac tissue (Colman et al., 2014; Varela et al., 2016), the presence of fibroblasts (Ashihara et al., 2012; Gomez et al., 2014a; Morgan et al., 2016), and that PSs increase in the zones with fibroblasts (Saha et al., 2018).
Waks and Josephson (2014) demonstrated that the rotation dynamics depends on the atrial tissue (RA or LA) and its electrophysiological characteristics, as did we in the present study, in which vulnerability to reentries and the dynamics of the rotation depended on the atrial region. LA presented slightly higher VWs, due to its shorter APD.
Gomez et al. (2014b) showed that the density of fibroblasts had a biphasic impact on the ventricular vulnerable window for reentry, while our results showed the same VW biphasic behavior for the first time in atrial tissue. Krul et al. (2015) found that local fibrosis is associated with reentrant activity, comparable to our results at low fibroblast density (10%), which can be considered as a region of local diffuse fibrosis, presented higher tissue vulnerability and resulted in multiple wavebreaks. When myofibroblast density was raised (20%), tissue vulnerability to reentry rose and conduction blocks occurred. However, at higher densities (40%) the conduction blocks also occurred but VW dropped, as was found by Campos et al. (2018). This suggests that myocyte-fibroblast coupling in PeAF plays an important role in AF (Morgan et al., 2014), with different effects in different atrial regions.
Limitations
Several limitations must be considered when drawing conclusions from the present study. Firstly, we did not consider the effect of electrical remodeling due to cytokines like TGF-B, which have been reported during structural remodeling and are known to affect myocyte AP (Burstein and Nattel, 2008; Nattel et al., 2008; Zahid et al., 2016).
Secondly, as fibroblast electrophysiology and differentiation into myofibroblasts still remain unclear, the models we used should be considered as an approximation (Nattel, 2018). Myofibroblast electrophysiological characteristics are still not well understood. Furthermore, it is important to highlight that although animal models and in vitro experiments prove the existence of electrotonic coupling between myocytes and fibroblasts (Kohl et al., 2005; Miragoli et al., 2006; Grand et al., 2014; Salvarani et al., 2017), this has not been reported in humans in vivo. However, human cardiac tissue with fibrosis presents altered electrical behavior (Krul et al., 2015; Fukumoto et al., 2016).
Thirdly, the pulmonary vein AP model was built on the basis of the reported experimental data for dogs and sheep APD90 and RMP. Electrophysiological data from human isolated myocytes from PV are still unavailable (Ehrlich et al., 2003).
Fourthly, while AF mechanisms are still unclear, there are different pathways by which reentrant drivers can be generated and maintained. Several studies have shown how non-myocyte cells like fibroblasts, macrophages (Hulsmans et al., 2017) and adipocytes (De Coster et al., 2018) alter the atrial substrate and promote arrhythmias. The effects of non-myocytes or endo-epicardial dissociation have not been considered in the present study, although these are factors which might alter the reentry dynamics and would be interesting to analyze in future studies (Verheule et al., 2014).
Fifthly, we are aware that fibrosis is a complex structure involving different actors such as collagen deposition, inflammatory cytokines and proteins which may alter the myocyte electrophysiology. There are different approaches to simulating the fibrotic regions such as non-conductive elements (Ten Tusscher et al., 2004), the paracrine effect changing the myocyte ion channel conductance (Zahid et al., 2016) and coupling to elements with a static RMP (Majumder et al., 2011) and using ionic models to describe the electrophysiology of fibroblasts (McDowell et al., 2013; Morgan et al., 2016; Saha et al., 2018).
Finally, we did not include the effect of different sizes or locations of fibrotic regions in the tissue. The fibrotic regions can in fact attract rotors (Roy et al., 2018), and future studies in this line would shed light into the mechanisms of chaotic rhythms when fibroblasts proliferate. Additionally, 3D simulations have shown how fibroblasts and non-conductive areas can modify the dynamics of the reentry (McDowell et al., 2013; Morgan et al., 2014). Zahid et al. (2016) introduced the paracrine effect, modifying myocyte electrophysiology in the fibrotic region in 3D simulations reducing the AP duration and the conductivity of the fibrotic region. We are also aware that 3D simulations include further structural details. The work by Ferrer et al. (2015), among others, presents a highly detailed atrial model. However, these simulations have a high computational cost while 2D simulations are computationally more efficient and provide a detailed insight into arrhythmia dynamics. Despite these limitations, we consider that our simulations represent a realistic heterogeneous PeAF remodeling (electrophysiological and structural) scenario in the different atrial regions and contribute a great deal to the understanding of AF mechanisms.
Conclusion
The results of the present simulation study show that fibroblast electrophysiology alters myocyte AP characteristics and leads to slower conduction velocity in atrial tissues affected by PeAF. These changes are heterogeneous within the different atrial regions.
Myocyte-fibroblast coupling creates a substrate in which the dynamics of arrhythmic reentries changes with the fibroblast density. By increasing the density of fibroblasts, reentries evolve from functional to anatomical around the obstacle formed by the fibrotic region. We also observed biphasic behavior of tissue vulnerability to reentries. Low myofibroblast densities (10 and 20%) increase the vulnerability to reentry, while a higher density (40%) reduces tissue vulnerability.
Statements
Author contributions
JoS, JaS, and BT contributed to the conception and design of the study. JG, LR, and LM-M contributed to the development of the model. JoS performed the computational simulations and JoS, JaS, and BT performed the analysis and interpretation of the results. All authors contributed to drafting of the manuscript and all revised, and approved the submitted version.
Funding
This work was supported by the “Plan Estatal de Investigación Científica y Técnica y de Innovación 2013–2016” from the Ministerio de Economía, Industria y Competitividad of Spain and Fondo Europeo de Desarrollo Regional (FEDER) DPI2016-75799-R (AEI/FEDER, UE), and by the Dirección General de Política Científica de la Generalitat Valenciana (PROMETEU 2016/088).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2019.00847/full#supplementary-material
References
1
AshiharaT.HaraguchiR.NakazawaK.NambaT.IkedaT.NakazawaY.et al (2012). The role of fibroblasts in complex fractionated electrograms during persistent/permanent atrial fibrillation: implications for electrogram-based catheter ablation.Circ. Res.110275–284. 10.1161/CIRCRESAHA.111.255026
2
BiernackaA.FrangogiannisN. G. (2011). Aging and cardiac fibrosis.Aging Dis.2158–173. 10.1016/j.bbi.2008.05.010
3
BrownT. R.Krogh-MadsenT.ChristiniD. J. (2015). Computational approaches to understanding the role of fibroblast-myocyte interactions in cardiac arrhythmogenesis.Biomed. Res. Int.2015:465714. 10.1155/2015/465714
4
BursteinB.NattelS. (2008). Atrial fibrosis: mechanisms and clinical relevance in atrial fibrillation.J. Am. Coll. Cardiol.51802–809. 10.1016/j.jacc.2007.09.064
5
BursteinB.QiX. Y.YehY. H.CalderoneA.NattelS. (2007). Atrial cardiomyocyte tachycardia alters cardiac fibroblast function: a novel consideration in atrial remodeling.Cardiovasc. Res.76442–452. 10.1016/j.cardiores.2007.07.013
6
CamellitiP.BorgT. K.KohlP. (2005). Structural and functional characterisation of cardiac fibroblasts.Cardiovasc. Res.6540–51. 10.1016/j.cardiores.2004.08.020
7
CamellitiP.DevlinG. P.MatthewsK. G.KohlP.GreenC. R. (2004). Spatially and temporally distinct expression of fibroblast connexins after sheep ventricular infarction.Cardiovasc. Res.62415–425. 10.1016/j.cardiores.2004.01.027
8
CamposF. O.ShiferawY.Weber dos SantosR.PlankG.BishopM. J. (2018). Microscopic isthmuses and fibrosis within the border zone of infarcted hearts promote calcium-mediated ectopy and conduction block.Front. Phys.6:57. 10.3389/fphy.2018.00057
9
ChacarS.FarèsN.BoisP.FaivreJ.-F. (2017). Basic signaling in cardiac fibroblasts.J. Cell. Physiol.232725–730. 10.1002/jcp.25624
10
ChatelierA.MercierA.TremblierB.ThériaultO.MoubarakM.BenamerN.et al (2012). A distinct de novo expression of Na v 1.5 sodium channels in human atrial fibroblasts differentiated into myofibroblasts.J. Physiol.5904307–4319. 10.1113/jphysiol.2012.233593
11
ColmanM. A.VarelaM.HancoxJ. C.ZhangH.AslanidiO. V. (2014). Evolution and pharmacological modulation of the arrhythmogenic wave dynamics in canine pulmonary vein model.Europace16416–423. 10.1093/europace/eut349
12
CourtemancheM.RamirezR. J.NattelS. (1998). Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model.Am. J. Physiol.275H301–H321. 10.1152/ajpheart.1998.275.1.H301
13
De CosterT.ClausP.KazbanovI. V.HaemersP.WillemsR.SipidoK. R.et al (2018). Arrhythmogenicity of fibro-fatty infiltrations.Sci. Rep.82050. 10.1038/s41598-018-20450-w
14
EhrlichJ. R.ChaT.-J.ZhangL.ChartierD.MelnykP.HohnloserS. H.et al (2003). Cellular electrophysiology of canine pulmonary vein cardiomyocytes: action potential and ionic current properties.J. Physiol.551801–813. 10.1113/jphysiol.2003.046417
15
EverettT. H.OlginJ. E. (2007). Atrial fibrosis and the mechanisms of atrial fibrillation.Hear. Rhythm4S24–S27. 10.1016/j.hrthm.2006.12.040
16
Fernandez-LozanoI.Toquero-RamosJ.Escudero-VelaC.Castedo-MejutoE.Escudier-VillaJ. M.Alonso-PulponL. (2006). [Left atrial posterior wall and pulmonary vein refractory periods are associated with atrial fibrillation inducibility in a swine model].Rev. Esp. Cardiol.59653–661.
17
FerrerA.SebastiánR.Sánchez-QuintanaD.RodríguezJ. F.GodoyE. J.MartínezL.et al (2015). Detailed anatomical and electrophysiological models of human atria and torso for the simulation of atrial activation.PLoS One10:e0141573. 10.1371/journal.pone.0141573
18
FukumotoK.HabibiM.IpekE. G.ZahidS.KhurramI. M.ZimmermanS. L.et al (2016). Association of left atrial local conduction velocity with late gadolinium enhancement on cardiac magnetic resonance in patients with atrial fibrillation.Circ. Arrhythm. Electrophysiol.9:e002897. 10.1161/CIRCEP.115.002897
19
GaspoR.BoschR. F.TalajicM.NattelS. (1997). Functional mechanisms underlying tachycardia-induced sustained atrial fibrillation in a chronic dog model.Circulation964027–4035. 10.1161/01.CIR.96.11.4027
20
GaudesiusG.MiragoliM.ThomasS. P.RohrS. (2003). Coupling of cardiac electrical activity over extended distances by fibroblasts of cardiac origin.Circ. Res.93421–428. 10.1161/01.RES.0000089258.40661.0C
21
GomezJ. F.CardonaK.MartinezL.SaizJ.TrenorB. (2014a). Electrophysiological and structural remodeling in heart failure modulate arrhythmogenesis. 2D simulation study.PLoS One10:e0117883. 10.1371/journal.pone.0103273
22
GomezJ. F.CardonaK.RomeroL.FerreroJ. M.TrenorB. (2014b). Electrophysiological and structural remodeling in heart failure modulate arrhythmogenesis. 1D simulation study.PLoS One9:e106602. 10.1371/journal.pone.0106602
23
GrandT.SalvaraniN.JoussetF.RohrS. (2014). Aggravation of cardiac myofibroblast arrhythmogeneicity by mechanical stress.Cardiovasc. Res.104489–500. 10.1093/cvr/cvu227
24
GrayR. A.PertsovA. M.JalifeJ. (1998). Spatial and temporal organization during cardiac fibrillation.Nature39275–78. 10.1038/32164
25
HarrildD. M.HenC. S. (2000). A computer model of normal conduction in the human atria.Circ. Res.87e25–e36. 10.1161/01.RES.87.7.e25
26
HulsmansM.ClaussS.XiaoL.AguirreA. D.KingK. R.HanleyA.et al (2017). Macrophages facilitate electrical conduction in the heart.Cell169510.e20–522.e20. 10.1016/j.cell.2017.03.050
27
IwasakiY. K.NishidaK.KatoT.NattelS. (2011). Atrial fibrillation pathophysiology: implications for management.Circulation1242264–2274. 10.1161/CIRCULATIONAHA.111.019893
28
JacquemetV.HenriquezC. S. (2007). Modelling cardiac fibroblasts: interactions with myocytes and their impact on impulse propagation.Europace 9(Suppl.6), vi29-vi37. 10.1093/europace/eum207
29
JacquemetV.HenriquezC. S. (2008). Loading effect of fibroblast-myocyte coupling on resting potential, impulse propagation, and repolarization: insights from a microstructure model.AJP Hear. Circ. Physiol.294H2040–H2052. 10.1152/ajpheart.01298.2007
30
JacquemetV.HenriquezC. S. (2009). Modulation of conduction velocity by nonmyocytes in the low coupling regime.IEEE Trans. Biomed. Eng.56893–896. 10.1109/TBME.2008.2006028
31
JalifeJ.KaurK. (2015). Atrial remodeling, fibrosis, and atrial fibrillation.Trends Cardiovasc. Med.25475–484. 10.1016/j.tcm.2014.12.015
32
JoussetF.MaguyA.RohrS.KuceraJ. P. (2016). Myofibroblasts electrotonically coupled to cardiomyocytes alter conduction: insights at the cellular level from a detailed in silico tissue structure model.Front. Physiol.7:496. 10.3389/fphys.2016.00496
33
KirchhofP.BenussiS.KotechaD.AhlssonA.AtarD.CasadeiB.et al (2016). 2016 ESC guidelines for the management of atrial fibrillation developed in collaboration with EACTS.Europace181609–1678. 10.1093/europace/euw295
34
KohlP.CamellitiP.BurtonF. L.SmithG. L. (2005). Electrical coupling of fibroblasts and myocytes: relevance for cardiac propagation.J. Electrocardiol.3845–50. 10.1016/j.jelectrocard.2005.06.096
35
KohlP.GourdieR. G. (2014). Fibroblast – myocyte electrotonic coupling: does it occur in native cardiac tissue?J. Mol. Cell. Cardiol.7037–46. 10.1016/j.yjmcc.2013.12.024
36
KohlP.KamkinA. G.KiselevaI. S.NobleD. (1994). Mechanosensitive fibroblasts in the sino-atrial node region of rat heart: interaction with cardiomyocytes and possible role.Exp. Physiol.79943–956. 10.1113/expphysiol.1994.sp003819
37
KoivumäkiJ.ClarkR. B.BelkeD.KondoC.FedakP.MaleckarM. M.et al (2014a). Na+ current expression in human atrial myofibroblasts: identity and functional roles.Front. Physiol.5:275. 10.3389/fphys.2014.00275
38
KoivumäkiJ. T.SeemannG.MaleckarM. M.TaviP. (2014b). In silico screening of the key cellular remodeling targets in chronic atrial fibrillation.PLoS Comput. Biol.10:e1003620. 10.1371/journal.pcbi.1003620
39
KoivumäkiJ. T.KorhonenT.TaviP. (2011). Impact of sarcoplasmic reticulum calcium release on calcium dynamics and action potential morphology in human atrial myocytes: a computational study.PLoS Comput. Biol.7:e1001067. 10.1371/journal.pcbi.1001067
40
Krogh-madsenT.AbbottG. W.ChristiniD. J. (2012). Effects of electrical and structural remodeling on atrial fibrillation maintenance?: a simulation study.PLoS Comput. Biol.8:e1002390. 10.1371/journal.pcbi.1002390
41
KruegerM.SeemannG.RhodeK.KellerD. U. J.SchillingC.ArujunaA.et al (2012). Personalization of atrial anatomy and elelectophysiology as a basis for clinical modeling of radio-frequency-ablation of atrial fibrillation.Med. Imaging IEEE Trans.321–12. 10.1109/TMI.2012.2201948
42
KruegerM. W.DornA.KellerD. U. J.HolmqvistF.CarlsonJ.PlatonovP. G.et al (2013). In-silico modeling of atrial repolarization in normal and atrial fibrillation remodeled state.Med. Biol. Eng. Comput.511105–1119. 10.1007/s11517-013-1090-1
43
KrulS. P. J.BergerW. R.SmitN. W.Van AmersfoorthS. C. M.DriessenA. H. G.Van BovenW. J.et al (2015). Atrial fibrosis and conduction slowing in the left atrial appendage of patients undergoing thoracoscopic surgical pulmonary vein isolation for atrial fibrillation.Circ. Arrhythm. Electrophysiol.8288–295. 10.1161/CIRCEP.114.001752
44
MacCannellK. A.BazzaziH.ChiltonL.ShibukawaY.ClarkR. B.GilesW. R. (2007). A mathematical model of electrotonic interactions between ventricular myocytes and fibroblasts.Biophys. J.924121–4132. 10.1529/biophysj.106.101410
45
MajumderR.NayakA. R.PanditR. (2011). Scroll-wave dynamics in human cardiac tissue: lessons from a mathematical model with inhomogeneities and fiber architecture.PLoS One6:e18052. 10.1371/journal.pone.0018052
46
MaleckarM. M.GreensteinJ. L.GilesW. R.TrayanovaN. A. (2009). Electrotonic coupling between human atrial myocytes and fibroblasts alters myocyte excitability and repolarization.Biophys. J.972179–2190. 10.1016/j.bpj.2009.07.054
47
Martinez-MateuL.RomeroL.Ferrer-AlberoA.SebastianR.Rodríguez MatasJ. F.JalifeJ.et al (2018). Factors affecting basket catheter detection of real and phantom rotors in the atria: a computational study.PLoS Comput. Biol.14:e1006017. 10.1371/journal.pcbi.1006017
48
McArthurL.ChiltonL.SmithG. L.NicklinS. A. (2015). Electrical consequences of cardiac myocyte: fibroblast coupling.Biochem. Soc. Trans.43513–518. 10.1042/BST20150035
49
McDowellK. S.VadakkumpadanF.BlakeR.BlauerJ.PlankG.MacLeodR. S.et al (2012). Methodology for patient-specific modeling of atrial fibrosis as a substrate for atrial fibrillation.J. Electrocardiol.45640–645. 10.1016/j.jelectrocard.2012.08.005
50
McDowellK. S.VadakkumpadanF.BlakeR.BlauerJ.PlankG.MacleodR. S.et al (2013). Mechanistic inquiry into the role of tissue remodeling in fibrotic lesions in human atrial fibrillation.Biophys. J.1042764–2773. 10.1016/j.bpj.2013.05.025
51
MiragoliM.GaudesiusG.RohrS. (2006). Electrotonic modulation of cardiac impulse conduction by myofibroblasts.Circ. Res.98801–810. 10.1161/01.RES.0000214537.44195.a3
52
MorganR.ColmanM.KrugerM.SeemannG.RhodeK.AslanidiO. (2014). “Evaluating effects of fibrosis in atrial arrhythmogenesis using 3D computational modelling,” in Proceedings of the international conference on Computing in Cardiology 2014, (Cambridge, MA), 765–768.
53
MorganR.ColmanM. A.ChubbH.SeemannG.AslanidiO. 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:474. 10.3389/fphys.2016.00474
54
NattelS. (2018). Molecular and cellular mechanisms of atrial fibrosis in atrial fibrillation.JACC Clin. Electrophysiol.3425–435. 10.1016/j.jacep.2017.03.002
55
NattelS.BursteinB.DobrevD. (2008). Atrial remodeling and atrial fibrillation: mechanisms and implications.Circ. Arrhythm. Electrophysiol.162–73. 10.1161/CIRCEP.107.754564
56
NattelS.DobrevD. (2017). Controversies about atrial fibrillation mechanisms: aiming for order in chaos and whether it matters.Circ. Res.1201396–1398. 10.1161/CIRCRESAHA.116.310489
57
NguyenT. P.XieY.GarfinkelA.QuZ.WeissJ. N. (2012). Arrhythmogenic consequences of myofibroblastmyocyte coupling.Cardiovasc. Res.93242–251. 10.1093/cvr/cvr292
58
NygrenA.FisetC.FirekL.ClarkJ. W.LindbladD. S.ClarkR. B.et al (1998). Mathematical model of an adult human atrial cell: the role of K+ currents in repolarization.Circ. Res.8263–81. 10.1161/01.RES.82.1.63
59
PouletC.KünzelS.BüttnerE.LindnerD.WestermannD.RavensU. (2016). Altered physiological functions and ion currents in atrial fibroblasts from patients with chronic atrial fibrillation.Physiol. Rep.4:e12681. 10.14814/phy2.12681
60
QuinnT. A.CamellitiP.Rog-ZielinskaE. A.SiedleckaU.PoggioliT.O’TooleE. T.et al (2016). Electrotonic coupling of excitable and nonexcitable cells in the heart revealed by optogenetics.Proc. Natl. Acad. Sci. U.S.A.11314852–14857. 10.1073/pnas.1611184114
61
RohrS. (2009). Myofibroblasts in diseased hearts: new players in cardiac arrhythmias?Hear. Rhythm6848–856. 10.1016/j.hrthm.2009.02.038
62
RohrS. (2011). Cardiac fibroblasts in cell culture systems: myofibroblasts all along?J. Cardiovasc. Pharmacol.57389–399. 10.1097/FJC.0b013e3182137e17
63
RoneyC. H.BayerJ. D.CochetH.MeoM.DuboisR.JaïsP.et al (2018). Variability in pulmonary vein electrophysiology and fibrosis determines arrhythmia susceptibility and dynamics.PLoS Comput. Biol.14:e1006166. 10.1371/journal.pcbi.1006166
64
RoneyC. H.BayerJ. D.ZahidS.MeoM.BoyleP. M. J.TrayanovaN. A.et al (2016). Modelling methodology of atrial fibrosis affects rotor dynamics and electrograms.Europace18(Suppl. 4), iv146–iv155. 10.1093/europace/euw365
65
RookM. B.Van GinnekenA. C. G.de JongeB.el AoumariA.GrosD.JongsmaH. J.et al (1992). Differences in gap junction channels between cardiac myocytes, fibroblasts, and heterologous pairs.Am. J. Physiol.263C959–C977. 10.1152/ajpcell.1992.263.5.C959
66
RoyA.VarelaM.AslanidiO. (2018). Image-based computational evaluation of the effects of atrial wall thickness and fibrosis on re-entrant drivers for atrial fibrillation.Front. Physiol.9:1352. 10.3389/fphys.2018.01352
67
SabirM. A.SosulskiF. W.FinlaysonA. J. (1974). Chlorogenic acid-protein interactions in sunflower.J. Agric. Food Chem.22575–578. 10.1021/jf60194a052
68
SachseF. B.MorenoA. P.AbildskovJ. A. (2008). Electrophysiological modeling of fibroblasts and their interaction with myocytes.Ann. Biomed. Eng.3641–56. 10.1007/s10439-007-9405-8
69
SahaM.RoneyC. H.BayerJ. D.MeoM.CochetH.DuboisR.et al (2018). Wavelength and fibrosis affect phase singularity locations during atrial fibrillation.Front. Physiol.9:1207. 10.3389/fphys.2018.01207
70
SalvaraniN.MaguyA.De SimoneS. A.MiragoliM.JoussetF.RohrS. (2017). TGF-β1(Transforming Growth Factor-β1) plays a pivotal role in cardiac myofibroblast arrhythmogenicity.Circ. Arrhythm. Electrophysiol.10:e004567. 10.1161/CIRCEP.116.004567
71
Sanchez-QuintanaD.Ramon Lopez-MínguezJ.PizarroG.MurilloM.Angel CabreraJ. (2012). Triggers and anatomical substrates in the genesis and perpetuation of atrial fibrillation.Curr. Cardiol. Rev.8310–326. 10.2174/157340312803760721
72
SatohH.DelbridgeL. M. D.BlatterL. A.BersD. M. (1996). Surface: volume relationship in cardiac myocytes studied with confocal microscopy and membrane capacitance measurements: species-dependence and developmental effects.Biophys. J.701494–1504. 10.1016/S0006-3495(96)79711-4
73
SridharS.VandersickelN.PanfilovA. V. (2017). Effect of myocyte-fibroblast coupling on the onset of pathological dynamics in a model of ventricular tissue.Sci. Rep.7:40985. 10.1038/srep40985
74
TanakaK.ZlochiverS.VikstromK. L.YamazakiM.MorenoJ.KlosM.et al (2007). Spatial distribution of fibrosis governs fibrillation wave dynamics in the posterior left atrium during heart failure.Circ. Res.101839–847. 10.1161/CIRCRESAHA.107.153858
75
Ten TusscherK.NobleD.NobleP. J.PanfilovA. V. (2004). A model for human ventricular tissue.Am. J. Physiol. Circ. Physiol.286H1573–H1589.
76
VarelaM.ColmanM. A.HancoxJ. C.AslanidiO. V. (2016). Atrial heterogeneity generates re-entrant substrate during atrial fibrillation and anti-arrhythmic drug action: mechanistic insights from canine atrial models.PLoS Comput. Biol.12:e1005245. 10.1371/journal.pcbi.1005245
77
VasquezC.MorenoA. P.BerbariE. J. (2004). Modeling fibroblast-mediated conduction in the ventricle.Comput. Cardiol.2004349–352. 10.1109/CIC.2004.1442944
78
VerheuleS.EcksteinJ.LinzD.MaesenB.BidarE.GharaviriA.et al (2014). Role of endo-epicardial dissociation of electrical activity and transmural conduction in the development of persistent atrial fibrillation.Prog. Biophys. Mol. Biol.115173–185. 10.1016/j.pbiomolbio.2014.07.007
79
WaksJ. W.JosephsonM. E. (2014). Mechanisms of atrial fibrillation – reentry, rotors and reality.Arrhythm. Electrophysiol. Rev.390–100. 10.15420/aer.2014.3.2.90
80
WarrenM.GuhaP. K.BerenfeldO.ZaitsevA.AnumonwoJ. M. B.DhamoonA. S.et al (2003). Blockade of the inward rectifying potassium current terminates ventricular fibrillation in the guinea pig heart.J. Cardiovasc. Electrophysiol.14621–631. 10.1046/j.1540-8167.2003.03006.x
81
WilhelmsM.HettmannH.MaleckarM. M.KoivumäkiJ. T.DösselO.SeemannG. (2013). Benchmarking electrophysiological models of human atrial myocytes.Front. Physiol.3:487. 10.3389/fphys.2012.00487
82
XieY.GarfinkelA.CamellitiP.KohlP.WeissJ. N.QuZ. (2009). Effects of fibroblast-myocyte coupling on cardiac conduction and vulnerability to reentry: a computational study.Hear. Rhythm61641–1649. 10.1016/j.hrthm.2009.08.003
83
YamazakiM.MironovS.TaravantC.BrecJ.VaqueroL. M.BandaruK.et al (2012). Heterogeneous atrial wall thickness and stretch promote scroll waves anchoring during atrial fibrillation.Cardiovasc. Res.9448–57. 10.1093/cvr/cvr357
84
YueL.XieJ.NattelS. (2011). Molecular determinants of cardiac fibroblast electrical function and therapeutic implications for atrial fibrillation.Cardiovasc. Res.89744–753. 10.1093/cvr/cvq329
85
ZahidS.CochetH.BoyleP. M.SchwarzE. L.WhyteK. N.VigmondE. J.et al (2016). Patient-derived models link re-entrant driver localization in atrial fibrillation to fibrosis spatial pattern.Cardiovasc. Res.110443–454. 10.1093/cvr/cvw073
86
ZhanH.XiaL.ShouG.ZangY.LiuF.CrozierS. (2014). Fibroblast proliferation alters cardiac excitation conduction and contraction: a computational study.J. Zhejiang Univ. Sci. B15225–242. 10.1631/jzus.B1300156
Summary
Keywords
atrial fibrillation, computer simulation, structural remodeling, myofibroblast, vulnerability
Citation
Sánchez J, Gomez JF, Martinez-Mateu L, Romero L, Saiz J and Trenor B (2019) 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
Received
16 December 2018
Accepted
19 June 2019
Published
04 July 2019
Volume
10 - 2019
Edited by
Natalia A. Trayanova, Johns Hopkins University, United States
Reviewed by
Sanjay Ram Kharche, University of Western Ontario, Canada; Craig Henriquez, Duke University, United States; Oleg Aslanidi, King’s College London, United Kingdom
Updates
Copyright
© 2019 Sánchez, Gomez, Martinez-Mateu, Romero, Saiz and Trenor.
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: Beatriz Trenor, btrenor@eln.upv.es
This article was submitted to Computational Physiology and Medicine, a section of the journal Frontiers in Physiology
Disclaimer
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.