Prone positioning redistributes gravitational stress in the lung in normal conditions and in simulations of oedema.

NEW FINDINGS: What is the central question of this study? How does the interaction between posture and gravity affect the stresses on the lung, particularly in highly inflated gravitationally non-dependent regions, which are potentially vulnerable to increased mechanical stress and injury? What is the main finding and its importance? Changes in stress attributable to gravity are not well characterized between postures. Using a new metric of gravitational stress, we show that regions of the lung near maximal inflation have the greatest gravitational stresses while supine, but not while prone. In simulations of increased lung weight consistent with severe pulmonary oedema, the prone lung has lower gravitational stress in vulnerable, non-dependent regions, potentially protecting them from overinflation and injury. ABSTRACT: Prone posture changes the gravitational vector, and potentially the stress induced by tissue deformation, because a larger lung volume is gravitationally dependent when supine, but non-dependent when prone. To evaluate this, 10 normal subjects (six male and four female; age, means ± SD = 27 ± 6 years; height, 171 ± 9 cm; weight, 69 ± 13 kg; forced expiratory volume in the first second/forced expiratory volume as a percentage of predicted, 93 ± 6%) were imaged at functional residual capacity, supine and prone, using magnetic resonance imaging, to quantify regional lung density. We defined regional gravitational stress as the cumulative weight, per unit area, of the column of lung tissue below each point. Gravitational stress was compared between regions of differing inflation to evaluate differences between highly stretched, and thus potentially vulnerable, regions and less stretched lung. Using reference density values for normal lungs at total lung capacity (0.10 ± 0.03 g/ml), regions were classified as highly inflated (density < 0.13 g/ml, i.e., close to total lung capacity), intermediate (0.13 ≤ density < 0.16 g/ml) or normally inflated (density ≥ 0.16 g/ml). Gravitational stress differed between inflation categories while supine (-1.6 ± 0.3 cmH2 O highly inflated; -1.4 ± 0.3 cmH2 O intermediate; -1.1 ± 0.1 cmH2 O normally inflated; P = 0.05) but not while prone (-1.4 ± 0.2 cmH2 O highly inflated; -1.3 ± 0.2 cmH2 O intermediate; -1.3 ± 0.1 cmH2 O normally inflated; P = 0.39), and increased more with height from dependent lung while supine (-0.24 ± 0.02 cmH2 O/cm supine; -0.18 ± 0.04 cmH2 O/cm prone; P = 0.05). In simulated severe pulmonary oedema, the gradient in gravitational stress increased in both postures (all P < 0.0001), was greater in the supine posture than when prone (-0.57 ± 0.21 cmH2 O/cm supine; -0.34 ± 0.16 cmH2 O/cm prone; P = 0.0004) and was similar to the gradient calculated from supine computed tomography images in a patient with acute respiratory distress syndrome (-0.51 cmH2 O/cm). The non-dependent lung has greater gravitational stress while supine and might be protected while prone, particularly in the presence of oedema.

Oxygen-Sensitive MRI: A Predictive Imaging Biomarker for Tumor Radiation Response?

PURPOSE: To develop a noninvasive prognostic imaging biomarker related to hypoxia to predict SABR tumor control. METHODS AND MATERIALS: A total of 145 subcutaneous syngeneic Dunning prostate R3327-AT1 rat tumors were focally irradiated once using cone beam computed tomography guidance on a small animal irradiator at 225 kV. Various doses in the range of 0 to 100 Gy were administered, while rats breathed air or oxygen, and tumor control was assessed up to 200 days. Oxygen-sensitive magnetic resonance imaging (MRI) (T1-weighted, ΔR1, ΔR2*) was applied to 79 of these tumors at 4.7 T to assess response to an oxygen gas breathing challenge on the day before irradiation as a probe of tumor hypoxia. RESULTS: Increasing radiation dose in the range of 0 to 90 Gy enhanced tumor control of air-breathing rats with a TCD50 estimated at 59.6 ± 1.5 Gy. Control was significantly improved at some doses when rats breathed oxygen during irradiation (eg, 40 Gy; P < .05), and overall there was a modest left shift in the control curve: TCD50(oxygen) = 53.1 ± 3.1 Gy (P < .05 vs air). Oxygen-sensitive MRI showed variable response to oxygen gas breathing challenge; the magnitude of T1-weighted signal response (%ΔSI) allowed stratification of tumors in terms of local control at 40 Gy. Tumors showing %ΔSI >0.922 with O2-gas breathing challenge showed significantly better control at 40 Gy during irradiation while breathing oxygen (75% vs 0%, P < .01). In addition, increased radiation dose (50 Gy) substantially overcame resistance, with 50% control for poorly oxygenated tumors. Stratification of dose-response curves based on %ΔSI >0.922 revealed different survival curves, with TCD50 = 36.2 ± 3.2 Gy for tumors responsive to oxygen gas breathing challenge; this was significantly less than the 54.7 ± 2.4 Gy for unresponsive tumors (P < .005), irrespective of the gas inhaled during tumor irradiation. CONCLUSIONS: Oxygen-sensitive MRI allowed stratification of tumors in terms of local control at 40 Gy, indicating its use as a potential predictive imaging biomarker. Increasing dose to 50 Gy overcame radiation resistance attributable to hypoxia in 50% of tumors.

Regional pulmonary perfusion patterns in humans are not significantly altered by inspiratory hypercapnia.

Pulmonary vascular tone is known to be sensitive to both local alveolar Po2 and Pco2. Although the effects of hypoxia are well studied, the hypercapnic response is relatively less understood. We assessed changes in regional pulmonary blood flow in humans in response to hypercapnia using previously developed MRI techniques. Dynamic measures of blood flow were made in a single slice of the right lung of seven healthy volunteers following a block-stimulus paradigm (baseline, challenge, recovery), with CO2 added to inspired gas during the challenge block to effect a 7-Torr increase in end-tidal CO2. Effects of hypercapnia on blood flow were evaluated based on changes in spatiotemporal variability (fluctuation dispersion, FD) and in regional perfusion patterns in comparison to hypoxic effects previously studied. Hypercapnia increased FD 2.5% from baseline (relative to control), which was not statistically significant (P = 0.07). Regional perfusion patterns were not significantly changed as a result of increased FICO2 (P = 0.90). Reanalysis of previously collected data using a similar protocol but with the physiological challenge replaced by decreased FIO2 (FIO2 = 0.125) showed marked flow redistribution (P = 0.01) with the suggestion of a gravitational pattern, demonstrating hypoxia has the ability to affect regional change with a global stimulus. Taken together, these data indicate that hypercapnia of this magnitude does not lead to appreciable changes in the distribution of pulmonary perfusion, and that this may represent an interesting distinction between the hypoxic and hypercapnic regulatory response.NEW & NOTEWORTHY Although it is well known that the pulmonary circulation responds to local alveolar hypoxia, and that this mechanism may facilitate ventilation-perfusion matching, the relative role of CO2 is not well appreciated. This study demonstrates that an inspiratory hypercapnic stimulus is significantly less effective at inducing changes in pulmonary perfusion patterns than inspiratory hypoxia, suggesting that in these circumstances hypercapnia is not sufficient to induce substantial integrated feedback control of ventilation-perfusion mismatch across the lung.

Oxygen-sensitive MRI assessment of tumor response to hypoxic gas breathing challenge.

Oxygen-sensitive MRI has been extensively used to investigate tumor oxygenation based on the response (R2 * and/or R1 ) to a gas breathing challenge. Most studies have reported response to hyperoxic gas indicating potential biomarkers of hypoxia. Few studies have examined hypoxic gas breathing and we have now evaluated acute dynamic changes in rat breast tumors. Rats bearing syngeneic subcutaneous (n = 15) or orthotopic (n = 7) 13762NF breast tumors were exposed to a 16% O2 gas breathing challenge and monitored using blood oxygen level dependent (BOLD) R2 * and tissue oxygen level dependent (TOLD) T1 -weighted measurements at 4.7 T. As a control, we used a traditional hyperoxic gas breathing challenge with 100% O2 on a subset of the subcutaneous tumor bearing rats (n = 6). Tumor subregions identified as responsive on the basis of R2 * dynamics coincided with the viable tumor area as judged by subsequent H&E staining. As expected, R2 * decreased and T1 -weighted signal increased in response to 100% O2 breathing challenge. Meanwhile, 16% O2 breathing elicited an increase in R2 *, but divergent response (increase or decrease) in T1 -weighted signal. The T1 -weighted signal increase may signify a dominating BOLD effect triggered by 16% O2 in the relatively more hypoxic tumors, whereby the influence of increased paramagnetic deoxyhemoglobin outweighs decreased pO2 . The results emphasize the importance of combined BOLD and TOLD measurements for the correct interpretation of tumor oxygenation properties.

Comparison of quantitative multiple-breath specific ventilation imaging using colocalized 2D oxygen-enhanced MRI and hyperpolarized 3He MRI.

Two magnetic resonance specific ventilation imaging (SVI) techniques, namely, oxygen-enhanced proton (OE-1H) and hyperpolarized 3He (HP-3He), were compared in eight healthy supine subjects [age 32 (6) yr]. An in-house radio frequency coil array for 1H configured with the 3He transmit-receive coil in situ enabled acquisition of SVI data from two nuclei from the same slice without repositioning the subjects. After 3 × 3 voxel downsampling to account for spatial registration errors between the two SV images, the voxel-by-voxel correlation coefficient of two SV maps ranged from 0.11 to 0.63 [0.46 mean (0.17 SD); P < 0.05]. Several indexes were analyzed and compared from the tidal volume-matched SV maps: the mean of SV log-normal distribution (SVmean), the standard deviation of the distribution as a measure of SV heterogeneity (SVwidth), and the gravitational gradient (SVslope). There were no significant differences in SVmean [OE-1H: 0.28 (0.08) and HP-3He: 0.32 (0.14)], SVwidths [OE-1H: 0.28 (0.08) and HP-3He: 0.27 (0.10)], and SVslopes [OE-1H: -0.016 (0.006) cm-1 and HP-3He: -0.013 (0.007) cm-1]. Despite the statistical similarities of the population averages, Bland-Altman analysis demonstrated large individual intertechnique variability. SDs of differences in these indexes were 42% (SVmean), 46% (SVwidths), and 62% (SVslopes) of their corresponding overall mean values. The present study showed that two independent, spatially coregistered, SVI techniques presented a moderate positive voxel-by-voxel correlation. Population averages of SVmean, SVwidth, and SVslope were in close agreement. However, the lack of agreement when the data sets were analyzed individually might indicate some fundamental mechanistic differences between the techniques. NEW & NOTEWORTHY To the best of our knowledge, this is the first cross-comparison of two different specific ventilation (SV) MRI techniques in the human lung (i.e., oxygen-enhanced proton and hyperpolarized 3He). The present study showed that two types of spatially coregistered SV images presented a modest positive correlation. The two techniques also yielded similar population averages of SV indexes such as log-normal mean, SV heterogeneity, and the gravitational slope, albeit with some intersubject variability.

MR-CBCT image-guided system for radiotherapy of orthotopic rat prostate tumors.

Multi-modality image-guided radiotherapy is the standard of care in contemporary cancer management; however, it is not common in preclinical settings due to both hardware and software limitations. Soft tissue lesions, such as orthotopic prostate tumors, are difficult to identify using cone beam computed tomography (CBCT) imaging alone. In this study, we characterized a research magnetic resonance (MR) scanner for preclinical studies and created a protocol for combined MR-CBCT image-guided small animal radiotherapy. Two in-house dual-modality, MR and CBCT compatible, phantoms were designed and manufactured using 3D printing technology. The phantoms were used for quality assurance tests and to facilitate end-to-end testing for combined preclinical MR and CBCT based treatment planning. MR and CBCT images of the phantoms were acquired utilizing a Varian 4.7 T scanner and XRad-225Cx irradiator, respectively. The geometry distortion was assessed by comparing MR images to phantom blueprints and CBCT. The corrected MR scans were co-registered with CBCT and subsequently used for treatment planning. The fidelity of 3D printed phantoms compared to the blueprint design yielded favorable agreement as verified with the CBCT measurements. The geometric distortion, which varied between -5% and 11% throughout the scanning volume, was substantially reduced to within 0.4% after correction. The distortion free MR images were co-registered with the corresponding CBCT images and imported into a commercial treatment planning software SmART Plan. The planning target volume (PTV) was on average 19% smaller when contoured on the corrected MR-CBCT images relative to raw images without distortion correction. An MR-CBCT based preclinical workflow was successfully designed and implemented for small animal radiotherapy. Combined MR-CBCT image-guided radiotherapy for preclinical research potentially delivers enhanced relevance to human radiotherapy for various disease sites. This novel protocol is wide-ranging and not limited to the orthotopic prostate tumor study presented in the study.

Measurement of the distribution of ventilation-perfusion ratios in the human lung with proton MRI: comparison with the multiple inert-gas elimination technique.

We have developed a novel functional proton magnetic resonance imaging (MRI) technique to measure regional ventilation-perfusion (V̇A/Q̇) ratio in the lung. We conducted a comparison study of this technique in healthy subjects (n = 7, age = 42 ± 16 yr, Forced expiratory volume in 1 s = 94% predicted), by comparing data measured using MRI to that obtained from the multiple inert gas elimination technique (MIGET). Regional ventilation measured in a sagittal lung slice using Specific Ventilation Imaging was combined with proton density measured using a fast gradient-echo sequence to calculate regional alveolar ventilation, registered with perfusion images acquired using arterial spin labeling, and divided on a voxel-by-voxel basis to obtain regional V̇A/Q̇ ratio. LogSDV̇ and LogSDQ̇, measures of heterogeneity derived from the standard deviation (log scale) of the ventilation and perfusion vs. V̇A/Q̇ ratio histograms respectively, were calculated. On a separate day, subjects underwent study with MIGET and LogSDV̇ and LogSDQ̇ were calculated from MIGET data using the 50-compartment model. MIGET LogSDV̇ and LogSDQ̇ were normal in all subjects. LogSDQ̇ was highly correlated between MRI and MIGET (R = 0.89, P = 0.007); the intercept was not significantly different from zero (-0.062, P = 0.65) and the slope did not significantly differ from identity (1.29, P = 0.34). MIGET and MRI measures of LogSDV̇ were well correlated (R = 0.83, P = 0.02); the intercept differed from zero (0.20, P = 0.04) and the slope deviated from the line of identity (0.52, P = 0.01). We conclude that in normal subjects, there is a reasonable agreement between MIGET measures of heterogeneity and those from proton MRI measured in a single slice of lung.NEW & NOTEWORTHY We report a comparison of a new proton MRI technique to measure regional V̇A/Q̇ ratio against the multiple inert gas elimination technique (MIGET). The study reports good relationships between measures of heterogeneity derived from MIGET and those derived from MRI. Although currently limited to a single slice acquisition, these data suggest that single sagittal slice measures of V̇A/Q̇ ratio provide an adequate means to assess heterogeneity in the normal lung.

Characterizing spatiotemporal information loss in sparse-sampling-based dynamic MRI for monitoring respiration-induced tumor motion in radiotherapy.

PURPOSE: Sparse-sampling and reconstruction techniques represent an attractive strategy to achieve faster image acquisition speeds, while maintaining adequate spatial resolution and signal-to-noise ratio in rapid magnetic resonance imaging (MRI). The authors investigate the use of one such sequence, broad-use linear acquisition speed-up technique (k-t BLAST) in monitoring tumor motion for thoracic and abdominal radiotherapy and examine the potential trade-off between increased sparsification (to increase imaging speed) and the potential loss of "true" information due to greater reliance on a priori information. METHODS: Lung tumor motion trajectories in the superior-inferior direction, previously recorded from ten lung cancer patients, were replayed using a motion phantom module driven by an MRI-compatible motion platform. Eppendorf test tubes filled with water which serve as fiducial markers were placed in the phantom. The modeled rigid and deformable motions were collected in a coronal image slice using balanced fast field echo in conjunction with k-t BLAST. Root mean square (RMS) error was used as a metric of spatial accuracy as measured trajectories were compared to input data. The loss of spatial information was characterized for progressively increasing acceleration factor from 1 to 16; the resultant sampling frequency was increased approximately from 2.5 to 19 Hz when the principal direction of the motion was set along frequency encoding direction. In addition to the phantom study, respiration-induced tumor motions were captured from two patients (kidney tumor and lung tumor) at 13 Hz over 49 s to demonstrate the impact of high speed motion monitoring over multiple breathing cycles. For each subject, the authors compared the tumor centroid trajectory as well as the deformable motion during free breathing. RESULTS: In the rigid and deformable phantom studies, the RMS error of target tracking at the acquisition speed of 19 Hz was approximately 0.3-0.4 mm, which was smaller than the reconstructed pixel resolution of 0.67 mm. In the patient study, the dynamic 2D MRI enabled the monitoring of cycle-to-cycle respiratory variability present in the tumor position. It was seen that the range of centroid motion as well as the area covered due to target motion during each individual respiratory cycle was underestimated compared to the entire motion range observed over multiple breathing cycles. CONCLUSIONS: The authors' initial results demonstrate that sparse-sampling- and reconstruction-based dynamic MRI can be used to achieve adequate image acquisition speeds without significant information loss for the task of radiotherapy guidance. Such monitoring can yield spatial and temporal information superior to conventional offline and online motion capture methods used in thoracic and abdominal radiotherapy.

The effect of lung deformation on the spatial distribution of pulmonary blood flow.

KEY POINTS: Pulmonary perfusion measurement using magnetic resonance imaging combined with deformable image registration enabled us to quantify the change in the spatial distribution of pulmonary perfusion at different lung volumes. The current study elucidated the effects of tidal volume lung inflation [functional residual capacity (FRC) + 500 ml and FRC + 1 litre] on the change in pulmonary perfusion distribution. Changes in hydrostatic pressure distribution as well as transmural pressure distribution due to the change in lung height with tidal volume inflation are probably bigger contributors to the redistribution of pulmonary perfusion than the changes in pulmonary vasculature resistance caused by lung tissue stretch. ABSTRACT: Tidal volume lung inflation results in structural changes in the pulmonary circulation, potentially affecting pulmonary perfusion. We hypothesized that perfusion is recruited to regions receiving the greatest deformation from a tidal breath, thus ensuring ventilation-perfusion matching. Density-normalized perfusion (DNP) magnetic resonance imaging data were obtained in healthy subjects (n = 7) in the right lung at functional residual capacity (FRC), FRC+500 ml, and FRC+1.0 l. Using deformable image registration, the displacement of a sagittal lung slice acquired at FRC to the larger volumes was calculated. Registered DNP images were normalized by the mean to estimate perfusion redistribution (nDNP). Data were evaluated across gravitational regions (dependent, middle, non-dependent) and by lobes (upper, RUL; middle, RML; lower, RLL). Lung inflation did not alter mean DNP within the slice (P = 0.10). The greatest expansion was seen in the dependent region (P < 0.0001: dependent vs non-dependent, P < 0.0001: dependent vs middle) and RLL (P = 0.0015: RLL vs RUL, P < 0.0001: RLL vs RML). Neither nDNP recruitment to RLL [+500 ml = -0.047(0.145), +1 litre = 0.018(0.096)] nor to dependent lung [+500 ml = -0.058(0.126), +1 litre = -0.023(0.106)] were found. Instead, redistribution was seen in decreased nDNP in the non-dependent [+500 ml = -0.075(0.152), +1 litre = -0.137(0.167)) and increased nDNP in the gravitational middle lung [+500 ml = 0.098(0.058), +1 litre = 0.093(0.081)] (P = 0.01). However, there was no significant lobar redistribution (P < 0.89). Contrary to our hypothesis, based on the comparison between gravitational and lobar perfusion data, perfusion was not redistributed to the regions of the most inflation. This suggests that either changes in hydrostatic pressure or transmural pressure distribution in the gravitational direction are implicated in the redistribution of perfusion away from the non-dependent lung.

A continuous surface reconstruction method on point cloud captured from a 3D surface photogrammetry system.

PURPOSE: To accurately and efficiently reconstruct a continuous surface from noisy point clouds captured by a surface photogrammetry system (VisionRT). METHODS: The authors have developed a level-set based surface reconstruction method on point clouds captured by a surface photogrammetry system (VisionRT). The proposed method reconstructs an implicit and continuous representation of the underlying patient surface by optimizing a regularized fitting energy, offering extra robustness to noise and missing measurements. By contrast to explicit/discrete meshing-type schemes, their continuous representation is particularly advantageous for subsequent surface registration and motion tracking by eliminating the need for maintaining explicit point correspondences as in discrete models. The authors solve the proposed method with an efficient narrowband evolving scheme. The authors evaluated the proposed method on both phantom and human subject data with two sets of complementary experiments. In the first set of experiment, the authors generated a series of surfaces each with different black patches placed on one chest phantom. The resulting VisionRT measurements from the patched area had different degree of noise and missing levels, since VisionRT has difficulties in detecting dark surfaces. The authors applied the proposed method to point clouds acquired under these different configurations, and quantitatively evaluated reconstructed surfaces by comparing against a high-quality reference surface with respect to root mean squared error (RMSE). In the second set of experiment, the authors applied their method to 100 clinical point clouds acquired from one human subject. In the absence of ground-truth, the authors qualitatively validated reconstructed surfaces by comparing the local geometry, specifically mean curvature distributions, against that of the surface extracted from a high-quality CT obtained from the same patient. RESULTS: On phantom point clouds, their method achieved submillimeter reconstruction RMSE under different configurations, demonstrating quantitatively the faith of the proposed method in preserving local structural properties of the underlying surface in the presence of noise and missing measurements, and its robustness toward variations of such characteristics. On point clouds from the human subject, the proposed method successfully reconstructed all patient surfaces, filling regions where raw point coordinate readings were missing. Within two comparable regions of interest in the chest area, similar mean curvature distributions were acquired from both their reconstructed surface and CT surface, with mean and standard deviation of (µrecon=-2.7×10(-3) mm(-1), σrecon=7.0×10(-3) mm(-1)) and (µCT=-2.5×10(-3) mm(-1), σCT=5.3×10(-3) mm(-1)), respectively. The agreement of local geometry properties between the reconstructed surfaces and the CT surface demonstrated the ability of the proposed method in faithfully representing the underlying patient surface. CONCLUSIONS: The authors have integrated and developed an accurate level-set based continuous surface reconstruction method on point clouds acquired by a 3D surface photogrammetry system. The proposed method has generated a continuous representation of the underlying phantom and patient surfaces with good robustness against noise and missing measurements. It serves as an important first step for further development of motion tracking methods during radiotherapy.

The heterogeneity of regional specific ventilation is unchanged following heavy exercise in athletes.

Heavy exercise increases ventilation-perfusion mismatch and decreases pulmonary gas exchange efficiency. Previous work using magnetic resonance imaging (MRI) arterial spin labeling in athletes has shown that, after 45 min of heavy exercise, the spatial heterogeneity of pulmonary blood flow was increased in recovery. We hypothesized that the heterogeneity of regional specific ventilation (SV, the local tidal volume over functional residual capacity ratio) would also be increased following sustained exercise, consistent with the previously documented changes in blood flow heterogeneity. Trained subjects (n = 6, maximal O2 consumption = 61 ± 7 ml·kg(-1)·min(-1)) cycled 45 min at their individually determined ventilatory threshold. Oxygen-enhanced MRI was used to quantify SV in a sagittal slice of the right lung in supine posture pre- (preexercise) and 15- and 60-min postexercise. Arterial spin labeling was used to measure pulmonary blood flow in the same slice bracketing the SV measures. Heterogeneity of SV and blood flow were quantified by relative dispersion (RD = SD/mean). The alveolar-arterial oxygen difference was increased during exercise, 23.3 ± 5.3 Torr, compared with rest, 6.3 ± 3.7 Torr, indicating a gas exchange impairment during exercise. No significant change in RD of SV was seen after exercise: preexercise 0.78 ± 0.15, 15 min postexercise 0.81 ± 0.13, 60 min postexercise 0.78 ± 0.08 (P = 0.5). The RD of blood flow increased significantly postexercise: preexercise 1.00 ± 0.12, 15 min postexercise 1.15 ± 0.10, 45 min postexercise 1.10 ± 0.10, 60 min postexercise 1.19 ± 0.11, 90 min postexercise 1.11 ± 0.12 (P < 0.005). The lack of a significant change in RD of SV postexercise, despite an increase in the RD of blood flow, suggests that airways may be less susceptible to the effects of exercise than blood vessels.

Affine transformation registers small scale lung deformation.

To evaluate the nature of small scale lung deformation between multiple pulmonary magnetic resonance images, two different kinematic intensity based image registration techniques: affine and bicubic Hermite interpolation were tested. The affine method estimates uniformly distributed deformation metrics throughout the lung. The bicubic Hermite method allows the expression of heterogeneously distributed deformation metrics such as Lagrangian strain. A cardiac triggered inversion recovery technique was used to obtain 10 sequential images of pulmonary vessel structure in a sagittal plane in the right lung at FRC in 4 healthy subjects (Age: 28.5(6.2)). One image was used as the reference image, and the remaining images (target images) were warped onto the reference image using both image registration techniques. The normalized correlation between the reference and the transformed target images within the lung domain was used as a cost function for optimization, and the root mean square (RMS) of image intensity difference was used to evaluate the quality of the registration. Both image registration techniques significantly improved the RMS compared with non-registered target images (p= 0.04). The spatial mean (µE) and standard deviation (σ(E)) of Lagrangian strain were computed based on the spatial distribution of lung deformation approximated by the bicubic Hermite method, and were measured on the order of 10(-3) or less, which is virtually negligible. As a result, small scale lung deformation between FRC lung volumes is spatially uniform, and can be simply characterized by affine deformation even though the bicubic Hermite method is capable of expressing complicated spatial patterns of lung deformation.

Measuring lung water: ex vivo validation of multi-image gradient echo MRI.

PURPOSE: To validate a fast gradient echo sequence for rapid (9 s) quantitative imaging of lung water. MATERIALS AND METHODS: Eleven excised pig lungs were imaged with a fast GRE sequence in triplicate, in the sagittal plane at 2 levels of inflation pressure (5 and 15 cm H(2) O), an intervention that alters T(2) *, but not total lung water. Images were acquired alternating between two closely-spaced echoes and data were fit (voxel-by-voxel) to a single exponential to determine T(2) * and water content, and compared with gravimetric measurements of total water. RESULTS: T(2) * averaged 1.08 ± 0.02 ms at 5 cm H(2) O and 1.02 ± 0.02 ms at 15 cm H(2) O (P < 0.05). The measure was reliable (R(2) = 0.99), with an average mean error of 1.8%. There was a significant linear relationship between the two measures of water content: The regression equations for the relationship were y = 0.92x + 19 (R(2) = 0.94), and y = 1.04x + 4 (R(2) = 0.96), for 5 and 15 cm H(2) O inflation pressure respectively. Y-intercepts were not statistically different from zero (P = 0.86). CONCLUSION: The multi-echo GRE sequence is a reliable and valid technique to assess water content in the lung. This technique enables rapid assessment of lung water, which is advantageous for in vivo studies.

Magnetic resonance imaging quantification of pulmonary perfusion using calibrated arterial spin labeling.

UNLABELLED: This demonstrates a MR imaging method to measure the spatial distribution of pulmonary blood flow in healthy subjects during normoxia (inspired O(2), fraction (F(I)O(2)) = 0.21) hypoxia (F(I)O(2) = 0.125), and hyperoxia (F(I)O(2) = 1.00). In addition, the physiological responses of the subject are monitored in the MR scan environment. MR images were obtained on a 1.5 T GE MRI scanner during a breath hold from a sagittal slice in the right lung at functional residual capacity. An arterial spin labeling sequence (ASL-FAIRER) was used to measure the spatial distribution of pulmonary blood flow and a multi-echo fast gradient echo (mGRE) sequence was used to quantify the regional proton (i.e. H(2)O) density, allowing the quantification of density-normalized perfusion for each voxel (milliliters blood per minute per gram lung tissue). With a pneumatic switching valve and facemask equipped with a 2-way non-rebreathing valve, different oxygen concentrations were introduced to the subject in the MR scanner through the inspired gas tubing. A metabolic cart collected expiratory gas via expiratory tubing. Mixed expiratory O(2) and CO(2) concentrations, oxygen consumption, carbon dioxide production, respiratory exchange ratio, respiratory frequency and tidal volume were measured. Heart rate and oxygen saturation were monitored using pulse-oximetry. Data obtained from a normal subject showed that, as expected, heart rate was higher in hypoxia (60 bpm) than during normoxia (51) or hyperoxia (50) and the arterial oxygen saturation (SpO(2)) was reduced during hypoxia to 86%. Mean ventilation was 8.31 L/min BTPS during hypoxia, 7.04 L/min during normoxia, and 6.64 L/min during hyperoxia. Tidal volume was 0.76 L during hypoxia, 0.69 L during normoxia, and 0.67 L during hyperoxia. Representative quantified ASL data showed that the mean density normalized perfusion was 8.86 ml/min/g during hypoxia, 8.26 ml/min/g during normoxia and 8.46 ml/min/g during hyperoxia, respectively. In this subject, the relative dispersion, an index of global heterogeneity, was increased in hypoxia (1.07 during hypoxia, 0.85 during normoxia, and 0.87 during hyperoxia) while the fractal dimension (Ds), another index of heterogeneity reflecting vascular branching structure, was unchanged (1.24 during hypoxia, 1.26 during normoxia, and 1.26 during hyperoxia). Overview. This protocol will demonstrate the acquisition of data to measure the distribution of pulmonary perfusion noninvasively under conditions of normoxia, hypoxia, and hyperoxia using a magnetic resonance imaging technique known as arterial spin labeling (ASL). RATIONALE: Measurement of pulmonary blood flow and lung proton density using MR technique offers high spatial resolution images which can be quantified and the ability to perform repeated measurements under several different physiological conditions. In human studies, PET, SPECT, and CT are commonly used as the alternative techniques. However, these techniques involve exposure to ionizing radiation, and thus are not suitable for repeated measurements in human subjects.

Lung volume does not alter the distribution of pulmonary perfusion in dependent lung in supine humans.

There is a gravitational influence on pulmonary perfusion, including in the most dependent lung, where perfusion is reduced, termed Zone 4. Studies using xenon-133 show Zone 4 behaviour, present in the dependent 4 cm at total lung capacity (TLC), affects the dependent 11 cm at functional residual capacity (FRC) and almost all the lung at residual volume (RV). These differences were ascribed to increased resistance in extra-alveolar vessels at low lung volumes although other mechanisms have been proposed. To further evaluate the behaviour of perfusion in dependent lung using a technique that directly measures pulmonary perfusion and corrects for tissue distribution by measuring regional proton density, seven healthy subjects (age = 38 ± 6 years, FEV1 = 104 ± 7% predicted) underwent magnetic resonance imaging in supine posture. Data were acquired in the right lung during breath-holds at RV, FRC and TLC. Arterial spin labelling quantified regional pulmonary perfusion, which was normalized for regional proton density measured using a fast low-angle shot technique. The height of the onset of Zone 4 behaviour was not different between lung volumes (P = 0.23). There were no significant differences in perfusion (expressed as ml min⁻¹ g⁻¹) between lung volumes in the gravitationally intermediate (RV = 8.9 ± 3.1, FRC = 8.1 ± 2.9, TLC = 7.4 ± 3.6; P = 0.26) and dependent lung (RV = 6.6 ± 2.4, FRC = 6.1 ± 2.1, TLC = 6.4 ± 2.6; P = 0.51). However, at TLC perfusion was significantly lower in non-dependent lung than at FRC or RV (3.6 ± 3.3, 7.7 ± 1.5, 7.9 ± 2.0, respectively; P < 0.001). These data suggest that the mechanism of the reduction in perfusion in dependent lung is unlikely to be a result of lung volume related increases in resistance in extra-alveolar vessels. In supine posture, the gravitational influence on perfusion is remarkably similar over most of the lung, irrespective of lung volume.

Vertical distribution of specific ventilation in normal supine humans measured by oxygen-enhanced proton MRI.

Specific ventilation (SV) is the ratio of fresh gas entering a lung region divided by its end-expiratory volume. To quantify the vertical (gravitationally dependent) gradient of SV in eight healthy supine subjects, we implemented a novel proton magnetic resonance imaging (MRI) method. Oxygen is used as a contrast agent, which in solution changes the longitudinal relaxation time (T1) in lung tissue. Thus alterations in the MR signal resulting from the regional rise in O(2) concentration following a sudden change in inspired O(2) reflect SV-lung units with higher SV reach a new equilibrium faster than those with lower SV. We acquired T1-weighted inversion recovery images of a sagittal slice of the supine right lung with a 1.5-T MRI system. Images were voluntarily respiratory gated at functional residual capacity; 20 images were acquired with the subject breathing air and 20 breathing 100% O(2), and this cycle was repeated five times. Expired tidal volume was measured simultaneously. The SV maps presented an average spatial fractal dimension of 1.13 ± 0.03. There was a vertical gradient in SV of 0.029 ± 0.012 cm(-1), with SV being highest in the dependent lung. Dividing the lung vertically into thirds showed a statistically significant difference in SV, with SV of 0.42 ± 0.14 (mean ± SD), 0.29 ± 0.10, and 0.24 ± 0.08 in the dependent, intermediate, and nondependent regions, respectively (all differences, P < 0.05). This vertical gradient in SV is consistent with the known gravitationally induced deformation of the lung resulting in greater lung expansion in the dependent lung with inspiration. This SV imaging technique can be used to quantify regional SV in the lung with proton MRI.

Rapid intravenous infusion of 20 ml/kg saline does not impair resting pulmonary gas exchange in the healthy human lung.

Rapid infusion of intravenous saline, a model of pulmonary interstitial edema, alters the distribution of pulmonary perfusion, raises pulmonary capillary blood volume, and increases bronchial wall thickness in humans. We hypothesized that infusion would disrupt pulmonary gas exchange by increasing ventilation/perfusion ((.)VA/(.)Q) inequality as opposed to a diffusive impairment in O2 exchange. Seven males (26 +/- 3 yr; FEV1: 110 +/- 16% predicted.) performed spirometry and had (.)VA/(.)Q mismatch measured using the multiple inert gas elimination technique, before and after 20 ml/kg iv of normal saline delivered in approximately 30 min. Infusion increased thoracic fluid content from transthoracic impedance by 12% (P < 0.0001) and left FVC unchanged but reduced expiratory flows (FEF(25-75) falling from 5.1 +/- 0.4 to 4.2 +/- 0.4 l/s, P < 0.05). However, (.)VA/(.)Q mismatch as measured by the log standard deviation of the ventilation (LogSD(.)V) and perfusion (LogSD(.)Q) distributions remained unchanged; LogSD(.)V: 0.40 +/- 0.03 pre, 0.38 +/- 0.04 post, NS; LogSD(.)Q: 0.38 +/- 0.03 pre, 0.37 +/- 0.03 post, NS. There was no significant change in arterial PO2 (99 +/- 2 pre, 99 +/- 3 mmHg post, NS) but arterial PCO2 was decreased (38.7 +/- 0.6 pre, 36.8 +/- 1.2 mmHg post, P < 0.05). Thus, infusion compressed small airways and caused a mild degree of hyperventilation. There was no evidence for a diffusive limitation to O2 exchange, with the measured-predicted alveolar-arterial oxygen partial pressure difference being unaltered by infusion at FIO2 = 0.125 (4.3 +/- 1.0 pre, 5.2 +/- 1.0 post, NS). After infusion, the fraction of perfusion going to areas with (.)VA/(.)Q < 1 was increased when a subject breathed a hyperoxic gas mixture [0.72 +/- 0.06 (FIO2 = 0.21), 0.80 +/- 0.06 (FIO2 = 0.30), P < 0.05] with similar effects on ventilation in the face of unchanged (.)VA and (.)Q. These results suggest active control of blood flow to regions of decreased ventilation during air breathing, thus minimizing the gas exchange consequences.

Quantitative MRI measurement of lung density must account for the change in T(2) (*) with lung inflation.

PURPOSE: To evaluate lung water density at three different levels of lung inflation in normal lungs using a fast gradient echo sequence developed for rapid imaging. MATERIALS AND METHODS: Ten healthy volunteers were imaged with a fast gradient echo sequence that collects 12 images alternating between two closely spaced echoes in a single 9-s breathhold. Data were fit to a single exponential to determine lung water density and T(2) (*). Data were evaluated in a single imaging slice at total lung capacity (TLC), functional residual capacity (FRC), and residual volume (RV). Analysis of variance for repeated measures was used to statistically evaluate changes in T(2) (*) and lung water density across lung volumes, imaging plane, and spatial locations in the lung. RESULTS: In normal subjects (n = 10), T(2) (*) (and [lung density/water density]) was 1.2 +/- 0.1 msec (0.10 +/- 0.02), 1.8 +/- 0.2 ms (0.25 +/- 0.04), and 2.0 +/- 0.2 msec (0.27 +/- 0.03) at TLC, FRC, and RV, respectively. Results also show that there is a considerable intersubject variability in the values of T(2) (*). CONCLUSION: Data show that T(2) (*) in the lung is very short, and varies considerably with lung volume. Thus, if quantitative assessment of lung density within a breathhold is to be measured accurately, then it is necessary to also determine T(2) (*).

Pulmonary perfusion in the prone and supine postures in the normal human lung.

Prone posture increases cardiac output and improves pulmonary gas exchange. We hypothesized that, in the supine posture, greater compression of dependent lung limits regional blood flow. To test this, MRI-based measures of regional lung density, MRI arterial spin labeling quantification of pulmonary perfusion, and density-normalized perfusion were made in six healthy subjects. Measurements were made in both the prone and supine posture at functional residual capacity. Data were acquired in three nonoverlapping 15-mm sagittal slices covering most of the right lung: central, middle, and lateral, which were further divided into vertical zones: anterior, intermediate, and posterior. The density of the entire lung was not different between prone and supine, but the increase in lung density in the anterior lung with prone posture was less than the decrease in the posterior lung (change: +0.07 g/cm(3) anterior, -0.11 posterior; P < 0.0001), indicating greater compression of dependent lung in supine posture, principally in the central lung slice (P < 0.0001). Overall, density-normalized perfusion was significantly greater in prone posture (7.9 +/- 3.6 ml.min(-1).g(-1) prone, 5.1 +/- 1.8 supine, a 55% increase; P < 0.05) and showed the largest increase in the posterior lung as it became nondependent (change: +71% posterior, +58% intermediate, +31% anterior; P = 0.08), most marked in the central lung slice (P < 0.05). These data indicate that central posterior portions of the lung are more compressed in the supine posture, likely by the heart and adjacent structures, than are central anterior portions in the prone and that this limits regional perfusion in the supine posture.