Biomechanics of Human Fetal Hearts with Critical Aortic Stenosis

Myocardial biomechanical effects of fetal aortic valvuloplasty

Fetal critical aortic stenosis with evolving hypoplastic left heart syndrome (CAS-eHLHS) can progress to a univentricular (UV) birth malformation. Catheter-based fetal aortic valvuloplasty (FAV) can resolve stenosis and reduce the likelihood of malformation progression. However, we have limited understanding of the biomechanical impact of FAV and subsequent LV responses. Therefore, we performed image-based finite element (FE) modeling of 4 CAS-eHLHS fetal hearts, by performing iterative simulations to match image-based characteristics and then back-computing physiological parameters. We used pre-FAV simulations to conduct virtual FAV (vFAV) and compared pre-FAV and post-FAV simulations. vFAV simulations generally enabled partial restoration of several physiological features toward healthy levels, including increased stroke volume and myocardial strains, reduced aortic valve (AV) and mitral valve regurgitation (MVr) velocities, reduced LV and LA pressures, and reduced peak myofiber stress. FAV often leads to aortic valve regurgitation (AVr). Our simulations showed that AVr could compromise LV and LA depressurization but it could also significantly increase stroke volume and myocardial deformational stimuli. Post-FAV scans and simulations showed FAV enabled only partial reduction of the AV dissipative coefficient. Furthermore, LV contractility and peripheral vascular resistance could change in response to FAV, preventing decreases in AV velocity and LV pressure, compared with what would be anticipated from stenosis relief. This suggested that case-specific post-FAV modeling is required to fully capture cardiac functionality. Overall, image-based FE modeling could provide mechanistic details of the effects of FAV, but computational prediction of acute outcomes was difficult due to a patient-dependent physiological response to FAV.

Biventricular finite element modeling of the fetal heart in health and during critical aortic stenosis

Finite Element simulations are a robust way of investigating cardiac biomechanics. To date, it has only been performed with the left ventricle (LV) alone for fetal hearts, even though results are likely different with biventricular (BiV) simulations. In this research, we conduct BiV simulations of the fetal heart based on 4D echocardiography images to show that it can capture the biomechanics of the normal healthy fetal heart, as well as those of fetal aortic stenosis better than the LV alone simulations. We found that performing LV alone simulations resulted in overestimation of LV stresses and pressures, compared to BiV simulations. Interestingly, inserting a compliance between the LV and right ventricle (RV) in the lumped parameter model of the LV only simulation effectively resolved these overestimations, demonstrating that the septum could be considered to play a LV-RV pressure communication role. However, stresses and strains spatial patterns remained altered from BiV simulations after the addition of the compliance. The BiV simulations corroborated previous studies in showing disease effects on the LV, where fetal aortic stenosis (AS) drastically elevated LV pressures and reduced strains and stroke volumes, which were moderated down with the addition of mitral regurgitation (MR). However, BiV simulations enabled an evaluation of the RV as well, where we observed that effects of the AS and MR on pressures and stroke volumes were generally much smaller and less consistent. The BiV simulations also enabled investigations of septal dynamics, which showed a rightward shift with AS, and partial restoration with MR. Interestingly, AS tended to enhance RV stroke volume, but MR moderated that down.

Computational approaches for mechanobiology in cardiovascular development and diseases

The cardiovascular development in vertebrates evolves in response to genetic and mechanical cues. The dynamic interplay among mechanics, cell biology, and anatomy continually shapes the hydraulic networks, characterized by complex, non-linear changes in anatomical structure and blood flow dynamics. To better understand this interplay, a diverse set of molecular and computational tools has been used to comprehensively study cardiovascular mechanobiology. With the continual advancement of computational capacity and numerical techniques, cardiovascular simulation is increasingly vital in both basic science research for understanding developmental mechanisms and disease etiologies, as well as in clinical studies aimed at enhancing treatment outcomes. This review provides an overview of computational cardiovascular modeling. Beginning with the fundamental concepts of computational cardiovascular modeling, it navigates through the applications of computational modeling in investigating mechanobiology during cardiac development. Second, the article illustrates the utility of computational hemodynamic modeling in the context of treatment planning for congenital heart diseases. It then delves into the predictive potential of computational models for elucidating tissue growth and remodeling processes. In closing, we outline prevailing challenges and future prospects, underscoring the transformative impact of computational cardiovascular modeling in reshaping cardiovascular science and clinical practice.

Role of tissue biomechanics in the formation and function of myocardial trabeculae in zebrafish embryos

Cardiac trabeculae are uneven ventricular muscular structures that develop during early embryonic heart development at the outer curvature of the ventricle. Their biomechanical function is not completely understood, and while their formation is known to be mechanosensitive, it is unclear whether ventricular tissue internal stresses play an important role in their formation. Here, we performed imaging and image-based cardiac biomechanics simulations on zebrafish embryonic ventricles to investigate these issues. Microscopy-based ventricular strain measurements show that the appearance of trabeculae coincided with enhanced deformability of the ventricular wall. Image-based biomechanical simulations reveal that the presence of trabeculae reduces ventricular tissue internal stresses, likely acting as structural support in response to the geometry of the ventricle. Passive ventricular pressure-loading experiments further reveal that the formation of trabeculae is associated with a spatial homogenization of ventricular tissue stiffnesses in healthy hearts, but gata1 morphants with a disrupted trabeculation process retain a spatial stiffness heterogeneity. Our findings thus suggest that modulating ventricular wall deformability, stresses, and stiffness are among the biomechanical functions of trabeculae. Further, experiments with gata1 morphants reveal that a reduction in fluid pressures and consequently ventricular tissue internal stresses can disrupt trabeculation, but a subsequent restoration of ventricular tissue internal stresses via vasopressin rescues trabeculation, demonstrating that tissue stresses are important to trabeculae formation. Overall, we find that tissue biomechanics is important to the formation and function of embryonic heart trabeculation. KEY POINTS: Trabeculations are fascinating and important cardiac structures and their abnormalities are linked to embryonic demise. However, their function in the heart and their mechanobiological formation processes are not completely understood. Our imaging and modelling show that tissue biomechanics is the key here. We find that trabeculations enhance cardiac wall deformability, reduce fluid pressure stresses, homogenize wall stiffness, and have alignments that are optimal for providing load-bearing structural support for the heart. We further discover that high ventricular tissue internal stresses consequent to high fluid pressures are needed for trabeculation formation through a rescue experiment, demonstrating that myocardial tissue stresses are as important as fluid flow wall shear stresses for trabeculation formation.

Pre-intervention myocardial stress is a good predictor of aortic valvoluplasty outcome for fetal critical aortic stenosis and evolving HLHS

Fetal critical aortic stenosis with evolving hypoplastic left heart syndrome (CAS-eHLHS) causes biomechanical and functional aberrations, leading to a high risk of progression to hypoplastic left heart syndrome (HLHS) at birth. Fetal aortic valvuloplasty (FAV) can resolve outflow obstruction and may reduce progression risk. However, it is currently difficult to accurately predict which patients will respond to the intervention and become functionally biventricular (BV) at birth, as opposed to becoming functionally univentricular (UV). This prediction is important for patient selection, parental counselling, and surgical planning. Therefore, we investigated whether biomechanics parameters from pre-FAV image-based computations could robustly distinguish between CAS-eHLHS cases with BV or UV outcomes in a retrospective cohort. To do so we performed image-based finite element biomechanics modelling of nine CAS-eHLHS cases undergoing intervention and six healthy fetal control hearts, and found that a biomechanical parameter, peak systolic myofibre stress, showed a uniquely large difference between BV and UV cases, which had a larger magnitude effect than echocardiography parameters. A simplified equation was derived for quick and easy estimation of myofibre stress from echo measurements via principal component analysis. When tested on a retrospective cohort of 37 CAS-eHLHS cases, the parameter outperformed other parameters in predicting UV versus BV outcomes, and thus has a high potential of improving outcome predictions, if incorporated into patient selection procedures. Physiologically, high myocardial stresses likely indicate a healthier myocardium that can withstand high stresses and resist pathological remodelling, which can explain why it is a good predictor of BV outcomes. KEY POINTS: Predicting the morphological birth outcomes (univentricular versus biventricular) of fetal aortic valvuloplasty for fetal aortic stenosis with evolving HLHS is important for accurate patient selection, parental counselling and management decisions. Computational simulations show that a biomechanics parameter, pre-intervention peak systolic myofibre stress, is uniquely robust in distinguishing between such outcomes, outperforming all echo parameters. An empirical equation was developed to quickly compute peak systolic myofibre stress from routine echo measurements and was the best predictor of outcomes among a wide range of parameters tested.

Effects of myocardial sheetlet sliding on left ventricular function

Left ventricle myocardium has a complex micro-architecture, which was revealed to consist of myocyte bundles arranged in a series of laminar sheetlets. Recent imaging studies demonstrated that these sheetlets re-orientated and likely slided over each other during the deformations between systole and diastole, and that sheetlet dynamics were altered during cardiomyopathy. However, the biomechanical effect of sheetlet sliding is not well-understood, which is the focus here. We conducted finite element simulations of the left ventricle (LV) coupled with a windkessel lumped parameter model to study sheetlet sliding, based on cardiac MRI of a healthy human subject, and modifications to account for hypertrophic and dilated geometric changes during cardiomyopathy remodeling. We modeled sheetlet sliding as a reduced shear stiffness in the sheet-normal direction and observed that (1) the diastolic sheetlet orientations must depart from alignment with the LV wall plane in order for sheetlet sliding to have an effect on cardiac function, that (2) sheetlet sliding modestly aided cardiac function of the healthy and dilated hearts, in terms of ejection fraction, stroke volume, and systolic pressure generation, but its effects were amplified during hypertrophic cardiomyopathy and diminished during dilated cardiomyopathy due to both sheetlet angle configuration and geometry, and that (3) where sheetlet sliding aided cardiac function, it increased tissue stresses, particularly in the myofibre direction. We speculate that sheetlet sliding is a tissue architectural adaptation to allow easier deformations of the LV walls so that LV wall stiffness will not hinder function, and to provide a balance between function and tissue stresses. A limitation here is that sheetlet sliding is modeled as a simple reduction in shear stiffness, without consideration of micro-scale sheetlet mechanics and dynamics.

Myocardial Biomechanics and the Consequent Differentially Expressed Genes of the Left Atrial Ligation Chick Embryonic Model of Hypoplastic Left Heart Syndrome

Left atrial ligation (LAL) of the chick embryonic heart is a model of the hypoplastic left heart syndrome (HLHS) where a purely mechanical intervention without genetic or pharmacological manipulation is employed to initiate cardiac malformation. It is thus a key model for understanding the biomechanical origins of HLHS. However, its myocardial mechanics and subsequent gene expressions are not well-understood. We performed finite element (FE) modeling and single-cell RNA sequencing to address this. 4D high-frequency ultrasound imaging of chick embryonic hearts at HH25 (ED 4.5) were obtained for both LAL and control. Motion tracking was performed to quantify strains. Image-based FE modeling was conducted, using the direction of the smallest strain eigenvector as the orientations of contractions, the Guccione active tension model and a Fung-type transversely isotropic passive stiffness model that was determined via micro-pipette aspiration. Single-cell RNA sequencing of left ventricle (LV) heart tissues was performed for normal and LAL embryos at HH30 (ED 6.5) and differentially expressed genes (DEG) were identified.After LAL, LV thickness increased by 33%, strains in the myofiber direction increased by 42%, while stresses in the myofiber direction decreased by 50%. These were likely related to the reduction in ventricular preload and underloading of the LV due to LAL. RNA-seq data revealed potentially related DEG in myocytes, including mechano-sensing genes (Cadherins, NOTCH1, etc.), myosin contractility genes (MLCK, MLCP, etc.), calcium signaling genes (PI3K, PMCA, etc.), and genes related to fibrosis and fibroelastosis (TGF-β, BMP, etc.). We elucidated the changes to the myocardial biomechanics brought by LAL and the corresponding changes to myocyte gene expressions. These data may be useful in identifying the mechanobiological pathways of HLHS.

Recent advances in quantifying the mechanobiology of cardiac development via computational modeling

Mechanical forces are essential for coordinating cardiac morphogenesis, but much remains to be discovered about the interactions between mechanical forces and the mechanotransduction pathways they activate. Due to the elaborate and fundamentally multi-physics and multi-scale nature of cardiac mechanobiology, a complete understanding requires multiple experimental and analytical techniques. We identify three fundamental tools used in the field to probe these interactions: high resolution imaging, genetic and molecular analysis, and computational modeling. In this review, we focus on computational modeling and present recent studies employing this tool to investigate the mechanobiological pathways involved with cardiac development. These works demonstrate that understanding the detailed spatial and temporal patterns of biomechanical forces is crucial to building a comprehensive understanding of mechanobiology during cardiac development, and that computational modeling is an effective and efficient tool for obtaining such detail. In this context, multidisciplinary studies combining all three tools present the most compelling results.

Contribution of Ventricular Motion and Sampling Location to Discrepancies in Two-Dimensional Versus Three-Dimensional Fetal Ventricular Strain Measures

BACKGROUND

Echocardiographic quantification of fetal cardiac strain is important to evaluate function and the need for intervention, with both two-dimensional (2D) and three-dimensional (3D) strain measurements currently feasible. However, discrepancies between 2D and 3D measurements have been reported, the etiologies of which are unclear. This study sought to determine the etiologies of the differences between 2D and 3D strain measurements.

METHODS

A validated cardiac motion-tracking algorithm was used on 3D cine ultrasound images acquired in 26 healthy fetuses. Both 2D and 3D myocardial strain quantifications were performed on each image set for controlled comparisons. Finite element modeling of 2 left ventricle (LV) models with minor geometrical differences were performed with various helix angle configurations for validating image processing results.

RESULTS

Three-dimensional longitudinal strain (LS) was significantly lower than 2D LS for the LV free wall and septum but not for the right ventricular (RV) free wall, while 3D circumferential strain (CS) was significantly higher than 2D CS for the LV, RV, and septum. The LS discrepancy was due to 2D long-axis imaging not capturing the out-of-plane motions associated with LV twist, while the CS discrepancy was due to the systolic motion of the heart toward the apex that caused out-of-plane motions in 2D short-axis imaging. A timing mismatch between the occurrences of peak longitudinal and circumferential dimensions caused a deviation in zero-strain referencing between 2D and 3D strain measurements, contributing to further discrepancies between the 2.

CONCLUSIONS

Mechanisms for discrepancies between 2D and 3D strain measurements in fetal echocardiography were identified, and inaccuracies associated with 2D strains were highlighted. Understanding of this mechanism is useful and important for future standardization of fetal cardiac strain measurements, which we propose to be important in view of large discrepancies in measured values in the literature.

The dependency of fetal left ventricular biomechanics function on myocardium helix angle configuration

The helix angle configuration of the myocardium is understood to contribute to the heart function, as finite element (FE) modeling of postnatal hearts showed that altered configurations affected cardiac function and biomechanics. However, similar investigations have not been done on the fetal heart. To address this, we performed image-based FE simulations of fetal left ventricles (LV) over a range of helix angle configurations, assuming a linear variation of helix angles from epicardium to endocardium. Results showed that helix angles have substantial influence on peak myofiber stress, cardiac stroke work, myocardial deformational burden, and spatial variability of myocardial strain. A good match between LV myocardial strains from FE simulations to those measured from 4D fetal echo images could only be obtained if the transmural variation of helix angle was generally between 110 and 130°, suggesting that this was the physiological range. Experimentally discovered helix angle configurations from the literature were found to produce high peak myofiber stress, high cardiac stroke work, and a low myocardial deformational burden, but did not coincide with configurations that would optimize these characteristics. This may suggest that the fetal development of myocyte orientations depends concurrently on several factors rather than a single factor. We further found that the shape, rather than the size of the LV, determined the manner at which helix angles influenced these characteristics, as this influence changed significantly when the LV shape was varied, but not when a heart was scaled from fetal to adult size while retaining the same shape. This may suggest that biomechanical optimality would be affected during diseases that altered the geometric shape of the LV.

Soft-Tissue Material Properties and Mechanogenetics during Cardiovascular Development

During embryonic development, changes in the cardiovascular microstructure and material properties are essential for an integrated biomechanical understanding. This knowledge also enables realistic predictive computational tools, specifically targeting the formation of congenital heart defects. Material characterization of cardiovascular embryonic tissue at consequent embryonic stages is critical to understand growth, remodeling, and hemodynamic functions. Two biomechanical loading modes, which are wall shear stress and blood pressure, are associated with distinct molecular pathways and govern vascular morphology through microstructural remodeling. Dynamic embryonic tissues have complex signaling networks integrated with mechanical factors such as stress, strain, and stiffness. While the multiscale interplay between the mechanical loading modes and microstructural changes has been studied in animal models, mechanical characterization of early embryonic cardiovascular tissue is challenging due to the miniature sample sizes and active/passive vascular components. Accordingly, this comparative review focuses on the embryonic material characterization of developing cardiovascular systems and attempts to classify it for different species and embryonic timepoints. Key cardiovascular components including the great vessels, ventricles, heart valves, and the umbilical cord arteries are covered. A state-of-the-art review of experimental techniques for embryonic material characterization is provided along with the two novel methods developed to measure the residual and von Mises stress distributions in avian embryonic vessels noninvasively, for the first time in the literature. As attempted in this review, the compilation of embryonic mechanical properties will also contribute to our understanding of the mature cardiovascular system and possibly lead to new microstructural and genetic interventions to correct abnormal development.

Deficient Myocardial Organization and Pathological Fibrosis in Fetal Aortic Stenosis-Association of Prenatal Ultrasound with Postmortem Histology

In fetal aortic stenosis (AS), it remains challenging to predict left ventricular development over the course of pregnancy. Myocardial organization, differentiation and fibrosis could be potential biomarkers relevant for biventricular outcome. We present four cases of fetal AS with varying degrees of severity and associate myocardial deformation on fetal ultrasound with postmortem histopathological characteristics. During routine fetal echocardiography, speckle tracking recordings of the cardiac four-chamber view were performed to assess myocardial strain as parameter for myocardial deformation. After pregnancy termination, postmortem cardiac specimens were examined using immunohistochemical labeling (IHC) of key markers for myocardial organization, differentiation and fibrosis and compared to normal fetal hearts. Two cases with critical AS presented extremely decreased left ventricular (LV) strain on fetal ultrasound. IHC showed overt endocardial fibro-elastosis, which correlated with pathological fibrosis patterns in the myocardium and extremely disturbed cardiomyocyte organization. The LV in severe AS showed mildly reduced myocardial strain and less severe disorganization of the cardiomyocytes. In conclusion, the degree of reduction in myocardial deformation corresponded with high extent to the amount of pathological fibrosis patterns and cardiomyocyte disorganization. Myocardial deformation on fetal ultrasound seems to hold promise as a potential biomarker for left ventricular structural damage in AS.

Mending a broken heart: In vitro, in vivo and in silico models of congenital heart disease

Birth defects contribute to ∼0.3% of global infant mortality in the first month of life, and congenital heart disease (CHD) is the most common birth defect among newborns worldwide. Despite the significant impact on human health, most treatments available for this heterogenous group of disorders are palliative at best. For this reason, the complex process of cardiogenesis, governed by multiple interlinked and dose-dependent pathways, is well investigated. Tissue, animal and, more recently, computerized models of the developing heart have facilitated important discoveries that are helping us to understand the genetic, epigenetic and mechanobiological contributors to CHD aetiology. In this Review, we discuss the strengths and limitations of different models of normal and abnormal cardiogenesis, ranging from single-cell systems and 3D cardiac organoids, to small and large animals and organ-level computational models. These investigative tools have revealed a diversity of pathogenic mechanisms that contribute to CHD, including genetic pathways, epigenetic regulators and shear wall stresses, paving the way for new strategies for screening and non-surgical treatment of CHD. As we discuss in this Review, one of the most-valuable advances in recent years has been the creation of highly personalized platforms with which to study individual diseases in clinically relevant settings.