Computational and Mathematical Methods in Medicine
Volume 2017 (2017), Article ID 9463010, 10 pages
https://doi.org/10.1155/2017/9463010

Research Article

Simulation of Atrial Fibrosis Using Coupled Myocyte-Fibroblast Cellular and Human Atrial Models

Yuan Gao, Yinglan Gong, and Ling Xia

Key Laboratory for Biomedical Engineering of Ministry of Education, Department of Biomedical Engineering, Zhejiang University, Hangzhou 310027, China

Correspondence should be addressed to Yinglan Gong; yinglangong@zju.edu.cn and Ling Xia; xialing@zju.edu.cn

Received 3 July 2017; Revised 10 November 2017; Accepted 3 December 2017; Published 26 December 2017

Academic Editor: Marta Parazzini

Copyright © 2017 Yuan Gao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Atrial fibrosis is characterized by expansion of extracellular matrix and increase in the number of fibroblasts which has been associated with the development and maintenance of atrial arrhythmias. However, the mechanisms how the fibrosis contributes to atrial arrhythmia remain incompletely understood. In this study, we used a proposed fibroblast model coupled with the human atrial myocyte to investigate the effects of fibrosis on atrial excitability and repolarization at both cellular and macroscopic levels. The 12-lead electrocardiogram (ECG) was also simulated to explore the index of clinical diagnosis for fibrosis. The simulation results showed that the fibrosis can modify action potential morphology of human atrial myocyte, slow down wave propagation, and have rate adaptation, thus causing the atrial electrical heterogeneity. The fibrosis alone was sufficient to cause arrhythmia, induce reentry wave, and result in low amplitude and wide P waves at normal heart rate and significant prolonged and inverse P waves at high heart rate. All these symptoms aggravated when the level of fibrosis increased. Our simulations demonstrated that fibrosis is the substrate of atrial arrhythmia and thereby may be a potential target in the treatment of atrial arrhythmias.

1. Introduction

Fibrosis is a critical aspect of cardiac dysfunction following myocardial infarction, hypertension, heart failure, and severe arrhythmia [1, 2]. It is characterized by expansion of myocardial extracellular matrix and increase in the number of cardiac fibroblasts (Fbs). It increases with age and can lead to slow propagation and contribute to unidirectional block [3]. Apart from electrical remodeling and contractile dysfunction, atrial fibrosis has been shown to increase the susceptibility to atrial fibrillation (AF) and may serve as a critical substrate in the formation of the arrhythmia [4, 5]. A study about the small mouse atria has indicated that the alterations in atrial conduction produced by atrial interstitial fibrosis alone were sufficient to produce a substrate for AF [6].

But the mechanisms by how the fibrosis contributes to atrial arrhythmia remain incompletely understood. Cardiac fibroblast was the essential cell type in heart that is responsible for the homeostasis of the extracellular matrix. But it may transform to a myofibroblast phenotype and contribute to cardiac fibrosis when under impairment [7]. Many studies have revealed that myocytes and fibroblasts are functionally coupled by gap junctions when fibrosis occurred and may interact and modify impulse conduction via two-way electronic signal [8–10]. These results suggest a potential role for myocyte-fibroblast coupling in atrial arrhythmia.

Although there has been significant progress in the study of fibrosis over the past few decades, there still remain many unanswered questions. Investigation on human atrial myocardium with fibrosis is still very limited due to the difficulties in accessing experimental data. Therefore, alternative methods such as computer simulations are of great importance [11]. The main purpose of this study was to investigate the effects of fibrosis on atrial excitability and repolarization. To achieve this goal, we used a proposed fibroblast model [12] coupled with the Courtemanche et al. model [13] for atrial myocyte to stimulate the changes of action potential (AP) at the cellular level. At the macroscopic level, the AP models were incorporated into a human atrial model to simulate cardiac excitation conduction. Since the 12-lead ECG has been widely accepted as the main noninvasive diagnostic method of cardiac disease, we simulated the fibrosis ECG patterns to serve as an assistant method to study the underlying relationship between morphology of ECG and fibrosis atrial electrical activity.

2. Material and Methods

2.1. Myocyte-Fibroblast Electrophysiological Coupling

In our simulation, the atrial myocyte was surrounded by a selected number of fibroblasts in normal tissue, which was regarded as one myocyte/Fbs complex unit. The paradigm is illustrated in Figure 1.

(a)

(b)

Figure 1: Illustration of myocyte/Fbs complex unit. The rectangle indicates a normal atrial myocyte, the ellipse indicates a myocyte/Fbs complex unit, and the circle around the myocyte indicates a fibroblast. (a) View of one myocyte/Fbs complex unit: a single human atrial myocyte is coupled a selected number of homogeneous fibroblasts via assigning an intercellular conductance (gap) [14]. In this study, the number of fibroblasts was 1 to 6. (b) Schematic representation of 2D atrial tissue: myocyte/Fbs complex unit that replaces the normal myocyte random existing in the tissue indicates the fibrosis proliferous area [15].

The electrophysiological model of atrial myocyte is based on the Courtemanche et al. [13], which can accurately represent the atrial AP dynamics. Using specific formulations of Na⁺, K⁺, and Ca²⁺ currents based on the experimental data recorded from human atrial myocytes, along with representations of pump, exchange, and background currents, the model is computationally efficient and correctly replicates physiologically repolarization process.

The atrial fibroblast model was based on the active 1 model described by Maleckar et al. [12]. The model contains four active membrane ionic currents: time and voltage dependent K⁺ current, inward-rectifying K⁺ current, Na⁺-K⁺ pump current, and background Na⁺ current. This model originates from a previous mathematical formulation of MacCannell et al. [14] and is modified based on the experimental data to well represent human atrial fibroblast properties. The membrane potential of the coupling model is governed by the following equation:

where and represent the membrane potential of the atrial myocyte and fibroblast, respectively, and are the membrane capacitance of the myocyte and fibroblast, respectively, is the net membrane current of the myocyte, and is the net membrane current of fibroblast. is the stimulus current applied to myocyte membrane, is the current that flows through the gap junction between the myocyte and each fibroblast, is the total number of fibroblasts, and represents the gap junction conductance. According to the experimental reports, ranges from 6.3 to 75 pF [16]; between a myocyte and a fibroblast ranges from 0.3 to 8 nS in cultured cells [17]. In this study, we used = 100 pF, = 6.3 pF, and = 3 nS.

2.2. Model of 3D Atria

The atrial specimen was obtained from a healthy male adult in Zhujiang Hospital, Southern Medical University, China. The use of the heart for research purpose was approved by the local Ethics Committee of the Southern Medical University. The National Rules and Regulations on Heart research were strictly followed. The specimen was scanned by spiral computer tomography (Philips/Brilliance 64) with a resolution of 512 pixels by 512 pixels and a spatial resolution of 0.3574 × 0.3574 × 0.33 mm. Details of the model were described in our previous study [18–20].

In the model, the conduction system included sinoatrial node (SAN), Bachmann’s bundle (BB), crista terminalis (CT), pectinate muscles (PM), slow pathway (SP), and fast pathway (FP). In order to simulate the anisotropy, the fiber orientation of atria was contained. During the propagation, each myocardial unit has specific electrophysiological parameters associated with the action potential of the cell unit and conduction velocity.

The propagation of action potential was based on the monodomain model [21]:

where is transmembrane voltage, is the membrane capacitance, is surface-to-volume ratio, and we used = 1000 cm⁻¹ and = 1 μF/cm² [22]. is the ratio of conductivity extracellular to intracellular, is cellular conductivity, and conductivity values were based on experimental data studied by Roth [23]. The transversal conductivity was set to 0.15 S/m for myocyte and CT, 0.11 S/m for PM and BB, and 0.05 S/m for SAN. The ratio for cross-axis to long-axis conductivity was set to 1 : 9. was set to 1.0. is the sum of ionic currents, and is the sum of applied stimulus currents.

The equation was solved numerically using explicit Euler method based on parallel computational techniques. The simulations were performed on a cluster of networked Dawning TC4000L system. The hardware architecture is symmetric multiprocessor shared memory that contains one management node and ten computation nodes. Each computation node consisted of two Intel Xeon 5335 processors and 4 GB of shared memory. The total theoretical computing capacity can be up to 184 Gflops. We used MPI to implement the communication between nodes.

The torso model in our study was taken from the virtual male subject of the United States. The body surface potentials generated by the cardiac sources satisfy the Poisson equation with Newman boundary conditions:

where is the tissue dependent conductivity tensor, is the quasi static potential, is the density of the equivalent dipole sources, is the normal vector, is the body surface, which encloses the volume conductor .

Using the Green second identity,

with

The differential equation for as (3) can be solved as the following integral equation:

where = is the vector between the field point and source point ; is the heart area, is the conductivity junction surface, and its inside and outside conductivities are ⁻ and ⁺, respectively. Further details of the model can be found in our previous studies [24, 25].

The 12-leads ECGs are calculated as [26]

where is left arm surface potential, is right arm surface potential, is left leg surface potential, is each precordial lead , and is each precordial surface potential.

3. Results

3.1. Action Potential of Atrial Myocyte Coupling Fibroblasts

In our study, we first simulated the APs of one human atrial myocyte coupled different number of fibroblasts with a of 3 nS at stimulation frequencies between 1 and 3 Hz. More than 5 to 10 stimuli were applied to achieve steady state conditions.

Figure 2 shows the APs of atrial myocyte when coupled with 1, 4, and 6 fibroblasts, respectively, in comparison with the noncoupled control for pacing at 1 Hz. In Table 1, we compare typical characteristics of action potentials for the coupled myocyte. By comparison, we can see that, with more coupled fibroblasts, the maximum overshoot potential is decreased and resting potential is increased. When the number of fibroblasts increased, the membrane potential during the plateau was less depolarized, and the AP duration to 60% repolarization (APD₆₀) decreased quickly. But the 90% repolarization (APD₉₀) was prolonged with the increasing of fibroblasts.

Table 1: Characteristics for an atrial myocyte coupled Fbs with a of 3 nS for 1 Hz pacing.

#FB
(mV)
(mV)
APD60 (ms)
APD90 (ms)

0
−81.2
24.4
230.8
315.7

1
−79.6
23.2
210.6
320.6

4
−74.1
18.1
114.0
349.1

6
−69.2
6.5
71.6
412.6

#FB indicates the number of fibroblasts coupled myocyte, indicates the resting potential, and is the maximum overshoot value.

Figure 2: Illustration of changes in the waveform of atrial myocyte action potential when coupling to Fbs with a of 3 nS for 1 Hz pacing.

Figure 3 shows the APD₉₀ of atrial myocyte coupled fibroblasts at stimulation frequencies between 1 and 3 Hz. Each of the cells displayed marked frequency dependent adaptation of the APD with a shortening at higher stimulation rates. The APD of noncoupled myocyte changed smoothly at low rates and sharply at high rates. However, the APD of coupled myocyte changed linearly with the rate increased. Remarkably, when coupled with 6 fibroblasts, the myocyte failed to achieve steady state condition at high frequency.

Figure 3: APD₉₀ of one atrial myocyte coupled Fbs at different stimulation frequencies.

3.2. Simulation of 3D Atrial Exciting Sequence Map and 12-Lead ECG

Although the ratio has not been conclusive, fibrosis is characterized by the increase in the number of fibroblasts [27]. In this study, we simulated one atrial myocyte coupled 4 fibroblasts as a myocyte/Fbs complex unit. Then such a unit was used to randomly replace one normal myocyte in atria to simulate fibrosis. 10% and 70% of total atrial myocytes were replaced by the units to simulate gentle and severe fibrosis.

Figure 4 shows the exciting sequence map of atria with 10% and 70% fibrosis at 1 Hz pacing (i.e., the sinus rhythm was 60 beats per minute). For comparison, the exciting sequence map of normal atria was also presented. The first column of Figure 4 shows the typical normal atrial exciting sequence maps. The depolarization duration of the right atrium was 86 ms. The last part to be activated was at the site below the inferior vena cava. The total depolarization duration of the atria was 102 ms at the place of posterior left atrial wall. The second column shows the exciting sequence maps of the atria with 10% fibrosis. The right atrial depolarization duration was 90 ms and the total depolarization duration was 106 ms. The last activated areas were nearly the same as the normal case. The third column was the exciting sequence maps of the atria with 70% fibrosis. The last place to be activated in the right atrium was the right lateral wall at the time of 104 ms. The total depolarization duration of the atria was 133 ms nearly at the same place as the normal case.

Figure 4: Simulated activation sequences at 1 Hz pacing. First column: normal atria; second column: atria with 10% fibrosis; third column: atria with 70% fibrosis. The first row is the anterior view of the atria; the second row is the posterior view of atria. The arrows indicate the direction of the wave propagation. The color bar on the right-hand side indicates the propagation time with the units in milliseconds.

Figure 5 shows the 12-lead ECG of atrium with 10% and 70% fibrosis in comparison with normal atria at 1 Hz pacing in two cardiac cycles. Because the fibrosis decreased the cellular maximum overshoot potential and reduced the atrial electrical excitation, this led to the smaller amplitude of P wave on each lead. As the fibrosis prolonged the atrial excitation time, the P wave duration was extended. The phenomenon has an increasing trend with the fibrosis degree growing. Since the time of atrial repolarization was sufficient, ECG patterns in the second cycle did not significantly change.

Figure 5: 12-lead ECG at 1 Hz in two cycles. The black lines are the normal case, red lines are atria with 10% fibrosis, and blue lines are the atria with 70% fibrosis.

Figure 6 shows the exciting sequence map of atria with 10% and 70% fibrosis at 3 Hz pacing in three cardiac cycles. The first cycle of the excitation patterns was the same as Figure 4. The conduction pattern of atria with 10% fibrosis within the second cycle was slightly different. Since the repolarization of atrial myocytes near the proximal ends of conductive bundles was not sufficient, the excitation potential was far below the normal level. This led to the slow excitation propagation in the atria and slightly changed the exciting sequence. So the early part of P wave was gentle and an obvious wide P wave was visible on each lead (see Figure 7). The conduction pattern of atria with 10% fibrosis in the third cycle was more unusual. A quick exciting wave was initiated at the right posterior wall propagating around inferior vena cava from right atrium to the left atrium and merged with the planner wave propagated from the left atrial roof. It also propagated retrograde to the right lateral wall and merged with the wave propagated from the right anterior wall. This led to the notches of P waves appearing on most leads and inverse P wave on lead I.

Figure 6: Simulated activation sequences at 3 Hz pacing. First column: normal atria; second column: atria with 10% fibrosis; third column: atria with 70% fibrosis. First two rows: anterior and posterior view of first cardiac cycle. Since the total activation duration is less than 133 ms, time window of the beginning 140 ms was given. Middle two rows: anterior and posterior view of second cardiac cycle. Last two rows: anterior and posterior view of third cardiac cycle. The arrows indicate the direction of the wave propagation. The color bars on the right-hand side indicate the propagation time with the units in milliseconds.

Figure 7: 12-lead ECG at 3 Hz in three cycles. The black lines are the normal case, red lines are atria with 10% fibrosis, and blue lines are the atria with 70% fibrosis.

As the fibrosis level increased to 70%, the myocytes near the sinoatrial node failed to be activated in the second cycle. The first exciting area was the right posterior wall below the inferior vena cava, initiating the reentry wave meandered in both the left atrium and the right atrium. The total activating time persisted nearly the whole cardiac cycle and merged with the next cycle. So the prolonged and inverse P waves were more significant in the ECG patterns.

4. Discussions

In this study, based on a coupled myocyte-fibroblast model, we investigated the effects of fibrosis on human atria at both cellular and macroscopic level. At the cellular level, the AP morphology of the atrial myocyte is modified remarkably during depolarization and repolarization. With the number of coupled fibroblasts increasing, the resting membrane potential elevated and leads to weakened excitability along with the maximum overshoot potential decreased. Curtailment of 60% repolarization was seen and plateau of the AP disappeared depending on the increased number of coupled fibroblasts. However, the APD₉₀ was prolonged with the increasing of fibroblasts. The atrial myocyte also displayed marked frequency dependent on the APD with a shortening at higher stimulation rate.

At the macroscopic level, the fibrosis changed the exciting wave pattern in the atria. When the sinus rhythm was normal, the velocity of wave propagation slowed down with the increased level of fibrosis, leading to the total duration of atrial depolarization and repolarization being prolonged. This resulted in the smaller amplitude and wider P wave on each lead. Since the time of atrial repolarization was sufficient, the excitation and ECG patterns at each cardiac cycle did not have obvious change. With the heart rate increasing, the electrical anomalous conduction in the atria was more obvious. Local electrical conduction block and reentry waves were initiated and meandered in both the right atrium and the left atrium. Notches of P waves appeared on most leads; significant prolonged and inverse P waves were observed. The reduced amplitude and prolonged duration of P waves have been observed in previous ECG recordings of animal and human with atrial fibrosis [28–31]. But there was no report about inverse P wave in patient with atrial fibrosis. However, this morphology has been shown to be common in patients with AF [32–35]. Notably, other morphologies of P waves in our simulation (e.g., wide duration and notches) also can be seen in these patients. This suggested the possibility that patients with atrial fibrosis can show similar ECG patterns characterized by AF.

Since atrial fibrosis is associated with a variety of cardiomyopathies and one of the main factors determining the relapse after therapies, the evaluating of existence and level of atrial fibrosis in atria has increasing importance for therapy strategies [5]. Nowadays late gadolinium-enhanced MRI (LGE-MRI) is the mainstream technique for detecting fibrosis [36]. But the image quality of LGE-MRI scans is frequently poor due to various reasons (e.g., residual respiratory motion, heart rate variability, low resolution limited by the thinness of atrial wall, and confounded enhancement from surrounding heart substructures). These caused a large number of false positives in the atrial fibrosis delineation. ECG may serve as a clinical method to guide the LGE-MRI evaluation. This could have implications for treatment as well as future clinical trials. Our simulations provide a promising starting point for assessing patients with atrial fibrosis by ECG, and this could provide some references for clinical diagnosis.

It should be point out that there are still some limitations in this study. In our simulation, the atrial fibrosis model was constructed by uniformly random distributed fibrosis but not patient-derived. At the macroscopic level, we only considered the percentage increasing to represent the enhancement of fibrotic level, and the different number of coupled fibroblasts was not taken into account. Finally, the model used in this study was a static heart model with electrophysiological properties, but the mechanical functions have not been considered. Cardiac motion should be considered in future studies to further improve the simulation accuracy.

5. Conclusions

In this study, a coupled myocyte-fibroblast model has been used to investigate the effects of atrial fibrosis at cellular and human atrial levels. The 12-lead ECG was also simulated to explore the index of clinical diagnosis of fibrosis. The results show that the fibrosis can modify AP morphology of human atrial myocyte, slow down wave propagation, and have rate adaptation. These caused the atrial electrical heterogeneity. The fibrosis alone was sufficient to cause arrhythmia, induce reentry wave, and result in low amplitude and wide P waves at normal heart rate and significant prolonged and inverse P waves at high heart rate. These symptoms will aggravate with the level of fibrosis increased. Our simulations demonstrated that fibrosis is the substrate of atrial arrhythmia and thereby may be a potential target in treatment of atrial arrhythmia. The understanding of the ECG changes could provide some references for clinical diagnosis.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National R&D Program for Major Research Instruments (Grant no. 61527811), the National Natural Science Foundation of China (Grant no. 61701435), and the Zhejiang Provincial Natural Science Foundation of China (LY17H180003).

References

A. Piek, R. A. de Boer, and H. H. W. Silljé, “The fibrosis-cell death axis in heart failure,” Heart Failure Reviews, vol. 21, no. 2, pp. 199–211, 2016.

S. D. Prabhu and N. G. Frangogiannis, “The biological basis for cardiac repair after myocardial infarction: from inflammation to fibrosis,” Circulation Research, vol. 119, no. 1, pp. 91–112, 2016.

B. J. Hansen, J. Zhao, T. A. Csepe et al., “Atrial fibrillation driven by micro-anatomic intramural re-entry revealed by simultaneous sub-epicardial and sub-endocardial optical mapping in explanted human hearts,” European Heart Journal, vol. 36, no. 35, pp. 2390–2401, 2015.

H. Cochet, A. Mouries, H. Nivet et al., “Age, atrial fibrillation, and structural heart disease are the main determinants of left atrial fibrosis detected by delayed-enhanced magnetic resonance imaging in a general cardiology population,” Journal of Cardiovascular Electrophysiology, vol. 26, no. 5, pp. 484–492, 2015.

N. F. Marrouche, D. Wilber, G. Hindricks et al., “Association of atrial tissue fibrosis identified by delayed enhancement MRI and atrial fibrillation catheter ablation: the DECAAF Study,” Journal of the American Medical Association, vol. 311, no. 5, pp. 498–506, 2014.

S. Verheule, T. Sato, T. Everett et al., “Increased vulnerability to atrial fibrillation in transgenic mice with selective atrial fibrosis caused by overexpression of TGF-β1,” Circulation Research, vol. 94, no. 11, pp. 1458–1465, 2004.

V. Nagpal, R. Rai, A. T. Place et al., “Response to letter regarding article, “MiR-125b is critical for fibroblast-to-myofibroblast transition and cardiac fibrosis”,” Circulation, vol. 133, no. 24, pp. e714–e714, 2016.

J. Pellman, J. Zhang, and F. Sheikh, “Myocyte-fibroblast communication in cardiac fibrosis and arrhythmias: mechanisms and model systems,” Journal of Molecular and Cellular Cardiology, vol. 94, pp. 22–31, 2016.

S. L. K. Bowers and T. A. Baudino, “Cardiac myocyte-fibroblast interactions and the coronary vasculature,” Journal of Cardiovascular Translational Research, vol. 5, no. 6, pp. 783–793, 2012.

M. Miragoli, G. Gaudesius, and S. Rohr, “Electrotonic modulation of cardiac impulse conduction by myofibroblasts,” Circulation Research, vol. 98, no. 6, pp. 801–810, 2006.

S. Zahid, H. Cochet, P. M. Boyle et al., “Patient-derived models link re-entrant driver localization in atrial fibrillation to fibrosis spatial pattern,” Cardiovascular Research, vol. 110, no. 3, pp. 443–454, 2016.

M. M. Maleckar, J. L. Greenstein, W. R. Giles, and N. A. Trayanova, “Electrotonic coupling between human atrial myocytes and fibroblasts alters myocyte excitability and repolarization,” Biophysical Journal, vol. 97, no. 8, pp. 2179–2190, 2009.

M. Courtemanche, R. J. Ramirez, and S. Nattel, “Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 275, no. 1, pp. 301–321, 1998.

K. A. MacCannell, H. Bazzazi, L. Chilton, Y. Shibukawa, R. B. Clark, and W. R. Giles, “A mathematical model of electrotonic interactions between ventricular myocytes and fibroblasts,” Biophysical Journal, vol. 92, no. 11, pp. 4121–4132, 2007.

H. Zhan and L. Xia, “Excitation-contraction coupling between human atrial myocytes with fibroblasts and stretch activated channel current: a simulation study,” Computational and Mathematical Methods in Medicine, vol. 2013, no. 3, Article ID 238676, 9 pages, 2013.

C. Vasquez, A. P. Moreno, and E. J. Berbari, “Modeling fibroblast-mediated conduction in the ventricle,” Computers in Cardiology, pp. 349–352, 2004.

M. B. Rook, A. C. G. Van Ginneken, B. de Jonge, A. El Aoumari, D. Gros, and H. J. Jongsma, “Differences in gap junction channels between cardiac myocytes, fibroblasts, and heterologous pairs,” American Journal of Physiology-Cell Physiology, vol. 263, no. 5, pp. C959–C977, 1992.

L. Xia, Y. Gao, Q. Lu, D. Zheng, D. D. Deng, and L. Xia, “Preliminary simulation study of atrial fibrillation treatment procedure based on a detailed human atrial model,” Journal of Clinical Trials in Cardiology, vol. 2, no. 4, pp. 1–9, 2015.

D.-D. Deng, Y.-L. Gong, G.-F. Shou et al., “Simulation of biatrial conduction via different pathways during sinus rhythm with a detailed human atrial model,” Journal of Zhejiang University Science B, vol. 13, no. 9, pp. 676–694, 2012.

D. Deng, P. Jiao, X. Ye, and L. Xia, “An image-based model of the whole human heart with detailed anatomical structure and fiber orientation,” Computational and Mathematical Methods in Medicine, vol. 2012, no. 3, Article ID 891070, pp. 449–461, 2012.

Y. Zhang, L. Xia, Y. Gong, L. Chen, G. Hou, and M. Tang, “Parallel solution in simulation of cardiac excitation anisotropic propagation,” in International Conference on Functional Imaging and Modeling of the Heart, Lecture Notes in Computer Science, pp. 170–179, 2007.

M. Potse, B. Dubé, A. Vinet, and R. Cardinal, “A comparison of monodomain and bidomain propagation models for the human heart,” in Proceedings of the 28th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, (EMBS '06), pp. 3895–3898, New York, NY, USA, September 2006.

B. J. Roth, “Electrical conductivity values used with the bidomain model of cardiac tissue,” IEEE Transactions on Biomedical Engineering, vol. 44, no. 4, pp. 326–328, 1997.

L. Xia, M. Huo, Q. Wei, F. Liu, and S. Crozier, “Electrodynamic heart model construction and ECG simulation,” Methods of Information in Medicine, vol. 45, no. 5, pp. 564–573, 2006.

G. Shou, L. Xia, M. Jiang, F. Liu, and S. Crozier, “Forward and inverse solutions of electrocardiography problem using an adaptive BEM method,” International Conference on Functional Imaging and Modeling of the Heart, pp. 290–299, 2007.

P. Kligfield, L. S. Gettes, J. J. Bailey et al., “Recommendations for the standardization and interpretation of the electrocardiogram,” Heart Rhythm the Official Journal of the Heart Rhythm Society, vol. 115, no. 10, pp. 1306–1324, 2007.

P. Camelliti, T. K. Borg, and P. Kohl, “Structural and functional characterisation of cardiac fibroblasts,” Cardiovascular Research, vol. 65, no. 1, pp. 40–51, 2005.

S. Verheule, T. Sat, T. IV. Everett et al., “Increased vulnerability to atrial fibrillation in transgenic mice with selective atrial fibrosis caused by overexpression of TGF-β1,” Circulation Research, vol. 94, no. 11, pp. 1458–1465, 2004.

F. Zhang, M. Tang, W. Wang, L. Chen, and D. O. Cardiology, “The effection of spironolactone on spontaneous hypertensive rat's left atrial fibrosis and it's atrial vulnerability,” Journal of Clinical Cardiology, 2015.

T. T. Win, B. A. Venkatesh, G. J. Volpe et al., “Associations of electrocardiographic P-wave characteristics with left atrial function, and diffuse left ventricular fibrosis defined by cardiac magnetic resonance: the PRIMERI study,” Heart Rhythm, vol. 12, no. 1, pp. 155–162, 2015.

A. Maheshwari, F. L. Norby, E. Z. Soliman et al., “Relation of prolonged p-wave duration to risk of sudden cardiac death in the general population (from the Atherosclerosis risk in communities study),” American Journal of Cardiology, vol. 119, no. 9, pp. 1302–1306, 2017.

S. Rosenheck, “Signal-averaged P wave in patients with paroxysmal atrial fibrillation,” Pacing and Clinical Electrophysiology Pace, vol. 20, no. 10, pp. 2577–2586, 1997.

P. G. Platonov, J. Carlson, M. P. Ingemansson et al., “Detection of inter-atrial conduction defects with unfiltered signal-averaged P-wave ECG in patients with lone atrial fibrillation,” Europace: European Pacing, Arrhythmias, and Cardiac Electrophysiology, vol. 2, no. 1, pp. 32–41, 2000.

J. Carlson, R. Johansson, and S. Olsson, “Classification of electrocardiographic P-wave morphology,” in Proceedings of the 39th IEEE Conference on Decision and Control, pp. 1804–1809, Sydney, NSW, Australia.

F. Sadiq Ali, A. Enriquez, D. Redfearn, and A. Baranchuk, “P-wave pseudonormalization after iatrogenic coronary sinus isolation,” Journal of Electrocardiology, vol. 49, no. 1, pp. 13-14, 2016.

P. Kellman and A. E. Arai, “Cardiac imaging techniques for physicians: late enhancement,” Journal of Magnetic Resonance Imaging, vol. 36, no. 3, pp. 529–542, 2012.