1
|
Mukherjee T, Keshavarzian M, Fugate EM, Naeini V, Darwish A, Ohayon J, Myers KJ, Shah DJ, Lindquist D, Sadayappan S, Pettigrew RI, Avazmohammadi R. Complete spatiotemporal quantification of cardiac motion in mice through enhanced acquisition and super-resolution reconstruction. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.31.596322. [PMID: 38895261 PMCID: PMC11185553 DOI: 10.1101/2024.05.31.596322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
The quantification of cardiac motion using cardiac magnetic resonance imaging (CMR) has shown promise as an early-stage marker for cardiovascular diseases. Despite the growing popularity of CMR-based myocardial strain calculations, measures of complete spatiotemporal strains (i.e., three-dimensional strains over the cardiac cycle) remain elusive. Complete spatiotemporal strain calculations are primarily hampered by poor spatial resolution, with the rapid motion of the cardiac wall also challenging the reproducibility of such strains. We hypothesize that a super-resolution reconstruction (SRR) framework that leverages combined image acquisitions at multiple orientations will enhance the reproducibility of complete spatiotemporal strain estimation. Two sets of CMR acquisitions were obtained for five wild-type mice, combining short-axis scans with radial and orthogonal long-axis scans. Super-resolution reconstruction, integrated with tissue classification, was performed to generate full four-dimensional (4D) images. The resulting enhanced and full 4D images enabled complete quantification of the motion in terms of 4D myocardial strains. Additionally, the effects of SRR in improving accurate strain measurements were evaluated using an in-silico heart phantom. The SRR framework revealed near isotropic spatial resolution, high structural similarity, and minimal loss of contrast, which led to overall improvements in strain accuracy. In essence, a comprehensive methodology was generated to quantify complete and reproducible myocardial deformation, aiding in the much-needed standardization of complete spatiotemporal strain calculations.
Collapse
Affiliation(s)
- Tanmay Mukherjee
- Department of Biomedical Engineering, Texas A&M University, College Station, TX 77843, USA
| | - Maziyar Keshavarzian
- Department of Biomedical Engineering, Texas A&M University, College Station, TX 77843, USA
| | - Elizabeth M. Fugate
- Department of Radiology, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH 45229, USA
| | - Vahid Naeini
- Department of Biomedical Engineering, Texas A&M University, College Station, TX 77843, USA
| | - Amr Darwish
- Houston Methodist DeBakey Heart & Vascular Center, Houston, TX 77030, USA
| | - Jacques Ohayon
- Savoie Mont-Blanc University, Polytech Annecy-Chambéry, Le Bourget du Lac, France
- Laboratory TIMC-CNRS, UMR 5525, Grenoble-Alpes University, Grenoble, France
| | - Kyle J. Myers
- Hagler Institute for Advanced Study, Texas A&M University, College Station, TX 77843, USA
| | - Dipan J. Shah
- Houston Methodist DeBakey Heart & Vascular Center, Houston, TX 77030, USA
| | - Diana Lindquist
- Department of Radiology, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH 45229, USA
| | - Sakthivel Sadayappan
- Department of Internal Medicine, Division of Cardiovascular Health and Disease, University of Cincinnati College of Medicine, Cincinnati, OH 45267, USA
| | - Roderic I. Pettigrew
- School of Engineering Medicine, Texas AM University, Houston, TX 77030, USA
- Department of Cardiovascular Sciences, Houston Methodist Academic Institute, Houston, TX 77030, USA
| | - Reza Avazmohammadi
- Department of Biomedical Engineering, Texas A&M University, College Station, TX 77843, USA
- Department of Cardiovascular Sciences, Houston Methodist Academic Institute, Houston, TX 77030, USA
- J. Mike Walker ’66 Department of Mechanical Engineering, Texas A&M University, College Station, TX 77843, USA
| |
Collapse
|
2
|
Jain S, Belkadi H, Michaut A, Sart S, Gros J, Genet M, Baroud CN. Using a micro-device with a deformable ceiling to probe stiffness heterogeneities within 3D cell aggregates. Biofabrication 2024; 16:035010. [PMID: 38447213 DOI: 10.1088/1758-5090/ad30c7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 03/06/2024] [Indexed: 03/08/2024]
Abstract
Recent advances in the field of mechanobiology have led to the development of methods to characterise single-cell or monolayer mechanical properties and link them to their functional behaviour. However, there remains a strong need to establish this link for three-dimensional (3D) multicellular aggregates, which better mimic tissue function. Here we present a platform to actuate and observe many such aggregates within one deformable micro-device. The platform consists of a single polydimethylsiloxane piece cast on a 3D-printed mould and bonded to a glass slide or coverslip. It consists of a chamber containing cell spheroids, which is adjacent to air cavities that are fluidically independent. Controlling the air pressure in these air cavities leads to a vertical displacement of the chamber's ceiling. The device can be used in static or dynamic modes over time scales of seconds to hours, with displacement amplitudes from a fewµm to several tens of microns. Further, we show how the compression protocols can be used to obtain measurements of stiffness heterogeneities within individual co-culture spheroids, by comparing image correlations of spheroids at different levels of compression with finite element simulations. The labelling of the cells and their cytoskeleton is combined with image correlation methods to relate the structure of the co-culture spheroid with its mechanical properties at different locations. The device is compatible with various microscopy techniques, including confocal microscopy, which can be used to observe the displacements and rearrangements of single cells and neighbourhoods within the aggregate. The complete experimental and imaging platform can now be used to provide multi-scale measurements that link single-cell behaviour with the global mechanical response of the aggregates.
Collapse
Affiliation(s)
- Shreyansh Jain
- Institut Pasteur, Université Paris Cité, Physical Microfluidics and Bioengineering, 25-28 Rue du Dr Roux, 75015 Paris, France
- Laboratoire d' Hydrodynamique (LadHyX), CNRS, Ecole Polytechnique, Institut Polytechnique de Paris, 91128 Palaiseau, France
| | - Hiba Belkadi
- Institut Pasteur, Université Paris Cité, Physical Microfluidics and Bioengineering, 25-28 Rue du Dr Roux, 75015 Paris, France
- Laboratoire d' Hydrodynamique (LadHyX), CNRS, Ecole Polytechnique, Institut Polytechnique de Paris, 91128 Palaiseau, France
| | - Arthur Michaut
- Institut Pasteur, Université Paris Cité, Dynamic Regulation of Morphogenesis, 25-28 Rue du Dr Roux, 75015 Paris, France
| | - Sébastien Sart
- Institut Pasteur, Université Paris Cité, Physical Microfluidics and Bioengineering, 25-28 Rue du Dr Roux, 75015 Paris, France
- Laboratoire d' Hydrodynamique (LadHyX), CNRS, Ecole Polytechnique, Institut Polytechnique de Paris, 91128 Palaiseau, France
| | - Jérôme Gros
- Institut Pasteur, Université Paris Cité, Dynamic Regulation of Morphogenesis, 25-28 Rue du Dr Roux, 75015 Paris, France
| | - Martin Genet
- Laboratoire de Mécanique des Solides, CNRS, École Polytechnique, Institut Polytechnique de Paris, 91128 Palaiseau, France
- Inria, Palaiseau, France
| | - Charles N Baroud
- Institut Pasteur, Université Paris Cité, Physical Microfluidics and Bioengineering, 25-28 Rue du Dr Roux, 75015 Paris, France
- Laboratoire d' Hydrodynamique (LadHyX), CNRS, Ecole Polytechnique, Institut Polytechnique de Paris, 91128 Palaiseau, France
| |
Collapse
|
3
|
Genet M, Diaz J, Chapelle D, Moireau P. Reduced left ventricular dynamics modeling based on a cylindrical assumption. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2023; 39:e3711. [PMID: 37203282 DOI: 10.1002/cnm.3711] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/16/2022] [Revised: 02/11/2023] [Accepted: 04/02/2023] [Indexed: 05/20/2023]
Abstract
Biomechanical modeling and simulation is expected to play a significant role in the development of the next generation tools in many fields of medicine. However, full-order finite element models of complex organs such as the heart can be computationally very expensive, thus limiting their practical usability. Therefore, reduced models are much valuable to be used, for example, for pre-calibration of full-order models, fast predictions, real-time applications, and so forth. In this work, focused on the left ventricle, we develop a reduced model by defining reduced geometry & kinematics while keeping general motion and behavior laws, allowing to derive a reduced model where all variables & parameters have a strong physical meaning. More specifically, we propose a reduced ventricular model based on cylindrical geometry & kinematics, which allows to describe the myofiber orientation through the ventricular wall and to represent contraction patterns such as ventricular twist, two important features of ventricular mechanics. Our model is based on the original cylindrical model of Guccione, McCulloch, & Waldman (1991); Guccione, Waldman, & McCulloch (1993), albeit with multiple differences: we propose a fully dynamical formulation, integrated into an open-loop lumped circulation model, and based on a material behavior that incorporates a fine description of contraction mechanisms; moreover, the issue of the cylinder closure has been completely reformulated; our numerical approach is novel aswell, with consistent spatial (finite element) and time discretizations. Finally, we analyze the sensitivity of the model response to various numerical and physical parameters, and study its physiological response.
Collapse
Affiliation(s)
- Martin Genet
- LMS, École Polytechnique/CNRS/Institut Polytechnique de Paris, Palaiseau, France
- Inria, MΞDISIM Team, Inria Saclay-Ile de France, Palaiseau, France
| | - Jérôme Diaz
- LMS, École Polytechnique/CNRS/Institut Polytechnique de Paris, Palaiseau, France
- Inria, MΞDISIM Team, Inria Saclay-Ile de France, Palaiseau, France
| | - Dominique Chapelle
- LMS, École Polytechnique/CNRS/Institut Polytechnique de Paris, Palaiseau, France
- Inria, MΞDISIM Team, Inria Saclay-Ile de France, Palaiseau, France
| | - Philippe Moireau
- LMS, École Polytechnique/CNRS/Institut Polytechnique de Paris, Palaiseau, France
- Inria, MΞDISIM Team, Inria Saclay-Ile de France, Palaiseau, France
| |
Collapse
|
4
|
Laville C, Fetita C, Gille T, Brillet PY, Nunes H, Bernaudin JF, Genet M. Comparison of optimization parametrizations for regional lung compliance estimation using personalized pulmonary poromechanical modeling. Biomech Model Mechanobiol 2023; 22:1541-1554. [PMID: 36913005 PMCID: PMC10009868 DOI: 10.1007/s10237-023-01691-9] [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: 09/30/2022] [Accepted: 01/09/2023] [Indexed: 03/14/2023]
Abstract
Interstitial lung diseases, such as idiopathic pulmonary fibrosis (IPF) or post-COVID-19 pulmonary fibrosis, are progressive and severe diseases characterized by an irreversible scarring of interstitial tissues that affects lung function. Despite many efforts, these diseases remain poorly understood and poorly treated. In this paper, we propose an automated method for the estimation of personalized regional lung compliances based on a poromechanical model of the lung. The model is personalized by integrating routine clinical imaging data - namely computed tomography images taken at two breathing levels in order to reproduce the breathing kinematic-notably through an inverse problem with fully personalized boundary conditions that is solved to estimate patient-specific regional lung compliances. A new parametrization of the inverse problem is introduced in this paper, based on the combined estimation of a personalized breathing pressure in addition to material parameters, improving the robustness and consistency of estimation results. The method is applied to three IPF patients and one post-COVID-19 patient. This personalized model could help better understand the role of mechanics in pulmonary remodeling due to fibrosis; moreover, patient-specific regional lung compliances could be used as an objective and quantitative biomarker for improved diagnosis and treatment follow up for various interstitial lung diseases.
Collapse
Affiliation(s)
- Colin Laville
- Laboratoire de Mécanique des Solides, École Polytechnique/CNRS/IPP, Palaiseau, France
- Inria, Palaiseau, France
| | | | - Thomas Gille
- Hypoxie et Poumon, Université Sorbonne Paris Nord/INSERM, Bobigny, France
- Hôpital Avicenne, APHP, Bobigny, France
| | - Pierre-Yves Brillet
- Hypoxie et Poumon, Université Sorbonne Paris Nord/INSERM, Bobigny, France
- Hôpital Avicenne, APHP, Bobigny, France
| | - Hilario Nunes
- Hypoxie et Poumon, Université Sorbonne Paris Nord/INSERM, Bobigny, France
- Hôpital Avicenne, APHP, Bobigny, France
| | | | - Martin Genet
- Laboratoire de Mécanique des Solides, École Polytechnique/CNRS/IPP, Palaiseau, France
- Inria, Palaiseau, France
| |
Collapse
|
5
|
Arratia López P, Mella H, Uribe S, Hurtado DE, Sahli Costabal F. WarpPINN: Cine-MR image registration with physics-informed neural networks. Med Image Anal 2023; 89:102925. [PMID: 37598608 DOI: 10.1016/j.media.2023.102925] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2022] [Revised: 07/18/2023] [Accepted: 08/01/2023] [Indexed: 08/22/2023]
Abstract
The diagnosis of heart failure usually includes a global functional assessment, such as ejection fraction measured by magnetic resonance imaging. However, these metrics have low discriminate power to distinguish different cardiomyopathies, which may not affect the global function of the heart. Quantifying local deformations in the form of cardiac strain can provide helpful information, but it remains a challenge. In this work, we introduce WarpPINN, a physics-informed neural network to perform image registration to obtain local metrics of heart deformation. We apply this method to cine magnetic resonance images to estimate the motion during the cardiac cycle. We inform our neural network of the near-incompressibility of cardiac tissue by penalizing the Jacobian of the deformation field. The loss function has two components: an intensity-based similarity term between the reference and the warped template images, and a regularizer that represents the hyperelastic behavior of the tissue. The architecture of the neural network allows us to easily compute the strain via automatic differentiation to assess cardiac activity. We use Fourier feature mappings to overcome the spectral bias of neural networks, allowing us to capture discontinuities in the strain field. The algorithm is tested on synthetic examples and on a cine SSFP MRI benchmark of 15 healthy volunteers, where it is trained to learn the deformation mapping of each case. We outperform current methodologies in landmark tracking and provide physiological strain estimations in the radial and circumferential directions. WarpPINN provides precise measurements of local cardiac deformations that can be used for a better diagnosis of heart failure and can be used for general image registration tasks. Source code is available at https://github.com/fsahli/WarpPINN.
Collapse
Affiliation(s)
| | - Hernán Mella
- School of Electrical Engineering, Pontificia Universidad Católica de Valparaíso, Valparaíso, Chile
| | - Sergio Uribe
- Millennium Institute for Intelligent Healthcare Engineering, iHEALTH, Chile; Biomedical Imaging Center, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Daniel E Hurtado
- Department of Structural and Geotechnical Engineering, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile; Institute for Biological and Medical Engineering, Schools of Engineering, Medicine and Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Francisco Sahli Costabal
- Millennium Institute for Intelligent Healthcare Engineering, iHEALTH, Chile; Institute for Biological and Medical Engineering, Schools of Engineering, Medicine and Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, Chile; Department of Mechanical and Metallurgical Engineering, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile.
| |
Collapse
|
6
|
Computational Analysis of Cardiac Contractile Function. Curr Cardiol Rep 2022; 24:1983-1994. [PMID: 36301405 PMCID: PMC10091868 DOI: 10.1007/s11886-022-01814-1] [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] [Accepted: 10/14/2022] [Indexed: 01/11/2023]
Abstract
PURPOSE OF REVIEW Heart failure results in the high incidence and mortality all over the world. Mechanical properties of myocardium are critical determinants of cardiac function, with regional variations in myocardial contractility demonstrated within infarcted ventricles. Quantitative assessment of cardiac contractile function is therefore critical to identify myocardial infarction for the early diagnosis and therapeutic intervention. RECENT FINDINGS Current advancement of cardiac functional assessments is in pace with the development of imaging techniques. The methods tailored to advanced imaging have been widely used in cardiac magnetic resonance, echocardiography, and optical microscopy. In this review, we introduce fundamental concepts and applications of representative methods for each imaging modality used in both fundamental research and clinical investigations. All these methods have been designed or developed to quantify time-dependent 2-dimensional (2D) or 3D cardiac mechanics, holding great potential to unravel global or regional myocardial deformation and contractile function from end-systole to end-diastole. Computational methods to assess cardiac contractile function provide a quantitative insight into the analysis of myocardial mechanics during cardiac development, injury, and remodeling.
Collapse
|
7
|
Pawar A, Li L, Gosain AK, Umulis DM, Tepole AB. PDE-constrained shape registration to characterize biological growth and morphogenesis from imaging data. ENGINEERING WITH COMPUTERS 2022; 38:3909-3924. [PMID: 38046797 PMCID: PMC10691863 DOI: 10.1007/s00366-022-01682-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/04/2022] [Accepted: 05/20/2022] [Indexed: 12/05/2023]
Abstract
We propose a PDE-constrained shape registration algorithm that captures the deformation and growth of biological tissue from imaging data. Shape registration is the process of evaluating optimum alignment between pairs of geometries through a spatial transformation function. We start from our previously reported work, which uses 3D tensor product B-spline basis functions to interpolate 3D space. Here, the movement of the B-spline control points, composed with an implicit function describing the shape of the tissue, yields the total deformation gradient field. The deformation gradient is then split into growth and elastic contributions. The growth tensor captures addition of mass, i.e. growth, and evolves according to a constitutive equation which is usually a function of the elastic deformation. Stress is generated in the material due to the elastic component of the deformation alone. The result of the registration is obtained by minimizing a total energy functional which includes: a distance measure reflecting similarity between the shapes, and the total elastic energy accounting for the growth of the tissue. We apply the proposed shape registration framework to study zebrafish embryo epiboly process and tissue expansion during skin reconstruction surgery. We anticipate that our PDE-constrained shape registration method will improve our understanding of biological and medical problems in which tissues undergo extreme deformations over time.
Collapse
Affiliation(s)
- Aishwarya Pawar
- School of Mechanical Engineering, Purdue University, 585 Purdue Mall, West Lafayette, 47907, Indiana, USA
| | - Linlin Li
- Weldon School of Biomedical Engineering, Purdue University, 206 S Martin Jischke Dr, West Lafayette, 47907, Indiana, USA
| | - Arun K. Gosain
- Lurie Children’s Hospital, Northwestern University, 225 East Chicago Ave, Chicago, 60611, Illinois, USA
| | - David M. Umulis
- Weldon School of Biomedical Engineering, Purdue University, 206 S Martin Jischke Dr, West Lafayette, 47907, Indiana, USA
| | - Adrian Buganza Tepole
- School of Mechanical Engineering, Purdue University, 585 Purdue Mall, West Lafayette, 47907, Indiana, USA
- Weldon School of Biomedical Engineering, Purdue University, 206 S Martin Jischke Dr, West Lafayette, 47907, Indiana, USA
| |
Collapse
|
8
|
Berberoğlu E, Stoeck CT, Kozerke S, Genet M. Quantification of left ventricular strain and torsion by joint analysis of 3D tagging and cine MR images. Med Image Anal 2022; 82:102598. [PMID: 36049451 DOI: 10.1016/j.media.2022.102598] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2021] [Revised: 06/30/2022] [Accepted: 08/11/2022] [Indexed: 11/19/2022]
Abstract
Cardiovascular magnetic resonance (CMR) imaging is the gold standard for the non-invasive assessment of left-ventricular (LV) function. Prognostic value of deformation metrics extracted directly from regular SSFP CMR images has been shown by numerous studies in the clinical setting, but with some limitations to detect torsion of the myocardium. Tagged CMR introduces trackable features in the myocardium that allow for the assessment of local myocardial deformation, including torsion; it is, however, limited in the quantification of radial strain, which is a decisive metric for assessing the contractility of the heart. In order to improve SSFP-only and tagged-only approaches, we propose to combine the advantages of both image types by fusing global shape motion obtained from SSFP images with the local deformation obtained from tagged images. To this end, tracking is first performed on SSFP images, and subsequently, the resulting motion is utilized to mask and track tagged data. Our implementation is based on a recent finite element-based motion tracking tool with mechanical regularization. Joint SSFP and tagged images registration performance is assessed based on deformation metrics including LV strain and twist using human and in-house porcine datasets. Results show that joint analysis of SSFP and 3DTAG images provides better quantification of LV strain and twist as either data source alone.
Collapse
Affiliation(s)
- Ezgi Berberoğlu
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland; Laboratoire de Mécanique des Solides (LMS), École Polytechnique/C.N.R.S./Institut Polytechnique de Paris, Palaiseau, France; MΞDISIM team, Inria, Palaiseau, France
| | - Christian T Stoeck
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Sebastian Kozerke
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Martin Genet
- Laboratoire de Mécanique des Solides (LMS), École Polytechnique/C.N.R.S./Institut Polytechnique de Paris, Palaiseau, France; MΞDISIM team, Inria, Palaiseau, France.
| |
Collapse
|
9
|
Patte C, Brillet PY, Fetita C, Bernaudin JF, Gille T, Nunes H, Chapelle D, Genet M. Estimation of Regional Pulmonary Compliance in Idiopathic Pulmonary Fibrosis Based On Personalized Lung Poromechanical Modeling. J Biomech Eng 2022; 144:1139545. [PMID: 35292805 DOI: 10.1115/1.4054106] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Indexed: 11/08/2022]
Abstract
Pulmonary function is tightly linked to the lung mechanical behavior, especially large deformation during breathing. Interstitial lung diseases, such as Idiopathic Pulmonary Fibrosis (IPF), have an impact on the pulmonary mechanics and consequently alter lung function. However, IPF remains poorly understood, poorly diagnosed and poorly treated. Currently, the mechanical impact of such diseases is assessed by pressure-volume curves, giving only global information. We developed a poromechanical model of the lung that can be personalized to a patient based on routine clinical data. The personalization pipeline uses clinical data, mainly CT-images at two time steps and involves the formulation of an inverse problem to estimate regional compliances. The estimation problem can be formulated both in terms of "effective", i.e., without considering the mixture porosity, or "rescaled", i.e., where the first-order effect of the porosity has been taken into account, compliances. Regional compliances are estimated for one control subject and three IPF patients, allowing to quantify the IPF-induced tissue stiffening. This personalized model could be used in the clinic as an objective and quantitative tool for IPF diagnosis.
Collapse
Affiliation(s)
- Cécile Patte
- Inria, Palaiseau, France, Laboratoire de Mécanique des Solides, École Polytechnique/CNRS/IPP, Palaiseau, France
| | - Pierre-Yves Brillet
- Hypoxie et Poumon, Universit é Sorbonne Paris Nord/INSERM, Bobigny, France; Hôpital Avicenne, APHP, Bobigny, France
| | - Catalin Fetita
- SAMOVAR, Telecom SudParis/Institut Mines-Télécom/IPP, Évry, France
| | | | - Thomas Gille
- Hypoxie et Poumon, Universit é Sorbonne Paris Nord/INSERM, Bobigny, France; Hôpital Avicenne, APHP, Bobigny, France
| | - Hilario Nunes
- Hypoxie et Poumon, Universit é Sorbonne Paris Nord/INSERM, Bobigny, France; Hôpital Avicenne, APHP, Bobigny, France
| | - Dominique Chapelle
- Inria, Palaiseau, France, Laboratoire de Mécanique des Solides, École Polytechnique/CNRS/IPP, Palaiseau, France
| | - Martin Genet
- Laboratoire de Mecanique des Solides, École Polytechnique/CNRS/IPP, Palaiseau, France; Inria, Palaiseau, France
| |
Collapse
|
10
|
Stimm J, Nordsletten DA, Jilberto J, Miller R, Berberoğlu E, Kozerke S, Stoeck CT. Personalization of biomechanical simulations of the left ventricle by in-vivo cardiac DTI data: Impact of fiber interpolation methods. Front Physiol 2022; 13:1042537. [PMID: 36518106 PMCID: PMC9742433 DOI: 10.3389/fphys.2022.1042537] [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: 09/12/2022] [Accepted: 11/14/2022] [Indexed: 11/29/2022] Open
Abstract
Simulations of cardiac electrophysiology and mechanics have been reported to be sensitive to the microstructural anisotropy of the myocardium. Consequently, a personalized representation of cardiac microstructure is a crucial component of accurate, personalized cardiac biomechanical models. In-vivo cardiac Diffusion Tensor Imaging (cDTI) is a non-invasive magnetic resonance imaging technique capable of probing the heart's microstructure. Being a rather novel technique, issues such as low resolution, signal-to noise ratio, and spatial coverage are currently limiting factors. We outline four interpolation techniques with varying degrees of data fidelity, different amounts of smoothing strength, and varying representation error to bridge the gap between the sparse in-vivo data and the model, requiring a 3D representation of microstructure across the myocardium. We provide a workflow to incorporate in-vivo myofiber orientation into a left ventricular model and demonstrate that personalized modelling based on fiber orientations from in-vivo cDTI data is feasible. The interpolation error is correlated with a trend in personalized parameters and simulated physiological parameters, strains, and ventricular twist. This trend in simulation results is consistent across material parameter settings and therefore corresponds to a bias introduced by the interpolation method. This study suggests that using a tensor interpolation approach to personalize microstructure with in-vivo cDTI data, reduces the fiber uncertainty and thereby the bias in the simulation results.
Collapse
Affiliation(s)
- Johanna Stimm
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - David A Nordsletten
- Department of Biomedical Engineering and Cardiac Surgery, University of Michigan, Ann Arbor, MI, United States.,School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Javiera Jilberto
- Department of Biomedical Engineering and Cardiac Surgery, University of Michigan, Ann Arbor, MI, United States
| | - Renee Miller
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Ezgi Berberoğlu
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Sebastian Kozerke
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Christian T Stoeck
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland.,Division of Surgical Research, University Hospital Zurich, University Zurich, Zurich, Switzerland
| |
Collapse
|
11
|
Berberoğlu E, Stoeck CT, Moireau P, Kozerke S, Genet M. In-silico study of accuracy and precision of left-ventricular strain quantification from 3D tagged MRI. PLoS One 2021; 16:e0258965. [PMID: 34739495 PMCID: PMC8570486 DOI: 10.1371/journal.pone.0258965] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2021] [Accepted: 10/08/2021] [Indexed: 11/18/2022] Open
Abstract
Cardiac Magnetic Resonance Imaging (MRI) allows quantifying myocardial tissue deformation and strain based on the tagging principle. In this work, we investigate accuracy and precision of strain quantification from synthetic 3D tagged MRI using equilibrated warping. To this end, synthetic biomechanical left-ventricular tagged MRI data with varying tag distance, spatial resolution and signal-to-noise ratio (SNR) were generated and processed to quantify errors in radial, circumferential and longitudinal strains relative to ground truth. Results reveal that radial strain is more sensitive to image resolution and noise than the other strain components. The study also shows robustness of quantifying circumferential and longitudinal strain in the presence of geometrical inconsistencies of 3D tagged data. In conclusion, our study points to the need for higher-resolution 3D tagged MRI than currently available in practice in order to achieve sufficient accuracy of radial strain quantification.
Collapse
Affiliation(s)
- Ezgi Berberoğlu
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Christian T. Stoeck
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Philippe Moireau
- MΞDISIM team, Inria, Palaiseau, France
- Laboratoire de Mécanique des Solides (LMS), École Polytechnique, C.N.R.S., Institut Polytechnique de Paris, Palaiseau, France
| | - Sebastian Kozerke
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Martin Genet
- MΞDISIM team, Inria, Palaiseau, France
- Laboratoire de Mécanique des Solides (LMS), École Polytechnique, C.N.R.S., Institut Polytechnique de Paris, Palaiseau, France
- * E-mail:
| |
Collapse
|
12
|
Gusseva M, Hussain T, Friesen CH, Moireau P, Tandon A, Patte C, Genet M, Hasbani K, Greil G, Chapelle D, Chabiniok R. Biomechanical Modeling to Inform Pulmonary Valve Replacement in Tetralogy of Fallot Patients After Complete Repair. Can J Cardiol 2021; 37:1798-1807. [PMID: 34216743 PMCID: PMC9810481 DOI: 10.1016/j.cjca.2021.06.018] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Revised: 06/05/2021] [Accepted: 06/26/2021] [Indexed: 01/07/2023] Open
Abstract
BACKGROUND A biomechanical model of the heart can be used to incorporate multiple data sources (electrocardiography, imaging, invasive hemodynamics). The purpose of this study was to use this approach in a cohort of patients with tetralogy of Fallot after complete repair (rTOF) to assess comparative influences of residual right ventricular outflow tract obstruction (RVOTO) and pulmonary regurgitation on ventricular health. METHODS Twenty patients with rTOF who underwent percutaneous pulmonary valve replacement (PVR) and cardiovascular magnetic resonance imaging were included in this retrospective study. Biomechanical models specific to individual patient and physiology (before and after PVR) were created and used to estimate the RV myocardial contractility. The ability of models to capture post-PVR changes of right ventricular (RV) end-diastolic volume (EDV) and effective flow in the pulmonary artery (Qeff) was also compared with expected values. RESULTS RV contractility before PVR (mean 66 ± 16 kPa, mean ± standard deviation) was increased in patients with rTOF compared with normal RV (38-48 kPa) (P < 0.05). The contractility decreased significantly in all patients after PVR (P < 0.05). Patients with predominantly RVOTO demonstrated greater reduction in contractility (median decrease 35%) after PVR than those with predominant pulmonary regurgitation (median decrease 11%). The model simulated post-PVR decreased EDV for the majority and suggested an increase of Qeff-both in line with published data. CONCLUSIONS This study used a biomechanical model to synthesize multiple clinical inputs and give an insight into RV health. Individualized modeling allows us to predict the RV response to PVR. Initial data suggest that residual RVOTO imposes greater ventricular work than isolated pulmonary regurgitation.
Collapse
Affiliation(s)
- Maria Gusseva
- Inria, Palaiseau, France,LMS, École Polytechnique, CNRS, Institut Polytechnique de Paris, Palaiseau, France
| | - Tarique Hussain
- Division of Pediatric Cardiology, Department of Pediatrics, UT Southwestern Medical Center, Dallas, Texas, USA
| | - Camille Hancock Friesen
- Division of Pediatric Cardiothoracic Surgery, Department of Pediatrics, UT Southwestern Medical Center, Dallas, Texas, USA
| | - Philippe Moireau
- Inria, Palaiseau, France,LMS, École Polytechnique, CNRS, Institut Polytechnique de Paris, Palaiseau, France
| | - Animesh Tandon
- Division of Pediatric Cardiology, Department of Pediatrics, UT Southwestern Medical Center, Dallas, Texas, USA
| | - Cécile Patte
- Inria, Palaiseau, France,LMS, École Polytechnique, CNRS, Institut Polytechnique de Paris, Palaiseau, France
| | - Martin Genet
- Inria, Palaiseau, France,LMS, École Polytechnique, CNRS, Institut Polytechnique de Paris, Palaiseau, France
| | - Keren Hasbani
- Division of Pediatric Cardiology, Department of Pediatrics, Dell Medical School, University of Texas, Austin, Texas, USA
| | - Gerald Greil
- Division of Pediatric Cardiology, Department of Pediatrics, UT Southwestern Medical Center, Dallas, Texas, USA
| | - Dominique Chapelle
- Inria, Palaiseau, France,LMS, École Polytechnique, CNRS, Institut Polytechnique de Paris, Palaiseau, France
| | - Radomír Chabiniok
- Inria, Palaiseau, France,LMS, École Polytechnique, CNRS, Institut Polytechnique de Paris, Palaiseau, France,Division of Pediatric Cardiology, Department of Pediatrics, UT Southwestern Medical Center, Dallas, Texas, USA,School of Biomedical Engineering & Imaging Sciences, St Thomas’ Hospital, King’s College London, London, United Kingdom,Department of Mathematics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, Prague, Czech Republic
| |
Collapse
|
13
|
Weissmann J, Charles CJ, Richards AM, Yap CH, Marom G. Cardiac mesh morphing method for finite element modeling of heart failure with preserved ejection fraction. J Mech Behav Biomed Mater 2021; 126:104937. [PMID: 34979481 DOI: 10.1016/j.jmbbm.2021.104937] [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: 09/19/2021] [Revised: 10/21/2021] [Accepted: 10/24/2021] [Indexed: 10/20/2022]
Abstract
Numerical modeling of heart biomechanics can realistically capture morphological variations in diseases and has been helpful in advancing our understanding of the physiology. Subject-specific models require anatomic representation of medical images, and it is desirable to have a consistently repeatable models for any given morphology. In this study, we propose a novel and easily adaptable cardiac reconstruction algorithm by morphing an existing discretized mesh of an advanced finite element (FE) model, to match anatomies acquired from porcine cardiac magnetic resonance imaging (cMRI) scans. The morphing algorithm involves iterative FE simulations with visco-hyperelastic material properties. The living heart porcine model (LHPM) was chosen as the input baseline FE mesh, in order to preserve detailed anatomical features that cannot be captured in routine scans such as myofiber orientations and conduction pathways. The algorithm was demonstrated for the recreation of porcine hearts of a healthy subject and of a subject induced with heart failure with preserved ejection fraction (HFpEF) conditions, where there were substantial hypertrophy and anatomical alterations. We further used the morphed meshes for FE modeling of cardiac contraction and relaxation, thus demonstrating the applicability of the proposed algorithm in producing viable meshes. The results show that our algorithm can recreate the characteristic anatomical changes of cardiac remodeling, including heart muscle thickening, as well as replicate the reduction in ventricular volume. This algorithm allows for the creation of subject-specific models with the same mesh connectivity, thus enabling spatial comparison and analysis of pathologic progress.
Collapse
Affiliation(s)
| | - Christopher J Charles
- Department of Surgery, Yong Loo Lin School of Medicine, National University of Singapore, Singapore; Cardiovascular Research Institute, National University of Singapore, Singapore; Christchurch Heart Institute, Department of Medicine, University of Otago, Christchurch, New Zealand
| | - A Mark Richards
- Cardiovascular Research Institute, National University of Singapore, Singapore; Christchurch Heart Institute, Department of Medicine, University of Otago, Christchurch, New Zealand
| | - Choon Hwai Yap
- Department of Bioengineering, Imperial College London, UK
| | - Gil Marom
- School of Mechanical Engineering, Tel Aviv University, Israel.
| |
Collapse
|
14
|
Cesarovic N, Weisskopf M, Kron M, Glaus L, Peper ES, Buoso S, Suendermann S, Canic M, Falk V, Kozerke S, Emmert MY, Stoeck CT. Septaly Oriented Mild Aortic Regurgitant Jets Negatively Influence Left Ventricular Blood Flow-Insights From 4D Flow MRI Animal Study. Front Cardiovasc Med 2021; 8:711099. [PMID: 34434980 PMCID: PMC8380779 DOI: 10.3389/fcvm.2021.711099] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2021] [Accepted: 07/06/2021] [Indexed: 11/23/2022] Open
Abstract
Objectives: Paravalvular leakage (PVL) and eccentric aortic regurgitation remain a major clinical concern in patients receiving transcatheter aortic valve replacement (TAVR), and regurgitant volume remains the main readout parameter in clinical assessment. In this work we investigate the effect of jet origin and trajectory of mild aortic regurgitation on left ventricular hemodynamics in a porcine model. Methods: A pig model of mild aortic regurgitation/PVL was established by transcatheter piercing and dilating the non-coronary (NCC) or right coronary cusp (RCC) of the aortic valve close to the valve annulus. The interaction between regurgitant blood and LV hemodynamics was assessed by 4D flow cardiovascular MRI. Results: Six RCC, six NCC, and two control animals were included in the study and with one dropout in the NCC group, the success rate of model creation was 93%. Regurgitant jets originating from NCC were directed along the ventricular side of the anterior mitral leaflet and integrated well into the diastolic vortex forming in the left ventricular outflow tract. However, jets from the RCC were orientated along the septum colliding with flow within the vortex, and progressing down to the apex. As a consequence, the presence as well as the area of the vortex was reduced at the site of impact compared to the NCC group. Impairment of vortex formation was localized to the area of impact and not the entire vortex ring. Blood from the NCC jet was largely ejected during the following systole, whereas ejection of large portion of RCC blood was protracted. Conclusions: Even for mild regurgitation, origin and trajectory of the regurgitant jet does cause a different effect on LV hemodynamics. Septaly oriented jets originating from RCC collide with the diastolic vortex, reduce its size, and reach the apical region of the left ventricle where blood resides extendedly. Hence, RCC jets display hemodynamic features which may have a potential negative impact on the long-term burden to the heart.
Collapse
Affiliation(s)
- Nikola Cesarovic
- Department of Health Sciences and Technology, Swiss Federal Institute of Technology, Zurich, Switzerland.,Department of Cardiothoracic and Vascular Surgery, German Heart Center Berlin, Berlin, Germany
| | - Miriam Weisskopf
- Division of Surgical Research, University Hospital Zurich, University of Zurich, Zurich, Switzerland
| | - Mareike Kron
- Division of Surgical Research, University Hospital Zurich, University of Zurich, Zurich, Switzerland
| | - Lukas Glaus
- Department of Health Sciences and Technology, Swiss Federal Institute of Technology, Zurich, Switzerland
| | - Eva S Peper
- Institute for Biomedical Engineering, University and ETH Zürich, Zurich, Switzerland
| | - Stefano Buoso
- Institute for Biomedical Engineering, University and ETH Zürich, Zurich, Switzerland
| | - Simon Suendermann
- Department of Cardiothoracic and Vascular Surgery, German Heart Center Berlin, Berlin, Germany.,Department of Cardiovascular Surgery, Charité-Universitätsmedizin Berlin, Berlin, Germany
| | - Marko Canic
- Division of Surgical Research, University Hospital Zurich, University of Zurich, Zurich, Switzerland
| | - Volkmar Falk
- Department of Health Sciences and Technology, Swiss Federal Institute of Technology, Zurich, Switzerland.,Department of Cardiothoracic and Vascular Surgery, German Heart Center Berlin, Berlin, Germany.,Department of Cardiovascular Surgery, Charité-Universitätsmedizin Berlin, Berlin, Germany
| | - Sebastian Kozerke
- Institute for Biomedical Engineering, University and ETH Zürich, Zurich, Switzerland
| | - Maximilian Y Emmert
- Department of Cardiothoracic and Vascular Surgery, German Heart Center Berlin, Berlin, Germany.,Department of Cardiovascular Surgery, Charité-Universitätsmedizin Berlin, Berlin, Germany.,Institute for Regenerative Medicine, University of Zurich, Zurich, Switzerland
| | - Christian T Stoeck
- Institute for Biomedical Engineering, University and ETH Zürich, Zurich, Switzerland
| |
Collapse
|
15
|
Left Ventricular Torsion Obtained Using Equilibrated Warping in Patients with Repaired Tetralogy of Fallot. Pediatr Cardiol 2021; 42:1275-1283. [PMID: 33900430 PMCID: PMC9753563 DOI: 10.1007/s00246-021-02608-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/06/2020] [Accepted: 04/07/2021] [Indexed: 10/21/2022]
Abstract
Patients after surgical repair of Tetralogy of Fallot (rTOF) may suffer a decrease in left ventricular (LV) function. The aim of our study is to evaluate a novel method of assessing LV torsion in patients with rTOF, as an early indicator of systolic LV dysfunction. Motion tracking based on image registration regularized by the equilibrium gap principle, known as equilibrated warping, was employed to assess LV torsion. Seventy-six cases of rTOF and ten normal controls were included. The group of controls was assessed for reproducibility using both equilibrated warping and standard clinical tissue tracking software (CVI42, version 5.10.1, Calgary, Canada). Patients were dichotomized into two groups: normal vs. loss of torsion. Torsion by equilibrated warping was successfully obtained in 68 of 76 (89%) patients and 9 of 10 (90%) controls. For equilibrated warping, the intra- and interobserver coefficients of variation were 0.095 and 0.117, respectively, compared to 0.260 and 0.831 for tissue tracking by standard clinical software. The intra- and inter-observer intraclass correlation coefficients for equilibrated warping were 0.862 and 0.831, respectively, compared to 0.992 and 0.648 for tissue tracking. Loss of torsion was noted in 32 of the 68 (47%) patients with rTOF. There was no difference in LV or RV volumes or ejection fraction between these groups. The assessment of LV torsion by equilibrated warping is feasible and shows good reliability. Loss of torsion is common in patients with rTOF and its robust assessment might contribute into uncovering heart failure in an earlier stage.
Collapse
|
16
|
Mella H, Mura J, Wang H, Taylor MD, Chabiniok R, Tintera J, Sotelo J, Uribe S. HARP-I: A Harmonic Phase Interpolation Method for the Estimation of Motion From Tagged MR Images. IEEE TRANSACTIONS ON MEDICAL IMAGING 2021; 40:1240-1252. [PMID: 33434127 DOI: 10.1109/tmi.2021.3051092] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We proposed a novel method called HARP-I, which enhances the estimation of motion from tagged Magnetic Resonance Imaging (MRI). The harmonic phase of the images is unwrapped and treated as noisy measurements of reference coordinates on a deformed domain, obtaining motion with high accuracy using Radial Basis Functions interpolations. Results were compared against Shortest Path HARP Refinement (SP-HR) and Sine-wave Modeling (SinMod), two harmonic image-based techniques for motion estimation from tagged images. HARP-I showed a favorable similarity with both methods under noise-free conditions, whereas a more robust performance was found in the presence of noise. Cardiac strain was better estimated using HARP-I at almost any motion level, giving strain maps with less artifacts. Additionally, HARP-I showed better temporal consistency as a new method was developed to fix phase jumps between frames. In conclusion, HARP-I showed to be a robust method for the estimation of motion and strain under ideal and non-ideal conditions.
Collapse
|
17
|
Perotti LE, Verzhbinsky IA, Moulin K, Cork TE, Loecher M, Balzani D, Ennis DB. Estimating cardiomyofiber strain in vivo by solving a computational model. Med Image Anal 2021; 68:101932. [PMID: 33383331 PMCID: PMC7956226 DOI: 10.1016/j.media.2020.101932] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2020] [Revised: 11/22/2020] [Accepted: 11/27/2020] [Indexed: 11/19/2022]
Abstract
Since heart contraction results from the electrically activated contraction of millions of cardiomyocytes, a measure of cardiomyocyte shortening mechanistically underlies cardiac contraction. In this work we aim to measure preferential aggregate cardiomyocyte ("myofiber") strains based on Magnetic Resonance Imaging (MRI) data acquired to measure both voxel-wise displacements through systole and myofiber orientation. In order to reduce the effect of experimental noise on the computed myofiber strains, we recast the strains calculation as the solution of a boundary value problem (BVP). This approach does not require a calibrated material model, and consequently is independent of specific myocardial material properties. The solution to this auxiliary BVP is the displacement field corresponding to assigned values of myofiber strains. The actual myofiber strains are then determined by minimizing the difference between computed and measured displacements. The approach is validated using an analytical phantom, for which the ground-truth solution is known. The method is applied to compute myofiber strains using in vivo displacement and myofiber MRI data acquired in a mid-ventricular left ventricle section in N=8 swine subjects. The proposed method shows a more physiological distribution of myofiber strains compared to standard approaches that directly differentiate the displacement field.
Collapse
Affiliation(s)
- Luigi E Perotti
- Department of Mechanical and Aerospace Engineering, University of Central Florida, Orlando, FL, USA.
| | - Ilya A Verzhbinsky
- Department of Radiology, Stanford University, Stanford, CA, USA; Medical Scientist Training Program, University of California, San Diego, La Jolla, USA
| | - Kévin Moulin
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Tyler E Cork
- Department of Radiology, Stanford University, Stanford, CA, USA; Department of Bioengineering, Stanford University, Stanford, CA, USA
| | - Michael Loecher
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Daniel Balzani
- Chair of Continuum Mechanics, Ruhr University Bochum, Bochum, Germany
| | - Daniel B Ennis
- Department of Radiology, Stanford University, Stanford, CA, USA
| |
Collapse
|
18
|
Genet M, Patte C, Fetita C, Brillet PY, Chapelle D. Personalized pulmonary poromechanics. Comput Methods Biomech Biomed Engin 2020. [DOI: 10.1080/10255842.2020.1812842] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- M. Genet
- Laboratoire de Mécanique des Solides (LMS), École Polytechnique, CNRS, Institut Polytechnique de Paris, Palaiseau, France
- Inria, Palaiseau, France
| | - C. Patte
- Laboratoire de Mécanique des Solides (LMS), École Polytechnique, CNRS, Institut Polytechnique de Paris, Palaiseau, France
- Inria, Palaiseau, France
| | - C. Fetita
- SAMOVAR, Telecom SudParis, Institut Mines-Telecom, Paris, France
- Institut Polytechnique de Paris, Évry-Courcouronnes, France
| | - P.-Y. Brillet
- Hypoxie et Poumon, Université Sorbonne Paris Nord, Villetaneuse, France
- INSERM, Bobigny, France
| | - D. Chapelle
- Laboratoire de Mécanique des Solides (LMS), École Polytechnique, CNRS, Institut Polytechnique de Paris, Palaiseau, France
- Inria, Palaiseau, France
| |
Collapse
|
19
|
Wiputra H, Chan WX, Foo YY, Ho S, Yap CH. Cardiac motion estimation from medical images: a regularisation framework applied on pairwise image registration displacement fields. Sci Rep 2020; 10:18510. [PMID: 33116206 PMCID: PMC7595231 DOI: 10.1038/s41598-020-75525-4] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2020] [Accepted: 10/06/2020] [Indexed: 11/09/2022] Open
Abstract
Accurate cardiac motion estimation from medical images such as ultrasound is important for clinical evaluation. We present a novel regularisation layer for cardiac motion estimation that will be applied after image registration and demonstrate its effectiveness. The regularisation utilises a spatio-temporal model of motion, b-splines of Fourier, to fit to displacement fields from pairwise image registration. In the process, it enforces spatial and temporal smoothness and consistency, cyclic nature of cardiac motion, and better adherence to the stroke volume of the heart. Flexibility is further given for inclusion of any set of registration displacement fields. The approach gave high accuracy. When applied to human adult Ultrasound data from a Cardiac Motion Analysis Challenge (CMAC), the proposed method is found to have 10% lower tracking error over CMAC participants. Satisfactory cardiac motion estimation is also demonstrated on other data sets, including human fetal echocardiography, chick embryonic heart ultrasound images, and zebrafish embryonic microscope images, with the average Dice coefficient between estimation motion and manual segmentation at 0.82-0.87. The approach of performing regularisation as an add-on layer after the completion of image registration is thus a viable option for cardiac motion estimation that can still have good accuracy. Since motion estimation algorithms are complex, dividing up regularisation and registration can simplify the process and provide flexibility. Further, owing to a large variety of existing registration algorithms, such an approach that is usable on any algorithm may be useful.
Collapse
Affiliation(s)
- Hadi Wiputra
- Department of Biomedical Engineering, National University of Singapore, Singapore, 117583, Singapore
| | - Wei Xuan Chan
- Department of Biomedical Engineering, National University of Singapore, Singapore, 117583, Singapore
| | - Yoke Yin Foo
- Department of Biomedical Engineering, National University of Singapore, Singapore, 117583, Singapore
| | - Sheldon Ho
- Department of Biomedical Engineering, National University of Singapore, Singapore, 117583, Singapore
| | - Choon Hwai Yap
- Department of Bioengineering, Imperial College London, London, SW7 2AZ, UK.
| |
Collapse
|
20
|
Zou H, Leng S, Xi C, Zhao X, Koh AS, Gao F, Tan JL, Tan RS, Allen JC, Lee LC, Genet M, Zhong L. Three-dimensional biventricular strains in pulmonary arterial hypertension patients using hyperelastic warping. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2020; 189:105345. [PMID: 31982668 PMCID: PMC7198336 DOI: 10.1016/j.cmpb.2020.105345] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/16/2019] [Revised: 01/16/2020] [Accepted: 01/16/2020] [Indexed: 06/10/2023]
Abstract
BACKGROUND AND OBJECTIVE Evaluation of biventricular function is an essential component of clinical management in pulmonary arterial hypertension (PAH). This study aims to examine the utility of biventricular strains derived from a model-to-image registration technique in PAH patients in comparison to age- and gender-matched normal controls. METHODS A three-dimensional (3D) model was reconstructed from cine short- and long-axis cardiac magnetic resonance (CMR) images and subsequently partitioned into right ventricle (RV), left ventricle (LV) and septum. The hyperelastic warping method was used to register the meshed biventricular finite element model throughout the cardiac cycle and obtain the corresponding biventricular circumferential, longitudinal and radial strains. RESULTS Intra- and inter-observer reproducibility of biventricular strains was excellent with all intra-class correlation coefficients > 0.84. 3D biventricular longitudinal, circumferential and radial strains for RV, LV and septum were significantly decreased in PAH patients compared with controls. Receiver operating characteristic (ROC) analysis showed that the 3D biventricular strains were better early markers (Area under the ROC curve = 0.96 for RV longitudinal strain) of ventricular dysfunction than conventional parameters such as two-dimensional strains and ejection fraction. CONCLUSIONS Our highly reproducible methodology holds potential for extending CMR imaging to characterize 3D biventricular strains, eventually leading to deeper understanding of biventricular mechanics in PAH.
Collapse
Affiliation(s)
- Hua Zou
- National Heart Research Institute Singapore, National Heart Centre Singapore, Singapore
| | - Shuang Leng
- National Heart Research Institute Singapore, National Heart Centre Singapore, Singapore
| | - Ce Xi
- Department of Mechanical Engineering, Michigan State University, MI, United States
| | - Xiaodan Zhao
- National Heart Research Institute Singapore, National Heart Centre Singapore, Singapore
| | - Angela S Koh
- National Heart Research Institute Singapore, National Heart Centre Singapore, Singapore; Duke-NUS Medical School, Singapore
| | - Fei Gao
- National Heart Research Institute Singapore, National Heart Centre Singapore, Singapore
| | - Ju Le Tan
- National Heart Research Institute Singapore, National Heart Centre Singapore, Singapore; Duke-NUS Medical School, Singapore
| | - Ru-San Tan
- National Heart Research Institute Singapore, National Heart Centre Singapore, Singapore; Duke-NUS Medical School, Singapore
| | | | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, MI, United States
| | - Martin Genet
- Mechanics Department & Solid Mechanics Laboratory, École Polytechnique (Paris-Saclay University), Palaiseau, France; M3DISIM research team, INRIA (Paris-Saclay University), Palaiseau, France
| | - Liang Zhong
- National Heart Research Institute Singapore, National Heart Centre Singapore, Singapore; Duke-NUS Medical School, Singapore.
| |
Collapse
|
21
|
Liu W, Wang Z. Current Understanding of the Biomechanics of Ventricular Tissues in Heart Failure. Bioengineering (Basel) 2019; 7:bioengineering7010002. [PMID: 31861916 PMCID: PMC7175293 DOI: 10.3390/bioengineering7010002] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2019] [Revised: 12/17/2019] [Accepted: 12/18/2019] [Indexed: 12/17/2022] Open
Abstract
Heart failure is the leading cause of death worldwide, and the most common cause of heart failure is ventricular dysfunction. It is well known that the ventricles are anisotropic and viscoelastic tissues and their mechanical properties change in diseased states. The tissue mechanical behavior is an important determinant of the function of ventricles. The aim of this paper is to review the current understanding of the biomechanics of ventricular tissues as well as the clinical significance. We present the common methods of the mechanical measurement of ventricles, the known ventricular mechanical properties including the viscoelasticity of the tissue, the existing computational models, and the clinical relevance of the ventricular mechanical properties. Lastly, we suggest some future research directions to elucidate the roles of the ventricular biomechanics in the ventricular dysfunction to inspire new therapies for heart failure patients.
Collapse
Affiliation(s)
- Wenqiang Liu
- School of Biomedical Engineering, Colorado State University, Fort Collins, CO 80523, USA;
| | - Zhijie Wang
- School of Biomedical Engineering, Colorado State University, Fort Collins, CO 80523, USA;
- Department of Mechanical Engineering, Colorado State University, Fort Collins, CO 80523, USA
- Correspondence:
| |
Collapse
|
22
|
Finsberg H, Xi C, Zhao X, Tan JL, Genet M, Sundnes J, Lee LC, Zhong L, Wall ST. Computational quantification of patient-specific changes in ventricular dynamics associated with pulmonary hypertension. Am J Physiol Heart Circ Physiol 2019; 317:H1363-H1375. [PMID: 31674809 DOI: 10.1152/ajpheart.00094.2019] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Pulmonary arterial hypertension (PAH) causes an increase in the mechanical loading imposed on the right ventricle (RV) that results in progressive changes to its mechanics and function. Here, we quantify the mechanical changes associated with PAH by assimilating clinical data consisting of reconstructed three-dimensional geometry, pressure, and volume waveforms, as well as regional strains measured in patients with PAH (n = 12) and controls (n = 6) within a computational modeling framework of the ventricles. Modeling parameters reflecting regional passive stiffness and load-independent contractility as indexed by the tissue active tension were optimized so that simulation results matched the measurements. The optimized parameters were compared with clinical metrics to find usable indicators associated with the underlying mechanical changes. Peak contractility of the RV free wall (RVFW) γRVFW,max was found to be strongly correlated and had an inverse relationship with the RV and left ventricle (LV) end-diastolic volume ratio (i.e., RVEDV/LVEDV) (RVEDV/LVEDV)+ 0.44, R2 = 0.77). Correlation with RV ejection fraction (R2 = 0.50) and end-diastolic volume index (R2 = 0.40) were comparatively weaker. Patients with with RVEDV/LVEDV > 1.5 had 25% lower γRVFW,max (P < 0.05) than that of the control. On average, RVFW passive stiffness progressively increased with the degree of remodeling as indexed by RVEDV/LVEDV. These results suggest a mechanical basis of using RVEDV/LVEDV as a clinical index for delineating disease severity and estimating RVFW contractility in patients with PAH.NEW & NOTEWORTHY This article presents patient-specific data assimilation of a patient cohort and physical description of clinical observations.
Collapse
Affiliation(s)
- Henrik Finsberg
- Simula Research Laboratory, Oslo, Norway.,Center for Cardiological Innovation, Oslo, Norway.,Department of Informatics, University of Oslo, Oslo, Norway
| | - Ce Xi
- Department of Mechanical Engineering, Michigan State University, East Lansing, Michigan
| | | | - Ju Le Tan
- National Heart Center Singapore, Singapore
| | - Martin Genet
- Mechanics Department and Solid Mechanics Laboratory, École Polytechnique/Le Centre national de la recherche scientifique/Paris-Saclay University, Palaiseau, France.,M3DISIM research team, Institut national de recherche en informatique et en automatique/Paris-Saclay University, Palaiseau, France
| | - Joakim Sundnes
- Simula Research Laboratory, Oslo, Norway.,Center for Cardiological Innovation, Oslo, Norway.,Department of Informatics, University of Oslo, Oslo, Norway
| | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, Michigan
| | - Liang Zhong
- National Heart Center Singapore, Singapore.,Duke-National University of Singapore Medical School, Singapore
| | - Samuel T Wall
- Simula Research Laboratory, Oslo, Norway.,Center for Cardiological Innovation, Oslo, Norway
| |
Collapse
|
23
|
Zou H, Xi C, Zhao X, Koh AS, Gao F, Su Y, Tan RS, Allen J, Lee LC, Genet M, Zhong L. Quantification of Biventricular Strains in Heart Failure With Preserved Ejection Fraction Patient Using Hyperelastic Warping Method. Front Physiol 2018; 9:1295. [PMID: 30283352 PMCID: PMC6156386 DOI: 10.3389/fphys.2018.01295] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2017] [Accepted: 08/28/2018] [Indexed: 12/19/2022] Open
Abstract
Heart failure (HF) imposes a major global health care burden on society and suffering on the individual. About 50% of HF patients have preserved ejection fraction (HFpEF). More intricate and comprehensive measurement-focused imaging of multiple strain components may aid in the diagnosis and elucidation of this disease. Here, we describe the development of a semi-automated hyperelastic warping method for rapid comprehensive assessment of biventricular circumferential, longitudinal, and radial strains that is physiological meaningful and reproducible. We recruited and performed cardiac magnetic resonance (CMR) imaging on 30 subjects [10 HFpEF, 10 HF with reduced ejection fraction patients (HFrEF) and 10 healthy controls]. In each subject, a three-dimensional heart model including left ventricle (LV), right ventricle (RV), and septum was reconstructed from CMR images. The hyperelastic warping method was used to reference the segmented model with the target images and biventricular circumferential, longitudinal, and radial strain-time curves were obtained. The peak systolic strains are then measured and analyzed in this study. Intra- and inter-observer reproducibility of the biventricular peak systolic strains was excellent with all ICCs > 0.92. LV peak systolic circumferential, longitudinal, and radial strain, respectively, exhibited a progressive decrease in magnitude from healthy control→HFpEF→HFrEF: control (-15.5 ± 1.90, -15.6 ± 2.06, 41.4 ± 12.2%); HFpEF (-9.37 ± 3.23, -11.3 ± 1.76, 22.8 ± 13.1%); HFrEF (-4.75 ± 2.74, -7.55 ± 1.75, 10.8 ± 4.61%). A similar progressive decrease in magnitude was observed for RV peak systolic circumferential, longitudinal and radial strain: control (-9.91 ± 2.25, -14.5 ± 2.63, 26.8 ± 7.16%); HFpEF (-7.38 ± 3.17, -12.0 ± 2.45, 21.5 ± 10.0%); HFrEF (-5.92 ± 3.13, -8.63 ± 2.79, 15.2 ± 6.33%). Furthermore, septum peak systolic circumferential, longitudinal, and radial strain magnitude decreased gradually from healthy control to HFrEF: control (-7.11 ± 1.81, 16.3 ± 3.23, 18.5 ± 8.64%); HFpEF (-6.11 ± 3.98, -13.4 ± 3.02, 12.5 ± 6.38%); HFrEF (-1.42 ± 1.36, -8.99 ± 2.96, 3.35 ± 2.95%). The ROC analysis indicated LV peak systolic circumferential strain to be the most sensitive marker for differentiating HFpEF from healthy controls. Our results suggest that the hyperelastic warping method with the CMR-derived strains may reveal subtle impairment in HF biventricular mechanics, in particular despite a "normal" ventricular ejection fraction in HFpEF.
Collapse
Affiliation(s)
- Hua Zou
- National Heart Centre Singapore, Singapore, Singapore
| | - Ce Xi
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - Xiaodan Zhao
- National Heart Centre Singapore, Singapore, Singapore
| | - Angela S Koh
- National Heart Centre Singapore, Singapore, Singapore.,Duke-NUS Medical School, National University of Singapore, Singapore, Singapore
| | - Fei Gao
- National Heart Centre Singapore, Singapore, Singapore
| | - Yi Su
- Institute of High Performance Computing, A∗STAR, Singapore, Singapore
| | - Ru-San Tan
- National Heart Centre Singapore, Singapore, Singapore.,Duke-NUS Medical School, National University of Singapore, Singapore, Singapore
| | - John Allen
- Duke-NUS Medical School, National University of Singapore, Singapore, Singapore
| | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - Martin Genet
- Mechanics Department and Solid Mechanics Laboratory, École Polytechnique, C.N.R.S., Université Paris-Saclay, Palaiseau, France.,M3DISIM Team, I.N.R.I.A, Université Paris-Saclay, Palaiseau, France
| | - Liang Zhong
- National Heart Centre Singapore, Singapore, Singapore.,Duke-NUS Medical School, National University of Singapore, Singapore, Singapore
| |
Collapse
|