|
|
|||||||||
Basic Science Investigations |
1 Department of Nuclear Medicine, University Hospital Münster, Münster, Germany
2 Institute of Heart and Circulation Physiology, University Hospital Düsseldorf, Düsseldorf, Germany
3 Department of Nuclear Medicine, Guys & St Thomas Hospitals, London, United Kingdom
4 Institute of Pathophysiology, Center for Internal Medicine, University Hospital Essen, Essen, Germany
5 Center for Transgene Technology and Gene Therapy, Flanders Interuniversity Institute for Biotechnology, University of Leuven, Leuven, Belgium
| ABSTRACT |
|---|
|
|
|---|
Key Words: PET heart failure segmentation mice
| INTRODUCTION |
|---|
|
|
|---|
High-resolution small-animal PET promises to extend PET methodology to small animals such as rats and mice. These animals are being increasingly used to study genetic defects and their phenotypic expression. Measurement of myocardial gene expression and measurement of their phenotypic expression in cardiac parameters such as perfusion and metabolism have been recently reported for rats and mice (11,12).
The evaluation of geometric and functional parameters such as ventricular volumes and ejection fraction with PET has been published for rats (13), but corresponding data for mice are largely absent. Their much smaller heart size makes noninvasive imaging by PET more challenging. It requires equipment capable of obtaining excellent resolution. The 32-module quadHIDAC PET scanner (Oxford Positron Systems) satisfies this condition with an achievable resolution of <1 mm using a resolution recovery reconstruction technique (1416).
The measurement of parameters of cardiac function such as left ventricular (LV) dilation in addition to biochemical properties of the myocardium within the same PET scan is desirable because it obviates the need for a second imaging modality such as echocardiography or cardiac MRI. This approach reduces study complexity and, therefore, facilitates time- and cost-effective serial measurements.
In this work, the feasibility and accuracy of the assessment of LV dilation using a semiautomatic contour-detection algorithm working on nongated 18F-FDG PET images were investigated. The contour-detection algorithm has previously been successfully applied to PET and SPECT in humans (17,18) and was adapted to small-animal PET. The assessment was based on the comparison with MRI in 20 mice (10 wild-type, 10 transgenic), which underwent both PET and MRI. In Figure 1, a short-axis slice through the left ventricle of a wild-type mouse (Fig. 1A) and a transgenic mouse (Fig. 1B), as obtained by PET and MRI, respectively, is shown as an illustration of the 2 imaging modalities.
|
| MATERIALS AND METHODS |
|---|
|
|
|---|
-myosin heavy chain (courtesy of Michael D. Schneider, Baylor College of Medicine, Houston, TX). On inspection, knockout mice were indistinguishable from their nonknockout littermates.
|
For data acquisition, the mouse was placed on a heating pad to maintain a body temperature within the normal range. Cardiac excitation and respiration were monitored throughout the scan to exclude anomalies. All vital parameters were measured with a dedicated system (BioVet; Spin Systems Pty Ltd.). The temperature was measured rectally, electrocardiogram (ECG) electrodes were placed on both forepaws and the left hindpaw, and respiration was measured with a small pressure detector attached to the mouse. The animals were anesthetized by inhalation of isoflurane (1.5%) and intravenously injected with approximately 10 MBq 18F-FDG in 100 µL of saline solution.
One hour after injection data acquisition was started. List-mode data were acquired for 15 min and subsequently reconstructed into an image volume of 110 x 60 x 20 mm3, voxel size of 0.4 x 0.4 x 0.4 mm3, using a resolution recovery reconstruction algorithm (15). A resolution of 0.7 mm is achieved (full width at half maximum).
MRI.
MRI was performed using a Bruker DRX 9.4-T wide-bore NMR spectrometer equipped with an actively shielded 40-mm gradient set (capable of 1 T/m maximum gradient strength and 110-µs rise time at 100% gradient switching) and a 30-mm birdcage resonator (1H). Mice were anesthetized by inhalation of isoflurane (1.5%). High-resolution images of mouse hearts were acquired using an ECG-triggered fast gradient echo cine sequence with a flip angle of <45°. To obtain a maximum of images over the entire cardiac cycle (
100 ms at a heart frequency of 600 beats/min), the shortest possible echo time (1.8 ms) was applied. Using a repetition rate of 3.8 ms allowed for the acquisition of 26 frames during one RR interval, which leads to a sufficient temporal resolution to accurately determine diastole and systole within the cardiac cycle. The resulting pixel size was 117 x 117 µm2 (field of view, 30 x 30 mm2; 256 x 256 matrix). The total acquisition time per slice for one cine sequence was
2 min (22). After acquisition of scout scans, one short-axis slice at the height of the papillary muscle was acquired in cine mode as the basis for volume analysis.
Data Processing
PET.
Data were transferred to a personal computer with a 900-MHz AMD Athlon central processing unit (Advanced Micro Devices Inc.) running the Linux operating system for further analysis. Processing of PET images to obtain a measure for LV volumes was performed using in-house software programmed in MATLAB (The Mathworks) and C programming languages. The interactive preprocessing steps as well as the automatic contour-detection algorithm are described elsewhere (17,18).
In brief, the left ventricle is interactively aligned along its long axis and the valve plane is specified perpendicularly. Image data are automatically transformed from the original cartesian matrix to heart coordinates, a combination of spheric polar coordinates toward the apical portion and circular cylindric coordinates in the midventricular and basal portions of the left ventricle. The midmyocardial contour is mathematically represented in 3 dimensions as an elastic membrane. Internal forces inherent in the membrane model ensure a continuous contour and allow bridging the gaps in the tracer distributionfor example, in the presence of a myocardial scar. The distribution of tracer activity within the PET image data is the source of an attractive field of force acting on the membrane. Because the tracer accumulation is high in the myocardium, the midmyocardial contour is locally attracted toward the accumulated tracer in the myocardium. The shape of the midmyocardial contour umid is modeled as the solution to a second-order partial differential equation, where (x, y) denotes the myocardial position in heart coordinates, umid denotes the radial displacement of the midmyocardial contour, and F denotes the external force exerted on the contour:
![]() |
![]() |
fact is the activity distribution and g is a nonlinear weighting function dependent on the distance between contour umid and radial position z, where z is relative to the origin of the spheric polar coordinate system in the apical portion and to the long axis in the midventricular and basal portions of the left ventricle. a and b are scaling parameters:
![]() |
The problem is discretized using 24 azimuth steps of equal angle, 16 polar/axial steps, and 96 radial steps and solved numerically using a successive overrelaxation (SOR) algorithm.
The volume enclosed by the calculated midmyocardial contour (Fig. 2) is taken as a measure for LV size and denoted as LV volume index (LV-VIPET). Because PET data were acquired in static mode and LV contraction was consequently ignored, the obtained values are only indicative of LV size and are not to be mistaken for the LV volume.
|
Because PET data were acquired in static mode, LV-VIPET is comparable to neither the end-diastolic nor the end-systolic LV volumes measured by MRI. Therefore, a new parameter LV-VIMRI is calculated in the following way:
![]() |
Comparison of PET and MRI.
The agreement of PET and MRI measurements was analyzed by calculating correlation and linear regression as well as by using a graphical analysis method according to Bland and Altman (24). The fractional contributions xED and xES were varied to determine whether the choice for xED and xES significantly influences the result.
Additionally, as an example of application, the development of LV size with age in both wild-type and transgenic mice was analyzed with PET and MRI. Differences in behavior between the 2 groups of animals (wild-type vs. transgenic) were demonstrated. All calculations were performed twice, xED set to 0.5 and 0.8 (xES set to 0.5 and 0.2), again to exclude a relevant dependency on the particular choice of xED/xES.
Statistical analysis was performed using standard statistical software (JMP 5.1; SAS Institute Inc.). The groups of animals are described by their mean and SDs. Comparisons between groups were performed by the 2-sided Student t test (level of significance, 0.95).
| RESULTS |
|---|
|
|
|---|
Correlation
Measures for LV dilation (LV-VIPET/MRI) obtained by PET and MRI correlated well. Variation of the fractional contributions xED and xES in the calculation of LV-VIMRI had only little influence on the result (Table 2) with R = 0.89 for a wide range of values for xED and xES. The intersect of the correlation line naturally differed; the slope remained nearly constant. For display in the figures, xED and xES were set to 0.5. A scatter plot of this correlation is shown in Figure 3A. The linear regression (Table 2) shows that the PET measure for LV dilation (LV-VIPET) underestimates the corresponding measure for MRI (LV-VIMRI) for all sensible estimates of xED/xES. Graphical analysis according to Bland and Altman (24) does not reveal a dependency on volume (Fig. 3B). The data support the assumption of a linear relationship between values obtained by PET and MRI.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
The results of this study show a good correlation of LV dilation measured by PET and MRI. Results were consistent over a large range of LV volumes; a systematic variation with LV volume was not observed. Whereas measurement of cardiac geometry and functional parameters have been reported for PET in rats (13), the results published for PET in mice are lacking, likely because of the fact that accurate imaging of mouse hearts requires a much higher resolution due to the difference of cardiac dimensions between rats and mice. In the future, one would expect to see more PET studies of murine hearts as more high-resolution PET systems are installed.
Despite a good correlation between PET and MRI, the values for LV dilation obtained by PET were consistently lower than the corresponding values for MRI (mean for end-diastolic or end-systolic states). Obviously, the PET methodology in this work yields instead an index of LV dilation rather than the real end-diastolic or end-systolic volumes. This fact is reflected in the choice of terminology (LV-VI) we have used in the study. Because PET in this study relied on static nongated images, whereas MRI did not, identical measurements between the 2 modalities were not expected. In static PET, myocardial activity over different contraction states is integrated, thus rendering the method unsuitable to assess LV contractionthat is, in the presence of severe myocardial ischemia. Therefore, the objective of this study was limited to the investigation of whether PET is indicative of the presence of LV dilation in groups of transgenic mice with nonischemic cardiomyopathies. This objective can be achieved largely independent of the fractional contributions xED and xES of the end-diastolic and end-systolic volumes assumed in the calculation of LV-VIMRI, as the results clearly illustrate. The possible future usage of gated PET acquisition could expand the versatility of the methodology to also cover the measurement of LV volumes during the cardiac cycle and the assessment of global and regional LV function, albeit at the expense of a higher methodologic complexity. Also, correction of respiratory motion in PETwhich, unlike in MRI, is difficult to avoid because of longer acquisition intervalsshould be a goal in future development.
In our study, the midmyocardial contour has been chosen for analysis rather than the endocardial or epicardial contours because it is better defined and is less dependent on resolution and wall thickening. It is worth mentioning that the described contour-detection algorithm can be easily extended to the additional calculation of the endocardial and epicardial contours (18), possibly rendering it suitable for the calculation of the LV volumes and the ejection fraction on the basis of ECG-gated data acquisition. However, the midmyocardial contour is also affected by the limited resolution of PET. Volumes within the midmyocardial contour tend to be smaller in images obtained by imaging modalities with a lesser resolution such as PET when compared with higher-resolution imaging modalities such as MRI. This is an effect of convolution of the "real" myocardial signal (which is spatially curved) with the scanners point-spread-function, resulting in a shift of the midmyocardial contour toward the convex side in the reconstructed PET images.
Another limitation of the study is the usage of the method of Teichholz et al. (23) applied to one short-axis slice instead of a 3-dimensional reconstruction for the calculation of LV volumes in MRI. This approach reduces acquisition time and allows easy comparison with echocardiography. This approach was considered for our study because the investigated animals developed a cardiomyopathy with global but not regional wall motion abnormalities. In the case of regional wall motion abnormalities, it is known that the Teichholz method is often not adequate. Our own studies have shown a very close agreement between volume calculations using the Teichholz method and a 3-dimensional analysis for mice, with LV volumes ranging from 50 to 100 µL (unpublished data). The methodologic approach was identical to that used in this study. However, we would expect an even better agreement between PET and MRI when using a 3-dimensional approach also for MRI, especially in a future application to ECG-gated PET with the aim of an exact PET-based volume measurement.
Many progressing heart diseases are associated with changes in both metabolism and heart size, and transgenic mouse models are successfully used to obtain insight into the genetic and molecular mechanisms involved. One of these is the mouse model used in this study: the cardiomyocyte-specific survivin knockout mouse (19,20). These mice progressively develop dilated cardiomyopathy and die prematurely. Specialized imaging techniques, among them PET and MRI, were specifically adapted to the imaging and data acquisition in these mice. In general, MRI has been used routinely to study geometric and functional properties of the heart such as LV volumes (22), whereas PET has generally been used previously to study biochemical properties (25). The comparison of groups between the cardiomyocyte-restricted survivin-deficient and wild-type mice presented in this study exemplifies this approach: The data indicative of LV dilation were obtained using either MRI or PET. However, the 2 imaging modalities may show a different power for detecting significant differences, which may also depend on the underlying disease process, making further studies desirable to elucidate these differences and their additive diagnostic potential. The value and significance of our study lie in its ability to allow the simultaneous analysis of LV metabolic properties such as glucose uptake together with a reliable estimation of LV size and, in diseased hearts, the extent of LV dilation.
| CONCLUSION |
|---|
|
|
|---|
State-of-the-art equipment for obtaining such results in murine hearts is currently available. Further development of this method should strive to include gated imaging and even correction for respiratory motion.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Lars Stegger, PhD, Department of Nuclear Medicine, University Hospital Münster, Albert-Schweitzer-Strasse 33, 48149 Münster, Germany.
E-mail: stegger{at}uni-muenster.de
| REFERENCES |
|---|
|
|
|---|
Related articles in JNM:
This article has been cited by other articles:
![]() |
C. M. Prado, E. J. Fine, W. Koba, D. Zhao, M. A. Rossi, H. B. Tanowitz, and L. A. Jelicks Micro-Positron Emission Tomography in the Evaluation of Trypanosoma cruzi-Induced Heart Disease: Comparison with Other Modalities Am J Trop Med Hyg, November 1, 2009; 81(5): 900 - 905. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Stegger, E. Heijman, K. P. Schafers, K. Nicolay, M. A. Schafers, and G. J. Strijkers Quantification of Left Ventricular Volumes and Ejection Fraction in Mice Using PET, Compared with MRI J. Nucl. Med., January 1, 2009; 50(1): 132 - 138. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Stegger, A.-N. Hoffmeier, K. P. Schafers, S. Hermann, O. Schober, M. A. Schafers, and G. Theilmeier Accurate Noninvasive Measurement of Infarct Size in Mice with High-Resolution PET J. Nucl. Med., November 1, 2006; 47(11): 1837 - 1844. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | RSS | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |