- Open Access
A high-resolution cardiovascular magnetic resonance diffusion tensor map from ex-vivo C57BL/6 murine hearts
Journal of Cardiovascular Magnetic Resonance volume 16, Article number: 77 (2014)
The complex cardiac fiber structural organization and spatial arrangement of cardiomyocytes in laminar sheetlets contributes greatly to cardiac functional and contractile ejection patterns. This study presents the first comprehensive, ultra-high resolution, fully quantitative statistical tensor map of the fixed murine heart at isotropic resolution of 43 μm using diffusion tensor (DT) cardiovascular magnetic resonance (CMR).
Imaging was completed in approximately 12 hours using a six-directional encoding scheme, in five ex vivo healthy C57BL/6 mouse hearts. The tensor map constructed from this data provides an average description of the murine fiber architecture visualized with fiber tractography, and its population variability, using the latest advances in image tensor analysis and statistics.
Results show that non-normalized cardiac tensor maps are associated with mean fractional anisotropy of 0.25 ± 0.07 and mean diffusivity of 8.9 ± 1.6 × 10−4 mm2/s. Moreover, average mid-ventricular helical angle distributions ranged between –41 ± 3° and +52 ± 5° and were highly correlated with transmural depth, in agreement with prior published results in humans and canines. Calculated variabilities of local myocyte orientations were 2.0° and 1.4°. Laminar sheet orientation variability was found to be less stable at 2.6°. Despite such variations, the murine heart seems to be highly structured, particularly when compared to canines and humans.
This tensor map has the potential to yield an accurate mean representation and identification of common or unique features of the cardiac myocyte architecture, to establish a baseline standard reference of DTI indices, and to improve detection of biomarkers, especially in pathological states or post-transgenetic modifications.
The unique myocardial architecture, comprised of laminar sheetlets with spatially varying orientation - and helical fiber tracts, account for the heart’s efficient contractile and torsional mechanical function ,. In fact, evidence supports the spatial arrangement of cardiac myocytes in minimal surfaces (referred to as generalized helicoids) as an innate compensatory mechanism to optimize organ function and contractile ejection .
The inherent cardiac tissue structural-functional associations also underline the potential significance of myocyte diffusion tensor CMR (DTCMR) tractography - in remodeling or cellular disarray following both early and late pathological states (e.g. ischemia, cardiomyopathies) ,, pre- or post-cellular regenerative interventions.
While in vivo cardiac DTCMR tractography is still in the early stages of development for widespread pre-clinical  or clinical use -, ex vivo DTCMR has been established as a useful method to identify local, transmural distributed eigenvector orientation, and to construct three-dimensional (3D) myocardial fiber maps ,,-. Parallel to visualization of distributed eigenvector orientation, numerous concerted efforts have been made to quantitatively assess transmural helical angle (HA) distributions, fractional anisotropy (FA), and mean diffusivity (MD) in normal and pathological states. Prior ex vivo attempts span studies on humans -, canines ,-, goats , sheep , rabbits ,, rats ,,, and mice ,,,.
Historically, some earlier studies validated the direct correlation of myocyte eigenvector orientations with conventional invasive histology, despite scale differences, in freshly excised , perfused , and fixed  myocardium. More recent high-resolution DTCMR  and Gd-DPTA  studies validated the localized orthotropic nature of sheep and rat myocardium, respectively, with histology and identified spatial variations in its micro-laminar composition and associated branching.
Nevertheless, to date, no previous studies have presented a comprehensive, high-resolution, quantitative eigenvector orientation tensor map of the ex vivo murine heart. Additionally, only one prior study  targeted the extensively studied and stable C57BL/6 mouse strain. That study, however, focused only on quantifying helical angles in normal and hypertrophic hearts.
This study employs 3D, microscopic, spin-echo, diffusion-weighted CMR , to construct a myofiber tensor map of five ex vivo fixed hearts, at an isotropic ultra-high resolution of 43 μm. Main contributions of this study include: a) construction of a 3D, high-resolution, quantitative tensor map of the C57BL/6 murine heart, b) quantitative assessment of the variability of FA, MD, and HA distributions, and c) identification of transmural tissue structure anisotropy and the myocardial laminar sheetlet dominance.
This tensor map will be useful in its potential to a) characterize electro-mechanical function, reflecting an unbiased mean of a test population, serving as a deformable template, and providing tensor-based normalization of study subjects/species; b) yield an accurate mean representation and identify common or unique features of myofiber architecture, c) establish a baseline standard reference of DTI indices, and d) improve the detection of FA, MD, and HA biomarkers, particularly in pathological states.
To the best of our knowledge, this is the first murine distributed eigenvector orientation map constructed at ultra-high resolution (exceeding isotropic recordings at 100 μm by  in mice and non-isotropic studies at resolutions ranging from 78–2000 μm in other species ,,,,,. We hope that this effort will prove useful in mouse phenotyping, targeting regional cardiac function, transgenesis, and molecular imaging. The constructed murine tensor map and DTI data from this work are directly available for public use at http://www.civm.duhs.duke.edu/cardiacMRdiffusion2014 and at http://lbi-cy.com/?page_id=80.
All experimental procedures on mice were designed to minimize or avoid unnecessary pain or discomfort inflicted on the animals. All protocols were approved by the Duke University Institutional Animal Care and Use Committee (IACUC), in accordance with the American Physiological Society’s Guiding Principles for the Care and Use of Animals, European Animal Research directives, and International Guidelines for Animal Research.
Heart excision and sample preparation-fixation
Male C57BL/6 mice (n = 5, ages 8–12 weeks) were anesthetized via injection of intra-peritoneal ketamine/xylazine. Endotracheal intubation was performed with trans-tracheal illumination using a 22 G angio-catheter as endotracheal tube. Mice were then maintained under general anesthesia using 1.5% isoflurane in mixed gas, delivered via a custom rodent ventilator with a tidal volume of 100 μL at a rate of 110 breaths per minute. The right internal jugular vein was surgically isolated and cannulated with a flexible catheter connected to a pre-calibrated infusion pump, loaded with both normal saline and saline mixed in 10% ProHance (Gadoteridol, Bracco Diagnostics Inc., Milan,)/formalin solutions. ProHance is a gadolinium chelate that significantly reduces T1 of cardiac tissue, thus allowing for short TR imaging with near complete T1 recovery between pulses. A 1.3% agarose/PBS solution was prepared and set aside. With the internal jugular venous cannula secured in place, a midline abdominal incision was made to expose the peritoneal cavity. The diaphragm was then carefully incised, dissected away from the cardiac apex, and removed. The descending thoracic aorta and inferior vena cava were exposed and isolated. The anesthetized, ventilated animals were then perfused with approximately 50 mL 0.9% NaCl over a period of 5 minutes. When perfusion began, the inferior vena cava (IVC) and descending thoracic aorta were incised at the level of diaphragmatic hiatus, allowing blood to rapidly clear from the thorax, upper extremities, and head. At the end of 5 minutes, the saline reservoir was exchanged for a 10% ProHance/formalin solution, and the animals were perfused with this solution for an additional 5 minutes, allowing for rapid fixation of the myocardium. Once perfusion fixation was complete, the animals were perfused with a heated 1.3% liquid agarose solution. Heating the solution decreases viscosity such that the solution crosses the pulmonary capillary bed and enters the left side of the heart. The animals were perfused until the gel solution could be seen emanating from the transected descending thoracic aorta. At this point, the infusion pump was switched off, and the aorta and the IVC were clamped to prevent outflow of liquid agarose from the heart. The liquid was then allowed to cool and form a semi-solid gel within the cardiac chambers, preserving anatomic atrial and ventricular shape and preventing chamber collapse. The heart was then carefully excised by transecting the great vessels and pulmonary veins, with attention to maintaining the integrity of the right and left atrial walls. The excised hearts were immersed in Fomblin® (perfluoro-polyether) (Solvay Solexis, Inc. West Deptford, NJ, USA) and immobilized in a specialized MR-compatible specimen tube.
Diffusion tensor CMR
Imaging was performed at the Duke Center of In Vivo Microscopy (Duke University Medical Center, Durham, NC, USA) on a 9.4 T vertical-bore Oxford magnet equipped with a gradient system able to achieve peak values of 2000 mT/m. A 3D DW-SE pulse sequence was used to scan the samples with TE (11.8 ms) and TR (100 ms). The acquisition matrix was 256x256x256, resulting in an isotropic voxel size of 43 μm. Diffusion was encoded using a pair of half sinusoidal gradient lobes of 1600 mT/m in amplitude, 1.3 ms in width, and 6.8 ms in separation. One b0 and six diffusion weighted image stacks were acquired (in a total scan time of 12.7 hours) with an empirically chosen b value of 1852 s/mm2. Other imaging parameters were: Bandwidth = 62.5 kHz, FOV = 11 × 11 × 11 mm3 full-Fourier encoding, one average, single echo. All specimens were scanned in a 12 × 25 mm2 (diameter × length) solenoid radiofrequency coil.
Image registration using LDDMM
Estimation-Registrations of Diffusion Tensors: The diffusion weighted MR images were imported into the Diffusion Toolkit (Massachusetts General Hospital, Boston, MA, USA; http://trackvis.org/dtk/) along with the corresponding gradient table and b–value, and the diffusion tensors for each dataset were calculated (Figure 1). All diffusion datasets were co-registered with the corresponding b0 images (Figure 2). Registration was performed using the Advanced Normalization Tools (ANTs) . The geodesic mean, group-wise methodology ,,,, was adopted. According to this methodology, a source heart image set was chosen (template) from the population of imaged (target) datasets mapped to the template heart ,,,,. This process defined a new template heart, and the methodology was then iteratively repeated until an average geometry was obtained that resided within the center of all target geometries (Figure 2, 3). Specifically, Large Diffeomorphic Metric Mapping (LDDMM)  with Geodesic Symmetric Normalization (SyN) transformation  was applied to register each target dataset to the current template dataset (Figure 2).
Mutual information was defined as the similarity metric of the registered b0 images. The quality of the process was assessed using the union overlap measure (Jaccard coefficient) , by summing the myocardial voxels in the target image that matched myocardial voxels in the template image, and dividing by the total number of myocardial voxels (Figure 2). The resulting transformation matrices and deformation fields obtained from the b0 image registration were subsequently used to register datasets to the reference template of the heart geometry, and to reorient the diffusion tensors using the Preservation of the Principle Direction (PPD) methodology .
Diffusion tensor averaging – construction of mean myocyte orientation tensor map
Diffusion Tensor Normalization: To account for factors responsible for diffusion tensor dispersion not associated with intrinsic variability (for example, temperature), the diffusion tensors were globally normalized. Such normalization was performed in the Log-Euclidean space of the diffusion tensors (Dlog), where the individual mouse datasets were divided by the mean value of the tensors’ norm , accounting only for voxels within the myocardium ,,.
Tensor Averaging – Variability Quantification: The registered/reoriented diffusion tensor fields were then averaged independently at each voxel using Log-Euclidean tensor averaging calculus in ANTs , resulting in the mean tensor () (Figure 4) and the corresponding covariance matrix (Σ), in accordance to:
The variability of the mean diffusion tensor was assessed by projecting the covariance matrix on the orthonormal basis Wii and Wij (Equations 3, 4) . Such orthonormal basis allowed the calculation of the eigenvalue (λi) and eigenvector (vi) variabilities (Equations 5, 6) .
The calculated mean Diffusion Tensor was then re-imported in Diffusion Toolkit where myocyte tracking was performed, employing the second-order Runge–Kutta method and using an empirically chosen angle threshold of 40° . FA, MD, and HA maps were subsequently calculated, constructed, visualized (Figure 6), and quantified. Specifically, a sectorial analysis approach was endorsed (Figure 6), where the short-axis myocardium was segmented by an expert into anterior, lateral, inferior, and septal regions (in accordance to prior AHA guidelines and other publications ,). Correspondingly, basal (27 slices), middle (58 slices), and apical (57 slices) areas were identified along the long-axis. While segmented short-axis areas were chosen to have similar areas, the number of basal slices analyzed was not proportional to the number of middle-myocardial and apical slices (Figure 6). Misregistration errors at the location of large vessels led to exclusion of a number of basal slices from subsequent analyses.
For HA assessment, five transmural myocardial layers were defined (each having approximately 20% of transmural thickness), that included epicardial (Epi), epicardial-middle (Epi-Mid), middle (Mid), middle-endocardial (Mid-Endo), and endocardial (Endo) regions (Figure 6). Multiple radial profiles were taken spanning sectorial regions for each of the four short-axis regions, in accordance to methodologies proposed and published previously ,. Within each wedge, helical angles were transmurally measured along radial samples in 1.5° increments, over 3 consecutive slices chosen in the middle of each respective area.
Differing secondary and tertiary eigenvalues and their associated ratio (λ3/λ2) distributions (Figure 3) are indicative of the existence of possible trans-myocyte structural anisotropy, regardless of eigenvalue statistical biases and susceptibility to sorting errors . Cluster sectorial analyses and modal-histogram plots were generated to verify existence of dominant layers in mice (based on computation of intersection angles (φ)  Figure 7), for basal, middle and apical areas, for all transmural myocardial layers and sectors. Identification of dominant angles was subject to histogram peak detection based on derivative zero-crossings or cumulative integration inflexion points.
All results are reported as mean ± SD. Constancy of FA, MD, and linearity of HA values within short axis areas for basal, middle, and apical myocardium for the constructed tensor map is reported as between (averaged estimates for all the mice studied) coefficient of variabilities (e.g. CV = 100xSD/Mean). Non-parametric Wilcoxon, two-tailed t-tests were also employed (XLSTAT, Paris, France) to statistically compare FA, MD, and HA differences (sectorial and globally averaged values) at the 1% significance level. Assessment for the existence of significant variability for all such recorded parameters was achieved using one-way ANOVA tests (myocardial location).
For SNR measurements, four circular ROIs (having the same area of 0.23 mm2) were placed in the Anterior, Inferior, Lateral and Septal myocardial wall areas in transversal apical/middle/basal slices. The mean signal intensity was measured in these regions (using ImageJ, http://imagej.nih.gov/ij/, Bethesda, MD, USA) along with the standard deviation of the signal in an additional ROI positioned in the background (outside the image object, having an area of 1.84 mm2). The ratio between the mean intensity of the myocardial ROIs and the background SD yielded the SNR for each myocardial region. Given the large estimated SNR values, no bias correction was imposed to such estimates.
The sequential flow and data processing for the methodology adopted for the tensor map construction is summarized in Figure 1. In addition to MR imaging, the remaining steps include diffusion tensor calculation, registration, tensor re-orientation, averaging and normalization, and quantitative and statistical analyses.
For all mice included in the study, measured dimensions of the heart ranged from 7.4-7.9 mm (mean ± SD, 7.7 ± 0.2 mm) in long axis (base-apex), and 4.5–5.2 ± 0.5 mm in short axis (basal diameter). Such results closely match murine cardiac dimensions reported in the in vivo and ex vivo state ,. These results also suggest the consistency in tissue fixation, allowing proper quantitative diffusion biomarker estimations and comparisons. Mean SNR values for non-weighted (b0) and fully-weighted scans (b6) scans were 58.9 ± 8.3 and 33.5 ± 7.5, respectively. Such values are dismissive of the likelihood of presence of significant noise biases , on quantitative eigenvalue and eigenvector estimations.
For the tensor map construction, iterative group-wise LDDMM registration was employed (Figures 2 and 3). Three specific diffusion weighted views, indicative of the image quality of the adopted six-encoding directional scheme are shown in Figure 2. The quantitative accuracy of registration was based on estimation of the union-overlap measure for registered-target images [Jaccard coefficient] that ranged between 88.5–93.6% (mean ± SD; 90.3 ± 2.1%).
The statistical tensor map presented here (visually displayed using tractography in long and short-axis views, as shown in Figure 4) includes reporting on both global and local population myocyte tensor variability of the principle diffusion direction (Figure 5). Disparity in the trace of the tensors is observed in the histograms shown in Figure 3. Therefore, to eliminate possible dependency on diffusion values from external factors (leading to varying scaling of the population tensors), tensor normalization was implemented.
Global diffusion tensor variability is embedded in the covariance matrix (Equation 2); analyses were thus based on the calculation of mean eigenvalues and eigenvectors, but reported as the mode of the standard deviation of the eigenvalue distributions with respect to the mean ,, found to be 1.2% (σλ1), 1.3% (σλ2), and 1.4% (σλ3). Overall, the eigenvalue variability of the constructed tensor map (based on normalized tensors) is within acceptable ranges, indicative of intra-cardiac and intra-species diffusion rate homogeneity. Figure 5 also depicts variability maps for both FA and MD.
Quantitative results (Figure 6) led to mean diffusion tensor eigenvalues of λ1 = 1.12 ± 0.17, λ2 = 0.86 ± 0.17, λ3 = 0.69 ± 0.17 (×10−3) mm2/s, and to global mean FA and MD values of 0.25 ± 0.07 and 8.9 ± 1.6 × 10−4 mm2/s, respectively, for all hearts (before normalization) and to corresponding FA and MD values of 0.021 ± 0.01 and 1.00 ± 0.26 (after tensor normalization). Sectorial myocardial analyses yielded non-normalized and normalized FA/MD values, as shown on Table 1.
Sectorial statistical quantitative analyses resulted in CV values that ranged between (non-normalized/normalized) 3.8–15.8%/2.4–11.1% for MD and 14.7–32.0%/16.6–37.4% for FA. Higher CV values resulted in lower statistical power in biomarker estimation.
Statistical t-tests showed insignificant differences of mean normalized FA and MD for different sectors at the three cardiac levels [Basal-Middle-Apical]. Additionally, ANOVA results indicate no significant variability for FA and MD in all myocardial areas, in justification of homogeneity of distribution of reported values. Such findings demonstrate that diffusion properties are similar among the normal population of studied mice.
Diffusion orientation is associated with differing variability along each of the three principle directions, manifested as a measure of diffusional coherence along the three orthonormal directions ,,. Local myocyte orientation variability in longitudinal and equatorial planes is reported by modes σθ1,3 and σθ1,2; σθ2,3 reports the laminar structure variability. Specifically, the standard deviations exemplifying variability of rotational differences of eigenvector pairs with the third orthonormal (based on the arctan of Equation 6) ,, are found to be 2.0° (σθ1,2), 1.4° (σθ1,3), 2.6° (σθ2,3), respectively. Local myocyte orientation variability has also been traditionally analyzed by the transmural distribution of the HA. These are reported to range between −57° and +60° with typical sectorial transmural trajectories depicted on Figure 6. We found a correlated dependency of HA with transmural depth, with Figure 6 demonstrating a counter-clockwise epicardial to endocardial sectorial HA variation in all short-axis planes. The globally estimated HA exhibited a mean transmural variation of 111 ± 23°, in close agreement with prior studies  (Table 1).
Significant and extensive evidence from previous work - supports the aggregation and organization of cardiac myofibers in sheetlets with spatially varying orientation. Theoretical and histological reconstructions have predicted and verified the existence of two distinct sheetlet populations in canine , sheep , and ovine ventricles . Initial evidence for existence of trans-myocyte anisotropy (Figure 3) was therefore reinforced by the results of Figure 7, showing analyses of intersection angles based on tertiary eigenvector orientation, indicative of a distribution exhibiting sheetlet dominance—in agreement with earlier studies on canines .
The importance of quantitative DTCMR microscopy atlases and maps to the study of cardiac structure and function has been defined previously . Nevertheless, only two comprehensive cardiac atlas attempts have been documented, on canines  and humans -. Thus, this effort introduces, to the best of our knowledge, the first integrated, quantitative cardiac tensor map of the C57BL/6 mouse, one of the most extensively studied and stable murine strains. Compared to previous publications on DTCMR on the same mouse strain by Jiang et al.  and Schmit et al. , this study is unique in that it presents isotropic data at 43-μm resolution, and extends quantitative analyses beyond global eigenvalue estimation and HA variability , to include FA, MD, myofiber and laminar sheetlet variability, as well as trans-myocyte anisotropy.
For these reasons, this population-based tensor map defines a crucial juncture in achieving our long-term goal of establishing a common reference for future phenotypic studies. The tensor map can be used in several ways to help identify (either on a voxel-based or sectorial-based approach) common local murine myocyte architectural features and population-based variabilities for quantitative (normalized and non-normalized) comparisons, in both pathologic and transgenetic states. Additionally, it can serve as a template reference for inter- and intra-species comparisons, addressing questions about structural-functional conservation and embryonic-adult development.
Ventricular myofiber geometry and orientations were initially studied based on histological reconstructions pioneered by Nielsen et al. , with a maximum attained measurement resolution of 500 μm. DTCMR emerged and was established as a non-invasive method to surpass Nielsen’s approach, attaining interspecies, in-plane spatial resolutions that ranged from 78–400 μm ,,,.
While limited success has been evidenced thus far in in vivo animal  and human DTCMR studies , the vast majority of prior attempts focused on ex vivo protocols, mostly conducted under formalin fixation of tissue. For example, using light microscopy, Grimm et al.  quantified average sarcomeric shrinkage of 4.2% in rat papillary muscle, upon formalin fixation. Such results have been consistently supported by recent DTCMR reports indicating minimal formalin fixation effects on primary, secondary and tertiary eigenvalues , and correlation of these findings with invasive histology .
Basic geometric LV measurements reported in this study also dismiss possible formalin-based fixation effects. They also support consistency in inducing fixation at the end-diastolic conformational state. Additionally, the adopted six-directional encoding scheme is in compliance with the realm of a proper acquisition -, yielding comparatively adequately high SNR values ,, dismissing possible biomarker noise bias effects, with anticipated myofiber orientation accuracy of approximately 6°, on the order of Jiang’s study . Possible errors arising due to the adopted encoding scheme are compensated somewhat by the high SNR and the increased spatial resolution of the ex vivo DT MRI ,,.
Methodologically the adopted data registration schemes, based on LDDMM techniques , allow precise voxel and sectorial analyses, preserving the spatial variability of myofiber architecture. Although evaluation of all possible tensor reorientation strategies is beyond the scope of this work, DTI analyses yielded consistent mappings of scalar measures.
Similarly, patterns of eigenvalue and λ3/λ2 distributions highly agree with prior canine  and human , atlases, indicating trans-myocyte anisotropy. However, the reported σλi values (1.2, 1.3, and 1.4%) are by far smaller and less variable than canine (5.35, 6.35, and 8.69%). Of course, care must be exercised in interpreting differences between human, canine, and murine data (especially angular dispersions of tensor eigenvectors), given the disparity of acquired resolutions, field-of-views and the relative cardiac scale sizes of the three species. Regardless, in both human and canine species the pattern seems to be conserved, with the primary eigenvalue exhibiting lesser variability compared to λ2 and λ3. It is highly likely that the primary source of observed increased variability of λ2 and λ3 is attributed to the intrinsically lower signal associated with the respective diffusional processes.
Global tensor variability of 4.98 ± 1.6% compares well with reported values of 10% and 13.2% in canine  and humans ,, respectively, supporting structural stability. Variabilities of myofiber orientation (assessed by angular dispersion patterns of the primary and secondary, and tertiary and primary eigenvector pairs over their orthotropic eigenvector counterparts, v3 and v2, respectively) are reported to be 2.0° and 1.4°, respectively, in comparison to 7.9° /7.7° and 13.0° /11.5° in canine and humans ,,.
Laminar sheetlet orientation variability (assessed through v2 and v3 orientation about v1) seems to be less stable at 2.6°. Such pattern is in agreement with findings in canine and humans, in support of the loosely organized laminar sheetlet structure. Nevertheless, the murine heart seems to be highly structured and stable, exceedingly so when compared to the canine (22.7°) and humans (31.1°). The smaller global tensor variabilities observed in the mouse may be physiologic, with the small-sized murine heart having evolved a less variable structure (compared to other species), but also may likely be attributed in part to the high-resolution of the acquired data. At an isotropic resolution of 43 μm3, and a typical cardiac myocyte phenotype of column-shaped cells of approximately 50–150 μm in length and approximately 20 μm in width, imaging is representative of only very few cells.
The short- and long-axis maps of the norm of the covariance trace of Figure 5 indicate that although certain septal and basal areas exhibit higher spatial variability and heterogeneity across the LV, the myocardium retains a highly organized structure. Similar findings are evident for MD. However, increased spatial dispersion is documented for FA.
In addition to spatial variability, global and sectorial values (normalized and non-normalized) for HA, FA and MD are shown on Table 1. While close agreement of global statistics for all such indices is recorded against prior studies ,,,, increased sectorial CV was associated with FA. Overall, mean quantitative (non-normalized) FA results concur with prior studies reporting values ranging from 0.1–0.6  for various species, and exhibiting a mean value of 0.27 ± 0.06 for the murine C57BL/6 myocardium .
Additionally, the HA distribution of local myocyte orientation has been studied and quantified across myocardial segments. Both the sectorial pattern of transmural helical angles (HA) range and absolute values match those of Healy et al. , as well as the reported values in other studies ,,. The HA spans from −57° – +23° on the epicardium to +15°– +60° on the endocardium, with a noted correlation with transmural depth. Sectorial plots of transmural HA exhibit, however, an increased variability as reflected by the increased SD values (that ranged between various transmural locations ranged from 2–15°). Such variation is a composite reflection of anatomical variability, measurement noise, and registration errors.
Sheetlet dominance is confirmed at mean angles values of 62.7 (±11.8) ° and 99.5 (±11.9) °, in comparison to reports by Helm et al.  reporting angles at 35° and 115° and a mean difference of 78° in canines. Dominant muscle layers were also identified previously by numerous studies, including Dokos et al. , Arts et al. , both reporting angles at 45° and 135°. Such results are consistent with the sliding sheetlet hypothesis, showing that the laminar structure is comprised of two primary orientation populations ,,.
Quantitative murine myocardial tractography at an isotropic resolution of 43-μm of the ex-vivo heart was accomplished in approximately 12 hours of data acquisition, sensitized to 6 diffusional directions. Fiber tractography yielded local myocyte orientation maps associated with mean fractional anisotropy (FA) of 0.25 and mean mid-ventricular HA distributions that ranged between –41° – +52°, in agreement with prior published results in humans and canines. This mouse statistical tensor map is the first of its kind in that it quantitatively estimates both normalized and non-normalized diffusion biomarkers. We believe that the tensor map will be of value in the process of phenotyping pathological and transgenetic cardiac states.
Glossary of terms
CMR: Cardiovascular magnetic resonance
DTI: Diffusion tensor imaging
Fractional anisotropy (FA): Scalar value between zero and one that describes the degree of anisotropy of a diffusion process (i.e. the degree of motional restriction along specific directions). A value of zero means that diffusion is isotropic, i.e. it is unrestricted (or equally restricted) in all directions. A value of one means that diffusion occurs only along one axis and it is fully restricted along all other directions.
Helix (Helical) angle (HA): The angle of obliquity of the principal eigenvector relative to the local short axis plane or relative to local circumferential direction, if necessary, projected radially onto the local wall tangent plane. Positive and negative values can be used for clockwise or anticlockwise rotation, for example as viewed from the lumen outward.
Intersection angle: The angle formed between the local wall plane unit normal () and the projection of the tertiary eigenvector onto the local wall plane (defined by vectors t2 and n).
Laminar structure: An inclusive alternative to `sheetlet and shear layer’ structure.
Mean Diffusivity (MD): The trace of the diffusion tensor; it is a scalar measure of the total diffusion within a voxel.
Myocyte: Basic contractual cell unit of myocardium composed of (among others) myofibers, serially connected sarcomeres via gap-junctions.
Shear layers: Fissures or interstices between adjacent sheetlets. They are thought to allow adjacent sheetlets to change imbrigation angle and shear relative to one another during contractile function.
Sheetlet: Aggregate arrangements of a few cardiac myocytes in planes having microscopic thicknesses and limited spatial extent.
Tractography: Visualization of distributed eigenvector orientations by plotting virtual tracts that, at all points along them, are aligned with local eigenvectors.
Transmural angle: The angle of transmural obliquity, for example of the 2nd eigenvector of diffusion (or of laminar structures visualized directly by microscopy) relative to the local wall tangent plane, measured in the plane perpendicular to the 1st eigenvector direction (or the local myocytes).
SA participated in most of the experimental acquisitions and drafting of the manuscript and carried out the tensor processing and performed computational analyses. NB, EC, and GAJ carried out the ex-vivo fixation preparation, conducted the image acquisitions. JMP, GAJ, and CC participated in the design of the study and execution of the experiments. JMP contributed to the coordination efforts (with CC) for tensorial registration, calculations, and tensor map formation. He was partly involved with tensorial computational processing. CC conceived, designed and coordinated the study and performed the statistical analysis. CC was instrumental to draft the manuscript and coordinate efforts towards its completion and submission. SA, JMP, EC, GAJ, and CC have been involved in the interpretation of data. All authors participated in the design and drafting of the manuscript and its critical revision for important intellectual input. All authors also read and approved the final manuscript. GAJ and CC agree to be accountable for all aspects of the work by ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
American Heart Association
Diffusion tensor imaging
Diffusion weighted spin echo
Helical (or Helix) angle
Inferior vena cava
Large diffeomorphic metric mapping
Magnetic resonance imaging
Phosphate buffered solution
Preservation of principal direction
Geodesic symmetric normalization
Streeter DD, Spotnitz HM, Patel DP, Ross J, Sonneblick EH: Fiber orientation in the canine left ventricle during diastole and systole. Circ Res. 1969, 24 (3): 339-347. 10.1161/01.RES.24.3.339.
LeGrice IJ, Hunter PJ, Smaill BH: Laminar structure of the heart: a mathematical model. Am J Physiol. 1997, 272 (5 Pt 2): H2466-H2476.
Young AA, LeGrice MA, Young IJ, Smaill BH: Extended confocal microscopy of myocardial laminae and collagen network. J Microsc. 1998, 192 (Pt 2): 139-150. 10.1046/j.1365-2818.1998.00414.x.
Jiang Y, Pandya K, Smithies O, Hsu EW: Three-dimensional diffusion tensor microscopy of fixed mouse hearts. Magn Reson Med. 2004, 52 (3): 453-460. 10.1002/mrm.20191.
Helm PA: A Novel Technique for Quantifying Variability of Cardiac Anatomy: Application to the Dyssynchronous Failing Heart. PhD Thesis. 2005, Johns Hopkins University, Biomedical Engineering Department, Johns Hopkins University, Biomedical Engineering Department, 51-76. Chapter 5
Savadjiev P, Strijkers GJ, Bakermans AJ, Pluze E, Zucker SW, Siddiqi K: Heart wall myofibers are arranged in minimal surfaces to optimize organ function. Proc Natl Acad Sci. 2012, 109 (24): 9248-9253. 10.1073/pnas.1120785109.
Pierpaoli C, Bassel PJ: Toward a quantitative assessment of diffusion anisotropy. Magn Reson in Med. 1996, 36 (6): 893-906. 10.1002/mrm.1910360612. Erratum in: Magn Reson Med 1997, 37(6):972
Rohmer D, Sitek A, Gullberg GT: Reconstruction and visualization of fiber and laminar structure in the normal human heart from ex vivo diffusion tensor magnetic resonance imaging (DTMRI) data. Invest Radiol. 2009, 42 (11): 777-789. 10.1097/RLI.0b013e3181238330.
Zhukov L, Barr AH: Heart-Muscle Fiber Reconstruction from Diffusion Tensor MRI [abstract]. 2003
Sosnovik DE, Wang R, Dai G, Reese TG, Wedeen VJ: Diffusion MR tractography of the heart. J Cardiovasc Magn Reson. 2009, 11: 47-61. 10.1186/1532-429X-11-47.
Mekkaoui C, Nielles-Vallespin S, Gatehouse PD, Jackwski MP, Firmin DN, Sosnovik DE: Diffusion MRI tractography of the human heart in vivo at end-diastole and end-systole. J Cardiovasc Magn Reson. 2012, 14 (Supple 1): O49-10.1186/1532-429X-14-S1-O49.
Chen J, Song SL, Liu W, McLean M, Allen SJ, Tan J, Wickline SA, Yu X: Remodeling of cardiac fiber structure after infarction in rats quantified with diffusion tensor MRI. Am J Phys Heart Circ Physiol. 2003, 285: H946-H954.
Giannakidis A, Rohmer D, Veress AI, Gullberg GT: Diffusion tensor MRI-derived myocardial fiber disarray in hypertensive left ventricular hypertrophy: Visualization, quantification and the effect on mechanical function. 2012, Blackwell Publishing Limited, Chichester, West Sussex, United Kingdom
Huang S, Sosnovick DE: Molecular and microstructural imaging of the myocardium. Curr Cardiovasc Imag Rep. 2010, 3 (1): 26-33. 10.1007/s12410-010-9007-y.
Tseng WY, Reese TG, Weisskoff RM, Brady TJ, Wedeen VJ: Myocardial fiber shortening in humans: initial results of MR imaging. Radiology. 2000, 216 (1): 128-139. 10.1148/radiology.216.1.r00jn39128.
Dou J, Tseng WY, Reese TG, Weisskoff RM, Wedeen VJ: Combined diffusion and strain MRI reveals structure and function of human myocardial laminar sheets in vivo. Magn Reson Med. 2003, 50: 107-113. 10.1002/mrm.10482.
Wu MT, Tseng WY, Su MY, Liu CP, Chiou KR, Wedeen VJ, Reese TG, Yang CF: Diffusion tensor magnetic resonance imaging mapping the fiber architecture remodeling in human myocardium after infraction: correlation with viability and wall motion. Circulation. 2006, 114: 1036-1045. 10.1161/CIRCULATIONAHA.105.545863.
Helm PA, Younes L, Beg MF, Ennis DB, Leclercq C, Faris OP, McVeigh E, Kass D, Miller MI, Winslow RL: Evidence of structural remodeling in the dyssynchronous failing heart. Circ Res. 2005, 98 (1): 125-132. 10.1161/01.RES.0000199396.30688.eb.
Helm P, Beg MF, Miller MI, Winslow RL: Measuring and mapping cardiac fiber and laminar architecture using diffusion tensor MRI imaging. Ann NY Acad Sci. 2005, 1047: 296-307. 10.1196/annals.1341.026.
Helm PA, Tseng HJ, Younes L, McVeigh ER, Winslow RL: Ex vivo 3D diffusion tensor imaging and quantification of cardiac laminar structure. Magn Reson Med. 2005, 54: 850-859. 10.1002/mrm.20622.
Healy LJ, Jiang Y, Hsu EW: Quantitative comparison of myocardial fiber structure between mice, rabbit, and sheep using diffusion tensor cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2011, 13: 74-81. 10.1186/1532-429X-13-74.
Lombaert H, Peyrat JM, Croisille P, Rapacchi S, Fanton L, Clarysse P, Delingette H, Ayache N: Statistical analysis of the human cardiac fiber architecture from DT-MRI. Lecture Notes Comput Sci. 2011, 6666: 171-179. 10.1007/978-3-642-21028-0_22.
Lombaert H, Peyrat JM, Croisille P, Rapacchi S, Fanton L, Cheriet F, Clarysse P, Magnin I, Delingette H, Ayache N: Human atlas of the cardiac fiber architecture: study on a healthy population. IEEE Trans Med Imag. 2012, 31 (7): 1436-1447. 10.1109/TMI.2012.2192743.
Lombaert H, Peyrat JM, Fanton L, Cheriet F, Delinget H, Ayache N, Clarysse P, Magnin I, Croisille P: Statistical atlas of human fibers: comparison with abnormal heart. Imag Modell Challenges Lecture Notes Comput Sci. 2012, 7085: 207-213. 10.1007/978-3-642-28326-0_21.
Hsu EW, Muzikant AL, Matulevicius SA, Penland RC, Henriquez CS: Magnetic resonance myocardial fiber-orientation mapping with direct histological correlation. Am J Physiol Heart Circ Physiol. 1998, 274: H1627-H1634.
Wu EX, Wu Y, Tang H, Wang J, Yang J, Ng MC, Yang ES, Chan CW, Zhu S, Lau CP, Tse HD: Study of myocardial fiber pathway using magnetic resonance diffusion tensor imaging. Magn Reson Imaging. 2007, 25 (7): 1048-1057. 10.1016/j.mri.2006.12.008.
Peyrat JM, Sermesant M, Pennec X, Delingette H, Xu C, McVeigh ER, Ayanche N: A computational framework for the statistical analysis of cardiac diffusion tensors: application to a small database of canine hearts. IEEE Trans Med Imag. 2007, 26 (11): 1500-1514. 10.1109/TMI.2007.907286.
Geerts L, Bovendeerd P, Nicolay K, Arts T: Characterization of the normal cardiac myofiber field in goat measured with MR-diffusion tensor imaging. Am J Physiol Heart Circ Physiol. 2002, 283: H139-H145.
Scollan D, Holmes A, Winslow R, Forder JR: Histological validation of myocardial microstructure obtained from diffusion tensor magnetic resonance imaging. Am J Physiol Heart Circ Physiol. 1998, 275: H2308-H2318.
Chen J, Liu W, Zhang H, Lacy L, Yang X, Song SK, Wickline SA, 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-H1907. 10.1152/ajpheart.00041.2005.
Schmitt B, Fedarava K, Falkenberg J, Rothaus K, Bodhey NK, Reischauer C, Kozerke S, Schnackenburg B, Westermann D, Lunkenheimer PP, Andrerson RH, Berger F, Kuehne TJ: Three-dimensional alignment of the aggregated myocytes in the normal and hypertrophic murine heart. Appl Physiol. 2009, 107 (3): 921-927. 10.1152/japplphysiol.00275.2009.
Strijkers GJ, Bouts A, Blankesteijn WM, Peeters TH, Vilanova A, van Prooijen MC, Sanders HM, Heijman E, Nicolay K: Diffusion tensor imaging of left ventricular remodelling in response to myocardial infarction in the mouse. NMR in Biomed. 2009, 22 (2): 182-190. 10.1002/nbm.1299.
Holmes AA, Scollan DF, Winslow RL: Direct histological validation of diffusion tensor MRI in formaldehyde-fixed myocardium. Magn Reson Med. 2000, 44: 157-161. 10.1002/1522-2594(200007)44:1<157::AID-MRM22>3.0.CO;2-F.
Kung G, Nguyen TC, Itoh A, Skare S, Ingels NB, Miller DC, Ennis DB: The presence of two myocardial sheet populations confirmed by diffusion tensor MRI and histological validation. J Magn Reson Imaging. 2011, 34 (5): 1080-1091. 10.1002/jmri.22725.
Gilbert SH, Benoist D, Benson AP, White AP, White E, Tanner SF, Holden AV, Dobrynski 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-H298. 10.1152/ajpheart.00824.2011.
Angeli S, Befera N, Cofer G, Johnson GA, Constantinides C: Construction of a Fiber Atlas of the Murine Heart [abstract]. Proceedings of the International Society of Magnetic Resonance in Medicine. 2013, 1007-
Jiang Y, Johnson GA: Microscopic diffusion tensor atlas of the mouse brain. Neuroimage. 2011, 56 (3): 1235-1243. 10.1016/j.neuroimage.2011.03.031.
Avants BB, Epstein CL, Grossman M, Gee JC: Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labelling of elderly and neurodegenerative brain. Med Image Anal. 2008, 12 (1): 26-41. 10.1016/j.media.2007.06.004.
Guimond A, Meunier J, Thirion JP: Average brain models: a convergence study. Comput Vis Image Und. 2000, 77 (2): 192-210. 10.1006/cviu.1999.0815.
Beg MF, Miller MI, Trouve A, Younes L: Computing metrics via geodisics on flows of diffeomorphisms. Int J Comput Vis. 2005, 61 (2): 139-157. 10.1023/B:VISI.0000043755.93987.aa.
Klein A, Andersson J, Ardekani BA, Ashburner J, Avants B, Chiang MC, Christensen GE, Collins DL, Gee J, Hellier P, Song JH, Jenkinson M, Lepage C, ruechert D, Thompson P, Vercauteren T, Woods RP, Mann JJ, Persey RV: Evaluation of 14 nonlinear deformation algorithms applied to human brain MRI registration. Neuroimage. 2009, 46 (3): 786-802. 10.1016/j.neuroimage.2008.12.037.
Alexander DC, Pierpaoli C, Basser PJ, Gee JC: Spatial transformation of diffusion tensor magnetic resonance images. IEEE Trans Med Imag. 2001, 20 (11): 1131-1139. 10.1109/42.963816.
Arsigny V, Commonwick O, Pennec X, Ayache N: A Log-Euclidean framework for statistics on diffeomorphisms. Proc MICCAI. 2006, 9 (Pt 1): 924-31.
Arsigny V, Fillard P, Pennec X, Ayache N: Log-Euclidean metrics for fast and simple calculus on diffusion tensors. Magn Reson Med. 2006, 56: 411-421. 10.1002/mrm.20965.
Cerqueira MD, Weissman NJ, Dilisizian V, Jacobs AK, Kaul S, Laskey WK, Pennell DJ, Rumberger JA, Ryan T, Verani MS: Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart: a statement for healthcare professionals from the cardiac imaging committee of the council on clinical cardiology of the American Heart Association. Circulation. 2002, 105: 539-542. 10.1161/hc0402.102975.
Bastin ME, Armitage PA, Marchall I: A theoretical study of the effect of experimental noise on the measurement of anisotropy in diffusion imaging. Magn Reson Imag. 1998, 16: 773-785. 10.1016/S0730-725X(98)00098-8.
Basser PJ, Pajevic S: Statistical artifacts in diffusion tensor MRI (DT-MRI) caused by background noise. Magn Reson Med. 2000, 44: 41-50. 10.1002/1522-2594(200007)44:1<41::AID-MRM8>3.0.CO;2-O.
Arts T, Costa KD, Covell JW, McCulloch AD: Relating myocardial architecture to shear strain and muscle fiber orientation. Am J Physiol Heart Circ Physiol. 2001, 280: H222-H2229.
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 (3): H1324-H1330. 10.1152/ajpheart.00813.2004.
Nielsen PM, LeGrice IJ, Smaill BH, Hunter PJ: Mathematical models of geometry and fibrous structure of the heart. Am J Physiol. 1991, 260 (4 Pt 2): H1365-H1378.
Grimm AF, Wohlfart B: Sarcomere lengths at the peak of the length-tension curve in living and fixed rat papillary muscle. Acta Physiol Scan. 1974, 92 (4): 575-577. 10.1111/j.1748-1716.1974.tb05782.x.
Hasan KM, Parker DL, Alexander AL: Comparison of gradient encoding schemes for diffusion-tensor MRI. Magn Reson Imag. 2001, 13 (5): 769-780. 10.1002/jmri.1107.
Ni H, Kavcic V, Zhu T, Ekholm S, Zhong J: Effects of number of diffusion of gradient directions on derived diffusion tensor imaging indices in human brain. AJNR Am J Neuroradiol. 2006, 27: 1776-1781.
Landman BA, Farrell JAD, Jones CK, Smith SA, Prince JL, Mori S: Effects of diffusion weighting schemes on the reproducibility of DTI-derived fractional anisotropy, mean diffusivity, and principal eigenvector measurements at 1.5T. Neuroimage. 2007, 36: 1123-1138. 10.1016/j.neuroimage.2007.02.056.
Anderson AW: Theoretical analysis of the effects of noise on diffusion tensor imaging. Magn Reson Med. 2001, 46: 1174-1188. 10.1002/mrm.1315.
Alexander AL, Hasan KM, Tsuruda M, Parker DL: Analysis of partial volume effects in diffusion tensor imaging indices in human brain. Magn Reson Med. 2001, 46: 1174-1188. 10.1002/mrm.1315.
Dokos S, Smaill BH, Young AA, LeGrice IJ: Shear properties of passive ventricular myocardium. Am J Physiol Circ Physiol. 2002, 283: H2650-H2659.
Spotnitz HM, Soitnitz 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-331. 10.1016/0022-2828(74)90074-1.
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 (2–3): 319-330. 10.1016/j.pbiomolbio.2012.07.014.
Support was received (PI: C. Constantinides) from grant IPE/TEXNOLOGIA/MHXAN/0609 (BE)/05 from the Research Promotion Foundation and personal funds (C. Constantinides). All imaging was performed at the Duke Center for In Vivo Microscopy, an NIH/NIBIB National Biomedical Technology Resource Center (P41 EB015897).
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Angeli, S., Befera, N., Peyrat, JM. et al. A high-resolution cardiovascular magnetic resonance diffusion tensor map from ex-vivo C57BL/6 murine hearts. J Cardiovasc Magn Reson 16, 77 (2014). https://doi.org/10.1186/s12968-014-0077-x
- Diffusion cardiovascular magnetic resonance
- Tensor map, MRI