Combined diffusion and strain MRI reveals structure and function of human myocardial laminar sheets in vivo

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.

Biomedical Imaging in Experimental Models of Cardiovascular Disease

Major advances in biomedical imaging have occurred over the last 2 decades and now allow many physiological, cellular, and molecular processes to be imaged noninvasively in small animal models of cardiovascular disease. Many of these techniques can be also used in humans, providing pathophysiological context and helping to define the clinical relevance of the model. Ultrasound remains the most widely used approach, and dedicated high-frequency systems can obtain extremely detailed images in mice. Likewise, dedicated small animal tomographic systems have been developed for magnetic resonance, positron emission tomography, fluorescence imaging, and computed tomography in mice. In this article, we review the use of ultrasound and positron emission tomography in small animal models, as well as emerging contrast mechanisms in magnetic resonance such as diffusion tensor imaging, hyperpolarized magnetic resonance, chemical exchange saturation transfer imaging, magnetic resonance elastography and strain, arterial spin labeling, and molecular imaging.

Estimating cardiomyofiber strain in vivo by solving a computational model

Since heart contraction results from the electrically activated contraction of millions of cardiomyocytes, a measure of cardiomyocyte shortening mechanistically underlies cardiac contraction. In this work we aim to measure preferential aggregate cardiomyocyte ("myofiber") strains based on Magnetic Resonance Imaging (MRI) data acquired to measure both voxel-wise displacements through systole and myofiber orientation. In order to reduce the effect of experimental noise on the computed myofiber strains, we recast the strains calculation as the solution of a boundary value problem (BVP). This approach does not require a calibrated material model, and consequently is independent of specific myocardial material properties. The solution to this auxiliary BVP is the displacement field corresponding to assigned values of myofiber strains. The actual myofiber strains are then determined by minimizing the difference between computed and measured displacements. The approach is validated using an analytical phantom, for which the ground-truth solution is known. The method is applied to compute myofiber strains using in vivo displacement and myofiber MRI data acquired in a mid-ventricular left ventricle section in N=8 swine subjects. The proposed method shows a more physiological distribution of myofiber strains compared to standard approaches that directly differentiate the displacement field.

Magnetic Resonance-Based Characterization of Myocardial Architecture

Advances in technology have made it possible to image the microstructure of the heart with diffusion-weighted magnetic resonance. The technique provides unique insights into the cellular architecture of the myocardium and how this is perturbed in a range of disease contexts. In this review, the physical basis of diffusion MRI and the challenges of implementing it in the beating heart are discussed. Cutting edge acquisition and analysis techniques, as well as the results of initial clinical studies, are reported.

Fluorescence microscopy tensor imaging representations for large-scale dataset analysis

Understanding complex biological systems requires the system-wide characterization of cellular and molecular features. Recent advances in optical imaging technologies and chemical tissue clearing have facilitated the acquisition of whole-organ imaging datasets, but automated tools for their quantitative analysis and visualization are still lacking. We have here developed a visualization technique capable of providing whole-organ tensor imaging representations of local regional descriptors based on fluorescence data acquisition. This method enables rapid, multiscale, analysis and virtualization of large-volume, high-resolution complex biological data while generating 3D tractographic representations. Using the murine heart as a model, our method allowed us to analyze and interrogate the cardiac microvasculature and the tissue resident macrophage distribution and better infer and delineate the underlying structural network in unprecedented detail.

Accelerating cardiac diffusion tensor imaging combining local low-rank and 3D TV constraint

OBJECTIVE

Diffusion tensor magnetic resonance imaging (DT-MRI, or DTI) is a promising technique for invasively probing biological tissue structures. However, DTI is known to suffer from much longer acquisition time with respect to conventional MRI and the problem is worsened when dealing with in vivo acquisitions. Therefore, faster DTI for both ex vivo and in vivo scans is highly desired.

MATERIALS AND METHODS

This paper proposes a new compressed sensing (CS) reconstruction method that employs local low-rank (LLR) model and three-dimensional (3D) total variation (TV) constraint to reconstruct cardiac diffusion-weighted (DW) images from highly undersampled k-space data. The LLR model takes the set of DW images corresponding to different diffusion gradient directions as a 3D image volume and decomposes the latter into overlapping 3D blocks. Then, the 3D blocks are stacked as two-dimensional (2D) matrix. Finally, low-rank property is applied to each block matrix and the 3D TV constraint to the 3D image volume. The underlying constrained optimization problem is finally solved using the first-order fast method. The proposed method is evaluated on real ex vivo cardiac DTI data as a prerequisite to in vivo cardiac DTI applications.

RESULTS

The results on real human ex vivo cardiac DTI images demonstrate that the proposed method exhibits lower reconstruction errors for DTI indices, including fractional anisotropy (FA), mean diffusivities (MD), transverse angle (TA), and helix angle (HA), compared to existing CS-based DTI image reconstruction techniques.

CONCLUSION

The proposed method provides better reconstruction quality and more accurate DTI indices in comparison with the state-of-the-art CS-based DW image reconstruction methods.

A rule-based method for predicting the electrical activation of the heart with cardiac resynchronization therapy from non-invasive clinical data

Background

Cardiac Resynchronization Therapy (CRT) is one of the few effective treatments for heart failure patients with ventricular dyssynchrony. The pacing location of the left ventricle is indicated as a determinant of CRT outcome.

Objective

Patient specific computational models allow the activation pattern following CRT implant to be predicted and this may be used to optimize CRT lead placement.

Methods

In this study, the effects of heterogeneous cardiac substrate (scar, fast endocardial conduction, slow septal conduction, functional block) on accurately predicting the electrical activation of the LV epicardium were tested to determine the minimal detail required to create a rule based model of cardiac electrophysiology. Non-invasive clinical data (CT or CMR images and 12 lead ECG) from eighteen patients from two centers were used to investigate the models.

Results

Validation with invasive electro-anatomical mapping data identified that computer models with fast endocardial conduction were able to predict the electrical activation with a mean distance errors of 9.2 ± 0.5 mm (CMR data) or (CT data) 7.5 ± 0.7 mm.

Conclusion

This study identified a simple rule-based fast endocardial conduction model, built using non-invasive clinical data that can be used to rapidly and robustly predict the electrical activation of the heart. Pre-procedural prediction of the latest electrically activating region to identify the optimal LV pacing site could potentially be a useful clinical planning tool for CRT procedures.

High resolution in-vivo DT-CMR using an interleaved variable density spiral STEAM sequence

PURPOSE

Diffusion tensor cardiovascular magnetic resonance (DT-CMR) has a limited spatial resolution. The purpose of this study was to demonstrate high-resolution DT-CMR using a segmented variable density spiral sequence with correction for motion, off-resonance, and T2*-related blurring.

METHODS

A single-shot stimulated echo acquisition mode (STEAM) echo-planar-imaging (EPI) DT-CMR sequence at 2.8 × 2.8 × 8 mm³ and 1.8 × 1.8 × 8 mm³ was compared to a single-shot spiral at 2.8 × 2.8 × 8 mm³ and an interleaved spiral sequence at 1.8 × 1.8 × 8 mm³ resolution in 10 healthy volunteers at peak systole and diastasis. Motion-induced phase was corrected using the densely sampled central k-space data of the spirals. STEAM field maps and T2* measures were obtained using a pair of stimulated echoes each with a double spiral readout, the first used to correct the motion-induced phase of the second.

RESULTS

The high-resolution spiral sequence produced similar DT-CMR results and quality measures to the standard-resolution sequence in both cardiac phases. Residual differences in fractional anisotropy and helix angle gradient between the resolutions could be attributed to spatial resolution and/or signal-to-noise ratio. Data quality increased after both motion-induced phase correction and off-resonance correction, and sharpness increased after T2* correction. The high-resolution EPI sequence failed to provide sufficient data quality for DT-CMR reconstruction.

CONCLUSION

In this study, an in vivo DT-CMR acquisition at 1.8 × 1.8 mm² in-plane resolution was demonstrated using a segmented spiral STEAM sequence. Motion-induced phase and off-resonance corrections are essential for high-resolution spiral DT-CMR. Segmented variable density spiral STEAM was found to be the optimal method for acquiring high-resolution DT-CMR data.

Novel insights into in-vivo diffusion tensor cardiovascular magnetic resonance using computational modeling and a histology-based virtual microstructure

Purpose

To develop histology‐informed simulations of diffusion tensor cardiovascular magnetic resonance (DT‐CMR) for typical in‐vivo pulse sequences and determine their sensitivity to changes in extra‐cellular space (ECS) and other microstructural parameters.

Methods

We synthesised the DT‐CMR signal from Monte Carlo random walk simulations. The virtual tissue was based on porcine histology. The cells were thickened and then shrunk to modify ECS. We also created idealised geometries using cuboids in regular arrangement, matching the extra‐cellular volume fraction (ECV) of 16–40%. The simulated voxel size was 2.8 × 2.8 × 8.0 mm³ for pulse sequences covering short and long diffusion times: Stejskal–Tanner pulsed‐gradient spin echo, second‐order motion‐compensated spin echo, and stimulated echo acquisition mode (STEAM), with clinically available gradient strengths.

Results

The primary diffusion tensor eigenvalue increases linearly with ECV at a similar rate for all simulated geometries. Mean diffusivity (MD) varies linearly, too, but is higher for the substrates with more uniformly distributed ECS. Fractional anisotropy (FA) for the histology‐based geometry is higher than the idealised geometry with low sensitivity to ECV, except for the long mixing time of the STEAM sequence. Varying the intra‐cellular diffusivity (D_IC) results in large changes of MD and FA. Varying extra‐cellular diffusivity or using stronger gradients has minor effects on FA. Uncertainties of the primary eigenvector orientation are reduced using STEAM.

Conclusions

We found that the distribution of ECS has a measurable impact on DT‐CMR parameters. The observed sensitivity of MD and FA to ECV and D_IC has potentially interesting applications for interpreting in‐vivo DT‐CMR parameters.

Computational Modeling for Cardiac Resynchronization Therapy

Cardiac resynchronization therapy (CRT) is an effective treatment for heart failure (HF) patients with an electrical substrate pathology causing ventricular dyssynchrony. However 40-50% of patients do not respond to treatment. Cardiac modeling of the electrophysiology, electromechanics, and hemodynamics of the heart has been used to study mechanisms behind HF pathology and CRT response. Recently, multi-scale dyssynchronous HF models have been used to study optimal device settings and optimal lead locations, investigate the underlying cardiac pathophysiology, as well as investigate emerging technologies proposed to treat cardiac dyssynchrony. However the breadth of patient and experimental data required to create and parameterize these models and the computational resources required currently limits the use of these models to small patient numbers. In the future, once these technical challenges are overcome, biophysically based models of the heart have the potential to become a clinical tool to aid in the diagnosis and treatment of HF.

Myocardial Scar Delineation Using Diffusion Tensor Magnetic Resonance Tractography

Background

Late gadolinium enhancement (LGE) is the current standard for myocardial scar delineation. In this study, we introduce the tractographic propagation angle (PA), a metric of myofiber curvature (degrees/unit distance) derived from diffusion tensor imaging (DTI), and compare its use to LGE and invasive scar assessment by endocardial voltage mapping.

Methods and Results

DTI was performed on 7 healthy human volunteers, 5 patients with myocardial infarction, 6 normal mice, and 7 mice with myocardial infarction. LGE to delineate the infarct and border zones was performed with a 2‐dimensional inversion recovery gradient‐echo sequence. Ex vivo DTI was performed on 5 normal human and 5 normal sheep hearts. Endocardial electroanatomic mapping and subsequent ex vivo DTI was performed on 5 infarcted sheep hearts. PA in the normal human hearts varied smoothly and was generally <4. The mean PA in the infarct zone was significantly elevated (10.34±1.02 versus 4.05±0.45, P<0.05). Regions with a PA ≤4 consistently had a bipolar voltage ≥1.5 mV, whereas those with PA values between 4 and 10 had voltages between 0.5 and 1.5 mV. A PA threshold >4 was the most accurate DTI‐derived measure of infarct size and demonstrated the greatest correlation with LGE (r=0.95).

Conclusions

We found a strong correlation between infarct size by PA and LGE in both mice and humans. There was also an inverse relationship between PA values and endocardial voltage. The use of PA may enable myocardial scar delineation and characterization of arrhythmogenic substrate without the need for exogenous contrast agents.

Diffusion tensor cardiovascular magnetic resonance with a spiral trajectory: An in vivo comparison of echo planar and spiral stimulated echo sequences

PURPOSE

Diffusion tensor cardiovascular MR (DT-CMR) using stimulated echo acquisition mode (STEAM) with echo-planar-imaging (EPI) readouts is a low signal-to-noise-ratio (SNR) technique and therefore typically has a low spatial resolution. Spiral trajectories are more efficient than EPI, and could increase the SNR. The purpose of this study was to compare the performance of a novel STEAM spiral DT-CMR sequence with an equivalent established EPI technique.

METHODS

A STEAM DT-CMR sequence was implemented with a spiral readout and a reduced field of view. An in vivo comparison of DT-CMR parameters and data quality between EPI and spiral was performed in 11 healthy volunteers imaged in peak systole and diastasis at 3 T. The SNR was compared in a phantom and in vivo.

RESULTS

There was a greater than 49% increase in the SNR in vivo and in the phantom measurements (in vivo septum, systole: SNR_EPI  = 8.0 ± 2.2, SNR_spiral  = 12.0 ± 2.7; diastasis: SNR_EPI  = 8.1 ± 1.6, SNR_spiral  = 12.0 ± 3.7). There were no significant differences in helix angle gradient (HAG) (systole: HAG_EPI  = -0.79 ± 0.07 °/%; HAG_spiral  = -0.74 ± 0.16 °/%; P = 0.11; diastasis: HAG_EPI  = -0.63 ± 0.05 °/%; HAG_spiral  = -0.56 ± 0.14 °/%; P = 0.20), mean diffusivity (MD) in systole (MD_EPI  = 0.99 ± 0.06 × 10^-3 mm² /s, MD_spiral  = 1.00 ± 0.09 × 10^-3 mm² /s, P = 0.23) and secondary eigenvector angulation (E2A) (systole: E2A_EPI  = 61 ± 10 °; E2A_spiral  = 63 ± 10 °; P = 0.77; diastasis: E2A_EPI  = 18 ± 11 °; E2A_spiral  = 15 ± 8 °; P = 0.20) between the sequences. There was a small difference (≈ 20%) in fractional anisotropy (FA) (systole: FA_EPI  = 0.49 ± 0.03, FA_spiral  = 0.41 ± 0.04; P < 0.01; diastasis: FA_EPI  = 0.66 ± 0.05, FA_spiral  = 0.55 ± 0.03; P < 0.01) and mean diffusivity in diastasis (10%; MD_EPI  = 1.00 ± 0.12 × 10^-3 mm² /s, MD_spiral  = 1.10 ± 0.09 × 10^-3 mm² /s, P = 0.02).

CONCLUSION

This is the first study to demonstrate DT-CMR STEAM using a spiral trajectory. The SNR was increased by using a spiral rather than the more established EPI readout, and the DT-CMR parameters were largely similar between the two sequences. Magn Reson Med 80:648-654, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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.

Diffusion MRI in the heart

Diffusion MRI provides unique information on the structure, organization, and integrity of the myocardium without the need for exogenous contrast agents. Diffusion MRI in the heart, however, has proven technically challenging because of the intrinsic non-rigid deformation during the cardiac cycle, displacement of the myocardium due to respiratory motion, signal inhomogeneity within the thorax, and short transverse relaxation times. Recently developed accelerated diffusion-weighted MR acquisition sequences combined with advanced post-processing techniques have improved the accuracy and efficiency of diffusion MRI in the myocardium. In this review, we describe the solutions and approaches that have been developed to enable diffusion MRI of the heart in vivo, including a dual-gated stimulated echo approach, a velocity- (M1 ) or an acceleration- (M2 ) compensated pulsed gradient spin echo approach, and the use of principal component analysis filtering. The structure of the myocardium and the application of these techniques in ischemic heart disease are also briefly reviewed. The advent of clinical MR systems with stronger gradients will likely facilitate the translation of cardiac diffusion MRI into clinical use. The addition of diffusion MRI to the well-established set of cardiovascular imaging techniques should lead to new and complementary approaches for the diagnosis and evaluation of patients with heart disease. © 2015 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.

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.

Image-Based Investigation of Human in Vivo Myofibre Strain

Cardiac myofibre deformation is an important determinant of the mechanical function of the heart. Quantification of myofibre strain relies on 3D measurements of ventricular wall motion interpreted with respect to the tissue microstructure. In this study, we estimated in vivo myofibre strain using 3D structural and functional atlases of the human heart. A finite element modelling framework was developed to incorporate myofibre orientations of the left ventricle (LV) extracted from 7 explanted normal human hearts imaged ex vivo with diffusion tensor magnetic resonance imaging (DTMRI) and kinematic measurements from 7 normal volunteers imaged in vivo with tagged MRI. Myofibre strain was extracted from the DTMRI and 3D strain from the tagged MRI. We investigated: i) the spatio-temporal variation of myofibre strain throughout the cardiac cycle; ii) the sensitivity of myofibre strain estimates to the variation in myofibre angle between individuals; and iii) the sensitivity of myofibre strain estimates to variations in wall motion between individuals. Our analysis results indicate that end systolic (ES) myofibre strain is approximately homogeneous throughout the entire LV, irrespective of the inter-individual variation in myofibre orientation. Additionally, inter-subject variability in myofibre orientations has greater effect on the variabilities in myofibre strain estimates than the ventricular wall motions. This study provided the first quantitative evidence of homogeneity of ES myofibre strain using minimally-invasive medical images of the human heart and demonstrated that image-based modelling framework can provide detailed insight to the mechanical behaviour of the myofibres, which may be used as a biomarker for cardiac diseases that affect cardiac mechanics.

Relationship between cardiac diffusion tensor imaging parameters and anthropometrics in healthy volunteers

BACKGROUND

In vivo cardiac diffusion tensor imaging (cDTI) is uniquely capable of interrogating laminar myocardial dynamics non-invasively. A comprehensive dataset of quantative parameters and comparison with subject anthropometrics is required.

METHODS

cDTI was performed at 3T with a diffusion weighted STEAM sequence. Data was acquired from the mid left ventricle in 43 subjects during the systolic and diastolic pauses. Global and regional values were determined for fractional anisotropy (FA), mean diffusivity (MD), helix angle gradient (HAg, degrees/%depth) and the secondary eigenvector angulation (E2A). Regression analysis was performed between global values and subject anthropometrics.

RESULTS

All cDTI parameters displayed regional heterogeneity. The RR interval had a significant, but clinically small effect on systolic values for FA, HAg and E2A. Male sex and increasing left ventricular end diastolic volume were associated with increased systolic HAg. Diastolic HAg and systolic E2A were both directly related to left ventricular mass and body surface area. There was an inverse relationship between E2A mobility and both age and ejection fraction.

CONCLUSIONS

Future interpretations of quantitative cDTI data should take into account anthropometric variations observed with patient age, body surface area and left ventricular measurements. Further work determining the impact of technical factors such as strain and SNR is required.

Higher-Order Motion-Compensation for In Vivo Cardiac Diffusion Tensor Imaging in Rats

Motion of the heart has complicated in vivo applications of cardiac diffusion MRI and diffusion tensor imaging (DTI), especially in small animals such as rats where ultra-high-performance gradient sets are currently not available. Even with velocity compensation via, for example, bipolar encoding pulses, the variable shot-to-shot residual motion-induced spin phase can still give rise to pronounced artifacts. This study presents diffusion-encoding schemes that are designed to compensate for higher-order motion components, including acceleration and jerk, which also have the desirable practical features of minimal TEs and high achievable b-values. The effectiveness of these schemes was verified numerically on a realistic beating heart phantom, and demonstrated empirically with in vivo cardiac diffusion MRI in rats. Compensation for acceleration, and lower motion components, was found to be both necessary and sufficient for obtaining diffusion-weighted images of acceptable quality and SNR, which yielded the first in vivo cardiac DTI demonstrated in the rat. These findings suggest that compensation for higher order motion, particularly acceleration, can be an effective alternative solution to high-performance gradient hardware for improving in vivo cardiac DTI.

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.

Second-order motion-compensated spin echo diffusion tensor imaging of the human heart

PURPOSE

Myocardial microstructure has been challenging to probe in vivo. Spin echo-based diffusion-weighted sequences allow for single-shot acquisitions but are highly sensitive to cardiac motion. In this study, the use of second-order motion-compensated diffusion encoding was compared with first-order motion-compensated diffusion-weighted imaging during systolic contraction of the heart.

METHODS

First- and second-order motion-compensated diffusion encoding gradients were incorporated into a triggered single-shot spin echo sequence. The effect of contractile motion on the apparent diffusion coefficients and tensor orientations was investigated in vivo from basal to apical level of the heart.

RESULTS

Second-order motion compensation was found to increase the range of systolic trigger delays from 30%-55% to 15%-77% peak systole at the apex and from 25%-50% to 15%-79% peak systole at the base. Diffusion tensor analysis yielded more physiological transmural distributions when using second-order motion-compensated diffusion tensor imaging.

CONCLUSION

Higher-order motion-compensated diffusion encoding decreases the sensitivity to cardiac motion, thereby enabling cardiac DTI over a wider range of time points during systolic contraction of the heart.

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.

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.

Arrhythmogenic right ventricular cardiomyopathy mutations alter shear response without changes in cell-cell adhesion

AIMS

The majority of patients diagnosed with arrhythmogenic right ventricular cardiomyopathy (ARVC) have mutations in genes encoding desmosomal proteins, raising the possibility that abnormal intercellular adhesion plays an important role in disease pathogenesis. We characterize cell mechanical properties and molecular responses to oscillatory shear stress in cardiac myocytes expressing mutant forms of the desmosomal proteins, plakoglobin and plakophilin, which are linked to ARVC in patients.

METHODS AND RESULTS

Cells expressing mutant plakoglobin or plakophilin showed no differences in cell-cell adhesion relative to controls, while knocking down these proteins weakened cell-cell adhesion. However, cells expressing mutant plakoglobin failed to increase the amount of immunoreactive signal for plakoglobin or N-cadherin at cell-cell junctions in response to shear stress, as seen in control cells. Cells expressing mutant plakophilin exhibited a similar attenuation in the shear-induced increase in junctional plakoglobin immunoreactive signal in response to shear stress, suggesting that the phenotype is independent of the type of mutant protein being expressed. Cells expressing mutant plakoglobin also showed greater myocyte apoptosis compared with controls. Apoptosis rates increased greatly in response to shear stress in cells expressing mutant plakoglobin, but not in controls. Abnormal responses to shear stress in cells expressing either mutant plakoglobin or plakophilin could be reversed by SB216763, a GSK3β inhibitor.

CONCLUSIONS

Desmosomal mutations linked to ARVC do not significantly affect cell mechanical properties, but cause myocytes to respond abnormally to mechanical stress through a mechanism involving GSK3β. These results may help explain why patients with ARVC experience disease exacerbations following strenuous exercise.

Optimal diffusion weighting for in vivo cardiac diffusion tensor imaging

PURPOSE

To investigate the influence of the diffusion weighting on in vivo cardiac diffusion tensor imaging (cDTI) and obtain optimal parameters.

METHODS

Ten subjects were scanned using stimulated echo acquisition mode echo planar imaging with six b-values, from 50 to 950 s·mm(-2) , plus b = 15 s·mm(-2) reference. The relationship between b-value and both signal loss and signal-to-noise ratio measures was investigated. Mean diffusivity, fractional anisotropy, and helical angle maps were calculated using all possible b-value pairs to investigate the effects of diffusion weighting on the main and reference data.

RESULTS

Signal decay at low b-values was dominated by processes with high apparent diffusion coefficients, most likely microvascular perfusion. This effect could be avoided by diffusion weighting of the reference images. Parameter maps were improved with increased b-value until the diffusion-weighted signal approached the noise floor. For the protocol used in this study, b = 750 s·mm(-2) combined with 150 s·mm(-2) diffusion weighting of the reference images proved optimal.

CONCLUSION

Mean diffusivity, fractional anisotropy, and helical angle from cDTI are influenced by the b-value of the main and reference data. Using optimal values improves parameter maps and avoids microvascular perfusion effects. This optimized protocol should provide greater sensitivity to pathological changes in parameter maps.

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.

Accelerated MR diffusion tensor imaging using distributed compressed sensing

PURPOSE

Diffusion tensor imaging (DTI) is known to suffer from long acquisition time in the orders of several minutes or even hours. Therefore, a feasible way to accelerate DTI data acquisition is highly desirable. In this article, the feasibility and efficacy of distributed compressed sensing to fast DTI is investigated by exploiting the joint sparsity prior in diffusion-weighted images.

METHODS

Fully sampled DTI datasets were obtained from both simulated phantom and experimental heart sample, with diffusion gradient applied in six directions. The k-space data were undersampled retrospectively with acceleration factors from 2 to 6. Diffusion-weighted images were reconstructed by solving an l2-l1 norm minimization problem. Reconstruction performance with varied signal-to-noise ratio and acceleration factors were evaluated by root-mean-square error and maps of reconstructed DTI indices.

RESULTS

Superiority of distributed compressed sensing over basic compressed sensing was confirmed with simulation, and the reconstruction accuracy was influenced by signal-to-noise ratio and acceleration factors. Experimental results demonstrate that DTI indices including fractional anisotropy, mean diffusivities, and orientation of primary eigenvector can be obtained with high accuracy at acceleration factors up to 4.

CONCLUSION

Distributed compressed sensing is shown to be able to accelerate DTI and may be used to reduce DTI acquisition time practically.

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.

The diagnosis of hypertrophic cardiomyopathy by cardiovascular magnetic resonance

Hypertrophic cardiomyopathy (HCM) is the most common genetic disease of the heart. HCM is characterized by a wide range of clinical expression, ranging from asymptomatic mutation carriers to sudden cardiac death as the first manifestation of the disease. Over 1000 mutations have been identified, classically in genes encoding sarcomeric proteins. Noninvasive imaging is central to the diagnosis of HCM and cardiovascular magnetic resonance (CMR) is increasingly used to characterize morphologic, functional and tissue abnormalities associated with HCM. The purpose of this review is to provide an overview of the clinical, pathological and imaging features relevant to understanding the diagnosis of HCM. The early and overt phenotypic expression of disease that may be identified by CMR is reviewed. Diastolic dysfunction may be an early marker of the disease, present in mutation carriers prior to the development of left ventricular hypertrophy (LVH). Late gadolinium enhancement by CMR is present in approximately 60% of HCM patients with LVH and may provide novel information regarding risk stratification in HCM. It is likely that integrating genetic advances with enhanced phenotypic characterization of HCM with novel CMR techniques will importantly improve our understanding of this complex disease.

Quantitative comparison of myocardial fiber structure between mice, rabbit, and sheep using diffusion tensor cardiovascular magnetic resonance

BACKGROUND

Accurate interpretations of cardiac functions require precise structural models of the myocardium, but the latter is not available always and for all species. Although scaling or substitution of myocardial fiber information from alternate species has been used in cardiac functional modeling, the validity of such practice has not been tested.

METHODS

Fixed mouse (n = 10), rabbit (n = 6), and sheep (n = 5) hearts underwent diffusion tensor imaging (DTI). The myocardial structures in terms of the left ventricular fiber orientation helix angle index were quantitatively compared between the mouse rabbit and sheep hearts.

RESULTS

The results show that significant fiber structural differences exist between any two of the three species. Specifically, the subepicardial fiber orientation, and the transmural range and linearity of fiber helix angles are significantly different between the mouse and either rabbit or sheep. Additionally, a significant difference was found between the transmural helix angle range between the rabbit and sheep. Across different circumferential regions of the heart, the fiber orientation was not found to be significantly different.

CONCLUSIONS

The current study indicates that myocardial structural differences exist between different size hearts. An immediate implication of the present findings for myocardial structural or functional modeling studies is that caution must be exercised when extrapolating myocardial structures from one species to another.

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.

Oblique 3D MRI tags for the estimation of true 3D cardiac motion parameters

Aim of this study is to demonstrate the advantages of oblique 3D tags in cardiac magnetic resonance imaging (MRI) and the potential to accurately describe the complex motion of the myocardial wall. 3D cardiac Cine data were densely tagged with 3D oblique tags. The latter were tracked using Gabor analysis and active geometries. From the tag intersections, common 2D parameters such as long axis shortening, radial shortening and rotation were evaluated on a global as well as detailed local level. Finally, the same data were used to estimate left ventricular volume change and myocardial stress/strain. We have successfully tracked dense 3D tags and evaluated common parameters on a detailed local level. In addition, inherently 3D parameters could be estimated. Global motion data are in accordance with previously published data. Oblique tags allow for unambiguous localization of the tag plane in all MRI slices and in any time frame. In contrast to HARP, our tag tracking methodology allows for tracking of the tags even when they are dense. Motion parameters can be extracted in greater detail. Moreover, the intersections of dense oblique 3D tags provide a natural basis for a finite element model of the heart. Straight forward access to the 3D characteristics of the cardiac motion is provided.

The correlation of 3D DT-MRI fiber disruption with structural and mechanical degeneration in porcine myocardium

Evaluation of structural parameters following a myocardial infarction (MI) is important to assess left ventricular function and remodeling. In this study, we assessed the capability of 3D diffusion tensor magnetic resonance imaging (DT-MRI) to assess tissue degeneration shortly after an MI using a porcine model of infarction. Two days after an induced infarction, hearts were explanted and immediately scanned by a 3T MRI scanner with a diffusion tensor imaging protocol. 3D fiber tracks and clustering models were generated from the diffusion-weighted imaging data. We found in a normal explanted heart that DT-MRI fibers showed a multilayered helical structure, with fiber architecture and fiber density reflecting the integrity of muscle fibers. For infarcted heart explants, we observed either a lack of fibers or disruption of fibers in the infarcted regions. Contours of the disrupted DT-MRI fibers were found to be consistent with the infarcted regions. Both histological and mechanical analysis of the infarcted hearts suggested DT-MRI fiber disruption correlated with altered microstructure and tissue mechanics. The ability of 3D DT-MRI to accurately distinguish viable myocardium from dead myocardium only 2 days post infarct without the use of radioisotopes or ionotropic agents makes it a promising approach to evaluate cardiac damage early post-MI.

Resolving myoarchitectural disarray in the mouse ventricular wall with diffusion spectrum magnetic resonance imaging

The myoarchitecture of the ventricular wall provides a structural template dictating tissue-scale patterns of mechanical function. We studied whether myofiber tract imaging performed with MR diffusion spectrum imaging (DSI) tractography has the capacity to resolve abnormalities of ventricular myoarchitecture in a model of congenital hypertrophic cardiomyopathy (HCM) associated with the ablation of myosin binding protein-C (MyBP-C). Homozygous MyBP-C knockout mice were generated by deletion of exons 3-10 from the endogenous MyBP-C gene. Fiber alignment in the left ventricular wall of wild type mice was depicted through DSI tractography (and confirmed by multi-slice two-photon microscopy) as a set of helical structures whose angles display a continuous transition from negative in the subepicardium to positive in the subendocardium. In contrast, the hearts obtained from the MyBP-C knockouts displayed substantial myoarchitectural disarray, characterized by a loss of voxel-to-voxel orientational coherence for fibers principally located in the mid-myocardium-subendocardium and impairment of the transmural progression of helix angles. These results substantiate the use of DSI tractography in determining myoarchitectural disarray in models of cardiomyopathy and suggest a biological association between myofilament expression, cardiac fiber alignment, and torsional rotation in the setting of congenital HCM.

In vivo human 3D cardiac fibre architecture: reconstruction using curvilinear interpolation of diffusion tensor images

In vivo imaging of the cardiac 3D fibre architecture is still a challenge, but it would have many clinical applications, for instance to better understand pathologies and to follow up remodelling after therapy. Recently, cardiac MRI enabled the acquisition of Diffusion Tensor images (DTI) of 2D slices. We propose a method for the complete 3D reconstruction of cardiac fibre architecture in the left ventricular myocardium from sparse in vivo DTI slices. This is achieved in two steps. First we map non-linearly the left ventricular geometry to a truncated ellipsoid. Second, we express coordinates and tensor components in Prolate Spheroidal System, where an anisotropic Gaussian kernel regression interpolation is performed. The framework is initially applied to a statistical cardiac DTI atlas in order to estimate the optimal anisotropic bandwidths. Then, it is applied to in vivo beating heart DTI data sparsely acquired on a healthy subject. Resulting in vivo tensor field shows good correlation with literature, especially the elevation (helix) angle transmural variation. To our knowledge, this is the first reconstruction of in vivo human 3D cardiac fibre structure. Such approach opens up possibilities in terms of analysis of the fibre architecture in patients.

Diffusion MR tractography of the heart

Histological studies have shown that the myocardium consists of an array of crossing helical fiber tracts. Changes in myocardial fiber architecture occur in ischemic heart disease and heart failure, and can be imaged non-destructively with diffusion-encoded MR. Several diffusion-encoding schemes have been developed, ranging from scalar measurements of mean diffusivity to a 6-dimensional imaging technique known as diffusion spectrum imaging or DSI. The properties of DSI make it particularly suited to the generation of 3-dimensional tractograms of myofiber architecture. In this article we review the physical basis of diffusion-tractography in the myocardium and the attributes of the available techniques, placing particular emphasis on DSI. The application of DSI in ischemic heart disease is reviewed, and the requisites for widespread clinical translation of diffusion MR tractography in the heart are discussed.

Heterogeneity of left ventricular wall thickening mechanisms

BACKGROUND

Myocardial fibers are grouped into lamina (or sheets) 3 to 4 cells thick. Fiber shortening produces systolic left ventricular (LV) wall thickening primarily by laminar extension, thickening, and shear, but the regional variability and transmural distribution of these 3 mechanisms are incompletely understood.

METHODS AND RESULTS

Nine sheep had transmural radiopaque markers inserted into the anterior basal and lateral equatorial LV. Four-dimensional marker dynamics were studied with biplane videofluoroscopy to measure circumferential, longitudinal, and radial systolic strains in the epicardium, midwall, and endocardium. Fiber and sheet angles from quantitative histology allowed transformation of these strains into transmural contributions of sheet extension, thickening, and shear to systolic wall thickening. At all depths, systolic wall thickening in the anterior basal region was 1.6 to 1.9 times that in the lateral equatorial region. Interestingly, however, systolic fiber shortening was identical at each transmural depth in these regions. Endocardial anterior basal sheet thickening was >2 times greater than in the lateral equatorial region (epicardium, 0.16+/-0.15 versus 0.03+/-0.06; endocardium, 0.45+/-0.40 versus 0.17+/-0.09). Midwall sheet extension was >2 times that in the lateral wall (0.22+/-0.12 versus 0.09+/-0.06). Epicardial and midwall sheet shears in the anterior wall were approximately 2 times higher than in the lateral wall (epicardium, 0.14+/-0.07 versus 0.05+/-0.03; midwall, 0.21+/-0.12 versus 0.12+/-0.06).

CONCLUSIONS

These data demonstrate fundamentally different regional contributions of laminar mechanisms for amplifying fiber shortening to systolic wall thickening. Systolic fiber shortening was identical at each transmural depth in both the anterior and lateral LV sites. However, systolic wall thickening of the anterior site was much greater than that of the lateral site. Fiber shortening drives systolic wall thickening, but sheet dynamics and orientations are of great importance to systolic wall thickening. LV wall thickening and its clinical implications pivot on different wall thickening mechanisms in various LV regions. Attempts to implant healthy contractile cells into diseased hearts or to surgically manipulate LV geometry need to take into account not only cardiomyocyte contraction but also transmural LV intercellular architecture and geometry.

Associating the mesoscale fiber organization of the tongue with local strain rate during swallowing

The tongue is an intricately configured muscular organ that undergoes a stereotypical set of deformations during the course of normal human swallowing. In order to demonstrate quantitatively the relationship between 3D aligned lingual fiber organization and mechanics during swallowing, the tissue's myoarchitecture and strain rate were imaged before and during the propulsive phase of a 3.0ml water bolus swallow. Mesoscale fiber organization was imaged with high-resolution diffusion tensor imaging (DTI) and multi-voxel myofiber tracts generated along maximum diffusion vectors. Tissue compression/expansion was obtained via lingual pressure-gated phase-contrast (PC) MRI, a method which determines local strain rate as a function of the phase shift occurring along an applied gradient vector. The co-alignment of myofiber tract direction and the localized principal strain rate vectors was obtained by translating the strain rate tensor into the reference frame with the primary axis parallel to the maximum diffusion vector using Mohr's circle, resulting in the generation of fiber-aligned strain rate (FASR). DTI tractography displayed the complete fiber anatomy of the tongue, consisting of a core region of orthogonally aligned fibers encased within a longitudinal sheath, which merge with the externally connected styloglossus, hyoglossus, and genioglossus fibers. FASR images obtained in the mid-sagittal plane demonstrated that bolus propulsion was associated with prominent compressive strain aligned with the genioglossus muscle combined with expansive strain aligned with the verticalis and geniohyoid muscles. These data demonstrate that lingual deformation during swallowing involves complex interactions involving intrinsic and extrinsic muscles, whose contractility is directed by the alignment of mesoscale fiber tracts.

MR diffusion tensor imaging study of postinfarct myocardium structural remodeling in a porcine model

This study aimed to investigate postinfarct left ventricular (LV) fiber structural alterations by ex vivo diffusion tensor imaging (DTI) in a porcine heart model. In vivo cardiac MR imaging was first performed to measure ventricular function in six adult pigs with septal infarction near apex induced by the LAD ligation 13 weeks earlier. Hearts were then excised from the infarct pigs (n = 6) and six intact controls (n = 6) and fixed in formalin. High-resolution DTI was employed to examine changes in fractional anisotropy (FA), apparent diffusion coefficient (ADC), and transmural helix angle distribution in the infarct, adjacent and remote regions as compared to the sham regions in the controls. FA values were found to decrease in the infarct and differ between the adjacent and remote regions. ADC increase in the infarct region was substantial, while changes in the adjacent and remote regions were insignificant. Structurally, the double-helix myocardial structure shifted toward more left-handed around the infarcted myocardium. Accordingly, the histological analysis revealed clear fiber structural degradation in the adjacent region. These findings confirmed the subtle alterations in the myocardial fiber quality and structure not only in the infarcted but also in the surrounding noninfarcted myocardium or borderzone.

Reconstruction and visualization of fiber and laminar structure in the normal human heart from ex vivo diffusion tensor magnetic resonance imaging (DTMRI) data

OBJECTIVE

The human heart is composed of a helical network of muscle fibers organized to form sheets that are separated by cleavage planes responsible for the orthotropic mechanical properties of cardiac muscle. The purpose of this study is the reconstruction and visualization of these structures in 3 dimensions.

METHODS

Anisotropic least square filtering followed by fiber and sheet tracking techniques were applied to diffusion tensor magnetic resonance imaging data of the excised human heart. Fibers were reconstructed using the first eigenvectors of the diffusion tensors. The sheets were reconstructed using the second and third eigenvectors and visualized as surfaces.

RESULTS

The fibers are shown to lie in sheets that have transmural structure, which correspond to histologic studies published in the literature. Quantitative measurements show that the sheets as appose to the fibers are organized into laminar orientations without dominant populations.

CONCLUSIONS

A visualization algorithm was developed to demonstrate the complex 3-dimensional orientation of the fibers and sheets in human myocardium.

Laminar Arrangement of Ventricular Myocytes Influences Electrical Behavior of the Heart

The response of the heart to electrical shock, electrical propagation in sinus rhythm, and the spatiotemporal dynamics of ventricular fibrillation all depend critically on the electrical anisotropy of cardiac tissue. A long-held view of cardiac electrical anisotropy is that electrical conductivity is greatest along the myocyte axis allowing most rapid propagation of electrical activation in this direction, and that conductivity is isotropic transverse to the myocyte axis supporting a slower uniform spread of activation in this plane. In this context, knowledge of conductivity in two directions, parallel and transverse to the myofiber axis, is sufficient to characterize the electrical action of the heart. Here we present new experimental data that challenge this view. We have used a novel combination of intramural electrical mapping, and experiment-specific computer modeling, to demonstrate that left ventricular myocardium has unique bulk conductivities associated with three microstructurally-defined axes. We show that voltage fields induced by intramural current injection are influenced by not only myofiber direction, but also the transmural arrangement of muscle layers or myolaminae. Computer models of these experiments, in which measured 3D tissue structure was reconstructed in-silico, best matched recorded voltages with conductivities in the myofiber direction, and parallel and normal to myolaminae, set in the ratio 4:2:1, respectively. These findings redefine cardiac tissue as an electrically orthotropic substrate and enhance our understanding of how external shocks may act to successfully reset the fibrillating heart into a uniform electrical state. More generally, the mechanisms governing the destabilization of coordinated electrical propagation into ventricular arrhythmia need to be evaluated in the light of this discovery.

Regional localisation of left ventricular sheet structure: integration with current models of cardiac fibre, sheet and band structure

The architecture of the heart remains controversial despite extensive effort and recent advances in imaging techniques. Several opposing and non-mutually compatible models have been proposed to explain cardiac structure, and these models, although limited, have advanced the study and understanding of heart structure, function and development. We describe key areas of similarity and difference, highlight areas of contention and point to the important limitations of these models. Recent research in animal models on the nature, geometry and interaction of cardiac sheet structure allows unification of some seemingly conflicting features of the structural models. Intriguingly, evidence points to significant inter-individual structural variability (within constrained limits) in the canine, leading to the concept of a continuum (or distribution) of cardiac structures. This variability in heart structure partly explains the ongoing debate on myocardial architecture. These developments are used to construct an integrated description of cardiac structure unifying features of fibre, sheet and band architecture that provides a basis for (i) explaining cardiac electromechanics, (ii) computational simulations of cardiac physiology and (iii) designing interventions.