Histological validation of myocardial microstructure obtained from diffusion tensor magnetic resonance imaging

Three-dimensional cardiac computational modelling: methods, features and applications

The combination of computational models and biophysical simulations can help to interpret an array of experimental data and contribute to the understanding, diagnosis and treatment of complex diseases such as cardiac arrhythmias. For this reason, three-dimensional (3D) cardiac computational modelling is currently a rising field of research. The advance of medical imaging technology over the last decades has allowed the evolution from generic to patient-specific 3D cardiac models that faithfully represent the anatomy and different cardiac features of a given alive subject. Here we analyse sixty representative 3D cardiac computational models developed and published during the last fifty years, describing their information sources, features, development methods and online availability. This paper also reviews the necessary components to build a 3D computational model of the heart aimed at biophysical simulation, paying especial attention to cardiac electrophysiology (EP), and the existing approaches to incorporate those components. We assess the challenges associated to the different steps of the building process, from the processing of raw clinical or biological data to the final application, including image segmentation, inclusion of substructures and meshing among others. We briefly outline the personalisation approaches that are currently available in 3D cardiac computational modelling. Finally, we present examples of several specific applications, mainly related to cardiac EP simulation and model-based image analysis, showing the potential usefulness of 3D cardiac computational modelling into clinical environments as a tool to aid in the prevention, diagnosis and treatment of cardiac diseases.

Quantification of Shear Deformations and Corresponding Stresses in the Biaxially Tested Human Myocardium

One goal of cardiac research is to perform numerical simulations to describe/reproduce the mechanoelectrical function of the human myocardium in health and disease. Such simulations are based on a complex combination of mathematical models describing the passive mechanical behavior of the myocardium and its electrophysiology, i.e., the activation of cardiac muscle cells. The problem in developing adequate constitutive models is the shortage of experimental data suitable for detailed parameter estimation in specific functional forms. A combination of shear and biaxial extension tests with different loading protocols on different specimen orientations is necessary to capture adequately the direction-dependent (orthotropic) response of the myocardium. In most experimental animal studies, where planar biaxial extension tests on the myocardium have been conducted, the generated shear stresses were neither considered nor discussed. Hence, in this study a method is presented which allows the quantification of shear deformations and related stresses. It demonstrates an approach for experimenters as to how the generation of these shear stresses can be minimized during mechanical testing. Experimental results on 14 passive human myocardial specimens, obtained from nine human hearts, show the efficiency of this newly developed method. Moreover, the influence of the clamping technique of the specimen, i.e., the load transmission between the testing device and the tissue, on the stress response is determined by testing an isotropic material (Latex). We identified that the force transmission between the testing device and the specimen by means of hooks and cords does not influence the performed experiments. We further showed that in-plane shear stresses definitely exist in biaxially tested human ventricular myocardium, but can be reduced to a minimum by preparing the specimens in an appropriate manner. Moreover, we showed whether shear stresses can be neglected when performing planar biaxial extension tests on fiber-reinforced materials. The used method appears to be robust to quantify normal and shear deformations and corresponding stresses in biaxially tested human myocardium. This method can be applied for the mechanical characterization of any fiber-reinforced material using planar biaxial extension tests.

Numerical evaluation of myofiber orientation and transmural contractile strength on left ventricular function

The left ventricle (LV) of the heart is composed of a complex organization of cardiac muscle fibers, which contract to generate force and pump blood into the body. It has been shown that both the orientation and contractile strength of these myofibers vary across the ventricular wall. The hypothesis of the current study is that the transmural distributions of myofiber orientation and contractile strength interdependently impact LV pump function. In order to quantify these interactions a finite element (FE) model of the LV was generated, which incorporated transmural variations. The influences of myofiber orientation and contractile strength on the Starling relationship and the end-systolic (ES) apex twist of the LV were assessed. The results suggest that reductions in contractile strength within a specific transmural layer amplified the effects of altered myofiber orientation in the same layer, causing greater changes in stroke volume (SV). Furthermore, when the epicardial myofibers contracted the strongest, the twist of the LV apex was greatest, regardless of myofiber orientation. These results demonstrate the important role of transmural distribution of myocardial contractile strength and its interplay with myofiber orientation. The coupling between these two physiologic parameters could play a critical role in the progression of heart failure.

Simulation Methods and Validation Criteria for Modeling Cardiac Ventricular Electrophysiology

We describe a sequence of methods to produce a partial differential equation model of the electrical activation of the ventricles. In our framework, we incorporate the anatomy and cardiac microstructure obtained from magnetic resonance imaging and diffusion tensor imaging of a New Zealand White rabbit, the Purkinje structure and the Purkinje-muscle junctions, and an electrophysiologically accurate model of the ventricular myocytes and tissue, which includes transmural and apex-to-base gradients of action potential characteristics. We solve the electrophysiology governing equations using the finite element method and compute both a 6-lead precordial electrocardiogram (ECG) and the activation wavefronts over time. We are particularly concerned with the validation of the various methods used in our model and, in this regard, propose a series of validation criteria that we consider essential. These include producing a physiologically accurate ECG, a correct ventricular activation sequence, and the inducibility of ventricular fibrillation. Among other components, we conclude that a Purkinje geometry with a high density of Purkinje muscle junctions covering the right and left ventricular endocardial surfaces as well as transmural and apex-to-base gradients in action potential characteristics are necessary to produce ECGs and time activation plots that agree with physiological observations.

Probing dynamic myocardial microstructure with cardiac magnetic resonance diffusion tensor imaging

This article is an invited editorial comment on the paper entitled "In vivo cardiovascular magnetic resonance diffusion tensor imaging shows evidence of abnormal myocardial laminar orientations and mobility in hypertrophic cardiomyopathy" by Ferreira et al., and published as Journal of Cardiovascular Magnetic Resonance 2014; 16:87.

In vivo cardiovascular magnetic resonance diffusion tensor imaging shows evidence of abnormal myocardial laminar orientations and mobility in hypertrophic cardiomyopathy

BACKGROUND

Cardiac diffusion tensor imaging (cDTI) measures the magnitudes and directions of intramyocardial water diffusion. Assuming the cross-myocyte components to be constrained by the laminar microstructures of myocardium, we hypothesized that cDTI at two cardiac phases might identify any abnormalities of laminar orientation and mobility in hypertrophic cardiomyopathy (HCM).

METHODS

We performed cDTI in vivo at 3 Tesla at end-systole and late diastole in 11 healthy controls and 11 patients with HCM, as well as late gadolinium enhancement (LGE) for detection of regional fibrosis.

RESULTS

Voxel-wise analysis of diffusion tensors relative to left ventricular coordinates showed expected transmural changes of myocardial helix-angle, with no significant differences between phases or between HCM and control groups. In controls, the angle of the second eigenvector of diffusion (E2A) relative to the local wall tangent plane was larger in systole than diastole, in accord with previously reported changes of laminar orientation. HCM hearts showed higher than normal global E2A in systole (63.9° vs 56.4° controls, p=0.026) and markedly raised E2A in diastole (46.8° vs 24.0° controls, p<0.001). In hypertrophic regions, E2A retained a high, systole-like angulation even in diastole, independent of LGE, while regions of normal wall thickness did not (LGE present 57.8°, p=0.0028, LGE absent 54.8°, p=0.0022 vs normal thickness 38.1°).

CONCLUSIONS

In healthy controls, the angles of cross-myocyte components of diffusion were consistent with previously reported transmural orientations of laminar microstructures and their changes with contraction. In HCM, especially in hypertrophic regions, they were consistent with hypercontraction in systole and failure of relaxation in diastole. Further investigation of this finding is required as previously postulated effects of strain might be a confounding factor.

Characterization of diffuse fibrosis in the failing human heart via diffusion tensor imaging and quantitative histological validation

Non-invasive imaging techniques are highly desirable as an alternative to conventional biopsy for the characterization of the remodeling of tissues associated with disease progression, including end-stage heart failure. Cardiac diffusion tensor imaging (DTI) has become an established method for the characterization of myocardial microstructure. However, the relationships between diffuse myocardial fibrosis, which is a key biomarker for staging and treatment planning of the failing heart, and measured DTI parameters have yet to be investigated systematically. In this study, DTI was performed on left ventricular specimens collected from patients with chronic end-stage heart failure as a result of idiopathic dilated cardiomyopathy (n = 14) and from normal donors (n = 5). Scalar DTI parameters, including fractional anisotropy (FA) and mean (MD), primary (D1 ), secondary (D2 ) and tertiary (D3 ) diffusivities, were correlated with collagen content measured by digital microscopy. Compared with hearts from normal subjects, the FA in failing hearts decreased by 22%, whereas the MD, D2 and D3 increased by 12%, 14% and 24%, respectively (P < 0.01). No significant change was detected for D1 between the two groups. Furthermore, significant correlation was observed between the DTI scalar indices and quantitative histological measurements of collagen (i.e. fibrosis). Pearson's correlation coefficients (r) between collagen content and FA, MD, D2 and D3 were -0.51, 0.59, 0.56 and 0.62 (P < 0.05), respectively. The correlation between D1 and collagen content was not significant (r = 0.46, P = 0.05). Computational modeling analysis indicated that the behaviors of the DTI parameters as a function of the degree of fibrosis were well explained by compartmental exchange between myocardial and collagenous tissues. Combined, these findings suggest that scalar DTI parameters can be used as metrics for the non-invasive assessment of diffuse fibrosis in failing hearts.

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.

Dual-phase cardiac diffusion tensor imaging with strain correction

Purpose

In this work we present a dual-phase diffusion tensor imaging (DTI) technique that incorporates a correction scheme for the cardiac material strain, based on 3D myocardial tagging.

Methods

In vivo dual-phase cardiac DTI with a stimulated echo approach and 3D tagging was performed in 10 healthy volunteers. The time course of material strain was estimated from the tagging data and used to correct for strain effects in the diffusion weighted acquisition. Mean diffusivity, fractional anisotropy, helix, transverse and sheet angles were calculated and compared between systole and diastole, with and without strain correction. Data acquired at the systolic sweet spot, where the effects of strain are eliminated, served as a reference.

Results

The impact of strain correction on helix angle was small. However, large differences were observed in the transverse and sheet angle values, with and without strain correction. The standard deviation of systolic transverse angles was significantly reduced from 35.9±3.9° to 27.8°±3.5° (p<0.001) upon strain-correction indicating more coherent fiber tracks after correction. Myocyte aggregate structure was aligned more longitudinally in systole compared to diastole as reflected by an increased transmural range of helix angles (71.8°±3.9° systole vs. 55.6°±5.6°, p<0.001 diastole). While diastolic sheet angle histograms had dominant counts at high sheet angle values, systolic histograms showed lower sheet angle values indicating a reorientation of myocyte sheets during contraction.

Conclusion

An approach for dual-phase cardiac DTI with correction for material strain has been successfully implemented. This technique allows assessing dynamic changes in myofiber architecture between systole and diastole, and emphasizes the need for strain correction when sheet architecture in the heart is imaged with a stimulated echo approach.

Images as drivers of progress in cardiac computational modelling

Computational models have become a fundamental tool in cardiac research. Models are evolving to cover multiple scales and physical mechanisms. They are moving towards mechanistic descriptions of personalised structure and function, including effects of natural variability. These developments are underpinned to a large extent by advances in imaging technologies. This article reviews how novel imaging technologies, or the innovative use and extension of established ones, integrate with computational models and drive novel insights into cardiac biophysics. In terms of structural characterization, we discuss how imaging is allowing a wide range of scales to be considered, from cellular levels to whole organs. We analyse how the evolution from structural to functional imaging is opening new avenues for computational models, and in this respect we review methods for measurement of electrical activity, mechanics and flow. Finally, we consider ways in which combined imaging and modelling research is likely to continue advancing cardiac research, and identify some of the main challenges that remain to be solved.

Interrogation of living myocardium in multiple static deformation states with diffusion tensor and diffusion spectrum imaging

Diffusion tensor magnetic resonance imaging (MRI) reveals valuable insights into tissue histo-anatomy and microstructure, and has steadily gained traction in the cardiac community. Its wider use in small animal cardiac imaging in vivo has been constrained by its extreme sensitivity to motion, exaggerated by the high heart rates usually seen in rodents. Imaging of the isolated heart eliminates respiratory motion and, if conducted on arrested hearts, cardiac pulsation. This serves as an important intermediate step for basic and translational studies. However, investigating the micro-structural basis of cardiac deformation in the same heart requires observations in different deformation states.

Here, we illustrate the imaging of isolated rat hearts in three mechanical states mimicking diastole (cardioplegic arrest), left-ventricular (LV) volume overload (cardioplegic arrest plus LV balloon inflation), and peak systole (lithium-induced contracture). An optimised MRI-compatible Langendorff perfusion setup with the radio-frequency (RF) coil integrated into the wet chamber was developed for use in a 9.4T horizontal bore scanner. Signal-to-noise ratio improved significantly, by 75% compared to a previous design with external RF coil, and stability tests showed no significant changes in mean T₁, T₂ or LV wall thickness over a 170 min period. In contracture, we observed a significant reduction in mean fractional anisotropy from 0.32 ± 0.02 to 0.28 ± 0.02, as well as a significant rightward shift in helix angles with a decrease in the proportion of left-handed fibres, as referring to the locally prevailing cell orientation in the heart, from 24.9% to 23.3%, and an increase in the proportion of right-handed fibres from 25.5% to 28.4%. LV overload, in contrast, gave rise to a decrease in the proportion of left-handed fibres from 24.9% to 21.4% and an increase in the proportion of right-handed fibres from 25.5% to 26.0%. The modified perfusion and coil setup offers better performance and control over cardiac contraction states.

We subsequently performed high-resolution diffusion spectrum imaging (DSI) and 3D whole heart fibre tracking in fixed ex vivo rat hearts in slack state and contracture. As a model-free method, DSI augmented the measurements of water diffusion by also informing on multiple intra-voxel diffusion orientations and non-Gaussian diffusion. This enabled us to identify the transition from right- to left-handed fibres from the subendocardium to the subepicardium, as well as voxels in apical regions that were traversed by multiple fibres. We observed that both the mean generalised fractional anisotropy and mean kurtosis were lower in hearts in contracture compared to the slack state, by 23% and 9.3%, respectively. While its heavy acquisition burden currently limits the application of DSI in vivo, ongoing work in acceleration techniques may enable its use in live animals and patients. This would provide access to the as yet unexplored dimension of non-Gaussian diffusion that could serve as a highly sensitive marker of cardiac micro-structural integrity.

Microstructural impact of ischemia and bone marrow-derived cell therapy revealed with diffusion tensor magnetic resonance imaging tractography of the heart in vivo

BACKGROUND

The arrangement of myofibers in the heart is highly complex and must be replicated by injected cells to produce functional myocardium. A novel approach to characterize the microstructural response of the myocardium to ischemia and cell therapy, with the use of serial diffusion tensor magnetic resonance imaging tractography of the heart in vivo, is presented.

METHODS AND RESULTS

Validation of the approach was performed in normal (n=6) and infarcted mice (n=6) as well as healthy human volunteers. Mice (n=12) were then injected with bone marrow mononuclear cells 3 weeks after coronary ligation. In half of the mice the donor and recipient strains were identical, and in half the strains were different. A positive response to cell injection was defined by a decrease in mean diffusivity, an increase in fractional anisotropy, and the appearance of new myofiber tracts with the correct orientation. A positive response to bone marrow mononuclear cell injection was seen in 1 mouse. The response of the majority of mice to bone marrow mononuclear cell injection was neutral (9/12) or negative (2/12). The in vivo tractography findings were confirmed with histology.

CONCLUSIONS

Diffusion tensor magnetic resonance imaging tractography was able to directly resolve the ability of injected cells to generate new myofiber tracts and provided a fundamental readout of their regenerative capacity. A highly novel and translatable approach to assess the efficacy of cell therapy in the heart is thus presented.

Quantitative study of the effect of tissue microstructure on contraction in a computational model of rat left ventricle

Tissue microstructure, in particular the alignment of myocytes (fibre direction) and their lateral organisation into sheets, is fundamental to cardiac function. We studied the effect of microstructure on contraction in a computational model of rat left ventricular electromechanics. Different fibre models, globally rule-based or locally optimised to DT-MRI data, were compared, in order to understand whether a subject-specific fibre model would enhance the predictive power of our model with respect to the global ones. We also studied the impact of sheets on ventricular deformation by comparing: (a) a transversely isotropic versus an orthotropic material law and (b) a linear model with a bimodal model of sheet transmural variation. We estimated ejection fraction, wall thickening and base-to-apex shortening and compared them with measures from cine-MRI. We also evaluated Lagrangian strains as local metrics of cardiac deformation. Our results show that the subject-specific fibre model provides little improvement in the metric predictions with respect to global fibre models while material orthotropy allows closer agreement with measures than transverse isotropy. Nonetheless, the impact of sheets in our model is smaller than that of fibres. We conclude that further investigation of the modelling of sheet dynamics is necessary to fully understand the impact of tissue structure on cardiac deformation.

The Structural Basis of Functional Improvement in Response to Human Umbilical Cord Blood Stem Cell Transplantation in Hearts With Postinfarct LV Remodeling

Cellular therapy for myocardial repair has been one of the most intensely investigated interventional strategies for acute myocardial infarction. Although the therapeutic potential of stem cells has been demonstrated in various studies, the underlying mechanisms for such improvements are poorly understood. In the present study, we investigated the long-term effects of stem cell therapy on both myocardial fiber organization and regional contractile function using a rat model of postinfarct remodeling. Human nonhematopoietic umbilical cord blood stem cells (nh-UCBSCs) were administered via tail vein to rats 2 days after infarct surgery. Animals were maintained without immunosuppressive therapy. In vivo and ex vivo MR imaging was performed on infarct hearts 10 months after cell transplantation. Compared to the age-matched rats exposed to the identical surgery, both global and regional cardiac functions of the nh-UCBSC-treated hearts, such as ejection fraction, ventricular strain, and torsion, were significantly improved. More importantly, the treated hearts exhibited preserved fiber orientation and water diffusivities that were similar to those in sham-operated control hearts. These data provide the first evidence that nh-UCBSC treatment may prevent/delay untoward structural remodeling in postinfarct hearts, which supports the improved LV function observed in vivo in the absence of immunosuppression, suggesting a beneficial paracrine effect occurred with the cellular therapy.

Diffusion MRI tractography of the developing human fetal heart

Objective

Human myocardium has a complex and anisotropic 3D fiber pattern. It remains unknown, however, when in fetal life this anisotropic pattern develops and whether the human heart is structurally fully mature at birth. We aimed here to use diffusion tensor MRI (DTI) tractography to characterize the evolution of fiber architecture in the developing human fetal heart.

Methods

Human fetal hearts (n = 5) between 10–19 weeks of gestation were studied. The heart from a 6-day old neonate and an adult human heart served as controls. The degree of myocardial anisotropy was measured by calculating the fractional anisotropy (FA) index. In addition, fiber tracts were created by numerically integrating the primary eigenvector field in the heart into coherent streamlines.

Results

At 10–14 weeks the fetal hearts were highly isotropic and few tracts could be resolved. Between 14–19 weeks the anisotropy seen in the adult heart began to develop. Coherent fiber tracts were well resolved by 19 weeks. The 19-week myocardium, however, remained weakly anisotropic with a low FA and no discernable sheet structure.

Conclusions

The human fetal heart remains highly isotropic until 14–19 weeks, at which time cardiomyocytes self-align into coherent tracts. This process lags 2–3 months behind the onset of cardiac contraction, which may be a prerequisite for cardiomyocyte maturation and alignment. No evidence of a connective tissue scaffold guiding this process could be identified by DTI. Maturation of the heart’s sheet structure occurs late in gestation and evolves further after birth.

Cardiovascular magnetic resonance: deeper insights through bioengineering

Heart disease is the main cause of morbidity and mortality worldwide, with coronary artery disease, diabetes, and obesity being major contributing factors. Cardiovascular magnetic resonance (CMR) can provide a wealth of quantitative information on the performance of the heart, without risk to the patient. Quantitative analyses of these data can substantially augment the diagnostic quality of CMR examinations and can lead to more effective characterization of disease and quantification of treatment benefit. This review provides an overview of the current state of the art in CMR with particular regard to the quantification of motion, both microscopic and macroscopic, and the application of bioengineering analysis for the evaluation of cardiac mechanics. We discuss the current clinical practice and the likely advances in the next 5-10 years, as well as the ways in which clinical examinations can be augmented by bioengineering analysis of strain, compliance, and stress.

Patient-specific modeling of the heart: estimation of ventricular fiber orientations

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

Cardiomyocyte architectural plasticity in fetal, neonatal, and adult pig hearts delineated with diffusion tensor MRI

Cardiomyocyte organization is a critical determinant of coordinated cardiac contractile function. Because of the acute opening of the pulmonary circulation, the relative workload of the left ventricle (LV) and right ventricle (RV) changes substantially immediately after birth. We hypothesized that three-dimensional cardiomyocyte architecture might be required to adapt rapidly to accommodate programmed perinatal changes of cardiac function. Isolated fixed hearts from pig fetuses or pigs at midgestation, preborn, postnatal day 1 (P1), postnatal day 5, postnatal day 14 (P14), and adulthood (n = 5 for each group) were acquired for diffusion-weighted magnetic resonance imaging. Cardiomyocyte architecture was visualized by three-dimensional fiber tracking and was quantitatively evaluated by the measured helix angle (α(h)). Upon the completion of MRI, hearts were sectioned and stained with hematoxylin/eosin (H&E) to evaluate cardiomyocyte alignment, with picrosirius red to evaluate collagen content, and with anti-Ki67 to evaluate postnatal cell proliferation. The helical architecture of cardiomyocyte was observed as early as the midgestational period. Postnatal changes of cardiomyocyte architecture were observed from P1 to P14, which primary occurred in the septum and RV free wall (RVFW). In the septum, the volume ratio of LV- vs. RV-associated cardiomyocytes rapidly changed from RV-LV balanced pattern at birth to LV dominant pattern by P14. In the RVFW, subendocardial α(h) decreased by ~30° from P1 to P14. These findings indicate that the helical architecture of cardiomyocyte is developed as early as the midgestation period. Substantial and rapid adaptive changes in cardiac microarchitecture suggested considerable developmental plasticity of cardiomyocyte form and function in the postnatal period in response to altered cardiac mechanical function.

Fiber architecture in remodeled myocardium revealed with a quantitative diffusion CMR tractography framework and histological validation

BACKGROUND

The study of myofiber reorganization in the remote zone after myocardial infarction has been performed in 2D. Microstructural reorganization in remodeled hearts, however, can only be fully appreciated by considering myofibers as continuous 3D entities. The aim of this study was therefore to develop a technique for quantitative 3D diffusion CMR tractography of the heart, and to apply this method to quantify fiber architecture in the remote zone of remodeled hearts.

METHODS

Diffusion Tensor CMR of normal human, sheep, and rat hearts, as well as infarcted sheep hearts was performed ex vivo. Fiber tracts were generated with a fourth-order Runge-Kutta integration technique and classified statistically by the median, mean, maximum, or minimum helix angle (HA) along the tract. An index of tract coherence was derived from the relationship between these HA statistics. Histological validation was performed using phase-contrast microscopy.

RESULTS

In normal hearts, the subendocardial and subepicardial myofibers had a positive and negative HA, respectively, forming a symmetric distribution around the midmyocardium. However, in the remote zone of the infarcted hearts, a significant positive shift in HA was observed. The ratio between negative and positive HA variance was reduced from 0.96 ± 0.16 in normal hearts to 0.22 ± 0.08 in the remote zone of the remodeled hearts (p < 0.05). This was confirmed histologically by the reduction of HA in the subepicardium from -52.03° ± 2.94° in normal hearts to -37.48° ± 4.05° in the remote zone of the remodeled hearts (p < 0.05).

CONCLUSIONS

A significant reorganization of the 3D fiber continuum is observed in the remote zone of remodeled hearts. The positive (rightward) shift in HA in the remote zone is greatest in the subepicardium, but involves all layers of the myocardium. Tractography-based quantification, performed here for the first time in remodeled hearts, may provide a framework for assessing regional changes in the left ventricle following infarction.

In vivo diffusion tensor MRI of the human heart: reproducibility of breath-hold and navigator-based approaches

The aim of this study was to implement a quantitative in vivo cardiac diffusion tensor imaging (DTI) technique that was robust, reproducible, and feasible to perform in patients with cardiovascular disease. A stimulated-echo single-shot echo-planar imaging (EPI) sequence with zonal excitation and parallel imaging was implemented, together with a novel modification of the prospective navigator (NAV) technique combined with a biofeedback mechanism. Ten volunteers were scanned on two different days, each time with both multiple breath-hold (MBH) and NAV multislice protocols. Fractional anisotropy (FA), mean diffusivity (MD), and helix angle (HA) fiber maps were created. Comparison of initial and repeat scans showed good reproducibility for both MBH and NAV techniques for FA (P > 0.22), MD (P > 0.15), and HA (P > 0.28). Comparison of MBH and NAV FA (FAMBHday1 = 0.60 ± 0.04, FANAVday1 = 0.60 ± 0.03, P = 0.57) and MD (MDMBHday1 = 0.8 ± 0.2 × 10(-3) mm(2) /s, MDNAVday1 = 0.9 ± 0.2 × 10(-3) mm(2) /s, P = 0.07) values showed no significant differences, while HA values (HAMBHday1Endo = 22 ± 10°, HAMBHday1Mid-Endo = 20 ± 6°, HAMBHday1Mid-Epi = -1 ± 6°, HAMBHday1Epi = -17 ± 6°, HANAVday1Endo = 7 ± 7°, HANAVday1Mid-Endo = 13 ± 8°, HANAVday1Mid-Epi = -2 ± 7°, HANAVday1Epi = -14 ± 6°) were significantly different. The scan duration was 20% longer with the NAV approach. Currently, the MBH approach is the more robust in normal volunteers. While the NAV technique still requires resolution of some bulk motion sensitivity issues, these preliminary experiments show its potential for in vivo clinical cardiac diffusion tensor imaging and for delivering high-resolution in vivo 3D DTI tractography of the heart.

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.

Focal but reversible diastolic sheet dysfunction reflects regional calcium mishandling in dystrophic mdx mouse hearts

Cardiac dysfunction is a primary cause of patient mortality in Duchenne muscular dystrophy, potentially related to elevated cytosolic calcium. However, the regional versus global functional consequences of cellular calcium mishandling have not been defined in the whole heart. Here we sought for the first time to elucidate potential regional dependencies between calcium mishandling and myocardial fiber/sheet function as a manifestation of dystrophin-deficient (mdx) cardiomyopathy. Isolated-perfused hearts from 16-mo-old mdx (N = 10) and wild-type (WT; N = 10) were arrested sequentially in diastole and systole for diffusion tensor MRI quantification of myocardial sheet architecture and function. When compared with WT hearts, mdx hearts exhibited normal systolic sheet architecture but a lower diastolic sheet angle magnitude (|β|) in the basal region. The regional diastolic sheet dysfunction was normalized by reducing perfusate calcium concentrations. Optical mapping of calcium transients in isolated hearts (3 mdx and 4 WT) revealed a stretch-inducible regional defect of intracellular calcium reuptake, reflected by a 25% increase of decay times (T(50)) and decay constants, at the base of mdx hearts. The basal region of mdx hearts also exhibited greater fibrosis than did the apex, which matched the regional sheet dysfunction. We conclude that myocardial diastolic sheet dysfunction is observed initially in basal segments along with calcium mishandling, ultimately culminating in increased fibrosis. The preservation of relatively normal calcium reuptake and diastolic/systolic sheet mechanics throughout the rest of the heart, together with the rapid reversibility of functional defects by reducing cytosolic calcium, points to the significance of regional mechanical factors in the progression of the disease.

Cardiac arrhythmia mechanisms in rats with heart failure induced by pulmonary hypertension

Pulmonary hypertension provokes right heart failure and arrhythmias. Better understanding of the mechanisms underlying these arrhythmias is needed to facilitate new therapeutic approaches for the hypertensive, failing right ventricle (RV). The aim of our study was to identify the mechanisms generating arrhythmias in a model of RV failure induced by pulmonary hypertension. Rats were injected with monocrotaline to induce either RV hypertrophy or failure or with saline (control). ECGs were measured in conscious, unrestrained animals by telemetry. In isolated hearts, electrical activity was measured by optical mapping and myofiber orientation by diffusion tensor-MRI. Sarcoplasmic reticular Ca(2+) handling was studied in single myocytes. Compared with control animals, the T-wave of the ECG was prolonged and in three of seven heart failure animals, prominent T-wave alternans occurred. Discordant action potential (AP) alternans occurred in isolated failing hearts and Ca(2+) transient alternans in failing myocytes. In failing hearts, AP duration and dispersion were increased; conduction velocity and AP restitution were steeper. The latter was intrinsic to failing single myocytes. Failing hearts had greater fiber angle disarray; this correlated with AP duration. Failing myocytes had reduced sarco(endo)plasmic reticular Ca(2+)-ATPase activity, increased sarcoplasmic reticular Ca(2+)-release fraction, and increased Ca(2+) spark leak. In hypertrophied hearts and myocytes, dysfunctional adaptation had begun, but alternans did not develop. We conclude that increased electrical and structural heterogeneity and dysfunctional sarcoplasmic reticular Ca(2+) handling increased the probability of alternans, a proarrhythmic predictor of sudden cardiac death. These mechanisms are potential therapeutic targets for the correction of arrhythmias in hypertensive, failing RVs.

Label-free photoacoustic microscopy of myocardial sheet architecture

Cardiac myofibers are organized into sheet architectures, which contribute to up to 40% of the heart wall thickening for ejection of blood for circulation. It is important to delineate the sheet architecture for a better understanding of cardiac mechanisms. However, current sheet imaging technologies are limited by fixation-induced dehydration/deformation and low spatial resolution. Here we implemented high-resolution label-free photoacoustic microscopy (PAM) of the myocardial sheet architecture. With high endogenous optical-absorption contrast originating mainly from cytochrome, myoglobin, and melanin, PAM can image the unfixed, unstained and unsliced heart without introducing deformation artifacts. A fresh blood-free mouse heart was imaged by PAM ex vivo. The three-dimensional branching sheets were clearly identified within 150 [micro sign]m depth. Various morphological parameters were derived from the PAM image. The sheet thickness (80 ± 10 μm) and the cleavage height (11 ± 1 μm) were derived from an undehydrated heart for the first time. Therefore, PAM has the potential for the functional imaging of sheet architecture in ex vivo perfused and viable hearts.

A novel rule-based algorithm for assigning myocardial fiber orientation to computational heart models

Electrical waves traveling throughout the myocardium elicit muscle contractions responsible for pumping blood throughout the body. The shape and direction of these waves depend on the spatial arrangement of ventricular myocytes, termed fiber orientation. In computational studies simulating electrical wave propagation or mechanical contraction in the heart, accurately representing fiber orientation is critical so that model predictions corroborate with experimental data. Typically, fiber orientation is assigned to heart models based on Diffusion Tensor Imaging (DTI) data, yet few alternative methodologies exist if DTI data is noisy or absent. Here we present a novel Laplace-Dirichlet Rule-Based (LDRB) algorithm to perform this task with speed, precision, and high usability. We demonstrate the application of the LDRB algorithm in an image-based computational model of the canine ventricles. Simulations of electrical activation in this model are compared to those in the same geometrical model but with DTI-derived fiber orientation. The results demonstrate that activation patterns from simulations with LDRB and DTI-derived fiber orientations are nearly indistinguishable, with relative differences ≤6%, absolute mean differences in activation times ≤3.15 ms, and positive correlations ≥0.99. These results convincingly show that the LDRB algorithm is a robust alternative to DTI for assigning fiber orientation to computational heart models.

Three-dimensional mechanisms of increased vulnerability to electric shocks in myocardial infarction: altered virtual electrode polarizations and conduction delay in the peri-infarct zone

Defibrillation efficacy is decreased in infarcted hearts, but the mechanisms by which infarcted hearts are more vulnerable to electric shocks than healthy hearts remain poorly understood. The goal of this study was to provide insight into the 3D mechanisms for the increased vulnerability to electric shocks in infarcted hearts. We hypothesized that changes in virtual electrode polarizations (VEPs) and propagation delay through the peri-infarct zone (PZ) were responsible. We developed a micro anatomically detailed rabbit ventricular model with chronic myocardial infarction from magnetic resonance imaging and enriched the model with data from optical mapping experiments. We further developed a control model without the infarct. The simulation protocol involved apical pacing followed by biphasic shocks. Simulation results from both models were compared.The upper limit of vulnerability(ULV) was 8 V cm(-1) in the infarction model and 4 V cm(-1) in the control model. VEPs were less pronounced in the infarction model, providing a larger excitable area for postshock propagation but smaller transmembrane potential gradients to initiate new wavefronts. Initial post-shock transmural activation occurred at a later time in the infarction model, and the PZ served to delay propagation in subsequent beats. The presence of the PZ was found to be responsible for the increased vulnerability.

Image-based estimation of ventricular fiber orientations for personalized modeling of cardiac electrophysiology

Technological limitations pose a major challenge to acquisition of myocardial fiber orientations for patient-specific modeling of cardiac (dys)function and assessment of therapy. The objective of this project was to develop a methodology to estimate cardiac fiber orientations from in vivo images of patient heart geometries. An accurate representation of ventricular geometry and fiber orientations was reconstructed, respectively, from high-resolution ex vivo structural magnetic resonance (MR) and diffusion tensor (DT) MR images of a normal human heart, referred to as the atlas. Ventricular geometry of a patient heart was extracted, via semiautomatic segmentation, from an in vivo computed tomography (CT) image. Using image transformation algorithms, the atlas ventricular geometry was deformed to match that of the patient. Finally, the deformation field was applied to the atlas fiber orientations to obtain an estimate of patient fiber orientations. The accuracy of the fiber estimates was assessed using six normal and three failing canine hearts. The mean absolute difference between inclination angles of acquired and estimated fiber orientations was 15.4°. Computational simulations of ventricular activation maps and pseudo-ECGs in sinus rhythm and ventricular tachycardia indicated that there are no significant differences between estimated and acquired fiber orientations at a clinically observable level.

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.

Electromechanical and structural alterations in the aging rabbit heart and aorta

Aging increases the risk for arrhythmias and sudden cardiac death (SCD). We aimed at elucidating aging-related electrical, functional, and structural changes in the heart and vasculature that account for this heightened arrhythmogenic risk. Young (5-9 mo) and old (3.5-6 yr) female New Zealand White (NZW) rabbits were subjected to in vivo hemodynamic, electrophysiological, and echocardiographic studies as well as ex vivo optical mapping, high-field magnetic resonance imaging (MRI), and histochemical experiments. Aging increased aortic stiffness (baseline pulse wave velocity: young, 3.54 ± 0.36 vs. old, 4.35 ± 0.28 m/s, P < 0.002) and diastolic (end diastolic pressure-volume relations: 3.28 ± 0.5 vs. 4.95 ± 1.5 mmHg/ml, P < 0.05) and systolic (end systolic pressure-volume relations: 20.56 ± 4.2 vs. 33.14 ± 8.4 mmHg/ml, P < 0.01) myocardial elastances in old rabbits. Electrophysiological and optical mapping studies revealed age-related slowing of ventricular and His-Purkinje conduction (His-to-ventricle interval: 23 ± 2.5 vs. 31.9 ± 2.9 ms, P < 0.0001), altered conduction anisotropy, and a greater inducibility of ventricular fibrillation (VF, 3/12 vs. 7/9, P < 0.05) in old rabbits. Histochemical studies confirmed an aging-related increased fibrosis in the ventricles. MRI showed a deterioration of the free-running Purkinje fiber network in ventricular and septal walls in old hearts as well as aging-related alterations of the myofibrillar orientation and myocardial sheet structure that may account for this slowed conduction velocity. Aging leads to parallel stiffening of the aorta and the heart, including an increase in systolic stiffness and contractility and diastolic stiffness. Increasingly, anisotropic conduction velocity due to fibrosis and altered myofibrillar orientation and myocardial sheet structure may contribute to the pathogenesis of VF in old hearts. The aging rabbit model represents a useful tool for elucidating age-related changes that predispose the aging heart to arrhythmias and SCD.

Image-based estimation of ventricular fiber orientations for patient-specific simulations

Patient-specific simulation of heart (dys)function aimed at personalizing cardiac therapy are hampered by the absence of in vivo imaging technology for clinically acquiring myocardial fiber orientations. In this research, we develop a methodology to predict ventricular fiber orientations of a patient heart, given the geometry of the heart and an atlas. We test the methodology by comparing the estimated fiber orientations with measured ones, and by quantifying the effect of the estimation error on outcomes of electrophysiological simulations, in normal and failing canine hearts. The new insights obtained from the project will pave the way for the development of patient-specific models of the heart that can aid physicians in personalized diagnosis and decisions regarding electrophysiological interventions.

Automated determination of axonal orientation in the deep white matter of the human brain

The wide-spread utilization of diffusion-weighted imaging in the clinical neurosciences to assess white-matter (WM) integrity and architecture calls for robust validation strategies applied to the data that are acquired with noninvasive imaging. However, the pathology and detailed fiber architecture of WM tissue can only be observed postmortem. With these considerations in mind, we designed an automated method for the determination of axonal orientation in high-resolution microscope images. The algorithm was tested on tissue that was stained using a silver impregnation technique that was optimized to resolve axonal fibers against very low levels of background. The orientation of individual nerve fibers was detected using spatial filtering and a template-matching algorithm, and the results are displayed as color-coded overlays. Quantitative models of WM fiber architecture at the microscopic level can lead to improved interpretation of low-resolution neuroimaging data and to more accurate mapping of fiber pathways in the human brain.

Diffusion-Tensor MRI Based Skeletal Muscle Fiber Tracking

A skeletal muscle's function is strongly influenced by the internal organization and geometric properties of its fibers, a property known as muscle architecture. Diffusion-tensor magnetic resonance imaging-based fiber tracking provides a powerful tool for non-invasive muscle architecture studies, has three-dimensional sensitivity, and uses a fixed frame of reference. Significant advances have been made in muscle fiber tracking technology, including defining seed points for fiber tracking, quantitatively characterizing muscle architecture, implementing denoising procedures, and testing validity and repeatability. Some examples exist of how these data can be integrated with those from other advanced MRI and computational methods to provide novel insights into muscle function. Perspectives are offered regarding future directions in muscle diffusion-tensor imaging, including needs to develop an improved understanding for the microstructural basis for reduced and anisotropic diffusion, establish the best practices for data acquisition and analysis, and integrate fiber tracking with other physiological data.

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.

Diffusion tensor imaging: exploring the motor networks and clinical applications

With the advances in diffusion magnetic resonance (MR) imaging techniques, diffusion tensor imaging (DTI) has been applied to a number of neurological conditions because DTI can demonstrate microstructures of the brain that are not assessable with conventional MR imaging. Tractography based on DTI offers gross visualization of the white matter fiber architecture in the human brain in vivo. Degradation of restrictive barriers and disruption of the cytoarchitecture result in changes in the diffusion of water molecules in various pathological conditions, and these conditions can also be assessed with DTI. Yet many factors may influence the ability to apply DTI clinically, so these techniques have to be used with a cautious hand.

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.

Perception-based visualization of manifold-valued medical images using distance-preserving dimensionality reduction

A method for visualizing manifold-valued medical image data is proposed. The method operates on images in which each pixel is assumed to be sampled from an underlying manifold. For example, each pixel may contain a high dimensional vector, such as the time activity curve (TAC) in a dynamic positron emission tomography (dPET) or a dynamic single photon emission computed tomography (dSPECT) image, or the positive semi-definite tensor in a diffusion tensor magnetic resonance image (DTMRI). A nonlinear mapping reduces the dimensionality of the pixel data to achieve two goals: distance preservation and embedding into a perceptual color space. We use multidimensional scaling distance-preserving mapping to render similar pixels (e.g., DT or TAC pixels) with perceptually similar colors. The 3D CIELAB perceptual color space is adopted as the range of the distance preserving mapping, with a final similarity transform mapping colors to a maximum gamut size. Similarity between pixels is either determined analytically as geodesics on the manifold of pixels or is approximated using manifold learning techniques. In particular, dissimilarity between DTMRI pixels is evaluated via a Log-Euclidean Riemannian metric respecting the manifold of the rank 3, second-order positive semi-definite DTs, whereas the dissimilarity between TACs is approximated via ISOMAP. We demonstrate our approach via artificial high-dimensional, manifold-valued data, as well as case studies of normal and pathological clinical brain and heart DTMRI, dPET, and dSPECT images. Our results demonstrate the effectiveness of our approach in capturing, in a perceptually meaningful way, important features in the data.

Models of cardiac electromechanics based on individual hearts imaging data: image-based electromechanical models of the heart

Current multi-scale computational models of ventricular electromechanics describe the full process of cardiac contraction on both the micro- and macro- scales including: the depolarization of cardiac cells, the release of calcium from intracellular stores, tension generation by cardiac myofilaments, and mechanical contraction of the whole heart. Such models are used to reveal basic mechanisms of cardiac contraction as well as the mechanisms of cardiac dysfunction in disease conditions. In this paper, we present a methodology to construct finite element electromechanical models of ventricular contraction with anatomically accurate ventricular geometry based on magnetic resonance and diffusion tensor magnetic resonance imaging of the heart. The electromechanical model couples detailed representations of the cardiac cell membrane, cardiac myofilament dynamics, electrical impulse propagation, ventricular contraction, and circulation to simulate the electrical and mechanical activity of the ventricles. The utility of the model is demonstrated in an example simulation of contraction during sinus rhythm using a model of the normal canine ventricles.

Transmural and apicobasal gradients in repolarization contribute to T-wave genesis in human surface ECG

The cellular basis of the T-wave morphology of surface ECG remains controversial in clinical cardiology. We examined the effect of action potential duration (APD) distribution on T-wave morphology using a realistic model of the human ventricle and torso. We developed a finite-element model of the ventricle consisting of ∼26 million elements, including the conduction system, each implemented with the ion current model of cardiomyocytes. This model was embedded in a torso model with distinct organ structures to obtain the standard ECG leads. The APD distribution was changed in the transmural direction by locating the M cells in either the endocardial or epicardial region. We also introduced apicobasal gradients by modifying the ion channel parameters. Both the transmural gradient (with M cells on the endocardial side) and the apicobasal gradient produced positive T waves, although a very large gradient was required for the apicobasal gradient. By contrast, T waves obtained with the transmural gradient were highly symmetric and, therefore, did not represent the true physiological state. Only combination of the transmural and the moderate apicobasal gradients produced physiological T waves in surface ECG. Positive T waves in surface ECG mainly originated from the transmural distribution of APD with M cells on the endocardial side, although the apicobasal gradient was also required to attain the physiological waveform.

Combined diffusion and strain tensor MRI reveals a heterogeneous, planar pattern of strain development during isometric muscle contraction

The purposes of this study were to create a three-dimensional representation of strain during isometric contraction in vivo and to interpret it with respect to the muscle fiber direction. Diffusion tensor MRI was used to measure the muscle fiber direction of the tibialis anterior (TA) muscle of seven healthy volunteers. Spatial-tagging MRI was used to measure linear strains in six directions during separate 50% maximal isometric contractions of the TA. The strain tensor (E) was computed in the TA's deep and superficial compartments and compared with the respective diffusion tensors. Diagonalization of E revealed a planar strain pattern, with one nonzero negative strain (ε(N)) and one nonzero positive strain (ε(P)); both strains were larger in magnitude (P < 0.05) in the deep compartment [ε(N) = -40.4 ± 4.3%, ε(P) = 35.1 ± 3.5% (means ± SE)] than in the superficial compartment (ε(N) = -24.3 ± 3.9%, ε(P) = 6.3 ± 4.9%). The principal shortening direction deviated from the fiber direction by 24.0 ± 1.3° and 39.8 ± 6.1° in the deep and superficial compartments, respectively (P < 0.05, deep vs. superficial). The deviation of the shortening direction from the fiber direction was due primarily to the lower angle of elevation of the shortening direction over the axial plane than that of the fiber direction. It is concluded that three-dimensional analyses of strain interpreted with respect to the fiber architecture are necessary to characterize skeletal muscle contraction in vivo. The deviation of the principal shortening direction from the fiber direction may relate to intramuscle variations in fiber length and pennation angle.

Image-based models of cardiac structure in health and disease

Computational approaches to investigating the electromechanics of healthy and diseased hearts are becoming essential for the comprehensive understanding of cardiac function. In this article, we first present a brief review of existing image-based computational models of cardiac structure. We then provide a detailed explanation of a processing pipeline which we have recently developed for constructing realistic computational models of the heart from high resolution structural and diffusion tensor (DT) magnetic resonance (MR) images acquired ex vivo. The presentation of the pipeline incorporates a review of the methodologies that can be used to reconstruct models of cardiac structure. In this pipeline, the structural image is segmented to reconstruct the ventricles, normal myocardium, and infarct. A finite element mesh is generated from the segmented structural image, and fiber orientations are assigned to the elements based on DTMR data. The methods were applied to construct seven different models of healthy and diseased hearts. These models contain millions of elements, with spatial resolutions in the order of hundreds of microns, providing unprecedented detail in the representation of cardiac structure for simulation studies.

Human soleus muscle architecture at different ankle joint angles from magnetic resonance diffusion tensor imaging

The orientation of muscle fibers influences the physiological cross-sectional area, the relationship between fiber shortening and aponeurosis shear, and the total force produced by the muscle. Such architectural parameters are challenging to determine particularly in vivo in multicompartment structures such as the human soleus with a complex arrangement of muscle fibers. The objective of this study was to map the fiber architecture of the human soleus in vivo at rest in both neutral and plantarflexed ankle positions using an MRI-based method of diffusion tensor imaging (DTI). Six subjects were imaged at 3 Tesla with the foot at rest in the two ankle positions. Eigenvalues, fractional anisotropy (FA), and eigenvector orientations of fibers in the different soleus subcompartments were evaluated after denoising of the diffusion tensor. The fiber architecture from DTI was similar to earlier studies based on a 3D fiber model from cadavers. The three eigenvalues of the diffusion tensor increased by ∼14% on increasing the joint plantarflexion angle in all of the soleus subcompartments, whereas FA showed a trend to decrease in the posterior and marginal soleus and to increase in the anterior soleus. The angle change in the lead eigenvector between the two foot positions was significant: ∼41° for the posterior soleus and ∼48° for the anterior soleus. Fibers tracked from the subcompartments support these changes seen in the eigenvector orientations. DTI-derived, subject-specific, muscle morphological data could potentially be used to model a more complete description of muscle performance and changes from disease.

Construction and validation of anisotropic and orthotropic ventricular geometries for quantitative predictive cardiac electrophysiology

Reaction-diffusion computational models of cardiac electrophysiology require both dynamic excitation models that reconstruct the action potentials of myocytes as well as datasets of cardiac geometry and architecture that provide the electrical diffusion tensor D, which determines how excitation spreads through the tissue. We illustrate an experimental pipeline we have developed in our laboratories for constructing and validating such datasets. The tensor D changes with location in the myocardium, and is determined by tissue architecture. Diffusion tensor magnetic resonance imaging (DT-MRI) provides three eigenvectors e(i) and eigenvalues λ(i) at each voxel throughout the tissue that can be used to reconstruct this architecture. The primary eigenvector e(1) is a histologically validated measure of myocyte orientation (responsible for anisotropic propagation). The secondary and tertiary eigenvectors (e(2) and e(3)) specify the directions of any orthotropic structure if λ(2) is significantly greater than λ(3)-this orthotropy has been identified with sheets or cleavage planes. For simulations, the components of D are scaled in the fibre and cross-fibre directions for anisotropic simulations (or fibre, sheet and sheet normal directions for orthotropic tissues) so that simulated conduction velocities match values from optical imaging or plunge electrode experiments. The simulated pattern of propagation of action potentials in the models is partially validated by optical recordings of spatio-temporal activity on the surfaces of hearts. We also describe several techniques that enhance components of the pipeline, or that allow the pipeline to be applied to different areas of research: Q ball imaging provides evidence for multi-modal orientation distributions within a fraction of voxels, infarcts can be identified by changes in the anisotropic structure-irregularity in myocyte orientation and a decrease in fractional anisotropy, clinical imaging provides human ventricular geometry and can identify ischaemic and infarcted regions, and simulations in human geometries examine the roles of anisotropic and orthotropic architecture in the initiation of arrhythmias.

Cardiovascular magnetic resonance imaging in experimental models

Cardiovascular magnetic resonance (CMR) imaging is the modality of choice for clinical studies of the heart and vasculature, offering detailed images of both structure and function with high temporal resolution.

Small animals are increasingly used for genetic and translational research, in conjunction with models of common pathologies such as myocardial infarction. In all cases, effective methods for characterising a wide range of functional and anatomical parameters are crucial for robust studies.

CMR is the gold-standard for the non-invasive examination of these models, although physiological differences, such as rapid heart rate, make this a greater challenge than conventional clinical imaging. However, with the help of specialised magnetic resonance (MR) systems, novel gating strategies and optimised pulse sequences, high-quality images can be obtained in these animals despite their small size.

In this review, we provide an overview of the principal CMR techniques for small animals for example cine, angiography and perfusion imaging, which can provide measures such as ejection fraction, vessel anatomy and local blood flow, respectively. In combination with MR contrast agents, regional dysfunction in the heart can also be identified and assessed. We also discuss optimal methods for analysing CMR data, particularly the use of semi-automated tools for parameter measurement to reduce analysis time. Finally, we describe current and emerging methods for imaging the developing heart, aiding characterisation of congenital cardiovascular defects.

Advanced small animal CMR now offers an unparalleled range of cardiovascular assessments. Employing these methods should allow new insights into the structural, functional and molecular basis of the cardiovascular system.

Molecular and Microstructural Imaging of the Myocardium

The past year has witnessed ongoing progress in the field of molecular MRI of the myocardium. In addition, several novel fluorescent agents have been introduced and used to image remodeling in the injured myocardium. New techniques to image myocardial microstructure, such as diffusion spectrum MRI, have also been introduced and have tremendous potential for integration and synergy with molecular MRI. In the current review we focus on these and other advances in the field that have occurred over the past year.