|
|
||||||||
HIGHLIGHTED TOPIC
Physiological Imaging of the Lung
1Department of Radiology and 2Division of Physiology, Department of Medicine, University of California, San Diego, La Jolla, California
Submitted 5 May 2006 ; accepted in final form 26 January 2007
| ABSTRACT |
|---|
|
|
|---|
0.1/decade until age 50–59 yr, and there was a significant positive relationship between RD and age (R = 0.48, P < 0.0005) and height (R = 0.39, P < 0.01), but not body mass index (R = 0.07, P = 0.67). Age and height combined in a multiple regression were significantly related to RD (R = 0.66, P < 0.0001). There was no significant relationship between RD and spirometry or arterial oxygen saturation. Ds was not related to age, height, spirometry, or arterial oxygen saturation. The lack of relationship between age and Ds argues against an intrinsic alteration in the pulmonary vascular branching with age as being responsible for the observed increase in global spatial perfusion heterogeneity measured by the RD. arterial spin labeling; pulmonary function; fractal dimension; relative dispersion
Arterial spin labeling (ASL) magnetic resonance imaging (MRI) techniques have been used to quantify regional perfusion in many organ systems (16, 35, 38, 40, 46). We have modified the general kinetic model originally developed for analyzing the ASL signal in the brain (46) to collect data in the lung within a single breath hold (3). This technique allows quantification of regional pulmonary perfusion using a specific MRI technique known as ASL-FAIRER (FAIRER stands for flow-sensitive alternating inversion recovery with an extra radio-frequency pulse, described below). While functional MRI techniques are more commonly used in the evaluation of perfusion in the brain (7), we have used ASL-FAIRER in the lung to identify changes in the global changes in spatial perfusion heterogeneity in response to physiological stimuli, such as acute hypoxia (30) and steep head-down tilt (28).
In this study, we used ASL-FAIRER to measure spatial perfusion heterogeneity in healthy, nonsmoking subjects with normal pulmonary function who were of different ages. We hypothesized that spatial perfusion heterogeneity would increase with increasing age, consistent with the changes in ventilation and gas exchange observed with normal aging (8, 17).
| METHODS |
|---|
|
|
|---|
MRI of Pulmonary Perfusion
Overview of ASL measurement of pulmonary perfusion. ASL inverts the proton magnetization of blood by applying a radio-frequency pulse, allowing these magnetically tagged protons in blood to act as an endogenous tracer for the evaluation of blood flow. During each measurement, two images are acquired of each section with the signal of blood prepared in a different way. Then the two images are subtracted, canceling the stationary signal within each image voxel and leaving just the signal due to blood delivered to that voxel during a specified time interval (3).
The ASL-FAIRER sequence used in this study was modified in house from one previously reported (33, 35) to allow for a single subtraction (3) and acquisition of data within a single breath hold. As stated, two images are acquired. The sequence is cardiac gated such that both images are acquired during diastole (while flow is minimal), and one complete systole interval occurs during the time between the labeling of the blood and the measurement of the delivery of the tagged blood to the imaged section. In the first image, termed the "tag" image, the magnetization of the blood is inverted at the beginning of the experiment with a preparatory inversion (180°) pulse applied to the whole lung (a spatially nonselective inversion). In the second "control" image, a preparatory inversion (180°) pulse is applied to only the section being imaged (a spatially selective inversion). For each image, a spatially selective 90° pulse is applied to the imaged section immediately after the preparatory inversion pulse. The effect of these radio-frequency pulses during the preparation phase is that the static magnetization within each voxel is reduced to zero in both experiments. The blood magnetization is fully inverted before the tag image, but it is fully relaxed before the control image. After a delay (TI), encompassing one R-R interval, the image is acquired. During this delay, blood flows into each voxel of the imaged section, and there is also relaxation of the magnetization. The static magnetization relaxes identically in both experiments, so when the two images are subtracted this signal is cancelled. In the tag image, the inverted magnetization of blood has relaxed part way to equilibrium, whereas for the control image the magnetization of blood remains fully relaxed. The difference signal (control – tag) measured for each voxel then reflects the amount of blood delivered, weighted with a decay factor due to the relaxation of the blood magnetization during the interval TI. Note that with TI chosen to approximate the R-R interval, the tagged blood signal is near its null point (zero signal) when the tag image is acquired, so a characteristic feature of these images is that vessels tend to be bright on the control image and dark on the tag image. In the ASL difference image (control – tag), the signal intensity for each voxel is proportional to the amount of blood delivered during one cardiac cycle (TI), and so is proportional to the local pulmonary blood flow expressed as volume of blood per volume of lung per minute.
The imaging scheme used for the ASL-FAIRER sequence was a half-Fourier acquisition single-shot turbo spin-echo (HASTE) acquisition (33). The HASTE imaging acquisition reduces susceptibility effects caused by inhomogeneous magnetic fields (10), maximizes the intrinsically low lung signal by allowing recovery of T2 signal loss (23), and reduces artifacts caused by respiratory and cardiac motion (27). The half-Fourier acquisition provides a twofold decrease in imaging time, allowing subsecond image acquisition. The application of this sequence to measuring pulmonary perfusion has been previously described in detail (3, 28, 30).
Data collection.
Subjects were imaged in the supine position using ASL-FAIRER along with a HASTE technique as described above. A phased-array torso coil was used. Imaging was performed at functional residual capacity, improving the signal-to-noise characteristics of the image by increasing the amount of tissue relative to air within the lungs (34). Following a minor amount of training, subjects were able to consistently breath hold at this level, as documented by the consistent position of the diaphragm on the nonsubtracted images. The imaging plane was a single, 15-mm-thick, coronal slice, chosen from an initial series of scout images, located within the posterior one-half of the lung in line with the posterior descending aorta. The position was chosen to avoid the heart and large central vessels. A TI of 700–800 ms was used. This was
80% of the R-R interval and was kept constant for each subject. Additional magnetic resonance parameters were the following: echo time = 36 ms, flip angle = 30°, number of excitations = 1, image matrix = 128 x 256, and field of view = 400 mm. Thus each voxel had a fixed dimension of
1.5 x 3 x 15 mm (
70 mm3). All sequence parameters were kept within the established Food and Drug Administration safety guidelines for routine clinical magnetic resonance examinations. Three to five perfusion-weighted image maps were obtained for each subject, and the resultant data were averaged.
Quantification of perfusion heterogeneity. MRI data were analyzed using in-house software developed in the MATLAB programming environment (The Mathworks, Natick, MA). For this study, image data from the right lung was used as a single, large, region of interest. The contour of the right lung was manually outlined on a nonsubtracted image, and this was used to define an image mask for the subtracted image. The mean signal intensity and SD were measured for all voxels within this region of interest, and the relative dispersion was calculated (19). The relative dispersion (also known as the coefficient of variation) of blood flow is defined as the SD of signal intensity divided by the mean signal. This provides a global index of spatial heterogeneity without considering the specific anatomic location of the heterogeneous flow. A second measure of perfusion heterogeneity, the fractal dimension, was calculated by calculating the relative dispersion of the image in progressively larger blocks centered on each voxel (19, 21). For blocks that fall partially outside the mask, the weights are adjusted to exclude the effect of the masked-out portions. The slope of the relationship, representing the fractal dimension, is obtained from the plot of log relative dispersion vs. log number of voxels averaged (19, 21). The fractal dimension is therefore a measure of the scale-independent perfusion heterogeneity (19) and indicative of pulmonary vascular branching structure (20). A fractal dimension of 1 represents a uniform distribution, whereas a fractal dimension of 1.5 represents a random distribution.
Data analysis. Data were first analyzed using analysis of covariance (Statview 4.1, SAS Institute, Cary, NC) to determine the effect of any potential gender differences, with age as a covariate for perfusion heterogeneity. Subjects were also grouped by 10-yr intervals (age 20–29, 30–39, 40–40, 50–59, and 60+ yr), and the data were analyzed using ANOVA (Statview 4.1, SAS Institute, Cary, NC) with one grouping variable (age, 5 levels). Where overall significance was established, post hoc testing was applied to determine where the differences occurred. Linear regression was used to relate the change in blood flow heterogeneity to age and physiological variables. Significance was accepted at P < 0.05 (2 tailed). All data shown are means (SD).
| RESULTS |
|---|
|
|
|---|
|
|
|
Effect of age, height, body mass index, pulmonary function, and arterial oxygen saturation.
The mean perfusion heterogeneity, as measured using relative dispersion, was 1.00 (SD 0.18) averaged over all subjects. It was lowest in the subjects aged 20–29 yr [0.85 (SD 0.10)] and increased significantly with increasing age (Fig. 3; P < 0.005) by
0.1 per decade of age until age 50–59 yr. Using linear regression, there was a significant positive relationship between relative dispersion and age (Fig. 4; R = 0.48, P < 0.0005). Interestingly, there was also a significant positive relationship between relative dispersion and height (Fig. 5; R = 0.39, P < 0.01), and when age and height were combined in a multiple regression these two variables were highly significant (R = 0.66, P < 0.0001). The regression equation is the following: relative dispersion = 0.007·age (yr) + 0.009·height (cm) –0.809. There was no significant relationship between relative dispersion and body mass index (R = 0. 07, P = 0.67). The mean value for fractal dimension was 1.18 (SD 0.03) averaged over all subjects, and like the value for relative dispersion it was not different between men and women (P = 0.78). Unlike relative dispersion, fractal dimension did not change with age across 10-yr age groups (P = 0.3), and there was no significant relationship between relative dispersion and age (Fig. 4; R = 0.05, P = 0.70). There was also no significant relationship between fractal dimension and height, (R =0.11 , P = 0.45 ) or body mass index (R = 0.22 , P = 0.2).
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Aging and changes in pulmonary perfusion and ventilation-perfusion inequality. Very little is known about the spatial changes in pulmonary blood flow that accompany normal aging. Holland and colleagues (29) evaluated the regional distribution of pulmonary blood flow in six male subjects between the ages of 65 and 75 yr using an intravenous injection of 133Xe dissolved in saline. They found a relative increase in blood flow to the upper lungs while the subjects were upright, compared with similar measures made in younger men. However, pulmonary arterial pressure increases with increasing age (36); thus these findings may be explained partly on the basis of increased pulmonary vascular pressure leading to increased perfusion of zone 1 in the lung (2, 45). In addition, four of the six subjects had significant smoking histories, averaging 35 pack-years, and thus the significance of these findings is uncertain.
More data are available describing the changes in ventilation-perfusion inequality associated with aging. Using various imaging techniques, ventilation-perfusion inequality has been variably reported to be increased in older subjects compared with younger ones (29, 37) or to be unchanged (25). Using the multiple inert gas technique, considered to be the gold standard for measuring ventilation-perfusion inequality, Cardus et al. (8) found an increase in ventilation-perfusion inequality of
5–6% per decade in 64 normal subjects aged 18–71 yr, similar to the decreases found in spirometric data with increasing age (13). This increase was not related to lung volumes measured with spirometry or to other physiological variables such as height or weight of the subjects.
In our study, we found a similar increase in perfusion heterogeneity as measured by the relative dispersion, which is a global index of the degree of spatial heterogeneity that provides an overall index of spatial heterogeneity without considering the specific anatomic location of the heterogeneous flow. The increase in relative dispersion with increasing age of
5% per decade of age is similar to the magnitude of the change in ventilation-perfusion heterogeneity with age (8). Also similar to the study of Cardus et al. (8), the increase in perfusion heterogeneity we observed with age was not related to lung volumes measured with spirometry. One possible explanation for this may be that spirometry measures variables related to pulmonary ventilation rather than perfusion per se. Surprisingly, we found a positive relationship between pulmonary blood flow heterogeneity and height, such that taller individuals had greater pulmonary blood flow heterogeneity, a result not found by Cardus et al. Although it is possible that this may be explained partly on the basis of lung size, this would appear unlikely because of a lack of relationship between spirometric data and perfusion heterogeneity. However, although an individual with larger lungs by spirometry would be expected to have a larger pulmonary blood volume, the scaling of the two variables with body size may not exactly correspond, and thus it may be possible that perfusion heterogeneity would correlate with measures of pulmonary blood volume or diffusing capacity, which were not measured in the present study.
In simple terms, a fractal structure is one that is composed of pieces that are similar to the structure as a whole. The pattern of blood flow in the lung has been described as fractal in nature, and the fractal dimension of pulmonary blood flow has been reported for several mammals, dogs, sheep, and horses, (9, 21, 41), ranging from 1.09 to 1.17. Our value of 1.18 (SD 0.03) for healthy normal humans is therefore very similar. Interestingly, we did not find a relationship between fractal dimension and age. The implication of this finding is that the overall anatomic vascular structure in the lung is retained with increasing age in the healthy human lung, despite the global increase in spatial heterogeneity of pulmonary blood flow.
The reason for the increase in perfusion heterogeneity with increasing age is uncertain. One possibility is that this represents an intrinsic alteration in the pulmonary vasculature with aging because diffusing capacity for carbon monoxide (12), pulmonary capillary blood volume (22), and capillary density are known to decrease with age (6). However, the finding of the unchanged fractal dimension with age argues against the possibility of an alteration in branching structure, at least at the level of resolution of our technique. Alternately, because closing volume increases with aging (5), this may represent local hypoxic pulmonary vasoconstriction in response to regional alveolar hypoxia. However, 100% oxygen breathing does not alter areas of low ventilation-perfusion ratio present in the lung bases in elderly subjects (25), arguing against this idea. The lack of a change in fractal dimension does not rule this possibility out because the level of resolution of our measurements is not at a level to image blood flow in very small arterioles. As mentioned earlier, pulmonary arterial pressure increases with increasing age (36). However, our findings cannot represent an effect of increased pulmonary arterial pressure causing increase perfusion to the upper lung zones for two reasons. First, the effect of increased pressure would be expected to make perfusion more uniform rather than less uniform. Second, the imaging in our study was a coronal section of the lung in a supine subject, and thus the data from the image are all in the same gravitational plane, so differences observed between subjects with age cannot reflect gravitationally based changes.
Another possibility that should be considered is that any increases in the heterogeneity of lung density with age would likely be manifest in our study as increases in perfusion heterogeneity, because our technique measures perfusion per unit volume (voxel) rather than per gram of tissue. For example, in a lung in which perfusion per gram of lung tissue was uniform, regions of high or low tissue density would appear as regions of high and low perfusion, respectively, merely because the number of capillaries within a voxel is a function any compression of that voxel. Although data from studies using both computed tomography scanning and positron emission tomography do not support any significant age-related effect on overall lung density (4, 18), the effects on heterogeneity of lung density are not well described. In addition, there is evidence that alveolar size increases with age and that this increase is not uniform across the lung (15); thus alterations lung density heterogeneity may well play a role. Because we do not have regional lung density measurements in our subjects we cannot address the potential contribution of this mechanism to our observations.
Measurement of Perfusion Heterogeneity With ASL-FAIRER
The techniques used to evaluate pulmonary blood flow heterogeneity are an important feature of this study and have been recently described (3). MRI using ASL techniques has been widely used to determine regional blood flow in other organ systems, such as brain (7). The technique has been validated in tube-flow models (1), heart (38), brain (44), and skeletal muscle (39). There are a number of potential methods to evaluate the spatial heterogeneity of blood flow (reviewed in Ref. 19). Each coronal image generates roughly 15,000 voxels of data using this technique. A significant issue, then, is how to evaluate and express this large data set to compare among subjects. The simplest measurement of heterogeneity is relative dispersion, the SD of the image signal divided by the mean. The advantages of the relative dispersion are that it is robust and reliable (
4% variation over repeated measures as seen in Fig. 1), it is independent of changes in absolute signal intensity that might occur between subjects, and it is largely unaffected by small changes in body position or depth of inspiration. Also, because of these factors, it has the advantage that data can be compared between different imaging centers and allows the establishment of normative values for healthy subjects. The obvious disadvantage of the relative dispersion is the loss of specific anatomic information. Nonetheless, the ability to quantitatively characterize one subject as having a different heterogeneity than another has been shown to be a useful tool in previous studies using this technique (28, 30). In addition, we reported an additional measure of perfusion heterogeneity, the fractal dimension, which was unchanged with increasing age. Like relative dispersion, it is also a highly reliable measure. The nature of data in the present study differs from that previously reported in animals, since they are derived from a single coronal slice of lung and not from a three-dimensional data set. This means that lung voxels in a three-dimensional data set that would normally contribute to the calculation of the fractal dimension at a particular distance contribute to the calculation of the two-dimensional data set only if they fall within the imaging plane. However, fractal analysis of data from dog lung suggests that the heterogeneity of a single slice gives similar results to the analysis of the lung as a whole (43).
Limitations
There are some limitations to the study that should be considered. We imaged a coronal slice of the right lung that encompassed a sample of the lung field, and therefore changes in pulmonary blood flow heterogeneity we observed with increasing age could have been more or less pronounced in regions of the lung that were not assessed. Also, ASL-FAIRER provides an image map of all tagged protons move into the imaging slice during the delay between tagging and image acquisition, and thus components of both pulmonary arterial and venous blood flow are likely present, the significance of which is unclear (3). Nevertheless, this limitation applies to all of our subjects irrespective of age and thus should not affect our conclusions. Finally, we did not perform absolute quantification of pulmonary perfusion. However, because in healthy subjects the lungs receive virtually the entire cardiac output, which is readily measured using a variety of well-established techniques, the relative distribution of perfusion in the lung is the much more interesting data.
In conclusion, understanding the physiological changes that occur with healthy aging is important in the detection of early pathology, in the evaluation of disease progression, and in assessing the effects of therapy. In particular, with the development of new diagnostic techniques such as functional MRI of the lung, it is important to document any age-related changes in healthy subjects, because many patients with diseases such as chronic obstructive lung disease are elderly. If these techniques are to be used to probe pulmonary function in disease, the effects of disease need to be distinguished from the effects of normal aging. We have demonstrated an increase in a global measure of spatial pulmonary blood flow heterogeneity, the relative dispersion with aging among healthy, nonsmoking subjects using ASL-FAIRER MRI and established normative data for this measure. The possibility of age related changes in lung density heterogeneity as an explanation for our findings cannot be ruled out. However, the lack of a change in the fractal dimension with aging suggests that the overall pulmonary vascular structure is retained in normal aging.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. J. Arai, A. C. Henderson, D. J. Dubowitz, D. L. Levin, P. J. Friedman, R. B. Buxton, G. K. Prisk, and S. R. Hopkins Hypoxic pulmonary vasoconstriction does not contribute to pulmonary blood flow heterogeneity in normoxia in normal supine humans J Appl Physiol, April 1, 2009; 106(4): 1057 - 1064. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. R. Hopkins, A. C. Henderson, D. L. Levin, K. Yamada, T. Arai, R. B. Buxton, and G. K. Prisk Vertical gradients in regional lung density and perfusion in the supine human lung: the Slinky effect J Appl Physiol, July 1, 2007; 103(1): 240 - 248. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |