Validating the distribution of specific ventilation in healthy humans measured using proton MR imaging

Assessing the pulmonary vascular responsiveness to oxygen with proton MRI

Ventilation-perfusion matching occurs passively and is also actively regulated through hypoxic pulmonary vasoconstriction (HPV). The extent of HPV activity in humans, particularly normal subjects, is uncertain. Current evaluation of HPV assesses changes in ventilation-perfusion relationships/pulmonary vascular resistance with hypoxia and is invasive, or unsuitable for patients because of safety concerns. We used a noninvasive imaging-based approach to quantify the pulmonary vascular response to oxygen as a metric of HPV by measuring perfusion changes between breathing 21% and 30%O2 using arterial spin labeling (ASL) MRI. We hypothesized that the differences between 21% and 30%O2 images reflecting HPV release would be 1) significantly greater than the differences without [Formula: see text] changes (e.g., 21-21% and 30-30%O2) and 2) negatively associated with ventilation-perfusion mismatch. Perfusion was quantified in the right lung in normoxia (baseline), after 15 min of 30% O2 breathing (hyperoxia) and 15 min normoxic recovery (recovery) in healthy subjects (7 M, 7 F; age = 41.4 ± 19.6 yr). Normalized, smoothed, and registered pairs of perfusion images were subtracted and the mean square difference (MSD) was calculated. Separately, regional alveolar ventilation and perfusion were quantified from specific ventilation, proton density, and ASL imaging; the spatial variance of ventilation-perfusion (σ2V̇a/Q̇) distributions was calculated. The O2-responsive MSD was reproducible (R2 = 0.94, P < 0.0001) and greater (0.16 ± 0.06, P < 0.0001) than that from subtracted images collected under the same [Formula: see text] (baseline = 0.09 ± 0.04, hyperoxia = 0.08 ± 0.04, recovery = 0.08 ± 0.03), which were not different from one another (P = 0.2). The O2-responsive MSD was correlated with σ2V̇a/Q̇ (R2 = 0.47, P = 0.007). These data suggest that active HPV optimizes ventilation-perfusion matching in normal subjects. This noninvasive approach could be applied to patients with different disease phenotypes to assess HPV and ventilation-perfusion mismatch.NEW & NOTEWORTHY We developed a new proton MRI method to noninvasively quantify the pulmonary vascular response to oxygen. Using a hyperoxic stimulus to release HPV, we quantified the resulting redistribution of perfusion. The differences between normoxic and hyperoxic images were greater than those between images without [Formula: see text] changes and negatively correlated with ventilation-perfusion mismatch. This suggests that active HPV optimizes ventilation-perfusion matching in normal subjects. This approach is suitable for assessing patients with different disease phenotypes.

Review of oxygen-enhanced lung mri: Pulse sequences for image acquisition and T1 measurement

Oxygen-enhanced MR imaging (OE-MRI) is a special proton imaging technique that can be performed without modifying the scanner hardware. Many fundamental studies have been conducted following the initial reporting of this technique in 1996, illustrating the high potential for its clinical application. This review aims to summarise and analyse current pulse sequences and T1 measurement methods for OE-MRI, including fundamental theories, existing pulse sequences applied to OE-MRI acquisition and T1 mapping. Wash-in and wash-out time identify lung function and are sensitive to ventilation; thus, dynamic OE-MRI is also discussed in this review. We compare OE-MRI with the primary competitive technique, hyperpolarised gas MRI. Finally, an overview of lower-field applications of OE-MRI is highlighted, as relatively recent publications demonstrated positive results. Lower-field OE-MRI, which is lower than 1.5 T, could be an alternative modality for detecting lung diseases. This educational review is aimed at researchers who want a quick summary of the steps needed to perform pulmonary OE-MRI with a particular focus on sequence design, settings, and quantification methods.

Magnetic Resonance Imaging of Aerosol Deposition

Nuclear magnetic resonance imaging (MRI) uses non-ionizing radiation and offers a host of contrast mechanisms with the potential to quantify aerosol deposition. This chapter introduces the physics of MRI, its use in lung imaging, and more specifically, the methods that are used for the detection of regional distributions of inhaled particles. The most common implementation of MRI is based on imaging of hydrogen atoms (1H) in water. The regional deposition of aerosol particles can be measured by the perturbation of the acquired 1H signals via labeling of the aerosol with contrast agents. Existing in vitro human and in vivo animal model measurements of regional aerosol deposition in the respiratory tract are described, demonstrating the capability of MRI to assess aerosol deposition in the lung.

Quantitative Imaging Metrics for the Assessment of Pulmonary Pathophysiology: An Official American Thoracic Society and Fleischner Society Joint Workshop Report

Multiple thoracic imaging modalities have been developed to link structure to function in the diagnosis and monitoring of lung disease. Volumetric computed tomography (CT) renders three-dimensional maps of lung structures and may be combined with positron emission tomography (PET) to obtain dynamic physiological data. Magnetic resonance imaging (MRI) using ultrashort-echo time (UTE) sequences has improved signal detection from lung parenchyma; contrast agents are used to deduce airway function, ventilation-perfusion-diffusion, and mechanics. Proton MRI can measure regional ventilation-perfusion ratio. Quantitative imaging (QI)-derived endpoints have been developed to identify structure-function phenotypes, including air-blood-tissue volume partition, bronchovascular remodeling, emphysema, fibrosis, and textural patterns indicating architectural alteration. Coregistered landmarks on paired images obtained at different lung volumes are used to infer airway caliber, air trapping, gas and blood transport, compliance, and deformation. This document summarizes fundamental "good practice" stereological principles in QI study design and analysis; evaluates technical capabilities and limitations of common imaging modalities; and assesses major QI endpoints regarding underlying assumptions and limitations, ability to detect and stratify heterogeneous, overlapping pathophysiology, and monitor disease progression and therapeutic response, correlated with and complementary to, functional indices. The goal is to promote unbiased quantification and interpretation of in vivo imaging data, compare metrics obtained using different QI modalities to ensure accurate and reproducible metric derivation, and avoid misrepresentation of inferred physiological processes. The role of imaging-based computational modeling in advancing these goals is emphasized. Fundamental principles outlined herein are critical for all forms of QI irrespective of acquisition modality or disease entity.

Phase-cycled balanced SSFP imaging for non-contrast-enhanced functional lung imaging

PURPOSE

To introduce phase-cycled balanced SSFP (bSSFP) acquisition as an alternative in Fourier decomposition MRI for improved robustness against field inhomogeneities.

METHODS

Series 2D dynamic lung images were acquired in 5 healthy volunteers at 1.5 T and 3 T using bSSFP sequence with multiple RF phase increments and compared with conventional single RF phase increment acquisitions. The approach was evaluated based on functional map homogeneity analysis, while ensuring image and functional map quality by means of SNR and contrast-to-noise ratio analyses.

RESULTS

At both field strengths, functional maps obtained with phase-cycled acquisitions displayed improved robustness against local signal losses compared with single-phase acquisitions. The coefficient of variation (mean ± SD, across volunteers) measured in the ventilation maps resulted in 29.7 ± 2.6 at 1.5 T and 37.5 ± 3.1 at 3 T for phase-cycled acquisitions, compared with 39.9 ± 5.2 at 1.5 T and 49.5 ± 3.7 at 3 T for single-phase acquisitions, indicating a significant improvement ( p < 0.05 $$ p<0.05 $$ ) in ventilation map homogeneity.

CONCLUSIONS

Phase-cycled bSSFP acquisitions improve robustness against field inhomogeneity artifacts and significantly improve ventilation map homogeneity at both field strengths. As such, phase-cycled bSSFP may serve as a robust alternative in lung function assessments.

Measuring short-term changes in specific ventilation using dynamic specific ventilation imaging

Specific ventilation imaging (SVI) measures the spatial distribution of specific ventilation (SV) in the lung with MRI by using inhaled oxygen as a contrast agent. Because of the inherently low signal to noise ratio in the technique, multiple switches between inspiring air and O₂ are utilized, and the high spatial resolution SV distribution determined as an average over the entire imaging period (~20 minutes). We hypothesized that a trade-off between spatial and temporal resolution could allow imaging at a higher temporal resolution, at the cost of a coarser, yet acceptable, spatial resolution. The appropriate window length and spatial resolution compromise was determined by generating synthetic data with signal- and contrast-to-noise characteristics reflective of that in previously published experimental data, with a known and unchanging distribution of SV, and showed that acceptable results could be obtained in an imaging period of ~7 minutes (80 breaths), with a spatial resolution of ~1cm³. Previously published data were then reanalyzed. The average heterogeneity of the temporally resolved maps of SV were not different to the previous overall analysis, however the temporally resolved maps were less effective at detecting the amount of bronchoconstriction resulting from methacholine administration. The results further indicated that the initial response to inhaled methacholine and subsequent inhalation of albuterol were largely complete within ~22 minutes and ~9 minutes respectively, although there was a tendency for an ongoing developing effect in both cases. These results suggest that it is feasible to use a shortened SVI protocol, with a modest sacrifice in spatial resolution, in order to measure temporally dynamic processes.

Simulation-based optimisation to quantify heterogeneity of specific ventilation and perfusion in the lung by the Inspired Sinewave Test

The degree of specific ventilatory heterogeneity (spatial unevenness of ventilation) of the lung is a useful marker of early structural lung changes which has the potential to detect early-onset disease. The Inspired Sinewave Test (IST) is an established noninvasive 'gas-distribution' type of respiratory test capable of measuring the cardiopulmonary parameters. We developed a simulation-based optimisation for the IST, with a simulation of a realistic heterogeneous lung, namely a lognormal distribution of spatial ventilation and perfusion. We tested this method in datasets from 13 anaesthetised pigs (pre and post-lung injury) and 104 human subjects (32 healthy and 72 COPD subjects). The 72 COPD subjects were classified into four COPD phenotypes based on 'GOLD' classification. This method allowed IST to identify and quantify heterogeneity of both ventilation and perfusion, permitting diagnostic distinction between health and disease states. In healthy volunteers, we show a linear relationship between the ventilatory heterogeneity versus age ([Formula: see text]). In a mechanically ventilated pig, IST ventilatory heterogeneity in noninjured and injured lungs was significantly different (p < 0.0001). Additionally, measured indices could accurately identify patients with COPD (area under the receiver operating characteristic curve is 0.76, p < 0.0001). The IST also could distinguish different phenotypes of COPD with 73% agreement with spirometry.

Utilization of 19F MRI for Identification of Iraq-Afghanistan War Lung Injury

INTRODUCTION

There is mounting evidence of respiratory problems related to military service in the Middle East in the past two decades due to environmental exposures during deployment (eg, sand storms and burn pits). This pilot study tests the hypothesis that regional lung function in subjects with prior deployment in Iraq and/or Afghanistan with suspected War Lung Injury (WLI) would be worse than subjects with normal lung function.

MATERIALS AND METHODS

Five subjects meeting the inclusion and exclusion criteria were recruited for this pilot study. All subjects underwent spirometry, high-resolution chest computed tomography imaging, and 19F MRI.

RESULTS

While the WLI subjects had normal pulmonary function tests and normal high-resolution chest computed tomography evaluations, their regional lung function from 19F MRI was abnormal with compartments with poor function showing slower filling time constants for ventilation. The scans of suspected WLI subjects show higher fractional lung volume with slow filling compartments similar to patients with chronic obstructive pulmonary disease in contrast to normal subjects.

CONCLUSIONS

This is consistent with our premise that WLI results in abnormal lung function and reflects small airways dysfunction and suggests that we may be able to provide a more sensitive tool for evaluation of WLI suspected cases.

Vaping disrupts ventilation-perfusion matching in asymptomatic users

Inhalation of e-cigarette's aerosols (vaping) has the potential to disrupt pulmonary gas exchange, but the effects in asymptomatic users are unknown. We assessed ventilation-perfusion (V̇_A/Q̇) mismatch in asymptomatic e-cigarette users, using magnetic resonance imaging (MRI). We hypothesized that vaping induces V̇_A/Q̇ mismatch through alterations in both ventilation and perfusion distributions. Nine young, asymptomatic "Vapers" with >1-yr vaping history, and no history of cardiopulmonary disease, were imaged supine using proton MRI, to assess the right lung at baseline and immediately after vaping. Seven young "Controls" were imaged at baseline only. Relative dispersion (SD/means) was used to quantify the heterogeneity of the individual ventilation and perfusion distributions. V̇_A/Q̇ mismatch was quantified using the second moments of the ventilation and perfusion versus V̇_A/Q̇ ratio distributions, log scale, LogSDV̇, and LogSDQ̇, respectively, analogous to the multiple inert gas elimination technique. Spirometry was normal in both groups. Ventilation heterogeneity was similar between groups at baseline (Vapers, 0.43 ± 0.13; Controls, 0.51 ± 0.11; P = 0.13) but increased after vaping (to 0.57 ± 0.17; P = 0.03). Perfusion heterogeneity was greater (P = 0.04) in Vapers at baseline (0.53 ± 0.06) compared with Controls (0.44 ± 0.10) but decreased after vaping (to 0.42 ± 0.07; P = 0.005). Vapers had greater (P = 0.01) V̇_A/Q̇ mismatch at baseline compared with Controls (LogSDQ̇ = 0.61 ± 0.12 vs. 0.43 ± 0.12), which was increased after vaping (LogSDQ̇ = 0.73 ± 0.16; P = 0.03). V̇_A/Q̇ mismatch is greater in Vapers and worsens after vaping. This suggests subclinical alterations in lung function not detected by spirometry.NEW & NOTEWORTHY This research provides evidence of vaping-induced disruptions in ventilation-perfusion matching in young, healthy, asymptomatic adults with normal spirometry who habitually vape. The changes in ventilation and perfusion distributions, both at baseline and acutely after vaping, and the potential implications on hypoxic vasoconstriction are particularly relevant in understanding the pathogenesis of vaping-induced dysfunction. Our imaging-based approach provides evidence of potential subclinical alterations in lung function below thresholds of detection using spirometry.

Ventilatory heterogeneity in the normal human lung is unchanged by controlled breathing

Measurement of ventilation heterogeneity with the multiple-breath nitrogen washout (MBW) is usually performed using controlled breathing with a fixed tidal volume and breathing frequency. However, it is unclear whether controlled breathing alters the underlying ventilatory heterogeneity. We hypothesized that the width of the specific ventilation distribution (a measure of heterogeneity) would be greater in tests performed during free breathing compared with those performed using controlled breathing. Eight normal subjects (age range = 23-50 yr, 5 female/3 male) twice underwent MRI-based specific ventilation imaging consisting of five repeated cycles with the inspired gas switching between 21% and 100% O₂ every ~2 min (total imaging time = ~20 min). In each session, tests were performed with free breathing (FB, no constraints) and controlled breathing (CB) at a respiratory rate of 12 breaths/min and no tidal volume control. The specific ventilation (SV) distribution in a mid-sagittal slice of the right lung was calculated, and the heterogeneity was calculated as the full width at half max of a Gaussian distribution fitted on a log scale (SV width). Free breathing resulted in a range of breathing frequencies from 8.7 to 15.9 breaths/min (mean = 11.5 ± 2.2, P = 0.62, compared with CB). Heterogeneity (SV width) was unchanged by controlled breathing (FB: 0.38 ± 0.12; CB: 0.34 ± 0.09, P = 0.18, repeated-measures ANOVA). The imposition of a controlled breathing frequency did not significantly affect the heterogeneity of ventilation in the normal lung, suggesting that MBW and specific ventilation imaging as typically performed provide an unperturbed measure of ventilatory heterogeneity.NEW & NOTEWORTHY By using MRI-based specific ventilation imaging (SVI), we showed that the heterogeneity of specific ventilation was not different comparing free breathing and breathing with the imposition of a fixed breathing frequency of 12 breaths/min. Thus, multiple-breath washout and SVI as typically performed provide an unperturbed measure of ventilatory heterogeneity.

Ventilation/Perfusion Relationships and Gas Exchange: Measurement Approaches

Ventilation-perfusion ( V ˙ A / Q ˙ ) matching, the regional matching of the flow of fresh gas to flow of deoxygenated capillary blood, is the most important mechanism affecting the efficiency of pulmonary gas exchange. This article discusses the measurement of V ˙ A / Q ˙ matching with three broad classes of techniques: (i) those based in gas exchange, such as the multiple inert gas elimination technique (MIGET); (ii) those derived from imaging techniques such as single-photon emission computed tomography (SPECT), positron emission tomography (PET), magnetic resonance imaging (MRI), computed tomography (CT), and electrical impedance tomography (EIT); and (iii) fluorescent and radiolabeled microspheres. The focus is on the physiological basis of these techniques that provide quantitative information for research purposes rather than qualitative measurements that are used clinically. The fundamental equations of pulmonary gas exchange are first reviewed to lay the foundation for the gas exchange techniques and some of the imaging applications. The physiological considerations for each of the techniques along with advantages and disadvantages are briefly discussed. © 2020 American Physiological Society. Compr Physiol 10:1155-1205, 2020.

Ventilation-perfusion heterogeneity measured by the multiple inert gas elimination technique is minimally affected by intermittent breathing of 100% O2

Proton magnetic resonance (MR) imaging to quantify regional ventilation-perfusion ( V ˙ A / Q ˙ ) ratios combines specific ventilation imaging (SVI) and separate proton density and perfusion measures into a composite map. Specific ventilation imaging exploits the paramagnetic properties of O2 , which alters the local MR signal intensity, in an FI O2 -dependent manner. Specific ventilation imaging data are acquired during five wash-in/wash-out cycles of breathing 21% O2 alternating with 100% O2 over ~20 min. This technique assumes that alternating FI O2 does not affect V ˙ A / Q ˙ heterogeneity, but this is unproven. We tested the hypothesis that alternating FI O2 exposure increases V ˙ A / Q ˙ mismatch in nine patients with abnormal pulmonary gas exchange and increased V ˙ A / Q ˙ mismatch using the multiple inert gas elimination technique (MIGET).The following data were acquired (a) breathing air (baseline), (b) breathing alternating air/100% O2 during an emulated-SVI protocol (eSVI), and (c) 20 min after ambient air breathing (recovery). MIGET heterogeneity indices of shunt, deadspace, ventilation versus V ˙ A / Q ˙ ratio, LogSD V ˙ , and perfusion versus V ˙ A / Q ˙ ratio, LogSD Q ˙ were calculated. LogSD V ˙ was not different between eSVI and baseline (1.04 ± 0.39 baseline, 1.05 ± 0.38 eSVI, p = .84); but was reduced compared to baseline during recovery (0.97 ± 0.39, p = .04). There was no significant difference in LogSD Q ˙ across conditions (0.81 ± 0.30 baseline, 0.79 ± 0.15 eSVI, 0.79 ± 0.20 recovery; p = .54); Deadspace was not significantly different (p = .54) but shunt showed a borderline increase during eSVI (1.0% ± 1.0 baseline, 2.6% ± 2.9 eSVI; p = .052) likely from altered hypoxic pulmonary vasoconstriction and/or absorption atelectasis. Intermittent breathing of 100% O2 does not substantially alter V ˙ A / Q ˙ matching and if SVI measurements are made after perfusion measurements, any potential effects will be minimized.

Proton MRI of the Lung: How to Tame Scarce Protons and Fast Signal Decay

Pulmonary proton MRI techniques offer the unique possibility of assessing lung function and structure without the requirement for hyperpolarization or dedicated hardware, which is mandatory for multinuclear acquisition. Five popular approaches are presented and discussed in this review: 1) oxygen enhanced (OE)-MRI; 2) arterial spin labeling (ASL); 3) Fourier decomposition (FD) MRI and other related methods including self-gated noncontrast-enhanced functional lung (SENCEFUL) MR and phase-resolved functional lung (PREFUL) imaging; 4) dynamic contrast-enhanced (DCE) MRI; and 5) ultrashort TE (UTE) MRI. While DCE MRI is the most established and well-studied perfusion measurement, FD MRI offers a free-breathing test without any contrast agent and is predestined for application in patients with renal failure or with low compliance. Additionally, FD MRI and related methods like PREFUL and SENCEFUL can act as an ionizing radiation-free V/Q scan, since ventilation and perfusion information is acquired simultaneously during one scan. For OE-MRI, different concentrations of oxygen are applied via a facemask to assess the regional change in T₁ , which is caused by the paramagnetic property of oxygen. Since this change is governed by a combination of ventilation, diffusion, and perfusion, a compound functional measurement can be achieved with OE-MRI. The known problem of fast T₂ * decay of the lung parenchyma leading to a low signal-to-noise ratio is bypassed by the UTE acquisition strategy. Computed tomography (CT)-like images allow the assessment of lung structure with high spatial resolution without ionizing radiation. Despite these different branches of proton MRI, common trends are evident among pulmonary proton MRI: 1) free-breathing acquisition with self-gating; 2) application of UTE to preserve a stronger parenchymal signal; and 3) transition from 2D to 3D acquisition. On that note, there is a visible convergence of the different methods and it is not difficult to imagine that future methods will combine different aspects of the presented methods.

Heavy upright exercise increases ventilation-perfusion mismatch in the basal lung: indirect evidence for interstitial pulmonary edema

Ventilation-perfusion (V̇a/Q̇) mismatch during exercise may result from interstitial pulmonary edema if increased pulmonary vascular pressure causes fluid efflux into the interstitium. If present, the increased fluid may compress small airways or blood vessels, disrupting V̇a/Q̇ matching, but this is unproven. We hypothesized that V̇a/Q̇ mismatch would be greatest in basal lung following heavy upright exercise, consistent with hydrostatic forces favoring edema accumulation in the gravitationally dependent lung. We applied new tools to reanalyze previously published magnetic resonance imaging data to determine regional V̇a/Q̇ mismatch following 45 min of heavy upright exercise in six athletes (V̇o_2max = 61 ± 7 mL·kg^-1·min^-1). In the supine posture, regional alveolar ventilation and local perfusion were quantified from specific ventilation imaging, proton density, and arterial spin labeling data in a single sagittal slice of the right lung before exercise (PRE), 15 min after exercise (POST), and in recovery 60 min after exercise (REC). Indices of V̇a/Q̇ mismatch [second moments (log scale) of ventilation (LogSD_V) and perfusion (LogSD_Q) vs. V̇a/Q̇ distributions] were calculated for apical, middle, and basal lung thirds, which represent gravitationally nondependent, middle, and dependent regions, respectively, during upright exercise. LogSD_V increased after exercise only in the basal lung (PRE 0.46 ± 0.06, POST 0.57 ± 0.14, REC 0.55 ±0.14, P = 0.01). Similarly, LogSD_Q increased only in the basal lung (PRE 0.40 ± 0.06, POST 0.51 ± 0.10, REC 0.44 ± 0.09, P = 0.04). Increased V̇a/Q̇ mismatch in the basal lung after exercise is potentially consistent with interstitial pulmonary edema accumulating in gravitationally dependent lung during exercise.NEW & NOTEWORTHY We reanalyzed previously published MRI data with new tools and found increased ventilation-perfusion mismatch only in the basal lung of athletes following 45 min of cycling exercise. This is consistent with the development of interstitial edema in the gravitationally dependent lung during heavy exercise.

Quantitative Mapping of Specific Ventilation in the Human Lung using Proton Magnetic Resonance Imaging and Oxygen as a Contrast Agent

Specific ventilation imaging (SVI) is a functional magnetic resonance imaging technique capable of quantifying specific ventilation - the ratio of the fresh gas entering a lung region divided by the region's end-expiratory volume - in the human lung, using only inhaled oxygen as a contrast agent. Regional quantification of specific ventilation has the potential to help identify areas of pathologic lung function. Oxygen in solution in tissue shortens the tissue's longitudinal relaxation time (T1), and thus a change in tissue oxygenation can be detected as a change in T1-weighted signal with an inversion recovery acquired image. Following an abrupt change between two concentrations of inspired oxygen, the rate at which lung tissue within a voxel equilibrates to a new steady-state reflects the rate at which resident gas is being replaced by inhaled gas. This rate is determined by specific ventilation. To elicit this sudden change in oxygenation, subjects alternately breathe 20-breath blocks of air (21% oxygen) and 100% oxygen while in the MRI scanner. A stepwise change in inspired oxygen fraction is achieved through use of a custom three-dimensional (3D)-printed flow bypass system with a manual switch during a short end-expiratory breath hold. To detect the corresponding change in T1, a global inversion pulse followed by a single shot fast spin echo sequence was used to acquire two-dimensional T1-weighted images in a 1.5 T MRI scanner, using an eight-element torso coil. Both single slice and multi-slice imaging are possible, with slightly different imaging parameters. Quantification of specific ventilation is achieved by correlating the time-course of signal intensity for each lung voxel with a library of simulated responses to the air/oxygen stimulus. SVI estimations of specific ventilation heterogeneity have been validated against multiple breath washout and proved to accurately determine the heterogeneity of the specific ventilation distribution.

Free-breathing quantification of regional ventilation derived by phase-resolved functional lung (PREFUL) MRI

PURPOSE

To test the feasibility of regional fully quantitative ventilation measurement in free breathing derived by phase-resolved functional lung (PREFUL) MRI in the supine and prone positions. In addition, the influence of T₂ * relaxation time on ventilation quantification is assessed.

METHODS

Twelve healthy volunteers underwent functional MRI at 1.5 T using a 2D triple-echo spoiled gradient echo sequence allowing for quantitative measurement of T₂ * relaxation time. Minute ventilation (ΔV) was quantified by conventional fractional ventilation (FV) and the newly introduced regional ventilation (VR), which corrects volume errors due to image registration. ΔV_FV versus ΔV_VR and ΔV_VR versus ΔV_VR with T₂ * correction were compared using Bland-Altman plots and correlation analysis. The repeatability and physiological plausibility of all measurements were tested in the supine and prone positions.

RESULTS

On global and regional scales a strong correlation was observed between ΔV_FV versus ΔV_VR and ΔV_VR versus ΔV_VRT2* (r > 0.93); however, regional Bland-Altman analysis showed systematic differences (p < 0.0001). Unlike ΔV_VRT2* , ΔV_VR and ΔV_FV showed expected physiologic anterior-posterior gradients, which decreased in the supine but not in the prone position at second measurement during 3 min in the same position. For all quantification methods a moderate repeatability (coefficient of variation <20%) of ventilation was found.

CONCLUSION

A fully quantified regional ventilation measurement using ΔV_VR in free breathing is feasible and shows physiologically plausible results. In contrast to conventional ΔV_FV, volume errors due to image registration are eliminated with the ΔV_VR approach. However, correction for the T₂ * effect remains challenging.

Regional airflow obstruction after bronchoconstriction and subsequent bronchodilation in subjects without pulmonary disease

Some subjects with asthma have ventilation defects that are resistant to bronchodilator therapy, and it is thought that these resistant defects may be due to ongoing inflammation or chronic airway remodeling. However, it is unclear whether regional obstruction due to bronchospasm alone persists after bronchodilator therapy. To investigate this, six young, healthy subjects, in whom inflammation and remodeling were assumed to be absent, were bronchoconstricted with a PC₂₀ [the concentration of methacholine that elicits a 20% drop in forced expiratory volume in 1 s (FEV₁)] dose of methacholine and subsequently bronchodilated with a standard dose of albuterol on three separate occasions. Specific ventilation imaging, a proton MRI technique, was used to spatially map specific ventilation across 80% of each subject's right lung in each condition. The ratio between regional specific ventilation at baseline and after intervention was used to classify areas that had constricted. After albuterol rescue from methacholine bronchoconstriction, 12% (SD 9) of the lung was classified as constricted. Of the 12% of lung units that were classified as constricted after albuterol, approximately half [7% (SD 7)] had constricted after methacholine and failed to recover, whereas half [6% (SD 4)] had remained open after methacholine but became constricted after albuterol. The incomplete regional recovery was not reflected in the subjects' FEV₁ measurements, which did not decrease from baseline (P = 0.97), nor was it detectable as an increase in specific ventilation heterogeneity (P = 0.78).NEW & NOTEWORTHY In normal subjects bronchoconstricted with methacholine and subsequently treated with albuterol, not all regions of the healthy lung returned to their prebronchoconstricted specific ventilation after albuterol, despite full recovery of integrative lung indexes (forced expiratory volume in 1 s and specific ventilation heterogeneity). The regions that remained bronchoconstricted following albuterol were those with the highest specific ventilation at baseline, which suggests that they may have received the highest methacholine dose.

T1 , T1 contrast, and Ernst-angle images of four rat-lung pathologies

PURPOSE

To initiate the archive of relaxation-weighted images that may help discriminate between pulmonary pathologies relevant to acute respiratory distress syndrome. MRI has the ability to distinguish pathologies by providing a variety of different contrast mechanisms. Lungs have historically been difficult to image with MRI but image quality is sufficient to begin cataloging the appearance of pathologies in T₁ - and T₂ -weighted images. This study documents T₁ and the use of T₁ contrast with four experimental rat lung pathologies.

METHODS

Inversion-recovery and spoiled steady state images were made at 1.89 T to measure T₁ and document contrast in rats with atelectasis, lipopolysaccharide-induced inflammation, ventilator-induced lung injury (VILI), and injury from saline lavage. Higher-resolution Ernst-angle images were made to see patterns of lung infiltrations.

RESULTS

T₁ -weighted images showed minimal contrast between pathologies, similar to T₁ -weighted images of other soft tissues. Images taken shortly after magnetization inversion and displayed with inverted contrast highlight lung pathologies. Ernst-angle images distinguish the effects of T₁ relaxation and spin density and display distinctive patterns. T₁ for pathologies were: atelectasis, 1.25 ± 0.046 s; inflammation from instillation of lipopolysaccharide, 1.24 ± 0.015 s; VILI, 1.55 ± 0.064 s (p = 0.0022 vs. normal lung); and injury from saline lavage, 1.90±0.080 s (p = 0.0022 vs. normal lung; p = 0.0079 vs. VILI). T₁ of normal lung and erector spinae muscle were 1.25 ± 0.028 s and 1.02 ± 0.027 s, respectively (p = 0.0022).

CONCLUSIONS

Traditional T₁ -weighting is subtle. However, images made with inverted magnetization and inverted contrast highlight the pathologies and Ernst-angle images aid in distinguishing pathologies.

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.

The spatial pattern of methacholine bronchoconstriction recurs when supine, independently of posture during provocation, but does not recur between postures

The location of lung regions with compromised ventilation (often called ventilation defects) during a bronchoconstriction event may be influenced by posture. We aimed to determine the effect of prone vs. supine posture on the spatial pattern of methacholine-induced bronchoconstriction in six healthy adults (ages 21-41, three females) using specific ventilation imaging. Three postural conditions were chosen to assign the effect of posture to the drug administration and/or imaging phase of the experiment - supine methacholine administration followed by supine imaging, prone methacholine administration followed by supine imaging, and prone methacholine administration followed by prone imaging. The two conditions in which imaging was performed supine had similar spatial patterns of bronchoconstriction despite a change in posture during methacholine administration; the odds ratio for recurrent constriction was mean (SD) = 7.4 (3.9). Conversely, dissimilar spatial patterns of bronchoconstriction emerged when posture during imaging was changed; the odds ratio for recurrent constriction between the prone methacholine/supine imaging condition and the prone methacholine/prone imaging condition was 1.2 (0.9). Logistic regression showed that height above the dependent lung border was a significant negative predictor of constriction in the two supine imaging conditions (p<0.001 for each), but not in the prone imaging condition (p=0.20). These results show that the spatial pattern of methacholine bronchoconstriction is recurrent in the supine posture, regardless of whether methacholine is given prone or supine, but that prone posture during imaging eliminates that recurrent pattern and reduces its dependence on gravitational height.

Spatial persistence of reduced specific ventilation following methacholine challenge in the healthy human lung

Specific ventilation imaging was used to identify regions of the healthy lung (6 supine subjects, ages 21-41 yr, 3 men) that experienced a fall in specific ventilation following inhalation of methacholine. This test was repeated 1 wk later and 3 mo later to test for spatial recurrence. Our data showed that 53% confidence interval (CI; 46%, 59%) of volume elements that constricted during one methacholine challenge did so again in another and that this quantity did not vary with time; 46% CI (28%, 64%) recurred 1 wk later, and 56% CI (51%, 61%) recurred 3 mo later. Previous constriction was a strong predictor for future constriction. Volume elements that constricted during one challenge were 7.7 CI (5.2, 10.2) times more likely than nonconstricted elements to constrict in a second challenge, regardless of whether the second episode was 1 wk [7.7 CI (2.9, 12.4)] or 3 mo [7.7 CI (4.6, 10.8)] later. Furthermore, posterior lung elements were more likely to constrict following methacholine than anterior lung elements (volume fraction 0.43 ± 0.22 posterior vs. 0.10 ± 0.03 anterior; P = 0.005), and basal elements that constricted were more likely than their apical counterparts to do so persistently through all three trials (volume fraction 0.14 ± 0.04 basal vs. 0.04 ± 0.04 apical; P = 0.003). Taken together, this evidence suggests a physiological predisposition toward constriction in some lung elements, especially those located in the posterior and basal lung when the subject is supine. NEW & NOTEWORTHY The spatial pattern of bronchoconstriction following methacholine is persistent over time in healthy individuals, in whom chronic inflammation and airway remodeling are assumed to be absent. This suggests that regional lung inflation and airway structure may play dominant roles in determining the spatial pattern of methacholine bronchoconstriction.

In silico modeling of oxygen-enhanced MRI of specific ventilation

Specific ventilation imaging (SVI) proposes that using oxygen-enhanced 1H MRI to capture signal change as subjects alternatively breathe room air and 100% O2 provides an estimate of specific ventilation distribution in the lung. How well this technique measures SV and the effect of currently adopted approaches of the technique on resulting SV measurement is open for further exploration. We investigated (1) How well does imaging a single sagittal lung slice represent whole lung SV? (2) What is the influence of pulmonary venous blood on the measured MRI signal and resultant SVI measure? and (3) How does inclusion of misaligned images affect SVI measurement? In this study, we utilized two patient-based in silico models of ventilation, perfusion, and gas exchange to address these questions for normal healthy lungs. Simulation results from the two healthy young subjects show that imaging a single slice is generally representative of whole lung SV distribution, with a calculated SV gradient within 90% of that calculated for whole lung distributions. Contribution of O2 from the venous circulation results in overestimation of SV at a regional level where major pulmonary veins cross the imaging plane, resulting in a 10% increase in SV gradient for the imaging slice. A worst-case scenario simulation of image misalignment increased the SV gradient by 11.4% for the imaged slice.

Free-breathing Pulmonary MR Imaging to Quantify Regional Ventilation

Purpose To measure regional specific ventilation with free-breathing hydrogen 1 (¹H) magnetic resonance (MR) imaging without exogenous contrast material and to investigate correlations with hyperpolarized helium 3 (³He) MR imaging and pulmonary function test measurements in healthy volunteers and patients with asthma. Materials and Methods Subjects underwent free-breathing ¹H and static breath-hold hyperpolarized ³He MR imaging as well as spirometry and plethysmography; participants were consecutively recruited between January and June 2017. Free-breathing ¹H MR imaging was performed with an optimized balanced steady-state free-precession sequence; images were retrospectively grouped into tidal inspiration or tidal expiration volumes with exponentially weighted phase interpolation. MR imaging volumes were coregistered by using optical flow deformable registration to generate ¹H MR imaging-derived specific ventilation maps. Hyperpolarized ³He MR imaging- and ¹H MR imaging-derived specific ventilation maps were coregistered to quantify regional specific ventilation within hyperpolarized ³He MR imaging ventilation masks. Differences between groups were determined with the Mann-Whitney test and relationships were determined with Spearman (ρ) correlation coefficients. Statistical analyses were performed with software. Results Thirty subjects (median age: 50 years; interquartile range [IQR]: 30 years), including 23 with asthma and seven healthy volunteers, were evaluated. Both ¹H MR imaging-derived specific ventilation and hyperpolarized ³He MR imaging-derived ventilation percentage were significantly greater in healthy volunteers than in patients with asthma (specific ventilation: 0.14 [IQR: 0.05] vs 0.08 [IQR: 0.06], respectively, P < .0001; ventilation percentage: 99% [IQR: 1%] vs 94% [IQR: 5%], P < .0001). For all subjects, ¹H MR imaging-derived specific ventilation correlated with plethysmography-derived specific ventilation (ρ = 0.54, P = .002) and hyperpolarized ³He MR imaging-derived ventilation percentage (ρ = 0.67, P < .0001) as well as with forced expiratory volume in 1 second (FEV₁) (ρ = 0.65, P = .0001), ratio of FEV₁ to forced vital capacity (ρ = 0.75, P < .0001), ratio of residual volume to total lung capacity (ρ = -0.68, P < .0001), and airway resistance (ρ = -0.51, P = .004). ¹H MR imaging-derived specific ventilation was significantly greater in the gravitational-dependent versus nondependent lung in healthy subjects (P = .02) but not in patients with asthma (P = .1). In patients with asthma, coregistered ¹H MR imaging specific ventilation and hyperpolarized ³He MR imaging maps showed that specific ventilation was diminished in corresponding ³He MR imaging ventilation defects (0.05 ± 0.04) compared with well-ventilated regions (0.09 ± 0.05) (P < .0001). Conclusion ¹H MR imaging-derived specific ventilation correlated with plethysmography-derived specific ventilation and ventilation defects seen by using hyperpolarized ³He MR imaging. ^© RSNA, 2018 Online supplemental material is available for this article.

Pulmonary Delivery of Therapeutic and Diagnostic Gases

The 21st Congress for the International Society for Aerosols in Medicine included, for the first time, a session on Pulmonary Delivery of Therapeutic and Diagnostic Gases. The rationale for such a session within ISAM is that the pulmonary delivery of gaseous drugs in many cases targets the same therapeutic areas as aerosol drug delivery, and is in many scientific and technical aspects similar to aerosol drug delivery. This article serves as a report on the recent ISAM congress session providing a synopsis of each of the presentations. The topics covered are the conception, testing, and development of the use of nitric oxide to treat pulmonary hypertension; the use of realistic adult nasal replicas to evaluate the performance of pulsed oxygen delivery devices; an overview of several diagnostic gas modalities; and the use of inhaled oxygen as a proton magnetic resonance imaging (MRI) contrast agent for imaging temporal changes in the distribution of specific ventilation during recovery from bronchoconstriction. Themes common to these diverse applications of inhaled gases in medicine are discussed, along with future perspectives on development of therapeutic and diagnostic gases.

Quantification of regional deformation of the lungs by non-rigid registration of three-dimensional contrast-enhanced magnetic resonance imaging

BACKGROUND

Assessment of lung function is vital for the diagnosis of a variety of pathological conditions. Research has been proposed to study pulmonary mechanics and kinematics using two-dimensional (2D) magnetic resonance imaging (MRI). This allows estimation of regional lung tissue mechanics but is limited to 2D information. An approach based on three-dimensional (3D) contrast-enhanced MR angiogram of pulmonary blood vessels and a non-rigid image registration technique is proposed for quantification of lung regional deformations, which can potentially be used for assessment of pulmonary parenchymal mechanics and regional ventilation for disease diagnosis without ionizing radiation.

METHODS

On three volunteers, an end-expiration scan and end-inspiration scan was acquired successively for each volunteer using a 3D breath-hold contrast-enhanced MRI sequence several minutes after gadolinium injection. Subsequently, a rectangle box lung mask is manually selected for each end-expiration scan, applying non-rigid registration algorithms using cubic B-splines as transformations to align each pair of images. This incorporates the Normalized Correlation Coefficient similarity with the bending energy term as cost function with a multi-resolution multi-grid approach. Finally, the lung regional 3D deformations were obtained using the transformations obtained by registration. The alignment accuracy after non-rigid registration was estimated by using a set of branch points of pulmonary blood vessels as anatomical landmarks for each pair of images.

RESULTS

With contrast enhancement, the pulmonary blood vessel signal was enhanced, which greatly facilitated the non-rigid registration in the lung parenchyma. The average landmarks distances in three pairs of datasets are reduced from 17.9, 20.3 and 16.3 mm, to 1.0, 1.6 and 1.2 mm, respectively, by non-rigid registration. After registration, the average distances error of each pair of datasets was less than 0.6 mm in the right-to-left (RL) direction, less than 0.9 mm in the inferior-to-superior (IS) direction, and less than 1.2 mm in the anterior-to-posterior (AP) direction.

CONCLUSIONS

Results demonstrated that the proposed method can accurately register lungs with large deformations to evaluate lung regional deformation. It may be used for quantitative assessment of 3D lung regional ventilation avoiding ionizing radiation.

Ventilation heterogeneity measured by multiple breath inert gas testing is not affected by inspired oxygen concentration in healthy humans

Multiple breath washout (MBW) and oxygen-enhanced MRI techniques use acute exposure to 100% oxygen to measure ventilation heterogeneity. Implicit is the assumption that breathing 100% oxygen does not induce changes in ventilation heterogeneity; however, this is untested. We hypothesized that ventilation heterogeneity decreases with increasing inspired oxygen concentration in healthy subjects. We performed MBW in 8 healthy subjects (4 women, 4 men; age = 43 ± 15 yr) with normal pulmonary function (FEV₁ = 98 ± 6% predicted) using 10% argon as a tracer gas and oxygen concentrations of 12.5%, 21%, or 90%. MBW was performed in accordance with ERS-ATS guidelines. Subjects initially inspired air followed by a wash-in of test gas. Tests were performed in balanced order in triplicate. Gas concentrations were measured at the mouth, and argon signals rescaled to mimic a N₂ washout, and analyzed to determine the distribution of specific ventilation (SV). Heterogeneity was characterized by the width of a log-Gaussian fit of the SV distribution and from S_acin and S_cond indexes derived from the phase III slope. There were no significant differences in the ventilation heterogeneity due to altered inspired oxygen: histogram width (hypoxia 0.57 ± 0.11, normoxia 0.60 ± 0.08, hyperoxia 0.59 ± 0.09, P = 0.51), S_cond (hypoxia 0.014 ± 0.011, normoxia 0.012 ± 0.015, hyperoxia 0.010 ± 0.011, P = 0.34), or S_acin (hypoxia 0.11 ± 0.04, normoxia 0.10 ± 0.03, hyperoxia 0.12 ± 0.03, P = 0.23). Functional residual capacity was increased in hypoxia (P = 0.04) and dead space increased in hyperoxia (P = 0.0001) compared with the other conditions. The acute use of 100% oxygen in MBW or MRI is unlikely to affect ventilation heterogeneity.NEW & NOTEWORTHY Hyperoxia is used to measure the distribution of ventilation in imaging and MBW but may alter the underlying ventilation distribution. We used MBW to evaluate the effect of inspired oxygen concentration on the ventilation distribution using 10% argon as a tracer. Short-duration exposure to hypoxia (12.5% oxygen) and hyperoxia (90% oxygen) during MBW had no significant effect on ventilation heterogeneity, suggesting that hyperoxia can be used to assess the ventilation distribution.

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.

Regional Ventilation Is the Main Determinant of Alveolar Deposition of Coarse Particles in the Supine Healthy Human Lung During Tidal Breathing

BACKGROUND

To quantify the relationship between regional lung ventilation and coarse aerosol deposition in the supine healthy human lung, we used oxygen-enhanced magnetic resonance imaging and planar gamma scintigraphy in seven subjects.

METHODS

Regional ventilation was measured in the supine posture in a 15 mm sagittal slice of the right lung. Deposition was measured by using planar gamma scintigraphy (coronal scans, 40 cm FOV) immediately postdeposition, 1 hour 30 minutes and 22 hours after deposition of ^99mTc-labeled particles (4.9 μm MMAD, GSD 2.5), inhaled in the supine posture (flow 0.5 L/s, 15 breaths/min). The distribution of retained particles at different times was used to infer deposition in different airway regions, with 22 hours representing alveolar deposition. The fraction of total slice ventilation per quartile of lung height from the lung apex to the dome of the diaphragm at functional residual capacity was computed, and co-registered with deposition data-apices aligned-using a transmission scan as reference. The ratio of fractional alveolar deposition to fractional ventilation of each quartile (r) was used to evaluate ventilation and deposition matching (r > 1, regional aerosol deposition fraction larger than regional ventilation fraction).

RESULTS

r was not significantly different from 1 for all regions (1.04 ± 0.25, 1.08 ± 0.22, 1.03 ± 0.17, 0.92 ± 0.13, apex to diaphragm, p > 0.40) at the alveolar level (r_22h). For retention times r_0h and r_1h30, only the diaphragmatic region at r_1h30 differed significantly from 1.

CONCLUSIONS

These results support the hypothesis that alveolar deposition is directly proportional to ventilation for ∼5 μm particles that are inhaled in the supine posture and are consistent with previous simulation predictions that show that convective flow is the main determinant of aerosol transport to the lung periphery.

Susceptibility to high-altitude pulmonary edema is associated with a more uniform distribution of regional specific ventilation

High-altitude pulmonary edema (HAPE) is a potentially fatal condition affecting high-altitude sojourners. The biggest predictor of HAPE development is a history of prior HAPE. Magnetic resonance imaging (MRI) shows that HAPE-susceptible (with a history of HAPE), but not HAPE-resistant (with a history of repeated ascents without illness) individuals develop greater heterogeneity of regional pulmonary perfusion breathing hypoxic gas (O₂ = 12.5%), consistent with uneven hypoxic pulmonary vasoconstriction (HPV). Why HPV is uneven in HAPE-susceptible individuals is unknown but may arise from regionally heterogeneous ventilation resulting in an uneven stimulus to HPV. We tested the hypothesis that ventilation is more heterogeneous in HAPE-susceptible subjects (n = 6) compared with HAPE-resistant controls (n = 7). MRI specific ventilation imaging (SVI) was used to measure regional specific ventilation and the relative dispersion (SD/mean) of SVI used to quantify baseline heterogeneity. Ventilation heterogeneity from conductive and respiratory airways was measured in normoxia and hypoxia (O₂ = 12.5%) using multiple-breath washout and heterogeneity quantified from the indexes S_cond and S_acin, respectively. Contrary to our hypothesis, HAPE-susceptible subjects had significantly lower relative dispersion of specific ventilation than the HAPE-resistant controls [susceptible = 1.33 ± 0.67 (SD), resistant = 2.36 ± 0.98, P = 0.05], and S_acin tended to be more uniform (susceptible = 0.085 ± 0.009, resistant = 0.113 ± 0.030, P = 0.07). S_cond was not significantly different between groups (susceptible = 0.019 ± 0.007, resistant = 0.020 ± 0.004, P = 0.67). S_acin and S_cond did not change significantly in hypoxia (P = 0.56 and 0.19, respectively). In conclusion, ventilation heterogeneity does not change with short-term hypoxia irrespective of HAPE susceptibility, and lesser rather than greater ventilation heterogeneity is observed in HAPE-susceptible subjects. This suggests that the basis for uneven HPV in HAPE involves vascular phenomena.NEW & NOTEWORTHY Uneven hypoxic pulmonary vasoconstriction (HPV) is thought to incite high-altitude pulmonary edema (HAPE). We evaluated whether greater heterogeneity of ventilation is also a feature of HAPE-susceptible subjects compared with HAPE-resistant subjects. Contrary to our hypothesis, ventilation heterogeneity was less in HAPE-susceptible subjects and unaffected by hypoxia, suggesting a vascular basis for uneven HPV.

Using Hyperpolarized 129Xe MRI to Quantify the Pulmonary Ventilation Distribution

RATIONALE AND OBJECTIVES

Ventilation heterogeneity is impossible to detect with spirometry. Alternatively, pulmonary ventilation can be imaged three-dimensionally using inhaled ¹²⁹Xe magnetic resonance imaging (MRI). To date, such images have been quantified primarily based on ventilation defects. Here, we introduce a robust means to transform ¹²⁹Xe MRI scans such that the underlying ventilation distribution and its heterogeneity can be quantified.

MATERIALS AND METHODS

Quantitative ¹²⁹Xe ventilation MRI was conducted in 12 younger (24.7 ± 5.2 years) and 10 older (62.2 ± 7.2 years) healthy individuals, as well as in 9 younger (25.9 ± 6.4 yrs) and 10 older (63.2 ± 6.1 years) asthmatics. The younger healthy population was used to establish a reference ventilation distribution and thresholds for six intensity bins. These bins were used to display and quantify the ventilation defect region (VDR), the low ventilation region (LVR), and the high ventilation region (HVR).

RESULTS

The ventilation distribution in young subjects was roughly Gaussian with a mean and standard deviation of 0.52 ± 0.18, resulting in VDR = 2.1 ± 1.3%, LVR = 15.6 ± 5.4%, and HVR = 17.4 ± 3.1%. Older healthy volunteers exhibited a significantly right-skewed distribution (0.46 ± 0.20, P = 0.034), resulting in significantly increased VDR (7.0 ± 4.8%, P = 0.008) and LVR (24.5 ± 11.5%, P = 0.025). In the asthmatics, VDR and LVR increased in the older population, and HVR was significantly reduced (13.5 ± 4.6% vs 18.9 ± 4.5%, P = 0.009). Quantitative ¹²⁹Xe MRI also revealed altered ventilation heterogeneity in response to albuterol in two asthmatics with normal spirometry.

CONCLUSIONS

Quantitative ¹²⁹Xe MRI provides a robust and objective means to display and quantify the pulmonary ventilation distribution, even in subjects who have airway function impairment not appreciated by spirometry.

Physiology for the pulmonary functional imager

As pulmonary functional imaging moves beyond the realm of the radiologist and physicist, it is important that imagers have a common language and understanding of the relevant physiology of the lung. This review will focus on key physiological concepts and pitfalls relevant to functional lung imaging.

Improving the Accuracy of Registration-Based Biomechanical Analysis: A Finite Element Approach to Lung Regional Strain Quantification

Tissue deformation plays an important role in lung physiology, as lung parenchyma largely deforms during spontaneous ventilation. However, excessive regional deformation may lead to lung injury, as observed in patients undergoing mechanical ventilation. Thus, the accurate estimation of regional strain has recently received great attention in the intensive care community. In this work, we assess the accuracy of regional strain maps computed from direct differentiation of B-Spline (BS) interpolations, a popular technique employed in non-rigid registration of lung computed tomography (CT) images. We show that, while BS-based registration methods give excellent results for the deformation transformation, the strain field directly computed from BS derivatives results in predictions that largely oscillate, thus introducing important errors that can even revert the sign of strain. To alleviate such spurious behavior, we present a novel finite-element (FE) method for the regional strain analysis of lung CT images. The method follows from a variational strain recovery formulation, and delivers a continuous approximation to the strain field in arbitrary domains. From analytical benchmarks, we show that the FE method results in errors that are a fraction of those found for the BS method, both in an average and pointwise sense. The application of the proposed FE method to human lung CT images results in 3D strain maps are heterogeneous and smooth, showing high consistency with specific ventilation maps reported in the literature. We envision that the proposed FE method will considerably improve the accuracy of image-based biomechanical analysis, making it reliable enough for routine medical applications.

Rapid Prototyping of Inspired Gas Delivery System for Pulmonary MRI Research

Specific ventilation imaging (SVI) is a noninvasive magnetic resonance imaging (MRI)-based method for determining the regional distribution of inspired air in the lungs, useful for the assessment of pulmonary function in medical research. This technique works by monitoring the rate of magnetic resonance signal change in response to a series of imposed step changes in inspired oxygen concentration. The current SVI technique requires a complex system of tubes, valves, and electronics that are used to supply and rapidly switch inspired gases while subjects are imaged, which makes the technique difficult to translate into the clinical setting. This report discusses the design and implementation of custom three-dimensional (3D) printed hardware that greatly simplifies SVI measurement of lung function. Several hardware prototypes were modeled using computer-aided design software and printed for evaluation. After finalization of the design, the new delivery system was evaluated based on O₂ and N₂ concentration step responses and validated against the current SVI protocol. The design performed rapid switching of supplied gas within 250 ms and consistently supplied the desired concentration of O₂ during operation. It features a reduction in the number of commercial hardware components, from five to one, and a reduction in the number of gas lines between the operator's room and the scanner room, from four to one, as well as a substantially reduced preparation time from 25 to 5 min. 3D printing is well suited to the design of inexpensive custom MRI compatible hardware, making it particularly useful in imaging-based research.