Computational hemodynamic optimization predicts dominant aortic arch selection is driven by embryonic outflow tract orientation in the chick embryo

In the early embryo, a series of symmetric, paired vessels, the aortic arches, surround the foregut and distribute cardiac output to the growing embryo and fetus. During embryonic development, the arch vessels undergo large-scale asymmetric morphogenesis to form species-specific adult great vessel patterns. These transformations occur within a dynamic biomechanical environment, which can play an important role in the development of normal arch configurations or the aberrant arch morphologies associated with congenital cardiac defects. Arrested migration and rotation of the embryonic outflow tract during late stages of cardiac looping has been shown to produce both outflow tract and several arch abnormalities. Here, we investigate how changes in flow distribution due to a perturbation in the angular orientation of the embryonic outflow tract impact the morphogenesis and growth of the aortic arches. Using a combination of in vivo arch morphometry with fluorescent dye injection and hemodynamics-driven bioengineering optimization-based vascular growth modeling, we demonstrate that outflow tract orientation significantly changes during development and that the associated changes in hemodynamic load can dramatically influence downstream aortic arch patterning. Optimization reveals that balancing energy expenditure with diffusive capacity leads to multiple arch vessel patterns as seen in the embryo, while minimizing energy alone led to the single arch configuration seen in the mature arch of aorta. Our model further shows the critical importance of the orientation of the outflow tract in dictating morphogenesis to the adult single arch and accurately predicts arch IV as the dominant mature arch of aorta. These results support the hypothesis that abnormal positioning of the outflow tract during early cardiac morphogenesis may lead to congenital defects of the great vessels due to altered hemodynamic loading.     

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.     

The total cavopulmonary connection resistance: a significant impact on single ventricle hemodynamics at rest and exercise

Little is known about the impact of the total cavopulmonary connection (TCPC) on resting and exercise hemodynamics in a single ventricle (SV) circulation. The aim of this study was to elucidate this mechanism using a lumped parameter model of the SV circulation. Pulmonary vascular resistance (1.96+/-0.80 WU) and systemic vascular resistances (18.4+/-7.2 WU) were obtained from catheterization data on 40 patients with a TCPC. TCPC resistances (0.39+/-0.26 WU) were established using computational fluid dynamic simulations conducted on anatomically accurate three-dimensional models reconstructed from MRI (n=16). These parameters were used in a lumped parameter model of the SV circulation to investigate the impact of TCPC resistance on SV hemodynamics under resting and exercise conditions. A biventricular model was used for comparison. For a biventricular circulation, the cardiac output (CO) dependence on TCPC resistance was negligible (sensitivity=-0.064 l.min(-1).WU(-1)) but not for the SV circulation (sensitivity=-0.88 l.min(-1).WU(-1)). The capacity to increase CO with heart rate was also severely reduced for the SV. At a simulated heart rate of 150 beats/min, the SV patient with the highest resistance (1.08 WU) had a significantly lower increase in CO (20.5%) compared with the SV patient with the lowest resistance (50%) and normal circulation (119%). This was due to the increased afterload (+35%) and decreased preload (-12%) associated with the SV circulation. In conclusion, TCPC resistance has a significant impact on resting hemodynamics and the exercise capacity of patients with a SV physiology.     

Analysis of early embryonic great-vessel microcirculation in zebrafish using high-speed confocal μPIV

In the developing cardiovascular system, hemodynamic vascular loading is critical for angiogenesis and cardiovascular adaptation. Normal zebrafish embryos with transgenically-labeled endothelial and red blood cells provide an excellent in vivo model for studying the fluid-flow induced vascular loading. To characterize the developmental hemodynamics of early embryonic great-vessel microcirculation in the zebrafish embryo, two complementary studies (experimental and numerical) are presented. Quantitative comparison of the wall shear stress (WSS) at the first aortic arch (AA1) of wild-type zebrafish embryos during two consecutive developmental stages is presented, using time-resolved confocal micro-particle image velocimetry (μPIV). Analysis showed that there was significant WSS difference between 32 and 48 h post-fertilization (hpf) wild-type embryos, which correlates with normal arch morphogenesis. The vascular distensibility of the arch wall at systole and the acceleration/deceleration rates of time-lapse phase-averaged streamwise blood flow curves were also analyzed. To estimate the influence of a novel intermittent red-blood cell (RBC) loading on the endothelium, a numerical two-phase, volume of fluid (VOF) flow model was further developed with realistic in vivo conditions. These studies showed that near-wall effects and cell clustering increased WSS augmentation at a minimum of 15% when the distance of RBC from arch vessel wall was less than 3 μm or when RBC cell-to-cell distance was less than 3 μm. When compared to a smooth wall, the WSS augmentation increased by a factor of ~1.4 due to the roughness of the wall created by the endothelial cell profile. These results quantitatively highlight the contribution of individual RBC flow patterns on endothelial WSS in great-vessel microcirculation and will benefit the quantitative understanding of mechanotransduction in embryonic great vessel biology, including arteriovenous malformations (AVM).     

Single-step stereolithography of complex anatomical models for optical flow measurements

Transparent stereolithographic rapid prototyping (RP) technology has already demonstrated in literature to be a practical model construction tool for optical flow measurements such as digital particle image velocimetry (DPIV), laser doppler velocimetry (LDV), and flow visualization. Here, we employ recently available transparent RP resins and eliminate time-consuming casting and chemical curing steps from the traditional approach. This note details our methodology with relevant material properties and highlights its advantages. Stereolithographic model printing with our procedure is now a direct single-step process, enabling faster geometric replication of complex computational fluid dynamics (CFD) models for exact experimental validation studies. This methodology is specifically applied to the in vitro flow modeling of patient-specific total cavopulmonary connection (TCPC) morphologies. The effect of RP machining grooves, surface quality, and hydrodynamic performance measurements as compared with the smooth glass models are also quantified.     

Interaction between alk1 and blood flow in the development of arteriovenous malformations

Arteriovenous malformations (AVMs) are fragile direct connections between arteries and veins that arise during times of active angiogenesis. To understand the etiology of AVMs and the role of blood flow in their development, we analyzed AVM development in zebrafish embryos harboring a mutation in activin receptor-like kinase I (alk1), which encodes a TGFβ family type I receptor implicated in the human vascular disorder hereditary hemorrhagic telangiectasia type 2 (HHT2). Our analyses demonstrate that increases in arterial caliber, which stem in part from increased cell number and in part from decreased cell density, precede AVM development, and that AVMs represent enlargement and stabilization of normally transient arteriovenous connections. Whereas initial increases in endothelial cell number are independent of blood flow, later increases, as well as AVMs, are dependent on flow. Furthermore, we demonstrate that alk1 expression requires blood flow, and despite normal levels of shear stress, some flow-responsive genes are dysregulated in alk1 mutant arterial endothelial cells. Taken together, our results suggest that Alk1 plays a role in transducing hemodynamic forces into a biochemical signal required to limit nascent vessel caliber, and support a novel two-step model for HHT-associated AVM development in which pathological arterial enlargement and consequent altered blood flow precipitate a flow-dependent adaptive response involving retention of normally transient arteriovenous connections, thereby generating AVMs.     

Spatiotemporal remodeling of embryonic aortic arch: stress distribution,microstructure, and vascular growth in silico

Biomechanics and Modeling in Mechanobiology - The microstructure for mature vessels has been investigated in detail, while there is limited information about the embryonic stages, in spite of their...     

Characterizaton of the Vessel Geometry,Flow Mechanics and Wall Shear Stress in the Great Arteries of Wildtype Prenatal Mouse

Introduction

Abnormal fluid mechanical environment in the pre-natal cardiovascular system is hypothesized to play a significant role in causing structural heart malformations. It is thus important to improve our understanding of the prenatal cardiovascular fluid mechanical environment at multiple developmental time-points and vascular morphologies. We present such a study on fetal great arteries on the wildtype mouse from embryonic day 14.5 (E14.5) to near-term (E18.5).

Methods

Ultrasound bio-microscopy (UBM) was used to measure blood velocity of the great arteries. Subsequently, specimens were cryo-embedded and sectioned using episcopic fluorescent image capture (EFIC) to obtain high-resolution 2D serial image stacks, which were used for 3D reconstructions and quantitative measurement of great artery and aortic arch dimensions. EFIC and UBM data were input into subject-specific computational fluid dynamics (CFD) for modeling hemodynamics.

Results

In normal mouse fetuses between E14.5–18.5, ultrasound imaging showed gradual but statistically significant increase in blood velocity in the aorta, pulmonary trunk (with the ductus arteriosus), and descending aorta. Measurement by EFIC imaging displayed a similar increase in cross sectional area of these vessels. However, CFD modeling showed great artery average wall shear stress and wall shear rate remain relatively constant with age and with vessel size, indicating that hemodynamic shear had a relative constancy over gestational period considered here.

Conclusion

Our EFIC-UBM-CFD method allowed reasonably detailed characterization of fetal mouse vascular geometry and fluid mechanics. Our results suggest that a homeostatic mechanism for restoring vascular wall shear magnitudes may exist during normal embryonic development. We speculate that this mechanism regulates the growth of the great vessels.     

Oscillating Couette flow for in vitro cell loading

Synovial joints are loaded by weight bearing, stretching, and fluid-driven shear. To simulate in vitro fluid-driven shear, we developed an "oscillating Couette flow mechanical shear loader". Oscillating Couette flow mimics relative motion of articular surfaces; hence, characterizing flow-induced shear by the loader enhances understanding of mechanotransduction in the joint tissue. Here, the analytical and computational models for an oscillating Couette flow were used to predict time-varying shear distribution on a plate surface, applying numerical simulation to evaluate the effects of finite plate dimension in a 2D flow. Shear stress on the plate was significantly different from that in simpler models (unbounded plates and viscous low-frequency flow). High-stress spots appeared near the leading and trailing edges of a moving plate, and a relatively uniform shear region was restricted to the interior area. Stress prediction in an example experimental geometry is presented, where the frequency and finite width effects are feasibly accounted.