- Open Access
In vivo cardiovascular magnetic resonance diffusion tensor imaging shows evidence of abnormal myocardial laminar orientations and mobility in hypertrophic cardiomyopathy
Journal of Cardiovascular Magnetic Resonance volume 16, Article number: 87 (2014)
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).
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.
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°).
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.
The compact myocardium of the left ventricle (LV) of humans and other mammals shows remarkable microstructural adaptation for global contractile function. Cardiomyocytes of the subendocardial and subepicardial layers contract almost orthogonally to each other, each angled obliquely, at helix-angles (HA) up to +60° and -60° relative to the circumferential alignments of myocytes of the mid layer ,. This means that myocytes contract in one layer in the same direction as the simultaneous counter-thickening of the orthogonally orientated myocytes of another layer. However, the potential conflict between these strains is resolved and enabled to contribute to wall thickening by the small scale laminar structures called sheetlets and shear layers (Table1) ,.
Populations of sheetlets are thought to collectively swivel around angles of approximately 45° to the local wall tangent plane, their transmural angles being smaller in diastole and larger in systole (Table1) . This allows the systolic shortening across the myocytes of one wall layer, imposed by the contraction of differently aligned myocytes of other wall layers, to be accommodated and translated to enhanced wall thickening . In this way, wall thickening can be attributed to sheetlet reorientation with shear layer slippage, in addition to slight systolic increase of myocyte diameter . Alternating, herringbone-like patterns of obliquely opposed sub-populations of sheetlets have been demonstrated histologically in ex vivo studies of several regions of the LVs of different mammalian species ,,,.
Ex vivo cardiac Diffusion Tensor Imaging (cDTI), based on measurements of diffusion distances of water molecules in myocardium, has been found capable of interrogating myocyte and sheetlet orientations ,,-. From the diffusion tensor of each voxel it is possible to extract a set of three eigenvectors. The direction of the principal eigenvector (E1, associated with the largest eigenvalue) has been shown to coincide with the orientations of the myocytes (frequently described as the fibre direction), when compared to histological results ,. The angle between E1 and the local radial direction is commonly referred to as the helix-angle. The secondary and tertiary eigenvectors (E2 and E3) are perpendicular to each other and E1. It has been argued that E2 and E3 align with the more and the less dominant local sheetlet populations respectively ,-. The direction and magnitude of E3 then indicates the orientation and prevalence of a secondary sheetlet sub population . However this can only, by definition, apply where two intravoxel sheetlet subpopulations are themselves perpendicular to each other. This ought not to be assumed, especially not at end systole and late diastole when the angles between adjoining subpopulations may be acute or obtuse ,,.
The integrity and mobility of the laminar microstructures is thought to be fundamental for LV function and could be altered in cardiovascular disease. Hypertrophic Cardiomyopathy (HCM) is an autosomal dominant genetic disease, which can be caused by mutations in the genes responsible for encoding proteins of the cardiac sarcomere ,. It is characterized by increased left ventricle (LV) wall thickness in the absence of increased LV afterload and is associated with diastolic dysfunction ,. Histological findings include myocyte hypertrophy and disarray, with an increase of interstitial connective tissue -. However, studies of myocardial strain in HCM have shown regional contractile heterogeneity that could not be fully explained by the distribution of regional fibrosis .
To date, most cDTI studies have been ex vivo, and evidence for the phasic reorientations of laminar structures has been provided by static, explanted hearts imaged separately in contracted and relaxed states ,. An exception was work by Dou et al. where 3D diffusion data was acquired in vivo at six different systolic trigger times .
We have recently described the development of a novel in vivo pulse sequence for cDTI which draws on technical improvements and also the benefits of imaging at 3 T. Measures derived from this technique were recently shown to be reproducible in both controls and in HCM subjects ,. In this study we used this in vivo sequence to analyze the cross-myocyte components of diffusion, which we hypothesized would provide information on mean intravoxel sheetlet orientations and their changes between end-systole and late-diastole in controls and show abnormalities in patients with HCM.
This study was approved by the National Research Ethics Service. Eleven healthy controls and 11 patients with HCM were recruited and gave written informed consent. Data in systole from 6 patients was previously reported for reproducibility . Imaging was performed using a 3 T scanner (Skyra, Siemens, Erlangen, Germany) with an 18 element anterior matrix coil and 8-12 elements of a matrix spine coil. After localization of the LV short and long axes, a retro-gated balanced Steady State Free Precession (bSSFP) cine sequence with a 40 ms temporal resolution was used to find the timing and duration of the subject specific end systolic and diastolic pauses, taking care to avoid the atrial kick in diastole.
In order to ensure the same myocardial region was imaged with cDTI in systole and diastole, a breath-hold spoiled gradient echo (GRE) sequence with a spatial modulation of magnetization (SPAMM) tagging pre-pulse and a 45 ms temporal resolution was used to acquire four-chamber and two-chamber views. The linear tags were separated by 16 mm and were perpendicular to the long-axis with an acquired spatial resolution of 2.1 × 1.7 mm in-plane and slice thickness of 6 mm. The displacement of the linear tag closest to the central mid-ventricular slice was manually tracked from the systolic to the diastolic phase. The positioning of the cDTI central slice was manually adjusted between the systolic and the diastolic acquisition accordingly, with the inter-slice distance kept constant.
The in vivo cardiac diffusion weighted stimulated echo acquisition mode (STEAM) single shot echo planar imaging (EPI) sequence used has been previously described . In short, this sequence runs over 2 heart beats and assumes that myocardium returns to the same position at the same encoding times in consecutive cycles, in this work, end systole and the late diastolic pause (diastasis). In order to minimize the length of the single shot EPI readout, parallel imaging and zonal excitation were implemented. Data were acquired in a series of breath-holds to minimize the effects of respiratory motion. Each breath-hold had a duration of 18 heartbeats. The following sequence parameters were used: reference plus 6 diffusion encoding directions at b =350 s/mm2 (at a heart rate of 60 beats per minute), fat saturation, TR = 2 RR intervals = 2000 ms (at a heart rate of 60 beats per minute), TE = 23 ms, BW = 2442Hz/pixel, GRAPPA  parallel imaging acceleration factor of 2, EPI echo train length = 22-28 readouts, depending on field of view, EPI echo-train readout duration = 11-14 ms, in plane acquired spatial resolution = 2.7 × 2.7 mm2, interpolated to 1.35 × 1.35 mm2 (k-space zero-filling by a factor of two), field of view 360 × 123-157 mm2 and a minimum of 8 averages (one average per breathold). Three slices with 8 mm thickness were acquired in the mid ventricular short axis with an inter-slice gap of 4 mm. No partial Fourier was used. The sequence trigger time was set to acquire cDTI data during the systolic and diastolic rest periods. Localized first and second-order shimming and frequency adjustment were performed over the extent of the LV within the imaging planes. It should be noted that the sequence runs a spoiler gradient in place of the diffusion encoding gradients for the reference (b0) acquisitions, which introduces a small but not negligible diffusion weighting to the "b0" images. This is taken into account during the tensor calculation.
Additional breath-hold CMR was performed to assess cardiac morphology and function in patients and controls. Retrospectively gated bSSFP cine sequences were acquired in 3 long-axis planes, followed by a contiguous stack of short-axis slices from the atrioventricular ring to the apex . Late gadolinium enhancement (LGE) imaging for the assessment of regional fibrosis was also performed, using an inversion recovery prepared spoiled GRE.
Ventricular volumes, function, mass, and ejection fraction was measured using a semi-automated threshold-based technique (CMRtools, Cardiovascular Imaging Solutions, London, UK), with indexing to body surface area . The LV epicardial and endocardial borders were defined by manual planimetry excluding papillary muscles. Each short-axis slice was divided into 12 segments. The presence of LGE was recorded, and when present its extent was quantified using CMR42 software (Circle Cardiovascular Imaging, Calgary, Canada) with a full-width half-maximum threshold and expressed as percent of LV mass.
Diffusion tensor analysis
All quantitative tensor information was post-processed with custom software, written in-house using MATLAB (Mathworks, Massachusetts, USA). In a first step, all the raw diffusion-weighted images were analyzed visually to reject frames corrupted by breathing motion or abnormal RR intervals. A minimum of two averages per direction, after frame rejection, was set as a requirement for inclusion in the study. The accepted frames were then co-registered to correct for intra-subject variations in breath-hold position (in-plane only) with a cross-correlation based algorithm for rigid displacement . The subject's RR-intervals during diffusion acquisition were then calculated from the acquisition times recorded in the DICOM files, and assumed to be half the time between frames (each frame runs over two heart-beats). This information was used to adjust the b-value for each image. A rank-2 diffusion tensor, and its eigensystem (set of three eigenvalues λ1, λ2, and λ3; and the respective eigenvectors: E1, E2, and E3) were then calculated for each voxel of the diffusion dataset. The eigenvalues and respective eigenvectors were sorted in a descending order according to the magnitude of the eigenvalues. A small percentage (3.2%) of the eigenvalues was negative due to motion, noise or mis-registration artifacts. This violates the assumption of a positive definite diffusion tensor, and therefore these values were averaged to the mean of neighbouring voxels.
Two diffusion tensor measurements were calculated from the eigensystem: Helix-angle (HA), and secondary eigenvector angle (E2A) ,. HA and E2A represent directional information from the first (E1) and second eigenvectors (E2) respectively. These were calculated by identifying the 3 orthogonal cardiac coordinate directions of each myocardial voxel: longitudinal, which is parallel to the left-ventricular long-axis and pointing towards the base; circumferential, which is tangential to the wall with a counter-clockwise direction when viewed from base to apex; and radial, which is given by the cross product of the previous two and pointing outward (Figure1). The primary eigenvector E1 was then projected radially to the local wall tangent plane. HA was defined as the angle in this plane between the E1 projection and the circumferential direction in the range -90 to 90 degrees, being positive (right-handed helix) if rotated counter-clockwise from the circumferential as viewed from the outside, and negative (left-handed helix) if rotated clockwise. The cross-myocyte plane, perpendicular to E1, was then calculated for each voxel. Then E2 was projected onto this plane and E2A calculated in it, between E2 and the illustrated cross-myocyte direction. This angle was measured in the range [-90, 90], being positive if rotated clockwise from the cross-myocyte direction when viewed in the more circumferential direction, and negative if rotated counter-clockwise (Figure1).
For 3D visualization and comparison of E1 distributions, principal eigenvector tractography was performed, color coded according to helix-angle (Figure2A,B,E,F). It must be emphasised that the obtained tracts are only used to aid on the visualisation of E1 orientation. No further analysis was performed on the tractograms, considering the limited number of slices acquired.
To display mean intra voxel sheetlet angulation in the septal and anterior wall regions, the full diffusion tensor was represented by superquadric glyphs, and superimposed with long thin white cylinders, which represent the orientation of E2. The tensors were color coded according to their absolute E2A value (Figure2C,D,G,H). 3D visualization of tensor information and tractography was performed using Python and PARAVIEW (Kitware, NM, USA).
Quantification and statistical analysis
The measured HA and E2A values were binned and plotted as histograms for all subjects in order to visualize the angle distribution. For further analysis of the effects of hypertrophy, E2A data was presented in a scatter plot. E2 mobility was defined as the change of mean absolute E2A values from diastole to systole. E2 mobility uses the absolute value of the angle so polarity is ignored to provide a measure of change of angulation relative to the local wall tangent plane from diastole to systole. Global E2A results included all LV segments in all 3 slices. For further insights in the HCM group, myocardial segments were categorized by consensus by 2 observers (12 segments per short-axis slice): segments with hypertrophy and LGE (H+LG+); segments with hypertrophy but no LGE (H+LG-); and segments with no hypertrophy or LGE (H-LG-). Segments were included in the H+ or LG+ categories even if hypertrophy or late-enhancement were only partially included. Mean absolute E2A values were compared between groups with the Mann-Whitney-Wilcoxon test. Statistical significance was considered present when p < 0.05 after Bonferroni correction for the number of tests performed.
The baseline demographics and clinical characteristics of control and patient groups showed no significant differences (Table2). The time spent in diffusion acquisition was 42 ± 7 min (mean ± SD) for controls, and 46 ± 7 min for HCM patients. The number of required breath-holds was 57 ± 4 for controls, and 63 ± 9 for HCM patients. The SD recorded on the RR interval, during diffusion acquisition, were 38 ± 9 ms for controls, and 41 ± 20 ms for HCM patients.
Average LV mass, wall thickness and ejection fraction were greater in the HCM group and end-systolic volume was smaller. The ventricular morphology, risk factors for sudden death and LGE patterns for the HCM patients are summarized in Table3. Ten out of 11 patients had asymmetrical septal hypertrophy and 1 patient had anterior wall hypertrophy. The median number of risk factors for sudden death was 0 (range 0-2). Six patients had LV outflow tract obstruction. All patients had regional fibrosis on LGE imaging usually affecting the septum and representing 9.1% of LV mass.
Diffusion tensor visualization
Figure2 shows 3-dimensional graphic representations of the orientations of E1 (mean intravoxel helix-angle) and absolute E2 (greater cross-myocyte direction of diffusion) rendered for 1 control and 1 HCM patient at both cardiac phases. The expected helix-angle distributions are observed, transitioning inwards transmurally from left-handed, through circumferential to right handed, without obvious changes of layer-wise orientation between phases, or between the control and the patient. In contrast, glyph representations of the diffusion tensor, and in particular orientation of E2, show unmistakable reorientation, reflecting phasic changes in absolute E2A. E2 tilts towards parallel to the wall-plane in diastole, and towards perpendicular in systole, particularly through mid-wall and subendocardial voxels, consistent with the predicted sheetlet reorientations during the cardiac cycle. In the HCM patient however, most obviously in the most hypertrophic wall region, the phasic changes are less pronounced, with the glyphs remaining relatively perpendicular to the local wall plane in diastole as well as in systole.
Diffusion tensor quantification
cDTI was successfully performed on all subjects. Only 3.2% of all eigenvalues were spuriously negative, for which we assigned the positive mean values of adjoining voxels. In total the number of myocardial voxels analyzed for each group was: controls-systole 28,836 voxels; controls-diastole 26,034 voxels; HCM-systole 45,791 voxels; HCM-diastole 44,661 voxels. The greater number of voxels in HCM patients occurred because of LV hypertrophy. The measured E2 angles were found to be meaningful (more details in Additional file 1: Appendix 1). An example of a mid-slice averaged magnitude image, together with the respective HA, E2A and absolute E2A maps at both cardiac phases is shown in Figure3 for a control and an HCM patient.
The distribution of HA is shown in Figure4. In general, HA maintained the same distribution between cardiac phases and also when comparing controls and HCM patients, all showing greatest frequency around the circumferential orientation (i.e. angles around 0 degrees). The main difference was in the controls, where diastolic distribution shows a narrower peak when compared to the systolic distribution. Contrasting with HA results, the E2A values did not follow an obvious spatial pattern within the myocardium, but showed substantial changes between the 2 cardiac phases, especially for controls (Figure3). A histogram of the measured myocardial E2A values for both groups is shown in Figure5. Systolic values peak at ±90 degrees for both groups, implying sheetlet orientations tilted in towards perpendicular to local wall planes, generally even more so in HCM than controls. However marked differences between the 2 groups were seen in diastole. In controls the diastolic angles peak at approximately 0 degrees, implying mean intravoxel sheetlet orientations towards parallel to local wall plane. However, in HCM, the diastolic distribution curves were flattened, showing a more homogeneous distribution across the range of angles without the tendency towards wall-parallel orientation seen in controls. The histograms were approximately symmetric, except for a slight predominance of negative systolic angles in controls.
Additional analysis of the E2A data is shown in the E2 mobility plots of Figure6. Figure6A shows the global mean absolute values for both cohorts at the 2 cardiac phases. There was a significant difference in global E2A in diastole between HCM and controls (p < 0.001). In systole, there was also a statistically significant difference between the 2 groups of lesser magnitude (p =0.026). Figure6B and Table4 show further analysis of E2 mobility in the HCM cohort, with the myocardial absolute E2A values divided into 3 categories: segments with hypertrophy and LGE (H+LG+); segments with hypertrophy but no LGE (H+LG-); and regions with no hypertrophy or LGE (H-LG-). There was a significant difference in E2A at diastole between non-hypertrophic segments (H-LG-) and the other two groups H+LG+ and H+LG- (p =0.0028 and p =0.0022 respectively, both after Bonferroni correction for 2 tests). There was also a significant difference in E2A in systole between hypertrophic regions and controls for H+LG+ segments and H+LG- segments (p =0.0060 and p =0.0030 respectively both after Bonferroni correction for 3 tests).
Impairment of E2 mobility in hypertrophic segments is shown in Figure7, where slices from 2 HCM hearts with asymmetric hypertrophy are compared with the corresponding slice from a control. The non-hypertrophic lateral wall in both HCM hearts approaches the absolute E2A angles measured in the control heart, while differences between the hypertrophic regions and the control values are clearly apparent in diastole. There was no obvious effect from the presence or absence of LGE on the E2A values in the hypertrophied anteroseptal wall, in accord with the results shown in Figure6B.
We believe this is the first report of in vivo systolic and diastolic measurements of the principal (E1) and second (E2) directional components of intramyocardial diffusion in HCM patients as measured by cDTI. The LV helix-angle (HA) distribution (Figures3 and 4), which is individually and collectively consistent with previous ex vivo studies, showed only slight changes between systole and diastole ,, or between control and HCM hearts . In contrast, we found marked changes between systolic and diastolic phases in the principal cross-myocyte direction of diffusion E2A, and these are consistent with previously reported changes of transmural angulation of the laminar sheetlet and shear layer structures ,,. We assume, as others have, that the measured cross-myocyte diffusions depend on the relative freedom of molecules to diffuse, particularly through the extracellular fluid of shear layers between sheetlets. The cyclic changes of E2 angle (E2A) measured in volunteers were in keeping with previously reported changes of transmural orientation of sheetlets and shear layers (Figure6). However, interpretation of them is complicated by the need to consider previously postulated effects of myocardial strain on diffusion. As microstructures and micro-strains are likely to be fundamental to interpretation of cDTI, we include histological sections from a pig heart prepared in our own institution (Figure8) and a schematic illustration (Figure9). Together, these are intended to help conceptualise the possible effects of myocardial microstructures and their expected patterns of deformation on diffusion.
The issue of strain
The microstructure of myocardium is not isotropic. However, with a view to its potential relevance to cDTI, Reese et al. considered the effects of cyclically applied strain on aqueous diffusion through a structurally isotropic material  They used an elastic gel phantom, periodically indented by a hemi-spherical driver. This caused radially directed compression of the gel with counter-stretching in the orthogonal circumferential directions, analogous with the radial thinning and bidirectional elongation of LV walls in diastole. On release, elastic restoration of the gel was associated with radially directed thickening and bidirectional shortening, analogous with systolic deformations of the LV wall, as imaged macroscopically. As predicted by the underlying theory of diffusion through an isotropic material during deformation, DTI in the `systolic' phase showed anisotropy, with increase of the radial and reduction of the circumferential components of diffusion. We have schematically illustrated this effect in Figure9A. The converse was found by DTI in the diastolic phase. Given normal, healthy myocardial strain, comparable strain-induced directions and amounts of anisotropy might be thought to account for the changes of E2A that we found in healthy volunteers. Strain might be taken as the underlying explanation, except for the following crucial considerations:
Firstly, the myocardium is not structurally isotropic, but orthotropic in most LV regions , meaning that it has different constituent properties in three orthogonal directions: the myocyte, sheetlet and sheetlet normal directions. Diffusion taking place through it is likely to be constrained by microstructural boundaries in ways that modify the postulated effects of strain on diffusion, particularly for the inherent long mixing-times (one heart-beat) of the STEAM sequence. Furthermore, the constituent microstructures undergo their own distinct small-scale deformations, including the shearing and cyclic reorientation of extracellular fluid in the interstitial shear layers. This is perhaps the most pertinent of several types of micro-deformation in contracting myocardium that may contribute to the measured anisotropies of diffusion. Others include the even smaller scale interdigitation of actin and myosin nanofilaments within sarcomeres, the lengthwise shortening and slight circumferential counter-thickening of whole myocytes and the reorientations of the sheetlets in which the myocytes are arrayed. These micro-deformations underlie but can, especially in the cross-myocyte directions, differ in magnitude and sign from deformations imaged macroscopically. The marked systolic thickening of myocardium is known to be caused more by laminar reorientations than by the thickening of individual myocytes . This amount of radial thickening, combined with systolic shortening in cross-myocyte as well as myocyte directions in wall tangent planes, is only achievable through the presence and mobility of laminar microstructures.
Secondly, the effects of strain cannot account for the exaggerated systolic anisotropy of diffusion that we found in the hypertrophied regions of HCM hearts. As is well known in HCM, strain is reduced, especially in hypertrophic regions . However, in systole, we measured above normal E2As implying marked predominance of radial over tangential components of diffusion. In diastole, on the other hand, in spite of the fact that even the hypertrophic myocardium had undergone radial thickening with slight longitudinal and circumferential shortening during systole, albeit of limited amount, there remained a predominance of radial over tangential components of cross-myocyte diffusion, with E2As greater than 45°. From the postulated effects of strain alone, the opposite would have been the case. To appreciate an alternative, more plausible interpretation of the HCM DTI results, the function of laminar microstructures should be considered further.
Sheetlet and shear layer function
The contraction of myocytes of one myocardial layer, for example the sub-epicardial layer, imposes cross-myocyte shortening on the simultaneously contracting myocytes of a different layer, in this case the subendocardial. As mentioned in the introduction, it is this potential conflict that requires the presence and function of laminae: the sheetlets and shear layers not only offer compliance that allows cross-myocyte shortening in spite of their individual thickening, but the changing angulations of laminae, like the folds of a compressed concertina , contribute to wall thickening. The reorientations of laminae of a particular region are driven not by their own force, but rather by the contraction of deeper and/or more superficial myocytes of the same region. Contraction of myocytes through all depths of any given segment causes laminae to adopt steeper, through-wall obliquity, with increase of E2A values and of wall thickness. This underlies our interpretation of the elevated E2As found in the hypertrophic regions of HCM patients, in diastole as well as systole.
We had hypothesized that sheetlet mobility in HCM would be impaired, which we found. Our results also show that the diastolic orientations appear markedly abnormal, with the E2A of HCM hearts retaining relatively systolic conformation in diastole. This abnormality was associated with increased wall thickness, independent of the presence of myocardial LGE (Figure6B). This suggests that the wall's state of contraction or hypertrophy (not presence of LGE) is the dominant determinant of E2A in HCM, a finding that is consistent with findings by Aletras et al., that regional contractile heterogeneity in HCM could not be explained by the distribution of fibrosis . It must however be noted that although capable of detecting regional fibrosis, LGE imaging is not suitable to image diffuse fibrosis . T1-mapping is an emerging technique more suitable to detect and quantify diffuse fibrosis .
Although less marked than the diastolic abnormality, we also found exaggerated E2A values in systole, consistent with a hypercontracted state.
Our in vivo data, which suggests relatively systolic sheetlet orientations and impaired relaxation is consistent with current hypotheses of the pathophysiological consequences of HCM-causing mutations from biochemical studies and transgenic models. A common property of HCM-causing mutations is increased myofilament or myosin ATPase calcium sensitivity -. These mutations increase cardiomyocyte tension for a given cytoplasmic calcium level, particularly in the diastolic range, slowing relaxation kinetics, limiting diastolic relaxation, and under some circumstances increasing the tension generated during systole. Exaggeration and persistence of systolic tone, and delayed relaxation, have implications for myocardial perfusion given that coronary flow predominantly occurs in diastole enabled by myocardial relaxation , and can be compromised in HCM .
The improved understanding of pathophysiological mechanisms has potential implications for novel approaches to therapeutic intervention ,. It may be feasible to use cDTI to identify abnormal myocardium, before the onset of fibrosis, and monitor its response to the effects a drug, known or novel, initially in animal experiments and, if apparently successful, ultimately in patients.
No macroscopic evidence of disarray in the HA maps of HCM patients was found. These findings contrast with the work of Tseng et al., which found more longitudinally oriented HAs in the hypertrophied epicardial septum when compared to the normal heart . We consider there may be 2 factors, which have contributed to this difference. In some of the HCM subjects, longitudinal oriented HAs that were considered to be hypertrophy from the septo-marginal trabeculations of the right-ventricle were excluded from analysis of the left ventricle. Additionally our cohort of HCM subjects was relatively old with relatively few risk factors and might therefore have had limited myocardial disarray.
The spatial-resolution achievable in vivo cDTI is limited by acquisition parameters chosen for suitability in patients with limited breath-holding ability. In controls, diastolic wall thinning reduces effective resolution of differences between layers, although this was less problematic in HCM. A zero-filling interpolation algorithm was used. This has the potential to increase spatial resolution by using the corners of k-space, although the spatial resolution improvement becomes angle dependent . The impact of this interpolation in measuring mean E2 angles was tested in 18 datasets acquired for another study; results showed minimal differences of less than 3 degrees when compared to no interpolation. A maximum difference lower than 8 degrees was found for the helical-angle. We therefore expect this interpolation to have minimal impact in the diffusion tensor calculation.
There are a number of limitations to this study. As discussed above, the sensitivity to myocardial strain of the cDTI technique remains unclear. It is likely to differ from that of a homogenous, elastic gel and further work is required to investigate whether or not it contributes to the results reported.
The relation between E2A and the orientations of sheetlets within a voxel is unlikely to be straightforward. Thousands of cardiomyocytes arrayed in hundreds of sheetlets, with potentially two or more opposing populations may occupy a single 2.7 × 2.7 × 8 mm3 imaging voxel ,. The sheetlets are expected to change orientation through the cycle relative to one another as well as to the local wall plane. Angles between adjacent sheetlet populations of the mid layer might, for example, if viewed like the folds of an accordion, be obtuse when extended longitudinally in diastole, but acute when brought together in systole. Were a voxel to contain 2 equal opposing populations, this might result in very small, almost wall parallel, E2As in diastole even though individual sheetlet populations remained oblique. Interestingly, however, our analysis of E2A supports the predominance of one orientation of sheetlet population in most voxels. We have not yet identified any consistent pattern of distribution of the `positive' and `negative' slopes of these between individuals. The different and unknown orientations of 2 or more sheetlet populations in a voxel cannot be elucidated using a single tensor. cDTI techniques with higher angular diffusion resolution (multi-tensor, diffusion spectrum imaging, and q-ball) would be needed for this , although the acquisition times required would currently be prohibitive in a clinical setting.
Myocardial coverage was limited to 3 slices. The lack of whole heart coverage made it challenging to define the longitudinal wall tangent direction. In this work, it was taken as perpendicular to the short-axis plane, which is a better approximation for the basal and mid than the more apical slice. Intra-subject variations in breath-hold are corrected with a rigid in-plane registration algorithm; through-plane variation remains uncorrected. The SPAMM imaging was used to track manually the imaged slices from systole to diastole, with the objective of reducing differences in the imaged volume between the two phases. This method is very approximate and it does not take into account the complex cardiac deformation present entirely, particularly in the HCM cohort. Nevertheless we believe this method to be sufficient for the DTI metrics shown, as we are only showing differences in global mean values.
The diffusion-encoding scheme used had only 6 directions. Recent DTI work in the brain showed comparable Mean Diffusivity (MD) and Fractional Anisotropy (FA) between 6 directions and 30 or 60 direction schemes (acquisition time kept constant) . DTI numerical simulations also show comparable tensor accuracy when trading directions for averages . It remains unknown if 6 direction schemes are comparable to higher order schemes for cDTI. A recent study from our group shows that both HA and E2A maps become smoother with higher b-values (>750 s/mm2). However E2A global mean values do not change substantially .
The averaging of non-positive definite tensors was performed for the eigenvalues directly. We appreciate that this is not a robust method, it maintains the neighbouring FA and MD, but it lacks sophisticated tensor interpolation. Nevertheless it affects only a small percentage (3.2%) of the eigenvalues measured, so we expect the impact to be small for the global means of E2A.
In vivo cDTI in normals delivered apparently meaningful measurements of cross-myocyte diffusion, in keeping with previously reported myocardial sheetlet and shear layer reorientations from more wall-parallel in diastole to more wall-perpendicular in systole. In HCM, impaired diastolic reorientation of sheetlet populations was identified which was related to increased wall thickness and not LGE. The apparent persistence in diastole of a relatively systolic conformation in HCM may represent a novel in vivo human pathophysiological finding and would be consistent with biochemical studies and transgenic models that point to increased calcium sensitivity and impaired diastolic relaxation associated with known causative mutations.
Although the extent of the strain contribution to these results is currently uncertain, we believe that the measurements described in this study may represent a step towards in vivo investigation of myocardial structural dynamics.
The potential of this novel non-invasive approach to assess HCM pathophysiology in-vivo at the tissue level, combined with advances in the understanding of underlying molecular mechanisms, has potential for the development and evaluation of novel therapeutic approaches for patients with HCM.
Joint first authors: Pedro F Ferreira, Philip J Kilner and Laura-Ann McGill.
Joint senior authors: David N Firmin and Dudley J Pennell.
Pettigrew JB: On the arrangement of the muscular fibres in the ventricles of the vertebrate, with physiological remarks. Philos Trans. 1864, 154: 445-500. 10.1098/rstl.1864.0014.
Smerup M, Nielsen E, Agger P, Frandsen J, Vestergaard-Poulsen P, Andersen J, Nyengaard J, Pedersen M, Ringgaard S, Hjortdal V, Lunkenheimer PP, Anderson RH: The three-dimensional arrangement of the myocytes aggregated together within the mammalian ventricular myocardium. Anat Rec. 2009, 292: 1-11. 10.1002/ar.20798.
Sands GB, Gerneke DA, Hooks DA, Green CR, Smaill BH, LeGrice IJ: Automated imaging of extended tissue volumes using confocal microscopy. Microsc Res Tech. 2005, 67: 227-39. 10.1002/jemt.20200.
Hales PW, Schneider JE, Burton RA, Wright BJ, Bollensdorff C, Kohl P: Histo-anatomical structure of the living isolated rat heart in two contraction states assessed by diffusion tensor MRI. Prog Biophys Mol Biol. 2012, 110: 319-30. 10.1016/j.pbiomolbio.2012.07.014.
Harrington KB, Rodriguez F, Cheng A, Langer F, Ashikaga H, Daughters GT, Criscione JC, Ingels NB, Miller DC: Direct measurement of transmural laminar architecture in the anterolateral wall of the ovine left ventricle: new implications for wall thickening mechanics. Am J Physiol Heart Circ Physiol. 2005, 288: H1324-30. 10.1152/ajpheart.00813.2004.
Spotnitz HM, Spotnitz WD, Cottrell TS, Spiro D, Sonnenblick EH: Cellular basis for volume related wall thickness changes in the rat left ventricle. J Mol Cell Cardiol. 1974, 6: 317-31. 10.1016/0022-2828(74)90074-1.
LeGrice IJ, Takayama Y, Covell JW: Transverse shear along myocardial cleavage planes provides a mechanism for normal systolic wall thickening. Circ Res. 1995, 77: 182-93. 10.1161/01.RES.77.1.182.
Kung GL, Nguyen TC, Itoh A, Skare S, Ingels NB, Miller DC, Ennis DB: The presence of two local myocardial sheet populations confirmed by diffusion tensor MRI and histological validation. J Magn Reson Imaging. 2011, 34: 1080-91. 10.1002/jmri.22725.
Gilbert SH, Benoist D, Benson AP, White E, Tanner SF, Holden AV, Dobrzynski H, Bernus O, Radjenovic A: Visualization and quantification of whole rat heart laminar structure using high- spatial resolution contrast-enhanced MRI. Am J Physiol Heart Circ Physiol. 2012, 302: H287-98. 10.1152/ajpheart.00824.2011.
Scollan DF, Holmes A, Winslow R, Forder J: Histological validation of myocardial microstructure obtained from diffusion tensor magnetic resonance imaging. Am J Physiol. 1998, 275 (6 Pt 2): H2308-18.
Scollan DF, Holmes A, Zhang J, Winslow RL: Reconstruction of cardiac ventricular geometry and fiber orientation using magnetic resonance imaging. Ann Biomed Eng. 2000, 28: 934-44. 10.1114/1.1312188.
Tseng WY, Wedeen VJ, Reese TG, Smith RN, Halpern EF: Diffusion tensor MRI of myocardial fibers and sheets: correspondence with visible cut-face texture. J Magn Reson Imaging. 2003, 17: 31-42. 10.1002/jmri.10223.
Maron MS: Clinical utility of cardiovascular magnetic resonance in hypertrophic cardiomyopathy. J Cardiovasc Magn Reson. 2012, 14: 13-10.1186/1532-429X-14-13.
Maron BJ, Maron MS: Hypertrophic cardiomyopathy. Lancet. 2013, 381: 242-55. 10.1016/S0140-6736(12)60397-3.
Ho CY, Sweitzer NK, McDonough B, Maron BJ, Casey SA, Seidman JG, Seidman CE, Solomon SD: Assessment of diastolic function with doppler tissue imaging to predict genotype in preclinical hypertrophic cardiomyopathy. Circulation. 2002, 105: 2992-7. 10.1161/01.CIR.0000019070.70491.6D.
Motoyasu M, Kurita T, Onishi K, Uemura S, Tanigawa T, Okinaka T, Takeda K, Nakano T, Ito M, Sakuma H: Correlation between late gadolinium enhancement and diastolic function in hypertrophic cardiomyopathy assessed by magnetic resonance imaging. Circ J. 2008, 72: 378-83. 10.1253/circj.72.378.
Factor SM, Butany J, Sole MJ, Wigle ED, Williams WC, Rojkind M: Pathologic fibrosis and matrix connective tissue in the subaortic myocardium of patients with hypertrophic cardiomyopathy. J Am Coll Cardiol. 1991, 17: 1343-51. 10.1016/S0735-1097(10)80145-7.
Shirani J, Pick R, Roberts WC, Maron BJ: Morphology and significance of the left ventricular collagen network in young patients with hypertrophic cardiomyopathy and sudden cardiac death. J Am Coll Cardiol. 2000, 35: 36-44. 10.1016/S0735-1097(99)00492-1.
Hughes SE: The pathology of hypertrophic cardiomyopathy. Histopathology. 2004, 44: 412-27. 10.1111/j.1365-2559.2004.01835.x.
Aletras AH, Tilak GS, Hsu LY, Arai AE: Heterogeneity of intramural function in hypertrophic cardiomyopathy: mechanistic insights from MRI late gadolinium enhancement and high-resolution displacement encoding with stimulated echoes strain maps. Circ Cardiovasc Imaging. 2011, 4: 425-34. 10.1161/CIRCIMAGING.110.958751.
Chen J, Liu W, Zhang H, Lacy L, Yang X, Song SK, Wickline S, Yu X: Regional ventricular wall thickening reflects changes in cardiac fiber and sheet structure during contraction: quantification with diffusion tensor MRI. Am J Physiol Heart Circ Physiol. 2005, 289: H1898-907. 10.1152/ajpheart.00041.2005.
Dou J, Reese TG, Tseng WY, Wedeen VJ: Cardiac diffusion MRI without motion effects. Magn Reson Med. 2002, 48: 105-14. 10.1002/mrm.10188.
Nielles-Vallespin S, Mekkaoui C, Gatehouse P, Reese TG, Keegan J, Ferreira PF, Collins S, Speier P, Feiweier T, de Silva R, Jackowski MP, Pennell DJ, Sosnovik DE, Firmin D: In vivo diffusion tensor MRI of the human heart: reproducibility of breath-hold and navigator-based approaches. Magn Reson Med. 2012, 70: 454-65. 10.1002/mrm.24488.
McGill LA, Ismail TF, Nielles-Vallespin S, Ferreira P, Scott AD, Roughton M, Kilner PJ, Ho SY, McCarthy KP, Gatehouse PD, de Silva R, Speier P, Feiweier T, Mekkaoui C, Sosnovik DE, Prasad SK, Firmin DN, Pennell DJ: Reproducibility of in-vivo diffusion tensor cardiovascular magnetic resonance in hypertrophic cardiomyopathy. J Cardiovasc Magn Reson. 2012, 14: 86-10.1186/1532-429X-14-86.
Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, Kiefer WB, Haase A: Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn Reson Med. 2002, 47: 1202-10. 10.1002/mrm.10171.
Kramer CM, Barkhausen J, Flamm SD, Kim RJ, Nagel E: Standardized cardiovascular magnetic resonance imaging (CMR) protocols, society for cardiovascular magnetic resonance: board of trustees task force on standardized protocols. J Cardiovasc Magn Reson. 2008, 10: 35-10.1186/1532-429X-10-35.
Du Bois D, Du Bois EF: A formula to estimate the approximate surface area if height and weight be known. Nutrition. 1916, 5: 303-11.
Guizar-Sicairos M, Thurman ST, Fienup JR: Efficient subpixel image registration algorithms. Opt Lett. 2008, 33: 156-8. 10.1364/OL.33.000156.
Basser PJ: Inferring microstructural features and the physiological state of tissues from diffusion- weighted images. NMR Biomed. 1995, 8: 333-44. 10.1002/nbm.1940080707.
Stoeck CT, Toussaint N, Boesiger P, Kozerke S: Dual Heart-Phase Cardiac DTI Using Local-Look STEAM. 2012, ISMRM, Milan
Reese TG, Wedeen VJ, Weisskoff RM: Measuring diffusion in the presence of material strain. J Magn Reson B. 1996, 112: 253-8. 10.1006/jmrb.1996.0139.
Ellims AH, Iles LM, L-h L, Hare JL, Kaye DM, Taylor AJ: Diffuse myocardial fibrosis in hypertrophic cardiomyopathy can be identified by cardiovascular magnetic resonance, and is associated with left ventricular diastolic dysfunction. J Cardiovasc Magn Reson. 2012, 14: 76-10.1186/1532-429X-14-76.
Moon JC, Messroghli DR, Kellman P, Piechnik SK, Robson MD, Ugander M, Gatehouse PD, Arai AE, Friedrich MG, Neubauer S, Schulz-Menger J, Schelbert EB: Society for cardiovascular magnetic resonance imaging, and cardiovascular magnetic resonance working group of the European society of cardiology. Myocardial T1 mapping and extracellular volume quantification: a society for cardiovascular magnetic resonance (SCMR) and CMR working group of the European society of cardiology consensus statement. J Cardiovasc Magn Reson. 2013, 15: 92-10.1186/1532-429X-15-92.
Gomes AV, Harada K, Potter JD: A mutation in the n-terminus of troponin I that is associated with hypertrophic cardiomyopathy affects the Ca(2+)-sensitivity, phosphorylation kinetics and proteolytic susceptibility of troponin. J Mol Cell Cardiol. 2005, 39: 754-65. 10.1016/j.yjmcc.2005.05.013.
Baudenbacher F, Schober T, Pinto JR, Sidorov VY, Hilliard F, Solaro RJ, Potter JD, Knollmann BC: Myofilament Ca2+ sensitization causes susceptibility to cardiac arrhythmia in mice. J Clin Invest. 2008, 118: 3893-903.
Pinto JR, Parvatiyar MS, Jones MA, Liang J, Ackerman MJ, Potter JD: A functional and structural study of troponin c mutations related to hypertrophic cardiomyopathy. J Biol Chem. 2009, 284: 19090-100. 10.1074/jbc.M109.007021.
Willott RH, Gomes AV, Chang AN, Parvatiyar MS, Pinto JR, Potter JD: Mutations in troponin that cause HCM, DCM and RCM: what can we learn about thin filament function?. J Mol Cell Cardiol. 2010, 48: 882-92. 10.1016/j.yjmcc.2009.10.031.
Davies JE, Whinnett ZI, Francis DP, Manisty CH, Aguado-Sierra J, Willson K, Foale RA, Malik IS, Hughes AD, Parker KH, Mayet J: Evidence of a dominant backward-propagating "suction" wave responsible for diastolic coronary filling in humans, attenuated in left ventricular hypertrophy. Circulation. 2006, 113: 1768-78. 10.1161/CIRCULATIONAHA.105.603050.
Cecchi F, Olivotto I, Gistri R, Lorenzoni R, Chiriatti G, Camici PG: Coronary microvascular dysfunction and prognosis in hypertrophic cardiomyopathy. N Engl J Med. 2003, 349: 1027-35. 10.1056/NEJMoa025050.
Maron MS, Olivotto I, Maron BJ, Prasad SK, Cecchi F, Udelson JE, Camici PG: The case for myocardial ischemia in hypertrophic cardiomyopathy. J Am Coll Cardiol. 2009, 54: 866-75. 10.1016/j.jacc.2009.04.072.
Ashrafian H, McKenna WJ, Watkins H: Disease pathways and novel therapeutic targets in hypertrophic cardiomyopathy. Circ Res. 2011, 109: 86-96. 10.1161/CIRCRESAHA.111.242974.
Tseng W-Y, Dou J, Reese TG, Wedeen VJ: Imaging myocardial fiber disarray and intramural strain hypokinesis in hypertrophic cardiomyopathy with MRI. J Magn Reson Imaging. 2006, 23 (1): 1-8. 10.1002/jmri.20473.
Bernstein MA, Fain SB, Riederer SJ: Effect of windowing and zero-filled reconstruction of MRI data on spatial resolution and acquisition strategy. J Magn Reson Imaging. 2001, 14 (3): 270-80. 10.1002/jmri.1183.
del Monte F, O'Gara P, Poole-Wilson PA, Yacoub M, Harding SE: Cell geometry and contractile abnormalities of myocytes from failing human left ventricle. Cardiovasc Res. 1995, 30: 281-90. 10.1016/S0008-6363(95)00040-2.
Sosnovik DE, Wang R, Dai G, Reese TG, Wedeen VJ: Diffusion MR tractography of the heart. J Cardiovasc Magn Reson. 2009, 11: 47-10.1186/1532-429X-11-47.
Lebel C, Benner T, Beaulieu C: Six is enough? Comparison of diffusion parameters measured using six or more diffusion-encoding gradient directions with deterministic tractography. Magn Reson Med. 2012, 68: 474-83. 10.1002/mrm.23254.
Gahm JK, Kindlmann G, Ennis DB: The effects of noise over the complete space of diffusion tensor shape. Med Image Anal. 2014, 18: 197-210. 10.1016/j.media.2013.10.009.
Scott AD, Ferreira PF, Nielles-Vallespin S, Gatehouse P, McGill LA, Kilner P, Pennell DJ, Firmin DN: Optimal diffusion weighting for in vivo cardiac diffusion tensor imaging. Magn Reson Med. 2014,in press
The authors are grateful to the staff of the CMR Unit and the Cardiovascular Biomedical Research Unit, Royal Brompton Hospital for their support with this work.
The Guest Editor for this paper was Leon Axel, Professor of Radiology, Medicine, and Neuroscience and Physiology at the New York University School of Medicine. firstname.lastname@example.org.
This work was supported by the National Institute of Health Research Cardiovascular Biomedical Research Unit at the Royal Brompton Hospital and Imperial College, London.
Prof Dudley Pennell is a consultant to Siemens and a director of Cardiovascular Imaging Solutions. Royal Brompton Hospital has research collaboration agreements with Siemens AG Medical Solutions.
All authors contributed in the design, and intellectual conception. PF, PK, LAM, SNV and AS participated in the acquisition, analysis and interpretation of the data. All authors participated in revision, read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Ferreira, P.F., Kilner, P.J., McGill, LA. et al. In vivo cardiovascular magnetic resonance diffusion tensor imaging shows evidence of abnormal myocardial laminar orientations and mobility in hypertrophic cardiomyopathy. J Cardiovasc Magn Reson 16, 87 (2014). https://doi.org/10.1186/s12968-014-0087-8
- Diffusion tensor imaging
- Hypertrophic cardiomyopathy
- Cardiovascular magnetic resonance
- Myocardial architecture
- Laminar structure
- Sheet and shear layers
- Diastolic dysfunction