Abstract

Purpose. To evaluate whether 3T clinical MRI with a small-animal coil and gradient-echo (GE) sequence could be used to characterize long-term left ventricular remodelling (LVR) following nonreperfused myocardial infarction (MI) using semi-automatic segmentation software (SASS) in a rat model. Materials and Methods. 5 healthy rats were used to validate left ventricular mass (LVM) measured by MRI with postmortem values. 5 sham and 7 infarcted rats were scanned at 2 and 4 weeks after surgery to allow for functional and structural analysis of the heart. Measurements included ejection fraction (EF), end-diastolic volume (EDV), end-systolic volume (ESV), and LVM. Changes in different regions of the heart were quantified using wall thickness analyses. Results. LVM validation in healthy rats demonstrated high correlation between MR and postmortem values. Functional assessment at 4 weeks after MI revealed considerable reduction in EF, increases in ESV, EDV, and LVM, and contractile dysfunction in infarcted and noninfarcted regions. Conclusion. Clinical 3T MRI with a small animal coil and GE sequence generated images in a rat heart with adequate signal-to-noise ratio (SNR) for successful semiautomatic segmentation to accurately and rapidly evaluate long-term LVR after MI.

1. Introduction

Assessment of ventricular function and mass is an important guide for predicting survival in patients suffering from different cardiovascular diseases, such as ischemic heart disease [1]. In addition, the efficacy of long-term drug therapy of compounds such as angiotensin-converting enzyme (ACE) inhibitor and beta-adrenergic blocking agents are evaluated by assessment of left ventricular (LV) mass and function [2]. Rat models have been widely used to study LV performance and remodelling after induction of myocardial infarction (MI) [3, 4]. In addition, the effects of ACE inhibitors and captopril have been examined in this way [1, 5]. Hence, noninvasive methods to accurately quantify ventricular mass and function are invaluable to clinical and preclinical research.

Cardiac magnetic resonance imaging has been considered the optimal tool for noninvasive evaluation of LV function and mass due to its ability to acquire images with high spatial resolution and superb soft tissue contrast without risk of radiation [68]. Also, it allows for noninvasive imaging on the same animal over a prolonged period of time, providing subject-specific longitudinal data and reducing the number of animals needed.

Rodent (rat and mouse) models are preferred over large animal models because of their small size, easy availability and low maintenance cost. Furthermore, transgenic and gene-knockout technologies allow design and generation of rodents that mimic human cardiovascular disease [9]. However, rat models of heart failure are widely and successfully used and are preferred over mouse models because of reasons explained by Patten and Hall-Porter [10]: it is easier to perform open-chest surgical procedures and invasive hemodynamic assessments, and a greater number of postmortem histological or molecular biological analyses can be done given the approximately 10-fold greater myocardial mass of rats compared with mice.

Imaging rats is challenging due to their small heart size and rapid heart rates [11]. Thus, most of these studies are conducted in high-field nonclinical magnetic resonance (MR) scanners with high performance gradient systems and small RF coils for rapid and high signal-to-noise ratio (SNR) imaging. Several articles have reported assessment of LV function in rats and mice using dedicated high-field MR scanners (≥4.7 Tesla) [1217]. However, these scanners are not available at many institutions [18]. Although cardiac imaging in rodents using clinical grade scanners has been implemented using optimized radio frequency (RF) coils, the images produced suffer from reduced signal and contrast due to lower field strength [17]. Consequently, high-resolution images (<250 μm in-plane resolution) acquired in clinical scanners with gradient echo (GE) sequence generally result in poor definition of the endocardial and epicardial boundaries, requiring manual intervention during the segmentation process [19, 20].

The purpose of this study was to assess whether 3T clinical MRI with a small animal coil and GE sequence could generate images with sufficient contrast and resolution to characterize using semi-automatic segmentation software (SASS), long-term LV remodelling (LVR), and changes in functional parameters—such as LV ejection fraction (EF), end-diastolic (ED), and end-systolic (ES) volume—following induction of nonreperfused MI in rats.

2. Materials and Methods

2.1. Magnet and Gradient System

All scans were performed on a MAGNETOM Allegra 3T MRI scanner (Siemens, Erlangen, Germany). This system is equipped with maximum gradient strength of 40 mT·m−1 and slew rate of 400 mT·m−1·ms−1. A custom built, circularly polarized, 70 mm (inner diameter) birdcage radiofrequency (RF) coil was used to transmit and receive the signal.

2.2. Animal Model

All animal studies were approved by the Faculty of Health Sciences Animal Research Ethics Committee (REC) and performed in accordance with the National Institutes of Health (NIH, Bethesda, MD, USA) guidelines. 17 male outbred Wistar albino rats were used for the study that weighed 1 9 5 . 3 ± 1 0 . 8  g (standard deviation) at implant. The animals were split into three groups, namely, healthy/control ( 𝑛 = 5 ), sham ( 𝑛 = 5 ), and MI ( 𝑛 = 7 ). In the MI group, a nonreperfused MI was induced as previously described by Dobner et al. [21]. Briefly, rats were placed in an induction chamber and anesthetized with air flow containing 5% isoflurane and 1.5% Oxygen. Animals were then intubated with a 16G catheter and mechanically ventilated with 1.5–2% isoflurane and 0.3% oxygen. A left intercostal thoracotomy was performed and MI was induced by permanent ligation of the left anterior descending (LAD) coronary artery. A 6-0 prolene suture was used to ligate and placed 2 to 3 mm below the left atrial appendix. Fibrillation and ischemia of the LV wall confirmed the infarction. Sham animals were subjected to the same procedure without ligation of LAD coronary artery. Following the surgical procedure, an analgesic dose of 100 μL Temgesic was administered by intramuscular injection to the hind leg. Animals were monitored daily for the duration of the study. Analgesic (100 μL Temgesic) was administered every morning for 48 hrs after surgery. All animals survived surgery.

3. MR Image Acquisition

3.1. Animal Preparation

Baseline imaging using healthy animals ( 𝑛 = 5 ) was performed with a GE fast low-angle shot (FLASH) sequence [4, 17]. Infarcted and sham rats were imaged 2 and 4 weeks after surgery. Before scanning, anaesthesia was induced and maintained through a nose cone containing an atmosphere of 1.5% isoflurane in air. Subdermal electrocardiogram (ECG) needle electrodes were attached on the right front and left rear legs. Respiration signal was superimposed on the measured ECG signal by closely coupling one of the ECG lead wires to the rat’s abdomen. A rectal temperature probe was used to measure the rat’s core body temperature, which was maintained by circulating hot water through a tube in the rat bed. ECG, respiration, and core body temperature were monitored by a dedicated small animal MR-compatible monitoring and gating system (SA Instruments, Stony Brook, NY, USA).

3.2. MR Parameters

MR parameters for short-axis images were field of view (FOV) = 60 mm, matrix size = 256 × 256, in-plane resolution 234 μm, slice thickness 1.0 mm, TR 7.7 ms, TE 3.5 ms, flip angle 25° to 40°, and 3 signal averages. For the four-week scan, FOV and matrix size were increased to 75 mm and 320 × 320, respectively. In all cases, the FLASH sequence was triggered by both ECG and respiration. One k-space line was acquired per frame resulting in 18 to 30 frames per cardiac cycle depending on the heart rate. Three averages were implemented to yield adequate SNR. Acquisition time was approximately 2.5 minutes per slice, and 15 to 16 slices were typically required to cover the entire LV. One to two additional slices were acquired to cover the entire LV of the MI group at the 4-week scan due to animal growth.

4. Histological Preparation and Infarct Size Measurement

A day after the completion of MRI scanning, all animals were euthanized with intraperitoneally injected sodium pentobarbitone whilst under isoflurane-induced anaesthesia. Hearts excised from healthy rats were trimmed free of right ventricle and atria to determine true LV mass. In the MI group, infarct size was determined histologically as described by Dobner et al. [21]. Briefly, the hearts were injected with 1 mL ice-cold saturated KCl solution after anaesthesia to arrest them in diastolic state and then excised for weighing. 3 μm sections were obtained from 4 equal divisions of every heart. Each section was stained with Masson’s trichrome to determine the extent of infarction and images captured with a Nikon E1000 M camera (Nikon, Tokyo, Japan) with a 0.5x magnification lens. Infarct size was determined using the midline method as previously described by Takagawa et al. [22]. Briefly, infarct size was measured as the sum of midline length of infarcted myocardium from all sections divided by the sum of midline LV circumferences from all sections multiplied by 100. Measurements were acquired with Visiopharm Integrator Systems (Visiopharm, Hørsholm, Denmark).

5. MR Image Analysis

The image stacks with largest LV cavity were used as the end-diastolic (ED) images and the smallest as end-systolic (ES) images. From ED and ES short-axis images, LV borders were segmented using SASS, called Segment (http://segment.heiberg.se/) without any long-axis correction. This software has been shown to yield low interobserver and intraobserver variability (<5%) in volumetric analyses [23]. We implemented the LV segmentation technique and definitions for the most apical and basal slice that were previously explained by Riegler et al. [23]. Briefly, the most apical slice is defined as the first slice with no blood pool while the most basal slice is the one that shows a remaining part of the thick myocardium and is below the aortic valve. Papillary muscles were segmented as part of the myocardium (Figures 1(a) and 1(b)).

The segmentation method employed in Segment software is based on the concept of deformable models, but extended with an enhanced and fast edge detection scheme [24]. Analyses using Segment software was initiated using the “left ventricular tool” algorithm. First, a marker was positioned in the center of the LV blood pool, from which the algorithm fits a deformable model that expands toward the endocardial border to include the LV blood pool across the various short-axis slices. Afterward, manual corrections were performed as necessary for the epicardial border, and for the endocardial border where the papillary muscles appeared separated from the myocardium. If the myocardium moved out of the plane due to long-axis contraction, a straight line was drawn to merge the two visible ends of the myocardium (Figure 1(c)), as described by Barbier et al. [25].

ED volume (EDV), ES volume (ESV), and LV mass (LVM) were measured by SASS. Other functional parameters were calculated as follows: stroke volume (SV) = EDV − ESV, EF = 100 × SV/EDV, and cardiac output (CO) = SV × Heart Rate.

LV wall thickness (WTh) was measured 4 weeks after surgery using the same software. WTh analyses were performed on midventricular slices as explained by Ross et al. [15] and Ojha et al. [26]. Briefly, following the delineation of epicardial and endocardial borders, LV wall was divided into 16 equal sectors using the right anterior ventricular insertion as a reference point. 4 contiguous measurements were averaged to form 4 equal sectors of 90° each. These four sectors were called anterior, lateral, posterior, and septal. Papillary muscles were not included in this analysis. Subsequent to wall thickness measurements, % systolic wall thickening (%WTh) was calculated as 100 × (systolic wall thickness − diastolic wall thickness)/diastolic wall thickness.

6. Statistical Analysis

If not indicated otherwise, all data are specified as mean ± standard error of the mean (SEM). Correlation and Bland-Altman analyses [27] were used to examine the relation and agreement between MR LVM and post-mortem measurements. An unpaired two-tailed t-test was used to compare LV function and mass between sham and MI groups. 𝑃 < 0 . 0 5 was considered significant.

7. Results

All 17 animals survived the duration of the study. LVM mass validation showed very close correlation between the two methods. Rats in the MI group underwent considerable LVR, such as a significant drop in EF and increase in LVM, EDV, and ESV as highlighted below.

7.1. LV Mass Authentication

Healthy rats were used to validate MR LV mass with the true wet LV mass found by measurement on an electronic weight balance a day after baseline imaging. The mean MR LV mass was 5 1 2 . 4 ± 3 1 . 8  mg, while the mean wet LV mass was 5 0 9 . 6 ± 2 8  mg, correlation 0.87. Bland-Altman analysis revealed a mean difference of 2 mg between the two methods with 95% confidence interval of 72 mg (Figure 2).

7.2. Structural and Global Functional Analysis

Short axis images were analysed with Segment. The midventricular images (Figure 3) at ES phase and histological images (Figure 4) acquired from sham and infarcted rats illustrate LV dilation and wall thinning in the anterior and lateral regions of the heart. Dilation of ESV is evident by 2 and 4 weeks after surgery demonstrating contractile dysfunction (Figure 3). From post-mortem measurements, infarct size was found to be 2 7 . 2 ± 2 . 6 %. Table 1 describes cardiac parameters of sham and MI animals after surgery.

At 2 weeks after surgery, EDV did not increase significantly (Figure 5(a)), while there was significant increase in ESV in the MI group (Table 1, Figure 5(b)) compared to the sham group ( 𝑃 < 0 . 0 5 ). There was no significant change in CO and SV at 2 and 4 weeks after surgery. At 4 weeks after surgery, there was significant enlargement of EDV and ESV in the MI rats compared to sham rats at 4 weeks aftert surgery ( 𝑃 < 0 . 0 5 , Table 1). For the MI group, EF dropped significantly at 2 weeks and was even lower at 4 weeks after surgery (Table 1, Figure 5(c)). LV mass was significantly greater in MI rats compared to sham rats both at 2 and 4 weeks after surgery (Table 1, Figure 5(d)).

7.3. Regional Functional Analysis

The MI group revealed dyskinesia and severe wall thinning in the infarcted area (Figure 3) at 4 weeks after surgery. In order to quantify this, we performed WTh analyses for a midventricular slice 4 weeks after surgery in the MI and sham rats. At ED in the MI group, there was a significant drop in wall thickness in the anterior region (from 1 . 7 ± 0 . 1  mm to 1 . 3 ± 0 . 1  mm, 𝑃 < 0 . 0 5 ), and a modest decrease in the lateral region (Figure 6(a)). In addition, in the posterior region there was significant increase in wall thickness (from 1 . 6 ± 0 . 0 4  mm to 2 . 0 ± 0 . 1  mm, 𝑃 < 0 . 0 5 ) while there was modest increase in the septal region. Increases in the posterior and septal regions indicate the presence of hypertrophy. At ES in the MI group, there was a significant drop in wall thickness in the anterior (from 2 . 8 ± 0 . 2  mm to 1 . 9 ± 0 . 2  mm, 𝑃 < 0 . 0 5 ) and lateral regions (from 3 . 3 ± 0 . 1  mm to 2 . 2 ± 0 . 3  mm, 𝑃 < 0 . 0 5 , Figure 6(b)), while there was a modest decrease in wall thickness in lateral and septal regions. Histological analysis of the MI group revealed complete absence of scar tissue in the septal region (Figure 4(a)). In contrast, most of the anterior and lateral quadrants of the heart demonstrate the presence of scar tissue and severe wall thinning (Figure 4(a)).These observations were in agreement with WTh analysis as these two quadrants suffered severe contractile dysfunction as a result of MI, as evident from the significantly reduced wall thickening during systole (Figure 6(b)).

8. Discussion

This study demonstrates serial assessment of rat heart, structurally and functionally, after nonreperfused MI using a 3T clinical MR scanner. The characteristics of LVR, such as LV dilation, hypertrophy, and wall thinning, were evident on these images. The implementation of the MR sequence parameters and a small animal coil yielded images with sufficient SNR and contrast to delineate major anatomical landmarks, such as papillary muscles. In addition, these parameters produced images with excellent spatial resolution within reasonable scan times. This level of the resolution made it possible to perform accurate quantitative WTh analyses, thereby providing extra information and better understanding of LVR in the rat heart.

Due to very low temporal resolution (39 ms), Franco et al. [20] performed multiple acquisitions to determine trigger delay to obtain ED and ES frames for functional analysis. However, our temporal resolution (7.7 ms) was higher than previously reported values with clinical and non-clinical scanners [26, 2830]. This resolution resulted in large number of frames per cardiac cycle which enabled us to obtain accurately ED and ES frames without the need of multiple acquisitions, resulting in reduction in scan time, and accurate global and regional functional analyses [20]. Furthermore, such high temporal resolution is necessary for depiction of LV dynamics in rats, especially with the heart rates of 300 beats per minute or above [31].

In our study, the results of LVM measured by MRI correlated well with post-mortem values. In addition, Bland-Altman analysis revealed very low mean difference and low variation (<15% of mean of true wet LVM) between the two methods. Moreover, these results are in close agreement with previous studies [16, 31]. MRI has been considered the gold standard for accurate measurements of LVM in small animal models in non-clinical scanners [17, 30, 32]. However, this study shows that accurate and reliable volumetric and mass measurements can also be performed on data acquired with a 3T clinical scanner.

The contrast between LV blood pool and myocardium was large enough to effectively use SASS for segmentation to evaluate cardiac parameters. We compared our EF values, a major determinant of cardiac function, with reported values in the manner explained by Riegler et al. [23] and Delattre at al. [28]. In our study, EF in healthy animals ( 7 1 . 8 ± 1 . 3 %) was equivalent to reported values [4, 33]. Similarly, EF in sham rats showed slight variation but was equivalent to values reported by Nahrendorf et al. [34]. EF in the MI group at 4 weeks was similar to the value reported by Carr et al. [35], but was higher than those reported by Stuckey et al. [36] and Nahrendorf et al. [4]. The variation in EF in infarcted rats most likely results from variation in infarct size, leading to differences in level of contractility [1].

During segmentation, some previous studies considered papillary muscles as part of the LV blood pool (not part of myocardium) [7, 37, 38]. This approach has been adopted due to its simplicity, as the shape of the LV blood pool is made much simpler, thus making the evaluation of LV volume faster and more time efficient [39]. However, omission of papillary muscles causes significant differences in measurements of EDV, ESV, and LVM, but not in EF [39, 40]. Nonetheless, possible effects of omitting papillary muscles should be taken into account when performing data analysis, especially during estimation of ESV [23]. In our study, papillary muscles were treated as part of myocardium because they contribute considerably to the blood pool at end diastole [40]. In addition, strong correlation between post-mortem and MRI LVM can be achieved with the inclusion of these muscles in the LV myocardium [31, 40, 41].

The analysis time in our study ( 2 5 ± 5 minutes) was longer than previously reported values (10 and 8 minutes) using other SASS [42, 43]. Nevertheless, these studies conducted cardiac analysis in mouse heart, thus the number of slices covering the whole heart with the same slice thickness was lower. Time savings for SASS are enormous, compared to manual segmentation. It can range between 20% and 70% depending on the number of cardiac frames segmented manually [23]. In our study, the segmentation time was significantly shorter than reported values using manual segmentation software (53 minutes and 2.5 hours) [23, 42], which led to fast and accurate volumetric analysis. The time savings allowed us to perform longitudinal studies and analyze cardiac parameters over time, such as EDV, ESV, EF, SV, and CO. Also, it granted us the opportunity to perform regional functional analysis. Similar to a previous study [23], a significant amount of time was spent manually adjusting the contour line defining the epicardial border. Therefore, contrast enhancement between myocardium and surrounding tissues may lead to greater time savings.

In the sham group there was no significant variation in SV, CO, and LVM at either 2 or 4 weeks (Table 1). In the MI group, rats underwent LVR, in terms of functional and structural changes (Table 1), after LAD occlusion. The major functional change was a significant drop in EF due to the large extent of the fibrous tissue (loss of myocardium) causing contractile dysfunction in the heart (Figures 4(a) and 6(c)) [3, 44]. The contractile dysfunction also caused significant increase in ESV (Figure 5(b)) at 2 weeks. Furthermore, there was a slight increase in EDV at 2 weeks, which would result in maintenance of SV [1]. The increase in ventricular volumes would likely amplify diastolic and systolic wall stress which can serve as a stimulus for further ventricular dilation [1]. Thus, further significant increase in ESV and EDV was observed 4 weeks after surgery. LV hypertrophy was evident in infarcted rats as indicated by severe increase in heart weight and LVM at both 2 and 4 weeks, probably due to progressive LV dilation (Table 1 and Figure 5(d)) [4]. In parallel with increase in LVM, LV hypertrophy was also evident in posterior and septal regions (Figure 6(a)) through WTh analysis. This was also observed by Litwin et al. [45] and Rubin et al. [46] due to increase in myocyte cell volume in attempts to stabilise contractile function and reduce progressive dilation [47].

In the MI group, collagen scarring was evident in anterior and lateral regions (blue-stained region) of the heart 4 weeks after surgery (Figure 4(a)). By contrast, rats in the sham group had no scarring (Figures 4(c) and 4(d)). This observation confirms that infarction was induced only by ligation of the LAD coronary artery in the MI group. In regions with infarction, fibroblasts are transformed into myofibroblasts which are responsible for collagen formation to stabilise distending forces and prevent further deformation [47]. By 3 to 4 weeks, the infarcted region is completely healed by formation of fibrous tissue (collagen scar) (Figure 4(b)) [46, 48]. Regional functional analysis revealed significant reduction in wall thickness at ED in the anterior region and a modest decrease in the lateral region (Figure 6(a)). This was a result of necrosis of myocytes in infarcted regions [1]. At ES, WTh analysis revealed significant reduction in anterior and lateral regions and modest reduction in the remaining regions of the heart (Figure 6(b)). Furthermore, %WTh revealed reduced wall thickening in all regions, which is significant only in lateral and posterior regions (Figure 6(c)). These reductions in ES wall thickness and wall thickening arise due to contractile dysfunction, not only, in infarcted regions but also in non-infarcted regions [16, 49]. However, anterior and lateral regions showed a greater deficit in function (Figures 6(a) and 6(b)) than other regions because of considerable loss of viable myocardium (Figure 4(b)). This loss also strongly affects functional parameters as explained above.

In summary, we have demonstrated assessment of a rat model of nonreperfused MI using a 3T clinical MRI and SASS, called Segment. Following the induction of MI, considerable changes in functional and structural parameters were found, such as a decrease in EF and increases in EDV, ESV, and LVM which were quantified effectively with SASS in a short period of time. These results are in good agreement with previous studies conducted with non-clinical scanners [4, 33, 34]. Furthermore, the high in-plane resolution of our images (234 μm) enabled us to perform WTh analysis revealing important relationships between LV function and LVR.

In this study, ECG and respiratory triggering with 3 signal averages were necessary to acquire high-quality images. Some images, occasionally, suffered when rats had arrhythmic heartbeats. Due to long TE, some flow artefacts were observed, however those artefacts were far away from ED and ES frames, allowing robust and reliable functional and structural assessment of the rat heart. Hardware limitations and dephasing gradient lobes in the FLASH sequence limit the achievable reduction in TE, but implementation of spiral imaging can eradicate flow artefacts [50]. A large flip angle was used to provide high contrast between endocardium and the LV blood pool. Unfortunately, the epicardium did not benefit from this contrast phenomenon causing ambiguity in delineating the epicardial border. Riegler et al. [23] and Ross et al. [15] reported similar limitations due to the very small contrast difference between epicardium and surrounding tissue, and susceptibility differences between lung and the heart due to magnetic field inhomogeneities.

9. Conclusion

Clinical 3T MRI with a small animal coil and GE sequence generated images of rat heart with adequate SNR for successful semi-automatic segmentation to accurately and rapidly evaluate long-term LVR following nonreperfused MI.

Acknowledgments

The authors would like to thank Fulufhelo Masithulela for assistance with handling of the animals, Dr. Vijay Dahya for advice and guidance on image segmentation, Charles Harris for construction of the rat bed, Mogammad Fakier for construction of the automatic temperature regulation system, and Dr. Lawrence Wald and Azma Mareyam from the Martinos Center for Biomedical Imaging at Massachusetts General Hospital for guidance in building the small animal coil. This study was supported by the South African Research Chairs Initiative of the Department of Science and Technology and National Research Foundation of South Africa, the Medical Research Council of South Africa, and NIH Grant no. R21AA017410.

Supplementary Materials

1. 4_weeks_infarcted_rat_4_chamber.avi: A movie showing long-axis view of infarcted rat heart 4 weeks after surgery. Wall thinning is observed at the anterior segment of the heart.

2. 4_weeks_infarcted_rat_shortaxis.avi: A movie showing short-axis view of infarcted rat heart 4 weeks after surgery. Wall thinning is observed at the anterior segment of the heart.

3. 4_weeks_sham_rat_4_chamber.avi: A movie showing long-axis view of sham rat.

4. 4_weeks_sham_rat_shortaxis.avi: A movie showing short-axis view of sham rat.

  1. Supplementary Videos