1
|
Borja MG, Martinez-Legazpi P, Nguyen C, Flores O, Kahn AM, Bermejo J, Del Álamo JC. Deriving phenotype-representative left ventricular flow patterns by reduced-order modeling and classification. Comput Biol Med 2024; 179:108760. [PMID: 38944903 DOI: 10.1016/j.compbiomed.2024.108760] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2023] [Revised: 06/12/2024] [Accepted: 06/13/2024] [Indexed: 07/02/2024]
Abstract
BACKGROUND Extracting phenotype-representative flow patterns and their associated numerical metrics is a bottleneck in the clinical translation of advanced cardiac flow imaging modalities. We hypothesized that reduced-order models (ROMs) are a suitable strategy for deriving simple and interpretable clinical metrics of intraventricular flow suitable for further assessments. Combined with machine learning (ML) flow-based ROMs could provide new insight to help diagnose and risk-stratify patients. METHODS We analyzed 2D color-Doppler echocardiograms of 81 non-ischemic dilated cardiomyopathy (DCM) patients, 51 hypertrophic cardiomyopathy (HCM) patients, and 77 normal volunteers (Control). We applied proper orthogonal decomposition (POD) to build patient-specific and cohort-specific ROMs of LV flow. Each ROM aggregates a low number of components representing a spatially dependent velocity map modulated along the cardiac cycle by a time-dependent coefficient. We tested three classifiers using deliberately simple ML analyses of these ROMs with varying supervision levels. In supervised models, hyperparameter grid search was used to derive the ROMs that maximize classification power. The classifiers were blinded to LV chamber geometry and function. We ran vector flow mapping on the color-Doppler sequences to help visualize flow patterns and interpret the ML results. RESULTS POD-based ROMs stably represented each cohort through 10-fold cross-validation. The principal POD mode captured >80 % of the flow kinetic energy (KE) in all cohorts and represented the LV filling/emptying jets. Mode 2 represented the diastolic vortex and its KE contribution ranged from <1 % (HCM) to 13 % (DCM). Semi-unsupervised classification using patient-specific ROMs revealed that the KE ratio of these two principal modes, the vortex-to-jet (V2J) energy ratio, is a simple, interpretable metric that discriminates DCM, HCM, and Control patients. Receiver operating characteristic curves using V2J as classifier had areas under the curve of 0.81, 0.91, and 0.95 for distinguishing HCM vs. Control, DCM vs. Control, and DCM vs. HCM, respectively. CONCLUSIONS Modal decomposition of cardiac flow can be used to create ROMs of normal and pathological flow patterns, uncovering simple interpretable flow metrics with power to discriminate disease states, and particularly suitable for further processing using ML.
Collapse
Affiliation(s)
- María Guadalupe Borja
- Department of Mechanical and Aerospace Engineering, University of California San Diego, La Jolla, CA, USA
| | - Pablo Martinez-Legazpi
- Department of Mathematical Physics and Fluids, Facultad de Ciencias, Universidad Nacional de Educación a Distancia, UNED and CIBERCV, Madrid, Spain
| | - Cathleen Nguyen
- Department of Mechanical Engineering, University of Washington, Seattle, WA, USA
| | - Oscar Flores
- Department of Aerospace Engineering, Universidad Carlos III de Madrid, Leganés, Spain
| | - Andrew M Kahn
- Division of Cardiovascular Medicine, University of California San Diego, La Jolla, CA, USA
| | - Javier Bermejo
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón and CIBERCV, Madrid, Spain
| | - Juan C Del Álamo
- Department of Mechanical and Aerospace Engineering, University of California San Diego, La Jolla, CA, USA; Department of Mechanical Engineering, University of Washington, Seattle, WA, USA; Center for Cardiovascular Biology, University of Washington, Seattle, WA, USA; Division of Cardiology, University of Washington, Seattle, WA, USA.
| |
Collapse
|
2
|
Bennati L, Puppini G, Giambruno V, Luciani GB, Vergara C. Image-Based Computational Fluid Dynamics to Compare Two Repair Techniques for Mitral Valve Prolapse. Ann Biomed Eng 2024:10.1007/s10439-024-03597-8. [PMID: 39120769 DOI: 10.1007/s10439-024-03597-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Accepted: 07/29/2024] [Indexed: 08/10/2024]
Abstract
OBJECTIVE : The treatment of mitral valve prolapse involves two distinct repair techniques: chordal replacement (Neochordae technique) and leaflet resection (Resection technique). However, there is still a debate in the literature about which is the optimal one. In this context, we performed an image-based computational fluid dynamic study to evaluate blood dynamics in the two surgical techniques. METHODS : We considered a healthy subject (H) and two patients (N and R) who underwent surgery for prolapse of the posterior leaflet and were operated with the Neochordae and Resection technique, respectively. Computational Fluid Dynamics (CFD) was employed with prescribed motion of the entire left heart coming from cine-MRI images, with a Large Eddy Simulation model to describe the transition to turbulence and a resistive method for managing valve dynamics. We created three different virtual scenarios where the operated mitral valves were inserted in the same left heart geometry of the healthy subject to study the differences attributed only to the two techniques. RESULTS : We compared the three scenarios by quantitatively analyzing ventricular velocity patterns and pressures, transition to turbulence, and the ventricle ability to prevent thrombi formation. From these results, we found that the operative techniques affected the ventricular blood dynamics in different ways, with variations attributed to the reduced mobility of the Resection posterior leaflet. Specifically, the Resection technique resulted in turbulent forces, related with the risk of hemolysis formation, up to 640 Pa, while the other two scenarios exhibited a maximum of 240 Pa. Moreover, in correspondence of the ventricular apex, the Resection technique reduced the areas with low velocity to 15%, whereas the healthy case and the Neochordae case maintained these areas at 30 and 48%, respectively. Our findings suggest that the Neochordae technique developed a more physiological flow with respect to the Resection technique. CONCLUSION: Resection technique gives rise to a different direction of the mitral jet during diastole increasing the ability to washout the ventricular apex preventing from thrombi formation, but at the same time it promotes turbulence formation that is associated with ventricular effort and risk of hemolysis.
Collapse
Affiliation(s)
- Lorenzo Bennati
- Department of Surgery, Dentistry, Pediatrics, and Obstetrics/Gynecology, University of Verona, Piazzale Ludovico Antonio Scuro 10, Verona, 37134, Italy
| | - Giovanni Puppini
- Department of Radiology, University of Verona, Piazzale Stefani 1, Verona, 37126, Italy
| | - Vincenzo Giambruno
- Division of Cardiac Surgery, Department of Surgery, Dentistry, Pediatrics, and Obstetrics/Gynecology, University of Verona, Piazzale Stefani 1, Verona, 37126, Italy
| | - Giovanni Battista Luciani
- Division of Cardiac Surgery, Department of Surgery, Dentistry, Pediatrics, and Obstetrics/Gynecology, University of Verona, Piazzale Stefani 1, Verona, 37126, Italy
| | - Christian Vergara
- LaBS, Dipartimento di Chimica, Materiali e Ingegneria Chimica "Giulio Natta", Politecnico di Milano, Piazza Leonardo da Vinci 32, Milan, 20133, Italy.
| |
Collapse
|
3
|
Rodríguez-González E, Martínez-Legazpi P, González-Mansilla A, Espinosa MÁ, Mombiela T, Guzmán De-Villoria JA, Borja MG, Díaz-Otero F, Gómez de Antonio R, Fernández-García P, Fernández-Ávila AI, Pascual-Izquierdo C, Del Álamo JC, Bermejo J. Cardiac stasis imaging, stroke, and silent brain infarcts in patients with nonischemic dilated cardiomyopathy. Am J Physiol Heart Circ Physiol 2024; 327:H446-H453. [PMID: 38847759 DOI: 10.1152/ajpheart.00245.2024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/18/2024] [Revised: 05/23/2024] [Accepted: 05/31/2024] [Indexed: 07/17/2024]
Abstract
Cardioembolic stroke is one of the most devastating complications of nonischemic dilated cardiomyopathy (NIDCM). However, in clinical trials of primary prevention, the benefits of anticoagulation are hampered by the risk of bleeding. Indices of cardiac blood stasis may account for the risk of stroke and be useful to individualize primary prevention treatments. We performed a cross-sectional study in patients with NIDCM and no history of atrial fibrillation (AF) from two sources: 1) a prospective enrollment of unselected patients with left ventricular (LV) ejection fraction <45% and 2) a retrospective identification of patients with a history of previous cardioembolic neurological event. The primary end point integrated a history of ischemic stroke or the presence intraventricular thrombus, or a silent brain infarction (SBI) by imaging. From echocardiography, we calculated blood flow inside the LV, its residence time (TR) maps, and its derived stasis indices. Of the 89 recruited patients, 18 showed a positive end point, 9 had a history of stroke or transient ischemic attack (TIA) and 9 were diagnosed with SBIs in the brain imaging. Averaged TR, [Formula: see text] performed well to identify the primary end point [AUC (95% CI) = 0.75 (0.61-0.89), P = 0.001]. When accounting only for identifying a history of stroke or TIA, AUC for [Formula: see text] was 0.92 (0.85-1.00) with odds ratio = 7.2 (2.3-22.3) per cycle, P < 0.001. These results suggest that in patients with NIDCM in sinus rhythm, stasis imaging derived from echocardiography may account for the burden of stroke.NEW & NOTEWORTHY Patients with nonischemic dilated cardiomyopathy (NIDCM) are at higher risk of stroke than their age-matched population. However, the risk of bleeding neutralizes the benefit of preventive oral anticoagulation. In this work, we show that in patients in sinus rhythm, the burden of stroke is related to intraventricular stasis metrics derived from echocardiography. Therefore, stasis metrics may be useful to personalize primary prevention anticoagulation in these patients.
Collapse
Affiliation(s)
- Elena Rodríguez-González
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Madrid, Spain
- Department of Medicine, Facultad de Medicina, Universidad Complutense, Madrid, Spain
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares, Madrid, Spain
| | - Pablo Martínez-Legazpi
- Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares, Madrid, Spain
- Department of Mathematical Physics and Fluids, Facultad de Ciencias, Universidad Nacional de Educación a Distancia, Madrid, Spain
| | - Ana González-Mansilla
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Madrid, Spain
- Department of Medicine, Facultad de Medicina, Universidad Complutense, Madrid, Spain
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares, Madrid, Spain
| | - M Ángeles Espinosa
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Madrid, Spain
- Department of Medicine, Facultad de Medicina, Universidad Complutense, Madrid, Spain
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares, Madrid, Spain
| | - Teresa Mombiela
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Madrid, Spain
- Department of Medicine, Facultad de Medicina, Universidad Complutense, Madrid, Spain
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares, Madrid, Spain
| | - Juan A Guzmán De-Villoria
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Department of Radiology, Hospital General Universitario Gregorio Marañón, Madrid, Spain
- Centro de Investigación Biomédica en Red de Salud Mental, Madrid, Spain
| | - Maria Guadalupe Borja
- Department of Mechanical and Aerospace Engineering, University of California San Diego, La Jolla, California, United States
| | - Fernando Díaz-Otero
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Department of Neurology, Hospital General Universitario Gregorio Marañón, Madrid, Spain
| | - Rubén Gómez de Antonio
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Department of Hematology, Hospital General Universitario Gregorio Marañón, Madrid, Spain
| | - Pilar Fernández-García
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Centro de Investigación Biomédica en Red de Salud Mental, Madrid, Spain
| | - Ana I Fernández-Ávila
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Madrid, Spain
- Department of Medicine, Facultad de Medicina, Universidad Complutense, Madrid, Spain
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares, Madrid, Spain
| | - Cristina Pascual-Izquierdo
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Department of Hematology, Hospital General Universitario Gregorio Marañón, Madrid, Spain
| | - Juan C Del Álamo
- Division of Cardiology, Department of Mechanical Engineering, Center for Cardiovascular Biology, University of Washington, Seattle, Washington, United States
| | - Javier Bermejo
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Madrid, Spain
- Department of Medicine, Facultad de Medicina, Universidad Complutense, Madrid, Spain
- Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares, Madrid, Spain
| |
Collapse
|
4
|
Bennati L, Giambruno V, Renzi F, Di Nicola V, Maffeis C, Puppini G, Luciani GB, Vergara C. Turbulent blood dynamics in the left heart in the presence of mitral regurgitation: a computational study based on multi-series cine-MRI. Biomech Model Mechanobiol 2023; 22:1829-1846. [PMID: 37400622 PMCID: PMC10613156 DOI: 10.1007/s10237-023-01735-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2023] [Accepted: 06/04/2023] [Indexed: 07/05/2023]
Abstract
In this work, we performed a computational image-based study of blood dynamics in the whole left heart, both in a healthy subject and in a patient with mitral valve regurgitation. We elaborated multi-series cine-MRI with the aim of reconstructing the geometry and the corresponding motion of left ventricle, left atrium, mitral and aortic valves, and aortic root of the subjects. This allowed us to prescribe such motion to computational blood dynamics simulations where, for the first time, the whole left heart motion of the subject is considered, allowing us to obtain reliable subject-specific information. The final aim is to investigate and compare between the subjects the occurrence of turbulence and the risk of hemolysis and of thrombi formation. In particular, we modeled blood with the Navier-Stokes equations in the arbitrary Lagrangian-Eulerian framework, with a large eddy simulation model to describe the transition to turbulence and a resistive method to manage the valve dynamics, and we used a finite element discretization implemented in an in-house code for the numerical solution.
Collapse
Affiliation(s)
- Lorenzo Bennati
- Department of Surgery, Dentistry, Pediatrics, and Obstetrics/Gynecology, University of Verona, Piazzale Ludovico Antonio Scuro 10, 37134, Verona, Italy
| | - Vincenzo Giambruno
- Division of Cardiac Surgery, Department of Surgery, Dentistry, Pediatrics, and Obstetrics/Gynecology, University of Verona, Piazzale Stefani 1, 37126, Verona, Italy
| | - Francesca Renzi
- Department of Surgery, Dentistry, Pediatrics, and Obstetrics/Gynecology, University of Verona, Piazzale Ludovico Antonio Scuro 10, 37134, Verona, Italy
| | - Venanzio Di Nicola
- Division of Cardiac Surgery, Department of Surgery, Dentistry, Pediatrics, and Obstetrics/Gynecology, University of Verona, Piazzale Stefani 1, 37126, Verona, Italy
| | - Caterina Maffeis
- Department of Surgery, Dentistry, Pediatrics, and Obstetrics/Gynecology, University of Verona, Piazzale Ludovico Antonio Scuro 10, 37134, Verona, Italy
| | - Giovanni Puppini
- Department of Radiology, University of Verona, Piazzale Stefani 1, 37126, Verona, Italy
| | - Giovanni Battista Luciani
- Division of Cardiac Surgery, Department of Surgery, Dentistry, Pediatrics, and Obstetrics/Gynecology, University of Verona, Piazzale Stefani 1, 37126, Verona, Italy
| | - Christian Vergara
- LaBS, Dipartimento di Chimica, Materiali e Ingegneria Chimica "Giulio Natta", Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milan, Italy.
| |
Collapse
|
5
|
Postigo A, Viola F, Chazo C, Martínez-Legazpi P, González-Mansilla A, Rodríguez-González E, Fernández-Avilés F, Del Álamo JC, Ebbers T, Bermejo J. Assessment of Blood Flow Transport in the Left Ventricle Using Ultrasound. Validation Against 4-D Flow Cardiac Magnetic Resonance. ULTRASOUND IN MEDICINE & BIOLOGY 2022; 48:1822-1832. [PMID: 35764455 PMCID: PMC10408642 DOI: 10.1016/j.ultrasmedbio.2022.05.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Revised: 05/04/2022] [Accepted: 05/08/2022] [Indexed: 06/15/2023]
Abstract
Four-dimensional flow cardiac magnetic resonance (CMR) is the reference technique for analyzing blood transport in the left ventricle (LV), but similar information can be obtained from ultrasound. We aimed to validate ultrasound-derived transport in a head-to-head comparison against 4D flow CMR. In five patients and two healthy volunteers, we obtained 2D + t and 3D + t (4D) flow fields in the LV using transthoracic echocardiography and CMR, respectively. We compartmentalized intraventricular blood flow into four fractions of end-diastolic volume: direct flow (DF), retained inflow (RI), delayed ejection flow (DEF) and residual volume (RV). Using ultrasound we also computed the properties of LV filling waves (percentage of LV penetration and percentage of LV volume carried by E/A waves) to determine their relationships with CMR transport. Agreement between both techniques for quantifying transport fractions was good for DF and RV (Ric [95% confidence interval]: 0.82 [0.33, 0.97] and 0.85 [0.41, 0.97], respectively) and moderate for RI and DEF (Ric= 0.47 [-0.29, 0.88] and 0.55 [-0.20, 0.90], respectively). Agreement between techniques to measure kinetic energy was variable. The amount of blood carried by the E-wave correlated with DF and RV (R = 0.75 and R = 0.63, respectively). Therefore, ultrasound is a suitable method for expanding the analysis of intraventricular flow transport in the clinical setting.
Collapse
Affiliation(s)
- Andrea Postigo
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón, CIBERCV, Madrid, Spain
| | - Federica Viola
- Department of Health, Medicine and Caring Sciences and Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden
| | - Christian Chazo
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón, CIBERCV, Madrid, Spain
| | - Pablo Martínez-Legazpi
- Department of Mathematical Physics and Fluids, Facultad de Ciencias, Universidad Nacional de Educación a Distancia, UNED and CIBERCV, Madrid, Spain
| | - Ana González-Mansilla
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón, CIBERCV, Madrid, Spain
| | - Elena Rodríguez-González
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón, CIBERCV, Madrid, Spain
| | - Francisco Fernández-Avilés
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón, CIBERCV, Madrid, Spain
| | - Juan C Del Álamo
- Mechanical Engineering Department, Center for Cardiovascular Biology, Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, Washington, USA
| | - Tino Ebbers
- Department of Health, Medicine and Caring Sciences and Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden
| | - Javier Bermejo
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón, CIBERCV, Madrid, Spain.
| |
Collapse
|
6
|
Sakakibara T, Suwa K, Ushio T, Wakayama T, Alley M, Saotome M, Satoh H, Maekawa Y. Intra-Left Ventricular Hemodynamics Assessed with 4D Flow Magnetic Resonance Imaging in Patients with Left Ventricular Thrombus. Int Heart J 2021; 62:1287-1296. [PMID: 34853222 DOI: 10.1536/ihj.20-792] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Left ventricular thrombus (LVT) has been identified to be crucial in patients with reduced ejection fraction (EF). Three-dimensional cine phase-contrast magnetic resonance imaging (4D flow MRI) can visualize the intra-LV vortex during diastole and quantify the maximum flow velocity (Vmax) at the apex. In this study, we investigated whether the change in the intra-LV vortex was associated with the presence of LVT in patients with cardiac disease.In total, 36 patients (63.5 ± 11.9 years, 28 men, 12/24 with/without LVT) with diffuse LV dysfunction underwent 4D flow MRI. The relative vortex area using streamline images and Vmax of blood flow toward the apex at the apical left ventricle were evaluated. The correlation between the relative vortex area and Vmax was assessed using Pearson's correlation coefficient. The ability to detect LVT was evaluated using the area under the curve (AUC) of the receiver operating characteristic.The relative vortex area was found to be smaller (27 ± 10% versus 45 ± 11%, P = 0.000026), whereas Vmax at the apical left ventricle was lower (19.1 ± 4.4 cm/second versus 27.4 ± 8.9 cm/second, P = 0.0006) in patients with LVT. Vmax at the apical left ventricle demonstrated significant correlations with the relative vortex area (r = 0.43, P = 0.01) and relative transverse length of the vortex (r = 0.45, P = 0.007). The AUC was 0.91 for the relative vortex area, whereas it was 0.80 for Vmax in the apical left ventricle.A smaller LV vortex and lower flow velocity at the LV apex were associated with LVT in patients with reduced EF.
Collapse
Affiliation(s)
- Tomoaki Sakakibara
- Division of Cardiology, Internal Medicine 3, Hamamatsu University School of Medicine
| | - Kenichiro Suwa
- Division of Cardiology, Internal Medicine 3, Hamamatsu University School of Medicine
| | - Takasuke Ushio
- Department of Radiology, Hamamatsu University School of Medicine
| | | | - Marcus Alley
- Division of Radiology, Stanford University School of Medicine, Stanford
| | - Masao Saotome
- Division of Cardiology, Internal Medicine 3, Hamamatsu University School of Medicine
| | - Hiroshi Satoh
- Department of Cardiology, Fujinomiya City General Hospital
| | - Yuichiro Maekawa
- Division of Cardiology, Internal Medicine 3, Hamamatsu University School of Medicine
| |
Collapse
|
7
|
Wojtas K, Kozłowski M, Orciuch W, Makowski Ł. Computational Fluid Dynamics Simulations of Mitral Paravalvular Leaks in Human Heart. MATERIALS 2021; 14:ma14237354. [PMID: 34885504 PMCID: PMC8658524 DOI: 10.3390/ma14237354] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Revised: 11/21/2021] [Accepted: 11/27/2021] [Indexed: 11/16/2022]
Abstract
In recent years, computational fluid dynamics (CFD) has been extensively used in biomedical research on heart diseases due to its non-invasiveness and relative ease of use in predicting flow patterns inside the cardiovascular system. In this study, a modeling approach involving CFD simulations was employed to study hemodynamics inside the left ventricle (LV) of a human heart affected by a mitral paravalvular leak (PVL). A simplified LV geometry with four PVL variants that varied in shape and size was studied. Predicted blood flow parameters, mainly velocity and shear stress distributions, were used as indicators of how presence of PVLs correlates with risk and severity of hemolysis. The calculations performed in the study showed a high risk of hemolysis in all analyzed cases, with the maximum shear stress values considerably exceeding the safe level of 300 Pa. Results of our study indicated that there was no simple relationship between PVL geometry and the risk of hemolysis. Two factors that potentially played a role in hemolysis severity, namely erythrocyte exposure time and the volume of fluid in which shear stress exceeded a critical value, were not directly proportional to any of the characteristic geometrical parameters (shape, diameters, circumference, area, volume) of the PVL channel. Potential limitations of the proposed simplified approach of flow analysis are discussed, and possible modifications to increase the accuracy and plausibility of the results are presented.
Collapse
Affiliation(s)
- Krzysztof Wojtas
- Faculty of Chemical and Process Engineering, Warsaw University of Technology, Warynskiego 1, 00-645 Warsaw, Poland; (W.O.); (Ł.M.)
- Correspondence: ; Tel.: +48-22-234-6275
| | - Michał Kozłowski
- Department of Cardiology and Structural Heart Diseases, Medical University of Silesia, Ziolowa 45/47, 40-635 Katowice, Poland;
| | - Wojciech Orciuch
- Faculty of Chemical and Process Engineering, Warsaw University of Technology, Warynskiego 1, 00-645 Warsaw, Poland; (W.O.); (Ł.M.)
| | - Łukasz Makowski
- Faculty of Chemical and Process Engineering, Warsaw University of Technology, Warynskiego 1, 00-645 Warsaw, Poland; (W.O.); (Ł.M.)
| |
Collapse
|
8
|
Pedrizzetti G, Faganello G, Croatto E, Di Lenarda A. The hemodynamic power of the heart differentiates normal from diseased right ventricles. J Biomech 2021; 119:110312. [PMID: 33609983 DOI: 10.1016/j.jbiomech.2021.110312] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2020] [Revised: 01/31/2021] [Accepted: 02/03/2021] [Indexed: 12/15/2022]
Abstract
Cardiac mechanics is primarily described by the pressure-volume relationship. The ventricular pressure-volume loop displays the instantaneous relationship between intraventricular pressure and volume throughout the cardiac cycle; however, it does not consider the shape of the ventricles, their spatiotemporal deformation patterns, and how these balance with the flowing blood. Our study demonstrates that the pressure-volume relationship represents a first level of approximation for the mechanical power of the ventricles, while, at a further level of approximation, the importance of hemodynamic power emerges through the balance between deformation patterns and fluid dynamics. The analysis is preliminarily tested in a healthy subject's right ventricle and two patients. Moreover, patients' geometry was then rescaled to present a normal volumetric profile to verify whether results were affected by volume size or by the spatiotemporal pattern of how that volume profile was achieved. Results show that alterations of hemodynamic power were found in the abnormal ventricles and that they were not directly caused by the ventricular size but by changes in the ability of intraventricular pressure gradient to generate blood flow. Therefore, hemodynamic power represents a physics-based measure that takes into account the dynamics of the space-time shape changes in combination with blood flow. Hemodynamic power is assessed non-invasively using cardiac imaging techniques and can be an early indicator of cardiac dysfunction before changes occur in volumetric measurements. These preliminary results provide a physical ground to evaluate its diagnostic or prognostic significance in future clinical studies.
Collapse
Affiliation(s)
- Gianni Pedrizzetti
- Department of Engineering and Architecture, University of Trieste, Italy.
| | - Giorgio Faganello
- Cardiovascular Department, Azienda Sanitaria Universitaria Giuliano Isontina, Trieste, Italy
| | - Elisa Croatto
- Cardiovascular Department, Azienda Sanitaria Universitaria Giuliano Isontina, Trieste, Italy
| | - Andrea Di Lenarda
- Cardiovascular Department, Azienda Sanitaria Universitaria Giuliano Isontina, Trieste, Italy
| |
Collapse
|
9
|
Reindl M, Reinstadler SJ. E-wave propagation index (EPI) - A promising echocardiographic marker to improve left ventricular thrombus detection after STEMI? Int J Cardiol 2021; 331:331-332. [PMID: 33508340 DOI: 10.1016/j.ijcard.2021.01.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/09/2021] [Accepted: 01/21/2021] [Indexed: 10/22/2022]
Affiliation(s)
- Martin Reindl
- University Clinic of Internal Medicine III, Cardiology and Angiology, Medical University of Innsbruck, Anichstrasse 35, A-6020 Innsbruck, Austria
| | - Sebastian J Reinstadler
- University Clinic of Internal Medicine III, Cardiology and Angiology, Medical University of Innsbruck, Anichstrasse 35, A-6020 Innsbruck, Austria.
| |
Collapse
|
10
|
Association between the E-wave propagation index and left ventricular thrombus formation after ST-elevation myocardial infarction. Int J Cardiol 2020; 326:213-219. [PMID: 33152416 DOI: 10.1016/j.ijcard.2020.10.078] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/31/2020] [Revised: 10/24/2020] [Accepted: 10/28/2020] [Indexed: 11/20/2022]
Abstract
OBJECTIVE To explore the association between E-wave propagation index (EPI) as a marker of apical washout and the risk of left ventricular thrombus (LVT) formation in patients with ST-elevation myocardial infarction (STEMI). METHODS We performed a post-hoc analysis on 364 prospectively enrolled STEMI patients from a single-center. Non-contrast transthoracic echocardiographic examinations were performed a median of 2 days (IQR:1-3 days) after PCI. The endpoint was LVT formation, identified retrospectively. Univariable and multivariable logistic regression was applied to assess the association between EPI and LVT formation. Multivariable adjustments included LVEF, LAD culprit, prior myocardial infarction, heart rate, and early myocardial relaxation velocity. Area under receiver operating characteristic curves (AUC) was used to assess the diagnostic ability. RESULTS AND CONCLUSIONS Among 364 patients, 31 (8.5%) developed LVT. The mean age was 62 years, 75% were men, and mean LVEF was 46%. Patients developing LVT had increased heart rate, lower LVEF, impaired GLS, and more frequently had prior myocardial infarction. Variables associated with low values of EPI included, among others, LVEF, LV aneurysm, and GLS. EPI and LVT formation were significantly associated in the univariable model (OR = 1.87 (1.53-2.28), p < 0.001), and EPI showed an AUC of 0.90. After multivariable adjustments, EPI and LVT formation remained significantly associated (OR = 1.79 (1.42-2.27), p < 0.001). Patients with an EPI < 1.0 had a 23 times higher likelihood of LVT formation (OR = 23.41 (10.06-54.49), p < 0.001). EPI and LVT formation are strongly associated in patients with STEMI, with low values of EPI indicating a markedly increased probability of LVT formation.
Collapse
|
11
|
Nguyen TD, Kadri OE, Voronov RS. An Introductory Overview of Image-Based Computational Modeling in Personalized Cardiovascular Medicine. Front Bioeng Biotechnol 2020; 8:529365. [PMID: 33102452 PMCID: PMC7546862 DOI: 10.3389/fbioe.2020.529365] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2020] [Accepted: 08/31/2020] [Indexed: 02/05/2023] Open
Abstract
Cardiovascular diseases account for the number one cause of deaths in the world. Part of the reason for such grim statistics is our limited understanding of the underlying mechanisms causing these devastating pathologies, which is made difficult by the invasiveness of the procedures associated with their diagnosis (e.g., inserting catheters into the coronal artery to measure blood flow to the heart). Likewise, it is also difficult to design and test assistive devices without implanting them in vivo. However, with the recent advancements made in biomedical scanning technologies and computer simulations, image-based modeling (IBM) has arisen as the next logical step in the evolution of non-invasive patient-specific cardiovascular medicine. Yet, due to its novelty, it is still relatively unknown outside of the niche field. Therefore, the goal of this manuscript is to review the current state-of-the-art and the limitations of the methods used in this area of research, as well as their applications to personalized cardiovascular investigations and treatments. Specifically, the modeling of three different physics – electrophysiology, biomechanics and hemodynamics – used in the cardiovascular IBM is discussed in the context of the physiology that each one of them describes and the mechanisms of the underlying cardiac diseases that they can provide insight into. Only the “bare-bones” of the modeling approaches are discussed in order to make this introductory material more accessible to an outside observer. Additionally, the imaging methods, the aspects of the unique cardiac anatomy derived from them, and their relation to the modeling algorithms are reviewed. Finally, conclusions are drawn about the future evolution of these methods and their potential toward revolutionizing the non-invasive diagnosis, virtual design of treatments/assistive devices, and increasing our understanding of these lethal cardiovascular diseases.
Collapse
Affiliation(s)
- Thanh Danh Nguyen
- Otto H. York Department of Chemical and Materials Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States
| | - Olufemi E Kadri
- Otto H. York Department of Chemical and Materials Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States.,UC-P&G Simulation Center, University of Cincinnati, Cincinnati, OH, United States
| | - Roman S Voronov
- Otto H. York Department of Chemical and Materials Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States.,Department of Biomedical Engineering, Newark College of Engineering, New Jersey Institute of Technology, Newark, NJ, United States
| |
Collapse
|
12
|
Delgado-Montero A, Martinez-Legazpi P, Desco MM, Rodríguez-Pérez D, Díaz-Otero F, Rossini L, Pérez del Villar C, Rodríguez-González E, Chazo C, Benito Y, Flores O, Antoranz JC, Fernández-Avilés F, del Álamo JC, Bermejo J. Blood Stasis Imaging Predicts Cerebral Microembolism during Acute Myocardial Infarction. J Am Soc Echocardiogr 2020; 33:389-398. [DOI: 10.1016/j.echo.2019.09.020] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/07/2019] [Revised: 09/25/2019] [Accepted: 09/25/2019] [Indexed: 10/25/2022]
|
13
|
Collia D, Zovatto L, Pedrizzetti G. Analysis of mitral valve regurgitation by computational fluid dynamics. APL Bioeng 2019; 3:036105. [PMID: 31893254 PMCID: PMC6932856 DOI: 10.1063/1.5097245] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2019] [Accepted: 07/31/2019] [Indexed: 02/05/2023] Open
Abstract
The clinical syndrome of mitral insufficiency is a common consequence of mitral valve (MV) prolapse, when the MV leaflets do not seal the closed orifice and blood regurgitates back to the atrium during ventricular contraction. There are different types of MV prolapse that may influence the degree of regurgitation also in relation to the left ventricle (LV) geometry. This study aims to provide some insight into the fluid dynamics of MV insufficiency in view of improving the different measurements available in the clinical setting. The analysis is performed by a computational fluid dynamics model coupled with an asymptotic model of the MV motion. The computational fluid dynamics solution uses the immersed boundary method that is efficiently integrated with clinical imaging technologies. Healthy and dilated LVs obtained by multislice cardiac MRI are combined with simplified models of healthy and pathological MVs deduced from computed tomography and 4D-transesophageal echocardiography. The results demonstrated the properties of false regurgitation, blood that did not cross the open MV orifice and returns into the atrium during the backward motion of the MV leaflets, whose entity should be accounted when evaluating small regurgitation. The regurgitating volume is found to be proportional to the effective orifice area, with the limited dependence of the LV geometry and type of prolapse. These affect the percentage of old blood returning to the atrium which may be associated with thrombogenic risk.
Collapse
Affiliation(s)
- Dario Collia
- Department of Engineering and Architecture, University of Trieste, P. Europa 1, 34127 Trieste, Italy
| | - Luigino Zovatto
- Department of Engineering and Architecture, University of Trieste, P. Europa 1, 34127 Trieste, Italy
| | - Gianni Pedrizzetti
- Department of Engineering and Architecture, University of Trieste, P. Europa 1, 34127 Trieste, Italy
| |
Collapse
|
14
|
Benito Y, Martinez-Legazpi P, Rossini L, Pérez del Villar C, Yotti R, Martín Peinador Y, Rodríguez-Pérez D, Desco MM, Medrano C, Antoranz JC, Fernández-Avilés F, del Álamo JC, Bermejo J. Age-Dependence of Flow Homeostasis in the Left Ventricle. Front Physiol 2019; 10:485. [PMID: 31105588 PMCID: PMC6498893 DOI: 10.3389/fphys.2019.00485] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2019] [Accepted: 04/08/2019] [Indexed: 11/29/2022] Open
Abstract
Background: Intracardiac flow homeostasis requires avoiding blood stasis and platelet activation during its transit through the cardiac chambers. However, the foundations of intraventricular blood washout and its exposure to shear stresses have been poorly addressed. We aimed to characterize and quantify these features in a wide population of healthy subjects and assess the relationships of these indices with age. Methods: We used color-Doppler echocardiography and custom post-processing methods to study 149 healthy volunteers from 26 days to 80 years old. From the intraventricular flow-velocity fields we obtained personalized maps of (1) the residence time of blood in the LV, and (2) the shear index, a metric accounting for the strongest occurrence of shear stresses inside the chamber. From these maps we derived quantitative indices of the overall intraventricular blood washout and shear exposure. We addressed the age-dependence of these indices and analyzed their relationship with age-related changes in filling-flow. Results: The entire intraventricular blood pool was replaced before 8 cycles. Average residence time of blood inside the LV was <3 cycles in all subjects and followed an inverse U-shape relationship with age, increasing from median (IQR) of 1.0 (0.7 to 1.2) cycles in the 1st year of life to 1.8 (1.4-2.2) cycles in young adults (17-30 years old), becoming shorter again thereafter. Shear index showed no relation with age and was bounded around 20 dyn·s/cm2. Regions with the longest residence time and highest shear index were identified near the apex. Differences in the degree of apical penetration of the filling waves and the duration of the late-filling phase explained the age-dependence of residence time (R adj 2 = 0.48, p < 0.001). Conclusions: In average, blood spends 1 to 3 beats inside the LV with very low shear stress rates. The apical region is the most prone to blood stasis, particularly in mid-aged adults. The washout of blood in the normal LV is age-dependent due to physiological changes in the degree of apical penetration of the filling waves.
Collapse
Affiliation(s)
- Yolanda Benito
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón and CIBERCV, Madrid, Spain
| | - Pablo Martinez-Legazpi
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón and CIBERCV, Madrid, Spain
| | - Lorenzo Rossini
- Department of Mechanical and Aerospace Engineering, University of California, San Diego, La Jolla, CA, United States
| | - Candelas Pérez del Villar
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón and CIBERCV, Madrid, Spain
| | - Raquel Yotti
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón and CIBERCV, Madrid, Spain
| | - Yolanda Martín Peinador
- Centro de Salud Goya, Dirección Asistencial Centro, Atención Primaria de Madrid, Madrid, Spain
| | - Daniel Rodríguez-Pérez
- Department of Mathematical Physics and Fluids, Facultad de Ciencias, Universidad Nacional de Educación a Distancia, Madrid, Spain
| | - M. Mar Desco
- Department of Mathematical Physics and Fluids, Facultad de Ciencias, Universidad Nacional de Educación a Distancia, Madrid, Spain
| | - Constancio Medrano
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón and CIBERCV, Madrid, Spain
| | - Jose Carlos Antoranz
- Department of Mathematical Physics and Fluids, Facultad de Ciencias, Universidad Nacional de Educación a Distancia, Madrid, Spain
| | - Francisco Fernández-Avilés
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón and CIBERCV, Madrid, Spain
| | - Juan C. del Álamo
- Department of Mechanical and Aerospace Engineering, University of California, San Diego, La Jolla, CA, United States
- Institute for Engineering in Medicine, University of California, San Diego, La Jolla, CA, United States
| | - Javier Bermejo
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Facultad de Medicina, Universidad Complutense de Madrid, Instituto de Investigación Sanitaria Gregorio Marañón and CIBERCV, Madrid, Spain
| |
Collapse
|
15
|
Sampath K, Harfi TT, George RT, Katz J. Optimized Time-Resolved Echo Particle Image Velocimetry– Particle Tracking Velocimetry Measurements Elucidate Blood Flow in Patients With Left Ventricular Thrombus. J Biomech Eng 2018; 140:2668583. [DOI: 10.1115/1.4038886] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2017] [Indexed: 02/04/2023]
Abstract
Contrast ultrasound is a widely used clinical tool to obtain real-time qualitative blood flow assessments in the heart, liver, etc. Echocardiographic particle image velocimetry (echo-PIV) is a technique for obtaining quantitative velocity maps from contrast ultrasound images. However, unlike optical particle image velocimetry (PIV), routine echo images are prone to nonuniform spatiotemporal variations in tracer distribution, making analysis difficult for standard PIV algorithms. This study introduces optimized procedures that integrate image enhancement, PIV, and particle tracking velocimetry (PTV) to obtain reliable time-resolved two-dimensional (2D) velocity distributions. During initial PIV analysis, multiple results are obtained by varying processing parameters. Optimization involving outlier removal and smoothing is used to select the correct vector. These results are used in a multiparameter PTV procedure. To demonstrate their clinical value, the procedures are implemented to obtain velocity and vorticity distributions over multiple cardiac cycles using images acquired from four left ventricular thrombus (LVT) patients. Phase-averaged data elucidate flow structure evolution over the cycle and are used to calculate penetration depth and strength of left ventricular (LV) vortices, as well as apical velocity induced by them. The present data are consistent with previous time-averaged results for the minimum vortex penetration depth associated with LVT occurrence. However, due to decay and fragmentation of LV vortices, as they migrate away from the mitral annulus, in two cases with high penetration, there is still poor washing near the resolved clot throughout the cycle. Hence, direct examination of entire flow evolution may be useful for assessing risk of LVT relapse before prescribing anticoagulants.
Collapse
Affiliation(s)
- Kaushik Sampath
- Department of Mechanical Engineering, Johns Hopkins University, 3400 North Charles Street, Latrobe 223, Baltimore, MD 21218 e-mail:
| | - Thura T. Harfi
- Division of Cardiology, Department of Medicine, Johns Hopkins University, 600 North Wolfe Street, Baltimore, MD 21287 e-mail:
| | - Richard T. George
- Division of Cardiology, Department of Medicine, Johns Hopkins University, 600 North Wolfe Street, Baltimore, MD 21287 e-mail:
| | - Joseph Katz
- Department of Mechanical Engineering, Johns Hopkins University, 3400 North Charles Street, Latrobe 122, Baltimore, MD 21218 e-mail:
| |
Collapse
|
16
|
Vedula V, Lee J, Xu H, Kuo CCJ, Hsiai TK, Marsden AL. A method to quantify mechanobiologic forces during zebrafish cardiac development using 4-D light sheet imaging and computational modeling. PLoS Comput Biol 2017; 13:e1005828. [PMID: 29084212 PMCID: PMC5679653 DOI: 10.1371/journal.pcbi.1005828] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2017] [Revised: 11/09/2017] [Accepted: 10/15/2017] [Indexed: 01/09/2023] Open
Abstract
Blood flow and mechanical forces in the ventricle are implicated in cardiac development and trabeculation. However, the mechanisms of mechanotransduction remain elusive. This is due in part to the challenges associated with accurately quantifying mechanical forces in the developing heart. We present a novel computational framework to simulate cardiac hemodynamics in developing zebrafish embryos by coupling 4-D light sheet imaging with a stabilized finite element flow solver, and extract time-dependent mechanical stimuli data. We employ deformable image registration methods to segment the motion of the ventricle from high resolution 4-D light sheet image data. This results in a robust and efficient workflow, as segmentation need only be performed at one cardiac phase, while wall position in the other cardiac phases is found by image registration. Ventricular hemodynamics are then quantified by numerically solving the Navier-Stokes equations in the moving wall domain with our validated flow solver. We demonstrate the applicability of the workflow in wild type zebrafish and three treated fish types that disrupt trabeculation: (a) chemical treatment using AG1478, an ErbB2 signaling inhibitor that inhibits proliferation and differentiation of cardiac trabeculation; (b) injection of gata1a morpholino oligomer (gata1aMO) suppressing hematopoiesis and resulting in attenuated trabeculation; (c) weak-atriumm58 mutant (wea) with inhibited atrial contraction leading to a highly undeveloped ventricle and poor cardiac function. Our simulations reveal elevated wall shear stress (WSS) in wild type and AG1478 compared to gata1aMO and wea. High oscillatory shear index (OSI) in the grooves between trabeculae, compared to lower values on the ridges, in the wild type suggest oscillatory forces as a possible regulatory mechanism of cardiac trabeculation development. The framework has broad applicability for future cardiac developmental studies focused on quantitatively investigating the role of hemodynamic forces and mechanotransduction during morphogenesis.
Collapse
Affiliation(s)
- Vijay Vedula
- Department of Pediatrics (Cardiology), Stanford University, Stanford, California, United States of America
| | - Juhyun Lee
- Department of Bioengineering, University of California, Los Angeles, Los Angeles, California, United States of America
| | - Hao Xu
- Department of Electrical Engineering, University of Southern California, Los Angeles, California, United States of America
| | - C.-C. Jay Kuo
- Department of Electrical Engineering, University of Southern California, Los Angeles, California, United States of America
| | - Tzung K. Hsiai
- Department of Bioengineering, University of California, Los Angeles, Los Angeles, California, United States of America
- Department of Medicine, Division of Cardiology, University of California, Los Angeles, Los Angeles, California, United States of America
| | - Alison L. Marsden
- Department of Pediatrics (Cardiology), Stanford University, Stanford, California, United States of America
- Department of Bioengineering, Stanford University, Stanford, California, United States of America
- Institute for Computational and Mathematical Engineering (ICME), Stanford University, Stanford, California, United States of America
| |
Collapse
|