|
|
||||||||
Technical Developments |
1 From the Department of Radiology, University of Pennsylvania Medical Center, 1 Silverstein, 3400 Spruce St, Philadelphia, PA 19104-4283. Received March 5, 1998; revision requested June 5; final revision received October 27; accepted March 16, 1999. Supported in part by a seed grant from the Society for Thoracic Radiology. Address reprint requests to D.A.R. (e-mail: roberts@oasis.rad.upenn.edu).
| Abstract |
|---|
|
|
|---|
Index terms: Lung, MR, 60.121412, 60.121417 Lung, perfusion, 60.91 Lung, ventilation, 60.91
| Introduction |
|---|
|
|
|---|
Magnetic resonance (MR) imaging of the lung parenchyma is intrinsically difficult due to the large degree of microscopic magnetic field inhomogeneity, which causes the observed transverse relaxation time (T2*) to be extremely short (6,7). However, the recent development of gradient-echo pulse sequences with extremely short echo time and projection reconstruction techniques has allowed investigators to image the pulmonary parenchyma (810). Recently, several MR methods of pulmonary perfusion imaging have been described (1116). To date, the MR methods described have been breath-hold techniques that are based on dynamic observation of the pulmonary magnetization after either injection of a bolus of paramagnetic contrast agent (13,15,16) or spatially selective inversion of the mediastinum (12,14). An alternative noninvasive technique exploits differences in signal intensity (SI) between different cardiac phases (11). Although promising, these methods necessitate breath-holding maneuvers, which are often impossible for patients with pulmonary disease. In addition, injection of a contrast agent for the perfusion study introduces additional risk and expense. The quality of subsequent MR angiograms is also reduced by contrast material injection owing to reduced background suppression, which causes a reduction in the contrast-to-noise ratio.
The arterial spin-tagging method of perfusion imaging is based on the steady-state magnetic labeling of the arterial water supplying an organ (17,18). This results in a measurable change in the steady-state magnetization from which values for tissue-specific perfusion can be calculated. This method has been applied successfully in both animal and human studies (1723). To our knowledge, the application of this method to the study of pulmonary perfusion has not been reported. In this study, we evaluated use of a three-dimensional, respiratory-triggered implementation of this technique to depict regional pulmonary perfusion in healthy subjects without the need for breath holding.
| Materials and Methods |
|---|
|
|
|---|
Image Acquisition
MR imaging data were acquired with a 1.5-T imaging system equipped with high-power gradients (Echospeed Signa, revision 5.6 software; GE Medical Systems, Milwaukee, Wis). The gradients had a peak amplitude of 2.2 G/cm. For monitoring purposes, cardiac electrodes were placed on the chest wall, and a respiratory excursion transducer was looped around the chest. A quadrature birdcage coil encompassing the entire body was used for signal transmission. To maximize the signal-to-noise ratio (SNR), a phased-array coil was placed around the thorax for signal reception that consisted of anterior and posterior two-element arrays. An automated gradient shim was performed, and the radio-frequency (RF) synthesizer frequency was centered on the proton water resonance frequency.
The pulse sequence used to acquire data appears in Figure 1. It is a three-dimensional, compact, gradient-echo pulse sequence that incorporates an asymmetric, truncated sinc RF slab excitation pulse of 1.28-msec duration. The following imaging parameters were used: repetition time msec/echo time msec of 50/716760 with 30° flip angle, 710-mm section thickness, 4048-cm field of view, and 256 x 64 x 16 matrix. The nominal in-plane resolution was approximately 2 x 8 mm. To minimize the echo time, a fractional echo was acquired with a 64-kHz receiver bandwidth. A sagittal volume was excited that included the right lung, as indicated in Figure 2. Trapezoidal gradient and RF pulses were placed in the interpulse interval to label blood flow according to the principle of flow-induced adiabatic fast passage (24,25). The offset frequency of the labeling pulse was adjusted in each case so that the location of the plane of arterial labeling intersected the right pulmonary artery (Fig 2). In all subjects, the right pulmonary artery was readily identified on the scout localizer image by both the radiologist (D.A.R.) and technologist. The amplitude of the labeling RF pulse, Bi, was 24 mG, and its duration was 40 msec, corresponding to a duty cycle of 80%. The amplitude of the inversion gradient, Gi, was 0.1 G/cm. Under these conditions, the labeling RF pulse created a band of saturation approximately 0.6 mm thick, as indicated by the lines in Figure 2. On the basis of measurements obtained on our localizer images, the average diameter of the right pulmonary artery was 1.6 cm. Assuming that the right pulmonary artery receives half the normal cardiac output of 5 L/min, then the estimated average velocity of flow is 8.3 cm/sec. This yields an adiabaticity factor (21) for the inversion of 2.8, yielding an estimated degree of inversion of 95%. As the duty cycle is 80%, this yields an estimated net inversion efficiency of 76%. This likely represents an overestimation of the degree of inversion as higher velocities occur during the cardiac cycle that yield lower inversion efficiencies. Although one would ideally like to achieve 100% inversion efficiency to maximize the observed SI change, this was the highest achievable value to remain within limitations of the specific absorption rate.
|
|
|
This method, as described, does not provide an absolute quantitation of perfusion. Instead, images are produced of the relative regional SI change induced by means of arterial labeling. For this reason, accurate knowledge of the degree of arterial inversion at the level of the pulmonary artery is not important for the detection of spatial defects in a given individual. However, this does become important when attempting to compare results obtained in the same individual at different points in time.
Exercise Studies
Subjects (n = 8) were instructed to exercise by performing a leg-raising maneuver until the point of exhaustion. The basal heart rate and the heart rate immediately after exercise were recorded. Perfusion data were acquired before and after exercise. Comparison of the heart rate and the pre- and postexercise data was performed by means of the paired Student t test. In one patient, images obtained after exercise were corrupted by motion artifact and were excluded. Therefore, postexercise data were available in only seven patients.
Data Analysis
All data analysis was performed off-line on a workstation (Silicon Graphics, Mountain View, Calif) by using the Interactive Data Language (Research Systems, Boulder, Colo). Calculation of the perfusion images typically required 10 minutes per patient, including both data transfer to the workstation and subsequent analysis. The total time necessary to produce perfusion images of the right lung was 30 minutes.
Prior to formation of magnitude image subtraction data, the inversion image data were corrected for off-resonance effects of the labeling RF. A region of interest was drawn interactively around the lung in each image in the control and magnetization transfer control data sets. Pixels composing pulmonary parenchyma were then identified by means of segmentation of the magnitude image data within these regions of interest by using an upper threshold of 150. This value was chosen empirically, as segmentation at this level eliminated contributions from the pulmonary arteries in all patients. Note that this form of segmentation does not eliminate contributions from the bronchi, which are expected to have zero SI as they contain air. This represents a small error in the analysis. The degree of magnetization transfer saturation, S, caused by the labeling pulse in pulmonary parenchyma was computed as the ratio of the average SI with the pulse divided by the average SI without the pulse. This was computed for each section in the magnetization transfer control acquisition as a function of distance, z', from the plane of the lateral labeling pulse. It is important to note that the physical location of the plane used to determine the degree of saturation was different from the corresponding plane in the inversion data set. Therefore, a point-by-point correction could not be performed. The computation yielded a correction factor, S'(z'), for magnetization transfer effect as a function of z'. Assuming that the lung is spatially homogeneous in its off-resonance properties, the appropriate correction factor for arterial spin-tagging images is S(z) = S'(z'), where z is the distance of an arterial spin-tagging image from the plane of arterial labeling. In some cases, the extreme medial or lateral images in the volume contained so little lung parenchyma that an accurate calculation could not be performed. In these cases, no correction was performed (S[z] = 1.0).
Images of the regional percentage SI change after correction for off-resonance effects (
CORR) induced upon arterial labeling in a section located a distance, z, from the plane of arterial labeling were calculated as
CORR was set to zero. The SI in the control image was determined by computing the average value in a small region of interest placed around the peripheral aspect of the lung in the central section. The noise was computed as the SD in a small region of interest outside the body, in the upper left corner of the same section. This location was chosen to minimize the effect of potential respiratory motion artifacts on the noise calculation. The SNR was then computed as their ratio.
Values for the mean SI change in pulmonary parenchyma were computed as follows: (a) a region of interest was interactively placed around the lung in each of the control images, and (b) segmentation within this region was then performed with use of an upper threshold of 150 to identify pixels composing lung parenchyma. Inspection of the segmented images showed that this eliminated contributions from the central pulmonary arteries.
| Results |
|---|
|
|
|---|
|
|
We also obtained preliminary data in a patient with chronic obstructive pulmonary disease (Fig 5). The reformatted sagittal computed tomographic (CT) images reveal bullous disease, which asymmetrically affects the right upper and middle lobes, with relative sparing of the right lower lobe (Fig 5, top row). Despite the presence of bullous disease, an SNR of approximately 11:1 was obtained in the affected parts of the lung on the control images (Fig 5, middle row). The arterial spin-tagging perfusion images demonstrated marked hypoperfusion in a distribution that closely paralleled the distribution of chronic obstructive pulmonary disease (Fig 5, bottom row). Hypoperfusion was observed to the extreme lateral aspect of the lung, where off-resonance effects were quite small.
|
| Discussion |
|---|
|
|
|---|
However, the presence of thromboembolic disease leads to decreases in pulmonary perfusion and blood volume while at the same time prolonging the transit time. All of these changes are expected to lead to a marked reduction in the observed SI change upon arterial labeling. Therefore, the arterial spin-tagging method should be especially sensitive to the presence of pulmonary thromboembolic disease. This may have important clinical implications as recent data suggest that as many as 30% of patients with pulmonary emboli have involvement at the subsegmental level or below (32). This area has been difficult to evaluate with cross-sectional imaging techniques (33,34).
The assumption of a spatially uniform off-resonance characteristic in the lung represents a potential source of error. It is indeed probable that in the presence of disease there are marked regional variations in the magnetization transfer properties of the lung that would invalidate this assumption. The use of amplitude-modulated control irradiation would eliminate this as a potential source of error; however, such pulses could not be implemented at present owing to limitations of the RF amplifier.
The spin-tagging pulse labels all blood flow occurring at the level of the pulmonary hilum. Therefore, the perfusion images reflect contributions from both the pulmonary and bronchial arteries. This is an important difference from radionuclide ventilation-perfusion scanning in which, in the absence of shunt formation, uptake reflects only pulmonary arterial flow. Results obtained with the arterial spin-tagging technique may differ from those obtained on ventilation-perfusion scans, especially in patients with chronic lung disease in whom prominent bronchial collateral vessels may exist. In addition, this may decrease the sensitivity of this method for the detection of pulmonary embolism as bronchial supply to the affected segment may persist.
The arterial spin-tagging method may be used to obtain noninvasive images of pulmonary blood flow without the need for patient breath holding. The arterial spin-tagging technique is entirely noninvasive and eliminates the risks and expense associated with contrast material administration. In comparison with previously described MR techniques for pulmonary perfusion imaging (11,12,14), this technique does not require breath holding, a task that is often impossible for patients with pulmonary disease. In fact, this method is ideally suited to the evaluation of disease in patients receiving artificial ventilation as nearly perfect respiratory triggering may be achieved in this setting. As a steady-state technique, it theoretically allows a greater effect of arterial labeling to build up during the acquisition. Alternative MR pulmonary perfusion imaging techniques that require the intravenous injection of paramagnetic contrast agents have also been described (13,16). These methods suffer from the need for breath holding, a task that is often impossible for patients with clinically important pulmonary disease. Although pulmonary MR angiography has been shown to be a highly accurate and safe method for the diagnosis of pulmonary embolism, the need to inject gadolinium for the perfusion study inevitably reduces the quality of subsequent MR angiograms owing to decreased background suppression. In addition, the accuracy of pulmonary MR angiography at the subsegmental level and below is not known (34). There is evidence that a major fraction of pulmonary emboli occur at the subsegmental level or below (32). As the spin-tagging method does not involve contrast material injection, this technique may become useful as a means of focusing the MR angiographic examination in a specific region of the lung, thereby potentially increasing the spatial resolution and diagnostic yield of pulmonary MR angiograms and allowing the diagnosis of subsegmental pulmonary embolism. In addition, in conjunction with recently developed noninvasive methods of MR ventilation imaging (35,36), it has the potential to constitute a component of an entirely noninvasive, high-spatial-resolution form of MR ventilation-perfusion imaging.
The clinical utility of the arterial spin-tagging method, although promising, may be limited by several technical factors. Specifically, its success rate and accuracy in a patient population are unknown. Techniques must be developed to allow perfusion imaging of the left lung, where effects due to cardiac motion and intracardiac labeling of blood flow must be addressed. Such techniques may lead to increases in acquisition time. The arterial spin-tagging method in isolation is likely to suffer from indeterminate results in the clinical setting of pulmonary embolism for the same reasons that radionuclide scanning does. Despite these limitations, it is reasonable to conjecture that this method, coupled with high-spatial-resolution MR angiography and ventilation imaging, may form an important part of an integrated functional and anatomic MR examination of the lung.
| Acknowledgments |
|---|
| Footnotes |
|---|
Author contributions: Guarantor of integrity of entire study, D.A.R.; study concepts, D.A.R., J.S.L.; study design, D.A.R., M.D.S., R.E.L.; definition of intellectual content, M.D.S., D.A.R., J.S.L.; literature research, D.A.R.; clinical studies, D.A.R., R.E.L., L.D., M.D.S.; data acquisition, J.A.H., D.A.R., R.R.R.; data analysis, D.A.R., R.R.R.; statistical analysis, D.A.R.; manuscript preparation and editing, D.A.R.; manuscript review, R.R.R., W.B.G., R.E.L., M.D.S.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. C. Miller, H. H. Pien, D. Sahani, A. G. Sorensen, and J. H. Thrall Imaging Angiogenesis: Applications and Potential for Drug Development J Natl Cancer Inst, February 2, 2005; 97(3): 172 - 187. [Abstract] [Full Text] [PDF] |
||||
![]() |
U. J. Schoepf and P. Costello CT Angiography for Diagnosis of Pulmonary Embolism: State of the Art Radiology, February 1, 2004; 230(2): 329 - 337. [Abstract] [Full Text] [PDF] |
||||
![]() |
U. J. Schoepf, R. Bruening, H. Konschitzky, C. R. Becker, A. Knez, J. Weber, O. Muehling, P. Herzog, A. Huber, R. Haberl, et al. Pulmonary Embolism: Comprehensive Diagnosis by Using Electron-Beam CT for Detection of Emboli and Assessment of Pulmonary Blood Flow Radiology, December 1, 2000; 217(3): 693 - 700. [Abstract] [Full Text] |
||||
![]() |
D. D. Maki, W. B. Gefter, and A. Alavi Recent Advances in Pulmonary Imaging Chest, November 1, 1999; 116(5): 1388 - 1402. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |