Journal of Applied Physiology Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 93: 1841-1851, 2002; doi:10.1152/japplphysiol.00223.2002
8750-7587/02 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (41)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Musch, G.
Right arrow Articles by Venegas, J. G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Musch, G.
Right arrow Articles by Venegas, J. G.
Vol. 93, Issue 5, 1841-1851, November 2002

Topographical distribution of pulmonary perfusion and ventilation, assessed by PET in supine and prone humans

Guido Musch1, J. Dominick H. Layfield1,4, R. Scott Harris2, Marcos F. Vidal Melo1, Tilo Winkler1,5, Ronald J. Callahan3, Alan J. Fischman3, and Jose G. Venegas1

Departments of 1 Anesthesia and Critical Care, 2 Medicine (Pulmonary and Critical Care Unit), and 3 Radiology (Division of Nuclear Medicine), Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts 02114; 4 Massachusetts Institute of Technology, Boston, Massachusetts 02139; and 5 Clinic of Anesthesiology and Intensive Care Medicine, University Clinic Carl Gustav Carus, Dresden University of Technology, Dresden 01307, Germany


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

Using positron emission tomography (PET) and intravenously injected 13N2, we assessed the topographical distribution of pulmonary perfusion (Q) and ventilation (V) in six healthy, spontaneously breathing subjects in the supine and prone position. In this technique, the intrapulmonary distribution of 13N2, measured during a short apnea, is proportional to regional Q. After resumption of breathing, regional specific alveolar V (sVA, ventilation per unit of alveolar gas volume) can be calculated from the tracer washout rate. The PET scanner imaged 15 contiguous, 6-mm-thick, slices of lung. Vertical gradients of Q and sVA were computed by linear regression, and spatial heterogeneity was assessed from the squared coefficient of variation (CV2). Both CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP> and CV<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP> were corrected for the estimated contribution of random imaging noise. We found that 1) both Q and V had vertical gradients favoring dependent lung regions, 2) vertical gradients were similar in the supine and prone position and explained, on average, 24% of Q heterogeneity and 8% of V heterogeneity, 3) CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP> was similar in the supine and prone position, and 4) CV<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP> was lower in the prone position. We conclude that, in recumbent, spontaneously breathing humans, 1) vertical gradients favoring dependent lung regions explain a significant fraction of heterogeneity, especially of Q, and 2) although Q does not seem to be systematically more homogeneous in the prone position, differences in individual behaviors may make the prone position advantageous, in terms of V-to-Q matching, in selected subjects.

functional lung imaging; positron emission tomography; gas exchange; heterogeneity; prone position


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

RECENTLY THERE HAS BEEN A renewed interest in the effects of body position changes on the topographical distribution of pulmonary perfusion (Q) and ventilation (V) for two reasons. First, the prone position is increasingly used in the treatment of patients with acute respiratory distress syndrome (ARDS), and a recent clinical trial suggested that prone positioning may improve survival in a subset of patients with ARDS (6). Whether the prone position is associated with a more uniform distribution of Q, V, or both is unclear (16). Although animal studies have shown Q, V, and V/Q to be more uniform in the prone position (15, 29), the results in normal humans are discordant. Whereas some studies reported a more uniform distribution of Q (17, 18) and V (1) in the prone position, other studies reported no difference in Q (10, 11) and V (4, 11) gradients between the supine and prone positions or even greater Q gradients in the prone position (2). Second, animal studies have employed position changes to draw inferences on the role of gravity in determining the distribution of Q. Some of these studies suggested that gravity is a minor determinant of Q distribution in recumbent postures (7), whereas other studies found a much greater contribution of gravity to the distribution of Q (9, 29). The extent to which these results in animals can be extrapolated to humans is unclear.

Positron emission tomography (PET) imaging of the lung during a constant rate intravenous infusion of [13N]nitrogen (13N2) has been used to measure regional V/Q in humans (24). We modified that technique into a bolus infusion of 13N2 during a short apnea. By measuring the concentration of 13N2 during apnea and the ensuing period of 13N2 washout, the topographical distributions of Q and V can be assessed separately. The potential advantages of this method include its ability to assess both Q and V with a single administration of tracer, its minimal invasiveness, and the low radiation exposure to the subject. These characteristics might make this method suitable for future clinical applications in a variety of lung diseases.

In the present study, we used the 13N2 bolus infusion technique and PET to pursue two aims: 1) to assess the topographical distribution of Q and V in healthy, spontaneously breathing humans in the supine and prone position and 2) to quantify the amount of spatial heterogeneity of Q and V explained by a vertical gradient, after accounting for the estimated contribution of random imaging noise to the measured heterogeneity of Q and V.


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

Experimental Subjects

We studied six healthy adults (3 men, 3 women; ages 19-32 yr). All subjects were nonsmokers and nonobese and had normal pulmonary function tests. Informed consent was obtained from each subject before the study.

13N2-Saline Bolus Infusion Technique

To assess regional lung function, we used the 13N2-saline bolus infusion technique (9, 14, 29, 32). Simultaneously with the start of a 13N2-saline bolus infusion, the collection of a PET scan of the thorax was initiated. Each PET scan consisted of a series of consecutive PET images, acquired during an apnea (40 s) and the ensuing period of tracer washout. As the injected tracer reached the lung, its concentration rose until it reached a plateau for the remainder of the apnea (~30 s). When the subject resumed breathing, the tracer concentration decreased as 13N2 was eliminated by ventilation (washout). From the tracer kinetics (Fig. 1), regional lung function was assessed for each volume element of lung in the tomogram (i.e., voxel) as follows.


View larger version (16K):
[in this window]
[in a new window]
 
Fig. 1.   Tracer kinetics plot showing the time course of the mean concentration of 13N2 in the imaged lung field. Within ~10 s from the start of the bolus infusion, the tracer concentration reached a plateau that was maintained through the remainder of the apnea. As the subject resumed breathing, 13N2 was eliminated by ventilation (washout).

Assessment of Q. The theoretical basis for the measurement of Q from the 13N2 concentration during apnea was provided in detail by Mijailovich et al. (14). Because of the low solubility of 13N2 in blood and tissues (partition coefficient lambda water/air = 0.015 at 37°C), on arrival into the pulmonary capillaries, virtually all the tracer diffuses into the alveolar air space at first pass (23, 24, 30), and regional tracer content during apnea is proportional to regional Q (14, 32). Mean-normalized regional Q could thus be expressed as the ratio between tracer concentration in each voxel and mean tracer concentration of all voxels in the imaged lung field, measured during the plateau phase of the apnea. This approach to the calculation of Q is justified by the fact that 13N2 reabsorption into pulmonary venous blood during the apnea is negligible in a normally aerated and perfused lung (see APPENDIX and Fig. 1).

Assessment of V. As the subject resumed spontaneous breathing, 13N2 was eliminated from the lung almost exclusively by ventilation (30). Specific alveolar ventilation [sVA, alveolar ventilation (VA) per unit of alveolar gas volume] was calculated as the reciprocal of the time constant (tau ) of the tracer washout curve (see APPENDIX). For a monoexponential washout, tau  corresponds to the mean residence time (MRT) of the tracer (31) defined as
MRT<IT>=</IT><LIM><OP>∫</OP><LL>0</LL><UL><IT>∞</IT></UL></LIM> <FR><NU>c(<IT>t</IT>)</NU><DE>C<SUB>w</SUB></DE></FR> d<IT>t</IT> (1)
where c(t) is the concentration of 13N2 at time t and Cw is the concentration of 13N2 at the beginning of the washout [i.e., c(0) = Cw]. Because the amount of tracer remaining in the lung at the end of the washout was negligible (0.2-2.1% of the initial concentration), the integral in Eq. 1 was approximated by the integral of the tracer concentration over the finite washout imaging time. For each voxel, sVA was calculated as the reciprocal of the MRT.

Assessment of Q/V. From the derivation above, it follows that the integral of the tracer concentration during the washout is proportional to the Q-to-sVA ratio (Q/sVA; see APPENDIX). This can be intuitively appreciated from Eq. 1, which shows that the integral of the tracer washout curve is the product of Cw, which is proportional to Q, and the MRT.

As shown in the APPENDIX, neglecting tracer reabsorption by the pulmonary venous blood during the washout results in a theoretical error in the estimation of sVA and Q/sVA in the order of 1.5% in a normally aerated and perfused lung.

PET Transmission Scan

A 10-min transmission scan was recorded for each subject in both body positions by using a uniform rotating pin-source of 68Ge (Fig. 2). The transmission scan was used to correct the emission scan for energy attenuation caused by body tissues and supporting structures, to demarcate the lung field, and to calculate fractional gas content (Fgas), as previously described (4, 9, 30). Fgas represents the fraction of the volume of a voxel occupied by gas.


View larger version (76K):
[in this window]
[in a new window]
 
Fig. 2.   Positron emission tomography (PET) transmission, perfusion (Q), specific alveolar ventilation (sVA), and perfusion-to-ventilation ratio (Q/sVA) images. In each image, slices were ordered from the most cranial (upper left) to the most caudal (lower right). In each slice, the left lung is on the left and the right lung on the right. For Q, sVA, and Q/sVA, 2 color-coded scales are shown: a relative scale (right), in which the measurements are normalized by their mean value, and an absolute scale (left), with the original unit of measure (u.m.). Whereas the u.m. of sVA corresponds to the physical u.m. of specific ventilation (i.e., s-1), our measurement of Q does not have the physical u.m. of flow, because Q was assessed from the regional tracer concentration during apnea, measured by PET in µCi/ml. Black voxels within the sVA image represent outliers excluded from the analysis of sVA (see METHODS).

Experimental Setup

We used a PC-4096 PET scanner (Scanditronix) that imaged 15 contiguous, 6-mm-thick slices of thorax with an in-plane spatial resolution of 6 mm full width half maximum (FWHM). Impedance plethysmography (Respitrace, Non-Invasive Monitoring Systems, Miami Beach, FL) was used to monitor tidal volume and respiratory rate and to ensure that the subject held his or her breath at mean lung volume for the duration of the apnea.

13N2 gas was dissolved into sterile saline solution. The specific activity of the injectate was 0.38 ± 0.05 (means ± SD) mCi/ml and a dose of 8.9 ± 2.3 mCi per scan was injected at a rate of 5 ml/s in a peripheral vein. This corresponded to a radiation exposure of 17.8 ± 4.6 mrem per scan.

Experimental Protocol

The study protocol was approved by the Human Research Committee of our Institutional Review Board. Subjects were randomized to start in the supine or prone position (subjects 1, 2, and 5 started supine, and subjects 3, 4, and 6 started prone). After placement of an intravenous catheter in one antecubital vein, subjects lay recumbent on the PET scanner table with arms abducted and thorax in the PET scanner field up to the level of the axilla. A 10-min transmission scan was then collected, after which the subject was instructed to take two large breaths to total lung capacity and then to resume normal breathing for five additional breaths. During the expiratory phase of the fifth breath, when mean lung volume was reached, as assessed by the Respitrace, the subject was instructed to stop exhalation and remain apneic for 40 s. The lung volume signal from the Respitrace was monitored continuously to ensure the absence of ventilation during the apnea. At the end of this period, spontaneous ventilation was resumed and the subject was encouraged to maintain a regular breathing pattern despite the increased respiratory drive that followed the apnea. The apnea was performed at mean lung volume to ensure that the lung volume during the apnea was the same as the mean lung volume during the washout. This should limit the effect of image registration artifacts on the measurement of sVA.

Simultaneously with the beginning of the apnea, a bolus of 13N2-saline was injected intravenously and a PET scan lasting 3 min and 40 s (40 s of apnea followed by 3 min of washout) was started. Each PET scan consisted of a series of sequential PET images. At the end of the scan, the subject turned to the opposite body position, and the imaging protocol was repeated. Care was taken to image the same cross section of thorax in both body positions by aligning the laser pointer of the PET scanner with skin marks on the chest of the subject.

Four subjects had an additional emission scan in the last of the two body positions (subjects 2 and 5 had two prone scans and 3 and 4 had two supine scans). This third scan, taken ~30 min after the second scan, was used to assess reproducibility of the measurements.

Image Processing and Analysis

The projection data were corrected for nonuniformity of detector response, dead time, random coincidences, attenuation, and scattered radiation. The PET scanner was cross-calibrated with a well scintillation counter by comparing the scanner response from a fluoride-18 solution in a 20-cm cylindrical phantom with the response of the well counter to an aliquot of the same solution. All emission scans were reconstructed by using a conventional filtered back-projection algorithm to an in-plane resolution of 7 mm FWHM.

Selection of voxels for analysis. Lung "masks" were defined by thresholding the transmission scan (4, 24). These masks were then refined to exclude regions corresponding to the main bronchi and large pulmonary vessels. In subject 1 in the supine position, small dependent areas of tracer retention were identified on the last PET image of the washout and were excluded from the lung mask that was used to obtain the V and Q/V measurements.

Because sVA is effectively calculated as the ratio of two variables (i.e., the tracer concentration at the beginning of the washout and the integral of the washout curve; see Eq. 1), the frequency distribution of sVA is particularly prone to have outliers. These outliers would profoundly affect the estimate of V heterogeneity, although they probably reflect the uncertainty of our method in estimating sVA in regions of very low Q/V rather than true extreme values of sVA. Outliers, defined as voxels with sVA < [Q1 - 2(Q3 - Q1)] or sVA > [Q3 + 2(Q3 - Q1)], where Q1 and Q3 are the 25th and 75th percentile of the sVA distribution (26), were excluded from the analysis of V.

PET scans. Emission scans, corrected for tracer radioactive decay (half-life of 13N2: 9.96 min), were low-pass filtered to a length scale of 12 mm and corrected for edge effects (32). After scans were filtered, the size of the effective resolution element corresponded to 12 × 12 × 6 mm. For each voxel within the lung mask, Q, sVA, and Q/sVA were calculated from the tracer kinetics of the PET scan and displayed on a color-coded scale to yield functional images (Fig. 2).

Assessment of vertical ("gravitational") gradients. For each voxel, Q, sVA, Q/sVA, and Fgas were regressed vs. the vertical distance from the most dependent point of the imaged lung field. Gradients, measured from the slope of the corresponding regression line, were expressed in percent per centimeter, relative to the mean (Fig. 3). Negative gradients indicate a decrease of Q, sVA, Q/sVA, or Fgas from the dependent to the nondependent regions (i.e., higher values of Q, sVA, Q/sVA, or Fgas in dependent than in nondependent regions, Fig. 3).


View larger version (30K):
[in this window]
[in a new window]
 
Fig. 3.   Vertical Q gradient. For each PET scan, mean-normalized Q of the voxels included in the imaged lung field was regressed vs. the vertical distance of the voxel from the most dependent point of the imaged lung field. The vertical gradient was measured from the slope of the regression line (black line in the figure) and expressed in %/cm, relative to the mean. Negative gradients indicate a decrease of Q from the dependent to the nondependent regions (i.e., higher Q in dependent than in nondependent regions). In this figure, a gray color scale was used to represent the frequency distribution of the voxels as a function of Q and vertical distance.

Assessment of spatial heterogeneity. Spatial heterogeneity of Q, sVA, and Q/sVA was assessed from the squared coefficient of variation [CV2 = (SD/mean)2] of the respective functional images (29, 32).

Q. To correct the measured Q heterogeneity for the contribution of random imaging noise caused by finite count statistics, an approach similar to the one reported by Venegas et al. (33) was used. The assumption behind this approach is that, during the plateau phase of the apnea, regional tracer concentration remains constant. Thus differences in the tracer concentration measured in a given voxel on the PET images collected during the plateau of the apnea are due to random imaging noise. This assumption is valid if tracer reabsorption into pulmonary venous blood is negligible, a condition that is met in the normal human lung (see APPENDIX and Fig. 1). For each injection of tracer, a number of sequential PET images corresponding to the plateau phase of the apnea were identified by inspection. In the analysis of the data, these images were averaged in different combinations by calculating a voxel-by-voxel duration-weighted mean of the tracer activities of the images included in each combination. For example, if ai,j denotes the tracer activity of the ith voxel (i = 1, ... , N, where N is the number of voxels in the lung mask) in the jth image (j = 1, ... , M, where M is the number of images corresponding to the plateau of the apnea) and tj is the duration of the jth image, then the tracer activity of the ith voxel in the combination image made of, for instance, the first, second and Mth image is
a<SUB>i</SUB>=<FR><NU><LIM><OP>∑</OP><LL>j=1,2,M</LL></LIM> t<SUB>j</SUB>·a<SUB>i,j</SUB></NU><DE><LIM><OP>∑</OP><LL>j=1,2,M</LL></LIM> t<SUB>j</SUB></DE></FR> (2)
In particular, if j in the summations in Eq. 2 ranges from 1 to M, then ai represents the average activity of the ith voxel during the entire plateau phase of the apnea and the corresponding combination image is the Q image (Fig. 2). The CV2 of each combination image (i.e., the CV2 of {ai}i=1, ... , N) was then calculated and regressed vs. the reciprocal of the duration-weighted sum of the mean voxel activity of the images included in each combination. For the same example as above, in which we considered the combination image made of the first, second and Mth image, the duration-weighted sum of the mean voxel activities is
S<SUB>1,2<IT>,M</IT></SUB><IT>=</IT><LIM><OP>∑</OP><LL><IT>j=</IT>1,2<IT>,M</IT></LL></LIM><IT> t<SUB>j</SUB>·</IT><FENCE><FR><NU>1</NU><DE><IT>N</IT></DE></FR>  <LIM><OP>∑</OP><LL><IT>i=</IT>1</LL><UL><IT>N</IT></UL></LIM><IT> a<SUB>i,j</SUB></IT></FENCE> (3)
This sum, which represents the mean activity-time integral of the combination image, is proportional to the number of counts recorded by the PET scanner over the acquisition time of the combination image. The intercept of the regression line (Fig. 4) represents the estimated heterogeneity of Q, had a Q image with an infinite number of counts been acquired. The intercept is therefore a measure of Q heterogeneity, corrected for the contribution of random imaging noise (33). We will denote the value of the intercept as CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP> and the CV2 of the Q image (i.e., the CV2 of {ai}i = 1, ... ,N when j = 1, ... , M in Eq. 2) as CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>,im</SUB><SUP>2</SUP></UP>.


View larger version (12K):
[in this window]
[in a new window]
 
Fig. 4.   "Squared coefficient of variation (CV2) vs. activity" regression line. For each injection of tracer, the PET images from the plateau phase of the apnea were averaged in all possible combinations. The CV2 of each combination image was then calculated and regressed vs. the reciprocal of the mean activity-time integral of the image. The CV2 due to random noise for any PET image obtained from a subject could thus be estimated by multiplying the reciprocal of the mean activity-time integral of the image by the slope of the line obtained from that subject.

The importance of the linear relationship between CV2 and the reciprocal of the mean activity-time integral of a PET image (the "CV2 vs. activity" line, Fig. 4) is that it allows estimation of the CV2 due to random noise for a PET image with arbitrary mean activity and arbitrary duration. Indeed, if the mean activity and the acquisition time of any PET image from a subject are known, the CV2 due to random noise can be estimated by multiplying the reciprocal of the activity-time integral of the image by the slope of the line obtained from that subject. The CV2 vs. activity line was calculated for each subject in each body position and was used to correct the correspondent heterogeneity measurement.

CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP> was subdivided in the amount explained by the vertical gradient (CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>,g</SUB><SUP>2</SUP></UP>) and the residual heterogeneity (CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>,r</SUB><SUP>2</SUP></UP> = CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP> - CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>,g</SUB><SUP>2</SUP></UP>). CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>,g</SUB><SUP>2</SUP></UP> was calculated, by linear regression, from the Q image. The ratio R<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP> = (CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>,g</SUB><SUP>2</SUP></UP>)/(CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP>) thus represents the fraction of Q heterogeneity explained by the vertical gradient.

Q/V and V. For each PET scan, the duration-weighted sum of the mean activity of the PET images taken during tracer washout was calculated as
S<IT>=</IT><LIM><OP>∑</OP><LL><IT>k=</IT>1</LL><UL><IT>L</IT></UL></LIM><IT> t<SUB>k</SUB>·</IT><FENCE><FR><NU>1</NU><DE><IT>N</IT></DE></FR>  <LIM><OP>∑</OP><LL><IT>i=</IT>1</LL><UL><IT>N</IT></UL></LIM><IT> a<SUB>i,k</SUB></IT></FENCE> (4)
where k = 1, ... , L indexes the washout images. This sum is proportional to the number of counts recorded by the PET scanner during the washout phase [note that this sum is conceptually equivalent to the duration-weighted sum of the mean activity of the apnea images (Eq. 3) that was used to estimate the CV2 vs. activity regression line]. The reciprocal of this sum was then multiplied by the slope of the CV2 vs. activity regression line obtained from the apnea images of the same scan (Fig. 4). This yielded an estimate of the CV2 caused by random imaging noise for the PET image corresponding to the duration-weighted sum of the washout images, i.e., the Q/sVA image. The CV2 due to random imaging noise was then subtracted from the CV2 of the Q/sVA image (CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>/s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>,im</SUB><SUP>2</SUP></UP>) to obtain a measure of Q/sVA heterogeneity (CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>/s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP>). As for CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP>, CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>/s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP> was subdivided in the amount explained by the vertical gradient (CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>/s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>,g</SUB><SUP>2</SUP></UP>) and the residual heterogeneity (CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>/s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>,r</SUB><SUP>2</SUP></UP>), and the ratio R<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>/s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP> = (CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>/s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>,g</SUB><SUP>2</SUP></UP>)/(CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>/s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP>) was calculated.

To estimate the contribution of random noise to the measurement of sVA, we used a Monte Carlo approach. Because sVA was computed as the ratio of the Q image and the Q/sVA image, a ratio that corresponds to the reciprocal of the MRT (see Eq. 1), the parameters describing the distribution of noise in the Q image and in the Q/sVA image were used for the Monte Carlo simulation as follows. A normally distributed random variable (r.v.), with mean equal to the mean activity of the Q image and CV2 equal to the product of the reciprocal of the activity-time integral of the Q image (Eq. 3 for j = 1, ... , M) by the slope of the CV2 vs. activity regression line, was generated. The distribution of this r.v. characterizes the distribution of random noise in the Q image, under the assumption that such noise can be modeled as a normally distributed stochastic process. Similarly, a second normally distributed r.v., with mean equal to the mean activity of the Q/sVA image and CV2 equal to the product of the reciprocal of the activity-time integral of the Q/sVA image (Eq. 4) by the slope of the CV2 vs. activity regression line, was generated. The distribution of this second r.v. represents the distribution of random noise in a PET image with the same mean activity as the Q/sVA image. Because sVA was computed as the ratio of Q and Q/sVA, a third r.v. was generated as the ratio of the first and the second r.v. The distribution of this third r.v. represents the effect that random noise in the Q and Q/sVA images has on the distribution of sVA. One hundred thousand empiric realizations of this r.v. were generated by computer simulation (Matlab, The MathWorks, Natick, MA), and their CV2 was calculated. This CV2 represents an estimate of the sVA heterogeneity that could be attributed to random imaging noise. This was subtracted from the CV2 of the sVA image (CV<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>,im</SUB><SUP>2</SUP></UP>) to obtain a measure of sVA heterogeneity (CV<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP>). As with the previous variables, CV<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP> was subdivided in CV<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>,g</SUB><SUP>2</SUP></UP> and CV<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>,r</SUB><SUP>2</SUP></UP>, and the ratio R<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP> = (CV<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>,g</SUB><SUP>2</SUP></UP>)/(CV<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP>) was calculated.

Statistical Analysis

Least squares linear regression was used to estimate all regression relationships. When calculating the regression of the supine-to-prone difference in the Fgas gradient vs. the supine-to-prone difference in the Q gradient (see below), the constraint of a zero intercept was imposed. This constraint reflected the assumption that, if the Q gradients were equal in the supine and prone position, then the Fgas gradients would also be equal. This assumption was justified by the consideration that, if the lung behaved as a passive structure under the influence of hydrostatic forces, then the gradients of Q and lung density (and therefore of Fgas) would not differ between the two body positions.

Two-tailed Student's t-test was used to assess any significant difference from zero and to compare the results in the supine and prone position (t-test for paired data). Statistical significance was set at P < 0.05. To control the "experiment-wise" type I error and limit the number of statistical comparisons, the following approach was used. It was decided a priori that the primary comparisons of interest were those involving the gradients and the total heterogeneity. Only if these were statistically significant, then comparisons involving the different components of the total heterogeneity were performed as secondary endpoints. Values are presented as means ± SD.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

Vertical Gradients

The Q gradient was different from zero in both body positions, without a significant difference between supine and prone (Table 1). There was, however, large intersubject variability in the Q gradients, ranging from +0.5%/cm to -10.6%/cm (Fig. 5). Only in one subject (subject 6 in prone position) was there a slightly positive Q gradient; otherwise, Q was consistently greater in dependent than in nondependent lung regions (i.e., Q had a negative gradient).

                              
View this table:
[in this window]
[in a new window]
 
Table 1.   Vertical gradients



View larger version (20K):
[in this window]
[in a new window]
 
Fig. 5.   Vertical gradients of Q, sVA, and Q/sVA in the supine (S) and prone (P) position for the 6 subjects studied. Gradients are expressed in %/cm, relative to the mean. Negative gradients indicate a decrease from the dependent to the nondependent regions (i.e., higher values of Q, sVA, or Q/sVA in dependent regions; see also Fig. 3).

The Fgas gradient was positive in both body positions (Table 1). The difference in Fgas gradient between the supine and prone positions (Delta grad Fgas) was negatively correlated (r = -0.778) with the difference in Q gradient (Delta grad Q), as shown in Fig. 6. The slope of the regression line was -0.162 (95% confidence interval: -0.313 to -0.012).


View larger version (11K):
[in this window]
[in a new window]
 
Fig. 6.   Regression of the supine-to-prone difference in fractional gas content gradient (Delta grad Fgas) vs. the supine-to-prone difference in Q gradient (Delta grad Q). Each point represents a subject. The 95% confidence interval of the slope of the regression line did not include -1 (dashed line).

There was a significant, negative sVA gradient in both supine and prone positions (Table 1), and all subjects had greater sVA in dependent regions (Fig. 5).

The average Q/sVA gradient was not significantly different from zero in either body position (Table 1). However, two subjects in both body positions (2 and 6) and one subject in the supine position (3) had positive Q/sVA gradients, consistent with the finding that, in these cases, the Q gradient was less pronounced than the sVA gradient (Fig. 5). Consequently, in these cases, the Q/sVA was higher in nondependent than in dependent regions. In contrast, negative Q/sVA gradients (subjects 1, 3 in prone position, 4, and 5) corresponded to cases in which the magnitude of the Q gradient was greater than the magnitude of the sVA gradient. This implies that, in these cases, the fraction of Q going to dependent lung regions was greater than the fraction of sVA in these regions (Fig. 5).

The difference in vertical gradients between the first and second scan obtained in the same body position (n = 4) was 0.4 ± 1.1%/cm for Q, 0.9 ± 1.2%/cm for sVA, and -0.5 ± 2.3%/cm for Q/sVA.

Heterogeneity

There was no statistically significant difference in Q heterogeneity, CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP>, between the supine and prone positions (Table 2). The fraction of Q heterogeneity explained by the vertical gradient, R<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP>, was significantly different from zero and averaged 0.22 and 0.26 in, respectively, the supine and prone position (Table 2). There was, however, large intersubject variability in R<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP>, which varied between 0.01 and 0.57 (Fig. 7). Despite this variability, there was internal consistency between the heterogeneity and the gradient data. Subjects with pronounced vertical Q gradients (subject 4 and subject 3 in prone position) had higher values of CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>,g</SUB><SUP>2</SUP></UP> and R<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP> than subjects with small gradients (Fig. 5, Fig. 7). The residual Q heterogeneity, CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>,r</SUB><SUP>2</SUP></UP>, was significantly different from zero (Table 2, Fig. 7).

                              
View this table:
[in this window]
[in a new window]
 
Table 2.   Spatial heterogeneity



View larger version (18K):
[in this window]
[in a new window]
 
Fig. 7.   Heterogeneity (CV2) of Q and sVA in the 6 subjects studied. For each subject, the bar on the left represents heterogeneity in the supine position (S), whereas the bar on the right corresponds to the prone position (P). The total heterogeneity (full height of the bar) was subdivided in the amount explained by the vertical gradient (white portion of the bar) and the residual heterogeneity (gray portion of the bar). The ratio of the white portion of the bar to the total height of the bar yields a visual interpretation of the fraction of heterogeneity explained by the vertical gradient (i.e., R<UP><SUB><A><AC>Q</AC><AC>˙</AC></A></SUB><SUP>2</SUP></UP> and R<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP>).

V heterogeneity, CV<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP>, was on average higher in the supine than in the prone position (Table 2). The fraction of V heterogeneity explained by the vertical gradient, R<UP><SUB>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP>, was different from zero in both body positions, but significant residual V heterogeneity persisted after subtraction of the component due to the vertical gradient (Table 2, Fig. 7).

There was no difference in Q/sVA heterogeneity, CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>/s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP>, between the supine and prone positions (Table 2). Because we did not find significant Q/sVA gradients, we did not test the components of CV<UP><SUB><A><AC>Q</AC><AC>˙</AC></A>/s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></SUB><SUP>2</SUP></UP>.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

Comparison to Other Methods and Critique of the Technique

Several methods have been employed to assess regional Q in animals and humans (5, 7, 8, 10, 15, 17-19, 27). Because 13N2 has a first-pass retention of virtually 100% in aerated lung regions (23, 30), the 13N2-saline bolus infusion technique is conceptually similar to methods that utilize tracers that are retained in the pulmonary microvasculature, such as microspheres or macroaggregated albumin. The advantage of using 13N2, though, is that N2, like O2 and CO2, is a low-molecular-weight gas that is dissolved in the blood. Its distribution through the pulmonary capillary network should therefore resemble that of the respiratory gases more closely than the bulkier (10-100 µm in diameter) microspheres or albumin macroaggregates. Furthermore, compared with the macroaggregated albumin method, the 13N2 infusion technique reduces radiation exposure for the subject because 13N2 has a shorter half-life and is eliminated more rapidly than 99mTc.

For the 13N2 content during apnea to accurately reflect regional Q, two basic conditions need to be met. First, the apnea itself needs not to influence the regional distribution of Q. Second, tracer removal by the bloodstream during the apnea needs to be negligible. Because tracer distribution into the lung is virtually complete within the first 10 s since the start of the infusion, the effect of the apnea should be minimal. To further limit the potential effect of the apnea on regional Q, the subjects took two inspirations to total lung capacity five breaths before the beginning of the apnea. This maneuver should prevent the development of alveolar atelectasis and reduce the level of hypercapnia reached during the apnea. Because of the low solubility of 13N2 in blood, tracer reabsorption into pulmonary venous blood during the apnea is minimal. In normally aerated and perfused lungs, tracer reabsorption can be expected to result in a <2% drop in tracer concentration during the apnea (see APPENDIX).

The MRT of a gas tracer has been previously used to characterize the regional distribution of V (31). We assessed V from the reciprocal of the MRT during the washout. This approach has the following implications. First, because 13N2 is removed from the lung also by the pulmonary venous blood, the MRT is shorter than if 13N2 had been eliminated exclusively by ventilation. However, in a normally aerated and perfused lung, this should lead to an overestimation of specific ventilation (i.e., ventilation per unit of compartment volume) by only 1.5% (see APPENDIX). Therefore, the reciprocal of the MRT appears to be a reasonable estimate of specific ventilation. Second, the reciprocal of the MRT represents the rate at which the concentration of 13N2 within the voxel decreases. If the voxel behaves as a well-mixed single compartment, the rate of decrease of the concentration of 13N2 within the voxel equals the rate of decrease of the alveolar concentration of 13N2 (although the alveolar and the voxel concentration of 13N2 may differ because lambda water/air = 0.015 for 13N2). Therefore, the reciprocal of the MRT effectively represents specific alveolar ventilation, sVA, and the fact that we do not measure the alveolar concentration of 13N2 but only its concentration within the voxel is immaterial. Third, because sVA reflects the net rate of regional gas transfer and intrinsically accounts for the effect of reinspiration of gas from the common dead space, the topographical distribution of sVA should be closely related to regional gas exchange.

Gradients and Heterogeneity

Q. In this study, we found vertical Q gradients favoring dependent lung regions in both supine and prone positions. These gradients are similar to those reported previously in spontaneously breathing humans with magnetic resonance imaging (28) and single-photon-emission computed tomography (19) and corroborate the results of earlier studies that employed low-resolution planar imaging techniques (11, 12). Whereas it is generally accepted that Q is highest in dependent lung regions in the supine position, the vertical dependence of Q in the prone position is more controversial. Animal data, mostly from quadrupeds, indicate that there may be structural factors (3) or biochemical modulators (20) of the pulmonary vasculature that favor perfusion in dorsal lung regions. This is consistent with the results of Treppo et al. (29), who found a marked vertical gradient in supine (where both structural/biochemical factors and gravity would tend to increase Q in dependent, dorsal, regions) but not in prone dogs. These factors, by contrasting the action of gravity in the prone position, would serve the purpose of rendering the distribution of Q more homogeneous in quadrupeds. Recent studies, employing single-photon-emission computed tomography, suggested that also in humans the vertical Q gradient is either less or nil in the prone position (17, 18). Our data do not corroborate these findings, and caution on the extent to which results of studies that showed greater dorsal Q in prone quadrupeds (7) can be extended to biped humans. Our results are in line with those of Amis et al. (2) and Jones at al. (10), which also showed higher Q in dependent (ventral) than in nondependent (dorsal) regions in prone humans. Our results do, however, show a wide variability in individual behaviors that may partly explain the discrepant findings in the literature. The fact that vertical gradients were reproducible on repeated measurements in the same body position argues against the possibility that this variability was a measurement artifact. The fact that the change in gradient between the supine and prone positions was not dictated by the order of randomization (i.e., gradients were not systematically higher, or lower, in the first than in the second position, Fig. 5) virtually excludes the possibility that the intersubject variability in the change in gradient was merely the result of a systematic difference between the gradients in the first and the second body position. Whereas supine gradients were very similar in five of the six subjects, prone gradients were more disperse (Fig. 5). This suggests that, on turning prone, Q redistributed toward ventral, dependent lung regions in all subjects, but the extent to which this occurred was much greater in some subjects (e.g., subject 4) than in others (e.g., subject 6). Future studies will be needed to address whether a lesser degree of Q redistribution toward dependent, ventral regions is associated with a favorable response to prone positioning in patients with ARDS, in whom lung densities have been shown to shift to ventral regions in the prone position (13).

Interestingly, the supine-to-prone difference in the Fgas gradient (Delta grad Fgas) was negatively correlated with the supine-to-prone difference in Q gradient (Delta grad Q), but the 95% confidence interval of the slope of the regression line did not include -1. If perfusion per unit of lung tissue, and therefore perfusion per alveolus, were constant and the redistribution of Q between the supine and prone positions were merely a consequence of redistribution of lung tissue, then the change in the mean-normalized Q gradient would equal the change in the mean-normalized lung density gradient. As a result, Delta grad Q would be equal in magnitude but opposite in sign to Delta grad Fgas, and the slope of the regression line in Fig. 6 would be -1. Our finding instead suggests that redistribution of lung tissue cannot be the only explanation for the redistribution of Q that is associated with body position changes and that perfusion per unit of lung tissue, and therefore perfusion per alveolus, increases when a nondependent region of lung becomes dependent as a result of the position change. This is consistent with a decreased vascular resistance of dependent lung regions compared with nondependent regions.

Over the past decade, animal data obtained with the intravenous injection of microspheres have suggested that vertical gradients explain only a minor fraction (i.e., ~5%) of Q heterogeneity in the supine or prone position (7). To avoid an underestimation of the relative importance of vertical gradients, because of an overestimation of the total heterogeneity, we applied a method that allows estimation of the heterogeneity due to finite count statistics, which can then be used to correct the total measured heterogeneity (33). It is important to acknowledge, though, that this method does not correct for systematic sources of imaging noise such as reconstruction and registration artifacts. On average, the vertical gradient explained ~24% of Q heterogeneity in both positions. Both the total heterogeneity and the residual heterogeneity were similar in the supine and prone position (Table 2), suggesting that, in humans, Q is not systematically more uniform in the prone position. There were marked differences, though, in individual behaviors, and in four subjects Q heterogeneity was lower in the prone position (Fig. 7). For subject 4, vertical gradients explained the majority of Q heterogeneity (57% in both positions), whereas for subject 6 the contribution of the vertical gradient to Q heterogeneity in the prone position was insignificant. These differences may reflect a variable contribution of structural factors and active regulation of the pulmonary vasculature to the topographical distribution of Q.

V. In both the supine and prone positions, we found vertical sVA gradients that favored dependent lung regions (Table 1, Fig. 5). These results are consistent with early reports (11, 22) and with more recent PET studies (4) but differ from the observations of other authors that reported V to be uniform (1) or even greater in nondependent (dorsal) than in dependent (ventral) regions in the prone position (17, 19, 21). This discrepancy of results in the literature is not surprising in view of the different factors that can affect the topographical distribution of V in awake, spontaneously breathing humans. These include the pleural pressure gradient and the pattern of contraction of the respiratory muscles. The vertical gradient of pleural pressure has been suggested to be less in the prone position (34). This would favor a more uniform distribution of V when prone. Although the difference was not statistically significant, we did find, on average, lower sVA gradients in prone subjects (Table 1), and sVA heterogeneity was significantly smaller in the prone than in the supine position (Table 2, Fig. 7). These findings are consistent with a more uniform distribution of V in the prone position.

In awake, spontaneously breathing humans, the pattern of respiratory muscle contraction, which may differ among subjects, can also influence regional V (25) and explain part of the discrepancies in the literature. This is consistent with our results, which show that only a small fraction of V heterogeneity is explained by the vertical gradient and that residual, "nongravitational," heterogeneity is significant.

Finally, our findings show that estimates of the regional distribution of V are profoundly affected by the occurrence of gas trapping. Dependent, crescent-shaped areas of very low V/Q were reported, in healthy humans, by Rhodes et al. (24) and attributed to low V. We detected areas of tracer retention during the washout (i.e., areas of gas trapping) in the dependent lung regions of subject 1 in the supine position. This finding suggests that airway closure may occur in dependent lung regions during tidal breathing. Despite the fact that these areas represented only 9% of the imaged lung field, the sVA gradient decreased from 0.9 to -1.5%/cm and the sVA heterogeneity decreased from 0.26 to 0.16 after these areas were removed from the lung mask of this subject. In contrast to methods like ours that use the intravenous infusion of inert insoluble gases, methods that image V by inhalation of tracer will not reveal areas of gas trapping, because inhaled tracer cannot reach the airspace distal to closed airways. These implications of different methods may represent a further explanation for the discordant data on the regional distribution of V reported in the literature.

Q/V. We did not find significant Q/sVA gradients in either body position (Table 1), and the spatial heterogeneity of Q/sVA was similar in the supine and prone position (Table 2). This is consistent with the rest of our results and suggests that, on average, the vertical gradient of Q match the gradient of V. These findings contrast with the results of animal studies that have shown both the vertical gradient and the heterogeneity of the V/Q ratio to be substantially reduced in the prone position (15, 29). Interspecies differences and the fact that animals were studied during anesthesia and mechanical ventilation may account for this discrepancy.

In conclusion, we have shown that regional Q and V can be assessed noninvasively, in humans, with the 13N2-saline bolus infusion technique and PET. Our results suggest that both Q and V are distributed preferentially to dependent lung regions in the supine as well as in the prone position and that, in awake spontaneously breathing humans, the prone position does not offer any systematic advantage, in terms of V-to-Q matching, compared with the supine position. However, individual differences, especially in the distribution of Q, may explain why the prone position is effective in improving gas exchange in some subjects. The next step is to identify patterns of distribution of Q and V that may be predictive of a favorable response to prone positioning.


    APPENDIX
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

The purpose of this appendix is to present a one-compartment model for a voxel of lung to 1) elucidate the rationale of the measurements of sVA and Q/sVA obtained with the 13N2-saline bolus infusion method, and 2) quantify the effect of 13N2 reabsorption into pulmonary venous blood on Q, sVA, and Q/sVA measurements.

Glossary


VA   Alveolar gas volume
 VA   Alveolar ventilation
sVA   Specific alveolar ventilation (sVA = VA/VA)
 Q   Pulmonary perfusion
Co   Alveolar concentration (activity) of 13N2 immediately after arrival of the bolus of tracer (proportional to Q; Ref. 14)
c(t)   Alveolar concentration (activity) of 13N2 at time t
Cw   Alveolar concentration (activity) of 13N2 at the beginning of the washout
 lambda water/air   partition coefficient for 13N2 (lambda  = 0.015 at 37°C)

The rate of change of c(t) is given by the following first-order differential equation (14)
<FR><NU>dc(<IT>t</IT>)</NU><DE>d<IT>t</IT></DE></FR><IT>=</IT>−c(<IT>t</IT>) <FR><NU><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></NU><DE>V<SC>a</SC></DE></FR><IT>−</IT>c(<IT>t</IT>)<IT>&lgr; </IT><FR><NU><A><AC>Q</AC><AC>˙</AC></A></NU><DE>V<SC>a</SC></DE></FR> (A1)
The solution is
c(<IT>t</IT>)<IT>=</IT>C<SUB>o</SUB><IT>e</IT><SUP><IT>−</IT>(<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>/V<SC>a</SC><IT>+&lgr;</IT>(<A><AC>Q</AC><AC>˙</AC></A>/V<SC>a</SC>))<IT>t</IT></SUP><IT>=</IT>C<SUB>o</SUB><IT>e</IT><SUP><IT>−</IT>(s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><IT>+&lgr;</IT>(<A><AC>Q</AC><AC>˙</AC></A>/V<SC>a</SC>))<IT>t</IT></SUP> (A2)
Integrating c(t) from the beginning of the washout [c(0) = Cw] to infinity yields
<LIM><OP>∫</OP><LL>0</LL><UL>∞</UL></LIM> C<SUB>w</SUB><IT>e</IT><SUP><IT>−</IT>(s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><IT>+&lgr;</IT>(<A><AC>Q</AC><AC>˙</AC></A>/V<SC>a</SC>))<IT>t</IT></SUP>d<IT>t=</IT><FR><NU>C<SUB>w</SUB></NU><DE>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><IT>+&lgr;</IT><FR><NU><A><AC>Q</AC><AC>˙</AC></A></NU><DE>V<SC>a</SC></DE></FR></DE></FR><IT> ∝ </IT><FR><NU><A><AC>Q</AC><AC>˙</AC></A></NU><DE>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><IT>+&lgr; </IT><FR><NU><A><AC>Q</AC><AC>˙</AC></A></NU><DE>V<SC>a</SC></DE></FR></DE></FR> (A3)
To quantify the effect of 13N2 reabsorption on the measurements of Q, sVA, and Q/sVA, we will use global values of sVA and Q/VA of normal human lungs in which Q = VA = 5 l/min and VA = 2.5 liters. This yields sVA = Q/VA = 0.033/s.

During the apnea, sVA = 0. According to Eq. A2, the tracer concentration is expected to decrease by 1.5% over the 30 s of apnea that follow the arrival of the bolus of tracer. This warrants the assumption that c(t) is virtually constant during the plateau phase of the apnea and that Co approx  Cw. Because Coproportional to Q (14) and Co approx  Cw, then Cwproportional to Q (Eq. A3).

During the washout, ventilation is resumed. Equation A2 shows that the reciprocal of the time constant of the tracer washout curve is related to sVA. Equation A3 shows that the integral of the tracer washout curve is related to Q/sVA.

The error in the estimation of sVA and Q/sVA due to the fact that we neglect tracer reabsorption by the pulmonary venous blood and estimate sVA from [sVA + lambda (Q/VA)] is
ϵ=<FR><NU>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><IT>+&lgr; </IT><FR><NU><A><AC>Q</AC><AC>˙</AC></A></NU><DE>V<SC>a</SC></DE></FR><IT>−</IT>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></NU><DE>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></DE></FR><IT>=&lgr;=</IT>1.5<IT>%</IT>


    ACKNOWLEDGEMENTS

We thank M. Martinez for contributing to the analysis of the data; S. Barrow, S. Weise, A. Bruce, and the "MGH Cyclotron Group" for technical assistance; Dr. T. Richter and K. O'Neill for insightful comments; and Dr. W. M. Zapol for continued support with this project.


    FOOTNOTES

This study was supported by National Heart, Lung, and Blood Institute Grants GM-07592, HL-38267, and HL-68011.

Address for reprint requests and other correspondence: G. Musch, Dept. of Anesthesia and Critical Care, CLN 309, Massachusetts General Hospital, 32 Fruit St., Boston, MA 02114 (E-mail: gmusch{at}partners.org).

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.

10.1152/japplphysiol.00223.2002

Received 15 March 2002; accepted in final form 3 July 2002.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

1.   Amis, TC, Jones HA, and Hughes JMB Effect of posture on inter-regional distribution of pulmonary ventilation in man. Respir Physiol 56: 145-167, 1984.

2.   Amis, TC, Jones HA, and Hughes JMB Effect of posture on inter-regional distribution of pulmonary perfusion and VA/Q ratios in man. Respir Physiol 56: 169-182, 1984.

3.   Beck, KC, and Rehder K. Differences in regional vascular conductances in isolated dog lungs. J Appl Physiol 61: 530-538, 1986.

4.   Brudin, LH, Rhodes CG, Valind SO, Jones T, Jonson B, and Hughes JMB Relationships between regional ventilation and vascular and extravascular volume in supine humans. J Appl Physiol 76: 1195-1204, 1994.

5.   Brudin, LH, Rhodes CG, Valind SO, Jones T, and Hughes JMB Interrelationships between regional blood flow, blood volume and ventilation in supine humans. J Appl Physiol 76: 1205-1210, 1994.

6.   Gattinoni, L, Tognoni G, Pesenti A, Taccone P, Mascheroni D, Labarta V, Malacrida R, Di Giulio P, Fumagalli R, Pelosi P, Brazzi L, and Latini R. Effect of prone positioning on the survival of patients with acute respiratory failure. N Engl J Med 345: 568-573, 2001.

7.   Glenny, RW, Lamm WJ, Albert RK, and Robertson HT. Gravity is a minor determinant of pulmonary blood flow distribution. J Appl Physiol 71: 620-629, 1991.

8.   Hakim, TS, Lisbona R, and Dean GW. Gravity-independent inequality in pulmonary blood flow in humans. J Appl Physiol 63: 1114-1121, 1987.

9.   Harris, RS, Willey-Courand DB, Head CA, Galletti GG, Call DM, and Venegas JG. Regional VA,Q, and VA/Q during PLV: effects of nitroprusside and inhaled nitric oxide. J Appl Physiol 92: 297-312, 2002.

10.   Jones, AT, Hansell DM, and Evans TW. Pulmonary perfusion in supine and prone positions: an electron-beam computed tomography study. J Appl Physiol 90: 1342-1348, 2001.

11.   Kaneko, K, Milic-Emili J, Dolovich MB, Downson A, and Bates DV. Regional distribution of ventilation and perfusion as a function of body position. J Appl Physiol 21: 767-777, 1966.

12.   Landmark, SJ, Knopp TJ, Rehder K, and Sessler AD. Regional pulmonary perfusion and V/Q in awake and anesthetized-paralyzed man. J Appl Physiol 43: 993-1000, 1977.

13.   Langer, M, Mascheroni D, Marcolin R, and Gattinoni L. The prone position in ARDS patients. A clinical study. Chest 94: 103-107, 1988.

14.   Mijailovich, SM, Treppo S, and Venegas JG. Effects of lung motion and tracer kinetics corrections on PET imaging of pulmonary function. J Appl Physiol 82: 1154-1162, 1997.

15.   Mure, M, Domino KB, Lindahl SGE, Hlastala MP, Altemeier WA, and Glenny RW. Regional ventilation-perfusion distribution is more uniform in the prone position. J Appl Physiol 88: 1076-1083, 2000.

16.   Mure, M, and Lindahl SGE Prone position improves gas exchange---but how? Acta Anaesthesiol Scand 45: 150-159, 2001.

17.   Mure, M, Nyren S, Jacobsson H, Larsson SA, and Lindahl SGE High continuous positive airway pressure level induces ventilation/perfusion mismatch in the prone position. Crit Care Med 29: 959-964, 2001.

18.   Nyren, S, Mure M, Jacobsson H, Larsson SA, and Lindahl SGE Pulmonary perfusion is more uniform in the prone than in the supine position: scintigraphy in healthy humans. J Appl Physiol 86: 1135-1141, 1999.

19.   Orphanidou, D, Hughes JMB, Myers MJ, Al-Suhali AR, and Henderson B. Tomography of regional ventilation and perfusion using krypton 81m in normal subjects and asthmatic patients. Thorax 41: 542-551, 1986.

20.   Pelletier, N, Robinson NE, Kaiser L, and Derksen FJ. Regional differences in endothelial function in horse lungs: possible role in blood flow distribution? J Appl Physiol 85: 537-542, 1998.

21.   Rehder, K, Knopp TJ, and Sessler AD. Regional intrapulmonary gas distribution in awake and anesthetized-paralyzed prone man. J Appl Physiol 45: 528-535, 1978.

22.   Rehder, K, Sessler AD, and Rodarte JR. Regional intrapulmonary gas distribution in awake and anesthetized-paralyzed man. J Appl Physiol 42: 391-402, 1977.

23.   Rhodes, CG, Valind SO, Brudin LH, Wollmer PE, Jones T, and Hughes JMB Quantification of regional V/Q ratios in humans by the use of PET. I. Theory. J Appl Physiol 66: 1896-1904, 1989.

24.   Rhodes, CG, Valind SO, Brudin LH, Wollmer PE, Jones T, Buckingham PD, and Hughes JMB Quantification of regional V/Q ratios in humans by the use of PET. II. Procedure and normal values. J Appl Physiol 66: 1905-1913, 1989.

25.   Roussos, CS, Fixley M, Genest J, Cosio M, Kelly S, Martin RR, and Engel LA. Voluntary factors influencing the distribution of inspired gas. Am Rev Respir Dis 116: 457-467, 1977.

26.   Sachs, L. Angewandte Statistik. Berlin: Springer-Verlag, 2002, p. 366.

27.   Schuster, DP, Kaplan JD, Gauvain K, Welch MJ, and Markham J. Measurement of regional pulmonary blood flow with PET. J Nucl Med 36: 371-377, 1995.

28.   Stock, KW, Chen Q, Levin D, Hatabu H, and Edelman RR. Demonstration of gravity-dependent lung perfusion with contrast-enhanced magnetic resonance imaging. J Magn Reson Imaging 9: 557-561, 1999.

29.   Treppo, S, Mijailovich SM, and Venegas JG. Contributions of pulmonary perfusion and ventilation to heterogeneity in V/Q measured by PET. J Appl Physiol 82: 1163-1176, 1997.

30.   Valind, SO, Wollmer PE, and Rhodes CG. Application of positron emission tomography in the lung. In: Positron Emission Tomography, edited by Reivich M, and Alavi A.. New York: Liss, 1985, p. 387-412.

31.   Venegas, JG, Custer J, Kamm RD, and Hales CA. Relationship for gas transport during high-frequency ventilation in dogs. J Appl Physiol 59: 1539-1547, 1985.

32.   Venegas, JG, and Galletti GG. Low-pass filtering, a new method of fractal analysis: application to PET images of pulmonary blood flow. J Appl Physiol 88: 1365-1373, 2000.

33.   Venegas, JG, Treppo S, and Mijailovich SM. Contributions of statistical noise to spatial heterogeneity of PET images of pulmonary function. In: Medical Imaging 1994: Physiology and Function from Multidimensional Images (Proceedings), edited by Hoffman EA, and Acharya RS.. Newport Beach, CA: SPIE-T Int. Soc. for Optical Engineering, 1994, vol. 2168, p. 361-368.

34.   Yang, QH, Kaplowitz MR, and Lai-Fook SJ. Regional variations in lung expansion in rabbits: prone vs. supine positions. J Appl Physiol 67: 1371-1376, 1989.


J APPL PHYSIOL 93(5):1841-1851
8750-7587/02 $5.00 Copyright © 2002 the American Physiological Society



This article has been cited by other articles:


Home page
J. Appl. Physiol.Home page
R. S. Harris, T. Winkler, G. Musch, M. F. Vidal Melo, T. Schroeder, N. Tgavalekos, and J. G. Venegas
The prone position results in smaller ventilation defects during bronchoconstriction in asthma
J Appl Physiol, July 1, 2009; 107(1): 266 - 274.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
S. Verbanck, Y. Kerckx, D. Schuermans, C. de Bisschop, H. Guenard, R. Naeije, W. Vincken, and A. Van Muylem
The effect of posture-induced changes in peripheral nitric oxide uptake on exhaled nitric oxide
J Appl Physiol, May 1, 2009; 106(5): 1494 - 1498.
[Abstract] [Full Text] [PDF]


Home page
RadiologyHome page
B. Driehuys, H. E. Moller, Z. I. Cleveland, J. Pollaro, and L. W. Hedlund
Pulmonary Perfusion and Xenon Gas Exchange in Rats: MR Imaging with Intravenous Injection of Hyperpolarized 129Xe
Radiology, April 29, 2009; (2009) 2513081550.
[Abstract] [Full Text]


Home page
J. Appl. Physiol.Home page
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]


Home page
J. Appl. Physiol.Home page
L. Campana, J. Kenyon, S. Zhalehdoust-Sani, Y.-S. Tzeng, Y. Sun, M. Albert, and K. R. Lutchen
Probing airway conditions governing ventilation defects in asthma via hyperpolarized MRI image functional modeling
J Appl Physiol, April 1, 2009; 106(4): 1293 - 1300.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
M. Hughes and J. B. West
Point:Counterpoint: Gravity is/is not the major factor determining the distribution of blood flow in the human lung
J Appl Physiol, May 1, 2008; 104(5): 1531 - 1533.
[Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
G. Musch, G. Bellani, M. F. Vidal Melo, R. S. Harris, T. Winkler, T. Schroeder, and J. G. Venegas
Relation between Shunt, Aeration, and Perfusion in Experimental Acute Lung Injury
Am. J. Respir. Crit. Care Med., February 1, 2008; 177(3): 292 - 300.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
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 page
Eur Respir JHome page
N. T. Tgavalekos, G. Musch, R. S. Harris, M. F. Vidal Melo, T. Winkler, T. Schroeder, R. Callahan, K. R. Lutchen, and J. G. Venegas
Relationship between airway narrowing, patchy ventilation and lung mechanics in asthmatics
Eur. Respir. J., June 1, 2007; 29(6): 1174 - 1181.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
R. S. Harris and D. P. Schuster
Visualizing lung function with positron emission tomography
J Appl Physiol, January 1, 2007; 102(1): 448 - 458.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
C. Dehnert, F. Risse, S. Ley, T. A. Kuder, R. Buhmann, M. Puderbach, E. Menold, D. Mereles, H.-U. Kauczor, P. Bartsch, et al.
Magnetic Resonance Imaging of Uneven Pulmonary Perfusion in Hypoxia in Humans
Am. J. Respir. Crit. Care Med., November 15, 2006; 174(10): 1132 - 1138.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
D. Chon, K. C. Beck, R. L. Larsen, H. Shikata, and E. A. Hoffman
Regional pulmonary blood flow in dogs by 4D-X-ray CT
J Appl Physiol, November 1, 2006; 101(5): 1451 - 1465.
[Abstract] [Full Text] [PDF]


Home page
Proc Am Thorac SocHome page
E. A. Hoffman, B. A. Simon, and G. McLennan
State of the Art. A Structural and Functional Assessment of the Lung via Multidetector-Row Computed Tomography: Phenotyping Chronic Obstructive Pulmonary Disease
Proceedings of the ATS, August 1, 2006; 3(6): 519 - 532.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
R. S. Harris, T. Winkler, N. Tgavalekos, G. Musch, M. F. V. Melo, T. Schroeder, Y. Chang, and J. G. Venegas
Regional Pulmonary Perfusion, Inflation, and Ventilation Defects in Bronchoconstricted Patients with Asthma
Am. J. Respir. Crit. Care Med., August 1, 2006; 174(3): 245 - 253.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Roentgenol.Home page
A. J. Yun, P. Y. Lee, and A. N. Gerber
Integrating systems biology and medical imaging: understanding disease distribution in the lung model.
Am. J. Roentgenol., April 1, 2006; 186(4): 925 - 930.
[Abstract] [Full Text] [PDF]


Home page
Exp PhysiolHome page
M. H. Tawhai, K. S. Burrowes, and E. A. Hoffman
Computational models of structure-function relationships in the pulmonary circulation and their validation
Exp Physiol, March 1, 2006; 91(2): 285 - 293.
[Abstract] [Full Text] [PDF]


Home page
Proc Am Thorac SocHome page
E. A. Hoffman and D. Chon
Computed Tomography Studies of Lung Ventilation and Perfusion
Proceedings of the ATS, December 1, 2005; 2(6): 492 - 498.
[Abstract] [Full Text] [PDF]


Home page
Proc Am Thorac SocHome page
G. Musch and J. G. Venegas
Positron Emission Tomography Imaging of Regional Pulmonary Perfusion and Ventilation
Proceedings of the ATS, December 1, 2005; 2(6): 522 - 527.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
N. T. Tgavalekos, M. Tawhai, R. S. Harris, G. Mush, M. Vidal-Melo, J. G. Venegas, and K. R. Lutchen
Identifying airways responsible for heterogeneous ventilation and mechanical dysfunction in asthma: an image functional modeling approach
J Appl Physiol, December 1, 2005; 99(6): 2388 - 2397.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
S. M. Walther, M. J. Johansson, T. Flatebo, A. Nicolaysen, and G. Nicolaysen
Marked differences between prone and supine sheep in effect of PEEP on perfusion distribution in zone II lung
J Appl Physiol, September 1, 2005; 99(3): 909 - 914.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
T. Richter, G. Bellani, R. S. Harris, M. F. V. Melo, T. Winkler, J. G. Venegas, and G. Musch
Effect of Prone Position on Regional Shunt, Aeration, and Perfusion in Experimental Acute Lung Injury
Am. J. Respir. Crit. Care Med., August 15, 2005; 172(4): 480 - 487.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
M. F. V. Melo, R. S. Harris, J. D. H. Layfield, and J. G. Venegas
Topographic Basis of Bimodal Ventilation-Perfusion Distributions during Bronchoconstriction in Sheep
Am. J. Respir. Crit. Care Med., April 1, 2005; 171(7): 714 - 721.
[Abstract] [Full Text] [PDF]


Home page
ChestHome page
F.-C. Lin, Y.-C. Chen, H.-I. Chang, and S.-C. Chang
Effect of Body Position on Gas Exchange in Patients With Idiopathic Pulmonary Alveolar Proteinosis: No Benefit of Prone Positioning
Chest, March 1, 2005; 127(3): 1058 - 1064.
[Abstract] [Full Text] [PDF]


Home page
Proc Am Thorac SocHome page
M. Dolovich and R. Labiris
Imaging Drug Delivery and Drug Responses in the Lung
Proceedings of the ATS, December 1, 2004; 1(4): 329 - 337.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
G. Peces-Barba, M. J. Rodriguez-Nieto, S. Verbanck, M. Paiva, and N. Gonzalez-Mangado
Lower pulmonary diffusing capacity in the prone vs. supine posture
J Appl Physiol, May 1, 2004; 96(5): 1937 - 1942.
[Abstract] [Full Text] [PDF]


Home page
RadiologyHome page
C. Fink, M. Puderbach, M. Bock, K.-P. Lodemann, I. Zuna, A. Schmahl, S. Delorme, and H.-U. Kauczor
Regional Lung Perfusion: Assessment with Partially Parallel Three-dimensional MR Imaging
Radiology, April 1, 2004; 231(1): 175 - 184.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
J. Petersson, A. Sanchez-Crespo, M. Rohdin, S. Montmerle, S. Nyren, H. Jacobsson, S. A. Larsson, S. G. E. Lindahl, D. Linnarsson, R. W. Glenny, et al.
Physiological evaluation of a new quantitative SPECT method measuring regional ventilation and perfusion
J Appl Physiol, March 1, 2004; 96(3): 1127 - 1136.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Cell Mol. Bio.Home page
D. P. Schuster, A. Kovacs, J. Garbow, and D. Piwnica-Worms
Recent Advances in Imaging the Lungs of Intact Small Animals
Am. J. Respir. Cell Mol. Biol., February 1, 2004; 30(2): 129 - 138.
[Abstract] [Full Text] [PDF]


Home page
J. Physiol.Home page
J. Gisolf, R. Wilders, R. V. Immink, J. J. van Lieshout, and J. M. Karemaker
Tidal volume, cardiac output and functional residual capacity determine end-tidal CO2 transient during standing up in humans
J. Physiol., January 15, 2004; 554(2): 579 - 590.
[Abstract] [Full Text] [PDF]


Home page
JNMHome page
M. F. Vidal Melo, D. Layfield, R. S. Harris, K. O'Neill, G. Musch, T. Richter, T. Winkler, A. J. Fischman, and J. G. Venegas
Quantification of Regional Ventilation-Perfusion Ratios with PET
J. Nucl. Med., December 1, 2003; 44(12): 1982 - 1991.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (41)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Musch, G.
Right arrow Articles by Venegas, J. G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Musch, G.
Right arrow Articles by Venegas, J. G.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online