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

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.

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.

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.

Airway closure is the predominant physiological mechanism of low ventilation seen on hyperpolarized helium-3 MRI lung scans

Hyperpolarized helium-3 MRI (³He MRI) provides detailed visualization of low- (hypo- and non-) ventilated lungs. Physiological measures of gas mixing may be assessed by multiple breath nitrogen washout (MBNW) and of airway closure by a forced oscillation technique (FOT). We hypothesize that in patients with asthma, areas of low-ventilated lung on ³He MRI are the result of airway closure. Ten control subjects, ten asthma subjects with normal spirometry (non-obstructed), and ten asthmatic subjects with reduced baseline lung function (obstructed) attended two testing sessions. On visit one, baseline plethysmography was performed followed by spirometry, MBNW, and FOT assessment pre and post methacholine challenge. On visit two, ³He MRI scans were conducted pre and post methacholine challenge. Post methacholine the volume of low-ventilated lung increased from 8.3% to 13.8% in the non-obstructed group (P = 0.012) and from 13.0% to 23.1% in the obstructed group (P = 0.001). For all subjects, the volume of low ventilation from ³He MRI correlated with a marker of airway closure in obstructive subjects, Xrs (6 Hz) and the marker of ventilation heterogeneity Scond with r² values of 0.61 (P < 0.001) and 0.56 (P < 0.001), respectively. The change in Xrs (6 Hz) correlated well (r² = 0.45, p < 0.001), whereas the change in Scond was largely independent of the change in low ventilation volume (r² = 0.13, P < 0.01). The only significant predictor of low ventilation volume from the multi-variate analysis was Xrs (6 Hz). This is consistent with the concept that regions of poor or absent ventilation seen on ³He MRI are primarily the result of airway closure.NEW & NOTEWORTHY This study introduces a novel technique of generating high-resolution 3D ventilation maps from hyperpolarized helium-3 MRI. It is the first study to demonstrate that regions of poor or absent ventilation seen on ³He MRI are primarily the result of airway closure.

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.

Obesity alters the topographical distribution of ventilation and the regional response to bronchoconstriction

Obesity is associated with reduced operating lung volumes that may contribute to increased airway closure during tidal breathing and abnormalities in ventilation distribution. We investigated the effect of obesity on the topographical distribution of ventilation before and after methacholine-induced bronchoconstriction using single-photon emission computed tomography (SPECT)-computed tomography (CT) in healthy subjects. Subjects with obesity ( n = 9) and subjects without obesity ( n = 10) underwent baseline and postbronchoprovocation SPECT-CT imaging, in which Technegas was inhaled upright and followed by supine scanning. Lung regions that were nonventilated (Ventnon), low ventilated (Ventlow), or well ventilated (Ventwell) were calculated using an adaptive threshold method and were expressed as a percentage of total lung volume. To determine regional ventilation, lungs were divided into upper, middle, and lower thirds of axial length, derived from CT. At baseline, Ventnon and Ventlow for the entire lung were similar in subjects with and without obesity. However, in the upper lung zone, Ventnon (17.5 ± 10.6% vs. 34.7 ± 7.8%, P < 0.001) and Ventlow (25.7 ± 6.3% vs. 33.6 ± 5.1%, P < 0.05) were decreased in subjects with obesity, with a consequent increase in Ventwell (56.8 ± 9.2% vs. 31.7 ± 10.1%, P < 0.001). The greater diversion of ventilation to the upper zone was correlated with body mass index ( rs = 0.74, P < 0.001), respiratory system resistance ( rs = 0.72, P < 0.001), and respiratory system reactance ( rs = −0.64, P = 0.003) but not with lung volumes or basal airway closure. Following bronchoprovocation, overall Ventnon increased similarly in both groups; however, in subjects without obesity, Ventnon only increased in the lower zone, whereas in subjects with obesity, Ventnon increased more evenly across all lung zones. In conclusion, obesity is associated with altered ventilation distribution during baseline and following bronchoprovocation, independent of reduced lung volumes.

NEW & NOTEWORTHY Using ventilation SPECT-computed tomography imaging in healthy subjects, we demonstrate that ventilation in obesity is diverted to the upper lung zone and that this is strongly correlated with body mass index but is independent of operating lung volumes and of airway closure. Furthermore, methacholine-induced bronchoconstriction only occurred in the lower lung zone in individuals who were not obese, whereas in subjects who were obese, it occurred more evenly across all lung zones. These findings show that obesity-associated factors alter the topographical distribution of ventilation.

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.

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.

Airway Transmural Pressures in an Airway Tree During Bronchoconstriction in Asthma

Airway transmural pressure in healthy homogeneous lungs with dilated airways is approximately equal to the difference between intraluminal and pleural pressure. However, bronchoconstriction causes airway narrowing, parenchymal distortion, dynamic hyperinflation, and the emergence of ventilation defects (VDefs) affecting transmural pressure. This study aimed to investigate the changes in transmural pressure caused by bronchoconstriction in a bronchial tree. Transmural pressures before and during bronchoconstriction were estimated using an integrative computational model of bronchoconstriction. Briefly, this model incorporates a 12-generation symmetric bronchial tree, and the Anafi and Wilson model for the individual airways of the tree. Bronchoconstriction lead to the emergence of VDefs and a relative increase in peak transmural pressures of up to 84% compared to baseline. The highest increase in peak transmural pressure occurred in a central airway outside of VDefs, and the lowest increase was 27% in an airway within VDefs illustrating the heterogeneity in peak transmural pressures within a bronchial tree. Mechanisms contributing to the increase in peak transmural pressures include increased regional ventilation and dynamic hyperinflation both leading to increased alveolar pressures compared to baseline. Pressure differences between intraluminal and alveolar pressure increased driven by the increased airway resistance and its contribution to total transmural pressure reached up to 24%. In conclusion, peak transmural pressure in lungs with VDefs during bronchoconstriction can be substantially increased compared to dilated airways in healthy homogeneous lungs and is highly heterogeneous. Further insights will depend on the experimental studies taking these conditions into account.

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.