Submillimeter diffusion tensor imaging and late gadolinium enhancement cardiovascular magnetic resonance of chronic myocardial infarction
© The Author(s). 2017
Received: 19 August 2016
Accepted: 20 December 2016
Published: 11 January 2017
Knowledge of the three-dimensional (3D) infarct structure and fiber orientation remodeling is essential for complete understanding of infarct pathophysiology and post-infarction electromechanical functioning of the heart. Accurate imaging of infarct microstructure necessitates imaging techniques that produce high image spatial resolution and high signal-to-noise ratio (SNR). The aim of this study is to provide detailed reconstruction of 3D chronic infarcts in order to characterize the infarct microstructural remodeling in porcine and human hearts.
We employed a customized diffusion tensor imaging (DTI) technique in conjunction with late gadolinium enhancement (LGE) cardiovascular magnetic resonance (CMR) on a 3T clinical scanner to image, at submillimeter resolution, myofiber orientation and scar structure in eight chronically infarcted porcine hearts ex vivo. Systematic quantification of local microstructure was performed and the chronic infarct remodeling was characterized at different levels of wall thickness and scar transmurality. Further, a human heart with myocardial infarction was imaged using the same DTI sequence.
The SNR of non-diffusion-weighted images was >100 in the infarcted and control hearts. Mean diffusivity and fractional anisotropy (FA) demonstrated a 43% increase, and a 35% decrease respectively, inside the scar tissue. Despite this, the majority of the scar showed anisotropic structure with FA higher than an isotropic liquid. The analysis revealed that the primary eigenvector orientation at the infarcted wall on average followed the pattern of original fiber orientation (imbrication angle mean: 1.96 ± 11.03° vs. 0.84 ± 1.47°, p = 0.61, and inclination angle range: 111.0 ± 10.7° vs. 112.5 ± 6.8°, p = 0.61, infarcted/control wall), but at a higher transmural gradient of inclination angle that increased with scar transmurality (r = 0.36) and the inverse of wall thickness (r = 0.59). Further, the infarcted wall exhibited a significant increase in both the proportion of left-handed epicardial eigenvectors, and in the angle incoherency. The infarcted human heart demonstrated preservation of primary eigenvector orientation at the thinned region of infarct, consistent with the findings in the porcine hearts.
The application of high-resolution DTI and LGE-CMR revealed the detailed organization of anisotropic infarct structure at a chronic state. This information enhances our understanding of chronic post-infarction remodeling in large animal and human hearts.
KeywordsMyocardial infarction Fiber structure Microstructural remodeling Diffusion tensor imaging Late gadolinium enhancement
Myocardial infarction (MI) is a major cause of death, affecting millions of people worldwide . The occurrence of MI initiates a complex time-dependent and dynamic process of cardiac remodeling that leads to changes in tissue composition, heart geometry, and organ function . While it is clear that MI can lead to heart failure and arrhythmias [3, 4], the exact link between post-MI structural remodeling and the electromechanical functioning of the heart is not completely understood. Accurate knowledge of infarct structure and fiber orientation remodeling in the intact heart is essential for understanding MI pathophysiology. The need for such data is underscored by the fact that infarct structural remodeling is complex and three-dimensional (3D) in nature, which is reflected in the associated changes in cardiac function. However, there is a paucity of data regarding the detailed three-dimensional scar geometry and MI fiber orientation remodeling in intact large animal and human hearts.
Early renditions of myocardial fiber structure were based on sectioning approaches . These methods produce excellent high-resolution data of local tissue structure ; however, combining these measurements together into a registered whole-organ data set, particularly in large animal and human hearts, is extremely difficult. Diffusion Tensor Imaging (DTI) is a non-destructive tool that utilizes the restricted diffusivity of water molecules to assess the tissue microstructure . DTI yields data on the mean diffusivity of water molecules, quantified by Mean Diffusivity (MD), as well as the directional variability of the water diffusion measured by Fractional Anisotropic (FA). Importantly, the principal diffusion eigenvector reflects the mean intravoxel orientation in the tissue. Ventricular fiber maps derived from DTI of formalin-fixed hearts correlate well with histological measurements [8, 9]. DTI has been utilized in infarcted animal and human hearts to characterize structural remodeling using diffusivity measures [10–12], and to assess the remodeling in fiber arrangement at the infarct and in remote locations [13–18]. However, reliable imaging of the fiber orientation at and near the infarct has proven challenging , because the chronically infarcted region is often associated with significant wall thinning, and hence a higher image spatial resolution is required to reliably track the fiber angles across the infarcted wall. In addition, the low diffusion anisotropy in the infarct [12, 14] increases the measurement uncertainties in determining the principal eigenvectors of the diffusion tensors , particularly in a low image signal-to-noise ratio (SNR) regime [19, 20]. Accurate reconstruction of fiber orientation in infarcted hearts thus necessitates DTI sequences that produce high image spatial resolution and high SNR. Our group has recently developed such methodology that has proven successful in imaging myofiber orientation in the thin atrial walls .
The goal of this study was to image the 3D chronic infarct structure at a submillimeter resolution and to provide a systematic analysis of the microstructural remodeling at the infarct. To do so, we employed our previously developed  3D DTI sequence and applied it to image infarcted porcine and human ventricles ex vivo. For the infarcted porcine hearts, this technique was used in conjunction with high-resolution T1-weighted late gadolinium enhancement (LGE) cardiovascular magnetic resonance (CMR), to reconstruct the ventricular fiber organization and scar geometry in the same heart at image resolution and SNR higher than previously achieved . The detailed knowledge of infarct microstructure and fiber orientation obtained in this study is expected to enhance our understanding of post-infarction remodeling that underlies rhythm and pump disorders, thus providing impetus to improvements in targeted therapies.
Specimen acquisition and preparation in porcine hearts
Anteroapical infarction was created by occluding the mid–left anterior descending (LAD) coronary artery in female Yorkshire porcine for 120 min using a balloon angioplasty catheter (n = 8). The hearts were excised at least 3.5 months after the induction of the MI (average MI age: 6.7 ± 2.9 months). As the gold standard for MI imaging, a double dose of Gd-DTPA (Magnevist®) was administered (0.2 mmol/kg) through Intravenous (IV) line 20 min before animal sacrifice. Under anesthesia, the animals were injected with heparin to prevent clot formation in the heart. Further, the heart was arrested using potassium chloride (KCl) to avoid contraction. Excision was performed under 5 min after the animal sacrifice, and the ventricles were filled with rubber (Task5™) to keep the heart in the natural unloaded shape. To avoid specimen dehydration and susceptibility artifacts resulting from the tissue-air interface, the hearts were submerged in perfluorocarbon (Fluorinert-77, 3 M) prior to subsequent imaging. To serve as controls, 4 normal porcine hearts were harvested from animals with no prior MI and prepared in a similar fashion. The weight of the animals at the time of the first procedure (MI induction for the infarcted hearts and the heart harvest for the control hearts) was 50 ± 18 kg, and at the time of harvest for the infarcted hearts was 121 ± 62 kg. Porcine age was approximated as 4 months at the time of the first procedure.
Late gadolinium enhancement CMR
For the infarcted porcine hearts, LGE was performed using a 3D T1-weighted Gradient Echo sequence with radiofrequency (RF) spoiling, approximately 1 h after the excision, to image the whole heart. Imaging was performed with the following parameters: acquired resolution = 0.25 × 0.25 × 0.50 mm3, typical field of view = 110 × 110 × 140 mm3, echo time (TE) = 2.3 ms, repetition time (TR) = 12 ms, flip-angle = 15°, number of averages = 3, scan duration = 1 h.
Diffusion tensor imaging
The specimens were fixed in 10% buffered formaldehyde more than 40 days prior to DTI acquisition. The time from harvest to fixation was about 2–3 h. Our 3D Fast Spin Echo DTI sequence developed previously  was used on a 3T clinical system (Achieva TX, Philips Healthcare, Best, The Netherlands) to image the whole hearts ex vivo (n = 4 control porcine, n = 8 infarcted porcine and one human heart). Imaging parameters were: TE = 63 ms, TR = 504 ms, bandwidth = 290.0 Hz/pixel, number of echoes = 2, diffusion gradients duration = 22.8 ms, time gap between diffusion pulses = 12.6 ms, maximum gradient strength = 60 mT/m, RF coil: Philips 8-channel head coil, number of diffusion encoding directions = 15, maximum b-value = 800 s/mm2, typical field of view = 110 × 115 × 130 mm3, acquired voxel dimension = 0.6 × 0.6 × 1.2 mm3, reconstructed voxel dimension (using zero-padding) = 0.4 mm3, and total scan duration ~ 42 h.
Tensor calculation and tissue segmentation
Fiber angle measurement and visualization
In all the hearts, the primary diffusion eigenvector angles were measured in a local coordinate system tangential to the LV endocardial surface (Fig. 1H). The unit vectors of the three orthogonal axes of this coordinate system, n, t and f, were defined as follows: n: the normal vector to the endocardial surface, t: circumferential vector, such that t = z × n (where z is the longitudinal unit vector directed from apex to base), and f = n × t. The definition of this coordinate system guarantees that t and f are tangential to the heart surface, even in the case that n does not lie within the short-axis plane of the heart due to the heart curvature. As shown in Fig. 1H, the orientation of the primary eigenvector was uniquely determined by two fiber angles measured in this coordinate system. The inclination angle was defined as the angle between the projection of the primary eigenvector onto the tangential plane (defined by t and f) and the local circumferential vector (t). Likewise, the imbrication angle was defined as the angle between the circumferential vector (t) and the projection of the eigenvector onto the plane defined by t and n. In a normal LV wall, the inclination angle changes smoothly from the epicardium with fibers having negative inclination angle (left-handed fibers) to the endocardium with fibers having positive inclination angle (right-handed fibers). The circumferentially running fibers at the midwall have an inclination angle close to zero. The imbrication angle is, on average, close to zero in a normal LV wall. Inclination and imbrication angles were measured at all myocardial voxels throughout each heart. In all the figures in this study, the vector field associated with the primary eigenvector was visualized in TrackVis . This was performed by representing each eigenvector by a small line/tube of fixed dimensions that is located at that voxel and aligns with the primary eigenvector. The fiber orientations were color-coded using the inclination angle (Fig. 1B).
Analysis of regional remodeling: definition of LV segments in porcine hearts
To characterize the 3D regional structural remodeling in MI, LVs of normal and infarcted porcine hearts were partitioned into small transmural segments, as shown in Fig. 1A, E. For each short-axis slice of 1.2 mm thickness (excluding apical slices), a polar coordinate system was defined with an origin at the center of the blood mass. In this coordinate system, each LV was divided circumferentially into 36 segments of 10 angular width per segment. This resulted in about N ~ 1000 segments per heart.
Infarcted segments were defined as the LV segments that contained fibrotic tissue on the LGE images (Fig. 1E, F). These segments were located at the anteroseptal wall, consistent with the LAD infarction protocol. The corresponding segments from the same anatomical locations in the normal hearts were selected as control segments (spanning 180° of the LV anteroseptal wall, as delineated by the green dashed lines in Fig. 1B). Due to the presence of viable tissue surrounding the scar and the complex scar geometry, the transmural infarcted segments could contain both fibrotic and non-fibrotic tissues (Fig. 1F).
Analysis of regional remodeling: quantification of local segment structure in porcine hearts
Subsequent to defining LV transmural segments, voxel data from individual infarcted and control segments were analyzed to determine local fiber angles, wall thickness, and fibrosis extent. The following metrics were defined for each segment: wall thickness, the wall thickness of the segment; scar transmurality, the ratio of number of fibrosis voxels (as obtained from the LGE data) in a segment to the total number of voxels in a segment; slope and intercept, the slope and the intercept of the linear regression model fitted to the inclination angle transmural profile data (angle vs. depth of the wall as measured from the epicardium); r 2 , the coefficient of determination (R-squared) of the regression model of the inclination angle profile; (epi-to-endo) inclination angle range, the range of inclination angle within each segment as calculated from the [1–99]% range of inclination angle distribution; imbrication angle mean, the average of the imbrication angles within a segment; left-handed ratio (LH ratio), the ratio of the number of voxels in which inclination angle < −15° to the total number of voxels within a segment; circumferential ratio, same ratio but calculated for the voxels with −15° < inclination angle < 15°; and right-handed ratio (RH ratio), same ratio but for voxels with inclination angle > 15°. We also quantified local intervoxel incoherency in eigenvector arrangement corresponding to both the inclination and imbrication angles (see the Additional file 1: Methods for the definition).
Statistics were performed on the quantities obtained from the infarcted and control segments throughout the porcine hearts to study the microstructural differences between these regions. Non-parametric Wilcoxon rank-sum test was performed between the averaged segment values from the infarcted (n = 8) and control (n = 4) hearts (per heart statistics). The same statistics was performed on the pooled voxel or segment data (per voxel/segment statistics). All the data is presented as mean ± standard deviation (STD). The statistics presented in the text of the manuscript are from the per heart analysis, unless stated otherwise. Further, Pearson correlations were calculated to study the pair-wise associations of slope, inclination angle range, wall thickness, and scar transmurality in the infarcted segments from porcine hearts.
Imaging and analysis of infarcted human heart specimen
In addition to the porcine hearts, an intact human heart was procured through the National Disease Research Interchange (NDRI, Philadelphia, PA). This heart was from a donor with a history of MI (93 years old female). The human heart specimen had already undergone fixation process at the time of acquisition and hence was not imaged using LGE-CMR (fixation time from harvest: 7 h). It was imaged using the same diffusion tensor imaging sequence at the resolution of 0.5 × 0.5 × 1.0 mm3. To examine infarct remodeling in the human heart, fiber angle measurement and eigenvector visualization were performed in the same way as in the porcine hearts.
LV characteristics of porcine hearts
The eight infarcted porcine hearts demonstrated significant global remodeling with changes in LV mass and LV blood volume (LV mass: 121 ± 47 g in MI vs. 70 ± 9 g in normal hearts, LV blood volume: 81 ± 26 cm3 in MI vs. 41 ± 25 cm3 in normal hearts). The normal hearts had an average wall thickness of 7.2 ± 1.2 mm; by comparison, the average wall thickness in the infarcted porcine hearts was 5.1 ± 1.2 mm for the infarcted regions and 9.5 ± 1.8 mm otherwise. The scar volume comprised 10.1 ± 7.6% of the LV myocardial volume.
Measures of diffusivity and anisotropy in porcine hearts
Diffusion scalar measurements in fibrotic and non-fibrotic tissues in the infarcted hearts, and the normal tissue
Mean ± STD
Fibrotic vs. Non-fibrotic
Fibrotic vs. Normal
Non-fibrotic vs. Normal
e1 (×10−4 mm2/s)
11.11 ± 0.64
8.87 ± 0.93
8.95 ± 0.42
e2 (×10−4 mm2/s)
9.03 ± 0.71
6.06 ± 0.92
5.66 ± 0.37
e3 (×10−4 mm2/s)
7.09 ± 0.88
4.77 ± 0.94
4.39 ± 0.54
MD (×10−4 mm2/s)
9.08 ± 0.72
6.57 ± 0.92
6.33 ± 0.41
0.24 ± 0.04
0.33 ± 0.06
0.37 ± 0.04
e1 (×10−4 mm2/s)
11.25 ± 2.16
9.01 ± 1.52
8.89 ± 1.14
e2 (×10−4 mm2/s)
9.23 ± 2.41
6.25 ± 1.53
5.63 ± 1.22
e3 (×10−4 mm2/s)
7.25 ± 2.47
4.95 ± 1.45
4.35 ± 1.19
MD (×10−4 mm2/s)
9.24 ± 2.26
6.74 ± 1.42
6.29 ± 1.09
0.23 ± 0.09
0.31 ± 0.09
0.37 ± 0.09
Diffusion eigenvector orientation in porcine hearts
Structural metrics in infarcted and control segements
Per heart statistics
Per segment statistics
Wall thickness (mm)
5.5 ± 1.1
7.0 ± 0.9
5.5 ± 2.9
7.0 ± 1.5
0.42 ± 0.12
0.44 ± 0.30
Inclination angle range (°)
111.0 ± 10.7
112.5 ± 6.8
109.9 ± 26.1
112.0 ± 15.4
21.6 ± 7.0
15.7 ± 1.1
21.4 ± 18.7
15.7 ± 3.7
−60.4 ± 9.1
−54.0 ± 5.3
−60.9 ± 28.6
−53.9 ± 13.4
0.69 ± 0.05
0.93 ± 0.03
0.69 ± 0.23
0.93 ± 0.07
Imbrication angle mean (°)
1.96 ± 11.03
0.84 ± 1.47
1.08 ± 16.63
0.86 ± 8.39
Inclination incoherency (°)
10.05 ± 1.42
5.63 ± 0.56
9.61 ± 3.19
5.70 ± 1.31
Imbrication incoherency (°)
8.62 ± 1.77
4.83 ± 0.58
8.14 ± 3.67
4.85 ± 1.85
0.51 ± 0.08
0.39 ± 0.05
0.52 ± 0.22
0.39 ± 0.14
0.26 ± 0.06
0.30 ± 0.04
0.26 ± 0.17
0.30 ± 0.10
0.23 ± 0.06
0.31 ± 0.03
0.21 ± 0.18
0.31 ± 0.11
Diffusion eigenvector orientation in the infarcted human heart
In this work, we studied chronic MI microstructural remodeling in eight porcine hearts and in a human heart non-destructively at submillimeter resolution. We examined the hearts’ fibrous structure using a customized DTI sequence on a clinical scanner (to accommodate the large specimen) that allowed acquisition of data over long scanning times with high image stability . This, in addition to the high SNR and low artifact characteristics of the 3D spin echo sequence used here, resulted in high image quality and high spatial resolution. By combining this technique with high-resolution LGE imaging, we were able to provide reconstructions of both fiber architecture and scar distribution in infarcted hearts with an unprecedented level of detail. The submillimeter voxel size of the data (voxel volume: 0.432 mm3) resulted in an average of 8 voxels across the infarcted wall, allowing for the characterization of the structural remodeling in the zone of infarct, a task that has previously been particularly challenging for DTI due to significant infarct wall thinning and limited image resolution . Importantly, it enabled us to systematically quantify the transmural pattern of diffusion eigenvector orientation in the porcine infarcts and to study the association between the level of structural remodeling and the extent of the infarct. Finally, we applied this technique to a human heart specimen to assess the remodeling at the thinned infarcted wall in the human heart.
The measurement of diffusion scalars in infarcted porcine hearts demonstrated an average of 43% increase in MD and a 35% decrease in FA inside the scar. Similar changes in MD and FA have been reported in DTI studies of myocardial infarcts in various species [11, 12, 14, 25, 26]. The increase in MD is indicative of a less restricted diffusivity and hence, of a larger diffusion volume for water molecules inside the scar. Myocyte death and subsequent collagen deposition following MI could explain this increase in diffusion volume. The same changes lead to alterations in the relative degree of diffusion anisotropy in the tissue, and therefore could be a factor in the reduction of FA. In addition, the dispersion of fiber angles (whether referring to the myofibers or collagen fibers in the infarct) within a voxel, either coherent or incoherent, could alone lead to a decrease in the measured FA due to an averaging effect. While coherent fiber dispersion exists in a normal LV due to the epi-to-endo change in inclination angle, an increase in this dispersion due to an increase in the transmural gradient of fiber angles (such as that shown in Fig. 5A) could, in principle, reduce FA. This would be particularly true for lower image resolution estimates of FA. In addition to this, incoherent dispersion of fibers at the infarct (fiber disarray) has been shown to play significant role in the reduction of FA. Using histological characterization of rat infarcts, Chen et al.  found a good correlation between the amount of fiber disarray and the decrease in FA value. The extent to which factors like fiber dispersion and changes in tissue composition contribute to the measured FA value is unknown and requires further investigation. In our study, despite the lower anisotropy inside the infarct, the measured FA in the majority of infarct regions was greater than that of the isotropic water, which revealed that the chronic scar was mostly comprised of anisotropic structures. This is in agreement with previous histological studies demonstrating that scarred tissue in the porcine heart has a high content and high alignment of collagen . We also found highly aligned collagen bundles from histological imaging of the infarcted region (Additional file 1: Figure S5), likely constituting the main source of anisotropy in the infarct.
The presence of diffusion anisotropy allowed us to employ the diffusion eigenvectors to investigate the anisotropic organization of the collagenous scar. We found that, on average, the orientation of the eigenvectors at the infarcted segments followed the pattern of the original fiber orientation, i.e. left-handed fibers at the epicardium to right-handed fibers at the endocardium. This result is consistent with previous findings in rats  and sheep  but is contrary to findings in another study in rats, which reported severe perturbation of eigenvector orientation in the infarct . In the current study, we provided further evidence for the preservation of primary eigenvector angles and the increase in the slope of transmural angle profile in chronic porcine infarcts and in a human infarct. Despite the preservation of average fiber angles at the infarcted segments and the local collagen fiber alignment in the histology images, we also found that the average intervoxel incoherency of the inclination and imbrication angles in the infarct is higher than the control segments. Monte-Carlo analysis (see Additional file 1: Figure S1B) showed that this increase in the incoherency could not completely stem from the increase in the uncertainty in the estimation of primary eigenvector resulting from reduction in FA in the infarcted region, suggesting the presence of a microstructural basis. Finally, the investigation of structural remodeling at different levels of wall thickness and of the degree of scar transmurality in our study indicates that the inclination angle range is primarily preserved in the infarcted regions and that the increase in the slope is mainly due to wall thinning.
In this study, we found an increase in the proportion of left-handed primary eigenvectors inside the infarcted segments of porcine hearts, which had also been observed in two previous studies of human hearts [13, 18]. We can speculate that this could be due to non-uniform wall thinning across the wall; the sub-endocardial wall and midwall are more likely to be affected by ischemia and undergo more thinning. In support of this, we indeed observed non-uniform distribution of scar across the wall, with fibrosis more concentrated at the sub-endocardial and midwall sub-layers of the infarcted wall (Fig. 7A). The distribution of primary eigenvector angle at the infarcted segments could have also been affected by the patterns of mechanical stress and strain during the process of remodeling. In general, despite slight differences in the shapes of the inclination angle profiles, there was a remarkable similarity between the angle profile trends in control and in infarcted segments. As a proof of concept, preservation of the trend of inclination angle profile across the region of infarct was also observed in the thinned wall of the single human heart in the study.
After MI, the organization of collagen fibers in the scar is influenced by structural and mechanical factors . The original extracellular matrix (ECM) acts as scaffold for the deposition of the new collagen . In addition, the passive stretch of the infarcted tissue by the surrounding myocardium could influence the collagen alignment and hence the anisotropic structure of the scar . The extent to which these factors contribute to scar structure is not known. Our observations in porcine and human ventricles suggest that the existing ECM might play a significant role in the alignment of collagen fibers, as the general transmural patterns of primary diffusion eigenvector angles are preserved inside the scar. This information could be important to therapeutic approaches in the treatment of the infarct, such as tissue engineering and regenerative medicine, as ECM orientation is likely to determine the orientation of regenerated myocytes .
Accurate structural data is essential for the construction and validation of whole-heart computational models [32–35]. These models have enabled advancements in the understanding and treatment of cardiac dysfunction. DTMRI data derived from normal hearts have been utilized by many groups to accurately represent the fiber orientation in cardiac models . However the paucity of data from infarcted hearts has prevented modeling research from fully evaluating the effect of structural remodeling on electrical and mechanical dysfunction in whole-heart cardiac models, particularly in large animal and human hearts . The data presented here provide unprecedented detail about myofiber orientation and scar geometry in infarcted hearts. These data could be employed to construct high-resolution image-based models to investigate the mechanistic links between infarct-related arrhythmias and structural remodeling, including fibrosis distribution and fiber arrangement . In addition, the 3D information regarding the collagen fiber orientation in the intact scar as provided here could improve the accuracy of the modeling approaches aimed at understanding the mechanical role of the passive scar on post-MI ventricular function [38–40]. This could ultimately lead to optimal design of therapies aimed at modifying the electromechanical properties of the infarct [31, 41, 42]. Furthermore, patient-specific models are being constructed from clinical images with clinical applications such as patient arrhythmia risk stratification and optimal treatment planning for rhythm disorders [43, 44]. A complete picture of the high-resolution structural detail in the infarct will assist in making the appropriate assumptions when creating patient-specific models from lower resolution clinical images . The integration of realistic fiber orientation in infarcted hearts, as obtained in this study, will enhance the accuracy of these clinical modeling efforts.
This study has several limitations. First, the LAD reperfusion infarct induction protocol in porcine hearts may render our findings not entirely applicable to all chronic infarctions in the human. Second, although the main findings from the analysis of the porcine hearts regarding the preservation of eigenvector orientation and the increase in the transmural slope were in agreement with the those of the human heart in this study, the single human specimen was not sufficient for a full characterization of infarct remodeling in human hearts and was only a simple existence demonstration. Future studies are needed to investigate this in a larger number of samples. Also since this specimen did not undergo LGE imaging, we were not able to accurately delineate the infarcted area and provide systematic analysis of the remodeling in this heart. Third, despite submillimeter resolution of the imaging with average of ~8 voxels across the wall in the infarct, we excluded the regions with wall thickness less than 2 mm to ensure reliable measurement of transmural angle profile with at least ~4 voxels across the wall; this excluded around 10% of the infarcted wall. Finally, while the fixation process could have affected the baseline values of diffusivity such as FA and MD in our study , it has been shown that it does not change the eigenvector orientation significantly after infarction . So we do not believe it could have affected the main findings of the study.
The application of 3D DTI and LGE-CMR revealed the fiber orientation and scar geometry in infarcted porcine hearts at an unparalleled resolution and SNR. The results demonstrated preservation of eigenvector orientation, with a higher transmural gradient of inclination angle at the thinned wall of infarct. Detailed information of post-infarction remodeling obtained in this study could pave the way for generation of accurate whole-heart models of infarcted hearts to investigate the mechanistic links between the structure and electromechanical function, and thus may lead to improvements in therapies after myocardial infarction.
Diffusion tensor imaging
Diffusion tensor magnetic resonance imaging
Left anterior descending
Late gadolinium enhancement
Magnetic resonance imaging
The authors acknowledge support by the National Institutes of Health (NIH) Pioneer Award (DP1HL123271) to NT, and the NIH intramural research program.
Availability of data and materials
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
All authors contributed to intellectual conception of the study. FP, DH, SM, EM designed the imaging. FP performed the image acquisition and data reconstruction. FP and EM analyzed the data. DH, MZ and HH helped in animal sample collection and LGE imaging. NG and DB performed the human heart specimen collection and contributed to the diffusion image acquisition. FP and DH performed the sample preparation. FP, DH and NT contributed to the histology. FP, NT and EM drafted and revised the manuscript. All authors read and approved the manuscript.
The authors declare that they have no financial or non-financial competing interests.
Consent for publication
Ethics approval and consent to participate
This study was performed with approval from the Johns Hopkins University Institutional Animal Care and Use Committee (IACUC). One intact human heart fixed in formalin was obtained from the National Disease Research Interchange (NDRI, Philadelphia, PA).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Mozaffarian D, Benjamin EJ, Go AS, Arnett DK, Blaha MJ, Cushman M, Das SR, de Ferranti S, Després J-P, Fullerton HJ, Howard VJ, Huffman MD, Isasi CR, Jiménez MC, Judd SE, Kissela BM, Lichtman JH, Lisabeth LD, Liu S, Mackey RH, Magid DJ, McGuire DK, Mohler ER, Moy CS, Muntner P, Mussolino ME, Nasir K, Neumar RW, Nichol G, Palaniappan L, Pandey DK, Reeves MJ, Rodriguez CJ, Rosamond W, Sorlie PD, Stein J, Towfighi A, Turan TN, Virani SS, Woo D, Yeh RW, Turner MB. Heart disease and stroke statistics—2016 update. Circulation. 2016;133:e38–e360.View ArticlePubMedGoogle Scholar
- Ertl G, Frantz S. Healing after myocardial infarction. Cardiovasc res. 2005;66:22–32.View ArticlePubMedGoogle Scholar
- Minicucci MF, Azevedo PS, Polegato BF, Paiva SAR, Zornoff LAM. Heart failure after myocardial infarction: clinical implications and treatment. Clin cardiol. 2011;34:410–4.View ArticlePubMedGoogle Scholar
- de Bakker JM, van Capelle FJ, Janse MJ, Wilde AA, Coronel R, Becker AE, Dingemans KP, van Hemel NM, Hauer RN. Reentry as a cause of ventricular tachycardia in patients with chronic ischemic heart disease: electrophysiologic and anatomic correlation. Circulation. 1988;77:589–606.View ArticlePubMedGoogle Scholar
- Streeter DD, Spotnitz HM, Patel DP, Ross J, Sonnenblick EH. Fiber orientation in the canine left ventricle during diastole and systole. Circ res. 1969;24:339–47.View ArticlePubMedGoogle Scholar
- Rutherford SL, Trew ML, Sands GB, Legrice IJ, Smaill BH. High-resolution 3-dimensional reconstruction of the infarct border zone impact of structural remodeling on electrical activation. Circ res. 2012;111:301–11.View ArticlePubMedGoogle Scholar
- Basser PJ, Mattiello J, Lebihan D. MR diffusion tensor spectroscopy and imaging. Biophys j. 1994;66:259.View ArticlePubMedPubMed CentralGoogle Scholar
- Hsu EW, Muzikant AL, Matulevicius SA, Penland RC, Henriquez CS. Magnetic resonance myocardial fiber-orientation mapping with direct histological correlation. Am j physiol circ physiol. 1998;274:H1627–34. http://ajpheart.physiology.org/content/274/5/H1627.Google Scholar
- Holmes AA, Scollan DF, Winslow RL. Direct histological validation of diffusion tensor MRI in formaldehyde-fixed myocardium. Magn reson med. 2000;44:157–61.View ArticlePubMedGoogle Scholar
- Strijkers GJ, Bouts A, Blankesteijn WM, Peeters THJM, Vilanova A, van Prooijen MC, Sanders HMHF, Heijman E, Nicolay K. Diffusion tensor imaging of left ventricular remodeling in response to myocardial infarction in the mouse. NMR biomed. 2009;22:182–90.View ArticlePubMedGoogle Scholar
- Pop M, Ghugre NR, Ramanan V, Morikawa L, Stanisz G, Dick AJ, Wright G a. Quantification of fibrosis in infarcted swine hearts by ex vivo late gadolinium-enhancement and diffusion-weighted MRI methods. Phys med biol. 2013;58:5009–28.View ArticlePubMedGoogle Scholar
- Kung GL, Ajijola O, Ramirez RJ, Gahm JK, Zhou W, Wisniewski N, Mahajan A, Garfinkel A, Shivkumar K, Ennis D. Microstructural remodeling in the post-infarct porcine heart measured by diffusion tensor MRI and T1-weighted late gadolinium enhancement MRI. J cardiovasc magn reson. 2012;14:2–4.View ArticleGoogle Scholar
- Wu M-T, Tseng W-YI, Su M-YM, Liu C-P, Chiou K-R, Wedeen VJ, Reese TG, Yang C-F. Diffusion tensor magnetic resonance imaging mapping the fiber architecture remodeling in human myocardium after infarction: correlation with viability and wall motion. Circulation. 2006;114:1036–45.View ArticlePubMedGoogle Scholar
- Wu EX, Wu Y, Nicholls JM, Wang J, Liao S, Zhu S, Lau C-P, Tse H-F. MR diffusion tensor imaging study of postinfarct myocardium structural remodeling in a porcine model. Magn reson med. 2007;58:687–95.View ArticlePubMedGoogle Scholar
- Zhang S, Crow JA, Yang X, Chen J, Borazjani A, Mullins KB, Chen W, Cooper RC, McLaughlin RM, Liao J. The correlation of 3D DT-MRI fiber disruption with structural and mechanical degeneration in porcine myocardium. Ann biomed eng. 2010;38:3084–95.View ArticlePubMedGoogle Scholar
- Walker JC, Guccione JM, Jiang Y, Zhang P, Wallace AW, Hsu EW, Ratcliffe MB. Helical myofiber orientation after myocardial infarction and left ventricular surgical restoration in sheep. J thorac cardiovasc surg. 2005;129:382–90.View ArticlePubMedGoogle Scholar
- Mekkaoui C, Huang S, Chen HH, Dai G, Reese TG, Kostis WJ, Thiagalingam A, Maurovich-Horvat P, Ruskin JN, Hoffmann U, Jackowski MP, Sosnovik DE. Fiber architecture in remodeled myocardium revealed with a quantitative diffusion CMR tractography framework and histological validation. J cardiovasc magn reson. 2012;14:70.View ArticlePubMedPubMed CentralGoogle Scholar
- Winklhofer S, Stoeck CT, Berger N, Thali M, Manka R, Kozerke S, Alkadhi H, Stolzmann P. Post-mortem cardiac diffusion tensor imaging: detection of myocardial infarction and remodeling of myofiber architecture. Eur radiol. 2014;24:2810–18.
- Farrell JAD, Landman BA, Jones CK, Smith SA, Prince JL, Van Zijl PCM, Mori S. Effects of signal-to-noise ratio on the accuracy and reproducibility of diffusion tensor imaging – derived fractional anisotropy, mean diffusivity, and principal eigenvector measurements at 1.5 T. J magn reson imaging. 2007;26:756–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Pierpaoli C, Basser PJ. Toward a quantitative assessment of diffusion anisotropy. Magn reson med. 1996;36:893–906.View ArticlePubMedGoogle Scholar
- Pashakhanloo F, Herzka DA, Ashikaga H, Mori S, Gai N, Bluemke DA, Trayanova NA, Mcveigh ER. Myofiber architecture of the human atria as revealed by submillimeter diffusion tensor imaging. Circ arrhythmia electrophysiol. 2016;9:e004133.View ArticleGoogle Scholar
- Jiang H, van Zijl PCM, Kim J, Pearlson GD, Mori S, Van Zijl PCM, Kim J. DtiStudio: resource program for diffusion tensor computation and fiber bundle tracking. Comput methods programs biomed. 2006;81:106–16.View ArticlePubMedGoogle Scholar
- Otsu N. A threshold selection method from gray-level histograms. IEEE trans syst man cybern. 1979;9:62–6.View ArticleGoogle Scholar
- Wang R, Benner T, Sorensen AG, Wedeen VJ. Diffusion toolkit : a software package for diffusion imaging data processing and tractography. Proc int soc magn reson med. 2007;15:3720. http://cds.ismrm.org/ismrm-2007/files/03720.pdf.Google Scholar
- Chen J, Song S-K, Liu W, Mclean M, Allen JS, Tan J, Wickline S a, Yu X. Remodeling of cardiac fiber structure after infarction in rats quantified with diffusion tensor MRI. Am j physiol heart circ physiol. 2003;285:H946–54.View ArticlePubMedGoogle Scholar
- Wu M-T, Su M-YM, Huang Y-L, Chiou K-R, Yang P, Pan H-B, Reese TG, Wedeen VJ, Tseng W-YI. Sequential changes of myocardial microstructure in patients postmyocardial infarction by diffusion-tensor cardiac MR: correlation with left ventricular structure and function. Circ cardiovasc imaging. 2009;2:32–40.View ArticlePubMedGoogle Scholar
- Holmes JW, Covell JW. Collagen fiber orientation in myocardial scar tissue. Cardiovasc pathol. 1996;1:15–22.Google Scholar
- Sosnovik DE, Wang R, Dai G, Wang T, Aikawa E, Novikov M, Rosenzweig A, Gilbert RJ, Wedeen VJ. Diffusion spectrum MRI tractography reveals the presence of a complex network of residual myofibers in infarcted myocardium. Circ cardiovasc imaging. 2009;2:206–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Zimmerman SD, Karlon WJ, Holmes JW, Omens JH, Covell JW. Structural and mechanical factors influencing infarct scar collagen organization. Am j physiol heart circ physiol. 2000;278:H194–200. http://ajpheart.physiology.org/content/278/1/H194.PubMedGoogle Scholar
- Fomovsky GM, Rouillard AD, Holmes JW. Regional mechanics determine collagen fiber structure in healing myocardial infarcts. J mol cell cardiol. 2012;52:1083–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Holmes JW, Laksman Z, Gepstein L. Making better scar: emerging approaches for modifying mechanical and electrical properties following infarction and ablation. Prog biophys mol biol. 2016;120:134–48.View ArticlePubMedGoogle Scholar
- Trayanova NA. Whole-heart modeling: applications to cardiac electrophysiology and electromechanics. Circ res. 2011;108:113–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Deng D, Arevalo H, Pashakhanloo F, Prakosa A, Ashikaga H, Mcveigh E, Halperin H, Trayanova N. Accuracy of prediction of infarct-related arrhythmic circuits from image-based models reconstructed from low and high resolution MRI. Front physiol. 2015;6:1–12.View ArticleGoogle Scholar
- Suinesiaputra A, McCulloch AD, Nash MP, Pontre B, Young AA. Cardiac image modelling: breadth and depth in heart disease. Med image anal. 2016;33:38–43.View ArticlePubMedGoogle Scholar
- Ukwatta E, Arevalo H, Rajchl M, White J, Pashakhanloo F, Prakosa A, Herzka DA, Mcveigh E, Lardo AC, Trayanova NA, Vadakkumpadan F, White J. Image-based reconstruction of three-dimensional myocardial infarct geometry for patient-specific modeling of cardiac electrophysiology. Med phys. 2015;42:4579–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang VY, Nielsen PMF, Nash MP. Image-based predictive modeling of heart mechanics. Annu rev biomed eng. 2015;17:351–83.
- Dillon SM, Allessie MA, Ursell PC, Wit AL, Dillon SM, Allessie MA, Ursell PC, Wit AL. Influences of anisotropic tissue structure on reentrant circuits in the epicardial border zone of subacute canine infarcts SM Dillon, MA allessie, PC Ursell and AL Wit Circulation research 1988, 63:182–206. Circ res. 1988;1988:182–206.View ArticleGoogle Scholar
- Voorhees AP, Han H-C. A model to determine the effect of collagen fiber alignment on heart function post myocardial infarction. Theor biol med model. 2014;11:6.View ArticlePubMedPubMed CentralGoogle Scholar
- Holmes JW, Nuñez J a, Covell JW. Functional implications of myocardial scar structure. Am j physiol. 1997;272:H2123–30. http://ajpheart.physiology.org/content/272/5/H2123.PubMedGoogle Scholar
- Mazhari R, McCulloch AD. Integrative models for understanding the structural basis of regional mechanical dysfunction in ischemic myocardium. Ann biomed eng. 2000;28:979–90.View ArticlePubMedGoogle Scholar
- Fomovsky GM, Macadangdang JR, Ailawadi G, Holmes JW. Model-based design of mechanical therapies for myocardial infarction. J cardiovasc transl res. 2011;4:82–91.View ArticlePubMedGoogle Scholar
- Gourdie RG, Dimmeler S, Kohl P. Novel therapeutic strategies targeting fibroblasts and fibrosis in heart disease. Nat rev drug discov. 2016;15:620–38. doi:10.1038/nrd.2016.89.
- Arevalo HJ, Vadakkumpadan F, Guallar E, Jebb A, Malamas P, Wu KC, Trayanova NA. Arrhythmia risk stratification of patients after myocardial infarction using personalized heart models. Nat commun. 2016;7:11437.View ArticlePubMedPubMed CentralGoogle Scholar
- Krishnamurthy A, Villongco CT, Chuang J, Frank LR, Nigam V, Belezzuoli E, Stark P, Krummen DE, Narayan S, Omens JH, McCulloch AD, Kerckhoffs RCP. Patient-specific models of cardiac biomechanics. J comput phys. 2013;244:4–21.View ArticlePubMedGoogle Scholar
- Mazumder R, Choi S, Clymer BD, White RD, Kolipaka A. Diffusion tensor imaging of healthy and infarcted porcine hearts: study on the impact of formalin fixation. J med imaging radiat sci. 2016;47:74–85.View ArticlePubMedGoogle Scholar