MRI-determined left ventricular "crescent effect": a consequence of the slight deviation of contents of the pericardial sack from the constant-volume state

During one cardiac cycle, the volume encompassed by the pericardial sack in healthy subjects remains nearly constant, with a transient +/-5% decrease in volume at end systole. This "constant-volume" attribute defines a constraint that the longitudinal versus radial pericardial contour dimension relationship must obey. Using cardiac MRI, we determined the extent to which the constant-volume attribute is valid from four-chamber slices (two-dimensional) compared with three-dimensional volumetric data. We also compared the relative percentage of longitudinal versus radial (short-axis) change in cross-sectional area (dimension) of the pericardial contour, thereby assessing the fate of the +/-5% end-systolic volume decrease. We analyzed images from 10 normal volunteers and 1 subject with congenital absence of the pericardium, obtained using a 1.5-T MR scanner. Short-axis cine loop stacks covering the entire heart were acquired, as were single four-chamber cine loops. In the short-axis and four-chamber slices, relative to midventricular end-diastolic location, end-systolic pericardial (left ventricular epicardial) displacement was observed to be radial and maximized at end systole. Longitudinal (apex to mediastinum) pericardial contour dimension change and pericardial area change on the four-chamber slice were negligible throughout the cardiac cycle. We conclude that the +/-5% end-systolic decrease in the volume encompassed by the pericardial sack is primarily accounted for by a "crescent effect" on short-axis views, manifesting as a nonisotropic radial diminution of the pericardial/epicardial contour of the left ventricle. This systolic drop in cardiac volume occurs primarily at the ventricular level and is made up during the subsequent diastole when blood crosses the pericardium in the pulmonary venous Doppler D wave during early rapid left ventricular filling.     

Dynamic geometry of the intact left ventricle

Knowledge of left ventricular chamber dynamics is central to our understanding of cardiac physiology. The complicated changes in left ventricular geometry observed in the dog during various phases of the cardiac cycle can be represented as distinct linear relationships between chamber eccentricity and intracavitary volume during diastole and ejection, and probably represent structural properties of the ventricular wall. Chamber geometry of the left ventricle is a major determinant of overall myocardial function. The slope of the radius of curvature (r) to wall thickness (h) relationship is a geometric constant that determines the mural force at any given transmural pressure. Chronic pressure and volume overload produce changes in this geometric relationship as a result of increased mural force resisting ejection. The adaptive mechanism of ventricular hypertrophy in this setting alters the r/h ratio and returns systolic mural force toward normal. Coronary occlusion induces acute changes in regional geometry characterized by holosystolic wall bulging and systolic wall thinning, which shift the r/h relationship upward and to the left. The geometric alteration during ischemia probably increases systolic mural force and could adversely affect myocardial function. Recent studies with patients have shown the r/h ratio to be of value in distinguishing between reversible and irreversible impairment of myocardial performance. Because most myocardial diseases produce major alterations in the structure of the ventricular wall, analysis of dynamic chamber geometry may prove of prognostic value in assessing patients with cardiac disorders.     

Fluid Dynamics of Ventricular Filling in the Embryonic Heart

The vertebrate embryonic heart first forms as a valveless tube that pumps blood using waves of contraction. As the heart develops, the atrium and ventricle bulge out from the heart tube, and valves begin to form through the expansion of the endocardial cushions. As a result of changes in geometry, conduction velocities, and material properties of the heart wall, the fluid dynamics and resulting spatial patterns of shear stress and transmural pressure change dramatically. Recent work suggests that these transitions are significant because fluid forces acting on the cardiac walls, as well as the activity of myocardial cells that drive the flow, are necessary for correct chamber and valve morphogenesis. In this article, computational fluid dynamics was used to explore how spatial distributions of the normal forces acting on the heart wall change as the endocardial cushions grow and as the cardiac wall increases in stiffness. The immersed boundary method was used to simulate the fluid-moving boundary problem of the cardiac wall driving the motion of the blood in a simplified model of a two-dimensional heart. The normal forces acting on the heart walls increased during the period of one atrial contraction because inertial forces are negligible and the ventricular walls must be stretched during filling. Furthermore, the force required to fill the ventricle increased as the stiffness of the ventricular wall was increased. Increased endocardial cushion height also drastically increased the force necessary to contract the ventricle. Finally, flow in the moving boundary model was compared to flow through immobile rigid chambers, and the forces acting normal to the walls were substantially different.     

Energetics of ventricular contraction as traced in the pressure-volume diagram

The pressure-volume (P-V) relationship of the canine left ventricle can reasonably be simulated by a time-varying elastance model. In this model the total mechanical energy generated by a contraction can be determined theoretically from the change in the elastance. Applying this theory to the actual left ventricle, we have found that the area in the P-V diagram circumscribed by the end-systolic P-V relation line, the end-diastolic P-V relation curve, and the systolic segment of the P-V trajectory is equivalent to the total mechanical energy generated by ventricular contraction. We call this area the systolic P-V area (PVA). We have studied experimentally the correlation between the PVA and myocardial oxygen consumption (VO2) in the canine left ventricle. VO2 was linearly correlated with PVA regardless of the contraction mode and loading conditions in a given left ventricle. The VO2-PVA relation parallel shifted upward with positive inotropic agents. This shift comprised a significant increase in VO2 component for the unloaded contraction. We therefore consider that further analyses of the VO2-PVA relationship will greatly promote our understanding of cardiac energetics.     

A Tissue-Level Electromechanical Model of the Left Ventricle: Application to the Analysis of Intraventricular Pressure

The ventricular pressure profile is characteristic of the cardiac contraction progress and is useful to evaluate the cardiac performance. In this contribution, a tissue-level electromechanical model of the left ventricle is proposed, to assist the interpretation of left ventricular pressure waveforms. The left ventricle has been modeled as an ellipsoid composed of twelve mechano-hydraulic sub-systems. The asynchronous contraction of these twelve myocardial segments has been represented in order to reproduce a realistic pressure profiles. To take into account the different energy domains involved, the tissue-level scale and to facilitate the building of a modular model, multiple formalisms have been used: Bond Graph formalism for the mechano-hydraulic aspects and cellular automata for the electrical activation. An experimental protocol has been defined to acquire ventricular pressure signals from three pigs, with different afterload conditions. Evolutionary Algorithms have been used to identify the model parameters in order to minimize the error between experimental and simulated ventricular pressure signals. Simulation results show that the model is able to reproduce experimental ventricular pressure. In addition, electro-mechanical activation times have been determined in the identification process. For example, the maximum electrical activation time is reached, respectively, 96.5, 139.3 and 131.5 ms for the first, second, and third pigs. These preliminary results are encouraging for the application of the model on non-invasive data like ECG, arterial pressure or myocardial strain.     

Mutation of weak atrium/atrial myosin heavy chain disrupts atrial function and influences ventricular morphogenesis in zebrafish

The embryonic vertebrate heart is composed of two major chambers, a ventricle and an atrium, each of which has a characteristic size, shape and functional capacity that contributes to efficient circulation. Chamber-specific gene expression programs are likely to regulate key aspects of chamber formation. Here, we demonstrate that epigenetic factors also have a significant influence on chamber morphogenesis. Specifically, we show that an atrium-specific contractility defect has a profound impact on ventricular development. We find that the zebrafish locus weak atrium encodes an atrium-specific myosin heavy chain that is required for atrial myofibrillar organization and contraction. Despite their atrial defects, weak atrium mutants can maintain circulation through ventricular contraction. However, the weak atrium mutant ventricle becomes unusually compact, exhibiting a thickened myocardial wall, a narrow lumen and changes in myocardial gene expression. As weak atrium/atrial myosin heavy chain is expressed only in the atrium, the ventricular phenotypes in weak atrium mutants represent a secondary response to atrial dysfunction. Thus, not only is cardiac form essential for cardiac function, but there also exists a reciprocal relationship in which function can influence form. These findings are relevant to our understanding of congenital defects in cardiac chamber morphogenesis.     

Possibility Study of Scale Invariant Feature Transform (SIFT) Algorithm Application to Spine Magnetic Resonance Imaging

The purpose of this study is an application of scale invariant feature transform (SIFT) algorithm to stitch the cervical-thoracic-lumbar (C-T-L) spine magnetic resonance (MR) images to provide a view of the entire spine in a single image. All MR images were acquired with fast spin echo (FSE) pulse sequence using two MR scanners (1.5 T and 3.0 T). The stitching procedures for each part of spine MR image were performed and implemented on a graphic user interface (GUI) configuration. Moreover, the stitching process is performed in two categories; manual point-to-point (mPTP) selection that performed by user specified corresponding matching points, and automated point-to-point (aPTP) selection that performed by SIFT algorithm. The stitched images using SIFT algorithm showed fine registered results and quantitatively acquired values also indicated little errors compared with commercially mounted stitching algorithm in MRI systems. Our study presented a preliminary validation of the SIFT algorithm application to MRI spine images, and the results indicated that the proposed approach can be performed well for the improvement of diagnosis. We believe that our approach can be helpful for the clinical application and extension of other medical imaging modalities for image stitching.     

Finite-element analysis of left ventricular myocardial stresses

A versatile method of finite-element analysis is presented for the determination of the stress distributions in the left ventricular myocardial wall. The instantaneous shapes of the left ventricular myocardial wall, measured at 0,5 mm intervals and at a rate 0f 60 images/sec during a cardiac cycle, are approximated by axisymmetric shells following the approach of Gould et al. and analysed by the method of incremental loadings to account for the changing transmural pressure. The ventricular wall is mathematically divided up into coaxial rings of triagular cross sections so that determination of the stresses at any point within the wall can be achieved by assigning increased number of nodes across the wall thickness in the regions of the left ventricular wall where particular attention is needed. Appropriate boundary conditions are defined at the base of the left ventricle so that it can be treated as a shell with an open end. The computer program, which implements all the stress calculations involved, depends on the dimensions of the left ventricular wall measured from an operator-interactive roengen videometry system. It carries out the sequential formation of the nodes and elements and includes a CALCOMP subroutine to plot the finite-element partitioning of the instantaneous shape. Illustrative results of the end-diastolic stress distributions within the myocardial wall of a metabolically-supported, isolated, working canine left ventricle are given. This technique predicts higher endocardial meridional and hoop wall stresses relative to the stresses in the middle and epicardial region than those obtained with previous models.     

Patient-specific Modeling of the Heart: Estimation of Ventricular Fiber Orientations

Patient-specific simulations of heart (dys)function aimed at personalizing cardiac therapy are hampered by the absence of in vivo imaging technology for clinically acquiring myocardial fiber orientations. The objective of this project was to develop a methodology to estimate cardiac fiber orientations from in vivo images of patient heart geometries. An accurate representation of ventricular geometry and fiber orientations was reconstructed, respectively, from high-resolution ex vivo structural magnetic resonance (MR) and diffusion tensor (DT) MR images of a normal human heart, referred to as the atlas. Ventricular geometry of a patient heart was extracted, via semiautomatic segmentation, from an in vivo computed tomography (CT) image. Using image transformation algorithms, the atlas ventricular geometry was deformed to match that of the patient. Finally, the deformation field was applied to the atlas fiber orientations to obtain an estimate of patient fiber orientations. The accuracy of the fiber estimates was assessed using six normal and three failing canine hearts. The mean absolute difference between inclination angles of acquired and estimated fiber orientations was 15.4 °. Computational simulations of ventricular activation maps and pseudo-ECGs in sinus rhythm and ventricular tachycardia indicated that there are no significant differences between estimated and acquired fiber orientations at a clinically observable level.The new insights obtained from the project will pave the way for the development of patient-specific models of the heart that can aid physicians in personalized diagnosis and decisions regarding electrophysiological interventions.     

Electrical Wave Propagation in an Anisotropic Model of the Left Ventricle Based on Analytical Description of Cardiac Architecture

We develop a numerical approach based on our recent analytical model of fiber structure in the left ventricle of the human heart. A special curvilinear coordinate system is proposed to analytically include realistic ventricular shape and myofiber directions. With this anatomical model, electrophysiological simulations can be performed on a rectangular coordinate grid. We apply our method to study the effect of fiber rotation and electrical anisotropy of cardiac tissue (i.e., the ratio of the conductivity coefficients along and across the myocardial fibers) on wave propagation using the ten Tusscher–Panfilov (2006) ionic model for human ventricular cells. We show that fiber rotation increases the speed of cardiac activation and attenuates the effects of anisotropy. Our results show that the fiber rotation in the heart is an important factor underlying cardiac excitation. We also study scroll wave dynamics in our model and show the drift of a scroll wave filament whose velocity depends non-monotonically on the fiber rotation angle; the period of scroll wave rotation decreases with an increase of the fiber rotation angle; an increase in anisotropy may cause the breakup of a scroll wave, similar to the mother rotor mechanism of ventricular fibrillation.     

A framework for biomechanics simulations using four-chamber cardiac models

Computational cardiac models have been extensively used to study different cardiac biomechanics; specifically, finite-element analysis has been one of the tools used to study the internal stresses and strains in the cardiac wall during the cardiac cycle. Cubic-Hermite finite element meshes have been used for simulating cardiac biomechanics due to their convergence characteristics and their ability to capture smooth geometries compactly–fewer elements are needed to build the cardiac geometry–compared to linear tetrahedral meshes. Such meshes have previously been used only with simple ventricular geometries with non-physiological boundary conditions due to challenges associated with creating cubic-Hermite meshes of the complex heart geometry. However, it is critical to accurately capture the different geometric characteristics of the heart and apply physiologically equivalent boundary conditions to replicate the in vivo heart motion. In this work, we created a four-chamber cardiac model utilizing cubic-Hermite elements and simulated a full cardiac cycle by coupling the 3D finite element model with a lumped circulation model. The myocardial fiber-orientations were interpolated within the mesh using the Log-Euclidean method to overcome the singularity associated with interpolation of orthogonal matrices. Physiologically equivalent rigid body constraints were applied to the nodes along the valve plane and the accuracy of the resulting simulations were validated using open source clinical data. We then simulated a complete cardiac cycle of a healthy heart and a heart with acute myocardial infarction. We compared the pumping functionality of the heart for both cases by calculating the ventricular work. We observed a 20% reduction in acute work done by the heart immediately after myocardial infarction. The myocardial wall displacements obtained from the four-chamber model are comparable to actual patient data, without requiring complicated non-physiological boundary conditions usually required in truncated ventricular heart models.     

Use of Echocardiography Reveals Reestablishment of Ventricular Pumping Efficiency and Partial Ventricular Wall Motion Recovery upon Ventricular Cryoinjury in the Zebrafish

Aims

While zebrafish embryos are amenable to in vivo imaging, allowing the study of morphogenetic processes during development, intravital imaging of adults is hampered by their small size and loss of transparency. The use of adult zebrafish as a vertebrate model of cardiac disease and regeneration is increasing at high speed. It is therefore of great importance to establish appropriate and robust methods to measure cardiac function parameters.

Methods and Results

Here we describe the use of 2D-echocardiography to study the fractional volume shortening and segmental wall motion of the ventricle. Our data show that 2D-echocardiography can be used to evaluate cardiac injury and also to study recovery of cardiac function. Interestingly, our results show that while global systolic function recovered following cardiac cryoinjury, ventricular wall motion was only partially restored.

Conclusion

Cryoinjury leads to long-lasting impairment of cardiac contraction, partially mimicking the consequences of myocardial infarction in humans. Functional assessment of heart regeneration by echocardiography allows a deeper understanding of the mechanisms of cardiac regeneration and has the advantage of being easily transferable to other cardiovascular zebrafish disease models.     

MR tagging demonstrates quantitative differences in regional ventricular wall motion in mice, rats, and men

Rats and genetically manipulated mouse models have played an important role in the exploration of molecular causes of cardiovascular diseases. However, it has not been fully investigated whether mice or rats and humans manifest similar patterns of ventricular wall motion. Although similarities in anatomy and myofiber architecture suggest that fundamental patterns of ventricular wall motion may be similar, the considerable differences in heart size, heart rate, and sarcomeric protein isoforms may yield quantitative differences in ventricular wall mechanics. To further our understanding of the basic mechanisms of myofiber contractile performance, we quantified regional and global indexes of ventricular wall motion in mice, rats, and men using magnetic resonance (MR) imaging. Both regular cine and tagged MR images at apical, midventricular, and basal levels were acquired from six male volunteers, six Fischer 344 rats, and seven C57BL/6 mice. Morphological parameters and ejection fraction were computed directly from cine images. Myocardial twist (rotation angle), torsion (net twist per unit length), circumferential strain, and normalized radial shortening were calculated by homogeneous strain analysis from tagged images. Our data show that ventricular twist was conserved among the three species, leading to a significantly smaller torsion, measured as net twist per unit length, in men. However, both circumferential strain and normalized radial shortening were the largest in male subjects. Although other parameters, such as circumferential-longitudinal shear strain, need to be evaluated, and the causes of these differences in contractile mechanics remain to be elucidated, the preservation of twist appears fundamental to cardiac function and should be considered in studies that extrapolate data from animals to humans.     

A modified mathematical model of the anatomy of the cardiac left ventricle

A modification of the mathematical model of the shape and fiber direction field of the left cardiac ventricle is presented. The model was developed based on the idea of nested spiral surfaces. The ventricle is composed of surfaces that model myocardial layers. Each layer is filled with curves corresponding to myocardial fibers. The tangents to these curves form the myofiber direction field. A modified spherical coordinate system is linked with the model left ventricle, where the ventricular boundaries are coordinate surfaces. The model is based on echocardiographic, computed-tomography, or magnetic-resonance-imaging data. For this purpose, four-chamber and two-chamber echocardiography views or sections along the long axis of the left ventricle from these tomographic data in several positions are approximated with a model profile. To construct a 3D model, we then interpolate model parameters by periodic cubic splines and the vector field of the tangents to the model fibers is calculated. For verification of the model, we used diffusion-tensor magneticresonance-imaging data of the human heart.     

Automatic Localization of the Left Ventricle from Cardiac Cine Magnetic Resonance Imaging: A New Spectrum-Based Computer-Aided Tool

Traditionally, cardiac image analysis is done manually. Automatic image processing can help with the repetitive tasks, and also deal with huge amounts of data, a task which would be humanly tedious. This study aims to develop a spectrum-based computer-aided tool to locate the left ventricle using images obtained via cardiac magnetic resonance imaging. Discrete Fourier Transform was conducted pixelwise on the image sequence. Harmonic images of all frequencies were analyzed visually and quantitatively to determine different patterns of the left and right ventricles on spectrum. The first and fifth harmonic images were selected to perform an anisotropic weighted circle Hough detection. This tool was then tested in ten volunteers. Our tool was able to locate the left ventricle in all cases and had a significantly higher cropping ratio of 0.165 than did earlier studies. In conclusion, a new spectrum-based computer aided tool has been proposed and developed for automatic left ventricle localization. The development of this technique, which will enable the automatic location and further segmentation of the left ventricle, will have a significant impact in research and in diagnostic settings. We envisage that this automated method could be used by radiographers and cardiologists to diagnose and assess ventricular function in patients with diverse heart diseases.     

Alteration of collagenous protein profile in congestive heart failure secondary to myocardial infarction

The rat model of myocardial infarction is characterized by progressive cardiac hypertrophy and failure. Rats with infarcts greater than 30% of the left ventricle exhibited early and moderate, stages of heart failure 4 and 8 weeks after the occlusion of the left coronary artery, respectively. As heart failure is usually associated with remodeling of the extracellular matrix, a histological and biochemical study of cardiac collagenous proteins was carried out using failing hearts. Total collagen content in the right ventricle increased at 2, 4, and 8 weeks following occlusion of the left coronary artery whereas such a change in viable left ventricle was seen after 4 and 8 weeks. Total cardiac hydroxyproline concentration was increased in both right and left ventricular samples from the infarcted animals when compared to those of control; this increase was due to elevation of pepsin-insoluble collagen fraction. The myocardial noncollagenous/collagenous protein ratio was decreased in experimental right and left ventricular samples when compared to control samples. These findings suggest that an increase in cross-linking of cardiac collagen as well as disparate synthesis of collagenous and noncollagenous proteins occurs in this model of congestive heart, failure.     

Theoretical models in mechanics of the left ventricle

Different rheological concepts and theoretical studies have been recently presented using models of myocardial mechanics. Complex analysis of the mechanical behavior of the left ventricular wall have been developed in order to estimate the local stresses and deformations that occur during the heart cycle as well as the ventricular stroke volume and pressure. Theoretical models have taken into account non-linear and viscoelastic passive properties of the myocardium tissue, when subjected to large deformations, through given strain energy functions or stress-strain relations. Different prolate spheroid geometries have been considered for such thick shell cardiac structure. During the active state of the contraction, the rheological behavior of the fibers has been described using different muscle models and relationships between fiber tension and strain, and activation degree. A forthcoming approach for bridging the gap between the knowledge of the muscle fiber microrheological properties and the study of the mechanical behavior of the entire ventricle, consists in including anisotropic and inhomogeneous effects through fiber direction field.     

Structural finite deformation model of the left ventricle during diastole and systole

A model of left ventricular function is developed based on morphological characteristics of the myocardial tissue. The passive response of the three-dimensional collagen network and the active contribution of the muscle fibers are integrated to yield the overall response of the left ventricle which is considered to be a thick wall cylinder. The deformation field and the distributions of stress and pressure are determined at each point in the cardiac cycle by numerically solving three equations of equilibrium. Simulated results in terms of the ventricular deformation during ejection and isovolumic cycles are shown to be in good qualitative agreement with experimental data. It is shown that the collagen network in the heart has considerable effect on the pressure-volume loops. The particular pattern of spatial orientation of the collagen determines the ventricular recoil properties in early diastole. The material properties (myocardial stiffness and contractility) are shown to affect both the pressure-volume loop and the deformation pattern of the ventricle. The results indicate that microstructural consideration offer a realistic representation of the left ventricle mechanics.     

Extensive scar formation and regression during heart regeneration after cryoinjury in zebrafish

The zebrafish heart has the capacity to regenerate after ventricular resection. Although this regeneration model has proved useful for the elucidation of certain regeneration mechanisms, it is based on the removal of heart tissue rather than its damage. Here, we characterize the cellular response and regenerative capacity of the zebrafish heart after cryoinjury, an alternative procedure that more closely models the pathophysiological process undergone by the human heart after myocardial infarction (MI). Localized damage was induced in 25% of the ventricle by cryocauterization (CC). During the first 24 hours post-injury, CC leads to cardiomyocyte death within the injured area and the near coronary vasculature. Cell death is followed by a rapid proliferative response in endocardium, epicardium and myocardium. During the first 3 weeks post-injury cell debris was cleared and the injured area replaced by a massive scar. The fibrotic tissue was subsequently degraded and replaced by cardiac tissue. Although animals survived CC, their hearts showed nonhomogeneous ventricular contraction and had a thickened ventricular wall, suggesting that regeneration is associated with processes resembling mammalian ventricular remodeling after acute MI. Our results provide the first evidence that, like mammalian hearts, teleost hearts undergo massive fibrosis after cardiac damage. Unlike mammals, however, the fish heart can progressively eliminate the scar and regenerate the lost myocardium, indicating that scar formation is compatible with myocardial regeneration and the existence of endogenous mechanisms of scar regression. This finding suggests that CC-induced damage in zebrafish could provide a valuable model for the study of the mechanisms of scar removal post-MI.