- 1School of Biological Science and Medical Engineering, Southeast University, Nanjing, China
- 2Division of Cardiology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China
- 3School of Mechanical, Medical and Process Engineering, Queensland University of Technology (QUT), Brisbane, QLD, Australia
Background: It is well known that the thrombus triggering stroke in patients with atrial fibrillation (AF) mainly comes from the left atrial appendage (LAA). This study aims to characterize the morphological and hemodynamic parameters and evaluate their differences between AF patients with and without a stroke history.
Methods: Cardiac CT images were obtained from AF patients with (n = 10) and without a history of stroke (n = 10). 3D models of the left atrium (LA) were reconstructed by processing the CT image, and the LA/LAA morphological parameters were measured. Computational fluid dynamics (CFD) simulations were performed to calculate the hemodynamic parameters in LA. The species-transport model and discrete phase model (DPM) were applied to analyze blood residual ratio and particle residual ratio, two qualitative parameters for thrombus formation and flow-out potential, respectively.
Results: There were significant differences in LAA actual depth (p = 0.002), and direct length (p = 0.049) between the non-stroke and stroke groups. Significant differences were also found in certain hemodynamic parameters. The blood residual ratio in LAA was significantly smaller in the stroke group than in the non-stroke group (p < 0.05). Moreover, the particle residual ratio within LAA was significantly smaller in the stroke groups than in the non-stroke group (p < 0.05).
Conclusion: There are significant differences in both morphological and hemodynamic parameters between AF patients with and without a stroke history. A high blood residual ratio in LAA confirms that thrombus is more likely to form in AF patients. A significantly smaller particle residual ratio in the stroke group may suggest the thrombus formed with LAA is more likely to flow out of LAA, leading to a higher risk of stroke. The proposed morphological and hemodynamic parameters may be additional risk factors that can be used to better risk stratify AF patients.
Introduction
Atrial fibrillation (AF) is the most common arrhythmia, and the prevalence is more than 33 million worldwide (1). Additionally, the occurrence of AF is predicted to increase rapidly as the average age of the population grows (2). Moreover, AF is well known as an independent risk factor for stroke (3, 4). It was reported that 15–18% of stroke is related to AF (5), and the autopsy study revealed that the thrombus inducing stroke in patients with AF mainly originates from the left atrial appendage (LAA) (6, 7), which is a tubular structure derived from the wall of the left atrium (LA) (8, 9).
LAA is normally categorized into four different morphologies (namely Cactus, Chicken Wing, Windsock and Cauliflower) (10–12), which are believed to be associated with different levels of stroke risk. However, it is a subjective classification with a high intro-observer difference. An objective morphological classification method is needed to quantify the association between LAA morphology and stroke risk.
Meanwhile, a non-invasive and repeatable method for simulating various fluid dynamic conditions, computational fluid dynamics (CFD) is widely applied in the cardiac hemodynamic analysis. Early studies on the CFD model of the LA and LAA were based on ideal models and subject-specific models reconstructed from CT data (13–15). A few hemodynamic parameters of patient-specific models, such as velocity, pressure, vorticity, and wall shear stress have been reported (16–18). The hemodynamic parameters in the LAA with the four different LAA morphologies have been discussed (19–21). Camara and colleagues (22, 23) investigated the change in flow field after LAA occlusion and performed a sensitivity analysis to identify the most relevant LA and LAA morphological parameters that are associated with blood flow dynamics. However, there were not intuitive descriptions of the thrombosis risk in previous studies.
To this end, this study aims to develop an objective morphological parameterization method and an intuitive description of the thrombogenesis risk. The species-transport model and discrete phase model (DPM) were applied to calculate the blood residual ratio and particle residual ratio, two qualitative parameters for thrombus formation and flow-out potential, respectively. Patients with persistent AF were divided into the non-stroke group and stroke group. Their LA/LAA morphological parameters and hemodynamic parameters were compared.
Materials and Methods
Study Population
Twenty patients with persistent AF from Jiangsu Province Hospital were included in this study. The patients were divided into two groups based on their stroke history: non-stroke group (n = 10) and stroke group (n = 10). The inclusion criteria were patients with persistent atrial fibrillation who underwent cardiac CT imaging before AF ablation. The exclusion criteria were: (1) valvular heart disease; (2) non cardioembolic stroke; (3) incomplete left atrium models due to insufficient cardiac CT imaging qualities. All the patients were on the blood thinner and the cardiac CT examination was taken timely for the stroke group after their stroke. The studies involving human participants were reviewed and approved by the Southeast University Human Research Ethics Committee with all participants providing written informed consent.
Images were acquired by using CT systems (Somatom Force and Definition; Siemens) with an ECG-gated dual-source single-energy protocol over the entire heart. The acquisition parameters were as follows: retrospective ECG-gated protocol; automatic tube voltage adjustment (Care kV; Siemens) from 80 to 120 kVp; automatic tube current modulation ranging between 150 and 400 mAs. Beta-blockers and nitroglycerine were used if necessary. Images was obtained following the administration of contrast medium (370 mg iodine/mL, Ultravist, Bayer Schering Pharma, Berlin, Germany) via a 20-G intravenous catheter, followed by 40 mL saline flush using a dual-head power injector. The injection duration was determined by scan time plus 8 s and the flow rate ranging from 4.0 to 6.0 mL/s was tailored to the body weight, heart rate (HR) and tube voltage.
Geometrical Model and Morphological Parameters
The LA geometry was segmented from the cardiac CT images using Mimics Medical 21.0 (Materialise, Inc., Leuven, Belgium) (see Figure 1), and then reconstructed with Geomagic studio (Geomagic, Inc., Morrisville, the United States). The morphological parameters were measured manually using Siemens NX (Siemens Digital Industries Software, Inc., Plano, the United States).
Figure 1. The modeling of LA and LAA. The red arrows point to PVs, and the yellow arrows point to the MV, which is cut manually. The purple area is the LA segmented using the region growing method. (A) The coronal image of cardiac CT. (B) The axial image of cardiac CT. (C) The sagittal image of cardiac CT. (D) The 3D model of LA and LV reconstructed from the cardiac CT.
In the measurement, the volume and surface area of the LA and LAA were recorded (see Figure 2A) and so were the cross section of the LAA orifice and its long and short axes (see Figure 2B). Meanwhile, to evaluate the tortuosity of the LAA, the direct length (the distance between the LAA tip and the orifice center) and the actual depth (the distance along the centerline from orifice center to the LAA tip) were measured (see Figure 2C). The tortuosity was calculated based on the LAA direct length and actual depth:
Figure 2. The morphology of the LA and LAA. (A) The LA structure. The LAA, MV and PVs (i.e., RSPV, RIPV, LSPV, and LIPV) are shown. (B) The long axis and short axis of the LAA orifice. (C) The direct length and actual depth of the LA.
Computational Fluid Dynamics Analysis
For each patient, based on patient-specific 3D geometrical model, the pulmonary veins (PVs) were set as the velocity-inlet, and the mitral valve (MV) was the pressure-outlet. For the inlets, the velocity profile was derived from the measurement of pulmonary venous profile for a persistent-AF patient, and detailed in our former study (24) as:
where V is the flow velocity for the inlet, and t is the time in the cardiac time.
The open and close conditions of MV were set by the high- and low-pressure, which were expressed as:
The diastole is from 0 to 0.26 s, and the systole is from 0.26 to 0.7 s.
Tetrahedral-element meshes were generated by Ansys Meshing using a maximum element size of 0.8 mm and a minimum size of 0.2 mm, which were confirmed to be sufficient by performing a denser mesh.
The CFD analysis were performed on the node (2*Intel Xeon E5-2680 v3) in the High-Performance Computing Center of Collaborative Innovation Centre of Advanced Microstructures using the commercial software Fluent 2019 R3 (Ansys, Inc., Canonsburg, the United States).
In the simulation, the blood was set as incompressible Newtonian fluid with density ρ = 1,060 kg/m3 and viscosity μ = 0.0035 Pa⋅s. The flow was simulated by the realizable k-ε turbulence model with enhanced wall treatment, and the walls were assumed as rigid and with no-slip conditions to represent the worst AF symptom. The coupled second order upwind method was used for the simulations, and the time step was set as 0.001 s.
Meanwhile, to evaluate the hemodynamics in a more intuitive way, the species-transport model was employed to study the mechanisms of thrombus formation and the DPM analysis was performed to stratify the thrombus dynamics.
For the species-transport model, the mass fraction of each species was estimated through the convective-diffusion equation:
where Yi is the mass fraction of the ith species, and Ji is the diffusive flux of the ith species, which is calculated as follows in turbulence:
where SCt is the turbulent Schmidt number and Di, m is the diffusion coefficient of the ith species.
For the DPM, the trajectory of discrete-phase particles is solved by integrating the particle force differential equation in the Laplace coordinate system. The force balance equation of particles in the Cartesian coordinate system is in calculated as follows:
where FD(u−up) is the drag force per unit mass of particles:
where u is fluid velocity, up is particle velocity, μ is fluid dynamic viscosity, ρ is fluid density, ρp is particle density, dp is particle diameter, Re is relative Reynolds number, and CD is drag coefficient.
For the species-transport model, the initial blood in the LA was defined as one phase, which would be replaced by the fresh blood (defined as the other phase) flowing from the PVs. The residual rate of the initial blood was calculated, and the higher blood residual ratio corresponds to the more probability of thrombus formation in the LAA. In order to study the mechanism of thrombus flowing-out from the LAA, particles were set on the LAA walls at the beginning of the third cardiac cycle, the stable flow condition. The risk of thrombus flowing-out was represented by the residual ratio of particles in the LAA, and a higher particle residual ratio corresponds to a higher ratio of particles flowing out of the LAA, and thus a higher risk of stroke.
Statistical Analysis
Statistical analysis was performed the open-source software Jamovi. The data were tested for normal distribution using Shapiro-Wilk test and the results were described as mean ± SD or median (IQR). Independent samples t-test or Mann-Whitney U test was used to compare the geometrical and hemodynamic parameters between the non-stroke group and stroke group. It was considered statistically significant when p-value < 0.05.
Results
Patient Characteristics
Patient characteristics are shown in Table 1. No significant (p > 0.05) differences existed in age, hypertension, diabetes and left ventricular ejection fraction (LVEF) between the two groups.
Morphological Parameters
Morphological parameters of the LA and LAA are presented in Table 2. Both the actual depth and direct length of LAA were significantly larger in the non-stroke group than in the stroke group (p < 0.05). The surface area and tortuosity of LAA showed certain differences between the two groups but they were not statistically significant. There were no significant differences in the LA volume, LA surface area, LAA volume, LAA orifice long radius, LAA orifice short radius and LAA orifice ovality (the ratio of long radius to short radius).
Blood Residual Ratio
The residual ratios of initial blood at the end of each cardiac cycle in the LA and LAA were recorded, which represent how efficient the blood clears out from LAA at each cardia cycle. An example of the simulation results of the blood residual ratio is shown in Figure 3A. The ratios were compared between the non-stroke and stroke groups in Table 3.
Figure 3. The results of species-transport and DPM models. (A) The contour of the initial blood. The initial blood show in blue, and the fresh blood show in red. (B) The particle distribution.
As shown in Table 3, from the third cycle, the residual ratios for the LA were significantly smaller in the stroke group than in the non-stroke group (p < 0.05). While for the LAA, the residual ratios were significantly smaller in the stroke group than in the non-stroke group (p < 0.05) from the first cycle.
Meanwhile, we also calculated the time when the residual ratios of the initial blood reduced to 20, 15, and 10%. No significant differences existed in the time parameters in the LA (p > 0.05) between the two groups, but the time parameters for LAA were significantly larger in the stroke group than in the non-stroke group (Table 3).
Particle Residual Ratio
DPM analysis was used to quantify the efficiency of thrombus particles within LAA flowing out of LAA. Figure 3B shows an example of the particle distributions within LAA in a patient case. The residual ratio of the particles in LAA at each cardiac cycle was recorded, which are shown in Table 4. The particles were released at the third cardiac cycle, so the results were recorded from the third cycle. No significant (p > 0.05) differences existed in the residual ratios of particles in LAA between the two groups at the beginning three cycles, but the particle residual ratios became significantly smaller in the stroke groups than in the non-stroke group (p < 0.05) from the fourth cardiac cycle.
Discussion
AF is a common heart disease worldwide and known as an independent risk factor for stroke. The CHA2DS2-VASc score is an effective method for diagnosing whether a patient with persistent AF needs to be anticoagulated (25). However, the patients with low scores are still at risk of stroke (26). In order to further explore the mechanism of thrombus formation and stroke risk, and develop additional risk factors to better assess the risk of stroke in AF patients, we developed a quantitative strategy to characterize the morphological and hemodynamic parameters in LA/LAA. We divided the patients with persistent AF into two groups: non-stroke group and stroke group, and measured the morphological parameters of the LA and LAA and analyzed the hemodynamic characteristics.
Our results confirmed that there were some differences in the morphological parameters of LA/LAA between the non-stroke and stroke groups. Compared with the non-stroke group, the LAAs of the stroke group had a smaller actual depth and direct length. The results are consistent with a previous study by (27), who also proposed the feasibility of using morphological parameters for clinical diagnosis of AF patients. A combination of the CHA2DS2-VASc score and morphological parameters may provide a more accurate guide to diagnose whether patients with persistent AF need to be anticoagulated.
For the first time, we simulated the blood flow renewal process in the LA by applying the species-transport model. The results showed that the blood flow in the LA and LAA was more likely to be renewed in the stroke group than in the non-stroke group. This seems to contradict our conjecture that stagnant blood results in blood clots and stroke group should have a lower renewal ratio (higher blood residual ratio) based on the traditional hypothesis. One possible explanation may be that although the stroke group has a lower blood residual ratio, the residual ratios in both groups were very high within their LAA and the risk of thrombus formation were both high within the LAA in patients with AF. The lower blood residual is correlated with the shorter LAA length in the stroke group. The results suggest that the LAA is an important factor of thrombosis and the morphological parameters of the LAA are more notable in the stroke groups than in the non-stroke group. It was reported that thrombus formation is likely to happen in areas with low flow rates (28), so the areas under certain velocity may be a better way to show the risk of thrombus formation. Further extensive patient studies are required to further explore the underlying mechanism.
The process of thrombus discharged from the LAA was also studied by employing the DPM analysis. Unsurprisingly, the results showed that the thrombus had a higher propensity to expel in the stroke group than in the non-stroke group. This result showed that there is a higher likelihood of thrombus flowing out of LAA in the stroke group, which may suggest that a higher discharge of thrombus in LAA may lead to a higher risk of stroke. The proposed particle analysis may act as a useful tool to assess the risk of cardioembolic stroke in patients with AF.
Some limitations may exist in our study: First, only twenty patients with persistent AF were enrolled; Second, the wall of LA/LAA was assumed to be rigid wall. It will provide more convincing results, if a larger cohort of patients is enrolled and fluid-structure-interaction method is employed.
Conclusion
There are significant differences in both morphological and hemodynamic parameters between patients with and without a stroke history. Particularly, the LAAs of the stroke group is shorter and less-tortuous than the non-stroke group. Both groups have a high blood residual ratio in LAA which may suggest the thrombus formation is likely to happen in all the AF patients. For the blood clots formed within LAA, the flowing out risk of the blood clots is significantly higher in the stroke group, suggesting the particle residual ratio could be used as a hemodynamic parameter to assess stroke risk in AF patients. The proposed morphological and hemodynamic parameters may be additional risk factors that can be used to better risk stratify the patients with AF.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the Southeast University Human Research Ethics Committee. Written informed consent was not required for this study, in accordance with the local legislation and institutional requirements.
Author Contributions
LW and ZW made major contribution to the writing of the article, model construction, and data analysis. Z-YL was responsible for overall design and execution of the research. RF worked on CT image segmentation and morphological parameters measurement. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by the National Nature Science Foundation of China (11972118, 61821002, and 11772093) and ARC (DP200103492 and DP200101970).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Wijesurendra RS, Casadei B. Mechanisms of atrial fibrillation. Heart. (2019) 105:1860–7. doi: 10.1136/heartjnl-2018-314267
2. Iwasaki Y, Nishida K, Kato T, Nattel S. Atrial fibrillation pathophysiology. Circulation. (2011) 124:2264–74. doi: 10.1161/CIRCULATIONAHA.111.019893
3. Warraich HJ, Gandhavadi M, Manning WJ. Mechanical discordance of the left atrium and appendage. Stroke. (2014) 45:1481–4. doi: 10.1161/STROKEAHA.114.004800
4. Yaghi S, Song C, Gray WA, Furie KL, Elkind MSV, Kamel H. Left atrial appendage function and stroke risk. Stroke. (2015) 46:3554–9. doi: 10.1161/STROKEAHA.115.011273
5. Wolf PA, Abbott RD, Kannel WB. Atrial fibrillation as an independent risk factor for stroke: the Framingham study. Stroke. (1991) 22:983–8. doi: 10.1161/01.STR.22.8.983
6. Blackshear JL, Odell JA. Appendage obliteration to reduce stroke in cardiac surgical patients with atrial fibrillation. Ann Thorac Surg. (1996) 61:755–9. doi: 10.1016/0003-4975(95)00887-X
7. Wang Y, Biase LD, Horton RP, Nguyen T, Morhanty P, Natale A. Left atrial appendage studied by computed tomography to help planning for appendage closure device placement. J Cardiovasc Electrophysiol. (2010) 21:973–82. doi: 10.1111/j.1540-8167.2010.01814.x
8. Al-Saady NM, Obel OA, Camm AJ. Left atrial appendage: structure, function, and role in thromboembolism. Heart. (1999) 82:547. doi: 10.1136/hrt.82.5.547
9. Di Biase L, Burkhardt JD, Mohanty P, Sanchez J, Mohanty S, Horton R, et al. Left atrial appendage. Circulation. (2010) 122:109–18. doi: 10.1161/CIRCULATIONAHA.109.928903
10. Ernst G, Stöllberger C, Abzieher F, Veit-Dirscherl W, Bonner E, Bibus B, et al. Morphology of the left atrial appendage. Anat Rec. (1995) 242:553–61. doi: 10.1002/ar.1092420411
11. Di Biase L, Santangeli P, Anselmino M, Mohanty P, Salvetti I, Gili S, et al. Does the left atrial appendage morphology correlate with the risk of stroke in patients with atrial fibrillation?: results from a multicenter study. J Am Coll Cardiol. (2012) 60:531–8. doi: 10.1016/j.jacc.2012.04.032
12. Beigel R, Wunderlich NC, Ho SY, Arsanjani R, Siegel RJ. The left atrial appendage: anatomy, function, and noninvasive evaluation. JACC Cardiovasc Imaging. (2014) 7:1251–65. doi: 10.1016/j.jcmg.2014.08.009
13. Zhang LT, Gay M. Characterizing left atrial appendage functions in sinus rhythm and atrial fibrillation using computational models. J Biomech. (2008) 41:2515–23. doi: 10.1016/j.jbiomech.2008.05.012
14. Di Martino ES, Bellini C, Schwartzman DS. In vivo porcine left atrial wall stress: computational model. J Biomech. (2011) 44:2589–94. doi: 10.1016/j.jbiomech.2011.08.023
15. Di Martino ES, Bellini C, Schwartzman DS. In vivo porcine left atrial wall stress: effect of ventricular tachypacing on spatial and temporal stress distribution. J Biomech. (2011) 44:2755–60. doi: 10.1016/j.jbiomech.2011.09.003
16. Koizumi R, Funamoto K, Hayase T, Kanke Y, Shibata M, Shiraishi Y, et al. Numerical analysis of hemodynamic changes in the left atrium due to atrial fibrillation. J Biomech. (2015) 48:472–8. doi: 10.1016/j.jbiomech.2014.12.025
17. Vedula V, George R, Younes L, Mittal R. Hemodynamics in the left atrium and its effect on ventricular flow patterns. J Biomech Eng. (2015) 137:111003. doi: 10.1115/1.4031487
18. Otani T, Al-Issa A, Pourmorteza A, McVeigh ER, Wada S, Ashikaga H. A computational framework for personalized blood flow analysis in the human left atrium. Ann Biomed Eng. (2016) 44:3284–94. doi: 10.1007/s10439-016-1590-x
19. Valinoti M, Fabbri C, Turco D, Mantovan R, Pasini A, Corsi C. 3D patient-specific models for left atrium characterization to support ablation in atrial fibrillation patients. Magn Reson Imaging. (2018) 45:51–7. doi: 10.1016/j.mri.2017.09.012
20. Masci A, Alessandrini M, Forti D, Menghini F, Dedé L, Tomasi C, et al. A proof of concept for computational fluid dynamic analysis of the left atrium in atrial fibrillation on a patient-specific basis. J Biomech Eng. (2019) 142:011002. doi: 10.1115/1.4044583
21. Masci A, Barone L, Dedè L, Fedele M, Tomasi C, Quarteroni A, et al. The impact of left atrium appendage morphology on stroke risk assessment in atrial fibrillation: a computational fluid dynamics study. Front Physiol. (2019) 9:1938. doi: 10.3389/fphys.2018.01938
22. García-Isla G, Olivares AL, Silva E, Nuñez-Garcia M, Butakoff C, Sanchez-Quintana D, et al. Sensitivity analysis of geometrical parameters to study haemodynamics and thrombus formation in the left atrial appendage. Int J Numer Method Biomed Eng. (2018) 34:e3100. doi: 10.1002/cnm.3100
23. Aguado AM, Olivares AL, Yagüe C, Silva E, Nuñez-García M, Fernandez-Quilez Á, et al. In silico optimization of left atrial appendage occluder implantation using interactive and modeling tools. Front Physiol. (2019) 10:237. doi: 10.3389/fphys.2019.00237
24. Fang R, Li Y, Zhang Y, Chen Q, Liu Q, Li Z. Impact of left atrial appendage location on risk of thrombus formation in patients with atrial fibrillation. Biomech Model Mechanobiol. (2021) 20:1431–43. doi: 10.1007/s10237-021-01454-4
25. Lip GYH, Nieuwlaat R, Pisters R, Lane DA, Crijns HJGM. Refining clinical risk stratification for predicting stroke and thromboembolism in atrial fibrillation using a novel risk factor-based approach: the euro heart survey on atrial fibrillation. Chest. (2010) 137:263–72. doi: 10.1378/chest.09-1584
26. Lee JM, Kim J-B, Uhm J-S, Pak H-N, Lee M-H, Joung B. Additional value of left atrial appendage geometry and hemodynamics when considering anticoagulation strategy in patients with atrial fibrillation with low CHA2DS2-VASc scores. Heart Rhythm. (2017) 14:1297–301. doi: 10.1016/j.hrthm.2017.05.034
27. Bieging ET, Morris A, Chang L, Dagher L, Marrouche NF, Cates J. Statistical shape analysis of the left atrial appendage predicts stroke in atrial fibrillation. Int J Cardiovasc Imaging. (2021) 37:2521–7. doi: 10.1007/s10554-021-02262-8
Keywords: atrial fibrillation, left atrial appendage, stroke, hemodynamics, computational fluid dynamics
Citation: Wang L, Wang Z, Fang R and Li Z-Y (2022) Evaluation of Stroke Risk in Patients With Atrial Fibrillation Using Morphological and Hemodynamic Characteristics. Front. Cardiovasc. Med. 9:842364. doi: 10.3389/fcvm.2022.842364
Received: 23 December 2021; Accepted: 11 April 2022;
Published: 29 April 2022.
Edited by:
Liza Thomas, University of New South Wales, AustraliaReviewed by:
Kelvin Kian Loong Wong, University of Saskatchewan, CanadaGeorgios Benetos, National and Kapodistrian University of Athens, Greece
Copyright © 2022 Wang, Wang, Fang and Li. 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: Zhi-Yong Li, enlsaWNhbUBnbWFpbC5jb20=
†These authors have contributed equally to this work