Flow in the early embryonic human heart: a numerical study

Hemodynamic and Structural Comparison of Human Fetal Heart Development Between Normally Growing and Hypoplastic Left Heart Syndrome-Diagnosed Hearts

Congenital heart defects (CHDs) affect a wide range of societies with an incidence rate of 1.0–1.2%. These defects initiate at the early developmental stage and result in critical health disorders. Although genetic factors play a role in the formation of CHDs, the occurrence of cases in families with no history of CHDs suggests that mechanobiological forces may also play a role in the initiation and progression of CHDs. Hypoplastic left heart syndrome (HLHS) is a critical CHD, which is responsible for 25–40% of all prenatal cardiac deaths. The comparison of healthy and HLHS hearts helps in understanding the main hemodynamic differences related to HLHS. Echocardiography is the most common imaging modality utilized for fetal cardiac assessment. In this study, we utilized echocardiographic images to compare healthy and HLHS human fetal hearts for determining the differences in terms of heart chamber dimensions, valvular flow rates, and hemodynamics. The cross-sectional areas of chamber dimensions are determined from 2D b-mode ultrasound images. Valvular flow rates are measured via Doppler echocardiography, and hemodynamic quantifications are performed with the use of computational fluid dynamics (CFD) simulations. The obtained results indicate that cross-sectional areas of the left and right sides of the heart are similar for healthy fetuses during gestational development. The left side of HLHS heart is underdeveloped, and as a result, the hemodynamic parameters such as flow velocity, pressure, and wall shear stress (WSS) are significantly altered compared to those of healthy hearts.

Retrograde flow in aortic isthmus in normal and fetal heart disease by principal component analysis and computational fluid dynamics

OBJECTIVES

Reverse flow Retrograde flow (RF) of blood in the aortic isthmus can be observed in different types of fetal heart disease (FHD), including abnormalities in heart structure and function. This study sought to investigate the relationship between RF and blood flow parameters, and develop a computational fluid dynamics (CFD) model to understand the mechanisms underlying this observation.

MATERIAL AND METHODS

A total of 281 fetuses (gestational age [GA] 26.6±.3 weeks) with FHD and 2803 normal fetuses (GA: 26.1±.1 weeks) by fetal echocardiography collected from May 2016 to December 2018. Principal component analysis (PCA) was performed to find the relationship and the CFD model reconstructed from 3D/4D spatio-temporal image correlation (STIC) images to simulate hemodynamics.

RESULTS

There was a significant difference in the percentages of RF between the study (80/201 (39%)) and control (29/2803 (1%)) groups (p < 0.05). The RF occur when the aorta flow rate (left heart) is reduced to 60% by CFD stimulation. Pearson correlation analysis showed significant correlations between flow rate and wall shear stress(WSS) (r = .883, p = 0.047) variables at the AI.

CONCLUSION

Volumetric flow rate of AO or left heart was the main component of the cause of RF. The hemodynamics of the cardiovascular system have highly complex behavior hinge on the turbulent nature of circulating blood flow.

Vortex Dynamics in Trabeculated Embryonic Ventricles

Proper heart morphogenesis requires a delicate balance between hemodynamic forces, myocardial activity, morphogen gradients, and epigenetic signaling, all of which are coupled with genetic regulatory networks. Recently both in vivo and in silico studies have tried to better understand hemodynamics at varying stages of veretebrate cardiogenesis. In particular, the intracardial hemodynamics during the onset of trabeculation is notably complex—the inertial and viscous fluid forces are approximately equal at this stage and small perturbations in morphology, scale, and steadiness of the flow can lead to significant changes in bulk flow structures, shear stress distributions, and chemical morphogen gradients. The immersed boundary method was used to numerically simulate fluid flow through simplified two-dimensional and stationary trabeculated ventricles of 72, 80, and 120 h post fertilization wild type zebrafish embryos and ErbB2-inhibited embryos at seven days post fertilization. A 2D idealized trabeculated ventricular model was also used to map the bifurcations in flow structure that occur as a result of the unsteadiness of flow, trabeculae height, and fluid scale (Re). Vortex formation occurred in intertrabecular regions for biologically relevant parameter spaces, wherein flow velocities increased. This indicates that trabecular morphology may alter intracardial flow patterns and hence ventricular shear stresses and morphogen gradients. A potential implication of this work is that the onset of vortical (disturbed) flows can upregulate Notch1 expression in endothelial cells in vivo and hence impacts chamber morphogenesis, valvulogenesis, and the formation of the trabeculae themselves. Our results also highlight the sensitivity of cardiac flow patterns to changes in morphology and blood rheology, motivating efforts to obtain spatially and temporally resolved chamber geometries and kinematics as well as the careful measurement of the embryonic blood rheology. The results also suggest that there may be significant changes in shear signalling due to morphological and mechanical variation across individuals and species.

Fluid dynamics in heart development: effects of hematocrit and trabeculation

Recent in vivo experiments have illustrated the importance of understanding the haemodynamics of heart morphogenesis. In particular, ventricular trabeculation is governed by a delicate interaction between haemodynamic forces, myocardial activity, and morphogen gradients, all of which are coupled to genetic regulatory networks. The underlying haemodynamics at the stage of development in which the trabeculae form is particularly complex, given the balance between inertial and viscous forces. Small perturbations in the geometry, scale, and steadiness of the flow can lead to changes in the overall flow structures and chemical morphogen gradients, including the local direction of flow, the transport of morphogens, and the formation of vortices. The immersed boundary method was used to solve the two-dimensional fluid-structure interaction problem of fluid flow moving through a two chambered heart of a zebrafish (Danio rerio), with a trabeculated ventricle, at 96 hours post fertilization (hpf). Trabeculae heights and hematocrit were varied, and simulations were conducted for two orders of magnitude of Womersley number, extending beyond the biologically relevant range (0.2-12.0). Both intracardial and intertrabecular vortices formed in the ventricle for biologically relevant parameter values. The bifurcation from smooth streaming flow to vortical flow depends upon the trabeculae geometry, hematocrit, and Womersley number, $Wo$. This work shows the importance of hematocrit and geometry in determining the bulk flow patterns in the heart at this stage of development.

Left atrial ligation alters intracardiac flow patterns and the biomechanical landscape in the chick embryo

BACKGROUND

Hypoplastic left heart syndrome (HLHS) is a major human congenital heart defect that results in single ventricle physiology and high mortality. Clinical data indicate that intracardiac blood flow patterns during cardiac morphogenesis are a significant etiology. We used the left atrial ligation (LAL) model in the chick embryo to test the hypothesis that LAL immediately alters intracardiac flow streams and the biomechanical environment, preceding morphologic and structural defects observed in HLHS.

RESULTS

Using fluorescent dye injections, we found that intracardiac flow patterns from the right common cardinal vein, right vitelline vein, and left vitelline vein were altered immediately following LAL. Furthermore, we quantified a significant ventral shift of the right common cardinal and right vitelline vein flow streams. We developed an in silico model of LAL, which revealed that wall shear stress was reduced at the left atrioventricular canal and left side of the common ventricle.

CONCLUSIONS

Our results demonstrate that intracardiac flow patterns change immediately following LAL, supporting the role of hemodynamics in the progression of HLHS. Sites of reduced WSS revealed by computational modeling are commonly affected in HLHS, suggesting that changes in the biomechanical environment may lead to abnormal growth and remodeling of left heart structures.

Biomechanics of early cardiac development

Biomechanics affect early cardiac development, from looping to the development of chambers and valves. Hemodynamic forces are essential for proper cardiac development, and their disruption leads to congenital heart defects. A wealth of information already exists on early cardiac adaptations to hemodynamic loading, and new technologies, including high-resolution imaging modalities and computational modeling, are enabling a more thorough understanding of relationships between hemodynamics and cardiac development. Imaging and modeling approaches, used in combination with biological data on cell behavior and adaptation, are paving the road for new discoveries on links between biomechanics and biology and their effect on cardiac development and fetal programming.

Computational fluid dynamics of developing avian outflow tract heart valves

Hemodynamic forces play an important role in sculpting the embryonic heart and its valves. Alteration of blood flow patterns through the hearts of embryonic animal models lead to malformations that resemble some clinical congenital heart defects, but the precise mechanisms are poorly understood. Quantitative understanding of the local fluid forces acting in the heart has been elusive because of the extremely small and rapidly changing anatomy. In this study, we combine multiple imaging modalities with computational simulation to rigorously quantify the hemodynamic environment within the developing outflow tract (OFT) and its eventual aortic and pulmonary valves. In vivo Doppler ultrasound generated velocity profiles were applied to Micro-Computed Tomography generated 3D OFT lumen geometries from Hamburger-Hamilton (HH) stage 16-30 chick embryos. Computational fluid dynamics simulation initial conditions were iterated until local flow profiles converged with in vivo Doppler flow measurements. Results suggested that flow in the early tubular OFT (HH16 and HH23) was best approximated by Poiseuille flow, while later embryonic OFT septation (HH27, HH30) was mimicked by plug flow conditions. Peak wall shear stress (WSS) values increased from 18.16 dynes/cm(2) at HH16 to 671.24 dynes/cm(2) at HH30. Spatiotemporally averaged WSS values also showed a monotonic increase from 3.03 dynes/cm(2) at HH16 to 136.50 dynes/cm(2) at HH30. Simulated velocity streamlines in the early heart suggest a lack of mixing, which differed from classical ink injections. Changes in local flow patterns preceded and correlated with key morphogenetic events such as OFT septation and valve formation. This novel method to quantify local dynamic hemodynamics parameters affords insight into sculpting role of blood flow in the embryonic heart and provides a quantitative baseline dataset for future research.

Fluid dynamics of heart development

The morphology, muscle mechanics, fluid dynamics, conduction properties, and molecular biology of the developing embryonic heart have received much attention in recent years due to the importance of both fluid and elastic forces in shaping the heart as well as the striking relationship between the heart's evolution and development. Although few studies have directly addressed the connection between fluid dynamics and heart development, a number of studies suggest that fluids may play a key role in morphogenic signaling. For example, fluid shear stress may trigger biochemical cascades within the endothelial cells of the developing heart that regulate chamber and valve morphogenesis. Myocardial activity generates forces on the intracardiac blood, creating pressure gradients across the cardiac wall. These pressures may also serve as epigenetic signals. In this article, the fluid dynamics of the early stages of heart development is reviewed. The relevant work in cardiac morphology, muscle mechanics, regulatory networks, and electrophysiology is also reviewed in the context of intracardial fluid dynamics.

Quantifying blood flow and wall shear stresses in the outflow tract of chick embryonic hearts

Wall shear stresses (WSS) exerted by blood flow on cardiac tissues modulate growth and development of the heart. To study the role of hemodynamic conditions on cardiac morphogenesis, here, we present a methodology that combines imaging and finite element modeling to quantify the in vivo blood flow dynamics and WSS in the cardiac outflow tract (OFT) of early chicken embryos (day 3 out of 21-day incubation period). We found a distinct blood flow field and heterogeneous distribution of WSS in the chicken embryonic heart OFT, which have physiological implications for OFT morphogenesis.

Hemodynamic patterning of the avian atrioventricular valve

In this study, we develop an innovative approach to rigorously quantify the evolving hemodynamic environment of the atrioventricular (AV) canal of avian embryos. Ultrasound generated velocity profiles were imported into Micro-Computed Tomography generated anatomically precise cardiac geometries between Hamburger-Hamilton (HH) stages 17 and 30. Computational fluid dynamic simulations were then conducted and iterated until results mimicked in vivo observations. Blood flow in tubular hearts (HH17) was laminar with parallel streamlines, but strong vortices developed simultaneous with expansion of the cushions and septal walls. For all investigated stages, highest wall shear stresses (WSS) are localized to AV canal valve-forming regions. Peak WSS increased from 19.34 dynes/cm(2) at HH17 to 287.18 dynes/cm(2) at HH30, but spatiotemporally averaged WSS became 3.62 dynes/cm(2) for HH17 to 9.11 dynes/cm(2) for HH30. Hemodynamic changes often preceded and correlated with morphological changes. These results establish a quantitative baseline supporting future hemodynamic analyses and interpretations.

The development of the heart and microcirculation: role of shear stress

It is evident that hemodynamic factors have a dominant function already during early cardiogenesis. Flow and ensuing shear stress are sensed by endothelial cells by, ciliary modified, cytoskeletal deformation which then activates a number of subcellular structures and molecules. Shear stress dependent changes mostly converge towards NF kappa B signaling and DNA binding, thereby altering metabolic paths and influencing differentiation of the cells. Geometry of the vascular system heavily affects the flow and shear patterns, as is the case in the adult vasculature where atheroprone areas nicely coincide with the frequency of the primary cilium as shear stress sensor.

In vitro hemodynamic investigation of the embryonic aortic arch at late gestation

This study focuses on the dynamic flow through the fetal aortic arch driven by the concurrent action of right and left ventricles. We created a parametric pulsatile computational fluid dynamics (CFD) model of the fetal aortic junction with physiologic vessel geometries. To gain a better biophysical understanding, an in vitro experimental fetal flow loop for flow visualization was constructed for identical CFD conditions. CFD and in vitro experimental results were comparable. Swirling flow during the acceleration phase of the cardiac cycle and unidirectional flow following mid-deceleration phase were observed in pulmonary arteries (PA), head-neck vessels, and descending aorta. Right-to-left (oxygenated) blood flowed through the ductus arteriosus (DA) posterior relative to the antegrade left ventricular outflow tract (LVOT) stream and resembled jet flow. LVOT and right ventricular outflow tract flow mixing had not completed until approximately 3.5 descending aorta diameters downstream of the DA insertion into the aortic arch. Normal arch model flow patterns were then compared to flow patterns of four common congenital heart malformations that include aortic arch anomalies. Weak oscillatory reversing flow through the DA junction was observed only for the Tetralogy of Fallot configuration. PA and hypoplastic left heart syndrome configurations demonstrated complex, abnormal flow patterns in the PAs and head-neck vessels. Aortic coarctation resulted in large-scale recirculating flow in the aortic arch proximal to the DA. Intravascular flow patterns spatially correlated with abnormal vascular structures consistent with the paradigm that abnormal intravascular flow patterns associated with congenital heart disease influence vascular growth and function.

Computerized three-dimensional analysis of the heart and great vessels in normal and holoprosencephalic human embryos

The developing heart and great vessels undergo drastic morphogenetic changes during the embryonic period. To analyze the normal and abnormal development of these organs, it is essential to visualize their structures in three and four dimensions, including the changes occurring with time. We have reconstructed the luminal structure of the hearts and great vessels of staged human embryos from serial histological sections to demonstrate their sequential morphological changes in three dimensions. The detailed structures of the embryonic heart and major arteries in normal and holoprosencephalic (HPE) human embryos could be reconstructed and visualized, and anatomical structures were analyzed using 3D images. By 3D analysis, cardiac anomalies such as double-outlet right ventricle and malrotation of the heart tube were identified in HPE embryos, which were not easily diagnosed by histological observation. Reconstruction and analysis of 3D images are useful for the study of anatomical structures of developing embryos and for identifying their abnormalities.

A driving force for change: interstitial flow as a morphoregulator

Dynamic stresses that are present in all living tissues drive small fluid flows, called interstitial flows, through the extracellular matrix. Interstitial flow not only helps to transport nutrients throughout the tissue, but also has important roles in tissue maintenance and pathobiology that have been, until recently, largely overlooked. Here, we present evidence for the various effects of interstitial flow on cell biology, including its roles in embryonic development, tissue morphogenesis and remodeling, inflammation and lymphedema, tumor biology and immune cell trafficking. We also discuss possible mechanisms by which interstitial flow can induce morphoregulation, including direct shear stress, matrix-cell transduction (as has been proposed in the endothelial glycocalyx) and the newly emerging concept of autologous gradient formation.

Physics-Driven CFD Modeling of Complex Anatomical Cardiovascular Flows?A TCPC Case Study

Recent developments in medical image acquisition combined with the latest advancements in numerical methods for solving the Navier-Stokes equations have created unprecedented opportunities for developing simple and reliable computational fluid dynamics (CFD) tools for meeting patient-specific surgical planning objectives. However, for CFD to reach its full potential and gain the trust and confidence of medical practitioners, physics-driven numerical modeling is required. This study reports on the experience gained from an ongoing integrated CFD modeling effort aimed at developing an advanced numerical simulation tool capable of accurately predicting flow characteristics in an anatomically correct total cavopulmonary connection (TCPC). An anatomical intra-atrial TCPC model is reconstructed from a stack of magnetic resonance (MR) images acquired in vivo. An exact replica of the computational geometry was built using transparent rapid prototyping. Following the same approach as in earlier studies on idealized models, flow structures, pressure drops, and energy losses were assessed both numerically and experimentally, then compared. Numerical studies were performed with both a first-order accurate commercial software and a recently developed, second-order accurate, in-house flow solver. The commercial CFD model could, with reasonable accuracy, capture global flow quantities of interest such as control volume power losses and pressure drops and time-averaged flow patterns. However, for steady inflow conditions, both flow visualization experiments and particle image velocimetry (PIV) measurements revealed unsteady, complex, and highly 3D flow structures, which could not be captured by this numerical model with the available computational resources and additional modeling efforts that are described. Preliminary time-accurate computations with the in-house flow solver were shown to capture for the first time these complex flow features and yielded solutions in good agreement with the experimental observations. Flow fields obtained were similar for the studied total cardiac output range (1-3 1/min); however hydrodynamic power loss increased dramatically with increasing cardiac output, suggesting significant energy demand at exercise conditions. The simulation of cardiovascular flows poses a formidable challenge to even the most advanced CFD tools currently available. A successful prediction requires a two-pronged, physics-based approach, which integrates high-resolution CFD tools and high-resolution laboratory measurements.