Three-dimensional alignment of the aggregated myocytes in the normal and hypertrophic murine heart

Diffusion Tensor Phenomapping of the Healthy and Pressure-Overloaded Human Heart

Current techniques to image the microstructure of the heart with diffusion tensor MRI (DTI) are highly under-resolved. We present a technique to improve the spatial resolution of cardiac DTI by almost 10-fold and leverage this to measure local gradients in cardiomyocyte alignment or helix angle (HA). We further introduce a phenomapping approach based on voxel-wise hierarchical clustering of these gradients to identify distinct microstructural microenvironments in the heart. Initial development was performed in healthy volunteers (n=8). Thereader, subjects with severe but well-compensated aortic stenosis (AS, n=10) were compared to age-matched controls (CTL, n=10). Radial HA gradient was significantly reduced in AS (8.0±0.8°/mm vs. 10.2±1.8°/mm, p=0.001) but the other HA gradients did not change significantly. Four distinct microstructural clusters could be idenJfied in both the CTL and AS subjects and did not differ significantly in their properties or distribution. Despite marked hypertrophy, our data suggest that the myocardium in well-compensated AS can maintain its microstructural coherence. The described phenomapping approach can be used to characterize microstructural plasticity and perturbation in any organ system and disease.

Cardiomyocyte orientation recovery at micrometer scale reveals long-axis fiber continuum in heart walls

Coordinated cardiomyocyte contraction drives the mammalian heart to beat and circulate blood. No consensus model of cardiomyocyte geometrical arrangement exists, due to the limited spatial resolution of whole heart imaging methods and the piecemeal nature of studies based on histological sections. By combining microscopy and computer vision, we produced the first-ever three-dimensional cardiomyocyte orientation reconstruction across mouse ventricular walls at the micrometer scale, representing a gain of three orders of magnitude in spatial resolution. We recovered a cardiomyocyte arrangement aligned to the long-axis direction of the outer ventricular walls. This cellular network lies in a thin shell and forms a continuum with longitudinally arranged cardiomyocytes in the inner walls, with a complex geometry at the apex. Our reconstruction methods can be applied at fine spatial scales to further understanding of heart wall electrical function and mechanics, and set the stage for the study of micron-scale fiber remodeling in heart disease.

Multiscale characterization of left ventricle active behavior in the mouse

The myocardium possesses an intricately designed microarchitecture to produce an optimal cardiac contraction. The contractile behavior of the heart is generated at the sarcomere level and travels across several length scales to manifest as the systolic function at the organ level. While passive myocardial behavior has been studied extensively, the translation of active tension produced at the fiber level to the organ-level function is not well understood. Alterations in cardiac systolic function are often key sequelae in structural heart diseases, such as myocardial infarction and systolic heart failure; thus, characterization of the contractile behavior of the heart across multiple length scales is essential to improve our understanding of mechanisms collectively leading to depressed systolic function. In this study, we present a methodology to characterize the active behavior of left ventricle free wall (LVFW) myocardial tissues in mice. Combined with active tests in papillary muscle fibers and conventional in vivo contractility measurement at the organ level in an animal-specific manner, we establish a multiscale active characterization of the heart from fiber to organ. In addition, we quantified myocardial architecture from histology to shed light on the directionality of the contractility at the tissue level. The LVFW tissue activation-relaxation behavior under isometric conditions was qualitatively similar to that of the papillary muscle fiber bundle. However, the maximum stress developed in the LVFW tissue was an order of magnitude lower than that developed by a fiber bundle, and the time taken for active forces to plateau was 2-3 orders of magnitude longer. Although the LVFW tissue exhibited a slightly stiffer passive response in the circumferential direction, the tissues produced significantly larger active stresses in the longitudinal direction during active testing. Also, contrary to passive viscoelastic stress relaxation, active stresses relaxed faster in the direction with larger peak stresses. The multiscale experimental pipeline presented in this work is expected to provide crucial insight into the contractile adaptation mechanisms of the heart with impaired systolic function. STATEMENT OF SIGNIFICANCE: Heart failure cause significant alterations to the contractile-relaxation behavior of the yocardium. Multiscale characterization of the contractile behavior of the myocardium is essential to advance our understanding of how contractility translates from fiber to organ and to identify the multiscale mechanisms leading to impaired cardiac function. While passive myocardial behavior has been studied extensively, the investigation of tissue-level contractile behavior remains critically scarce in the literature. To the best of our knowledge, our study here is the first to investigate the contractile behavior of the left ventricle at multiple length scales in small animals. Our results indicate that the active myocardial wall is a function of transmural depth and relaxes faster in the direction with larger peak stresses.

Detailed quantification of cardiac ventricular myocardial architecture in the embryonic and fetal mouse heart by application of structure tensor analysis to high resolution episcopic microscopic data

The mammalian heart, which is one of the first organs to form and function during embryogenesis, develops from a simple tube into a complex organ able to efficiently pump blood towards the rest of the body. The progressive growth of the compact myocardium during embryonic development is accompanied by changes in its structural complexity and organisation. However, how myocardial myoarchitecture develops during embryogenesis remain poorly understood. To date, analysis of heart development has focused mainly on qualitative descriptions using selected 2D histological sections. High resolution episcopic microscopy (HREM) is a novel microscopic imaging technique that enables to obtain high-resolution three-dimensional images of the heart and perform detailed quantitative analyses of heart development. In this work, we performed a detailed characterization of the development of myocardial architecture in wildtype mice, from E14.5 to E18.5, by means of structure tensor analysis applied to HREM images of the heart. Our results shows that even at E14.5, myocytes are already aligned, showing a gradual change in their helical angle from positive angulation in the endocardium towards negative angulation in the epicardium. Moreover, there is gradual increase in the degree of myocardial organisation concomitant with myocardial growth. However, the development of the myoarchitecture is heterogeneous showing regional differences between ventricles, ventricular walls as well as between myocardial layers, with different growth patterning between the endocardium and epicardium. We also found that the percentage of circumferentially arranged myocytes within the LV significantly increases with gestational age. Finally, we found that fractional anisotropy (FA) within the LV gradually increases with gestational age, while the FA within RV remains unchanged.

Development of Helical Myofiber Tracts in the Human Fetal Heart: Analysis of Myocardial Fiber Formation in the Left Ventricle From the Late Human Embryonic Period Using Diffusion Tensor Magnetic Resonance Imaging

Background

Detection of the fiber orientation pattern of the myocardium using diffusion tensor magnetic resonance imaging lags ≈12 weeks of gestational age (WGA) behind fetal myocardial remodeling with invasion by the developing coronary vasculature (8 WGA). We aimed to use diffusion tensor magnetic resonance imaging tractography to characterize the evolution of fiber architecture in the developing human heart from the later embryonic period.

Methods and Results

Twenty human specimens (8–24 WGA) from the Kyoto Collection of Human Embryos and Fetuses, including specimens from the embryonic period (Carnegie stages 20–23), were used. Diffusion tensor magnetic resonance imaging data were acquired with a 7T magnetic resonance system. Fractional anisotropy and helix angle were calculated using standard definitions. In all samples, the fibers ran helically in an organized pattern in both the left and right ventricles. A smooth transmural change in helix angle values (from positive to negative) was detected in all 16 directions of the ventricles. This feature was observed in almost all small (Carnegie stage 23) and large samples. A higher fractional anisotropy value was detected at the outer side of the anterior wall and septum at Carnegie stage 20 to 22, which spread around the ventricular wall at Carnegie stage 23 and in the early fetal samples (11–12 WGA). The fractional anisotropy value of the left ventricular walls decreased in samples with ≥13 WGA, which remained low (≈0.09) in larger samples.

Conclusions

From the human late embryonic period (from 8 WGA), the helix angle arrangement of the myocardium is comparable to that of the adult, indicating that the myocardial structure blueprint, organization, and integrity are already formed.

Myoarchitectural disarray of hypertrophic cardiomyopathy begins pre-birth

Myoarchitectural disarray – the multiscalar disorganisation of myocytes, is a recognised histopathological hallmark of adult human hypertrophic cardiomyopathy (HCM). It occurs before the establishment of left ventricular hypertrophy (LVH) but its early origins and evolution around the time of birth are unknown. Our aim is to investigate whether myoarchitectural abnormalities in HCM are present in the fetal heart. We used wild‐type, heterozygous and homozygous hearts (n = 56) from a Mybpc3‐targeted knock‐out HCM mouse model and imaged the 3D micro‐structure by high‐resolution episcopic microscopy. We developed a novel structure tensor approach to extract, display and quantify myocyte orientation and its local angular uniformity by helical angle, angle of intrusion and myoarchitectural disarray index, respectively, immediately before and after birth. In wild‐type, we demonstrate uniformity of orientation of cardiomyocytes with smooth transitions of helical angle transmurally both before and after birth but with traces of disarray at the septal insertion points of the right ventricle. In comparison, heterozygous mice free of LVH, and homozygous mice showed not only loss of the normal linear helical angulation transmural profiles observed in wild‐type but also fewer circumferentially arranged myocytes at birth. Heterozygous and homozygous showed more disarray with a wider distribution than in wild‐type before birth. In heterozygous mice, disarray was seen in the anterior, septal and inferior walls irrespective of stage, whereas in homozygous mice it extended to the whole LV circumference including the lateral wall. In conclusion, myoarchitectural disarray is detectable in the fetal heart of an HCM mouse model before the development of LVH.

High resolution 3-Dimensional imaging of the human cardiac conduction system from microanatomy to mathematical modeling

Cardiac arrhythmias and conduction disturbances are accompanied by structural remodelling of the specialised cardiomyocytes known collectively as the cardiac conduction system. Here, using contrast enhanced micro-computed tomography, we present, in attitudinally appropriate fashion, the first 3-dimensional representations of the cardiac conduction system within the intact human heart. We show that cardiomyocyte orientation can be extracted from these datasets at spatial resolutions approaching the single cell. These data show that commonly accepted anatomical representations are oversimplified. We have incorporated the high-resolution anatomical data into mathematical simulations of cardiac electrical depolarisation. The data presented should have multidisciplinary impact. Since the rate of depolarisation is dictated by cardiac microstructure, and the precise orientation of the cardiomyocytes, our data should improve the fidelity of mathematical models. By showing the precise 3-dimensional relationships between the cardiac conduction system and surrounding structures, we provide new insights relevant to valvar replacement surgery and ablation therapies. We also offer a practical method for investigation of remodelling in disease, and thus, virtual pathology and archiving. Such data presented as 3D images or 3D printed models, will inform discussions between medical teams and their patients, and aid the education of medical and surgical trainees.

Transmural gradients of myocardial structure and mechanics: Implications for fiber stress and strain in pressure overload

Although a truly complete understanding of whole heart activation, contraction, and deformation is well beyond our current reach, a significant amount of effort has been devoted to discovering and understanding the mechanisms by which myocardial structure determines cardiac function to better treat patients with cardiac disease. Several experimental studies have shown that transmural fiber strain is relatively uniform in both diastole and systole, in contrast to predictions from traditional mechanical theory. Similarly, mathematical models have largely predicted uniform fiber stress across the wall. The development of this uniform pattern of fiber stress and strain during filling and ejection is due to heterogeneous transmural distributions of several myocardial structures. This review summarizes these transmural gradients, their contributions to fiber mechanics, and the potential functional effects of their remodeling during pressure overload hypertrophy.

Parametric Modeling of the Mouse Left Ventricular Myocardial Fiber Structure

Magnetic resonance diffusion tensor imaging (DTI) has greatly facilitated detailed quantifications of myocardial structures. However, structural patterns, such as the distinctive transmural rotation of the fibers, remain incompletely described. To investigate the validity and practicality of pattern-based analysis, 3D DTI was performed on 13 fixed mouse hearts and fiber angles in the left ventricle were transformed and fitted to parametric expressions constructed from elementary functions of the prolate spheroidal spatial variables. It was found that, on average, the myocardial fiber helix angle could be represented to 6.5° accuracy by the equivalence of a product of 10th-order polynomials of the radial and longitudinal variables, and 17th-order Fourier series of the circumferential variable. Similarly, the fiber imbrication angle could be described by 10th-order polynomials and 24th-order Fourier series, to 5.6° accuracy. The representations, while relatively concise, did not adversely affect the information commonly derived from DTI datasets including the whole-ventricle mean fiber helix angle transmural span and atlases constructed for the group. The unique ability of parametric models for predicting the 3D myocardial fiber structure from finite number of 2D slices was also demonstrated. These findings strongly support the principle of parametric modeling for characterizing myocardial structures in the mouse and beyond.

Finite-Element Extrapolation of Myocardial Structure Alterations Across the Cardiac Cycle in Rats

Myocardial microstructures are responsible for key aspects of cardiac mechanical function. Natural myocardial deformation across the cardiac cycle induces measurable structural alteration, which varies across disease states. Diffusion tensor magnetic resonance imaging (DT-MRI) has become the tool of choice for myocardial structural analysis. Yet, obtaining the comprehensive structural information of the whole organ, in 3D and time, for subject-specific examination is fundamentally limited by scan time. Therefore, subject-specific finite-element (FE) analysis of a group of rat hearts was implemented for extrapolating a set of initial DT-MRI to the rest of the cardiac cycle. The effect of material symmetry (isotropy, transverse isotropy, and orthotropy), structural input, and warping approach was observed by comparing simulated predictions against in vivo MRI displacement measurements and DT-MRI of an isolated heart preparation at relaxed, inflated, and contracture states. Overall, the results indicate that, while ventricular volume and circumferential strain are largely independent of the simulation strategy, structural alteration predictions are generally improved with the sophistication of the material model, which also enhances torsion and radial strain predictions. Moreover, whereas subject-specific transversely isotropic models produced the most accurate descriptions of fiber structural alterations, the orthotropic models best captured changes in sheet structure. These findings underscore the need for subject-specific input data, including structure, to extrapolate DT-MRI measurements across the cardiac cycle.

Small animal cardiovascular MR imaging and spectroscopy

The use of MR imaging and spectroscopy for studying cardiovascular disease processes in small animals has increased tremendously over the past decade. This is the result of the remarkable advances in MR technologies and the increased availability of genetically modified mice. MR techniques provide a window on the entire timeline of cardiovascular disease development, ranging from subtle early changes in myocardial metabolism that often mark disease onset to severe myocardial dysfunction associated with end-stage heart failure. MR imaging and spectroscopy techniques play an important role in basic cardiovascular research and in cardiovascular disease diagnosis and therapy follow-up. This is due to the broad range of functional, structural and metabolic parameters that can be quantified by MR under in vivo conditions non-invasively. This review describes the spectrum of MR techniques that are employed in small animal cardiovascular disease research and how the technological challenges resulting from the small dimensions of heart and blood vessels as well as high heart and respiratory rates, particularly in mice, are tackled.

Effects of incremental beta-blocker dosing on myocardial mechanics of the human left ventricle: MRI 3D-tagging insight into pharmacodynamics supports theory of inner antagonism

Beta-blockers contribute to treatment of heart failure. Their mechanism of action, however, is incompletely understood. Gradients in beta-blocker sensitivity of helically aligned cardiomyocytes compared with counteracting transversely intruding cardiomyocytes seem crucial. We hypothesize that selective blockade of transversely intruding cardiomyocytes by low-dose beta-blockade unloads ventricular performance. Cardiac magnetic resonance imaging (MRI) 3D tagging delivers parameters of myocardial performance. We studied 13 healthy volunteers by MRI 3D tagging during escalated intravenous administration of esmolol. The circumferential, longitudinal, and radial myocardial shortening was determined for each dose. The curves were analyzed for peak value, time-to-peak, upslope, and area-under-the-curve. At low doses, from 5 to 25 μg·kg−1·min−1, peak contraction increased while time-to-peak decreased yielding a steeper upslope. Combining the values revealed a left shift of the curves at low doses compared with baseline without esmolol. At doses of 50 to 150 μg·kg−1·min−1, a right shift with flattening occurred. In healthy volunteers we found more pronounced myocardial shortening at low compared with clinical dosage of beta-blockers. In patients with ventricular hypertrophy and higher prevalence of transversely intruding cardiomyocytes selective low-dose beta-blockade could be even more effective. MRI 3D tagging could help to determine optimal individual beta-blocker dosing avoiding undesirable side effects.

A high-resolution cardiovascular magnetic resonance diffusion tensor map from ex-vivo C57BL/6 murine hearts

BACKGROUND

The complex cardiac fiber structural organization and spatial arrangement of cardiomyocytes in laminar sheetlets contributes greatly to cardiac functional and contractile ejection patterns. This study presents the first comprehensive, ultra-high resolution, fully quantitative statistical tensor map of the fixed murine heart at isotropic resolution of 43 μm using diffusion tensor (DT) cardiovascular magnetic resonance (CMR).

METHODS

Imaging was completed in approximately 12 hours using a six-directional encoding scheme, in five ex vivo healthy C57BL/6 mouse hearts. The tensor map constructed from this data provides an average description of the murine fiber architecture visualized with fiber tractography, and its population variability, using the latest advances in image tensor analysis and statistics.

RESULTS

Results show that non-normalized cardiac tensor maps are associated with mean fractional anisotropy of 0.25 ± 0.07 and mean diffusivity of 8.9 ± 1.6 × 10⁻⁴mm²/s. Moreover, average mid-ventricular helical angle distributions ranged between -41 ± 3° and +52 ± 5° and were highly correlated with transmural depth, in agreement with prior published results in humans and canines. Calculated variabilities of local myocyte orientations were 2.0° and 1.4°. Laminar sheet orientation variability was found to be less stable at 2.6°. Despite such variations, the murine heart seems to be highly structured, particularly when compared to canines and humans.

CONCLUSIONS

This tensor map has the potential to yield an accurate mean representation and identification of common or unique features of the cardiac myocyte architecture, to establish a baseline standard reference of DTI indices, and to improve detection of biomarkers, especially in pathological states or post-transgenetic modifications.

Generating fibre orientation maps in human heart models using Poisson interpolation

Smoothly varying muscle fibre orientations in the heart are critical to its electrical and mechanical function. From detailed histological studies and diffusion tensor imaging, we now know that fibre orientations in humans vary gradually from approximately -70° in the outer wall to +80° in the inner wall. However, the creation of fibre orientation maps for computational analyses remains one of the most challenging problems in cardiac electrophysiology and cardiac mechanics. Here, we show that Poisson interpolation generates smoothly varying vector fields that satisfy a set of user-defined constraints in arbitrary domains. Specifically, we enforce the Poisson interpolation in the weak sense using a standard linear finite element algorithm for scalar-valued second-order boundary value problems and introduce the feature to be interpolated as a global unknown. User-defined constraints are then simply enforced in the strong sense as Dirichlet boundary conditions. We demonstrate that the proposed concept is capable of generating smoothly varying fibre orientations, quickly, efficiently and robustly, both in a generic bi-ventricular model and in a patient-specific human heart. Sensitivity analyses demonstrate that the underlying algorithm is conceptually able to handle both arbitrarily and uniformly distributed user-defined constraints; however, the quality of the interpolation is best for uniformly distributed constraints. We anticipate our algorithm to be immediately transformative to experimental and clinical settings, in which it will allow us to quickly and reliably create smooth interpolations of arbitrary fields from data-sets, which are sparse but uniformly distributed.

Progressive changes in T₁, T₂ and left-ventricular histo-architecture in the fixed and embedded rat heart

Chemical tissue fixation, followed by embedding in either agarose or Fomblin, is common practice in time-intensive MRI studies of ex vivo biological samples, and is required to prevent tissue autolysis and sample motion. However, the combined effect of fixation and sample embedding may alter tissue structure and MRI properties. We investigated the progressive changes in T(1) and T(2) relaxation times, and the arrangement of locally prevailing cardiomyocyte orientation determined using diffusion tensor imaging, in embedded ex vivo rat hearts fixed using Karnovsky's solution (glutaraldehyde-formaldehyde mix). Three embedding media were investigated: (i) standard agarose (n = 3 hearts); (ii) Fomblin (n = 4 hearts); and (iii) iso-osmotic agarose (n = 3 hearts); in the latter, the osmolarity of the fixative and embedding medium was adjusted to 300 mOsm to match more closely that of native tissue. The T(1) relaxation time in the myocardium showed a pronounced decrease over a 48-h period following embedding in Fomblin (-11.3 ± 6.2%; mean ± standard deviation), but was stable in standard agarose- and iso-osmotic agarose-embedded hearts. The mean myocardial T(2) relaxation time increased in all embedded hearts: by 35.1 ± 14.7% with standard agarose embedding, 13.1 ± 5.6% with Fomblin and 13.3 ± 1.4% with iso-osmotic agarose. Deviation in the orientation of the primary eigenvector of the diffusion tensor occurred in all hearts (mean angular changes of 6.6°, 3.2° and 1.9° per voxel after 48 h in agarose-, Fomblin- and iso-osmotic agarose-embedded hearts, respectively), indicative of progressive structural changes in myocardial histo-architecture, in spite of previous exposure to fast-acting tissue fixation. Our results suggest that progressive structural changes occur in chemically fixed myocardium, and that the extent of these changes is modulated by the embedding medium, and by osmotic gradients between the fixative in the tissue and the surrounding medium.