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


     


J Appl Physiol 82: 1163-1176, 1997;
8750-7587/97 $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 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 Google Scholar
Google Scholar
Right arrow Articles by Treppo, S.
Right arrow Articles by Venegas, J. G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Treppo, S.
Right arrow Articles by Venegas, J. G.

Journal of Applied Physiology
Vol. 82, No. 4, pp. 1163-1176, April 1997
GAS EXCHANGE, MECHANICS, AND AIRWAYS

Contributions of pulmonary perfusion and ventilation to heterogeneity in VA/Q measured by PET

Steven Treppo, Srboljub M. Mijailovich, and José G. Venegas

Department of Anesthesia, Massachusetts General Hospital, and Harvard Medical School, Boston, Massachusetts 02114

ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
ACKNOWLEDGEMENTS
FOOTNOTES
REFERENCES


ABSTRACT

Treppo, Steven, Srboljub M. Mijailovich, and José G. Venegas. Contributions of pulmonary perfusion and ventilation to heterogeneity in VA/Q measured by PET. J. Appl. Physiol. 82(4): 1163-1176, 1997. To estimate the contributions of the heterogeneity in regional perfusion (Q) and alveolar ventilation (VA) to that of ventilation-perfusion ratio (VA/Q), we have refined positron emission tomography (PET) techniques to image local distributions of Q and VA per unit of gas volume content (sQ and sVA, respectively) and VA/Q in dogs. sVA was assessed in two ways: 1) the washout of 13NN tracer after equilibration by rebreathing (sVAi), and 2) the ratio of an apneic image after a bolus intravenous infusion of 13NN-saline solution to an image collected during a steady-state intravenous infusion of the same solution (sVAp). sVAp was systematically higher than sVAi in all animals, and there was a high spatial correlation between sQ and sVAp in both body positions (mean correlation was 0.69 prone and 0.81 supine) suggesting that ventilation to well-perfused units was higher than to those poorly perfused. In the prone position, the spatial distributions of sQ, sVAp, and VA/Q were fairly uniform with no significant gravitational gradients; however, in the supine position, these variables were significantly more heterogeneous, mostly because of significant gravitational gradients (15, 5.5, and -10%/cm, respectively) accounting for 73, 33, and 66% of the corresponding coefficient of variation (CV)2 values. We conclude that, in the prone position, gravitational forces in blood and lung tissues are largely balanced out by dorsoventral differences in lung structure. In the supine position, effects of gravity and structure become additive, resulting in substantial gravitational gradients in sQ and sVAp, with the higher heterogeneity in VA/Q caused by a gravitational gradient in sQ, only partially compensated by that in sVA.

positron emission tomography; body position; gas exchange; regional ventilation-perfusion ratio; dog; pulmonary heterogeneity; functional imaging


INTRODUCTION

WILSON AND BECK (30) have recently outlined a theoretical approach to assess the contributions of heterogeneities in alveolar ventilation (VA) and regional perfusion (Q) to the heterogeneity of the ventilation-perfusion ratio (VA/Q). Such an approach was illustrated with experimental data extrapolated from a compilation of reports by different investigators using different methodologies, leaving the quantitative validity of some of their conclusions in question. A number of recent studies have measured and characterized the effect of body position on the spatial heterogeneity of VA (2, 15, 27) and Q (3, 11, 12), but the individual contributions of these variables to the heterogeneity of VA/Q cannot be reliably assessed unless all variables are measured in the same individual. We have supplemented the positron imaging technique to measure VA/Q described by Rhodes and co-workers (21, 22) with independent measurements of Q and VA per unit of gas content (sQ and sVA, respectively) (17) using 13NN as the tracer gas. With the use of these techniques, we have imaged the distributions of sQ, sVA, and VA/Q and analyzed their spatial correlation and heterogeneity, including the contribution of gravitational gradients, in prone and supine dogs. The imaging data were also analyzed to yield and characterize distributions of VA/Q comparable to those generated by the multiple inert-gas-elimination technique (MIGET). This study provides the experimental data required by the method of Wilson and Beck (30) to assess the contributions of heterogeneities in VA and Q to the heterogeneity of VA/Q.


METHODS

The experimental protocol for these animal experiments was approved by the Massachusetts General Hospital Committee on Animal Care.

Experimental Setup

The experimental apparatus included a single-ring positron emission tomography (PET) camera, PCR-1, a mechanical ventilator, a rebreathing circuit, and an infusion system (Fig. 1). The PET camera, described in detail elsewhere (9), is a high-sensitivity stationary camera that is able to trigger the beginning of image collection to occur in synchrony with a signal from a mechanical ventilator. A closed-loop breathing circuit, including a CO2 absorber and supplemental oxygen, allowed ventilation with 13NN-labeled gas or with unlabeled gas while a solenoid-controlled valving system allowed rapid switching between these respiratory gas sources. 13NN-labeled gas, produced by a cyclotron, was introduced into the rebreathing circuit for inhaled tracer studies or forced into solution with previously degassed saline and temporarily stored in a chamber before intravenous infusion to the animal. In the latter case, a total of 200 ml of 13NN-labeled saline were produced, with specific activity ranging from 0.1 to 0.2 mCi/ml. The infusion system included a peristaltic pump and a remotely controlled solenoid valving system that allowed flushing of the tubing between the storing chamber and infusion catheter with 13NN-labeled saline. 13NN-labeled saline could be infused as a rapid bolus at a rate of 54 ml/min, or delivered by a peristaltic pump at a constant flow rate of 16 ml/min.
Fig. 1. Schematic representation of experimental apparatus. PET, positron emission tomography; ACg, specific activity of labeled gas in rebreathing circuit; ACv, specific activity of systemic blood; ACpa, specific activity of pulmonary artery blood; ACi, specific activity of intravenous infusate.
[View Larger Version of this Image (13K GIF file)]

Animal Preparation

Ten mongrel dogs weighing 14.4 ± 2.9 (SD) kg (range 8.5-19 kg) were intubated and mechanically ventilated under general anesthesia (pentobarbital sodium 30 mg/kg). The ventilator (Harvard Apparatus, Millis, MA) was set at a breathing frequency of 12 breaths/min, and the inspiratory time was set 30% of the breathing period. Tidal volume (VT) was set to maintain normocapnic arterial blood gases (mean VT = 22 ± 4.2 ml/kg, PCO2 = 35.5 ± 6 Torr). Femoral artery and vein were cannulated for blood sampling and periodic administration of anesthetic. A Swan-Ganz catheter (model 93A-131H-7F, Edwards Laboratory, Santa Ana, CA) was inserted in the right external jugular vein and advanced into the pulmonary artery, leaving its proximal port to drain into the superior vena cava. The proximal port of the catheter was used to deliver the 13NN-labeled saline solution, and the distal port was used for monitoring of pulmonary arterial pressure (Ppa) and sampling of pulmonary arterial blood. Before each imaging run, the lungs were inflated and sustained at a pressure between 20 and 30 cmH2O to minimize the occurrence of microatelectasis and the loss of compliance.

Imaging Scans

Image collection during mechanical ventilation was gated by using an electrical signal from the mechanical ventilator at the start of inspiration. The gating scheme consisted of a collection of two consecutive images of equal duration during the breathing cycle. Because inspiratory time was set at 30% of the total breathing period, the first image included all of inspiration and the initial part of exhalation, whereas the second image included the latter part of exhalation, in which expiratory flow was small and the lungs remained almost stationary at a lung volume close to resting functional residual capacity (FRC). Five animals were studied in the prone position and five in the supine position. A transverse cross section of the lungs intersecting the apex of the heart was selected for imaging. Such a position was determined after equilibrating the lungs with inhaled 13NN-labeled gas and advancing the lungs into the field of view of the camera until the highest count rate was recorded by the camera. This position corresponded to the slice of the lungs having the greatest cross-sectional area and provided consistency among animals. After positioning, the following scan sequences were performed: inhaled tracer, bolus infusion, constant infusion, blood pool scan, and transmission and uniformity scans. These scan sequences are described in detail below, and the average counts per voxel of the resulting PET images are given in Table 1.

Table 1. Average counts per voxel of collected PET images


PET Image (Tracer) Mean Counts/Voxel Total Counts of Image

Inhaled tracer equilibration (13NN) 1,542 3,000,000
Inhaled tracer washout (13NN) 125 250,000
Bolus infusion (13NN-saline) 755 1,510,000
Constant infusion (13NN-saline) 1,513 3,020,000
Blood pool scan (11CO) 1,300 2,600,000
Transmission scan 10,500,000
Uniformity scan 62,700,000

PET, positron emission tomography.

Inhaled tracer sequence. The lungs were ventilated with 13NN-labeled gas from the closed breathing circuit with 100% oxygen. Once equilibration of the tracer was achieved between lungs and breathing circuit (~5 min), a gated imaging sequence was collected for a total of 240 s. Because of the low solubility of nitrogen in body fluids and tissues, 13NN remained mostly confined to the air spaces within the lung. This equilibration image was, therefore, proportional to gas volume content. During the time of imaging, a 1-ml gas sample was obtained from the rebreathing circuit to assess its specific activity (ACg). At the end of collection of the equilibration image, the inspiratory gas was switched to room air, and three gated images of 60-s duration were collected consecutively as the tracer was washed out from the lungs.

Bolus-infusion sequence. Starting with a tracer-free lung, mechanical ventilation was interrupted at end exhalation, and immediately a bolus of 13NN-labeled saline solution was infused at a flow rate of 54 ml/min into the superior vena cava through the proximal port of the Swan-Ganz catheter. Infusion time ranged between 4 and 15 s (median = 10 s) depending on the specific activity of the infusate to produce images with a consistent number of counts per voxel. Simultaneously with the beginning of infusion, a collection of five consecutive images was initiated. The first three images had a scanning time of 5 s and the last two of 30 s. When imaging was completed, ventilation was resumed, and the tracer was washed out with no further imaging after the apneic period. A sample of the intravenous infusate was collected to assess its specific activity (ACi).

Constant rate-infusion sequence. During steady-state mechanical ventilation, infusion of 13NN-labeled saline was initiated into the proximal port of the Swan-Ganz catheter at a constant flow rate of 16 ml/min. Simultaneously with the start of infusion, a collection sequence of four 120-s duration gated images was initiated. Once a steady-state activity was reached (generally after the 1st 2 min; Fig. 2), samples of infusate solution and pulmonary arterial and systemic blood were obtained to assess their respective specific activities (ACi, ACpa, and ACv).
Fig. 2. Average time course of normalized count rates during constant-infusion protocol including all runs reported. Data are normalized by average count rate of last 3 images of the series. Note that after 1st 2 min of infusion regional counts reach a steady state.
[View Larger Version of this Image (13K GIF file)]

Blood pool scan. To estimate the contribution of counts originating from the pulmonary arterial blood per voxel during the constant-infusion (CI) imaging sequence, an additional blood pool scan was conducted. This was done by labeling the red blood cells with a temporary inhalation of 11C-labeled CO until the steady-state activity of the blood reached an adequate level. Two sequential images of 20-min duration were then collected. Blood samples were obtained at the beginning and the end of the imaging sequence to assess their respective specific activity needed to normalize the image and to correct for ventilatory losses of the CO tracer over the imaging time.

Transmission and uniformity scans. To correct for gamma ray energy attenuation caused by the supporting structures and body tissues of the animal, a tubular ring, concentric to the PET camera's field, was filled with 18F-labeled water, and a gated transmission scan was collected during breathing. At the end of this scan, the animal and supporting structures were removed from the camera's field, and a final uniformity field scan was conducted.

Data Analysis

Image processing. PET images were initially corrected for camera sensitivity and for tissue attenuation. Image reconstruction was then performed with a convolution back-projection algorithm by using a Hanning filter yielding an effective spatial resolution of 10 mm determined from the width at one-half height of a point source image. This degradation of resolution length (from 4.5 mm of the camera) was needed to attenuate random noise to levels <5% of the measured coefficient of variation (CV)2 in most images processed. Resulting images consisted of an interpolated matrix of 159 × 159 voxels of 0.2 × 0.2 cm × 5 mm, or 57% of the resolution length. These reconstructed images of local counts per voxel were then processed following the methodology described in detail in the accompanying paper (17) and briefly discussed below, to yield functional images with voxel values in physical units (Table 2).

Table 2. Mean of local values and heterogeneity parameters of measured variables averaged for prone and supine animals


Variable Prone (n = 5) Supine (n = 5) P Value

 VA/Q
  Mean 1.380 ± 0.290  1.060 ± 0.270  0.23
  CV2 0.019 ± 0.003  0.122 ± 0.036* 0.023
  CV2 residual 0.017 ± 0.002  0.031 ± 0.004* 0.0095
  Gradient, %/cm  -0.63 ± 1.021   -10.26 ± 1.498dagger * 0.0005
  Grad, %CV2 10.0 ± 4.0  66.0 ± 7.3* 0.0002
 Q
  Mean, ml · s-1 · cm-3 0.031 ± 0.006  0.061 ± 0.006* 0.0041
  CV2 0.066 ± 0.011  0.170 ± 0.023* 0.0034
  CV2 residual 0.059 ± 0.010  0.103 ± 0.013* 0.0161
  Gradient, %/cm 0.70 ± 1.649  9.55 ± 0.732dagger * 0.0017
  Grad, %CV2 8.8 ± 4.0  38.3 ± 6.2* 0.0027
sQ
  Mean, s-1 0.056 ± 0.011  0.103 ± 0.012* 0.0096
  CV2 0.034 ± 0.007  0.221 ± 0.037* 0.0031
  CV2 residual 0.027 ± 0.007  0.057 ± 0.009* 0.015
  Gradient, %/cm 0.91 ± 1.171  15.1 ± 0.49dagger * 0.0004
  Grad, %CV2 17.3 ± 11.0  72.8 ± 3.2* 0.0028
VA
  Mean, ml/cm3 0.570 ± 0.042  0.640 ± 0.029  0.09
  CV2 0.044 ± 0.004  0.080 ± 0.013* 0.022
  CV2 residual 0.044 ± 0.004  0.060 ± 0.008  0.064
  Gradient, %/cm  -0.37 ± 0.480   -5.11 ± 0.67dagger * 0.0003
  Grad, % CV2 1.2 ± 0.4  23.7 ± 5.0* 0.050
sVAi
  Mean, s-1 0.042 ± 0.005  0.042 ± 0.005  0.48
  CV2 0.034 ± 0.013  0.052 ± 0.010  0.16
  CV2 residual 0.024 ± 0.008  0.045 ± 0.009* 0.019
  Gradient, %/cm 3.78 ± 1.213dagger 0.83 ± 1.095  0.054
  Grad, %CV2 19.7 ± 4.1  10.2 ± 3.0* 0.05
sVAp
  Mean, s-1 0.070 ± 0.014  0.086 ± 0.016  0.23
  CV2 0.026 ± 0.007  0.063 ± 0.010* 0.0091
  CV2 residual 0.024 ± 0.007  0.041 ± 0.008  0.068
  Gradient, %/cm 0.343 ± 0.879  5.517 ± 1.076dagger * 0.0031
  Grad, %CV2 8.4 ± 5.7  33.2 ± 8.3* 0.022

Values are means ± SE. CV, coefficient of variation (SD/mean); CV residual, heterogeneity that remains after removal of any gradient in dorsoventral direction by linear regression; Gradient, magnitude of this gradient in % of average value per cm ( dagger significantly different from zero, P < 0.05); Grad (%CV2), regression coefficient of linear regression fit and, to the 1st order, contribution of gradient to total heterogeneity. VA/Q, ventilation-perfusion ratio; Q, perfusion; sQ, specific perfusion; VA, gas content; sVAp, specific ventilation from perfusion tracer; sVAi, specific ventilation from inhaled tracer. * significantly different from prone, P < 0.05).

SELECTION OF VOXELS FOR ANALYSIS AND VA/Q IMAGE. A steady-state CI image was created by adding, on a voxel-by-voxel basis, the last three images of the protocol sequence. An initial mask was created by thresholding the CI image to exclude areas outside of the lung field. A threshold of 30% was used initially and then refined in increments until no areas outside the lung field were included in the mask (mean threshold used was 33 ± 6%). A second mask was then created by thresholding the blood pool (Vb) image to define the heart and largest vessels. (This was also done in an iterative process where mean threshold used was 55 ± 5%.) This second mask was subtracted from the first one to exclude heart and vessels from further analysis. The masked CI image was then corrected for the contribution of pulmonary arterial blood activity with the algorithm described by Rhodes and co-workers (21, 22), where the volume of radiolabeled arterial blood is assumed to be 40% of the total pulmonary blood volume assessed from the 11CO scan. Thus a bloodcorrection image was subtracted voxel-by-voxel from the CI image. A temporary image was formed from the ratio of uncorrected to blood-corrected CI images, and the mask was further refined to exclude from the analysis additional areas with very high correction values. These areas, covering ~20 voxels out of 2,000 voxels on average, were typically located in proximity to the heart and large blood vessels and corresponded to overcorrected voxels by partial volume effects from the Vb scan. These final masks were then applied to all functional images analyzed. Cardiac output (QT) was calculated by using a mass balance from the infusion flow rate and the 13NN specific activities of the infusate and the pulmonary arterial blood.

ALVEOLAR GAS CONTENT. An alveolar gas content per voxel image (VA) was obtained by decay-correcting the voxel values of the equilibrated inhaled tracer scan, normalized by the specific activity of a gas sample to create an image in units of milliliters of gas content per cubic centimeter of voxel.

Q. Because of diffusion to neighboring regions and/or readsorption of the tracer into capillary blood, the voxel tracer concentration may have changed during the apneic period following the bolus intravenous infusion of 13NN-labeled saline. An estimation of, and correction for, these tracer kinetics effects was conducted on a voxel-by-voxel basis by assessing the differences in tracer content between the last two 30-s images of the sequence and then extrapolating the activity level expected at the time of the tracer's arrival to the alveoli (17). The tracer-kinetics-corrected image was normalized by the ratio of total infused activity to QI to yield an image of local Q in units of milliliter per minute of blood flow per cubic centimeter of voxel volume. Finally, the Q image was divided in a voxel-by-voxel manner by the VA image to yield an image of sQ (in units of s-1).

SVA. We derived images of regional sVA by using two independent methods. One method directly assessed the kinetics of inhaled NN2 tracer during a washout maneuver following an equilibration scan (sVAi), as described in the accompanying paper (17). A second method indirectly assessed ventilation (sVAp) as the ratio of the local concentration of the NN2 tracer, infused during apnea (distribution of Q), divided by the local concentration of the tracer during CI of the tracer in saline solution in steady-state breathing, after subtraction of activity from the pulmonary arterial blood (distribution of Q/sVA) (17). The resulting sVAp image represented exclusively sVA of perfused units, since unperfused units would not receive tracer during either of the two imaging protocols.

Assessment of spatial heterogeneity. The spatial heterogeneity of the functional images was assessed from the CV of the voxel data within the lung field defined as the SD normalized by the mean value of the data
CV = <FR><NU>SD</NU><DE>mean</DE></FR>
To assess the true regional heterogeneity of the images, it was necessary to estimate and correct for the contributions to heterogeneity caused by statistical noise and imaging artifacts including imperfect corrections for tissue attenuation, nonuniformity of the camera sensitivity, or by edge blurring created by finite spatial resolution of PET. Theoretical models to characterize these contributions exist (1, 4, 8), but the actual effect is dependent on camera design, reconstruction algorithm, and dimensions of the object. We used a recently described experimental method to evaluate such artifacts by using a lung phantom (25).

Briefly, for functional images derived from individual PET images (such as Q or VA), the contribution of noise to the CV2 was calculated as the random noise (g2) caused by the finite counts from the image. Thus, the noise-corrected CV (CVcr) was
CV<SUB>cr</SUB> = <RAD><RCD>CV<SUP>2</SUP> − <IT>g</IT><SUP>2</SUP></RCD></RAD>

The value of g2 was found to be inversely proportional to the average number of counts per voxel (<OVL><IT>n</IT></OVL>) of the original PET image or
<IT>g</IT><SUP>2</SUP> = <FR><NU><IT>k</IT><SUB>g</SUB></NU><DE><OVL><IT>n</IT></OVL></DE></FR>
where kg = 2.76 was a constant experimentally derived for our camera from phantom images (25) reconstructed with the same filter and convolution back-projection algorithm and the same masking thresholding as the ones used in this experiment.

For functional images derived from the ratio of two PET images (such as VA/Q, sVA, or sQ), the CV of the ratio image [CV2(x/y)] was corrected by the sum of random-noise contributions to the original images, yielding
CV<SUB>cr</SUB>(<IT>x</IT>/<IT>y</IT>) = <RAD><RCD>CV<SUP>2</SUP>(<IT>x</IT>/<IT>y</IT>) − ( <IT>g</IT><SUP>2</SUP><SUB><IT>x</IT></SUB> + <IT>g</IT><SUP>2</SUP><SUB><IT>y</IT></SUB>)</RCD></RAD>

VERTICAL (GRAVITATIONAL) HETEROGENEITY. To assess the contribution of any significant vertical gradient to the total heterogeneity of an image, the regional data were fitted by a linear-regression model using the vertical coordinate as independent variable. The CV2cr was divided into the part due to the residuals from the regression (nongravitational heterogeneity) and the part accounted for by the regression (gravitational contribution).

SPATIAL CORRELATIONS. The spatial correlation coefficient (RS) between two regional variables was calculated from the corresponding log-transformed functional images on a voxel-by-voxel basis. RS was calculated for the following pairs of variables: VA vs. Q and vs. sQ; and sQ vs. sVAi and vs. sVAp.

In those cases where the pairs of functional images were originally derived by using a common image, i.e., VA vs. sQ, Rs was corrected to eliminate the pseudocorrelation caused by imaging noise in the common image VA.1

FRACTIONAL DISTRIBUTIONS. Mean-normalized distribution histograms for VA, Q, sQ, sVA, and VA/Q and their corresponding log-transformed versions were generated for each of these functional images. Regional data were grouped by either fraction of total voxels (lung fraction), VA, VA and/or Q fractions. These distributions were then characterized by evaluating the corresponding Pearson coefficient of skewness (Skx) and coefficient of kurtosis (kappa ).

BIVARIATE DISTRIBUTIONS. Mean-normalized bivariate-distribution histograms for log-transformed mean-normalized sQ vs. sVA were also generated, in which sVA was calculated from either inhaled or perfused tracer. These bivariate distributions were plotted as three-dimensional surfaces with the z-axis corresponding to the fraction of total voxels having the corresponding relative values of sVA and sQ. Distributions were then averaged within each group of animals.

Shunt Fraction

Because of the low solubility of nitrogen in blood and tissues, an index of overall lung venous admixture can be calculated from the fraction of tracer recirculating back into the lungs during the steady-state period of CI. Such a recirculation fraction (FR) can be estimated from the ratio between the peripheral CV and Cpa simultaneously measured during the steady-state part of the CI protocol
<IT>F</IT><SUB>R</SUB> = <FR><NU>C<SUB>V</SUB></NU><DE>C<SUB>pa</SUB></DE></FR>

Statistical Analysis

Comparisons between supine and prone positions were made by using single-tailed Student's t-test for independent samples. Comparisons between average and CV values of sVAi and sVAp were made using multivariate analysis of variance with body position and method as factors. Statistical significance was taken at P < 0.05 level.


RESULTS

Mean ± SE values for supine and prone positions of the average voxel value within the imaged section for Q, sQ, VA, sVAp, sVAi, and VA/Q are presented in Table 2. Of these parameters, only Q and sQ were significantly greater in the supine compared with the prone position.

Measurement of Heterogeneity

The contributions from different factors to heterogeneity, such as noise or vertical gradients, are additive only when expressed in terms of CV2; thus the findings of this study are presented in Table 2 and Fig. 3 as CV2. Because the definition of CV (SD/mean) gives a more intuitive impression of the degree of heterogeneity, in the text we present the results in terms of the CV.
Fig. 3. Individual values of coefficient of variation (CV)2 for pulmonary perfusion per unit of gas volume (sQ), ventilation-perfusion ratio (VA/Q), and specific alveolar ventilation measured from kinetics of perfused tracer (sVAp) images in 5 supine (open circle ) and 5 prone animals (square ). Note that CV2 values of supine dogs are higher that those in prone dogs and that in all animals CV2 of VA/Q is lower than that of sQ.
[View Larger Version of this Image (18K GIF file)]

VA/Q. Although there were substantial differences in spatial heterogeneity among the different dogs (Fig. 3), the heterogeneity of VA/Q was significantly lower in the prone animals (CV = 0.14) compared with the supine animals (CV = 0.34) (Table 2). This difference in heterogeneity was partially accounted for by the presence of a systematic vertical gradient in the supine position that accounted for 66% of the CV2, whereby VA/Q decreased by 10%/cm distance in the direction of gravity. In contrast, the vertical gradient in the prone position was not significantly different from zero (-0.63%/cm). After we removed the vertical gradient from the images by linear regression, the differences in residual heterogeneity between supine and prone positions were still statistically significant but of a much lesser magnitude [CVr = 0.13 for prone and 0.17 for supine positions (Table 2)].

Q. QT, measured from the specific tracer activities of the pulmonary artery and the saline infusate (ACpa and ACi) was 1.01 ± 0.180 l/min for prone and 1.27 ± 0.340 l/min for supine animals. Mean values of average regional Q (<OVL><A><AC>Q</AC><AC>˙</AC></A></OVL>) were 0.031 ml · s-1 · cm-3 for prone and 0.061 ml · s-1 · cm-3 for supine dogs (Table 2). Regional distributions of Q either normalized by voxel volume Q or by alveloar gas volume sQ, were more heterogeneous in the supine than in the prone position. CVQ and CVsQ were 0.41 and 0.46 in supine and 0.25 and 0.18 in prone position, respectively (Table 2). The greater heterogeneity of Q and sQ in the supine position was due, in part, to consistent gravitational gradients, whereby the respective variables increased by 9.5 and 15.1%/cm distance in the direction of gravity. These gradients contributed to 38 and 73% of the total CV Q2 and CVsQ2 , respectively. In contrast, there were no consistent vertical gradients in the prone position.

When the regional distributions of sQ and Q for each position were compared, sQ was significantly less heterogeneous than Q in the prone position while the contrary was true in the supine position (Table 2). Differences in the width of the corresponding fractional distribution are consistent with this finding (Fig. 4, B and C, right) by showing that the average distribution of sQ in the prone position is wider than the corresponding distribution of Q, with the opposite happening in the supine position.
Fig. 4. Representative images of alveolar gas volume (VA; A); pulmonary perfusion (Q; B), and (sQ; C) in prone (left) and supine (middle) positions. Height of surfaces over the x-y plane represents relative value of corresponding variable. Plots on right are fractional distributions of corresponding variables averaged for all dogs studied. Overbars designate average values.
[View Larger Version of this Image (62K GIF file)]

VA. As illustrated by a wider fractional distribution of VA (Fig. 4A, right), the spatial heterogeneity of local gas content was greater in the supine position compared with the prone position (CV was 0.28 for supine and 0.21 for prone; Table 2). The higher heterogeneity of the supine position was partially accounted for by a systematic vertical gradient, whereby gas content decreased by 5.1%/cm distance in the direction of gravity. This gradient contributed to 23.7% of the total CV2. In contrast, the prone position had smaller and not significant vertical gradients (-0.370%/cm), without significant contributions to the total CV2 (1.2%).

sVA. sVA was directly assessed with the 13NN tracer delivered by inhalation (sVAi), and, indirectly, from the ratio of images with the 13NN delivered by intravenous infusion (sVAp). Two-way analysis of variance with repeated measures on the values of average s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC> (s<OVL><A><AC>V</AC><AC>˙</AC></A></OVL><SC>a</SC>) for both methods, including the fixed effect of body position, showed a significant effect of method but not of body position. Student's t-tests confirmed no statistical differences for s<OVL><A><AC>V</AC><AC>˙</AC></A></OVL><SC>a</SC><SUB>i</SUB> between supine (0.042 s-1; Table 2) and prone (0.042 s-1) positions and for s<OVL><A><AC>V</AC><AC>˙</AC></A></OVL><SC>a</SC><SUB>p</SUB> between supine (0.086 s-1) and prone (0.070 s-1) positions but demonstrated s<OVL><A><AC>V</AC><AC>˙</AC></A></OVL><SC>a</SC><SUB>p</SUB> to be significantly greater than s<OVL><A><AC>V</AC><AC>˙</AC></A></OVL><SC>a</SC><SUB>i</SUB> for each position. There was, however, a significant correlation between the individual values of s<OVL><A><AC>V</AC><AC>˙</AC></A></OVL><SC>a</SC><SUB>i</SUB> and s<OVL><A><AC>V</AC><AC>˙</AC></A></OVL><SC>a</SC><SUB>p</SUB>, and a linear fit between the two independent estimates of sVA had a slope of 1.66 and R2 = 0.59 (Fig. 5).
Fig. 5. Overall lung specific ventilation measured from washout of inhaled tracer (<OVL><SC>s</SC><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL><SUB>i</SUB>) plotted against corresponding <OVL>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL><SUB>p</SUB> for each animal studied. <OVL>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL><SUB>i</SUB> was consistently lower than <OVL>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL><SUB>p</SUB>.
[View Larger Version of this Image (18K GIF file)]

Local values of sVAp on a voxel-by-voxel basis were poorly correlated with, and consistently higher than, those of sVAi (Fig. 6). CV of sVAi of prone (0.18 ± 0.09) was not significantly different from that of supine position (0.22 ± 0.1), whereas the CV of sVAp of prone animals (0.16) was significantly lower than that of supine dogs (0.25). Variations in the width of the corresponding fractional distribution histograms of sVAp and sVAi illustrate these findings (Fig. 7).
Fig. 6. Representative distributions on a voxel-by-voxel basis of sVAp vs. sVAi for a representative dog in supine (A) and one in prone (B) position.
[View Larger Version of this Image (17K GIF file)]


Fig. 7. Representative images of sVAi (A) and sVAp (B) in prone (left) and supine positions (middle). Height of surfaces over the x-y plane represents relative value of corresponding variable. Plots on right are fractional distributions of corresponding variables averaged for all dogs studied.
[View Larger Version of this Image (42K GIF file)]

A consistent difference between the spatial distributions of sVAp and sVAi occurred in the supine position where sVAp presented a vertical gradient of 5.5%/cm length in the direction of gravity accounting for 33.2% of the total CV2, whereas sVAi had no significant gradient (Table 2).

Spatial Correlations

Q was found to have a high and significant spatial correlation with VA in the prone position (RS = 0.74 ± 0.05) but none in the supine position (RS = 0.23 ± 0.01) (Table 3). sQ was not spatially correlated with sVAi in either position (RS = -0.036 ± 0.13 for prone, and RS = 0.11 ± 0.23 for supine positions). In contrast, sQ was highly correlated to sVAp in both prone and supine positions (RS = 0.69 and 0.81, respectively). This correlation is illustrated in the bivariate distributions (Fig. 8 and Fig. 9), where most voxels contain combinations of log(s<A><AC>Q</AC><AC>˙</AC></A>/s<OVL><A><AC>Q</AC><AC>˙</AC></A></OVL>), where s<OVL><A><AC>Q</AC><AC>˙</AC></A></OVL> is average regional sQ, and log(sVAp/s<OVL><A><AC>V</AC><AC>˙</AC></A></OVL><SC>a</SC><SUB>p</SUB>) that fall along a 45° projection (constant VA/Q).

Table 3. Coefficients of spatial correlation between functional images of supine and prone dogs


Variables Spatial Correlations
Prone Supine

 Q vs. VA 0.743 ± 0.047dagger 0.229 ± 0.0107*dagger
sVAi vs. sQ  -0.036 ± 0.135  0.105 ± 0.226 
sVAp vs. sQ 0.693 ± 0.043dagger 0.809 ± 0.08dagger

Values are means ± SE; n = 5 dogs/group. * Significantly different from prone, P < 0.05; dagger significantly different from 0, P < 0.05.


Fig. 8. Bivariate distribution of log(s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><SUB>p</SUB>/<OVL>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL><SUB>p</SUB>) vs. log(sQ/s<OVL><A><AC>Q</AC><AC>˙</AC></A></OVL>) averaged over all dogs studied in prone position. Height of surface over the x-y plane in surface plot (bottom left) represents average fraction of voxels containing respective combination of relative sVAp and sQ. A contour plot (top right) illustrates correlation between these variables (direction of constant VA/Q isopleths is 45°).
[View Larger Version of this Image (20K GIF file)]


Fig. 9. Bivariate distribution of log(s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><SUB>p</SUB>/<OVL>s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL><SUB>p</SUB>) vs. log(s<A><AC>Q</AC><AC>˙</AC></A>/s<OVL><A><AC>Q</AC><AC>˙</AC></A></OVL> ) averaged over all dogs studied in supine position. Height of surface over the x-y plane in surface plot (bottom left) represents average fraction of voxels containing respective combination of relative sVAp and sQ. A contour plot (top right) illustrates correlation between these variables (direction of constant VA/Q isopleths is 45°).
[View Larger Version of this Image (22K GIF file)]

Fractional Distributions

Mean-normalized fractional distributions of functional images originating from the ratio of two PET images (sVAp, sQ, and VA/Q) were all skewed to the right, as shown by the mean Skx significantly greater than zero at both body positions (Table 4). Logarithmic transformation of the data results in unskewed fractional distributions with mean Skx values not different from zero. The log-transformed distributions were also mesokurtic, i.e., they had coefficient of kurtosis kappa  that does not appreciably deviate from 3. This suggests that the regional distributions closely resemble the shape of a log-normal distribution.

Table 4. Average values of skewness and kurtosis coefficients of mean-normalized distributions of lung function variables, and corresponding log-transformed variables for prone and supine dogs


Distribution Variable Grouping Variable Skewness Coefficient
Kurtosis Coefficient
Prone (n = 5) Supine (n = 5) Prone (n = 5) Supine (n = 5)

ln(V<SC>a</SC>/<OVL><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL>)
VL  -1.163 ± 0.083   -1.134 ± 0.073  4.546 ± 0.288  4.491 ± 0.200 
V<SC>a</SC>/<OVL><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL>
VL  -0.452 ± 0.054   -0.303 ± 0.070  2.805 ± 0.075  2.752 ± 0.081 
ln(s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><SUB>i</SUB>/s<OVL><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL><SUB>i</SUB>)
VL 0.136 ± 0.033   -0.137 ± 0.243  3.336 ± 0.116  5.447 ± 0.326*
s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><SUB>i</SUB>/s<OVL><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL><SUB>i</SUB>
VL 0.810 ± 0.080  1.202 ± 0.365  4.481 ± 0.311  11.227 ± 1.639 
ln(s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><SUB>p</SUB>/s<OVL><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL><SUB>p</SUB>)
VL  -0.309 ± 0.060   -0.533 ± 0.064  3.684 ± 0.187  3.711 ± 0.155 
s<A><AC>V</AC><AC>˙</AC></A><SC>a</SC><SUB>p</SUB>/s<OVL><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL><SUB>p</SUB>
VL 0.217 ± 0.016  0.267 ± 0.065  3.246 ± 0.055  3.173 ± 0.067 
ln(<A><AC>Q</AC><AC>˙</AC></A>/<OVL><A><AC>Q</AC><AC>˙</AC></A></OVL>)
VL  -0.743 ± 0.067   -0.538 ± 0.032  3.391 ± 0.102  2.861 ± 0.083 
<A><AC>Q</AC><AC>˙</AC></A>/<OVL><A><AC>Q</AC><AC>˙</AC></A></OVL>
VL  -0.085 ± 0.076  0.324 ± 0.024* 2.575 ± 0.026  2.381 ± 0.044 
ln(s<A><AC>Q</AC><AC>˙</AC></A>/<OVL>s<A><AC>Q</AC><AC>˙</AC></A></OVL>
VL  -0.051 ± 0.093   -0.278 ± 0.058  3.744 ± 0.283  2.597 ± 0.151 
s<A><AC>Q</AC><AC>˙</AC></A>/s<OVL><A><AC>Q</AC><AC>˙</AC></A></OVL>
VL 0.504 ± 0.077  0.644 ± 0.073  3.984 ± 0.359  2.963 ± 0.199 
VA 0.179 ± 0.070   -0.202 ± 0.123  4.035 ± 0.188  2.869 ± 0.143*
ln[(<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>/<A><AC>Q</AC><AC>˙</AC></A>)/(<OVL><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL>/<OVL><A><AC>Q</AC><AC>˙</AC></A></OVL>)]
 Q 0.110 ± 0.055  0.127 ± 0.106  4.208 ± 0.229  2.769 ± 0.123*
 VA 0.228 ± 0.069   -0.209 ± 0.129  4.099 ± 0.176  3.009 ± 0.163*
VA 0.545 ± 0.039  0.483 ± 0.134  3.504 ± 0.107  3.237 ± 0.239 
(<A><AC>V</AC><AC>˙</AC></A><SC>a</SC>/<A><AC>Q</AC><AC>˙</AC></A>)/(<OVL><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></OVL> /<OVL><A><AC>Q</AC><AC>˙</AC></A></OVL>)
 Q 0.534 ± 0.025  0.829 ± 0.122  3.572 ± 0.104  3.984 ± 0.324 
 VA 0.589 ± 0.031  0.503 ± 0.135  3.555 ± 0.091  3.373 ± 0.240

Values are means ± SE. Overbars designate average regional values. VL, fraction of lung volume. * Significantly different from prone, P < 0.05.

In contrast, fractional distributions of functional images originating from single PET images (VA, Q, and VA) were not skewed (Skx not different from zero) and were mesokurtic while logarithmic transformation of the voxel data made the distributions significantly skewed to the left (Skx < 0).

Shunt Fraction

Reflecting the higher degree of venous admixture expected from a less uniform VA/Q distribution, the FR of the infused tracer recirculating back into the lungs during the CI protocol was significantly higher in the supine position (FR = 0.018 ± 0.005) compared with the prone position (FR = 0.007 ± 0.005).


DISCUSSION

The most significant findings of this study were as follows. 1) The regional distribution of VA/Q was more heterogeneous in the supine position compared with the prone position. The higher heterogeneity supine was due to a significant vertical gradient that contributed more than one-half of the total CV2. 2) Differences in the heterogeneities between supine and prone in both Q and VA contributed to the differences in VA/Q heterogeneity. 3) In both body positions, there was a high spatial correlation between sQ and sVAp while there was no correlation between sQ and sVAi. Regional Q and VA were positively correlated in the prone position but not in the supine position.

Methodology

Most issues related to the PET imaging methodology have either been discussed by the original proponents of the CI technique (7, 21, 22) or have been discussed in detail in the accompanying paper (17). A major methodological departure from the original method to assess VA/Q was in the estimation of VA, where we imaged the lungs after equilibration with inhaled 13NN gas instead of deriving VA from the transmission scan as described by Rhodes et al. (22). This modification improved the quality and signal-to-noise ratio of the VA/Q images because it largely canceled out systematic imaging artifacts and because transmission scans with an equal number of events have greater inherent noise than emission scans. As pointed out by Brudin et al. (7), this improvement in image quality is only realized when imaging normal lungs, since delayed equilibration in areas of very low ventilation introduces errors in VA with the inhaled 13NN technique. Given the small size of our animals, it was crucial to use this modification to obtain the highest possible spatial resolution of our instrument, compatible with an appropriate signal-to-noise ratio.

Aside from the contributions of random noise to image heterogeneity, there are systematic distortions introduced by PET imaging. For example, the filtering involved in the convolution-back-projection algorithm not only affects the spatial resolution of the camera but also results in smearing of the lung edges that artifactually increases the estimated heterogeneity of an image that includes them. Voxel size of our images was 0.2 cm × 0.2 cm × 5 mm with the average amount of 37 ml of lung studied for the 10 animals (from an average of 2,000 voxels/animal). Maximal resolution of the instrument was 4.5 mm × 4.5 mm × 5 mm. The resolution length, defined as the width at half maximum of a point source, was increased to 1 cm after filtering. Other effects, such as imperfect uniformity calibration and transmission corrections, also contribute to an increase in the CV2 of images derived from single PET scans. Fortunately, these systematic distortions cancel out in images obtained from the voxel-by-voxel ratio of two independently acquired PET images such as in sVA, sQ, and VA/Q (25). The contribution of systematic imaging artifacts to the CV2 of single PET images such as VA or Q was also estimated by imaging a uniform lung-like phantom and reconstructing, masking, and thresholding the resulting images in the same way as our images. The CV2 measured from those images was 0.015, which corresponds to 34 and 25% of the CV2 in the prone position and to 18 and 12% in the supine position, measured for VA and Q, respectively.

Heterogeneity of VA/Q

The spatial heterogeneity in VA/Q was found to be higher in supine (CV = 0.34) than in prone animals (CV = 0.14). In the supine position, 66.0% of the CV2 was attributed to a significant vertical gradient (10.3%/cm). The prone position, in contrast, had no systematic vertical gradients in VA/Q, although there was substantial interanimal variability with gradients ranging from -3.3 to +2.0 %/cm. The residual heterogeneity of the supine position, after subtraction of the vertical effect, was still significantly greater than that of the prone position. This might be attributed to the inadequacy of the linear-regression model in describing nonlinear vertical gradients. The CI technique to assess the distribution of VA/Q was originally described by Rhodes et al. (21). Although the authors did not study the effect of body position, they reported a CV = 0.21 for supine, spontaneously breathing normal humans. More recently, the same group of investigators (5) reported a very modest vertical gradient in VA/Q, whereby VA/Q decreased, on average, by 2%/cm in the ventral-to-dorsal direction, with this gradient explaining only 20% of the CV2. Remarkably, the single patient studied in the prone position showed a more substantial gradient in the direction of gravity than the group of supine patients. The lower CV2 and gradients in humans compared with our dogs could be due to a difference in distribution of ventilation between spontaneously breathing and mechanically ventilated subjects (20). Also, partial volume effects, exaggerated by a lack of respiratory gating and poorer spatial resolution of their PET instrument (1.7 cm) compared with ours (1 cm), must have also accounted for lower CV in comparison with our study.

Using MIGET, Beck and co-workers (3) reported values of lnSD for the main VA/Q distribution peak of 0.45 and 0.35 for the supine and prone positions, respectively. For narrow distributions, the lnSD approximates the CV (30), and the results from MIGET appear somewhat greater than our direct measurement of CV for VA/Q. Intrinsic differences between PET and MIGET need to be considered before discussing these results. MIGET distributions are understood to reflect the overall heterogeneity of VA/Q of the whole lung and at all physiologically relevant length scales. In normal experimental conditions, however, the capability of MIGET to resolve narrow distributions of VA/Q is limited to distributions with SDlog 0.2 to 0.3 (19, 29). In contrast, the heterogeneity in VA/Q measured by PET in this study is based on actual topographical distribution but is limited to detect heterogeneities with length scales greater than the spatial resolution of our PET imaging method (1 cm) and to sampling a single transverse cross section of the lung.

Thus part of the higher CV seen by MIGET, compared with PET, could possibly be attributed to the limited sampling and spatial resolution of our study. Although for the supine position the CV values measured by MIGET (0.45) and PET (0.34) do not appear to be too different, for the prone position the CV for MIGET (0.35) was more than double that measured by PET (0.14). Thus, at first glance, one could attribute the greater CV recovered from MIGET to a limitation of the technique to resolve narrow distributions