Models of ventricular structure and function reviewed for clinical cardiologists

A meshless fragile points method for rule-based definition of myocardial fiber orientation

BACKGROUND AND OBJECTIVE

Rule-based methods are commonly used to estimate the arrangement of myocardial fibers by solving the Laplace problem with appropriate Dirichlet boundary conditions. Existing algorithms are using the Finite Element Method (FEM) to solve the Laplace-Dirichlet problem. However, meshless methods are under development for cardiac electrophysiology simulation. The objective of this work is to propose a meshless rule based method for the determination of myocardial fiber arrangement without requiring a mesh discretization as it is required by FEM.

METHODS

The proposed method employs the Fragile Points Method (FPM) for the solution of the Laplace-Dirichlet problem. FPM uses simple discontinuous trial functions and single-point exact integration for linear trial functions that set it as a promising alternative to the Finite Element Method. We derive the FPM formulation of the Laplace-Dirichlet and we estimate ventricular and atrial fiber arrangements according to rules based on histology findings for four different geometries. The obtained fiber arrangements from FPM are compared with the ones obtained from FEM by calculating the angle between the fiber vector fields of the two methods for three different directions (i.e., longitudinal, sheet, transverse).

RESULTS

The fiber arrangements that were generated with FPM were in close agreement with the generated arrangements from FEM for all three directions. The mean angle difference between the FPM and FEM vector fields were lower than 0.030^∘ for the ventricular fiber arrangements and lower than 0.036^∘ for the atrial fiber arrangements.

DISCUSSION

The proposed meshless rule-based method was proven to generate myocardial fiber arrangements with very close agreement with FEM while alleviates the requirement for a mesh of the latter. This is of great value for cardiac electrophysiology solvers that are based on meshless methods since they require a well defined myocardial fiber arrangement to simulate accurately the propagation of electrical signals in the heart. Combining a meshless solution for both the determination of the fibers and the electrical signal propagation can allow for solution that do not require the definition of a mesh. To our knowledge, this work is the first one to propose a meshless rule-based method for myocardial fiber arrangement determination.

An Appreciation of Anatomy in the Molecular World

Robert H. Anderson is one of the most important and accomplished cardiac anatomists of the last decades, having made major contributions to our understanding of the anatomy of normal hearts and the pathologies of acquired and congenital heart diseases. While cardiac anatomy as a research discipline has become largely subservient to molecular biology, anatomists like Professor Anderson demonstrate anatomy has much to offer. Here, we provide cases of early anatomical insights on the heart that were rediscovered, and expanded on, by molecular techniques: migration of neural crest cells to the heart was deduced from histological observations (1908) and independently shown again with experimental interventions; pharyngeal mesoderm is added to the embryonic heart (1973) in what is now defined as the molecularly distinguishable second heart field; chambers develop from the heart tube as regional pouches in what is now considered the ballooning model by the molecular identification of regional differentiation and proliferation. The anatomical discovery of the conduction system by Purkinje, His, Tawara, Keith, and Flack is a special case because the main findings were never neglected in later molecular studies. Professor Anderson has successfully demonstrated that sound knowledge of anatomy is indispensable for proper understanding of cardiac development.

A rule-based method to model myocardial fiber orientation in cardiac biventricular geometries with outflow tracts

Rule-based methods are often used for assigning fiber orientation to cardiac anatomical models. However, existing methods have been developed using data mostly from the left ventricle. As a consequence, fiber information obtained from rule-based methods often does not match histological data in other areas of the heart such as the right ventricle, having a negative impact in cardiac simulations beyond the left ventricle. In this work, we present a rule-based method where fiber orientation is separately modeled in each ventricle following observations from histology. This allows to create detailed fiber orientation in specific regions such as the endocardium of the right ventricle, the interventricular septum, and the outflow tracts. We also carried out electrophysiological simulations involving these structures and with different fiber configurations. In particular, we built a modeling pipeline for creating patient-specific volumetric meshes of biventricular geometries, including the outflow tracts, and subsequently simulate the electrical wavefront propagation in outflow tract ventricular arrhythmias with different origins for the ectopic focus. The resulting simulations with the proposed rule-based method showed a very good agreement with clinical parameters such as the 10 ms isochrone ratio in a cohort of nine patients suffering from this type of arrhythmia. The developed modeling pipeline confirms its potential for an in silico identification of the site of origin in outflow tract ventricular arrhythmias before clinical intervention.

The functional architecture of skeletal compared to cardiac musculature: Myocyte orientation, lamellar unit morphology, and the helical ventricular myocardial band

How the cardiomyocytes are aggregated within the heart walls remains contentious. We still do not fully understand how the end-to-end longitudinal myocytic chains are arranged, nor the true extent and shape of the lamellar units they aggregate to form. In this article, we show that an understanding of the complex arrangement of cardiac musculature requires knowledge of three-dimensional myocyte orientation (helical and intrusion angle), and appreciation of myocyte packing within the connective tissue matrix. We show how visualization and segmentation of high-resolution three-dimensional image data can accurately identify the morphology and orientation of the myocytic chains, and the lamellar units. Some maintain that the ventricles can be unwrapped in the form of a "helical ventricular myocardial band," that is, as a compartmentalized band with selective regional innervation and deformation, and a defined origin and insertion like most skeletal muscles. In contrast to the simpler interpretation of the helical ventricular myocardial band, we provide insight as to how the complex myocytic chains, the heterogeneous lamellar units, and connective tissue matrix form an interconnected meshwork, which facilitates the complex internal deformations of the ventricular wall. We highlight the dangers of disregarding the intruding cardiomyocytes. Preparation of the band destroys intruding myocytic chains, and thus disregards the functional implications of the antagonistic auxotonic forces they produce. We conclude that the ventricular myocardium is not analogous to skeletal muscle, but is a complex three-dimensional meshwork, with a heterogeneous branching lamellar architecture.

Linking Genes to Cardiovascular Diseases: Gene Action and Gene-Environment Interactions

A unique myocardial characteristic is its ability to grow/remodel in order to adapt; this is determined partly by genes and partly by the environment and the milieu intérieur. In the "post-genomic" era, a need is emerging to elucidate the physiologic functions of myocardial genes, as well as potential adaptive and maladaptive modulations induced by environmental/epigenetic factors. Genome sequencing and analysis advances have become exponential lately, with escalation of our knowledge concerning sometimes controversial genetic underpinnings of cardiovascular diseases. Current technologies can identify candidate genes variously involved in diverse normal/abnormal morphomechanical phenotypes, and offer insights into multiple genetic factors implicated in complex cardiovascular syndromes. The expression profiles of thousands of genes are regularly ascertained under diverse conditions. Global analyses of gene expression levels are useful for cataloging genes and correlated phenotypes, and for elucidating the role of genes in maladies. Comparative expression of gene networks coupled to complex disorders can contribute insights as to how "modifier genes" influence the expressed phenotypes. Increasingly, a more comprehensive and detailed systematic understanding of genetic abnormalities underlying, for example, various genetic cardiomyopathies is emerging. Implementing genomic findings in cardiology practice may well lead directly to better diagnosing and therapeutics. There is currently evolving a strong appreciation for the value of studying gene anomalies, and doing so in a non-disjointed, cohesive manner. However, it is challenging for many-practitioners and investigators-to comprehend, interpret, and utilize the clinically increasingly accessible and affordable cardiovascular genomics studies. This survey addresses the need for fundamental understanding in this vital area.

Comparison of diffusion tensor imaging by cardiovascular magnetic resonance and gadolinium enhanced 3D image intensity approaches to investigation of structural anisotropy in explanted rat hearts

BACKGROUND

Cardiovascular magnetic resonance (CMR) can through the two methods 3D FLASH and diffusion tensor imaging (DTI) give complementary information on the local orientations of cardiomyocytes and their laminar arrays.

METHODS

Eight explanted rat hearts were perfused with Gd-DTPA contrast agent and fixative and imaged in a 9.4T magnet by two types of acquisition: 3D fast low angle shot (FLASH) imaging, voxels 50 × 50 × 50 μm, and 3D spin echo DTI with monopolar diffusion gradients of 3.6 ms duration at 11.5 ms separation, voxels 200 × 200 × 200 μm. The sensitivity of each approach to imaging parameters was explored.

RESULTS

The FLASH data showed laminar alignments of voxels with high signal, in keeping with the presumed predominance of contrast in the interstices between sheetlets. It was analysed, using structure-tensor (ST) analysis, to determine the most (v1(ST)), intermediate (v2(ST)) and least (v3(ST)) extended orthogonal directions of signal continuity. The DTI data was analysed to determine the most (e1(DTI)), intermediate (e2(DTI)) and least (e3(DTI)) orthogonal eigenvectors of extent of diffusion. The correspondence between the FLASH and DTI methods was measured and appraised. The most extended direction of FLASH signal (v1(ST)) agreed well with that of diffusion (e1(DTI)) throughout the left ventricle (representative discrepancy in the septum of 13.3 ± 6.7°: median ± absolute deviation) and both were in keeping with the expected local orientations of the long-axis of cardiomyocytes. However, the orientation of the least directions of FLASH signal continuity (v3(ST)) and diffusion (e3(ST)) showed greater discrepancies of up to 27.9 ± 17.4°. Both FLASH (v3(ST)) and DTI (e3(DTI)) where compared to directly measured laminar arrays in the FLASH images. For FLASH the discrepancy between the structure-tensor calculated v3(ST) and the directly measured FLASH laminar array normal was of 9 ± 7° for the lateral wall and 7 ± 9° for the septum (median ± inter quartile range), and for DTI the discrepancy between the calculated v3(DTI) and the directly measured FLASH laminar array normal was 22 ± 14° and 61 ± 53.4°. DTI was relatively insensitive to the number of diffusion directions and to time up to 72 hours post fixation, but was moderately affected by b-value (which was scaled by modifying diffusion gradient pulse strength with fixed gradient pulse separation). Optimal DTI parameters were b = 1000 mm/s(2) and 12 diffusion directions. FLASH acquisitions were relatively insensitive to the image processing parameters explored.

CONCLUSIONS

We show that ST analysis of FLASH is a useful and accurate tool in the measurement of cardiac microstructure. While both FLASH and the DTI approaches appear promising for mapping of the alignments of myocytes throughout myocardium, marked discrepancies between the cross myocyte anisotropies deduced from each method call for consideration of their respective limitations.

Advances in induced pluripotent stem cells, genomics, biomarkers, and antiplatelet therapy highlights of the year in JCTR 2013

The Journal provides the clinician and scientist with the latest advances in discovery research, emerging technologies, preclinical research design and testing, and clinical trials. We highlight advances in areas of induced pluripotent stem cells, genomics, biomarkers, multimodality imaging, and antiplatelet biology and therapy. The top publications are critically discussed and presented along with anatomical reviews and FDA insight to provide context.

Mechanotransduction mechanisms for intraventricular diastolic vortex forces and myocardial deformations: part 1

Epigenetic mechanisms are fundamental in cardiac adaptations, remodeling, reverse remodeling, and disease. This two-article series proposes that variable forces associated with diastolic RV/LV rotatory intraventricular flows can exert physiologically and clinically important, albeit still unappreciated, epigenetic actions influencing functional and morphological cardiac adaptations and/or maladaptations. Taken in toto, the two-part survey formulates a new paradigm in which intraventricular diastolic filling vortex-associated forces play a fundamental epigenetic role, and examines how heart cells react to these forces. The objectives are to provide a perspective on vortical epigenetic effects, to introduce emerging ideas, and to suggest directions of multidisciplinary translational research. The main goal is to make pertinent biophysics and cytomechanical dynamic systems concepts accessible to interested translational and clinical cardiologists. I recognize that the diversity of the epigenetic problems can give rise to a diversity of approaches and multifaceted specialized research undertakings. Specificity may dominate the picture. However, I take a contrasting approach. Are there concepts that are central enough that they should be developed in some detail? Broadness competes with specificity. Would, however, this viewpoint allow for a more encompassing view that may otherwise be lost by generation of fragmented results? Part 1 serves as a general introduction, focusing on background concepts, on intracardiac vortex imaging methods, and on diastolic filling vortex-associated forces acting epigenetically on RV/LV endocardium and myocardium. Part 2 will describe pertinent available pluridisciplinary knowledge/research relating to mechanotransduction mechanisms for intraventricular diastolic vortex forces and myocardial deformations and to their epigenetic actions on myocardial and ventricular function and adaptations.

Studying the microanatomy of the heart in three dimensions: a practical update

The structure and function of the heart needs to be understood in three dimensions. We give a brief historical summary of the methods by which such an understanding has been sought, and some practical details of the relatively new technique of micro-CT with iodine contrast enhancement in samples from rat and rabbit. We discuss how the improved anatomical detail available in fixed cadaveric hearts will enhance our ability to model and to understand the integrated function of the cardiomyocytes, conducting tissues, and fibrous supporting structures that generate the pumping function of the heart.