A framework for myoarchitecture analysis of high resolution cardiac MRI and comparison with diffusion tensor MRI

Assessing Myocardial Architecture: The Challenges and Controversies

In recent decades, investigators have strived to describe and quantify the orientation of the cardiac myocytes in an attempt to classify their arrangement in healthy and diseased hearts. There are, however, striking differences between the investigations from both a technical and methodological standpoint, thus limiting their comparability and impeding the drawing of appropriate physiological conclusions from the structural assessments. This review aims to elucidate these differences, and to propose guidance to establish methodological consensus in the field. The review outlines the theory behind myocyte orientation analysis, and importantly has identified pronounced differences in the definitions of otherwise widely accepted concepts of myocytic orientation. Based on the findings, recommendations are made for the future design of studies in the field of myocardial morphology. It is emphasised that projection of myocyte orientations, before quantification of their angulation, introduces considerable bias, and that angles should be assessed relative to the epicardial curvature. The transmural orientation of the cardiomyocytes should also not be neglected, as it is an important determinant of cardiac function. Finally, there is considerable disagreement in the literature as to how the orientation of myocardial aggregates should be assessed, but to do so in a mathematically meaningful way, the normal vector of the aggregate plane should be utilised.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

Electrical wave propagation in an anisotropic model of the left ventricle based on analytical description of cardiac architecture

We develop a numerical approach based on our recent analytical model of fiber structure in the left ventricle of the human heart. A special curvilinear coordinate system is proposed to analytically include realistic ventricular shape and myofiber directions. With this anatomical model, electrophysiological simulations can be performed on a rectangular coordinate grid. We apply our method to study the effect of fiber rotation and electrical anisotropy of cardiac tissue (i.e., the ratio of the conductivity coefficients along and across the myocardial fibers) on wave propagation using the ten Tusscher–Panfilov (2006) ionic model for human ventricular cells. We show that fiber rotation increases the speed of cardiac activation and attenuates the effects of anisotropy. Our results show that the fiber rotation in the heart is an important factor underlying cardiac excitation. We also study scroll wave dynamics in our model and show the drift of a scroll wave filament whose velocity depends non-monotonically on the fiber rotation angle; the period of scroll wave rotation decreases with an increase of the fiber rotation angle; an increase in anisotropy may cause the breakup of a scroll wave, similar to the mother rotor mechanism of ventricular fibrillation.

DT-MRI measurement of myolaminar structure: accuracy and sensitivity to time post-fixation, b-value and number of directions

DT-MRI has been widely used to quantify myocardial fiber and laminar orientations. These structural orientations influence both the spread of excitation and the reorganization of the myocardium during contraction and are altered in disease states. Studies have sought to validate DT-MRI but questions remain about the accuracy of the method and its sensitivity to the time post-fixation and imaging parameters, including b-value, number of diffusion directions and image voxel size. The advent of high-spatial resolution ex vivo MRI and structure tensor (ST) analysis provides a means of direct validation of DT-MRI and assessment of sensitivity to the b-value, the number of diffusion directions and the image voxel size. We find that, with the fixation method we used, structure does not change with time (up to 72 hours). We show that DT-MRI and ST/HR-MRI are markedly similar measures of fiber orientation but DT-MRI and ST are much less similar measures of laminar orientation. DT-MRI performance is not sensitive to the number of directions, with similar structural orientations measured with 6 or 12 directions. Likewise, DT-MRI performance is generally insensitive to b-value, but laminar measurement is moderately more accurate at b = 500 than for higher b-values.

Mathematical model of the anatomy and fibre orientation field of the left ventricle of the heart

Background

One of the main factors affecting propagation of electrical waves and contraction in ventricles of the heart is anisotropy of cardiac tissue. Anisotropy is determined by orientation of myocardial fibres. Determining fibre orientation field and shape of the heart is important for anatomically accurate modelling of electrical and mechanical function of the heart. The aim of this paper is to introduce a theoretical rule-based model for anatomy and fibre orientation of the left ventricle (LV) of the heart and to compare it with experimental data. We suggest explicit analytical formulae that allow us to obtain the left ventricle form and its fibre direction field. The ventricle band concept of cardiac architecture given by Torrent-Guasp is chosen as the model postulate.

Methods

In our approach, anisotropy of the heart is derived from some general principles. The LV is considered as a set of identical spiral surfaces, each of which can be produced from the other by rotation around one vertical axis. Each spiral surface is filled with non-intersecting curves which represent myocardial fibres.

For model verification, we use experimental data on fibre orientation in human and canine hearts.

Results

LV shape and anisotropy are represented by explicit analytical expressions in a curvilinear 3-D coordinate system. The derived fibre orientation field shows good qualitative agreement with experimental data. The model reveals the most thorough quantitative simulation of fibre angles at the LV middle zone.

Conclusions

Our analysis shows that the band concept can generate realistic anisotropy of the LV. Our model shows good qualitative agreement between the simulated fibre orientation field and the experimental data on LV anisotropy, and the model can be used for various numerical simulations to study the effects of anisotropy on cardiac excitation and mechanical function.