Relating myocardial laminar architecture to shear strain and muscle fiber orientation

Investigating the three-dimensional myocardial micro-architecture in the laminar structure using X-ray phase-contrast microtomography

A comprehensive grasp of the myocardial micro-architecture is essential for understanding diverse heart functions. This study aimed to investigate three-dimensional (3D) cardiomyocyte arrangement in the laminar structure using X-ray phase-contrast microtomography. Using the ID-19 beamline at the European Synchrotron Radiation Facility, we imaged human left ventricular (LV) wall transparietal samples and reconstructed them with an isotropic voxel edge length of 3.5 μm. From the reconstructed volumes, we extracted different regions to analyze the orientation distribution of local cardiomyocyte aggregates, presenting findings in terms of helix and intrusion angles. In regions containing one sheetlet population, we observed cardiomyocyte aggregates running along the local LV wall's radial direction at the border of sheetlets, branching and merging into a complex network around connecting points of different sheetlets, and bending to accommodate vessel passages. In regions with two sheetlet populations, the helix angle of local cardiomyocyte aggregates experiences a nonmonotonic change, and some cardiomyocyte aggregates run along the local radial direction. X-ray phase-contrast microtomography is a valuable technique for investigating the 3D local myocardial architecture at microscopic level. The arrangement of local cardiomyocyte aggregates in the LV wall proves to be both regional and complex, intricately linked to the local laminar structure.

Full-Field Strain Measurements of the Muscle-Tendon Junction Using X-ray Computed Tomography and Digital Volume Correlation

We report, for the first time, the full-field 3D strain distribution of the muscle-tendon junction (MTJ). Understanding the strain distribution at the junction is crucial for the treatment of injuries and to predict tear formation at this location. Three-dimensional full-field strain distribution of mouse MTJ was measured using X-ray computer tomography (XCT) combined with digital volume correlation (DVC) with the aim of understanding the mechanical behavior of the junction under tensile loading. The interface between the Achilles tendon and the gastrocnemius muscle was harvested from adult mice and stained using 1% phosphotungstic acid in 70% ethanol. In situ XCT combined with DVC was used to image and compute strain distribution at the MTJ under a tensile load (2.4 N). High strain measuring 120,000 µε, 160,000 µε, and 120,000 µε for the first principal stain (εp1), shear strain (γ), and von Mises strain (εVM), respectively, was measured at the MTJ and these values reduced into the body of the muscle or into the tendon. Strain is concentrated at the MTJ, which is at risk of being damaged in activities associated with excessive physical activity.

Oscillation of the orientation of cardiomyocyte aggregates in human left ventricle free wall

Orientation of local cardiomyocyte aggregates in the human left ventricle free wall experiences an oscillation in the laminar structure regions, besides its gradual change trend. We described this oscillation using five transmural samples imaged at the European Synchrotron Radiation Facility with an isotropic voxel size of 3.5 × 3.5 × 3.5 μm³ . In the reconstructed volume of each sample, we manually selected a region containing a regular laminar structure as the region of interest and measured the distribution of the orientation of local cardiomyocyte aggregates inside using a Fourier-based method. Then, we extracted the gradual change part of the orientation of cardiomyocyte aggregates with a three-dimensional centered Gaussian filter and measured the angle between the original orientation vector of local cardiomyocyte aggregates and its gradual change part. Further, we assessed the measured angles in different local coordinates. The results indicate that the oscillation amplitude of the orientation of cardiomyocyte aggregates is regional in the left ventricle wall, which may promote our understanding of the rearrangement mechanism of the cardiomyocyte aggregates and provide a new biomarker to study the heart physiological status.

Myocardial mesostructure and mesofunction

The complex and highly organized structural arrangement of some five billion cardiomyocytes directs the coordinated electrical activity and mechanical contraction of the human heart. The characteristic transmural change in cardiomyocyte orientation underlies base-to-apex shortening, circumferential shortening, and left ventricular torsion during contraction. Individual cardiomyocytes shorten ∼15% and increase in diameter ∼8%. Remarkably, however, the left ventricular wall thickens by up to 30-40%. To accommodate this, the myocardium must undergo significant structural rearrangement during contraction. At the mesoscale, collections of cardiomyocytes are organized into sheetlets, and sheetlet shear is the fundamental mechanism of rearrangement that produces wall thickening. Herein, we review the histological and physiological studies of myocardial mesostructure that have established the sheetlet shear model of wall thickening. Recent developments in tissue clearing techniques allow for imaging of whole hearts at the cellular scale, whereas magnetic resonance imaging (MRI) and computed tomography (CT) can image the myocardium at the mesoscale (100 µm to 1 mm) to resolve cardiomyocyte orientation and organization. Through histology, cardiac diffusion tensor imaging (DTI), and other modalities, mesostructural sheetlets have been confirmed in both animal and human hearts. Recent in vivo cardiac DTI methods have measured reorientation of sheetlets during the cardiac cycle. We also examine the role of pathological cardiac remodeling on sheetlet organization and reorientation, and the impact this has on ventricular function and dysfunction. We also review the unresolved mesostructural questions and challenges that may direct future work in the field.

Multiscale Anisotropic Tissue Biofabrication via Bulk Acoustic Patterning of Cells and Functional Additives in Hybrid Bioinks

Recapitulation of the microstructural organization of cellular and extracellular components found in natural tissues is an important but challenging feat for tissue engineering, which demands innovation across both process and material fronts. In this work, a highly versatile ultrasound-assisted biofabrication (UAB) approach is demonstrated that utilizes radiation forces generated by superimposing ultrasonic bulk acoustic waves to rapidly organize arrays of cells and other biomaterial additives within single and multilayered hydrogel constructs. UAB is used in conjunction with a novel hybrid bioink system, comprising of cartilage-forming cells (human adipose-derived stem cells or chondrocytes) and additives to promote cell adhesion (collagen microaggregates or polycaprolactone microfibers) encapsulated within gelatin methacryloyl (GelMA) hydrogels, to fabricate cartilaginous tissue constructs featuring bulk anisotropy. The hybrid matrices fabricated under the appropriate synergistic thermo-reversible and photocrosslinking conditions demonstrate enhanced mechanical stiffness, stretchability, strength, construct shape fidelity and aligned encapsulated cell morphology and collagen II secretion in long-term culture. Hybridization of UAB is also shown with extrusion and stereolithography printing to fabricate constructs featuring 3D perfusable channels for vasculature combined with a crisscross or circumferential organization of cells and adhesive bioadditives, which is relevant for further translation of UAB toward complex physiological-scale biomimetic tissue fabrication.

Measurement of local orientation of cardiomyocyte aggregates in human left ventricle free wall samples using X-ray phase-contrast microtomography

Most cardiomyocytes in the left ventricle wall are grouped in aggregates of four to five units that are quasi-parallel to each other. When one or more "cardiomyocyte aggregates" are delimited by two cleavage planes, this defines a "sheetlet" that can be considered as a "work unit" that contributes to the thickening of the wall during the cardiac cycle. In this paper, we introduce the skeleton method to measure the local three-dimensional (3D) orientation of cardiomyocyte aggregates in the sheetlets in three steps: data segmentation; extraction of the skeleton of the sheetlets; and calculation of the local orientation of the cardiomyocyte aggregates inside the sheetlets. These data include a series of virtual tissue volumes and five transmural human left ventricle free wall samples, imaged with 3D synchrotron radiation phase-contrast microtomography, and reconstructed with a 3.5×3.5×3.5μm³ voxel size. We computed the local orientation of the cardiomyocyte aggregates inside the sheetlets with a working window of 112×112×112μm³ in size. These data demonstrate that the skeleton method can provide accurate 3D measurements and reliable screening of the 3D evolution of the orientation of cardiomyocyte aggregates within the sheetlets. We showed that in regions that contain one population of quasi-parallel sheetlets, the orientation of the cardiomyocyte aggregates undergo "oscillations" along the perpendicular direction of the sheetlets. In regions that contain two populations of sheetlets with a different angular range, we demonstrate some discontinuity of the helix angle of the cardiomyocyte aggregates at the interface between the two populations.

Architected fibrous scaffolds for engineering anisotropic tissues

Mimicking the native three-dimensional microenvironment is of crucial importance when biofabricating a new healthcare material. One aspect of the native tissue that is often omitted when designing a suitable scaffold is its anisotropy. Not only is matching native mechanical properties important when designing implantable scaffolds or healthcare materials, but matching physiological structure is also important as many cell populations respond differently to fiber orientation. Therefore, novel aligned electrospun scaffolds with varying fiber angles and spacing of bundles were created and mechanically characterized. Through controlling the angle between the fibers in each layer of the scaffold, a range of different physiological anisotropic mechanical properties were achieved that encompasses values found in native tissues. Extrapolation of this mechanical data allowed for any native tissue's anisotropic Young's modulus to be mimicked by electrospinning fibers at a particular angle. These electrospun scaffolds were then incorporated with cell-laden hydrogels to create hybrid structures that contain the benefits of both scaffolding techniques with the ability to encapsulate cells in the hydrogel. To conclude, this study develops a novel bundled fiber scaffold that was architected to yield anisotropic properties matching native tissues.

Does layer-specific strain using speckle tracking echocardiography improve the assessment of left ventricular myocardial deformation? A review

An increasing number of studies of left ventricular myocardial deformation have been published. Layer-specific strain using speckle tracking echocardiography to evaluate left ventricular function is not recommended in clinical practice. However, evaluation of myocardial mechanics using longitudinal and circumferential layer-specific strain enables the detection of subclinical impairment of myocardial deformation in various diseases. Unfortunately, normal values for longitudinal and circumferential strain have not been clearly defined. In normal subjects, layer-specific strain decreases from the endocardial to the epicardial layer, and from the apex to the base of the left ventricle. Although various studies have tried to define normal values for each layer in healthy subjects, studies with more subjects are needed. This tool has good reproducibility in terms of intraobserver and interobserver variability, but, as with monolayer strain, it has poor intervendor variability. Efforts that aim for standardization between vendors will be required before widespread use of this technique can be advocated.

Infarct Collagen Topography Regulates Fibroblast Fate via p38-Yes-Associated Protein Transcriptional Enhanced Associate Domain Signals

RATIONALE

Myocardial infarction causes spatial variation in collagen organization and phenotypic diversity in fibroblasts, which regulate the heart's ECM (extracellular matrix). The relationship between collagen structure and fibroblast phenotype is poorly understood but could provide insights regarding the mechanistic basis for myofibroblast heterogeneity in the injured heart.

OBJECTIVE

To investigate the role of collagen organization in cardiac fibroblast fate determination.

METHODS AND RESULTS

Biomimetic topographies were nanofabricated to recapitulate differential collagen organization in the infarcted mouse heart. Here, adult cardiac fibroblasts were freshly isolated and cultured on ECM topographical mimetics for 72 hours. Aligned mimetics caused cardiac fibroblasts to elongate while randomly organized topographies induced circular morphology similar to the disparate myofibroblast morphologies measured in vivo. Alignment cues also induced myofibroblast differentiation, as >60% of fibroblasts formed αSMA (α-smooth muscle actin) stress fibers and expressed myofibroblast-specific ECM genes like Postn (periostin). By contrast, random organization caused 38% of cardiac fibroblasts to express αSMA albeit with downregulated myofibroblast-specific ECM genes. Coupling topographical cues with the profibrotic agonist, TGFβ (transforming growth factor beta), additively upregulated myofibroblast-specific ECM genes independent of topography, but only fibroblasts on flat and randomly oriented mimetics had increased percentages of fibroblasts with αSMA stress fibers. Increased tension sensation at focal adhesions induced myofibroblast differentiation on aligned mimetics. These signals were transduced by p38-YAP (yes-associated protein)-TEAD (transcriptional enhanced associate domain) interactions, in which both p38 and YAP-TEAD (yes-associated protein transcriptional enhanced associate domain) binding were required for myofibroblast differentiation. By contrast, randomly oriented mimetics did not change focal adhesion tension sensation or enrich for p38-YAP-TEAD interactions, which explains the topography-dependent diversity in fibroblast phenotypes observed here.

CONCLUSIONS

Spatial variations in collagen organization regulate cardiac fibroblast phenotype through mechanical activation of p38-YAP-TEAD signaling, which likely contribute to myofibroblast heterogeneity in the infarcted myocardium.

Strain Curve Classification Using Supervised Machine Learning Algorithm with Physiologic Constraints

Speckle tracking echocardiography (STE) enables quantification of myocardial deformation by a generation of spatiotemporal strain curves or time-strain curves (TSCs). Currently, only assessment of peak global longitudinal strain is employed in clinical practice because of the uncertainty in the accuracy of STE. We describe a supervised machine learning, physiologically constrained, fully automatic algorithm, trained with labeled data, for classification of TSCs into physiologic or artifactual classes. The data set of 415 healthy patients, with three cine loops per patient, corresponding to the three standard 2-D longitudinal views, was processed using a previously published, in-house STE software termed K-SAD. We report an accuracy of 86.4% for classifying TSCs as physiologic, artifactual and undetermined curves. The positive predictive value for a physiologic strain curve is 89%. This is as a necessary step for a similar separation of pathologic conditions, to allow full utilization of the temporal information concealed in layer-specific segmental TSCs.

Anatomically correct assessment of the orientation of the cardiomyocytes using diffusion tensor imaging

Diffusion tensor imaging has been used for assessing the orientation of cardiac myocytes for decades. Striking methodological differences exist between studies when quantifying these orientations. This limits the comparability between studies, and impedes collaboration and the drawing of appropriate physiological conclusions. We have sought to elucidate these differences, permitting us to propose a standardised "tool set" that might better establish consensus in future studies. We fixed hearts from seven 25 kg pigs in formalin, and scanned them using diffusion tensor imaging. Using various angle definitions as found in literature, we assessed the orientations of cardiomyocytes, comparing them in terms of helical and intrusion angles, along with the orientation of their aggregations. The difference between assessment of the helical angle with and without relation to the epicardial curvature was 25.2° (SD: 7.9) at the base, 5.8° (1.9) at the equatorial level, and 28.0° (7.0) at the apex, ANOVA P = 0.001. In comparable fashion, the intrusion angle differed by 25.9° (12.9), 7.6° (0.98) and 17.5° (4.7), P = 0.01, and the angle of the aggregates (E3-angle) differed by 25.0° (13.5) at the base, 9.4° (1.7) at the equator, and 23.1° (6.2) apically, P = 0.003. When assessing 14 definitions used in literature to calculate the orientation of aggregates, only 4 rendered identical results. The findings show that any attempt to use projection of eigenvectors introduces considerable bias. The epicardial curvature of the ventricular cone needs to be taken into account when seeking to provide accurate quantification of the orientation of the aggregated cardiomyocytes, especially in the apical and basal regions. This means that projection of eigenvectors should be avoided prior to quantifying myocyte orientation, especially when assessing radial orientation. Based on our results, we suggest appropriate methods for valid assessment of myocyte orientation using diffusion tensor imaging.

Engineering hiPSC cardiomyocyte in vitro model systems for functional and structural assessment

The study of human cardiomyopathies and the development and testing of new therapies has long been limited by the availability of appropriate in vitro model systems. Cardiomyocytes are highly specialized cells whose internal structure and contractile function are sensitive to the local microenvironment and the combination of mechanical and biochemical cues they receive. The complementary technologies of human induced pluripotent stem cell (hiPSC) derived cardiomyocytes (CMs) and microphysiological systems (MPS) allow for precise control of the genetics and microenvironment of human cells in in vitro contexts. These combined systems also enable quantitative measurement of mechanical function and intracellular organization. This review describes relevant factors in the myocardium microenvironment that affect CM structure and mechanical function and demonstrates the application of several engineered microphysiological systems for studying development, disease, and drug discovery.

A Contemporary Look at Biomechanical Models of Myocardium

Understanding and predicting the mechanical behavior of myocardium under healthy and pathophysiological conditions are vital to developing novel cardiac therapies and promoting personalized interventions. Within the past 30 years, various constitutive models have been proposed for the passive mechanical behavior of myocardium. These models cover a broad range of mathematical forms, microstructural observations, and specific test conditions to which they are fitted. We present a critical review of these models, covering both phenomenological and structural approaches, and their relations to the underlying structure and function of myocardium. We further explore the experimental and numerical techniques used to identify the model parameters. Next, we provide a brief overview of continuum-level electromechanical models of myocardium, with a focus on the methods used to integrate the active and passive components of myocardial behavior. We conclude by pointing to future directions in the areas of optimal form as well as new approaches for constitutive modeling of myocardium.

Imaging technologies for cardiac fiber and heart failure: a review

There has been an increasing interest in studying cardiac fibers in order to improve the current knowledge regarding the mechanical and physiological properties of the heart during heart failure (HF), particularly early HF. Having a thorough understanding of the changes in cardiac fiber orientation may provide new insight into the mechanisms behind the progression of left ventricular (LV) remodeling and HF. We conducted a systematic review on various technologies for imaging cardiac fibers and its link to HF. This review covers literature reports from 1900 to 2017. PubMed and Google Scholar databases were searched using the keywords "cardiac fiber" and "heart failure" or "myofiber" and "heart failure." This review highlights imaging methodologies, including magnetic resonance diffusion tensor imaging (MR-DTI), ultrasound, and other imaging technologies as well as their potential applications in basic and translational research on the development and progression of HF. MR-DTI and ultrasound have been most useful and significant in evaluating cardiac fibers and HF. New imaging technologies that have the ability to measure cardiac fiber orientations and identify structural and functional information of the heart will advance basic research and clinical diagnoses of HF.

Kinematic boundary conditions substantially impact in silico ventricular function

Computational cardiac mechanical models, individualized to the patient, have the potential to elucidate the fundamentals of cardiac (patho-)physiology, enable non-invasive quantification of clinically significant metrics (eg, stiffness, active contraction, work), and anticipate the potential efficacy of therapeutic cardiovascular intervention. In a clinical setting, however, the available imaging resolution is often limited, which limits cardiac models to focus on the ventricles, without including the atria, valves, and proximal arteries and veins. In such models, the absence of surrounding structures needs to be accounted for by imposing realistic kinematic boundary conditions, which, for prognostic purposes, are preferably generic and thus non-image derived. Unfortunately, the literature on cardiac models shows no consistent approach to kinematically constrain the myocardium. The impact of different approaches (eg, fully constrained base, constrained epi-ring) on the predictive capacity of cardiac mechanical models has not been thoroughly studied. For that reason, this study first gives an overview of current approaches to kinematically constrain (bi) ventricular models. Next, we developed a patient-specific in silico biventricular model that compares well with literature and in vivo recorded strains. Alternative constraints were introduced to assess the influence of commonly used mechanical boundary conditions on both the predicted global functional behavior of the in-silico heart (cavity volumes, stroke volume, ejection fraction) and local strain distributions. Meaningful differences in global functioning were found between different kinematic anchoring strategies, which brought forward the importance of selecting appropriate boundary conditions for biventricular models that, in the near future, may inform clinical intervention. However, whilst statistically significant differences were also found in local strain distributions, these differences were minor and mostly confined to the region close to the applied boundary conditions.

A tissue-engineered scale model of the heart ventricle

Laboratory studies of the heart use cell and tissue cultures to dissect heart function yet rely on animal models to measure pressure and volume dynamics. Here, we report tissue-engineered scale models of the human left ventricle, made of nanofibrous scaffolds that promote native-like anisotropic myocardial tissue genesis and chamber-level contractile function. Incorporating neonatal rat ventricular myocytes or cardiomyocytes derived from human induced pluripotent stem cells, the tissue-engineered ventricles have a diastolic chamber volume of ~500 µl (comparable to that of the native rat ventricle and approximately 1/250 the size of the human ventricle), and ejection fractions and contractile work 50-250 times smaller and 104-108 times smaller than the corresponding values for rodent and human ventricles, respectively. We also measured tissue coverage and alignment, calcium-transient propagation and pressure-volume loops in the presence or absence of test compounds. Moreover, we describe an instrumented bioreactor with ventricular-assist capabilities, and provide a proof-of-concept disease model of structural arrhythmia. The model ventricles can be evaluated with the same assays used in animal models and in clinical settings.

In vivo estimation of passive biomechanical properties of human myocardium

Identification of in vivo passive biomechanical properties of healthy human myocardium from regular clinical data is essential for subject-specific modelling of left ventricle (LV). In this work, myocardium was defined by Holzapfel-Ogden constitutive law. Therefore, the objectives of the study were (a) to estimate the ranges of the constitutive parameters for healthy human myocardium using non-invasive routine clinical data, and (b) to investigate the effect of geometry, LV end-diastolic pressure (EDP) and fibre orientations on estimated values. In order to avoid invasive measurements and additional scans, LV cavity volume, measured from routine MRI, and empirical pressure-normalised-volume relation (Klotz-curve) were used as clinical data. Finite element modelling, response surface method and genetic algorithm were used to inversely estimate the constitutive parameters. Due to the ill-posed nature of the inverse optimisation problem, the myocardial properties was extracted by identifying the ranges of the parameters, instead of finding unique values. Additional sensitivity studies were carried out to identify the effect of LV EDP, fibre orientation and geometry on estimated parameters. Although uniqueness of the solution cannot be achieved, the normal ranges of the parameters produced similar mechanical responses within the physiological ranges. These information could be used in future computational studies for designing heart failure treatments. Graphical abstract.

In vivo estimation of passive biomechanical properties of human myocardium

Identification of in vivo passive biomechanical properties of healthy human myocardium from regular clinical data is essential for subject-specific modelling of left ventricle (LV). In this work, myocardium was defined by Holzapfel-Ogden constitutive law. Therefore, the objectives of the study were (a) to estimate the ranges of the constitutive parameters for healthy human myocardium using non-invasive routine clinical data, and (b) to investigate the effect of geometry, LV end-diastolic pressure (EDP) and fibre orientations on estimated values. In order to avoid invasive measurements and additional scans, LV cavity volume, measured from routine MRI, and empirical pressure-normalised-volume relation (Klotz-curve) were used as clinical data. Finite element modelling, response surface method and genetic algorithm were used to inversely estimate the constitutive parameters. Due to the ill-posed nature of the inverse optimisation problem, the myocardial properties was extracted by identifying the ranges of the parameters, instead of finding unique values. Additional sensitivity studies were carried out to identify the effect of LV EDP, fibre orientation and geometry on estimated parameters. Although uniqueness of the solution cannot be achieved, the normal ranges of the parameters produced similar mechanical responses within the physiological ranges. These information could be used in future computational studies for designing heart failure treatments.

Assessment of Myocardial Microstructural Dynamics by In Vivo Diffusion Tensor Cardiac Magnetic Resonance

BACKGROUND

Cardiomyocytes are organized in microstructures termed sheetlets that reorientate during left ventricular thickening. Diffusion tensor cardiac magnetic resonance (DT-CMR) may enable noninvasive interrogation of in vivo cardiac microstructural dynamics. Dilated cardiomyopathy (DCM) is a condition of abnormal myocardium with unknown sheetlet function.

OBJECTIVES

This study sought to validate in vivo DT-CMR measures of cardiac microstructure against histology, characterize microstructural dynamics during left ventricular wall thickening, and apply the technique in hypertrophic cardiomyopathy (HCM) and DCM.

METHODS

In vivo DT-CMR was acquired throughout the cardiac cycle in healthy swine, followed by in situ and ex vivo DT-CMR, then validated against histology. In vivo DT-CMR was performed in 19 control subjects, 19 DCM, and 13 HCM patients.

RESULTS

In swine, a DT-CMR index of sheetlet reorientation (E2A) changed substantially (E2A mobility ~46°). E2A changes correlated with wall thickness changes (in vivo r² = 0.75; in situ r² = 0.89), were consistently observed under all experimental conditions, and accorded closely with histological analyses in both relaxed and contracted states. The potential contribution of cyclical strain effects to in vivo E2A was ~17%. In healthy human control subjects, E2A increased from diastole (18°) to systole (65°; p < 0.001; E2A mobility = 45°). HCM patients showed significantly greater E2A in diastole than control subjects did (48 ; p < 0.001) with impaired E2A mobility (23°; p < 0.001). In DCM, E2A was similar to control subjects in diastole, but systolic values were markedly lower (40° ; p < 0.001) with impaired E2A mobility (20°; p < 0.001).

CONCLUSIONS

Myocardial microstructure dynamics can be characterized by in vivo DT-CMR. Sheetlet function was abnormal in DCM with altered systolic conformation and reduced mobility, contrasting with HCM, which showed reduced mobility with altered diastolic conformation. These novel insights significantly improve understanding of contractile dysfunction at a level of noninvasive interrogation not previously available in humans. (J Am Coll Cardiol 2017;69:661–76) Published by Elsevier on behalf of the American College of Cardiology Foundation.

Imaging the dynamics of cardiac fiber orientation in vivo using 3D Ultrasound Backscatter Tensor Imaging

The assessment of myocardial fiber disarray is of major interest for the study of the progression of myocardial disease. However, time-resolved imaging of the myocardial structure remains unavailable in clinical practice. In this study, we introduce 3D Backscatter Tensor Imaging (3D-BTI), an entirely novel ultrasound-based imaging technique that can map the myocardial fibers orientation and its dynamics with a temporal resolution of 10 ms during a single cardiac cycle, non-invasively and in vivo in entire volumes. 3D-BTI is based on ultrafast volumetric ultrasound acquisitions, which are used to quantify the spatial coherence of backscattered echoes at each point of the volume. The capability of 3D-BTI to map the fibers orientation was evaluated in vitro in 5 myocardial samples. The helicoidal transmural variation of fiber angles was in good agreement with the one obtained by histological analysis. 3D-BTI was then performed to map the fiber orientation dynamics in vivo in the beating heart of an open-chest sheep at a volume rate of 90 volumes/s. Finally, the clinical feasibility of 3D-BTI was shown on a healthy volunteer. These initial results indicate that 3D-BTI could become a fully non-invasive technique to assess myocardial disarray at the bedside of patients.

Passive diastolic modelling of human ventricles: Effects of base movement and geometrical heterogeneity

Left-ventricular (LV) remodelling, associated with diastolic heart failure, is driven by an increase in myocardial stress. Therefore, normalisation of LV wall stress is the cornerstone of many therapeutic treatments. However, information regarding such regional stress-strain for human LV is still limited. Thus, the objectives of our study were to determine local diastolic stress-strain field in healthy LVs, and consequently, to identify the regional variations amongst them due to geometric heterogeneity. Effects of LV base movement on diastolic model predictions, which were ignored in the literature, were further explored. Personalised finite-element modelling of five normal human bi-ventricles was carried out using subject-specific myocardium properties. Model prediction was validated individually through comparison with end-diastolic volume and a new shape-volume based measurement of LV cavity, extracted from magnetic resonance imaging. Results indicated that incorporation of LV base movement improved the model predictions (shape-volume relevancy of LV cavity), and therefore, it should be considered in future studies. The LV endocardium always experienced higher fibre stress compared to the epicardium for all five subjects. The LV wall near base experienced higher stress compared to equatorial and apical locations. The lateral LV wall underwent greater stress distribution (fibre and sheet stress) compared to other three regions. In addition, normal ranges of different stress-strain components in different regions of LV wall were reported for five healthy ventricles. This information could be used as targets for future computational studies to optimise diastolic heart failure treatments or design new therapeutic interventions/devices.

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.

Quantification of Coupled Stiffness and Fiber Orientation Remodeling in Hypertensive Rat Right-Ventricular Myocardium Using 3D Ultrasound Speckle Tracking with Biaxial Testing

Mechanical and structural changes of right ventricular (RV) in response to pulmonary hypertension (PH) are inadequately understood. While current standard biaxial testing provides information on the mechanical behavior of RV tissues using surface markers, it is unable to fully assess structural and mechanical properties across the full tissue thickness. In this study, the mechanical and structural properties of normotensive and pulmonary hypertension right ventricular (PHRV) myocardium through its full thickness were examined using mechanical testing combined with 3D ultrasound speckle tracking (3D-UST). RV pressure overload was induced in Sprague–Dawley rats by pulmonary artery (PA) banding. The second Piola–Kirchhoff stress tensors and Green-Lagrangian strain tensors were computed in the RV myocardium using the biaxial testing combined with 3D-UST. A previously established non-linear curve-fitting algorithm was applied to fit experimental data to a Strain Energy Function (SEF) for computation of myofiber orientation. The fiber orientations obtained by the biaxial testing with 3D-UST compared well with the fiber orientations computed from the histology. In addition, the re-orientation of myofiber in the right ventricular free wall (RVFW) along longitudinal direction (apex-to-outflow-tract direction) was noticeable in response to PH. For normotensive RVFW samples, the average fiber orientation angles obtained by 3D-UST with biaxial test spiraled from 20° at the endo-cardium to -42° at the epi-cardium (Δ = 62°). For PHRV samples, the average fiber orientation angles obtained by 3D-UST with biaxial test had much less spiral across tissue thickness: 3° at endo-cardium to -7° at epi-cardium (Δ = 10°, P<0.005 compared to normotensive).

Resolving Fine Cardiac Structures in Rats with High-Resolution Diffusion Tensor Imaging

Cardiac architecture is fundamental to cardiac function and can be assessed non-invasively with diffusion tensor imaging (DTI). Here, we aimed to overcome technical challenges in ex vivo DTI in order to extract fine anatomical details and to provide novel insights in the 3D structure of the heart. An integrated set of methods was implemented in ex vivo rat hearts, including dynamic receiver gain adjustment, gradient system scaling calibration, prospective adjustment of diffusion gradients, and interleaving of diffusion-weighted and non-diffusion-weighted scans. Together, these methods enhanced SNR and spatial resolution, minimised orientation bias in diffusion-weighting, and reduced temperature variation, enabling detection of tissue structures such as cell alignment in atria, valves and vessels at an unprecedented level of detail. Improved confidence in eigenvector reproducibility enabled tracking of myolaminar structures as a basis for segmentation of functional groups of cardiomyocytes. Ex vivo DTI facilitates acquisition of high quality structural data that complements readily available in vivo cardiac functional and anatomical MRI. The improvements presented here will facilitate next generation virtual models integrating micro-structural and electro-mechanical properties of the heart.

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.

Heterogeneity of Fractional Anisotropy and Mean Diffusivity Measurements by In Vivo Diffusion Tensor Imaging in Normal Human Hearts

Background

Cardiac diffusion tensor imaging (cDTI) by cardiovascular magnetic resonance has the potential to assess microstructural changes through measures of fractional anisotropy (FA) and mean diffusivity (MD). However, normal variation in regional and transmural FA and MD is not well described.

Methods

Twenty normal subjects were scanned using an optimised cDTI sequence at 3T in systole. FA and MD were quantified in 3 transmural layers and 4 regional myocardial walls.

Results

FA was higher in the mesocardium (0.46 ±0.04) than the endocardium (0.40 ±0.04, p≤0.001) and epicardium (0.39 ±0.04, p≤0.001). On regional analysis, the FA in the septum was greater than the lateral wall (0.44 ±0.03 vs 0.40 ±0.05 p = 0.04). There was a transmural gradient in MD increasing towards the endocardium (epicardium 0.87 ±0.07 vs endocardium 0.91 ±0.08×10^-3 mm²/s, p = 0.04). With the lateral wall (0.87 ± 0.08×10^-3 mm²/s) as the reference, the MD was higher in the anterior wall (0.92 ±0.08×10^-3 mm²/s, p = 0.016) and septum (0.92 ±0.07×10^-3 mm²/s, p = 0.028). Transmurally the signal to noise ratio (SNR) was greatest in the mesocardium (14.5 ±2.5 vs endocardium 13.1 ±2.2, p<0.001; vs epicardium 12.0 ± 2.4, p<0.001) and regionally in the septum (16.0 ±3.4 vs lateral wall 11.5 ± 1.5, p<0.001). Transmural analysis suggested a relative reduction in the rate of change in helical angle (HA) within the mesocardium.

Conclusions

In vivo FA and MD measurements in normal human heart are heterogeneous, varying significantly transmurally and regionally. Contributors to this heterogeneity are many, complex and interactive, but include SNR, variations in cardiac microstructure, partial volume effects and strain. These data indicate that the potential clinical use of FA and MD would require measurement standardisation by myocardial region and layer, unless pathological changes substantially exceed the normal variation identified.

Methods for using 3-D ultrasound speckle tracking in biaxial mechanical testing of biological tissue samples

Being multilayered and anisotropic, biological tissues such as cardiac and arterial walls are structurally complex, making the full assessment and understanding of their mechanical behavior challenging. Current standard mechanical testing uses surface markers to track tissue deformations and does not provide deformation data below the surface. In the study described here, we found that combining mechanical testing with 3-D ultrasound speckle tracking could overcome this limitation. Rat myocardium was tested with a biaxial tester and was concurrently scanned with high-frequency ultrasound in three dimensions. The strain energy function was computed from stresses and strains using an iterative non-linear curve-fitting algorithm. Because the strain energy function consists of terms for the base matrix and for embedded fibers, spatially varying fiber orientation was also computed by curve fitting. Using finite-element simulations, we first validated the accuracy of the non-linear curve-fitting algorithm. Next, we compared experimentally measured rat myocardium strain energy function values with those in the literature and found a matching order of magnitude. Finally, we retained samples after the experiments for fiber orientation quantification using histology and found that the results satisfactorily matched those computed in the experiments. We conclude that 3-D ultrasound speckle tracking can be a useful addition to traditional mechanical testing of biological tissues and may provide the benefit of enabling fiber orientation computation.

Assessment and comparison of left ventricular shear in normal and situs inversus totalis hearts by means of magnetic resonance tagging

Situs inversus totalis (SIT) is characterized by complete mirroring of gross cardiac anatomy and position combined with an incompletely mirrored myofiber arrangement, being normal at the apex but inverted at the base of the left ventricle (LV). This study relates myocardial structure to mechanical function by analyzing and comparing myocardial deformation patterns of normal and SIT subjects, focusing especially on circumferential-radial shear. In nine control and nine SIT normotensive human subjects, myocardial deformation was assessed from magnetic resonance tagging (MRT) image sequences of five LV short-axis slices. During ejection, no significant difference in either circumferential shortening (εcc) or its axial gradient (Δεcc) is found between corresponding LV levels in control and SIT hearts. Circumferential-radial shear (εcr) has a clear linear trend from apex-to-base in controls, while in SIT it hovers close to zero at all levels. Torsion as well as axial change in εcr (Δεcr) is as in controls in apical sections of SIT hearts but deviates significantly towards the base, changing sign close to the LV equator. Interindividual variability in torsion and Δεcr values is higher in SIT than in controls. Apex-to-base trends of torsion and Δεcr in SIT, changing sign near the LV equator, further substantiate a structural transition in myofiber arrangement close to the LV equator itself. Invariance of εcc and Δεcc patterns between controls and SIT subjects shows that normal LV pump function is achieved in SIT despite partial mirroring of myocardial structure leading to torsional and shear patterns that are far from normality.

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.

Ultrasound backscatter tensor imaging (BTI): analysis of the spatial coherence of ultrasonic speckle in anisotropic soft tissues

The assessment of fiber architecture is of major interest in the progression of myocardial disease. Recent techniques such as magnetic resonance diffusion tensor imaging (MR-DTI) or ultrasound elastic tensor imaging (ETI) can derive the fiber directions by measuring the anisotropy of water diffusion or tissue elasticity, but these techniques present severe limitations in a clinical setting. In this study, we propose a new technique, backscatter tensor imaging (BTI), which enables determination of the fiber directions in skeletal muscles and myocardial tissues, by measuring the spatial coherence of ultrasonic speckle. We compare the results to ultrasound ETI. Acquisitions were performed using a linear transducer array connected to an ultrasonic scanner mounted on a motorized rotation device with angles from 0° to 355° by 5° increments to image ex vivo bovine skeletal muscle and porcine left ventricular myocardial samples. At each angle, multiple plane waves were transmitted and the backscattered echoes recorded. The coherence factor was measured as the ratio of coherent intensity over incoherent intensity of backscattered echoes. In skeletal muscle, maximal/minimal coherence factor was found for the probe parallel/perpendicular to the fibers. In myocardium, the coherence was assessed across the entire myocardial thickness, and the position of maxima and minima varied transmurally because of the complex fibers distribution. In ETI, the shear wave speed variation with the probe angle was found to follow the coherence variation. Spatial coherence can thus reveal the anisotropy of the ultrasonic speckle in skeletal muscle and myocardium. BTI could be used on any type of ultrasonic scanner with rotating phased-array probes or 2-D matrix probes for noninvasive evaluation of myocardial fibers.

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.

Progression of myocardial remodeling and mechanical dysfunction in the spontaneously hypertensive rat

The progression of hypertensive heart disease (HHD) to heart failure (HF) is associated with myocardial remodeling. Corresponding changes in three-dimensional organization of cardiac extracellular matrix have not been quantified or related fully to the development of HF. Spontaneously hypertensive rats (SHRs) and Wistar-Kyoto controls were studied at 3, 12, 18, and 24 mo. Hemodynamic and morphological data, brain natriuretic peptide levels, and echocardiography demonstrate four distinct disease stages: systemic hypertension, diastolic dysfunction, early systolic failure, and decompensated HF. Passive left ventricular (LV) pressure-volume relationships were determined in vitro. Transmural specimens from the anterior LV free wall were imaged using extended-volume confocal microscopy, and three-dimensional myocardial architecture was quantified. In SHRs, LV compliance was reduced at 12 mo and increased progressively thereafter. However, it was less than in controls for filling pressures <10 mmHg and not significantly different at ≥10 mmHg. Myocyte cross section was enlarged, with increased variability from 12 mo, while collagen fraction increased progressively. Perimysial collagen fraction remained unchanged with age, although endomysial collagen increased from 12 mo. Perimysial collagen between adjacent muscle layers fused at 12 mo and continued to thicken subsequently, while muscle layers became more dispersed and disordered. We conclude that LV dilatation, which accompanies decompensated HF in this model of HHD, is not due to LV "softening." While perimysial (and endomysial) collagen networks are substantially remodeled, they are not dissolved, as has been proposed. We argue that progressive disruption of the laminar organization of LV myocardium may contribute to impaired systolic function in HHD.

Histo-anatomical structure of the living isolated rat heart in two contraction states assessed by diffusion tensor MRI

Deformation and wall-thickening of ventricular myocardium are essential for cardiac pump function. However, insight into the histo-anatomical basis for cardiac tissue re-arrangement during contraction is limited. In this report, we describe dynamic changes in regionally prevailing cardiomyocyte (fibre) and myolaminar (sheet) orientations, using Diffusion Tensor Imaging (DTI) of ventricles in the same living heart in two different mechanical states. Hearts, isolated from Sprague–Dawley rats, were Langendorff-perfused and imaged, initially in their slack state during cardioplegic arrest, then during lithium-induced contracture. Regional fibre- and sheet-orientations were derived from DTI-data on a voxel-wise basis. Contraction was accompanied with a decrease in left-handed helical fibres (handedness relative to the baso-apical direction) in basal, equatorial, and apical sub-epicardium (by 14.0%, 17.3%, 15.8% respectively; p < 0.001), and an increase in right-handed helical fibres of the sub-endocardium (by 11.0%, 12.1% and 16.1%, respectively; p < 0.001). Two predominant sheet-populations were observed, with sheet-angles of either positive (β+) or negative (β−) polarity relative to a ‘chamber-horizontal plane’ (defined as normal to the left ventricular long-axis). In contracture, mean ‘intersection’-angle (geometrically quantifiable intersection of sheet-angle projections) between β+ and β− sheet-populations increased from 86.2 ± 5.5° (slack) to 108.3 ± 5.4° (p < 0.001). Subsequent high-resolution DTI of fixed myocardium, and histological sectioning, reconfirmed the existence of alternating sheet-plane populations. Our results suggest that myocardial tissue layers in alternating sheet-populations align into a more chamber-horizontal orientation during contraction. This re-arrangement occurs via an accordion-like mechanism that, combined with inter-sheet slippage, can significantly contribute to ventricular deformation, including wall-thickening in a predominantly centripetal direction and baso-apical shortening.

Why SIT works: normal function despite typical myofiber pattern in Situs Inversus Totalis (SIT) hearts derived by shear-induced myofiber reorientation

The left ventricle (LV) of mammals with Situs Solitus (SS, normal organ arrangement) displays hardly any interindividual variation in myofiber pattern and experimentally determined torsion. SS LV myofiber pattern has been suggested to result from adaptive myofiber reorientation, in turn leading to efficient pump and myofiber function. Limited data from the Situs Inversus Totalis (SIT, a complete mirror image of organ anatomy and position) LV demonstrated an essential different myofiber pattern, being normal at the apex but mirrored at the base. Considerable differences in torsion patterns in between human SIT LVs even suggest variation in myofiber pattern among SIT LVs themselves. We addressed whether different myofiber patterns in the SIT LV can be predicted by adaptive myofiber reorientation and whether they yield similar pump and myofiber function as in the SS LV. With a mathematical model of LV mechanics including shear induced myofiber reorientation, we predicted myofiber patterns of one SS and three different SIT LVs. Initial conditions for SIT were based on scarce information on the helix angle. The transverse angle was set to zero. During reorientation, a non-zero transverse angle developed, pump function increased, and myofiber function increased and became more homogeneous. Three continuous SIT structures emerged with a different location of transition between normal and mirrored myofiber orientation pattern. Predicted SIT torsion patterns matched experimentally determined ones. Pump and myofiber function in SIT and SS LVs are similar, despite essential differences in myocardial structure. SS and SIT LV structure and function may originate from same processes of adaptive myofiber reorientation.

Deciphering the structure-function relation in healthy hearts is important to understand cardiac pathologies. In the structure-function relation, the myofiber orientation patterns play a central role. Between people with normal organ arrangement (Situs Solitus, SS) this pattern is strikingly similar. Such consistency in myocardial structure might be the result of an adaptation process to accommodate for homogeneous distribution of myofiber strain across the wall and for optimal pump function. The heart of people with a mirror-imaged position of their organs (Situs Inversus Totalis, SIT) has a modified myofiber orientation pattern with respect to SS: normal at the LV apex, but mirrored at the base. Hence, studying SIT hearts provides a unique possibility 1) for understanding adaptation mechanisms related to myofiber orientation and mechanical load, and 2) to gain additional insights into the structure-function relations of the LV. Through mathematical modeling of LV mechanics, we have found that myofiber orientation pattern in both SS and SIT may originate from same processes of adaptive myofiber reorientation. After reorientation, pump and local myofiber function were found to be similar between SS and SIT as well: a remarkable finding when considering the large difference in myofiber orientation pattern.

Anisotropic reinforcement of acute anteroapical infarcts improves pump function

BACKGROUND

We hypothesize that a therapy that improves left ventricular (LV) pump function early after infarction should decrease the need for compensation through sympathetic activation and dilation, thereby reducing the risk of developing heart failure. The mechanical properties of healing myocardial infarcts are an important determinant of LV function, yet improving function by altering infarct properties has proven unexpectedly difficult. Using a computational model, we recently predicted that stiffening a large anterior infarct anisotropically (in only one direction) would improve LV function, whereas isotropic stiffening, the focus of previous studies and therapies, would not. The goal of this study was to test the novel strategy of anisotropic infarct reinforcement.

METHODS AND RESULTS

We tested the effects of anisotropic infarct reinforcement in 10 open-chest dogs with large anteroapical infarcts that depressed LV pump function. We measured regional mechanics, LV volumes, and cardiac output at a range of preloads at baseline, 45 minutes after coronary ligation (ischemia), and 30 minutes later, after surgical reinforcement in the longitudinal direction (anisotropic). Ischemia shifted the end-systolic pressure-volume relationship and cardiac output curves rightward, decreasing cardiac output at matched end-diastolic pressure by 44%. Anisotropic reinforcement significantly improved systolic function without impairing diastolic function, recovering half the deficit in overall LV function.

CONCLUSIONS

We conclude that anisotropic reinforcement is a promising new approach to improving LV function after a large myocardial infarction.

Orthotropic active strain models for the numerical simulation of cardiac biomechanics

A model for the active deformation of cardiac tissue considering orthotropic constitutive laws is introduced and studied. In particular, the passive mechanical properties of the myocardium are described by the Holzapfel-Ogden relation, whereas the activation model is based on the concept of active strain. There, an incompatible intermediate configuration is considered, which entails a multiplicative decomposition between active and passive deformation gradients. The underlying Euler-Lagrange equations for minimizing the total energy are written in terms of these deformation factors, where the active part is assumed to depend, at the cell level, on the electrodynamics and on the specific orientation of the cardiomyocytes. The active strain formulation is compared with the classical active stress model from both numerical and modeling perspectives. The well-posedness of the linear system derived from a generic Newton iteration of the original problem is analyzed, and different mechanical activation functions are considered. Taylor-Hood and MINI finite elements are used in the discretization of the overall mechanical problem. The results of several numerical experiments show that the proposed formulation is mathematically consistent and is able to represent the main features of the phenomenon, while allowing savings in computational costs.

Kinematics of cardiac growth: in vivo characterization of growth tensors and strains

Progressive growth and remodeling of the left ventricle are part of the natural history of chronic heart failure and strong clinical indicators for survival. Accompanied by changes in cardiac form and function, they manifest themselves in alterations of cardiac strains, fiber stretches, and muscle volume. Recent attempts to shed light on the mechanistic origin of heart failure utilize continuum theories of growth to predict the maladaptation of the heart in response to pressure or volume overload. However, despite a general consensus on the representation of growth through a second order tensor, the precise format of this growth tensor remains unknown. Here we show that infarct-induced cardiac dilation is associated with a chronic longitudinal growth, accompanied by a chronic thinning of the ventricular wall. In controlled in vivo experiments throughout a period of seven weeks, we found that the lateral left ventricular wall adjacent to the infarct grows longitudinally by more than 10%, thins by more than 25%, lengthens in fiber direction by more than 5%, and decreases its volume by more than 15%. Our results illustrate how a local loss of blood supply induces chronic alterations in structure and function in adjacent regions of the ventricular wall. We anticipate our findings to be the starting point for a series of in vivo studies to calibrate and validate constitutive models for cardiac growth. Ultimately, these models could be useful to guide the design of novel therapies, which allow us to control the progression of heart failure.

Mapping myocardial fiber orientation using echocardiography-based shear wave imaging

The assessment of disrupted myocardial fiber arrangement may help to understand and diagnose hypertrophic or ischemic cardiomyopathy. We hereby proposed and developed shear wave imaging (SWI), which is an echocardiography-based, noninvasive, real-time, and easy-to-use technique, to map myofiber orientation. Five in vitro porcine and three in vivo open-chest ovine hearts were studied. Known in physics, shear wave propagates faster along than across the fiber direction. SWI is a technique that can generate shear waves travelling in different directions with respect to each myocardial layer. SWI further analyzed the shear wave velocity across the entire left-ventricular (LV) myocardial thickness, ranging between 10 (diastole) and 25 mm (systole), with a resolution of 0.2 mm in the middle segment of the LV anterior wall region. The fiber angle at each myocardial layer was thus estimated by finding the maximum shear wave speed. In the in vitro porcine myocardium (n=5) , the SWI-estimated fiber angles gradually changed from +80° ± 7° (endocardium) to +30° ± 13° (midwall) and -40° ± 10° (epicardium) with 0° aligning with the circumference of the heart. This transmural fiber orientation was well correlated with histology findings. SWI further succeeded in mapping the transmural fiber orientation in three beating ovine hearts in vivo. At midsystole, the average fiber orientation exhibited 71° ± 13° (endocardium), 27° ± 8° (midwall), and -26° ± 30° (epicardium). We demonstrated the capability of SWI in mapping myocardial fiber orientation in vitro and in vivo. SWI may serve as a new tool for the noninvasive characterization of myocardial fiber structure.

Clinical benefits of integrating cardiac and vascular models

INTRODUCTION

Recent advances in computational methods and medical imaging techniques have enabled non-invasive exploration of cardiovascular pathologies, from cardiac level to complex arterial networks. The potential of cardiac and vascular modeling in guiding and monitoring therapies could be further extended through the integration of the two systems.

AREAS COVERED

This review includes advances in methods for cardiac electromechanics and vascular flow simulations. The results of a literature search depicting the state of the art in cardiac and vascular modeling are reviewed. The paper goes on to address the benefits and challenges of combined cardiovascular modeling, highlighting the relevance of specific cardiovascular features and implementation. Various alternative approaches and insights on future directions are presented and analyzed with respect to their applicability to clinical practice.

EXPERT OPINION

The article has emerged from the exploration of currently available cardiac and vascular mathematical tools and their corresponding clinical application. The summarized analysis suggests that future efforts should be aimed at developing more accurate and patient-specific mathematical models integrating cardiac and vascular functions to enhance the knowledge of cardiovascular pathologies.

Visualization and quantification of whole rat heart laminar structure using high-spatial resolution contrast-enhanced MRI

It has been shown by histology that cardiac myocytes are organized into laminae and this structure is important in function, both influencing the spread of electrical activation and enabling myocardial thickening in systole by laminar sliding. We have carried out high-spatial resolution three-dimensional MRI of the ventricular myolaminae of the entire volume of the isolated rat heart after contrast perfusion [dimeglumine gadopentate (Gd-DTPA)]. Four ex vivo rat hearts were perfused with Gd-DTPA and fixative and high-spatial resolution MRI was performed on a 9.4T MRI system. After MRI, cryosectioning followed by histology was performed. Images from MRI and histology were aligned, described, and quantitatively compared. In the three-dimensional MR images we directly show the presence of laminae and demonstrate that these are highly branching and are absent from much of the subepicardium. We visualized these MRI volumes to demonstrate laminar architecture and quantitatively demonstrated that the structural features observed are similar to those imaged in histology. We showed qualitatively and quantitatively that laminar architecture is similar in the four hearts. MRI can be used to image the laminar architecture of ex vivo hearts in three dimensions, and the images produced are qualitatively and quantitatively comparable with histology. We have demonstrated in the rat that: 1) laminar architecture is consistent between hearts; 2) myolaminae are absent from much of the subepicardium; and 3) although localized orthotropy is present throughout the myocardium, tracked myolaminae are branching structures and do not have a discrete identity.

The presence of two local myocardial sheet populations confirmed by diffusion tensor MRI and histological validation

PURPOSE

To establish the correspondence between the two histologically observable and diffusion tensor MRI (DTMRI) measurements of myolaminae orientation for the first time and show that single myolaminar orientations observed in local histology may result from histological artifact.

MATERIALS AND METHODS

DTMRI was performed on six sheep left ventricles (LV), then corresponding direct histological transmural measurements were made within the anterobasal and lateral-equatorial LV. Secondary and tertiary eigenvectors of the diffusion tensor were compared with each of the two locally observable sheet orientations from histology. Diffusion tensor invariants were calculated to compare differences in microstructural diffusive properties between histological locations with one observable sheet population and two observable sheet populations.

RESULTS

Mean difference ± 1SD between DTMRI and histology measured sheet angles was 8° ± 27°. Diffusion tensor invariants showed no significant differences between histological locations with one observable sheet population and locations with two observable sheet populations.

CONCLUSION

DTMRI measurements of myolaminae orientations derived from the secondary and tertiary eigenvectors correspond to each of the two local myolaminae orientations observed in histology. Two local sheet populations may exist throughout LV myocardium, and one local sheet population observed in histology may be a result of preparation artifact.

Composite scaffold provides a cell delivery platform for cardiovascular repair

Control over cell engraftment, survival, and function remains critical for heart repair. We have established a tissue engineering platform for the delivery of human mesenchymal progenitor cells (MPCs) by a fully biological composite scaffold. Specifically, we developed a method for complete decellularization of human myocardium that leaves intact most elements of the extracellular matrix, as well as the underlying mechanical properties. A cell-matrix composite was constructed by applying fibrin hydrogel with suspended cells onto decellularized sheets of human myocardium. We then implanted this composite onto the infarct bed in a nude rat model of cardiac infarction. We next characterized the myogenic and vasculogenic potential of immunoselected human MPCs and demonstrated that in vitro conditioning with a low concentration of TGF-β promoted an arteriogenic profile of gene expression. When implanted by composite scaffold, preconditioned MPCs greatly enhanced vascular network formation in the infarct bed by mechanisms involving the secretion of paracrine factors, such as SDF-1, and the migration of MPCs into ischemic myocardium, but not normal myocardium. Echocardiography demonstrated the recovery of baseline levels of left ventricular systolic dimensions and contractility when MPCs were delivered via composite scaffold. This adaptable platform could be readily extended to the delivery of other reparative cells of interest and used in quantitative studies of heart repair.