# A computational bilayer surface model of human atria

Simon Labarthe, Edward Vigmond, Yves Coudière, Jacques Henry, Hubert Cochet, Pierre Jaïs

► To cite this version:

Simon Labarthe, Edward Vigmond, Yves Coudière, Jacques Henry, Hubert Cochet, et al. A computational bilayer surface model of human atria. FIMH 2013 - 7th International Conference on Functional Imaging and Modeling of the Heart, Jun 2013, London, United Kingdom. pp.27-34, 10.1007/978-3-642-38899-6_4. hal-01606999

HAL Id: hal-01606999

https://hal.science/hal-01606999

Submitted on 6 Jun 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# A Computational Bilayer Surface Model of Human Atria

Simon Labarthe, Edward Vigmond, Yves Coudière, Jacques Henry, Hubert Cochet, Pierre Jais

► To cite this version:

Simon Labarthe, Edward Vigmond, Yves Coudière, Jacques Henry, Hubert Cochet, et al. A Computational Bilayer Surface Model of Human Atria. FIMH 2013 - 7th International Conference on Functional Imaging and Modeling of the Heart, Jun 2013, London, United Kingdom. Springer, 2013, Lecture Notes In Computer Sciences.

HAL Id: hal-00802104

https://hal.inria.fr/hal-00802104v1

Submitted on 19 Mar 2013 (v1), last revised 28 Jan 2013 (v2)

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
# A Computational Bilayer Surface Model of Human Atria

Simon Labarthe<sup>1,2,3</sup>, Edward Vigmond<sup>1,2,3</sup>, Yves Coudière<sup>1,2,3</sup>,
Jacques Henry<sup>1,2,3</sup>, Hubert Cochet<sup>3,4</sup>, Pierre Jaïs<sup>3,4</sup>

1 - IMB - Université Bordeaux Segalen - Université Bordeaux 1,
Bordeaux, France

2 - INRIA, Bordeaux, France

3 - Liryc Institute - CHU / Université de Bordeaux / INSERM
U1045, Bordeaux, France

4 - Hôpital Haut-Lévêque CHU Bordeaux Université Bordeaux I

March 19, 2013

## Abstract

A bilayer surface model of human atria is presented. This model allows inclusion of transmural heterogeneities while maintaining small computational costs associated with surface models.

A physiological fibre structure is constructed on an imaged atrial geometry which includes the main structures of fibres in the right and left atria, including the transmural heterogeneities, and the main transsepital connexions. The validity of the corresponding bilayer mathematical model is assessed by comparing the case of two sheets of orthogonal fibres to a three dimensional slab of tissue. A comparison of bilayer and monolayer simulations of sinus propagation is finally proposed.

This atrial bilayer model is a good tool for performing complex atrial simulations in pathological cases.

## 1 Introduction

Atrial computational models have proven their utility in fundamental or clinical research by giving insight into atrial fibrillation (AF) initiation and perpetuation, radio-frequency ablation, and the development of signal analysis techniques [1].
As the electrical activity spreads across the atria, propagation patterns and P-wave duration are strongly influenced by atrial fibre structure [2, 3]. Modeling the anisotropy of the atria is, therefore, a major step in the construction of atrial models. This construction can result from an image-based process on ex-vivo hearts [3] or from a semi-automatic rule-based approach for acquired geometries of in-vivo atria [2, 4].

Atrial three dimensional models are usually divided into two main categories: surface and volumetric models. Whereas the former propose a homogeneous transmural fibre distribution and computationally light simulations, the latter capture important transmural heterogeneities or layer dissociations as observed in animal models [5, 6, 7] in exchange for increased computational costs [1].

As a compromise between the surface two dimensional and volumetric three dimensional models, we propose a two layer computational model of the atria that can reproduce complex three dimensional patterns while keeping surface computational load. Such a two layer approach has also been used in [8] on idealized geometries, but has not been demonstrated to reproduce proper behavior.

We first describe the construction of the model. We then give some mathematical arguments to the validity of the bilayer approach. We finally compare simulations of a normal action potential (AP) propagation in a bilayer and a monolayer model.

# 2 Method

## 2.1 Imaging the Atria

Bi-atrial geometry was acquired in vivo using ECG-gated contrast-enhanced multi-detector computed tomography (MDCT) in a patient with no history of atrial disorder. Imaging was performed using a 64-slice CT scanner (SO-MATOM Definition, Siemens Medical Solutions, Forchheim, Germany) during intravenous injection of a 120 mL bolus of iodinated contrast agent at the rate of 4 mL s⁻¹. ECG-gating was set for the acquisition window to occur at end-systole, when atrial motion is minimal. Acquisition parameters were: slice thickness 0.5 mm, tube Voltage 120 kV, maximum tube current 850 mAs, and gantry rotation time 330 ms. An imaging volume comprising the whole atria was reconstructed with a voxel size of 0.5x0.4x0.4 mm. Transaxial images were imported in DICOM format in a local database of the software OsiriX 3.6.1 (OsiriX foundation, Geneva, Switzerland). The endocardial surface was segmented automatically using region growth segmentation.
Volume rendering reconstructions were used to remove non atrial structures
from the segmentation. Pulmonary veins were cut several centimetres from
their ostia because myocardium is known to extend to the first centimeter
beyond the pulmonary veins ostia. The resulting segmentation was used
to compute a 3D mesh in vtk format using the software CardioviZ3D (IN-
RIA Asclepios, Sophia Antipolis, France). The mesh was smoothed with
weighted Laplacian smoothing and optimized with the software freeYams
to obtain a mesh suitable for calculation that contained 361483 nodes and
715034 elements with a mean diameter of 0.492 mm.

## 2.2 Construction of the Geometry of the Bilayer Model

### 2.2.1 Algorithms.

Two methods were used for fibres definition and mesh completion.

*Fibre definition.* The fibre structure was defined with a rule-based semi-automatic method [4] that was founded on two atrial characteristics : (1) atrial geometry can be divided in basic cylinder-like structures with circumferential fibre orientation. (2) Fibre arrangement is continuous. The atrial mesh was first manually decomposed in cylinder-like elements. The local circumferential structure was semi-automatically obtained by determining the cylindrical axis as an eigen vector of an inertia tensor. Finally, an in-painting algorithm extended continuously the fibre direction on the whole geometry.

*Manual mesh completion.* Small structures that were partially acquired, such as transseptal connexions, were manually constructed with the following workflow : the points of insertion of the new structure in the mesh were identified. Plane surfaces connecting those points were defined and meshed with the software Gmsh. The structure was then smoothed by a Laplacian method with the software MeshLab (3D-CoForm Project, Italy).

### 2.2.2 Main modeling options.

In order to describe the atria as a bilayer structure, we need to define some physiological a-priori modeling options, specially focusing on transmural heterogeneities.

*Right Atrium (RA).* According to histological descriptions of the atria [10], the RA can be seen as a monolayer structure, except in the appendage (RAA) and in the roof where the crista terminalis (CT) and the pectinate muscles (PM) bring thickness and fibres inhomogeneities. A first layer modeled the epicardial part of the RA including its main fibre structures : the
circular fibres that surrounded the superior veina cava (SVC) crossed the CT and covered the RAA to blend with the circular bundle of the vestibule. A second layer was applied on the location of PM and the CT to include the longitudinal fibre direction along those structures in the endocardial part of the RA. A poor conduction zone was introduced along the left side of the sino-atrial node (SAN) as assessed in [12].

*Left Atrium (LA)*. The left atrium (LA) fibres structure can be described as a two layer geometry. The main transmural heterogeneities are concentrated in the anterior wall and on the pulmonary veins (PV). The circumferential and oblique fibres of the anterior epicardium blend with the Bachmann's bundle whereas the endocardial fibres of the septoatrial bundle run from the vestibule to the dome of the LA. The fibre structure of the PV can show strong heterogeneities with non uniform arrangement and transmural discontinuities of fibres directions [11].

*Transseptal connexions*. The chambers of this atrial model were electrically insulated and connected by the three main interatrial connexions: the BB, the fossa ovalis (FO) and the coronary sinus (CS). The BB was segmented manually as a three-dimensional structure and then reduced to two surfaces: the first one corresponded to the posterior face of the BB and inserted into the septal wall of LA while the second one was constructed as a mean surface of the volumetric BB and blended into the septopulmonary bundle of the LA. The FO and the CS were identifiable after the segmentation step and connected manually.

## 2.3 Modeling the AP Propagation

The monodomain equation was used to describe the propagation of electrical conductivity

$$ \mathrm{div}(\bar{\sigma}_i \nabla V_m) = \beta (C_m \partial_t V_m + I_{\mathrm{ion}}(V_m, t, \bar{\zeta})) \quad (1) $$

where $V_m$ is the transmembrane voltage, $\bar{\sigma}_i$ is the conductivity tensor, $\beta$ is the ratio of membrane surface to volume in which it is contained, $C_m$ is the membrane capacitance (1F cm$^{-2}$), and $I_{\mathrm{ion}}$ is the transmembrane current which is a function of voltage, time, and $\bar{\zeta}$, a set of state variables.

*Mathematical assessment of the bilayer model validity.* A bilayer expression of this mathematical model has been asymptotically derived in the limit of vanishing thickness in [9] from a bilayer volumetric monodomain model. That reads:

$$ \beta(C_m \partial_t V_m^{(k)} + I_{\mathrm{ion}}(V_m^{(k)}, \bar{\zeta}^{(k)})) = \mathrm{div}(\sigma'^{(k)} \nabla V_m^{(k)}) + (-1)^k \sigma_c (V_m^{(1)} - V_m^{(2)}) \quad (2) $$
<figure><img src="image_2.png" /></figure>

Figure 1: Left:Inferior view of the atria. Right: Superior view. All distinct surfaces are colored uniquely. SVC-superior vena cava. IVC-inferior vena cava. BB-bachmanns bundle. CS-coronary sinus. BZ-block zone. SAN-sinoatrial node. RSPV-right superior pulmonary vein. RIPV-right inferior pulmonary vein.LSPV-left superior pulmonary vein. LIPV-left inferior pulmonary vein. RAA-right atrial appendage. LAA-left atrial appendage. CT-crist terminalis. PM- pectinate muscles.

in the layer $k = 1, 2$, where $V_m^{(k)}$ is the transmembrane potential in the respective layer and $\sigma_c = \frac{3}{h^2} \frac{\sigma_3^{(1)} \sigma_3^{(2)}}{\sigma_3^{(1)} + \sigma_3^{(2)}}$. $\sigma'^{(k)}$ and $\sigma_3^{(k)}$ are respectively the diffusion coefficients in the layer tangent plane and in the transmural direction. $h$ is the thickness of the layers.

Three dimensional, bilayer and monolayer finite element $P_1$-Lagrange simulations have been performed on, respectively, a 3D slab and its 2D equivalent for various $h$, with the Beeler-Reuter model and a semi-implicit Euler scheme. The activation maps of the mean potential in the thickness of each simulation have been compared. Variations of the coupling coefficient $\sigma_c$ were uniformly applied to check the robustness of the model under physiological thickness, far from the asymptotic regime.

*Bilayer atrial model.* The ionic model was described by the Courtemanche human atrial action potential model [13]. The sodium conductance was increased by a factor of four from the published value to achieve realistic propagation speeds. $\sigma_i$ and $\beta$ were set to (0.3, 0.04) S m$^{-1}$ and 0.14 m$^{-2}$. The structural block zone close to the SAN [12] was implemented with a passive ionic model with a resting level -80 mV and an intracellular conductivity reduced by a factor of 100. The SAN is located between the block zone and the crista terminalis. Accordingly, a small region was stimulated
<figure><img src="image_3.png" /><figcaption>Figure 2: Top: Anterior view of the epicardial fibres arrangement of both atria and zoom on the left inferior part of the septopulmonary bundle. Bottom: Map of the angle in radians formed by the endo and epicardial fibre direction. Dark blue (negative values): only the epicardial layer is modeled. Light blue: 0 rad. Same fibre direction in the endo and epicardial layer. Red : $\pi/2$ rad. Orthogonal directions.</figcaption></figure>

there to simulate the site of first excitation.

The equations were solved using the finite element method. Each surface was composed of linear triangular elements while connections between surfaces were implemented as one dimensional finite elements. 134510 line elements were added to couple the surfaces together. The system was solved using the CARP cardiac simulator [14] with a fixed time step of 25 s using a Crank-Nicolson approach to solve eq. 1. On a desktop workstation with Intel(R) Xeon X5650 CPUs running at 2.67GHz, it took seven seconds to simulate one ms on ten cores, while on a desktop with an NVidia Quadro 4000 (256 cores) it took 18 seconds per millisecond running it on the GPU.
# 3 Results

## 3.1 Fibres architecture

The main fibre structures as described in [10] are correctly reproduced by our model. In fig. 2, we can see that the differences between the endo and epicardial fibres orientations are mainly localized in the CT and PM in the RA, and in the anterior and posterior wall, together with the PVs in the LA, where the fibres in both layers can be orthogonal.

## 3.2 Simulations

*Validity of the bilayer approach.* The bilayer mathematical model has been shown to have second-order convergence towards the volumetric model when $h \to 0$. For physiological $h$, a correction can be brought to the coupling value $\sigma_c$ so that the accuracy of the bilayer surface model is preserved : complex depolarisation patterns can be obtained with a relative error of 3% respectively to 3D simulations (figure 3). We refer to [9] for a deeper study of this model.

*Bi-atrial simulations.* The propagation of electrical activity showed a normal pathway: it first propagated from SAN to the BB along the CT. LA was first activated by the BB, and then by the FO and the CS. Both atria were totally activated after 196 ms - fig. 4. The sinus wave ended between the left PVs and in the LAA. The effect of transmural heterogeneities are visible along the PM on the activation map. The comparison of the monolayer and bilayer activation maps showed that the main differences of propagation are in accordance with the map of fibres differences - fig. 2. The bilayer wave is faster in the PM direction, around the mitral valve and on the left PVs, and slower on the right PVs and on the anterior and posterior wall of the LA.

# 4 Discussion and conclusions

Clinical uses of atrial modeling or accurate atrial models calibrated for fundamental research need personalized fibres orientations, and then efficient methods of construction. The work-flow presented here is still time-consuming and has to be improved in view of such a systematic use. Other methods that are being developed for 3D models [2] could be adapted to bilayer surface models. Statistical methods of registration to transport this model to personalized geometries could be another way to shorten this process[15].
<figure><img src="image_4.png" /><figcaption>Figure 3: Left : Snapshots of transmurally averaged transmembrane poten-
tial maps at t = 16ms for a 1mm thick slab of tissue with two layers of
orthogonal fibres. Top: Volumetric model. Bottom: bilayer surface model
(corrected coefficient). Right : Relative L^1 error on activation map versus
the thickness of the tissue in log-log scale for the monolayer (blue, stars),
the bilayer (green, x) and the corrected bilayer model (red, +)</figcaption></figure>

Modeling assumptions were set prior to model construction. The bilayer
representation of some zones (LA,CT,PM) and the monolayer representation
of others (RA) were guided by histological descriptions. If more complex
transmural heterogeneities are observed, the model can be locally improved
by adding layers or by simulating small entities with a volumetric model.

A normal propagation was obtained. The main qualitative character-
istics of an atrial normal AP propagation were correctly reproduced : the
transseptal pathway and the endpoint of the sinus wave are in a physiolog-
ical range. The total activation time is larger than those reported by the
literature [16] but this can be tuned by further adjusting a combination of
sodium channel conductances and tissue conductivities.

The bilayer model demonstrated its capacity to capture some perturba-
tions of action potential propagation induced by transmural heterogeneities;
that represents an important improvement over monolayer surface models.
We could show that those perturbations were reliable on idealized geome-
tries, indicating that this model could be efficiently used on physiological
geometries. The demonstration of the potential of the bilayer approach
could be completed with a comparison between a volumetric model and its
equivalent bilayer surface model.

This bi-atrial bilayer model proved to be a good trade-off between surface
and volumetric models; it reproduced transmural heterogeneities that were
specific to volumetric models with a small computational load. It is an
<figure><img src="image_5.png" /><figcaption>Figure 4: Left : posterior view of the activation map with 5 ms isolines. Dark blue : 0 ms. Red : 196 ms. The deformation of the front is clearly visible in the PM. Right : posterior view of the difference between activation map of the bilayer (epicardial layer, CT, PM, endocardial layer of LA) and monolayer (epicardial layer) surface models. Colors are ranging from -10 ms (dark blue) to +10ms (red). 2 ms isodifference lines are shown.</figcaption></figure>

efficient tool to investigate the impact of the transmural heterogeneities in
pathological situations.

References

[1] Olaf Dössel, Martin W. Krueger, Frank M. Weber, et al. Computational modeling of the human atrial anatomy and electrophysiology. *Medical & Biological Engineering & Computing*, 50:773–799, 2012.

[2] Martin W. Krueger, Viktor Schmidt, Catalina Tobón et al. Modeling atrial fiber orientation in patient-specific geometries: A semi-automatic rule-based approach. In *FIMH'11*, pages 223–232, 2011.

[3] Jichao Zhao, Timothy D. Butters, Henggui Zhang et al. An image-based model of atrial muscular architecture / clinical perspective. *Circulation: Arrhythmia and Electrophysiology*, 5(2):361–370, 2012.

[4] Simon Labarthe, Yves Coudière, Jacques Henry, and Hubert Cochet. A Semi-Automatic Method To Construct Atrial Fibre Structures : a
Tool for Atrial Simulations. In *CinC 2012 - Computing in Cardiology*, volume 39, pages 881-884, Krakow, Pologne, November 2012.

[5] Katayoun Derakhchan, Danshi Li, Marc Courtemanche et al. Method for simultaneous epicardial and endocardial mapping of in vivo canine heart: Application to atrial conduction properties and arrhythmia mechanisms. *Journal of Cardiovascular Electrophysiology*, 12(5):548-555, 2001.

[6] Jens Eckstein, Bart Maesen, Dominik Linz et al. Time course and mechanisms of endo-epicardial electrical dissociation during atrial fibrillation in the goat. *Cardiovascular Research*, 89(4):816-824, 2011.

[7] Yoav Michowitz, Shiro Nakahara, Tara Bourke et al. Electrophysiolog-ical differences between the epicardium and the endocardium of the left atrium. *Pacing and Clinical Electrophysiology*, 34(1):37-46, 2011.

[8] Ali Gharaviri, Sander Verheule, Jens Eckstein et al. A computer model of endo-epicardial electrical dissociation and transmural conduction during atrial fibrillation. *Europace*, 14(suppl 5):v10-v16, 2012.

[9] Mostafa Bendahmane, Yves Coudière, Jacques Henry, and Simon Labarthe. An asymptotic two-layers monodomain model of cardiac elec-trophysiology in the atria. *To be published*

[10] Siew Yen Ho, Robert H. Anderson, and Damián Sánchez-Quintana. Atrial structure and fibres: morphologic bases of atrial conduction. *Cardiovascular Research*, 54(2):325-336, 2002.

[11] Tsukasa Saito, Kenji Waki and Anton E. Becker Left atrial myocardial extension onto pulmonary veins in humans: anatomic observations relevant for atrial arrhythmias. *Journal of cardiovascular electrophysiology*, 11:888-894, 2000.

[12] Vadim V. Fedorov, Richard B. Schuessler, Matthew Hemphill et al. Structural and functional evidence for discrete exit pathways that connect the canine sinoatrial node and atria. *Circulation Research*, 104(7):915-923, 2009.

[13] Marc Courtemanche, Rafael J. Ramirez, and Stanley Nattel. Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model. *American Journal of Physiology - Heart and Circulatory Physiology*, 275(1):H301-H321, 1998.
[14] Edward J Vigmond, Matt Hughes, G Plank, and L. Joshua Leon. Computational tools for modeling electrical activity in cardiac tissue. *Journal of Electrocardiology*, 36, Supplement 1(0):69–74, 2003.

[15] Nicolas Toussaint, Maxime Sermesant, Christian T. Stoeck et al. In vivo human 3d cardiac fibre architecture: Reconstruction using curvilinear interpolation of diffusion tensor images. *MICCAI 2010, volume 6361 of Lecture Notes in Computer Science*, pages 418–425. Springer Berlin Heidelberg, 2010.

[16] David M. Harrild and Craig S. Henriquez. A computer model of normal conduction in the human atria. *Circ Res*, 87(7):e25–36, 2000.