Cardiac motion estimation from medical images: a regularisation framework applied on pairwise image registration displacement fields

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.

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.

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.

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.

Deep Learning-Based Image Registration in Dynamic Myocardial Perfusion CT Imaging

Registration of dynamic CT image sequences is a crucial preprocessing step for clinical evaluation of multiple physiological determinants in the heart such as global and regional myocardial perfusion. In this work, we present a deformable deep learning-based image registration method for quantitative myocardial perfusion CT examinations, which in contrast to previous approaches, takes into account some unique challenges such as low image quality with less accurate anatomical landmarks, dynamic changes of contrast agent concentration in the heart chambers and tissue, and misalignment caused by cardiac stress, respiration, and patient motion. The introduced method uses a recursive cascade network with a ventricle segmentation module, and a novel loss function that accounts for local contrast changes over time. It was trained and validated on a dataset of n = 118 patients with known or suspected coronary artery disease and/or aortic valve insufficiency. Our results demonstrate that the proposed method is capable of registering dynamic cardiac perfusion sequences by reducing local tissue displacements of the left ventricle (LV), whereas contrast changes do not affect the registration and image quality, in particular the absolute CT (HU) values of the entire CT sequence. In addition, the deep learning-based approach presented reveals a short processing time of a few seconds compared to conventional image registration methods, demonstrating its application potential for quantitative CT myocardial perfusion measurements in daily clinical routine.

Fluid Mechanical Effects of Fetal Aortic Valvuloplasty for Cases of Critical Aortic Stenosis with Evolving Hypoplastic Left Heart Syndrome

Fetuses with critical aortic stenosis (FAS) are at high risk of progression to HLHS by the time of birth (and are thus termed "evolving HLHS"). An in-utero catheter-based intervention, fetal aortic valvuloplasty (FAV), has shown promise as an intervention strategy to circumvent the progression, but its impact on the heart's biomechanics is not well understood. We performed patient-specific computational fluid dynamic (CFD) simulations based on 4D fetal echocardiography to assess the changes in the fluid mechanical environment in the FAS left ventricle (LV) directly before and 2 days after FAV. Echocardiograms of five FAS cases with technically successful FAV were retrospectively analysed. FAS compromised LV stroke volume and ejection fraction, but FAV rescued it significantly. Calculations to match simulations to clinical measurements showed that FAV approximately doubled aortic valve orifice area, but it remained much smaller than in healthy hearts. Diseased LVs had mildly stenotic mitral valves, which generated fast and narrow diastolic mitral inflow jet and vortex rings that remained unresolved directly after FAV. FAV further caused aortic valve damage and high-velocity regurgitation. The high-velocity aortic regurgitation jet and vortex ring caused a chaotic flow field upon impinging the apex, which drastically exacerbated the already high energy losses and poor flow energy efficiency of FAS LVs. Two days after the procedure, FAV did not alter wall shear stress (WSS) spatial patterns of diseased LV but elevated WSS magnitudes, and the poor blood turnover in pre-FAV LVs did not significantly improve directly after FAV. FAV improved FAS LV's flow function, but it also led to highly chaotic flow patterns and excessively high energy losses due to the introduction of aortic regurgitation directly after the intervention. Further studies analysing the effects several weeks after FAV are needed to understand the effects of such biomechanics on morphological development.

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.

Cardio-respiratory motion-corrected 3D cardiac water-fat MRI using model-based image reconstruction

PURPOSE

Myocardial fat infiltrations are associated with a range of cardiomyopathies. The purpose of this study was to perform cardio-respiratory motion-correction for model-based water-fat separation to image fatty infiltrations of the heart in a free-breathing, non-cardiac-triggered high-resolution 3D MRI acquisition.

METHODS

Data were acquired in nine patients using a free-breathing, non-cardiac-triggered high-resolution 3D Dixon gradient-echo sequence and radial phase encoding trajectory. Motion correction was combined with a model-based water-fat reconstruction approach. Respiratory and cardiac motion models were estimated using a dual-mode registration algorithm incorporating both motion-resolved water and fat information. Qualitative comparisons of fat structures were made between 2D clinical routine reference scans and reformatted 3D motion-corrected images. To evaluate the effect of motion correction the local sharpness of epicardial fat structures was analyzed for motion-averaged and motion-corrected fat images.

RESULTS

The reformatted 3D motion-corrected reconstructions yielded qualitatively comparable fat structures and fat structure sharpness in the heart as the standard 2D breath-hold. Respiratory motion correction improved the local sharpness on average by 32% ± 24% with maximum improvements of 81% and cardiac motion correction increased the sharpness further by another 15% ± 11% with maximum increases of 31%. One patient showed a fat infiltration in the myocardium and cardio-respiratory motion correction was able to improve its visualization in 3D.

CONCLUSION

The 3D water-fat separated cardiac images were acquired during free-breathing and in a clinically feasible and predictable scan time. Compared to a motion-averaged reconstruction an increase in sharpness of fat structures by 51% ± 27% using the presented motion correction approach was observed for nine patients.

Fluid Mechanics of Fetal Left Ventricle During Aortic Stenosis with Evolving Hypoplastic Left Heart Syndrome

In cases of fetal aortic stenosis and evolving Hypoplastic Left Heart Syndrome (feHLHS), aortic stenosis is associated with specific abnormalities such as retrograde or bidirectional systolic transverse arch flow. Many cases progressed to hypoplastic left heart syndrome (HLHS) malformation at birth, but fetal aortic valvuloplasty can prevent the progression in many cases. Since both disease and intervention involve drastic changes to the biomechanical environment, in-vivo biomechanics likely play a role in inducing and preventing disease progression. However, the fluid mechanics of feHLHS is not well-characterized. Here, we conduct patient-specific echocardiography-based flow simulations of normal and feHLHS left ventricles (LV), to understand the essential fluid dynamics distinction between the two cohorts. We found high variability across feHLHS cases, but also the following unifying features. Firstly, feHLHS diastole mitral inflow was in the form of a narrowed and fast jet that impinged onto the apical region, rather than a wide and gentle inflow in normal LVs. This was likely due to a malformed mitral valve with impaired opening dynamics. This altered inflow caused elevated vorticity dynamics and wall shear stresses (WSS) and reduced oscillatory shear index at the apical zone rather than mid-ventricle. Secondly, feHLHS LV also featured elevated systolic and diastolic energy losses, intraventricular pressure gradients, and vortex formation numbers, suggesting energy inefficiency of flow and additional burden on the LV. Thirdly, feHLHS LV had poor blood turnover, suggesting a hypoxic environment, which could be associated with endocardial fibroelastosis that is often observed in these patients.

Fluid mechanics of the zebrafish embryonic heart trabeculation

Embryonic heart development is a mechanosensitive process, where specific fluid forces are needed for the correct development, and abnormal mechanical stimuli can lead to malformations. It is thus important to understand the nature of embryonic heart fluid forces. However, the fluid dynamical behaviour close to the embryonic endocardial surface is very sensitive to the geometry and motion dynamics of fine-scale cardiac trabecular surface structures. Here, we conducted image-based computational fluid dynamics (CFD) simulations to quantify the fluid mechanics associated with the zebrafish embryonic heart trabeculae. To capture trabecular geometric and motion details, we used a fish line that expresses fluorescence at the endocardial cell membrane, and high resolution 3D confocal microscopy. Our endocardial wall shear stress (WSS) results were found to exceed those reported in existing literature, which were estimated using myocardial rather than endocardial boundaries. By conducting simulations of single intra-trabecular spaces under varied scenarios, where the translational or deformational motions (caused by contraction) were removed, we found that a squeeze flow effect was responsible for most of the WSS magnitude in the intra-trabecular spaces, rather than the shear interaction with the flow in the main ventricular chamber. We found that trabecular structures were responsible for the high spatial variability of the magnitude and oscillatory nature of WSS, and for reducing the endocardial deformational burden. We further found cells attached to the endocardium within the intra-trabecular spaces, which were likely embryonic hemogenic cells, whose presence increased endocardial WSS. Overall, our results suggested that a complex multi-component consideration of both anatomic features and motion dynamics were needed to quantify the trabeculated embryonic heart fluid mechanics.

In the embryonic heart, the mechanical forces that blood fluid imposes on the cardiac tissues are known to be important biological stimuli that affect the proper heart development. We thus perform careful quantification of these forces, using the zebrafish embryo as a model. To do this, we perform high resolution imaging of zebrafish embryonic hearts and image-based flow simulations. We find that the use of a particular fish line that expresses fluorescence at the exact boundary between heart tissue and blood, that is the endocardial cell membrane boundary, is important to give high quality results. The heart’s inner surface has uneven trabeculation structures. We find that they cause fluid forces to have spatial variability and an oscillatory nature. We also find that there is a squeezing motion of cardiac tissues on the trabeculation fluid spaces, which is the main mechanism that generated fluid forces. Fluid forces are also affected by a number of cardiac cells that were developing into blood cells, lodged in the trabeculation fluid spaces. Our investigations provide an understanding of the complexity of the fluid forces on the inner surface of the embryonic heart, and our quantifications will be useful to future studies on the biology elicited by these fluid forces.

Spatiotemporal Bayesian Regularization for Cardiac Strain Imaging: Simulation and In Vivo Results

Cardiac strain imaging (CSI) plays a critical role in the detection of myocardial motion abnormalities. Displacement estimation is an important processing step to ensure the accuracy and precision of derived strain tensors. In this paper, we propose and implement Spatiotemporal Bayesian regularization (STBR) algorithms for two-dimensional (2-D) normalized cross-correlation (NCC) based multi-level block matching along with incorporation into a Lagrangian cardiac strain estimation framework. Assuming smooth temporal variation over a short span of time, the proposed STBR algorithm performs displacement estimation using at least four consecutive ultrasound radio-frequency (RF) frames by iteratively regularizing 2-D NCC matrices using information from a local spatiotemporal neighborhood in a Bayesian sense. Two STBR schemes are proposed to construct Bayesian likelihood functions termed as Spatial then Temporal Bayesian (STBR-1) and simultaneous Spatiotemporal Bayesian (STBR-2). Radial and longitudinal strain estimated from a finite-element-analysis (FEA) model of realistic canine myocardial deformation were utilized to quantify strain bias, normalized strain error and total temporal relative error (TTR). Statistical analysis with one-way analysis of variance (ANOVA) showed that all Bayesian regularization methods significantly outperform NCC with lower bias and errors (p < 0.001). However, there was no significant difference among Bayesian methods. For example, mean longitudinal TTR for NCC, SBR, STBR-1 and STBR-2 were 25.41%, 9.27%, 10.38% and 10.13% respectively An in vivo feasibility study using RF data from ten healthy mice hearts were used to compare the elastographic signal-to-noise ratio (SNR_e) calculated using stochastic analysis. STBR-2 had the highest expected SNR_e both for radial and longitudinal strain. The mean expected SNR_e values for accumulated radial strain for NCC, SBR, STBR-1 and STBR-2 were 5.03, 9.43, 9.42 and 10.58, respectively. Overall results suggest that STBR improves CSI in vivo.

Cardiac forces regulate zebrafish heart valve delamination by modulating Nfat signaling

In the clinic, most cases of congenital heart valve defects are thought to arise through errors that occur after the endothelial–mesenchymal transition (EndoMT) stage of valve development. Although mechanical forces caused by heartbeat are essential modulators of cardiovascular development, their role in these later developmental events is poorly understood. To address this question, we used the zebrafish superior atrioventricular valve (AV) as a model. We found that cellularized cushions of the superior atrioventricular canal (AVC) morph into valve leaflets via mesenchymal–endothelial transition (MEndoT) and tissue sheet delamination. Defects in delamination result in thickened, hyperplastic valves, and reduced heart function. Mechanical, chemical, and genetic perturbation of cardiac forces showed that mechanical stimuli are important regulators of valve delamination. Mechanistically, we show that forces modulate Nfatc activity to control delamination. Together, our results establish the cellular and molecular signature of cardiac valve delamination in vivo and demonstrate the continuous regulatory role of mechanical forces and blood flow during valve formation.

Why do developing zebrafish atrioventricular heart valves become hyperplastic under certain hemodynamic conditions? This study suggests that part of the answer lies in how the mechanosensitive Nfat pathway regulates the valve mesenchymal-to-endothelial transition.

Full cardiac cycle asynchronous temporal compounding of 3D echocardiography images

It is important to improve echocardiography image quality, because the accuracy of echocardiographic assessment and diagnosis relies on image quality. Previous work on 2D temporal image compounding for image frames with matching cardiac phases (synchronous), and for temporally neighbouring image frames (asynchronous) over small ranges of time frames showed good improvement to image quality. Here, we extend this by performing asynchronous temporal compounding to echocardiographic images in 3D, involving all frames within a cardiac cycle, via a robust 3D cardiac motion estimation algorithm to describe the large image deformations. After compounding, the images can be reanimated via the motion model. Various methods of fusing image frames together are tested, including mean, max, and wavelet methods, and outlier rejection algorithms. The compounding algorithm is applied on 3D human adult, porcine adolescent, and human fetal echocardiography images. Results show significant improvements to contrast-to-noise ratio (CNR) and boundary clarity, and significantly decreased variability in manual quantification of cardiac chamber volumes after compounding. Interestingly, compounding can extend the field of view of the echo images, by reconstructing cardiac structures that momentarily exceeded the field of view, using the motion estimation algorithm to calculate their locations outside the field of view during these time periods. Although all compounding methods provide general improvements, the mean method led to blurred boundaries, while the max methods led to high variability of CNR. Outlier rejection algorithms were found to be useful in addressing these weaknesses.

Fluid mechanics of the left atrial ligation chick embryonic model of hypoplastic left heart syndrome

Left atrial ligation (LAL) of the chick embryonic heart at HH21 is a model of the hypoplastic left heart syndrome (HLHS) disease, demonstrating morphological and hemodynamic features similar to human HLHS cases. Since it relies on mechanical intervention without genetic or pharmacological manipulations, it is a good model for understanding the biomechanics origins of such HLHS malformations. To date, however, the fluid mechanical environment of this model is poorly understood. In the current study, we performed 4D ultrasound imaging of LAL and normal chick embryonic hearts and 4D cardiac flow simulations to help shed light on the mechanical environment that may lead to the HLHS morphology. Results showed that the HH25 LAL atrial function was compromised, and velocities in the ventricle were reduced. The HH25 LAL ventricles developed a more triangular shape with a sharper apex, and in some cases, the atrioventricular junction shifted medially. These changes led to more sluggish flow near the ventricular free wall and apex, where more fluid particles moved in an oscillatory manner with the motion of the ventricular wall, while slowly being washed out, resulting in lower wall shear stresses and higher oscillatory indices. Consequent to these flow conditions, at HH28, even before septation is complete, the left ventricle was found to be hypoplastic while the right ventricle was found to be larger in compensation. Our results suggest that the low and oscillatory flow near the left side of the heart may play a role in causing the HLHS morphology in the LAL model.