Skip to main content

Novel insight into the detailed myocardial motion and deformation of the rodent heart using high-resolution phase contrast cardiovascular magnetic resonance

Abstract

Background

Phase contrast velocimetry cardiovascular magnetic resonance (PC-CMR) is a powerful and versatile tool allowing assessment of in vivo motion of the myocardium. However, PC-CMR is sensitive to motion related artifacts causing errors that are geometrically systematic, rendering regional analysis of myocardial function challenging. The objective of this study was to establish an optimized PC-CMR method able to provide novel insight in the complex regional motion and strain of the rodent myocardium, and provide a proof-of-concept in normal and diseased rat hearts with higher temporal and spatial resolution than previously reported.

Methods

A PC-CMR protocol optimized for assessing the motion and deformation of the myocardium in rats with high spatiotemporal resolution was established, and ten animals with different degree of cardiac dysfunction underwent examination and served as proof-of-concept. Global and regional myocardial velocities and circumferential strain were calculated, and the results were compared to five control animals. Furthermore, the global strain measurements were validated against speckle-tracking echocardiography, and inter- and intrastudy variability of the protocol were evaluated.

Results

The presented method allows assessment of regional myocardial function in rats with high level of detail; temporal resolution was 3.2 ms, and analysis was done using 32 circumferential segments. In the dysfunctional hearts, global and regional function were distinctly altered, including reduced global peak values, increased regional heterogeneity and increased index of dyssynchrony. Strain derived from the PC-CMR data was in excellent agreement with echocardiography (r = 0.95, p < 0.001; limits-of-agreement −0.02 ± 3.92%strain), and intra- and interstudy variability were low for both velocity and strain (limits-of-agreement, radial motion: 0.01 ± 0.32 cm/s and −0.06 ± 0.75 cm/s; circumferential strain: -0.16 ± 0.89%strain and −0.71 ± 1.67%strain, for intra- and interstudy, respectively).

Conclusion

We demonstrate, for the first time, that PC-CMR enables high-resolution evaluation of in vivo circumferential strain in addition to myocardial motion of the rat heart. In combination with the superior geometric robustness of CMR, this ultimately provides a tool for longitudinal studies of regional function in rodents with high level of detail.

Background

The intricate motion and deformation of the heart can be assessed in vivo with varying degree of detail using several different techniques, including sonomicrometry, echocardiography employing Tissue Doppler Imaging or speckle-tracking strain analysis, as well as various cardiovascular magnetic resonance (CMR) techniques. Compared to other modalities, CMR offers measurements of true 3D function with practically no limitations in visualization geometry, and thus provides complete freedom in choosing regions for examination. Different techniques for the evaluation of myocardial function are available, including myocardial tagging [1], strain-encoded CMR (SENC) [2], displacement-encoded imaging with stimulated echoes (DENSE) [3] and phase contrast imaging (PC-CMR) [4]. The latter two offer pixel-wise measurement of displacement and velocity, respectively, allowing for high-resolution evaluation of tissue function. However, PC-CMR is the only CMR technique that been shown to allow assessment of velocity with high spatial and temporal resolution [5] and, subsequently, displacement [6], strain rate [7] and strain [8] concurrently, throughout the entire cardiac cycle. While PC-CMR velocimetry is emerging as a powerful and versatile tool for assessment of tissue motion both in humans and in rodents [5, 9], it has only been reported so far for the assessment of average myocardial velocities within a slice, or in a few (<=8) circumferential segments. Also, it might be challenging to achieve the optimal temporal resolution to capture the fine details of cardiac motion in small animals. Thus, to accurately investigate the regional function e.g. in hearts with infarctions with various sizes, and in order to derive parameters such as subtle dyssynchrony, transmural functional gradients or longitudinal spread of dysfunction, data with higher spatial and temporal resolution along with appropriate post-processing procedures are essential.

In PC-CMR, the displacement of spins between the centers of the bipolar encoding gradient lobes is encoded into the phase of the MR signal, producing datasets with near-instantaneous velocities with temporal resolution equal to the TR [10, 11]. Although PC-CMR is prone to errors from several sources, including concomitant gradient [12] and eddy-current induced [13] artifacts, several approaches has been proposed to minimize these errors [12, 14, 15]. However, PC-CMR encoded acquisitions are intrinsically non-motion compensated and thus particularly sensitive to flow- and motion related artifacts, such as ghosting due to beat-to-beat variation in blood flow [16, 17]. In Cartesian imaging, motion-related artifacts manifest in the phase-encoding direction, and may (in a short-axis view) affect the measurements non-uniformly over the circumference of the myocardium. This could lead to systematic errors when regional myocardial function is to be assessed, reducing the effective level of detail available for analysis. Black-blood contrast is therefore essential in PC-CMR of the myocardium, reducing this effect [5]. In addition, the impact of directionally dependent artifacts in studies employing signal averaging can be reduced by an in-plane rotation of the field-of-view (FOV) between the acquisitions, referred to as rotating FOV (Figure 1).

Figure 1
figure 1

Geometrically systematic artifacts. Even after black-blood preparation, some flow artifacts remain. Although not obvious in the original end-systolic (ES) magnitude images (a,b), the difference between the magnitude images of end-diastole (ED) and end-systole (c,d) reveal geometrically systematic artifacts. A rotation of the field-of-view alters the characteristics of these artifacts. Both acquisitions exhibit unique artifacts in the phase-encoding direction (solid ellipses), not present in the other (dashed ellipse). In the average image (e), these artifacts are reduced. For illustrative purposes, a median filter has been applied for noise reduction in this figure. In c-e the grayscale has been exaggerated for clarity.

In this study, we present a PC-CMR approach for assessing left ventricular (LV) myocardial motion in rats, employing a rotating FOV along with optimized acquisition parameters and post-processing protocols. We aim to demonstrate the feasibility of PC-CMR for accurately describing both myocardial motion and deformation in rodents by deriving global parameters such as peak velocities and maximum circumferential strain, as well as describing the spatial variability in essential parameters with higher resolution than have been previously reported. By applying this method on rats with myocardial infarction, we were able to describe distinct alterations in regional myocardial function, compared to sham-operated controls. These findings were validated against speckle-tracking echocardiography, confirming our protocol with PC-CMR as an accurate tool for investigating regional myocardial function in rats.

Methods

Experimental animals

Myocardial infarction was induced in male Wistar rats (10–12 weeks old, ~300 g) as described previously [18], where the left coronary artery was occluded by a silk suture. Six weeks later, cardiac imaging was performed. Inclusion criterion was visible myocardial infarction on echocardiography, and we selected both rats with small and large infarctions (infarct size range 21%-46%). Sham-operated rats went through the same procedure, except no coronary artery ligation was performed. Animals were 16–18 weeks old at the time of cardiac imaging.

The animal weights during examination were in the range 400–450 g. For evaluation of the PC-CMR method, the animals (total N = 15) were divided into two groups, one cohort to validate the method relative to echocardiography (post-MI (N = 6), sham (N = 3)), and one cohort to perform inter- and intrastudy variability analysis (post-MI (N = 4), sham (N = 2)). All animals were cared for according to the Norwegian Animal Welfare Act. The use of animals was approved by the Norwegian Animal Research Authority (ID 3284), and conformed to the Guide for the Care and Use of Laboratory Animals published by the US National Institutes of Health and the European Convention for the Protection of Vertebrate Animals used for Experimental and Other Scientific Purposes (ETS no. 123).

Anesthesia was induced in a chamber with a mixture of O2 and 4.5% isoflurane, and maintained during experiments by administration of a mixture of O2 and 1.5% isoflurane in freely breathing animals. During the CMR experiments, body temperature was maintained using heated air; and ECG, respiration rate and animal temperature were constantly monitored. Respiration was registered by an air cushion. The heart rate was kept as constant as practically possible during experiments by minor adjustments of the level of anesthesia.

Echocardiography

Echocardiography examinations were performed on a Vevo 2100 (Visual Sonics Inc., Ontario, Canada) scanner with a 24 MHz transducer approximately one day prior to CMR scan. For the animals included in the validation study, global circumferential strain was calculated off-line by 2D speckle-tracking [19] in mid-ventricular short-axis slices. The temporal resolution of the echocardiography data was 152 frames per second.

CMR hardware and acquisition

CMR experiments were performed on a 9.4 T/210 mm/ASR horizontal bore magnet (Agilent Technologies, Inc., USA) with a high-performance actively shielded gradient coil (inner diameter 120 mm, rise time 180 μs, max strength 600 mT/m). A quadrature volume transmit coil (inner diameter 72 mm) was used in combination with a four-channel surface receive coil array dedicated to rat heart imaging. Bipolar motion encoding gradients were incorporated into an RF-spoiled gradient echo cine sequence [9]. The acquisition employed a nine-point balanced scheme [15] encoding motion in three orthogonal directions, and was non-interleaved, that is, each encoding step was recorded separately with full temporal resolution [20]. Black-blood contrast was achieved by placing saturation slices above and below the imaging slice, applied at the end of the cine train [5]. To allow for some decay of eddy currents in the system, a short delay (τ = 250 μs) was introduced between the second encoding gradient lobe and the readout gradient.

In each animal, a mid-ventricular LV short-axis slice was planned as described by Schneider et al. [21]. Acquisition was prospectively triggered by ECG R-peak and gated for respiration by pausing acquisition during respiratory motion. 70-80 time frames were recorded covering >130% of the r-r-interval, with a temporal resolution of 3.2 ms. Overshooting the r-r-interval permitted complete coverage of the diastole in combination with black-blood saturation; and provided time for decay of gradient-induced disturbances in the ECG signal.

To reduce the impact of directional-dependent artifacts, each slice was acquired twice (corresponding to 2x signal averaging), where the second acquisition was rotated (in-plane) at least 30 degrees with respect to the first.

Key imaging parameters were as follows: TE/TR = 2.2/3.2 ms; FOV = 50x50 mm; matrix 128x128, slice thickness 1.5 mm, flip angle 7°, receiver bandwidth = 156.25 kHz; venc = 13.9 cm/s. Acquisition time for a complete slice was 10-15 minutes, depending on heart rate.

CMR post-processing

The phase contrast data was extracted from the multi-receiver array coil data as described by Bernstein et al. [22], including complex-conjugated multiplication of the encoded scans with the reference scans individually for each coil element. Spatially-specific ECC was performed as previously described [13], excluding areas in the FOV subject to fold-over artifacts from analysis [15].

The data sets were subsequently semi-automatically segmented and analyzed using a purpose-written Matlab software (The MathWorks, Natick, USA). The only human inputs required during post-processing were 1) tracking of subendo- and subepicardial border at key time frames (segmentation of intermediate time frames were automatically interpolated), 2) definition of position of papillary muscles in the images, 3) definition of regions in the image subject to fold-over artifacts and 4) definition of the last diastolic time-point. Bulk cardiac motion was corrected for by subtraction of the average in-plane motion. Data points outside the myocardial mask were discarded, and the accepted data points were automatically divided into 32 segments [15]. Both the myocardial mask and segments followed the motion of the LV throughout the cycle. Finally, the velocity vector in each pixel was decomposed into a cardiopolar coordinate system, constituting of in-plane radial and tangential components and a through-plane longitudinal component [5]. The time points corresponding to peak and end-systole was automatically determined from peak global radial motion and minimum LV lumen area, respectively.

The data from the two individual acquisitions were independently processed, including division into segments using the papillary muscles as reference points, and segment-wise averaged as the last step of post-processing. To reduce low-pass filtering of the velocities following signal averaging (due to potential slight variation in heart rate), the data was normalized prior to combination. This was done by temporal stretching the data (using cubic spline interpolation) from one acquisition to the point where maximum correlation in global radial velocity between the acquisitions was achieved.

Myocardial trajectories and strain calculations

Pixel-by-pixel motion paths were calculated from the velocity data through forward-backward-integration-based Fourier tracking [6, 23], with nearest-neighbor interpolation estimating velocities at non-grid locations. Trajectories travelling out of the user-defined myocardial mask were automatically discarded, but no other signal filtering was employed. To include more data points into the analysis (since the in-slice myocardial area is larger in end-systole), the motion tracking was performed twice. This was done by extending the forward-backward motion tracking protocol to calculate motion paths with temporal origin of integration in both end-diastole and end-systole, resulting in two separate descriptions of the displacement field of the myocardium. In both datasets, circumferential strain was calculated in each of the 32 myocardial segments from the trajectories of the two adjacent segments, before the strain waveforms from the two individual trajectory tracings were averaged segment-wise.

The circumferential (Lagrangian) strain in segment s at time t was given by [24]

S c s t = x s 1 t x s + 1 t x s 1 1 x s + 1 1 1
(1)

where x s (t) is the mean in-plane position vector for all pixels in segment s at time t, and s-1 and s+1 refer to the two adjacent segments. As the motion paths were closed, it follows from Eq. 1 that Sc s (1) = Sc s (t ED ) = 0, where t ED is the time point corresponding to end-diastole.

Evaluation of global and regional cardiac function

To investigate global function, peak global radial velocities and circumferential strain (Sc) were calculated in each animal. Furthermore, in order to evaluate regional function, the following parameters were determined:

  • Dispersion of peak motion; the in-slice standard deviation over the 32 segments of the regional radial velocities at peak systole.

  • Dispersion of peak strain; the segment-wise standard deviation of the regional Sc at end-systole.

  • Coherence of motion waveforms: the mean temporal correlation coefficient of regional vs. global radial velocity tracings (as described by Markl et al. [25]).

  • Dispersion of motion waveforms: the standard deviation of the above, over the 32 segments.

  • Index of dyssynchrony: evaluated from cross-correlation delay analysis where the temporal shift in the regional velocity waveforms that maximized the correlation relative to the global motion was calculated (as described by Delfino et al. [26]). The standard deviation of the 32 delays in each animal was used as a single index of myocardial dyssynchrony.

Validation of strain calculations

In order to validate the PC-CMR-derived circumferential strain, the global Sc was compared to echocardiography-derived global Sc. The number of temporal sampling points for a complete cardiac cycle varied between the data sets, due to differences between the techniques and animal heart rate. The dataset with lowest number of sampling points had 23 data points covering the cardiac cycle. To allow temporal paired analysis between the methods, all datasets were re-sampled using cubic spline interpolation to 23 equally spaced time points, and synchronized to peak global Sc.

In addition, as an internal control, the PC-CMR-derived Sc was compared with Sc estimated directly from the segmentation polygons following the borders of the subepi- and subendocardium. Here, the mean global Sc was estimated from

S c g t = 1 2 L epi t L epi 1 L epi 1 + L endo t L endo 1 L endo 1
(2)

where Lepi(t) and Lendo(t) are the lengths, at time t, of the polygons delineating the subepi- and subendocardium, respectively.

Inter- and intrastudy variability

To evaluate inter- and intrastudy variability of the protocol, six animals underwent two PC-CMR examinations on separate days, one of which included two full acquisitions of the same mid-ventricular short-axis slice. Inter- and intrastudy variability in global myocardial velocities and Sc were analyzed using limits-of-agreement. To avoid temporal jitter, all data sets were normalized to end-systole [27].

Statistical analysis

Student’s t-test was used for statistical analysis when comparing the dysfunctional heart to the controls, and p-values < =0.05 were considered statistically significant. Statistical analysis was performed using Matlab. Measurements are presented as mean with standard deviation in parentheses, and correlation coefficients are Pearson’s r.

Results

Animal characteristics

Animal characteristics for the rats included in the validation study are listed in Table 1. During echocardiography, the mean heart rate for the 9 animals in the validation study was 354 (35) bpm. Mean heart rate in all 15 animals during CMR experiments was 371 (31) bmp. Heart rate was not significantly different between groups. On average, the standard deviation of the heart rate in individual animals throughout the CMR examination was 10 bpm.

Table 1 Body and organ weights for the animals included in the validation study

Analysis of myocardial motion

Examples of radial velocities in a representative post-MI heart are shown in Figure 2, and compared to a representative control heart. Distinct alterations in both global (Figure 2a) and regional (Figure 2b-c) motion are evident in the images. The latter also demonstrate the spatiotemporal resolution of the data. Likewise, global and regional Sc are depicted in Figure 3, comparing the same post-MI and control hearts. In Figure 3c and e, the dispersions of regional Sc at peak global Sc are illustrated, that is, the profile of the line marked in Figure 3b and d. In both hearts heterogeneity in regional Sc is evident; however major alterations in the post-MI hearts are clearly visible.

Figure 2
figure 2

Example of radial motion waveforms. The global radial (a) in-slice velocities for two representative animals are shown, one post-MI and one sham. Note the distinct reduction in peak velocities in the diseased heart. Also, spatiotemporally resolved motion maps are displayed as colored plots (b-c) where the y-axis is circumferential position (i.e. segment; direction anteroseptal-anterior-lateral-posterior-posterioseptal), and x-axis is time after r-peak. Green color is positive radial motion (i.e. contraction), red is negative. The altered motion, especially in the anterolateral wall where the infarction is located, is clearly visible. Corresponding CMR magnitude images are shown (d,e), illustrating the location of myocardial thinning in the infarcted heart. The line denotes the first segment and counting direction.

Figure 3
figure 3

Global and regional circumferential strain. For the same animals shown in Figure 2, the global circumferential strain (Sc) is shown here as a function of time (a). Furthermore, regional variation in temporally-resolved Sc is illustrated as color plots (b,d) where y-axis is circumferential position and x-axis is time after r-peak. Finally, the dispersion of Sc at end-systole (peak global Sc) is shown as a function of circumferential position (c,e). Mean and standard deviation of the global Sc at that time point is shown as vertical lines.

Central parameters on myocardial function are listed in Table 2. Peak systolic radial velocity was reduced in the post-MI animals compared to the control (p < 0.001), as was global Sc (p < 0.001). Intragroup variation in peak diastolic radial velocity was larger in the post-MI animals compared to controls, but the mean value was not significantly different between the groups. The dispersion (i.e., the standard deviation over the segments) of regional velocities and Sc in peak systole and end-systole, respectively, was increased in the post-MI animals (p = 0.02 and < 0.001). The regional radial velocity waveforms in the post-MI hearts exhibited lower correlation to the respective global velocities than in the control hearts, and the standard deviation (the spread) of the correlation coefficients was likewise increased in the post-MI hearts (both p < 0.001). The post-MI hearts also exhibited increased index of dyssynchrony, identified from the spread of cross-correlation delays over the circumference (p = 0.006).

Table 2 Selected parameters from PC-CMR acquisition

Validation of strain calculations

Analysis of temporally resolved global Sc demonstrated excellent correlation between PC-CMR and echocardiography data (r = 0.95, p < 0.001; N = 207). A linear fit revealed a close relationship between temporally resolved CMR- and echocardiography-derived data, with a slope not significantly different from 1.00 (95% confidence bounds: [0.97, 1.07], R2 = 0.90), see Figure 4a. Bland-Altman limits-of-agreement was −0.02 ± 3.92%strain (Figure 4b). Intra-animal analysis exhibited likewise strong correlation in all animals (mean r = 0.94 (0.08), p < 0.001 and N = 23 in all nine animals).

Figure 4
figure 4

Comparison of circumferential strain measurements from CMR and echocardiography. Scatterplot (a) and Bland-Altman plot (b) of the data from the validation of PC-CMR-derived global Sc against 2D speckle-tracking echocardiography. Each animal (N = 9) had 23 equally spaced time points covering the complete cardiac cycle, producing a total of N = 207 data points. The data demonstrate small limits-of-agreement with no significant bias.

The PC-CMR-derived Sc also correlated well with Sc estimated directly from the segmentation polygons (r = 0.95, p < 0.001; N = 207), with limits-of-agreement −1.4 ± 4.8%. However, linear fit revealed a relationship whose slope significantly different from 1.00 (95% confidence bounds: [0.76,0.82]). Intra-animal analysis also demonstrated good correlation (mean r = 0.95 (0.05), p < 0.001 and N = 23 in all nine animals).

Intra- and interstudy variability

The resulting intra- and interstudy limits-of-agreement are listed in Table 3. Intra- and interstudy variability in the actual velocity and Sc waveforms from a single animal are illustrated in Figure 5.

Table 3 Intra- and interstudy limits-of-agreement
Figure 5
figure 5

Inter- and intrastudy variability. Inter- and intrastudy variability in two animals (one post-MI and one sham) are shown, illustrating the correspondence between velocity waveforms and magnitude (a,b), and Sc (c,d). All slices consisted of two rotations, and for intrastudy analysis the two individual acquisitions of the two slices were interleaved, attempting to reduce any physiological effects of different time spent under anesthesia.

Discussion

In this study, we have presented a PC-CMR protocol for assessing myocardial motion in rats. Several steps were introduced in data acquisition and post-processing to optimize the protocol. PC-CMR-derived circumferential strain was validated against echocardiography, and we demonstrated that PC-CMR is capable of capturing fine details in the intricate motion of the rodent heart. The velocity and strain data exhibited distinct alterations, both globally and regionally, in the post-MI hearts vs. sham.

Our findings on global Sc agree well with an MR tagging study by Liu and colleagues [28], which reported mid-ventricular Sc in the normal rat heart as −19 (1)%. A study by the same group based on harmonic phase MR tagging [29] investigated the reduction in strain in infarcted hearts. Both studies had 15 temporal frames per cardiac cycle. In a study employing displacement-encoded CMR of infarcted mouse hearts, the circumferential strain one day post-surgery was found to be reduced from −16.4 (1.3)% in controls to −11.6 (1.8) and +4.2 (2.4)% in non-infarcted and infarcted regions, respectively [30]. They also demonstrated good correlation with MR tagging.

Distinct regional alterations in infarcted rat hearts have previously been demonstrated using speckle-tracking echocardiography [31], and our results agree well with their reported circumferential strain. That study also reported, in accordance with our findings, prominent heterogeneity in the Sc in healthy rat hearts. Our results are also in excellent agreement with other studies on global myocardial Sc in rat hearts [32, 33].

A recent study by Dall’Armellina et al. [5] employed PC-CMR in studying myocardial velocities in mice, reporting corresponding findings on myocardial motion. The waveform of the global mid-ventricular radial motion is very similar between the species; however we found peak radial velocities (both in systole and diastole) in normal rats to be roughly the double of what they found in mice. Also, the regional heterogeneity of the motion around the circumference seems more pronounced in normal rats compared to mice. While this might reflect an actual variance between the species, differences could also be attributed to a fewer number of circumferential segments employed in [5].

Validation

PC-CMR has been previously validated as a technique capable of accurately measuring velocity [4, 15], displacement [16, 34] and the deformation gradient of the myocardium [8]. It has been used for calculation of myocardial strain [4, 7] and has been compared to MR tagging (by use of “virtual tagging”) [35]. In our study, two-dimensional speckle-tracking echocardiography (2D-STE) serves as a method for comparison. 2D-STE has been validated against tagged CMR, sonomicrometry and tissue Doppler echocardiography [36], and previously compared to DENSE CMR in mice [37, 38]. Our findings demonstrate that global Sc in a mid-ventricular slice derived from PC-CMR correlates well with echocardiography-derived Sc, with narrow limits-of-agreement.

When comparing to mask-derived Sc, both PC-CMR and echo yielded a linear slope significantly different from 1.00, both results suggesting that the mask-derived Sc systematically overestimated the strain (this is also supported by the fact that the 95% confidence bounds for the linear slope of mask vs. PC-CMR and mask vs echo overlapped (data not shown), suggesting that their slopes were not significantly different). This is not surprising, as the circumferential change-of-length of the subendo- and subepicardium of the LV is expected to be a result of a combination of actual shortening of the LV both circumferentially and longitudinally, which both echo and PC-CMR would account for while the masks would not. Also, the blood volume in the myocardium itself varies throughout the cardiac cycle, contributing to the change in in-slice area of the short-axis LV images, thus not to be attributed to the actual deformation of the cardiomyocytes.

Inter- and intrastudy variability

We found that both intra- and interstudy variability were low, the latter being comparable to previously reported values for human myocardial PC-CMR [39, 40]. A source of variation in our data is suspected to be instability of animal physiology, as the data sets were acquired at different times after induction of anesthesia. Although all datasets were temporally normalized to account for varying lengths of the cardiac cycle, the actual velocities and peak Sc are not independent on heart rate and thus remain uncorrected.

Limitations

Geometrically systematic artifacts as they appear in this study, along with the rotating FOV approach, only apply to Cartesian imaging.

The choice to cover more than one r-r-interval doubles the scanning time since only every second r-peak was used as trigger point. However, this allowed complete coverage of the diastolic phase in combination with black-blood preparation, and improved the reliability of the triggering by allowing decay of currents in the ECG wires. Although the acquisition time per slice in our study was quite long, it is comparable to previous reports on MR-based strain assessment [29, 30]. Compared to echocardiography, the MR examination is considerably more time-consuming. Depending on study design and needs, the trade-off between data yield and acquisition time must be considered.

Although beyond the scope of this study, a direct comparison of PC-CMR-derived strain with MR tagging should be considered, the latter usually being considered as the reference standard for measurement of myocardial strain. While our study validated global strain measurements, future studies should also address comparison of evaluation of regional strain from different methods.

The algorithm for calculating circumferential strain from myocardial trajectories in this paper (Eq. 1) is rather simple compared to more complex approaches, such as spline-based deformation analysis [41]. However, the presented results suggest that our approach is appropriate, yielding accurate and reproducible results.

Since slice selection was done in the laboratory system and the heart moves longitudinally during contraction, slightly different parts of the myocardium may be imaged in different time points throughout the cardiac cycle. This motivated the choice of including end-systole as a temporal origin for estimating tissue trajectories. Compared to conventional forward-backward motion tracking, intra- and interstudy limits-of-agreements were reduced using this extension (Table 3). However, to accurately capture complex three-directional motion and thus true 3D strain, volumetric data is required [42], and should be addressed by future studies. Volumetric PC-CMR might be achieved by embedding velocity encoding gradients into conventional or accelerated 3D imaging protocols, and has been demonstrated to allow comprehensive evaluation of both blood flow and myocardial motion in humans [4345], but not, to our knowledge, in small animals.

Conclusion

In this study, we have presented an optimized PC-CMR protocol allowing assessment of the motion of the myocardium in rats with high detail, and provided a robust method for calculation of regional circumferential strain from the velocity data. By combining optimized slice planning, acquisition parameters and post-processing, exploration of the complex spatiotemporal pattern of in vivo motion and circumferential strain in the healthy and dysfunctional rat heart is feasible. We present, to our knowledge, the first study in small animals using PC-CMR to calculate strain.

Abbreviations

2D-STE:

Two-dimensional speckle tracking echocardiography

DENSE:

Displacement encoded imaging with stimulated echoes

FOV:

Field-of-view

LV:

Left ventricle

MI:

Myocardial infarction

PC-CMR:

Phase contrast cardiovascular magnetic resonance

Sc:

Circumferential strain

SENC:

Strain-encoded CMR.

References

  1. Shehata ML, Cheng S, Osman NF, Bluemke DA, Lima JA: Myocardial tissue tagging with cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2009, 11: 55-10.1186/1532-429X-11-55.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Osman NF, Sampath S, Atalar E, Prince JL: Imaging longitudinal cardiac strain on short-axis images using strain-encoded MRI. Magn Reson Med. 2001, 46: 324-34. 10.1002/mrm.1195.

    Article  CAS  PubMed  Google Scholar 

  3. Aletras AH, Ding S, Balaban RS, Wen H: DENSE: displacement encoding with stimulated echoes in cardiac functional MRI. J Magn Reson. 1999, 137: 247-52. 10.1006/jmre.1998.1676.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Arai AE, Gaither CC, Epstein FH, Balaban RS, Wolff SD: Myocardial velocity gradient imaging by phase contrast MRI with application to regional function in myocardial ischemia. Magn Reson Med. 1999, 42: 98-109. 10.1002/(SICI)1522-2594(199907)42:1<98::AID-MRM14>3.0.CO;2-H.

    Article  CAS  PubMed  Google Scholar 

  5. Dall'armellina E, Jung BA, Lygate CA, Neubauer S, Markl M, Schneider JE: Improved method for quantification of regional cardiac function in mice using phase-contrast MRI. Magn Reson Med. 2012, 67: 541-51. 10.1002/mrm.23022.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Zhu Y, Drangova M, Pelc NJ: Fourier tracking of myocardial motion using cine-PC data. Magn Reson Med. 1996, 35: 471-80. 10.1002/mrm.1910350405.

    Article  CAS  PubMed  Google Scholar 

  7. Delfino JG, Fornwalt BK, Eisner RL, Leon AR, Oshinski JN: Determination of transmural, endocardial, and epicardial radial strain and strain rate from phase contrast MR velocity data. J Magn Reson Imaging. 2008, 27: 522-8. 10.1002/jmri.21211.

    Article  PubMed  Google Scholar 

  8. Zhu Y, Drangova M, Pelc NJ: Estimation of deformation gradient and strain from cine-PC velocity data. IEEE Trans Med Imaging. 1997, 16: 840-51. 10.1109/42.650880.

    Article  CAS  PubMed  Google Scholar 

  9. Herold V, Morchel P, Faber C, Rommel E, Haase A, Jakob PM: In vivo quantitative three-dimensional motion mapping of the murine myocardium with PC-MRI at 17.6 T. Magn Reson Med. 2006, 55: 1058-64. 10.1002/mrm.20866.

    Article  PubMed  Google Scholar 

  10. Pelc NJ, Bernstein MA, Shimakawa A, Glover GH: Encoding strategies for three-direction phase-contrast MR imaging of flow. J Magn Reson Imaging. 1991, 1: 405-13. 10.1002/jmri.1880010404.

    Article  CAS  PubMed  Google Scholar 

  11. Johnson KM, Markl M: Improved SNR in phase contrast velocimetry with five-point balanced flow encoding. Magn Reson Med. 2010, 63: 349-55. 10.1002/mrm.22202.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Bernstein MA, Zhou XJ, Polzin JA, King KF, Ganin A, Pelc NJ, Glover GH: Concomitant gradient terms in phase contrast MR: analysis and correction. Magn Reson Med. 1998, 39: 300-8. 10.1002/mrm.1910390218.

    Article  CAS  PubMed  Google Scholar 

  13. Walker PG, Cranney GB, Scheidegger MB, Waseleski G, Pohost GM, Yoganathan AP: Semiautomated method for noise reduction and background phase error correction in MR phase velocity data. J Magn Reson Imaging. 1993, 3: 521-30. 10.1002/jmri.1880030315.

    Article  CAS  PubMed  Google Scholar 

  14. Rolf MP, Hofman MB, Gatehouse PD, Markenroth-Bloch K, Heymans MW, Ebbers T, Graves MJ, Totman JJ, Werner B, van Rossum AC: Sequence optimization to reduce velocity offsets in cardiovascular magnetic resonance volume flow quantification–a multi-vendor study. J Cardiovasc Magn Reson. 2011, 13: 18-10.1186/1532-429X-13-18.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Espe EK, Aronsen JM, Skrbic B, Skulberg VM, Schneider JE, Sejersted OM, Zhang L, Sjaastad I: Improved MR phase-contrast velocimetry using a novel nine-point balanced motion-encoding scheme with increased robustness to eddy current effects. Magn Reson Med. 2013, 69: 48-61. 10.1002/mrm.24226.

    Article  PubMed  Google Scholar 

  16. Lingamneni A, Hardy PA, Powell KA, Pelc NJ, White RD: Validation of cine phase-contrast MR imaging for motion analysis. J Magn Reson Imaging. 1995, 5: 331-8. 10.1002/jmri.1880050318.

    Article  CAS  PubMed  Google Scholar 

  17. McVeigh ER: MRI of myocardial function: motion tracking techniques. Magn Reson Imaging. 1996, 14: 137-50. 10.1016/0730-725X(95)02009-I.

    Article  CAS  PubMed  Google Scholar 

  18. Sjaastad I, Sejersted OM, Ilebekk A, Bjornerheim R: Echocardiographic criteria for detection of postinfarction congestive heart failure in rats. J Appl Physiol. 2000, 89: 1445-54.

    CAS  PubMed  Google Scholar 

  19. Bauer M, Cheng S, Jain M, Ngoy S, Theodoropoulos C, Trujillo A, Lin FC, Liao R: Echocardiographic speckle-tracking based strain imaging for rapid cardiovascular phenotyping in mice. Circ Res. 2011, 108: 908-16. 10.1161/CIRCRESAHA.110.239574.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Hamilton CA, Jordan JH, Kraft RA, Hundley WG: Noninterleaved velocity encodings for improved temporal and spatial resolution in phase-contrast magnetic resonance imaging. J Comput Assist Tomogr. 2010, 34: 570-4. 10.1097/RCT.0b013e3181d564e8.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Schneider JE, Wiesmann F, Lygate CA, Neubauer S: How to perform an accurate assessment of cardiac function in mice using high-resolution magnetic resonance imaging. J Cardiovasc Magn Reson. 2006, 8: 693-701. 10.1080/10976640600723664.

    Article  PubMed  Google Scholar 

  22. Bernstein MA, Grgic M, Brosnan TJ, Pelc NJ: Reconstructions of phase contrast, phased array multicoil data. Magn Reson Med. 1994, 32: 330-4. 10.1002/mrm.1910320308.

    Article  CAS  PubMed  Google Scholar 

  23. Pelc NJ, Drangova M, Pelc LR, Zhu Y, Noll DC, Bowman BS, Herfkens RJ: Tracking of cyclic motion with phase-contrast cine MR velocity data. J Magn Reson Imaging. 1995, 5: 339-45. 10.1002/jmri.1880050319.

    Article  CAS  PubMed  Google Scholar 

  24. Simpson RM, Keegan J, Firmin DN: MR assessment of regional myocardial mechanics. J Magn Reson Imaging. 2012, 37 (3): 576-99.

    Article  PubMed  Google Scholar 

  25. Markl M, Schneider B, Hennig J: Fast phase contrast cardiac magnetic resonance imaging: improved assessment and analysis of left ventricular wall motion. J Magn Reson Imaging. 2002, 15: 642-53. 10.1002/jmri.10114.

    Article  PubMed  Google Scholar 

  26. Delfino JG, Fornwalt BK, Eisner RL, Leon AR, Oshinski JN: Cross-correlation delay to quantify myocardial dyssynchrony from phase contrast magnetic resonance (PCMR) velocity data. J Magn Reson Imaging. 2008, 28: 1086-91. 10.1002/jmri.21566.

    Article  PubMed  Google Scholar 

  27. Jung B, Foll D, Bottler P, Petersen S, Hennig J, Markl M: Detailed analysis of myocardial motion in volunteers and patients using high-temporal-resolution MR tissue phase mapping. J Magn Reson Imaging. 2006, 24: 1033-9. 10.1002/jmri.20703.

    Article  PubMed  Google Scholar 

  28. Liu W, Ashford MW, Chen J, Watkins MP, Williams TA, Wickline SA, Yu X: MR tagging demonstrates quantitative differences in regional ventricular wall motion in mice, rats, and men. Am J Physiol Heart Circ Physiol. 2006, 291: H2515-21. 10.1152/ajpheart.01016.2005.

    Article  CAS  PubMed  Google Scholar 

  29. Liu W, Chen J, Ji S, Allen JS, Bayly PV, Wickline SA, Yu X: Harmonic phase MR tagging for direct quantification of Lagrangian strain in rat hearts after myocardial infarction. Magn Reson Med. 2004, 52: 1282-90. 10.1002/mrm.20276.

    Article  PubMed  Google Scholar 

  30. Gilson WD, Yang Z, French BA, Epstein FH: Complementary displacement-encoded MRI for contrast-enhanced infarct detection and quantification of myocardial function in mice. Magn Reson Med. 2004, 51: 744-52. 10.1002/mrm.20003.

    Article  PubMed  Google Scholar 

  31. Popovic ZB, Benejam C, Bian J, Mal N, Drinko J, Lee K, Forudi F, Reeg R, Greenberg NL, Thomas JD: Speckle-tracking echocardiography correctly identifies segmental left ventricular dysfunction induced by scarring in a rat model of myocardial infarction. Am J Physiol Heart Circ Physiol. 2007, 292: H2809-16. 10.1152/ajpheart.01176.2006.

    Article  CAS  PubMed  Google Scholar 

  32. Daire JL, Jacob JP, Hyacinthe JN, Croisille P, Montet-Abou K, Richter S, Botsikas D, Lepetit-Coiffe M, Morel D, Vallee JP: Cine and tagged cardiovascular magnetic resonance imaging in normal rat at 1.5 T: a rest and stress study. J Cardiovasc Magn Reson. 2008, 10: 48-10.1186/1532-429X-10-48.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Hyacinthe JN, Ivancevic MK, Daire JL, Vallee JP: Feasibility of complementary spatial modulation of magnetization tagging in the rat heart after manganese injection. NMR Biomed. 2008, 21: 15-21. 10.1002/nbm.1144.

    Article  CAS  PubMed  Google Scholar 

  34. Drangova M, Zhu Y, Bowman B, Pelc NJ: In vitro verification of myocardial motion tracking from phase-contrast velocity data. Magn Reson Imaging. 1998, 16: 863-70. 10.1016/S0730-725X(98)00100-3.

    Article  CAS  PubMed  Google Scholar 

  35. Masood S, Gao J, Yang GZ: Virtual tagging: numerical considerations and phantom validation. IEEE Trans Med Imaging. 2002, 21: 1123-31. 10.1109/TMI.2002.804429.

    Article  PubMed  Google Scholar 

  36. Biswas M, Sudhakar S, Nanda NC, Buckberg G, Pradhan M, Roomi AU, Gorissen W, Houle H: Two- and three-dimensional speckle tracking echocardiography: clinical applications and future directions. Echocardiography. 2013, 30: 88-105. 10.1111/echo.12079.

    Article  PubMed  Google Scholar 

  37. Li Y, Garson CD, Xu Y, Beyers RJ, Epstein FH, French BA, Hossack JA: Quantification and MRI validation of regional contractile dysfunction in mice post myocardial infarction using high resolution ultrasound. Ultrasound Med Biol. 2007, 33: 894-904. 10.1016/j.ultrasmedbio.2006.12.008.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Azam S, Desjardins CL, Schluchter M, Liner A, Stelzer JE, Yu X, Hoit BD: Comparison of velocity vector imaging echocardiography with magnetic resonance imaging in mouse models of cardiomyopathy. Circ Cardiovasc Imaging. 2012, 5: 776-81. 10.1161/CIRCIMAGING.111.972406.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Petersen SE, Jung BA, Wiesmann F, Selvanayagam JB, Francis JM, Hennig J, Neubauer S, Robson MD: Myocardial tissue phase mapping with cine phase-contrast mr imaging: regional wall motion analysis in healthy volunteers. Radiology. 2006, 238: 816-26. 10.1148/radiol.2383041992.

    Article  PubMed  Google Scholar 

  40. Simpson R, Keegan J, Firmin D: Efficient and reproducible high resolution spiral myocardial phase velocity mapping of the entire cardiac cycle. J Cardiovasc Magn Reson. 2013, 15: 34-10.1186/1532-429X-15-34.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Bergvall E, Hedstrom E, Bloch KM, Arheden H, Sparr G: Spline-based cardiac motion tracking using velocity-encoded magnetic resonance imaging. IEEE Trans Med Imaging. 2008, 27: 1045-53.

    Article  PubMed  Google Scholar 

  42. Zhu Y, Pelc NJ: Three-dimensional motion tracking with volumetric phase contrast MR velocity imaging. J Magn Reson Imaging. 1999, 9: 111-8. 10.1002/(SICI)1522-2586(199901)9:1<111::AID-JMRI15>3.0.CO;2-5.

    Article  CAS  PubMed  Google Scholar 

  43. Wigstrom L, Sjoqvist L, Wranne B: Temporally resolved 3D phase-contrast imaging. Magn Reson Med. 1996, 36: 800-3. 10.1002/mrm.1910360521.

    Article  CAS  PubMed  Google Scholar 

  44. Johnson KM, Lum DP, Turski PA, Block WF, Mistretta CA, Wieben O: Improved 3D phase contrast MRI with off-resonance corrected dual echo VIPR. Magn Reson Med. 2008, 60: 1329-36. 10.1002/mrm.21763.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Kvitting JP, Ebbers T, Engvall J, Sutherland GR, Wranne B, Wigstrom L: Three-directional myocardial motion assessed using 3D phase contrast MRI. J Cardiovasc Magn Reson. 2004, 6: 627-36. 10.1081/JCMR-120038692.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by the South-Eastern Norway Regional Health Authority, Anders Jahre’s fund for the Promotion of Science and Familien Blix’ Fond Til Fremme Av Medisinsk Forskning. JES is a Senior British Heart Foundation Basic Science Research Fellow.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Emil KS Espe.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

EKSE was involved in designing the study, collected, analyzed and interpreted CMR data, analyzed and interpreted echocardiography data and drafted the manuscript; JMA was involved in designing the study, provided and performed surgery on experimental animals, collected, analyzed and interpreted echocardiography data and provided critical review of the manuscript; KS and LZ were involved in designing the study, assisted in collection of CMR data, interpreted CMR data and provided critical review of the manuscript; JES was involved in designing the study, interpreted CMR data, and provided critical review of the manuscript; IS conceived and designed the study, collected, analyzed and interpreted echocardiography data, interpreted CMR data, and provided critical review of the manuscript. All authors have read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Espe, E.K., Aronsen, J.M., Skårdal, K. et al. Novel insight into the detailed myocardial motion and deformation of the rodent heart using high-resolution phase contrast cardiovascular magnetic resonance. J Cardiovasc Magn Reson 15, 82 (2013). https://doi.org/10.1186/1532-429X-15-82

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1532-429X-15-82

Keywords