 Research
 Open Access
 Published:
A highresolution cardiovascular magnetic resonance diffusion tensor map from exvivo C57BL/6 murine hearts
Journal of Cardiovascular Magnetic Resonance volume 16, Article number: 77 (2014)
Abstract
Background
The complex cardiac fiber structural organization and spatial arrangement of cardiomyocytes in laminar sheetlets contributes greatly to cardiac functional and contractile ejection patterns. This study presents the first comprehensive, ultrahigh resolution, fully quantitative statistical tensor map of the fixed murine heart at isotropic resolution of 43 μm using diffusion tensor (DT) cardiovascular magnetic resonance (CMR).
Methods
Imaging was completed in approximately 12 hours using a sixdirectional encoding scheme, in five ex vivo healthy C57BL/6 mouse hearts. The tensor map constructed from this data provides an average description of the murine fiber architecture visualized with fiber tractography, and its population variability, using the latest advances in image tensor analysis and statistics.
Results
Results show that nonnormalized cardiac tensor maps are associated with mean fractional anisotropy of 0.25 ± 0.07 and mean diffusivity of 8.9 ± 1.6 × 10^{−4} mm^{2}/s. Moreover, average midventricular helical angle distributions ranged between –41 ± 3° and +52 ± 5° and were highly correlated with transmural depth, in agreement with prior published results in humans and canines. Calculated variabilities of local myocyte orientations were 2.0° and 1.4°. Laminar sheet orientation variability was found to be less stable at 2.6°. Despite such variations, the murine heart seems to be highly structured, particularly when compared to canines and humans.
Conclusions
This tensor map has the potential to yield an accurate mean representation and identification of common or unique features of the cardiac myocyte architecture, to establish a baseline standard reference of DTI indices, and to improve detection of biomarkers, especially in pathological states or posttransgenetic modifications.
Introduction
The unique myocardial architecture, comprised of laminar sheetlets with spatially varying orientation [1][3] and helical fiber tracts, account for the heart’s efficient contractile and torsional mechanical function [4],[5]. 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 [6].
The inherent cardiac tissue structuralfunctional associations also underline the potential significance of myocyte diffusion tensor CMR (DTCMR) tractography [7][11] in remodeling or cellular disarray following both early and late pathological states (e.g. ischemia, cardiomyopathies) [12],[13], pre or postcellular regenerative interventions.
While in vivo cardiac DTCMR tractography is still in the early stages of development for widespread preclinical [14] or clinical use [15][17], ex vivo DTCMR has been established as a useful method to identify local, transmural distributed eigenvector orientation, and to construct threedimensional (3D) myocardial fiber maps [4],[5],[18][21]. 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 [22][24], canines [19],[25][27], goats [28], sheep [21], rabbits [21],[29], rats [12],[13],[30], and mice [4],[21],[31],[32].
Historically, some earlier studies validated the direct correlation of myocyte eigenvector orientations with conventional invasive histology, despite scale differences, in freshly excised [25], perfused [29], and fixed [33] myocardium. More recent highresolution DTCMR [34] and GdDPTA [35] studies validated the localized orthotropic nature of sheep and rat myocardium, respectively, with histology and identified spatial variations in its microlaminar composition and associated branching.
Nevertheless, to date, no previous studies have presented a comprehensive, highresolution, quantitative eigenvector orientation tensor map of the ex vivo murine heart. Additionally, only one prior study [31] 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, spinecho, diffusionweighted CMR [36],[37] to construct a myofiber tensor map of five ex vivo fixed hearts, at an isotropic ultrahigh resolution of 43 μm. Main contributions of this study include: a) construction of a 3D, highresolution, 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 electromechanical function, reflecting an unbiased mean of a test population, serving as a deformable template, and providing tensorbased 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 ultrahigh resolution (exceeding isotropic recordings at 100 μm by [4] in mice and nonisotropic studies at resolutions ranging from 78–2000 μm in other species [12],[19],[23],[24],[26],[33]. 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://lbicy.com/?page_id=80.
Methods
Ethics approval
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 preparationfixation
Male C57BL/6 mice (n = 5, ages 8–12 weeks) were anesthetized via injection of intraperitoneal ketamine/xylazine. Endotracheal intubation was performed with transtracheal illumination using a 22 G angiocatheter 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 precalibrated 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 T_{1} of cardiac tissue, thus allowing for short TR imaging with near complete T_{1} 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 semisolid 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® (perfluoropolyether) (Solvay Solexis, Inc. West Deptford, NJ, USA) and immobilized in a specialized MRcompatible 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 verticalbore Oxford magnet equipped with a gradient system able to achieve peak values of 2000 mT/m. A 3D DWSE 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 b_{0} 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/mm^{2}. Other imaging parameters were: Bandwidth = 62.5 kHz, FOV = 11 × 11 × 11 mm^{3} fullFourier encoding, one average, single echo. All specimens were scanned in a 12 × 25 mm^{2} (diameter × length) solenoid radiofrequency coil.
Image registration using LDDMM
EstimationRegistrations 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 coregistered with the corresponding b_{0} images (Figure 2). Registration was performed using the Advanced Normalization Tools (ANTs) [38]. The geodesic mean, groupwise methodology [5],[18],[23],[24],[39] 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 [5],[18],[39],[23],[24]. 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) [40] with Geodesic Symmetric Normalization (SyN) transformation [38] 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 b_{0} images. The quality of the process was assessed using the union overlap measure (Jaccard coefficient) [41], 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 b_{0} 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 [42].
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 LogEuclidean space of the diffusion tensors (D_{log}), where the individual mouse datasets were divided by the mean value of the tensors’ norm $\sqrt{\mathit{Tr}\left({D}_{\text{log}}^{2}\right)}$, accounting only for voxels within the myocardium [27],[43],[44].
Tensor Averaging – Variability Quantification: The registered/reoriented diffusion tensor fields were then averaged independently at each voxel using LogEuclidean tensor averaging calculus in ANTs [38], resulting in the mean tensor (${\overleftarrow{D}}_{\text{log}}$) (Figure 4) and the corresponding covariance matrix (Σ), in accordance to:
where N is the number of studied hearts, X is the voxel position, vec (ΔD_{i}) is the minimal representation of ΔD_{i}= $\mathit{log}\left({D}_{i}\right){\overleftarrow{D}}_{\mathit{log}}$, and Σ is a representation of the covariance matrix in the logspace [23],[24],[27].
The variability of the mean diffusion tensor was assessed by projecting the covariance matrix on the orthonormal basis W_{ii} and W_{ij} (Equations 3, 4) [27]. Such orthonormal basis allowed the calculation of the eigenvalue (λ_{i}) and eigenvector (v_{i}) variabilities (Equations 5, 6) [27].
The norm $\sqrt{\phantom{\rule{0.12em}{0ex}}\mathit{Tr}\left[\Sigma \left({\rm X}\right)\right]}$ of such covariance matrix subsequently allowed estimation of the relative variability of the diffusion tensor (Figure 5), in accordance to methodology published earlier [27].
Tractography
The calculated mean Diffusion Tensor was then reimported in Diffusion Toolkit where myocyte tracking was performed, employing the secondorder Runge–Kutta method and using an empirically chosen angle threshold of 40° [26]. 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 shortaxis myocardium was segmented by an expert into anterior, lateral, inferior, and septal regions (in accordance to prior AHA guidelines and other publications [21],[45]). Correspondingly, basal (27 slices), middle (58 slices), and apical (57 slices) areas were identified along the longaxis. While segmented shortaxis areas were chosen to have similar areas, the number of basal slices analyzed was not proportional to the number of middlemyocardial 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), epicardialmiddle (EpiMid), middle (Mid), middleendocardial (MidEndo), and endocardial (Endo) regions (Figure 6). Multiple radial profiles were taken spanning sectorial regions for each of the four shortaxis regions, in accordance to methodologies proposed and published previously [21],[30]. 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.
Transmyocyte anisotropy
Differing secondary and tertiary eigenvalues and their associated ratio (λ_{3}/λ_{2}) distributions (Figure 3) are indicative of the existence of possible transmyocyte structural anisotropy, regardless of eigenvalue statistical biases and susceptibility to sorting errors [7]. Cluster sectorial analyses and modalhistogram plots were generated to verify existence of dominant layers in mice (based on computation of intersection angles (φ) [20] 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 zerocrossings or cumulative integration inflexion points.
Statistics
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). Nonparametric Wilcoxon, twotailed ttests 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 oneway ANOVA tests (myocardial location).
For SNR measurements, four circular ROIs (having the same area of 0.23 mm^{2}) 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 mm^{2}). 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.
Results
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 reorientation, averaging and normalization, and quantitative and statistical analyses.
For all mice included in the study, measured dimensions of the heart ranged from 7.47.9 mm (mean ± SD, 7.7 ± 0.2 mm) in long axis (baseapex), 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 [21],[31]. These results also suggest the consistency in tissue fixation, allowing proper quantitative diffusion biomarker estimations and comparisons. Mean SNR values for nonweighted (b0) and fullyweighted 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 [46],[47] on quantitative eigenvalue and eigenvector estimations.
For the tensor map construction, iterative groupwise LDDMM registration was employed (Figures 2 and 3). Three specific diffusion weighted views, indicative of the image quality of the adopted sixencoding directional scheme are shown in Figure 2. The quantitative accuracy of registration was based on estimation of the unionoverlap measure for registeredtarget 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 shortaxis 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 [23],[24], 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 intracardiac and intraspecies 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}) mm^{2}/s, and to global mean FA and MD values of 0.25 ± 0.07 and 8.9 ± 1.6 × 10^{−4} mm^{2}/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 nonnormalized and normalized FA/MD values, as shown on Table 1.
Sectorial statistical quantitative analyses resulted in CV values that ranged between (nonnormalized/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 ttests showed insignificant differences of mean normalized FA and MD for different sectors at the three cardiac levels [BasalMiddleApical]. 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 [23],[24],[27]. 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) [23],[24], 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 counterclockwise epicardial to endocardial sectorial HA variation in all shortaxis planes. The globally estimated HA exhibited a mean transmural variation of 111 ± 23°, in close agreement with prior studies [4] (Table 1).
Significant and extensive evidence from previous work [1][3] 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 [48], sheep [34], and ovine ventricles [49]. Initial evidence for existence of transmyocyte 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 [19].
Discussion
The importance of quantitative DTCMR microscopy atlases and maps to the study of cardiac structure and function has been defined previously [4]. Nevertheless, only two comprehensive cardiac atlas attempts have been documented, on canines [27] and humans [22][24]. 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. [4] and Schmit et al. [31], 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 [4],[31] to include FA, MD, myofiber and laminar sheetlet variability, as well as transmyocyte anisotropy.
For these reasons, this populationbased tensor map defines a crucial juncture in achieving our longterm 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 voxelbased or sectorialbased approach) common local murine myocyte architectural features and populationbased variabilities for quantitative (normalized and nonnormalized) comparisons, in both pathologic and transgenetic states. Additionally, it can serve as a template reference for inter and intraspecies comparisons, addressing questions about structuralfunctional conservation and embryonicadult development.
Ventricular myofiber geometry and orientations were initially studied based on histological reconstructions pioneered by Nielsen et al. [50], with a maximum attained measurement resolution of 500 μm. DTCMR emerged and was established as a noninvasive method to surpass Nielsen’s approach, attaining interspecies, inplane spatial resolutions that ranged from 78–400 μm [4],[12],[19],[33].
While limited success has been evidenced thus far in in vivo animal [14] and human DTCMR studies [16], 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. [51] 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 [4],[12] and correlation of these findings with invasive histology [33].
Basic geometric LV measurements reported in this study also dismiss possible formalinbased fixation effects. They also support consistency in inducing fixation at the enddiastolic conformational state. Additionally, the adopted sixdirectional encoding scheme is in compliance with the realm of a proper acquisition [52][54], yielding comparatively adequately high SNR values [4],[21], dismissing possible biomarker noise bias effects, with anticipated myofiber orientation accuracy of approximately 6°, on the order of Jiang’s study [4]. 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 [42],[55],[56].
Methodologically the adopted data registration schemes, based on LDDMM techniques [19], 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 [27] and human [22],[24] atlases, indicating transmyocyte 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, fieldofviews 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 [27] and humans [23],[24], 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, v_{3} and v_{2}, 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 [23],[24],[27].
Laminar sheetlet orientation variability (assessed through v_{2} and v_{3} orientation about v_{1}) 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 smallsized murine heart having evolved a less variable structure (compared to other species), but also may likely be attributed in part to the highresolution of the acquired data. At an isotropic resolution of 43 μm^{3}, and a typical cardiac myocyte phenotype of columnshaped 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 longaxis 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 nonnormalized) 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 [4],[21],[31],[32], increased sectorial CV was associated with FA. Overall, mean quantitative (nonnormalized) FA results concur with prior studies reporting values ranging from 0.1–0.6 [29] for various species, and exhibiting a mean value of 0.27 ± 0.06 for the murine C57BL/6 myocardium [4].
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. [21], as well as the reported values in other studies [4],[12],[31]. 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. [20] 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. [57], Arts et al. [48], 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 [35],[58],[59].
Conclusions
Quantitative murine myocardial tractography at an isotropic resolution of 43μm of the exvivo 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 midventricular 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 nonnormalized 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 ($\widehat{n}$) 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 gapjunctions.
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).
Authors’ contributions
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 exvivo 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.
Abbreviations
 AHA:

American Heart Association
 ANOVA:

Analysisofvariance
 DTI:

Diffusion tensor imaging
 DWSE:

Diffusion weighted spin echo
 FA:

Fractional anisotropy
 FS:

Finite strain
 HA:

Helical (or Helix) angle
 IVC:

Inferior vena cava
 LDDMM:

Large diffeomorphic metric mapping
 MD:

Mean diffusivity
 MRI:

Magnetic resonance imaging
 PBS:

Phosphate buffered solution
 PPD:

Preservation of principal direction
 SD:

Standard deviation
 SNR:

Signaltonoise ratio
 SyN:

Geodesic symmetric normalization
References
 1.
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): 339347. 10.1161/01.RES.24.3.339.
 2.
LeGrice IJ, Hunter PJ, Smaill BH: Laminar structure of the heart: a mathematical model. Am J Physiol. 1997, 272 (5 Pt 2): H2466H2476.
 3.
Young AA, LeGrice MA, Young IJ, Smaill BH: Extended confocal microscopy of myocardial laminae and collagen network. J Microsc. 1998, 192 (Pt 2): 139150. 10.1046/j.13652818.1998.00414.x.
 4.
Jiang Y, Pandya K, Smithies O, Hsu EW: Threedimensional diffusion tensor microscopy of fixed mouse hearts. Magn Reson Med. 2004, 52 (3): 453460. 10.1002/mrm.20191.
 5.
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, 5176. Chapter 5
 6.
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): 92489253. 10.1073/pnas.1120785109.
 7.
Pierpaoli C, Bassel PJ: Toward a quantitative assessment of diffusion anisotropy. Magn Reson in Med. 1996, 36 (6): 893906. 10.1002/mrm.1910360612. Erratum in: Magn Reson Med 1997, 37(6):972
 8.
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): 777789. 10.1097/RLI.0b013e3181238330.
 9.
Zhukov L, Barr AH: HeartMuscle Fiber Reconstruction from Diffusion Tensor MRI [abstract]. 2003
 10.
Sosnovik DE, Wang R, Dai G, Reese TG, Wedeen VJ: Diffusion MR tractography of the heart. J Cardiovasc Magn Reson. 2009, 11: 4761. 10.1186/1532429X1147.
 11.
Mekkaoui C, NiellesVallespin S, Gatehouse PD, Jackwski MP, Firmin DN, Sosnovik DE: Diffusion MRI tractography of the human heart in vivo at enddiastole and endsystole. J Cardiovasc Magn Reson. 2012, 14 (Supple 1): O4910.1186/1532429X14S1O49.
 12.
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: H946H954.
 13.
Giannakidis A, Rohmer D, Veress AI, Gullberg GT: Diffusion tensor MRIderived 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
 14.
Huang S, Sosnovick DE: Molecular and microstructural imaging of the myocardium. Curr Cardiovasc Imag Rep. 2010, 3 (1): 2633. 10.1007/s124100109007y.
 15.
Tseng WY, Reese TG, Weisskoff RM, Brady TJ, Wedeen VJ: Myocardial fiber shortening in humans: initial results of MR imaging. Radiology. 2000, 216 (1): 128139. 10.1148/radiology.216.1.r00jn39128.
 16.
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: 107113. 10.1002/mrm.10482.
 17.
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: 10361045. 10.1161/CIRCULATIONAHA.105.545863.
 18.
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): 125132. 10.1161/01.RES.0000199396.30688.eb.
 19.
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: 296307. 10.1196/annals.1341.026.
 20.
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: 850859. 10.1002/mrm.20622.
 21.
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: 7481. 10.1186/1532429X1374.
 22.
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 DTMRI. Lecture Notes Comput Sci. 2011, 6666: 171179. 10.1007/9783642210280_22.
 23.
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): 14361447. 10.1109/TMI.2012.2192743.
 24.
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: 207213. 10.1007/9783642283260_21.
 25.
Hsu EW, Muzikant AL, Matulevicius SA, Penland RC, Henriquez CS: Magnetic resonance myocardial fiberorientation mapping with direct histological correlation. Am J Physiol Heart Circ Physiol. 1998, 274: H1627H1634.
 26.
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): 10481057. 10.1016/j.mri.2006.12.008.
 27.
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): 15001514. 10.1109/TMI.2007.907286.
 28.
Geerts L, Bovendeerd P, Nicolay K, Arts T: Characterization of the normal cardiac myofiber field in goat measured with MRdiffusion tensor imaging. Am J Physiol Heart Circ Physiol. 2002, 283: H139H145.
 29.
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: H2308H2318.
 30.
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: H1898H1907. 10.1152/ajpheart.00041.2005.
 31.
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: Threedimensional alignment of the aggregated myocytes in the normal and hypertrophic murine heart. Appl Physiol. 2009, 107 (3): 921927. 10.1152/japplphysiol.00275.2009.
 32.
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): 182190. 10.1002/nbm.1299.
 33.
Holmes AA, Scollan DF, Winslow RL: Direct histological validation of diffusion tensor MRI in formaldehydefixed myocardium. Magn Reson Med. 2000, 44: 157161. 10.1002/15222594(200007)44:1<157::AIDMRM22>3.0.CO;2F.
 34.
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): 10801091. 10.1002/jmri.22725.
 35.
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 highspatial resolution contrastenhanced MRI. Am J Physiol Heart Circ Physiol. 2012, 302: H287H298. 10.1152/ajpheart.00824.2011.
 36.
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
 37.
Jiang Y, Johnson GA: Microscopic diffusion tensor atlas of the mouse brain. Neuroimage. 2011, 56 (3): 12351243. 10.1016/j.neuroimage.2011.03.031.
 38.
Avants BB, Epstein CL, Grossman M, Gee JC: Symmetric diffeomorphic image registration with crosscorrelation: evaluating automated labelling of elderly and neurodegenerative brain. Med Image Anal. 2008, 12 (1): 2641. 10.1016/j.media.2007.06.004.
 39.
Guimond A, Meunier J, Thirion JP: Average brain models: a convergence study. Comput Vis Image Und. 2000, 77 (2): 192210. 10.1006/cviu.1999.0815.
 40.
Beg MF, Miller MI, Trouve A, Younes L: Computing metrics via geodisics on flows of diffeomorphisms. Int J Comput Vis. 2005, 61 (2): 139157. 10.1023/B:VISI.0000043755.93987.aa.
 41.
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): 786802. 10.1016/j.neuroimage.2008.12.037.
 42.
Alexander DC, Pierpaoli C, Basser PJ, Gee JC: Spatial transformation of diffusion tensor magnetic resonance images. IEEE Trans Med Imag. 2001, 20 (11): 11311139. 10.1109/42.963816.
 43.
Arsigny V, Commonwick O, Pennec X, Ayache N: A LogEuclidean framework for statistics on diffeomorphisms. Proc MICCAI. 2006, 9 (Pt 1): 92431.
 44.
Arsigny V, Fillard P, Pennec X, Ayache N: LogEuclidean metrics for fast and simple calculus on diffusion tensors. Magn Reson Med. 2006, 56: 411421. 10.1002/mrm.20965.
 45.
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: 539542. 10.1161/hc0402.102975.
 46.
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: 773785. 10.1016/S0730725X(98)000988.
 47.
Basser PJ, Pajevic S: Statistical artifacts in diffusion tensor MRI (DTMRI) caused by background noise. Magn Reson Med. 2000, 44: 4150. 10.1002/15222594(200007)44:1<41::AIDMRM8>3.0.CO;2O.
 48.
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: H222H2229.
 49.
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): H1324H1330. 10.1152/ajpheart.00813.2004.
 50.
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): H1365H1378.
 51.
Grimm AF, Wohlfart B: Sarcomere lengths at the peak of the lengthtension curve in living and fixed rat papillary muscle. Acta Physiol Scan. 1974, 92 (4): 575577. 10.1111/j.17481716.1974.tb05782.x.
 52.
Hasan KM, Parker DL, Alexander AL: Comparison of gradient encoding schemes for diffusiontensor MRI. Magn Reson Imag. 2001, 13 (5): 769780. 10.1002/jmri.1107.
 53.
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: 17761781.
 54.
Landman BA, Farrell JAD, Jones CK, Smith SA, Prince JL, Mori S: Effects of diffusion weighting schemes on the reproducibility of DTIderived fractional anisotropy, mean diffusivity, and principal eigenvector measurements at 1.5T. Neuroimage. 2007, 36: 11231138. 10.1016/j.neuroimage.2007.02.056.
 55.
Anderson AW: Theoretical analysis of the effects of noise on diffusion tensor imaging. Magn Reson Med. 2001, 46: 11741188. 10.1002/mrm.1315.
 56.
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: 11741188. 10.1002/mrm.1315.
 57.
Dokos S, Smaill BH, Young AA, LeGrice IJ: Shear properties of passive ventricular myocardium. Am J Physiol Circ Physiol. 2002, 283: H2650H2659.
 58.
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: 317331. 10.1016/00222828(74)900741.
 59.
Hales PW, Schneider JE, Burton RA, Wright BJ, Bollensdorff C, Kohl P: Histoanatomical structure of the living isolated rat heart in two contraction states assessed by diffusion tensor MRI. Prog Biophys Mol Biol. 2012, 110 (2–3): 319330. 10.1016/j.pbiomolbio.2012.07.014.
Acknowledgments
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).
Author information
Affiliations
Corresponding author
Additional information
Competing interests
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.
Rights and permissions
About this article
Cite this article
Angeli, S., Befera, N., Peyrat, J. et al. A highresolution cardiovascular magnetic resonance diffusion tensor map from exvivo C57BL/6 murine hearts. J Cardiovasc Magn Reson 16, 77 (2014). https://doi.org/10.1186/s129680140077x
Received:
Accepted:
Published:
Keywords
 Diffusion cardiovascular magnetic resonance
 Cardiac
 Myocyte
 Tractography
 Mouse
 Tensor map, MRI