Assessing potential errors of MRI-based measurements of pulmonary blood flow using a detailed network flow model

Vaping causes an acute BMI-dependent change in pulmonary blood flow

Vaping use has skyrocketed especially among young adults, however there is no consensus on how vaping impacts the lungs. We aimed to determine whether there were changes in lung function acutely after a standard vaping session or if there were differences in lung function metrics between a healthy never-vaping cohort (N = 6; 27.3 ± 3.0 years) and a young asymptomatic vaping cohort (N = 14; 26.4 ± 8.0 years) indicating chronic changes. Pulmonary function measurements and impulse oscillometry were obtained on all participants. Oxygen-enhanced and Arterial Spin Labelling MRI were used to measure specific ventilation and perfusion, respectively, before and after vaping, and in the control cohort at baseline. MRI metrics did not show any significant differences in specific ventilation or perfusion after vaping. Heart rate increased post-vaping (68.1 ± 10.5 to 71.3 ± 8.7, p = 0.020); however, this and other metrics did not show a nicotine dose-dependent effect. There was a significant negative correlation between BMI and change in mean perfusion post-vaping (p = 0.003); those with normal/low BMI showing an increase in perfusion and vice versa for high BMI. This may be due to subjects lying supine during vaping inhalation. Pulmonary function metrics indicative of airways resistance showed significant differences between the vaping and control cohorts indicating early airway changes.

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.

The spatial-temporal dynamics of pulmonary blood flow are altered in pulmonary arterial hypertension

Global fluctuation dispersion (FDglobal), a spatial-temporal metric derived from serial images of the pulmonary perfusion obtained with MRI-arterial spin labeling, describes temporal fluctuations in the spatial distribution of perfusion. In healthy subjects, FDglobal is increased by hyperoxia, hypoxia, and inhaled nitric oxide. We evaluated patients with pulmonary arterial hypertension (PAH, 4F, aged 47 ± 15, mean pulmonary artery pressure 48 ± 7 mmHg) and healthy controls (CON, 7F, aged 47 ± 12) to test the hypothesis that FDglobal is increased in PAH. Images were acquired at ∼4-5 s intervals during voluntary respiratory gating, inspected for quality, registered using a deformable registration algorithm, and normalized. Spatial relative dispersion (RD = SD/mean) and the percent of the lung image with no measurable perfusion signal (%NMP) were also assessed. FDglobal was significantly increased in PAH (PAH = 0.40 ± 0.17, CON = 0.17 ± 0.02, P = 0.006, a 135% increase) with no overlap in values between the two groups, consistent with altered vascular regulation. Both spatial RD and %NMP were also markedly greater in PAH vs. CON (PAH RD = 1.46 ± 0.24, CON = 0.90 ± 0.10, P = 0.0004; PAH NMP = 13.4 ± 6.1%; CON = 2.3 ± 1.4%, P = 0.001 respectively) consistent with vascular remodeling resulting in poorly perfused regions of lung and increased spatial heterogeneity. The difference in FDglobal between normal subjects and patients with PAH in this small cohort suggests that spatial-temporal imaging of perfusion may be useful in the evaluation of patients with PAH. Since this MR imaging technique uses no injected contrast agents and has no ionizing radiation it may be suitable for use in diverse patient populations.NEW & NOTEWORTHY Using proton MRI-arterial spin labeling to obtain serial images of pulmonary perfusion, we show that global fluctuation dispersion (FDglobal), a metric of temporal fluctuations in the spatial distribution of perfusion, was significantly increased in female patients with pulmonary arterial hypertension (PAH) compared with healthy controls. This potentially indicates pulmonary vascular dysregulation. Dynamic measures using proton MRI may provide new tools for evaluating individuals at risk of PAH or for monitoring therapy in patients with PAH.

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.

Abnormal pulmonary perfusion heterogeneity in patients with Fontan circulation and pulmonary arterial hypertension

KEY POINTS

The distribution of pulmonary perfusion is affected by gravity, vascular branching structure and active regulatory mechanisms, which may be disrupted by cardiopulmonary disease, but this is not well studied, particularly in rare conditions. We evaluated pulmonary perfusion in patients who had undergone Fontan procedure, patients with pulmonary arterial hypertension (PAH) and two groups of controls using a proton magnetic resonance imaging technique, arterial spin labelling to measure perfusion. Heterogeneity was assessed by the relative dispersion (SD/mean) and gravitational gradients. Gravitational gradients were similar between all groups, but heterogeneity was significantly increased in both patient groups compared to controls and persisted after removing contributions from large blood vessels and gravitational gradients. Patients with Fontan physiology and patients with PAH have increased pulmonary perfusion heterogeneity that is not explainable by differences in mean perfusion, gravitational gradients, or large vessel anatomy. This probably reflects vascular remodelling in PAH and possibly in Fontan physiology.

ABSTRACT

Many factors affect the distribution of pulmonary perfusion, which may be disrupted by cardiopulmonary disease, but this is not well studied, particularly in rare conditions. An example is following the Fontan procedure, where pulmonary perfusion is passive, and heterogeneity may be increased because of the underlying pathophysiology leading to Fontan palliation, remodelling, or increased gravitational gradients from low flow. Another is pulmonary arterial hypertension (PAH), where gravitational gradients may be reduced secondary to high pressures, but remodelling may increase perfusion heterogeneity. We evaluated regional pulmonary perfusion in Fontan patients (n = 5), healthy young controls (Fontan control, n = 5), patients with PAH (n = 6) and healthy older controls (PAH control) using proton magnetic resonance imaging. Regional perfusion was measured using arterial spin labelling. Heterogeneity was assessed by the relative dispersion (SD/mean) and gravitational gradients. Mean perfusion was similar (Fontan = 2.50 ± 1.02 ml min^-1  ml^-1 ; Fontan control = 3.09 ± 0.58, PAH = 3.63 ± 1.95; PAH control = 3.98 ± 0.91, P = 0.26), and the slopes of gravitational gradients were not different (Fontan = -0.23 ± 0.09 ml min^-1  ml^-1  cm^-1 ; Fontan control = -0.29 ± 0.23, PAH = -0.27 ± 0.09, PAH control = -0.25 ± 0.18, P = 0.91) between groups. Perfusion relative dispersion was greater in both Fontan and PAH than controls (Fontan = 1.46 ± 0.18; Fontan control = 0.99 ± 0.21, P = 0.005; PAH = 1.22 ± 0.27, PAH control = 0.91 ± 0.12, P = 0.02) but similar between patient groups (P = 0.13). These findings persisted after removing contributions from large blood vessels and gravitational gradients (all P < 0.05). We conclude that patients with Fontan physiology and PAH have increased pulmonary perfusion heterogeneity that is not explained by differences in mean perfusion, gravitational gradients, or large vessel anatomy. This probably reflects the effects of remodelling in PAH and possibly in Fontan physiology.

Computational modeling of the interactions between the maternal and fetal circulations in human pregnancy

In pregnancy, fetal growth is supported by its placenta. In turn, the placenta is nourished by maternal blood, delivered from the uterus, in which the vasculature is dramatically transformed to deliver this blood an ever increasing volume throughout gestation. A healthy pregnancy is thus dependent on the development of both the placental and maternal circulations, but also the interface where these physically separate circulations come in close proximity to exchange gases and nutrients between mum and baby. As the system continually evolves during pregnancy, our understanding of normal vascular anatomy, and how this impacts placental exchange function is limited. Understanding this is key to improve our ability to understand, predict, and detect pregnancy pathologies, but presents a number of challenges, due to the inaccessibility of the pregnant uterus to invasive measurements, and limitations in the resolution of imaging modalities safe for use in pregnancy. Computational approaches provide an opportunity to gain new insights into normal and abnormal pregnancy, by connecting observed anatomical changes from high-resolution imaging to function, and providing metrics that can be observed by routine clinical ultrasound. Such advanced modeling brings with it challenges to scale detailed anatomical models to reflect organ level function. This suggests pathways for future research to provide models that provide both physiological insights into pregnancy health, but also are simple enough to guide clinical focus. We the review evolution of computational approaches to understanding the physiology and pathophysiology of pregnancy in the uterus, placenta, and beyond focusing on both opportunities and challenges. This article is categorized under: Reproductive System Diseases >Computational Models.

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.

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 Po₂ and Pco₂. 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 CO₂ added to inspired gas during the challenge block to effect a 7-Torr increase in end-tidal CO₂. 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 CO₂ 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.

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̇o2max = 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 (LogSDV) and perfusion (LogSDQ) 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. LogSDV 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, LogSDQ 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.

Optimizing human pulmonary perfusion measurement using an in silico model of arterial spin labeling magnetic resonance imaging

Arterial spin labeling (ASL) magnetic resonance imaging (MRI) is an imaging methodology that uses blood as an endogenous contrast agent to quantify flow. One limitation of this method of capillary blood quantification when applied in the lung is the contribution of signals from non-capillary blood. Intensity thresholding is one approach that has been proposed for minimizing the non-capillary blood signal. This method has been tested in previous in silico modeling studies; however, it has only been tested under a restricted set of physiological conditions (supine posture and a cardiac output of 5 L/min). This study presents an in silico approach that extends previous intensity thresholding analysis to estimate the optimal "per-slice" intensity threshold value using the individual components of the simulated ASL signal (signal arising independently from capillary blood as well as pulmonary arterial and pulmonary venous blood). The aim of this study was to assess whether the threshold value should vary with slice location, posture, or cardiac output. We applied an in silico modeling approach to predict the blood flow distribution and the corresponding ASL quantification of pulmonary perfusion in multiple sagittal imaging slices. There was a significant increase in ASL signal and heterogeneity (COV = 0.90 to COV = 1.65) of ASL signals when slice location changed from lateral to medial. Heterogeneity of the ASL signal within a slice was significantly lower (P = 0.03) in prone (COV = 1.08) compared to in the supine posture (COV = 1.17). Increasing stroke volume resulted in an increase in ASL signal and conversely an increase in heart rate resulted in a decrease in ASL signal. However, when cardiac output was increased via an increase in both stroke volume and heart rate, ASL signal remained relatively constant. Despite these differences, we conclude that a threshold value of 35% provides optimal removal of large vessel signal independent of slice location, posture, and cardiac output.

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.

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.

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.

A statistical clustering approach to discriminating perfusion from conduit vessel signal contributions in a pulmonary ASL MR image

The measurement of pulmonary perfusion (blood delivered to the capillary bed within a voxel) using arterial spin labeling (ASL) magnetic resonance imaging is often complicated by signal artifacts from conduit vessels that carry blood destined for voxels at a distant location in the lung. One approach to dealing with conduit vessel contributions involves the application of an absolute threshold on the ASL signal. While useful for identifying a subset of the most dominant high signal conduit image features, signal thresholding cannot discriminate between perfusion and conduit vessel contributions at intermediate and low signal. As an alternative, this article discusses a data-driven statistical approach based on statistical clustering for characterizing and discriminating between capillary perfusion and conduit vessel contributions over the full signal spectrum. An ASL flow image is constructed from the difference between a pair of tagged magnetic resonance images. However, when viewed as a bivariate projection that treats the image pair as independent measures (rather than the univariate quantity that results from the subtraction of the two images), the signal associated with capillary perfusion contributions is observed to cluster independently of the signal associated with conduit vessel contributions. Analyzing the observed clusters using a Gaussian mixture model makes it possible to discriminate between conduit vessel and capillary-perfusion-dominated signal contributions over the full signal spectrum of the ASL image. As a demonstration of feasibility, this study compares the proposed clustering approach with the standard absolute signal threshold strategy in a small number of test images.

Advances in functional and structural imaging of the human lung using proton MRI

The field of proton lung MRI is advancing on a variety of fronts. In the realm of functional imaging, it is now possible to use arterial spin labeling (ASL) and oxygen-enhanced imaging techniques to quantify regional perfusion and ventilation, respectively, in standard units of measurement. By combining these techniques into a single scan, it is also possible to quantify the local ventilation-perfusion ratio, which is the most important determinant of gas-exchange efficiency in the lung. To demonstrate potential for accurate and meaningful measurements of lung function, this technique was used to study gravitational gradients of ventilation, perfusion, and ventilation-perfusion ratio in healthy subjects, yielding quantitative results consistent with expected regional variations. Such techniques can also be applied in the time domain, providing new tools for studying temporal dynamics of lung function. Temporal ASL measurements showed increased spatial-temporal heterogeneity of pulmonary blood flow in healthy subjects exposed to hypoxia, suggesting sensitivity to active control mechanisms such as hypoxic pulmonary vasoconstriction, and illustrating that to fully examine the factors that govern lung function it is necessary to consider temporal as well as spatial variability. Further development to increase spatial coverage and improve robustness would enhance the clinical applicability of these new functional imaging tools. In the realm of structural imaging, pulse sequence techniques such as ultrashort echo-time radial k-space acquisition, ultrafast steady-state free precession, and imaging-based diaphragm triggering can be combined to overcome the significant challenges associated with proton MRI in the lung, enabling high-quality three-dimensional imaging of the whole lung in a clinically reasonable scan time. Images of healthy and cystic fibrosis subjects using these techniques demonstrate substantial promise for non-contrast pulmonary angiography and detailed depiction of airway disease. Although there is opportunity for further optimization, such approaches to structural lung imaging are ready for clinical testing.

Inhaled nitric oxide alters the distribution of blood flow in the healthy human lung, suggesting active hypoxic pulmonary vasoconstriction in normoxia

Hypoxic pulmonary vasoconstriction (HPV) is thought to actively regulate ventilation-perfusion (V̇a/Q̇) matching, reducing perfusion in regions of alveolar hypoxia. We assessed the extent of HPV in the healthy human lung using inhaled nitric oxide (iNO) under inspired oxygen fractions (FiO2 ) of 0.125, 0.21, and 0.30 (a hyperoxic stimulus designed to abolish HPV without the development of atelectasis). Dynamic measures of blood flow were made in a single sagittal slice of the right lung of five healthy male subjects using an arterial spin labeling (ASL) MRI sequence, following a block stimulus pattern (3 × 60 breaths) with 40 ppm iNO administered in the central block. The overall spatial heterogeneity, spatiotemporal variability, and regional pattern of pulmonary blood flow was quantified as a function of condition (FiO2 × iNO state). While spatial heterogeneity did not change significantly with iNO administration or FiO2 , there were statistically significant increases in Global Fluctuation Dispersion, (a marker of spatiotemporal flow variability) when iNO was administered during hypoxia (5.4 percentage point increase, P = 0.003). iNO had an effect on regional blood flow that was FiO2 dependent (P = 0.02), with regional changes in the pattern of blood flow occurring in hypoxia (P = 0.007) and normoxia (P = 0.008) tending to increase flow to dependent lung at the expense of nondependent lung. These findings indicate that inhaled nitric oxide significantly alters the distribution of blood flow in both hypoxic and normoxic healthy subjects, and suggests that some baseline HPV may indeed be present in the normoxic lung.

Multi-scale computational models of the airways to unravel the pathophysiological mechanisms in asthma and chronic obstructive pulmonary disease (AirPROM)

THE RESPIRATORY SYSTEM COMPRISES SEVERAL SCALES OF BIOLOGICAL COMPLEXITY: the genes, cells and tissues that work in concert to generate resultant function. Malfunctions of the structure or function of components at any spatial scale can result in diseases, to the detriment of gas exchange, right heart function and patient quality of life. Vast amounts of data emerge from studies across each of the biological scales; however, the question remains: how can we integrate and interpret these data in a meaningful way? Respiratory disease presents a huge health and economic burden, with the diseases asthma and chronic obstructive pulmonary disease (COPD) affecting over 500 million people worldwide. Current therapies are inadequate owing to our incomplete understanding of the disease pathophysiology and our lack of recognition of the enormous disease heterogeneity: we need to characterize this heterogeneity on a patient-specific basis to advance healthcare. In an effort to achieve this goal, the AirPROM consortium (Airway disease Predicting Outcomes through patient-specific computational Modelling) brings together a multi-disciplinary team and a wealth of clinical data. Together we are developing an integrated multi-scale model of the airways in order to unravel the complex pathophysiological mechanisms occurring in the diseases asthma and COPD.

The effect of supine exercise on the distribution of regional pulmonary blood flow measured using proton MRI

The Zone model of pulmonary perfusion predicts that exercise reduces perfusion heterogeneity because increased vascular pressure redistributes flow to gravitationally nondependent lung, and causes dilation and recruitment of blood vessels. However, during exercise in animals, perfusion heterogeneity as measured by the relative dispersion (RD, SD/mean) is not significantly decreased. We evaluated the effect of exercise on pulmonary perfusion in six healthy supine humans using magnetic resonance imaging (MRI). Data were acquired at rest, while exercising (∼27% of maximal oxygen consumption) using a MRI-compatible ergometer, and in recovery. Images were acquired in most of the right lung in the sagittal plane at functional residual capacity, using a 1.5-T MR scanner equipped with a torso coil. Perfusion was measured using arterial spin labeling (ASL-FAIRER) and regional proton density using a fast multiecho gradient-echo sequence. Perfusion images were corrected for coil-based signal heterogeneity, large conduit vessels removed and quantified (in ml·min(-1)·ml(-1)) (perfusion), and also normalized for density and quantified (in ml·min(-1)·g(-1)) (density-normalized perfusion, DNP) accounting for tissue redistribution. DNP increased during exercise (11.1 ± 3.5 rest, 18.8 ± 2.3 exercise, 13.2 ± 2.2 recovery, ml·min(-1)·g(-1), P < 0.0001), and the increase was largest in nondependent lung (110 ± 61% increase in nondependent, 63 ± 35% in mid, 70 ± 33% in dependent, P < 0.005). The RD of perfusion decreased with exercise (0.93 ± 0.21 rest, 0.73 ± 0.13 exercise, 0.94 ± 0.18 recovery, P < 0.005). The RD of DNP showed a similar trend (0.82 ± 0.14 rest, 0.75 ± 0.09 exercise, 0.81 ± 0.10 recovery, P = 0.13). In conclusion, in contrast to animal studies, in supine humans, mild exercise decreased perfusion heterogeneity, consistent with Zone model predictions.

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.

The gravitational distribution of ventilation-perfusion ratio is more uniform in prone than supine posture in the normal human lung

The gravitational gradient of intrapleural pressure is suggested to be less in prone posture than supine. Thus the gravitational distribution of ventilation is expected to be more uniform prone, potentially affecting regional ventilation-perfusion (Va/Q) ratio. Using a novel functional lung magnetic resonance imaging technique to measure regional Va/Q ratio, the gravitational gradients in proton density, ventilation, perfusion, and Va/Q ratio were measured in prone and supine posture. Data were acquired in seven healthy subjects in a single sagittal slice of the right lung at functional residual capacity. Regional specific ventilation images quantified using specific ventilation imaging and proton density images obtained using a fast gradient-echo sequence were registered and smoothed to calculate regional alveolar ventilation. Perfusion was measured using arterial spin labeling. Ventilation (ml·min(-1)·ml(-1)) images were combined on a voxel-by-voxel basis with smoothed perfusion (ml·min(-1)·ml(-1)) images to obtain regional Va/Q ratio. Data were averaged for voxels within 1-cm gravitational planes, starting from the most gravitationally dependent lung. The slope of the relationship between alveolar ventilation and vertical height was less prone than supine (-0.17 ± 0.10 ml·min(-1)·ml(-1)·cm(-1) supine, -0.040 ± 0.03 prone ml·min(-1)·ml(-1)·cm(-1), P = 0.02) as was the slope of the perfusion-height relationship (-0.14 ± 0.05 ml·min(-1)·ml(-1)·cm(-1) supine, -0.08 ± 0.09 prone ml·min(-1)·ml(-1)·cm(-1), P = 0.02). There was a significant gravitational gradient in Va/Q ratio in both postures (P < 0.05) that was less in prone (0.09 ± 0.08 cm(-1) supine, 0.04 ± 0.03 cm(-1) prone, P = 0.04). The gravitational gradients in ventilation, perfusion, and regional Va/Q ratio were greater supine than prone, suggesting an interplay between thoracic cavity configuration, airway and vascular tree anatomy, and the effects of gravity on Va/Q matching.

Spatial-temporal dynamics of pulmonary blood flow in the healthy human lung in response to altered FI(O2)

The temporal dynamics of blood flow in the human lung have been largely unexplored due to the lack of appropriate technology. Using the magnetic resonance imaging method of arterial spin labeling (ASL) with subject-gated breathing, we produced a dynamic series of flow-weighted images in a single sagittal slice of the right lung with a spatial resolution of ~1 cm(3) and a temporal resolution of ~10 s. The mean flow pattern determined from a set of reference images was removed to produce a time series of blood flow fluctuations. The fluctuation dispersion (FD), defined as the spatial standard deviation of each flow fluctuation map, was used to quantify the changes in distribution of flow in six healthy subjects in response to 100 breaths of hypoxia (FI(O(2)) = 0.125) or hyperoxia (FI(O(2)) = 1.0). Two reference frames were used in calculation, one determined from the initial set of images (FD(global)), and one determined from the mean of each corresponding baseline or challenge period (FD(local)). FD(local) thus represented changes in temporal variability as a result of intervention, whereas FD(global) encompasses both FD(local) and any generalized redistribution of flow associated with switching between two steady-state patterns. Hypoxic challenge resulted in a significant increase (96%, P < 0.001) in FD(global) from the normoxic control period and in FD(local) (46%, P = 0.0048), but there was no corresponding increase in spatial relative dispersion (spatial standard deviation of the images divided by the mean; 8%, not significant). There was a smaller increase in FD(global) in response to hyperoxia (47%, P = 0.0015) for the single slice, suggestive of a more general response of the pulmonary circulation to a change from normoxia to hyperoxia. These results clearly demonstrate a temporal change in the sampled distribution of pulmonary blood flow in response to hypoxia, which is not observed when considering only the relative dispersion of the spatial distribution.

Imaging lung perfusion

From the first measurements of the distribution of pulmonary blood flow using radioactive tracers by West and colleagues (J Clin Invest 40: 1-12, 1961) allowing gravitational differences in pulmonary blood flow to be described, the imaging of pulmonary blood flow has made considerable progress. The researcher employing modern imaging techniques now has the choice of several techniques, including magnetic resonance imaging (MRI), computerized tomography (CT), positron emission tomography (PET), and single photon emission computed tomography (SPECT). These techniques differ in several important ways: the resolution of the measurement, the type of contrast or tag used to image flow, and the amount of ionizing radiation associated with each measurement. In addition, the techniques vary in what is actually measured, whether it is capillary perfusion such as with PET and SPECT, or larger vessel information in addition to capillary perfusion such as with MRI and CT. Combined, these issues affect quantification and interpretation of data as well as the type of experiments possible using different techniques. The goal of this review is to give an overview of the techniques most commonly in use for physiological experiments along with the issues unique to each technique.

Imaging for lung physiology: what do we wish we could measure?

The role of imaging as a tool for investigating lung physiology is growing at an accelerating pace. Looking forward, we wished to identify unresolved issues in lung physiology that might realistically be addressed by imaging methods in development or imaging approaches that could be considered. The role of imaging is framed in terms of the importance of good spatial and temporal resolution and the types of questions that could be addressed as these technical capabilities improve. Recognizing that physiology is fundamentally a quantitative science, a recurring emphasis is on the need for imaging methods that provide reliable measurements of specific physiological parameters. The topics included necessarily reflect our perspective on what are interesting questions and are not meant to be a comprehensive review. Nevertheless, we hope that this essay will be a spur to physiologists to think about how imaging could usefully be applied in their research and to physical scientists developing new imaging methods to attack challenging questions imaging could potentially answer.