In vivo cardiovascular magnetic resonance diffusion tensor imaging shows evidence of abnormal myocardial laminar orientations and mobility in hypertrophic cardiomyopathy
© Ferreira et al.; licensee BioMed Central Ltd. 2014
Received: 3 March 2014
Accepted: 24 September 2014
Published: 12 November 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.
Glossary of terms used
Used here for a heart muscle cell, capable of contracting lengthwise, along its principal axis. We have not used `fiber', which we do not consider appropriate in relation to myocardial structure.
Directions perpendicular to the principal axis of local myocytes.
Laminar arrays of myocytes, as visualized in cross-myocyte histological sections of mammalian myocardium. We use `sheetlets' rather than `sheets' because each is typically only a few myocytes thick and is of limited individual extent before either branching or merging with adjoining sheetlets. Sheetlets are generally orientated oblique to the local wall tangent plane, and can occur in sub-populations with opposing orientations that give herringbone-like appearances in some cross-myocyte sections.
Collagen lined fissures or interstices between adjacent sheetlets which may be important in relation to the extents and preferential directions of aqueous diffusion through myocardium over the course of a heart cycle. The fissures are thought to allow adjacent sheetlets to slide relative to one another as they collectively swivel, increasing their obliquity (transmural angle) relative to the local wall tangent plane during contraction.
An inclusive adjective referring to myocardial `sheetlet and shear layer' structure.
The principal eigenvector of aqueous diffusion. Previous studies have shown its direction to be aligned with local myocytes.
E2 and E3:
Represent the greater and the least eigenvectors, respectively, of cross-myocyte diffusion in a voxel. They are orthogonal to E1 and to each other. Although E2 and E3 may represent the sheetlet and sheet-normal directions perpendicular to myocytes, it should be remembered that DTI can only measure the mean of all diffusions in any single voxel and that a voxel of 3×3×8 mm3 may include two or more sub-populations of sheetlets with different orientations.
Helix angle or E1 angle (E1A):
Angle of obliquity of E1 in the local wall tangent plane relative to the local circumferential direction.
Transmural angle or E2 angle (E2A):
Angle of through-wall obliquity of E2 in the local cross-myocyte plane relative to the line of intersection of this plane with the local wall tangent plane (see Figure1).
Diffusion that takes place with equal a freedom in all directions, i.e. spherically, on average, from any given starting point.
Diffusion that is constrained by anisotropic structures so that it proceeds preferentially along certain directions or planes.
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.
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.
Baseline patient characteristics
HCM Patients (n =11)
Controls (n =11)
62 ± 9
56 ± 10
Body surface area (m2)
1.89 ± 0.2
1.83 ± 0.2
Body mass index (kg/m2)
26.1 ± 3.6
24.2 ± 2.3
RR interval (ms)
995 ± 150
984 ± 108
Mid systolic pause (ms)
373 ± 70
353 ± 41
Mid diastolic pause (ms)
798 ± 96
731 ± 79
Distance of mid image slice from apex (measured in systole as a proportion of total ventricular length)
0.53 ± 0.06
0.53 ± 0.04
Indexed LV EDV (mL/m2)
65 ± 15
78 ± 16
Indexed LV ESV (mL/m2)
14 ± 5
23 ± 5
LV ejection fraction (%)
79 ± 6
71 ± 3
LV mass index (g/m2)
104 ± 41
67 ± 15
Maximum end-diastolic wall thickness (mm)
23 ± 4
9 ± 1
Morphology and late gadolinium enhancement in the 11 HCM patients
LV morphology & risk factors
Asymmetrical septal hypertrophy n (%)
Anterior hypertrophy n (%)
LV outflow tract obstruction at rest n (%)
Number of risk factors for sudden death: median (range)
Presence of late gadolinium enhancement
Extent of late gadolinium enhancement (% of LV mass)
9.1 ± 3.6%
Location of gadolinium enhancement:
> 1 region
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
E2A diastole median (°)
E2A systole median (°)
E2A Mobility (o)
H+LG+ vs H-LG- p-values
H+LG- vs H-LG- p-values
Controls vs H+LG+
Controls vs H+LG-
Controls vs H-LG-
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.
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. email@example.com.
This work was supported by the National Institute of Health Research Cardiovascular Biomedical Research Unit at the Royal Brompton Hospital and Imperial College, London.
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Maron MS: Clinical utility of cardiovascular magnetic resonance in hypertrophic cardiomyopathy. J Cardiovasc Magn Reson. 2012, 14: 13-10.1186/1532-429X-14-13.PubMed CentralView ArticlePubMedGoogle Scholar
- Maron BJ, Maron MS: Hypertrophic cardiomyopathy. Lancet. 2013, 381: 242-55. 10.1016/S0140-6736(12)60397-3.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Hughes SE: The pathology of hypertrophic cardiomyopathy. Histopathology. 2004, 44: 412-27. 10.1111/j.1365-2559.2004.01835.x.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.Google Scholar
- Guizar-Sicairos M, Thurman ST, Fienup JR: Efficient subpixel image registration algorithms. Opt Lett. 2008, 33: 156-8. 10.1364/OL.33.000156.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Stoeck CT, Toussaint N, Boesiger P, Kozerke S: Dual Heart-Phase Cardiac DTI Using Local-Look STEAM. 2012, ISMRM, MilanGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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 pressGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.