Hypoxic pulmonary vasoconstriction does not contribute to pulmonary blood flow heterogeneity in normoxia in normal supine humans

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 Pulmonary Vasculature

The pulmonary circulation is a low-pressure, low-resistance circuit whose primary function is to deliver deoxygenated blood to, and oxygenated blood from, the pulmonary capillary bed enabling gas exchange. The distribution of pulmonary blood flow is regulated by several factors including effects of vascular branching structure, large-scale forces related to gravity, and finer scale factors related to local control. Hypoxic pulmonary vasoconstriction is one such important regulatory mechanism. In the face of local hypoxia, vascular smooth muscle constriction of precapillary arterioles increases local resistance by up to 250%. This has the effect of diverting blood toward better oxygenated regions of the lung and optimizing ventilation-perfusion matching. However, in the face of global hypoxia, the net effect is an increase in pulmonary arterial pressure and vascular resistance. Pulmonary vascular resistance describes the flow-resistive properties of the pulmonary circulation and arises from both precapillary and postcapillary resistances. The pulmonary circulation is also distensible in response to an increase in transmural pressure and this distention, in addition to recruitment, moderates pulmonary arterial pressure and vascular resistance. This article reviews the physiology of the pulmonary vasculature and briefly discusses how this physiology is altered by common circumstances.

Gas Exchange in the Lung

Gas exchange in the lung depends on tidal breathing, which brings new oxygen to and removes carbon dioxide from alveolar gas. This maintains alveolar partial pressures that promote passive diffusion to add oxygen and remove carbon dioxide from blood in alveolar capillaries. In a lung model without ventilation and perfusion (V̇AQ̇) mismatch, alveolar partial pressures of oxygen and carbon dioxide are primarily determined by inspiratory pressures and alveolar ventilation. Regions with shunt or low ratios worsen arterial oxygenation while alveolar dead space and high lung units lessen CO2 elimination efficiency. Although less common, diffusion limitation might cause hypoxemia in some situations. This review covers the principles of lung gas exchange and therefore mechanisms of hypoxemia or hypercapnia. In addition, we discuss different metrics that quantify the deviation from ideal gas exchange.

Impact of Vasodilation on Oxygen-Enhanced Functional Lung MRI at 0.55 T

BACKGROUND

Oxygen-enhanced magnetic resonance imaging (OE-MRI) can be used to assess regional lung function without ionizing radiation. Inhaled oxygen acts as a T1-shortening contrast agent to increase signal in T1-weighted (T1w) images. However, increase in proton density from pulmonary hyperoxic vasodilation may also contribute to the measured signal enhancement. Our aim was to quantify the relative contributions of the T1-shortening and vasodilatory effects of oxygen to signal enhancement in OE-MRI in both swine and healthy volunteers.

METHODS

We imaged 14 anesthetized female swine (47 ± 8 kg) using a prototype 0.55 T high-performance MRI system while experimentally manipulating oxygenation and blood volume independently through oxygen titration, partial occlusion of the vena cava for volume reduction, and infusion of colloid fluid (6% hydroxyethyl starch) for volume increase. Ten healthy volunteers were imaged before, during, and after hyperoxia. Two proton density-weighted (PDw) and 2 T1w ultrashort echo time images were acquired per experimental state. The median PDw and T1w percent signal enhancement (PSE), compared with baseline room air, was calculated after image registration and correction for lung volume changes. Differences in median PSE were compared using Wilcoxon signed rank test.

RESULTS

The PSE in PDw images after 100% oxygen was similar in swine (1.66% ± 1.41%, P = 0.01) and in healthy volunteers (1.99% ± 1.79%, P = 0.02), indicating that oxygen-induced pulmonary vasodilation causes ~2% lung proton density increase. The PSE in T1w images after 100% oxygen was also similar (swine, 9.20% ± 1.68%, P < 0.001; healthy volunteers, 10.10% ± 3.05%, P < 0.001). The PSE in T1w enhancement was oxygen dose-dependent in anesthetized swine, and we measured a dose-dependent PDw image signal increase from infused fluids.

CONCLUSIONS

The contribution of oxygen-induced vasodilation to T1w OE-MRI signal was measurable using PDw imaging and was found to be ~2% in both anesthetized swine and in healthy volunteers. This finding may have implications for patients with regional or global hypoxia or vascular dysfunction undergoing OE-MRI and suggest that PDw imaging may be useful to account for oxygen-induced vasodilation in OE-MRI.

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.

MR imaging of the airways

The need for airway imaging is defined by the limited sensitivity of common clinical tests like spirometry, lung diffusion (DLCO) and blood gas analysis to early changes of peripheral airways and to inhomogeneous regional distribution of lung function deficits. Therefore, X-ray and computed tomography (CT) are frequently used to complement the standard tests.As an alternative, magnetic resonance imaging (MRI) offers radiation-free lung imaging, but at lower spatial resolution. Non-contrast enhanced MRI shows healthy airways down to the first subsegmental level/4^th order (CT: eighth). Bronchiectasis can be identified by wall thickening and fluid accumulation. Smaller airways become visible, when altered by peribronchiolar inflammation or mucus retention (tree-in-bud sign).The strength of MRI is functional imaging. Dynamic, time-resolved MRI directly visualizes expiratory airway collapse down to the lobar level (CT: segmental level). Obstruction of even smaller airways becomes visible as air trapping on the expiratory scans. MRI with hyperpolarized noble gases (³He, ¹²⁹Xe) directly shows the large airways and peripheral lung ventilation. Dynamic contrast-enhanced MRI (DCE MRI) indirectly shows airway dysfunction as perfusion deficits resulting from hypoxic vasoconstriction of the dependent lung volumes. Further promising scientific approaches such as non-contrast enhanced, ventilation-/perfusion-weighted MRI from periodic signal changes of respiration and blood flow are in development.In summary, MRI of the lungs and airways excels with its unique combination of morphologic and functional imaging capacities for research (e.g., in chronic obstructive lung disease or asthma) as well as for clinical imaging (e.g., in cystic fibrosis).

Perfusion imaging heterogeneity during NO inhalation distinguishes pulmonary arterial hypertension (PAH) from healthy subjects and has potential as an imaging biomarker

BACKGROUND

Without aggressive treatment, pulmonary arterial hypertension (PAH) has a 5-year mortality of approximately 40%. A patient's response to vasodilators at diagnosis impacts the therapeutic options and prognosis. We hypothesized that analyzing perfusion images acquired before and during vasodilation could identify characteristic differences between PAH and control subjects.

METHODS

We studied 5 controls and 4 subjects with PAH using HRCT and ¹³NN PET imaging of pulmonary perfusion and ventilation. The total spatial heterogeneity of perfusion (CV²_Qtotal) and its components in the vertical (CV²_Qvgrad) and cranio-caudal (CV²_Qzgrad) directions, and the residual heterogeneity (CV²_Qr), were assessed at baseline and while breathing oxygen and nitric oxide (O₂ + iNO). The length scale spectrum of CV²_Qr was determined from 10 to 110 mm, and the response of regional perfusion to O₂ + iNO was calculated as the mean of absolute differences. Vertical gradients in perfusion (Q_vgrad) were derived from perfusion images, and ventilation-perfusion distributions from images of ¹³NN washout kinetics.

RESULTS

O₂ + iNO significantly enhanced perfusion distribution differences between PAH and controls, allowing differentiation of PAH subjects from controls. During O₂ + iNO, CV²_Qvgrad was significantly higher in controls than in PAH (0.08 (0.055-0.10) vs. 6.7 × 10^-3 (2 × 10^-4-0.02), p < 0.001) with a considerable gap between groups. Q_vgrad and CV²_Qtotal showed smaller differences: - 7.3 vs. - 2.5, p = 0.002, and 0.12 vs. 0.06, p = 0.01. CV²_Qvgrad had the largest effect size among the primary parameters during O₂ + iNO. CV²_Qr, and its length scale spectrum were similar in PAH and controls. Ventilation-perfusion distributions showed a trend towards a difference between PAH and controls at baseline, but it was not statistically significant.

CONCLUSIONS

Perfusion imaging during O2 + iNO showed a significant difference in the heterogeneity associated with the vertical gradient in perfusion, distinguishing in this small cohort study PAH subjects from controls.

The hypoxic pulmonary vasoconstriction: From physiology to clinical application in thoracic surgery

More than 70 years after its original report, the hypoxic pulmonary vasoconstriction (HPV) response continues to spark scientific interest on its mechanisms and clinical implications, particularly for anesthesiologists involved in thoracic surgery. Selective airway intubation and one-lung ventilation (OLV) facilitates the surgical intervention on a collapsed lung while the HPV redirects blood flow from the "upper" non-ventilated hypoxic lung to the "dependent" ventilated lung. Therefore, by limiting intrapulmonary shunting and optimizing ventilation-to-perfusion (V/Q) ratio, the fall in arterial oxygen pressure (PaO₂) is attenuated during OLV. The HPV involves a biphasic response mobilizing calcium within pulmonary vascular smooth muscles, which is activated within seconds after exposure to low alveolar oxygen pressure and that gradually disappears upon re-oxygenation. Many factors including acid-base balance, the degree of lung expansion, circulatory volemia as well as lung diseases and patient age affect HPV. Anesthetic agents, analgesics and cardiovascular medications may also interfer with HPV during the perioperative period. Since HPV represents the homeostatic mechanism for regional ventilation-to-perfusion matching and in turn, for optimal pulmonary oxygen uptake, a clear understanding of HPV is clinically relevant for all anesthesiologists.

Hypoxic pulmonary vasoconstriction as a regulator of alveolar-capillary oxygen flux: A computational model of ventilation-perfusion matching

The relationship between regional variabilities in airflow (ventilation) and blood flow (perfusion) is a critical determinant of gas exchange efficiency in the lungs. Hypoxic pulmonary vasoconstriction is understood to be the primary active regulator of ventilation-perfusion matching, where upstream arterioles constrict to direct blood flow away from areas that have low oxygen supply. However, it is not understood how the integrated action of hypoxic pulmonary vasoconstriction affects oxygen transport at the system level. In this study we develop, and make functional predictions with a multi-scale multi-physics model of ventilation-perfusion matching governed by the mechanism of hypoxic pulmonary vasoconstriction. Our model consists of (a) morphometrically realistic 2D pulmonary vascular networks to the level of large arterioles and venules; (b) a tileable lumped-parameter model of vascular fluid and wall mechanics that accounts for the influence of alveolar pressure; (c) oxygen transport accounting for oxygen bound to hemoglobin and dissolved in plasma; and (d) a novel empirical model of hypoxic pulmonary vasoconstriction. Our model simulations predict that under the artificial test condition of a uniform ventilation distribution (1) hypoxic pulmonary vasoconstriction matches perfusion to ventilation; (2) hypoxic pulmonary vasoconstriction homogenizes regional alveolar-capillary oxygen flux; and (3) hypoxic pulmonary vasoconstriction increases whole-lobe oxygen uptake by improving ventilation-perfusion matching.

The relationship between regional ventilation (airflow) and perfusion (blood flow) is a major determinant of gas exchange efficiency. Atelactasis and pulmonary vascular occlusive diseases, such as acute pulmonary embolism, are characterized by ventilation-perfusion mismatching and decreased oxygen in the bloodstream. Despite the physiological and medical importance of ventilation-perfusion matching, there are gaps in our knowledge of the regulatory mechanisms that maintain adequate gas exchange under pathological and normal conditions. Hypoxic pulmonary vasoconstriction is understood to be the primary regulator of ventilation-perfusion matching, where upstream arterioles constrict to direct blood flow away from areas that have low oxygen supply, yet it is not understood how this mechanism affects oxygen transport at the system level. In this study we present a computational model of the ventilation-perfusion matching and hypoxic pulmonary vasoconstriction to better understand how physiological regulation at the regional level scales to affect oxygen transport at the system level. Our model simulations predict that this regulatory mechanism improves the spatial overlap of airflow and blood flow, which serves to increase the uptake of oxygen into the bloodstream. This improved understanding of ventilation-perfusion matching may offer insights into the etiology of, and therapeutic interventions for diseases characterized by ventilation-perfusion mismatching.

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.

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.

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.

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.

Magnetic resonance angiography for the primary diagnosis of pulmonary embolism: A review from the international workshop for pulmonary functional imaging

Pulmonary contrast enhanced magnetic resonance angiography (CE-MRA) is useful for the primary diagnosis of pulmonary embolism (PE). Many sites have chosen not to use CE-MRA as a first line of diagnostic tool for PE because of the speed and higher efficacy of computerized tomographic angiography (CTA). In this review, we discuss the strengths and weaknesses of CE-MRA and the appropriate imaging scenarios for the primary diagnosis of PE derived from our unique multi-institutional experience in this area. The optimal patient for this test has a low to intermediate suspicion for PE based on clinical decision rules. Patients in extremis are not candidates for this test. Younger women (< 35 years of age) and patients with iodinated contrast allergies are best served by using this modality We discuss the history of the use of this test, recent technical innovations, artifacts, direct and indirect findings for PE, ancillary findings, and the effectiveness (patient outcomes) of CE-MRA for the exclusion of PE. Current outcomes data shows that CE-MRA and NM V/Q scans are effective alternative tests to CTA for the primary diagnosis of PE.

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.

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.

Iron, oxygen, and the pulmonary circulation

The human pulmonary vasculature vasoconstricts in response to a reduction in alveolar oxygen tension, a phenomenon termed hypoxic pulmonary vasoconstriction (HPV). This review describes the time course of this behavior, which occurs in distinct phases, and then explores the importance for HPV of the hypoxia-inducible factor (HIF) pathway. Next, the HIF-hydroxylase enzymes that act as molecular oxygen sensors within the HIF pathway are discussed. These enzymes are particularly sensitive to intracellular iron availability, which confers iron-sensing properties on the HIF pathway. Human studies of iron chelation and supplementation are then reviewed. These demonstrate that the iron sensitivity of the HIF pathway evident from in vitro experiments is relevant to human pulmonary vascular physiology. Next, the importance of iron status in high-altitude illness and chronic cardiopulmonary disease is explored, and the therapeutic potential of intravenous iron discussed. The review concludes by highlighting some further complexities that arise from interactions between the HIF pathway and other intracellular iron-sensing mechanisms.

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.

Gas exchange and ventilation-perfusion relationships in the lung

This review provides an overview of the relationship between ventilation/perfusion ratios and gas exchange in the lung, emphasising basic concepts and relating them to clinical scenarios. For each gas exchanging unit, the alveolar and effluent blood partial pressures of oxygen and carbon dioxide (PO2 and PCO2) are determined by the ratio of alveolar ventilation to blood flow (V'A/Q') for each unit. Shunt and low V'A/Q' regions are two examples of V'A/Q' mismatch and are the most frequent causes of hypoxaemia. Diffusion limitation, hypoventilation and low inspired PO2 cause hypoxaemia, even in the absence of V'A/Q' mismatch. In contrast to other causes, hypoxaemia due to shunt responds poorly to supplemental oxygen. Gas exchanging units with little or no blood flow (high V'A/Q' regions) result in alveolar dead space and increased wasted ventilation, i.e. less efficient carbon dioxide removal. Because of the respiratory drive to maintain a normal arterial PCO2, the most frequent result of wasted ventilation is increased minute ventilation and work of breathing, not hypercapnia. Calculations of alveolar-arterial oxygen tension difference, venous admixture and wasted ventilation provide quantitative estimates of the effect of V'A/Q' mismatch on gas exchange. The types of V'A/Q' mismatch causing impaired gas exchange vary characteristically with different lung diseases.

Hypoxic pulmonary vasoconstriction as a contributor to response in acute pulmonary embolism

Hypoxic pulmonary vasoconstriction (HPV) is an adaptive response unique to the lung whereby blood flow is diverted away from areas of low alveolar oxygen to improve ventilation-perfusion matching and resultant gas exchange. Some previous experimental studies have suggested that the HPV response to hypoxia is blunted in acute pulmonary embolism (APE), while others have concluded that HPV contributes to elevated pulmonary blood pressures in APE. To understand these contradictory observations, we have used a structure-based computational model of integrated lung function in 10 subjects to study the impact of HPV on pulmonary hemodynamics and gas exchange in the presence of regional arterial occlusion. The integrated model includes an experimentally-derived model for HPV. Its function is validated against measurements of pulmonary vascular resistance in normal subjects at four levels of inspired oxygen. Our results show that the apparently disparate observations of previous studies can be explained within a single model: the model predicts that HPV increases mean pulmonary artery pressure in APE (by 8.2 ± 7.0% in these subjects), and concurrently shows a reduction in response to hypoxia in the subjects who have high levels of occlusion and therefore maximal HPV in normoxia.

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.

Gas exchange consequences of left heart failure

This review explores the pathophysiology of gas exchange abnormalities arising consequent to either acute or chronic elevation of pulmonary venous pressures. The initial experimental studies of acute pulmonary edema outlined the sequence of events from lymphatic congestion with edema fluid to frank alveolar flooding and its resultant hypoxemia. Clinical studies of acute heart failure (HF) suggested that hypoxemia was associated only with the final stage of alveolar flooding. However, in patients with chronic heart failure and normal oxygenation, hypoxemia could be produced by the administration of potent pulmonary vasodilators, suggesting that hypoxic pulmonary vasoconstriction is an important reflex for these patients. Patients with chronic left HF commonly manifest a reduced diffusing capacity, an abnormality that appears to be a consequence of chronic elevation of left atrial pressure. That reduction in diffusing capacity does not appear to be primarily attributable to increases in lung water but is improved by any sustained treatment that improves overall cardiac function. Patients with heart failure may also manifest an abnormally elevated VE/VCO2 during exercise, and that exercise ventilation abnormality arises as a consequence of both alveolar hyperventilation and elevated physiologic dead space. That elevated exercise VE/VCO2 in an HF patient has proven to be a powerful predictor of an adverse outcome and hence it has received sustained attention in the HF literature. At least three of the classes of drugs used to treat HF will normalize the exercise VE/VCO2, suggesting that the excessive ventilation response may be linked to elevated sympathetic activity.

Determinants of pulmonary blood flow distribution

The primary function of the pulmonary circulation is to deliver blood to the alveolar capillaries to exchange gases. Distributing blood over a vast surface area facilitates gas exchange, yet the pulmonary vascular tree must be constrained to fit within the thoracic cavity. In addition, pressures must remain low within the circulatory system to protect the thin alveolar capillary membranes that allow efficient gas exchange. The pulmonary circulation is engineered for these unique requirements and in turn these special attributes affect the spatial distribution of blood flow. As the largest organ in the body, the physical characteristics of the lung vary regionally, influencing the spatial distribution on large-, moderate-, and small-scale levels.

Distribution of perfusion

Local driving pressures and resistances within the pulmonary vascular tree determine the distribution of perfusion in the lung. Unlike other organs, these local determinants are significantly influenced by regional hydrostatic and alveolar pressures. Those effects on blood flow distribution are further magnified by the large vertical height of the human lung and the relatively low intravascular pressures in the pulmonary circulation. While the distribution of perfusion is largely due to passive determinants such as vascular geometry and hydrostatic pressures, active mechanisms such as vasoconstriction induced by local hypoxia can also redistribute blood flow. This chapter reviews the determinants of regional lung perfusion with a focus on vascular tree geometry, vertical gradients induced by gravity, the interactions between vascular and surrounding alveolar pressures, and hypoxic pulmonary vasoconstriction. While each of these determinants of perfusion distribution can be examined in isolation, the distribution of blood flow is dynamically determined and each component interacts with the others so that a change in one region of the lung influences the distribution of blood flow in other lung regions.

Fractal analysis in radiological and nuclear medicine perfusion imaging: a systematic review

OBJECTIVES

To provide an overview of recent research in fractal analysis of tissue perfusion imaging, using standard radiological and nuclear medicine imaging techniques including computed tomography (CT), magnetic resonance imaging (MRI), ultrasound, positron emission tomography (PET) and single-photon emission computed tomography (SPECT) and to discuss implications for different fields of application.

METHODS

A systematic review of fractal analysis for tissue perfusion imaging was performed by searching the databases MEDLINE (via PubMed), EMBASE (via Ovid) and ISI Web of Science.

RESULTS

Thirty-seven eligible studies were identified. Fractal analysis was performed on perfusion imaging of tumours, lung, myocardium, kidney, skeletal muscle and cerebral diseases. Clinically, different aspects of tumour perfusion and cerebral diseases were successfully evaluated including detection and classification. In physiological settings, it was shown that perfusion under different conditions and in various organs can be properly described using fractal analysis.

CONCLUSIONS

Fractal analysis is a suitable method for quantifying heterogeneity from radiological and nuclear medicine perfusion images under a variety of conditions and in different organs. Further research is required to exploit physiologically proven fractal behaviour in the clinical setting.

KEY POINTS

• Fractal analysis of perfusion images can be successfully performed. • Tumour, pulmonary, myocardial, renal, skeletal muscle and cerebral perfusion have already been examined. • Clinical applications of fractal analysis include tumour and brain perfusion assessment. • Fractal analysis is a suitable method for quantifying perfusion heterogeneity. • Fractal analysis requires further research concerning the development of clinical applications.

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.

Pulmonary blood flow and its distribution in highly trained endurance athletes and healthy control subjects

Pulmonary blood flow (PBF) is an important determinant of endurance sports performance, yet studies investigating adaptations of the pulmonary circulation in athletes are scarce. In the present study, we investigated PBF, its distribution, and heterogeneity at baseline and during intravenous systemic adenosine infusion in 10 highly trained male endurance athletes and 10 untrained but fit healthy controls, using positron emission tomography and [15O]water at rest and during adenosine infusion at supine body posture. Our results indicate that PBF at rest and during adenosine stimulation was similar in both groups (213 ± 55 and 563 ± 138 ml·100 ml−1·min−1 in athletes and 206 ± 83 and 473 ± 212 ml·100 ml−1·min−1 in controls, respectively). Although the PBF response to adenosine was thus unchanged in athletes, overall PBF heterogeneity was reduced from rest to adenosine infusion (from 84 ± 18 to 70 ± 19%, P < 0.05), while remaining unchanged in healthy controls (77 ± 16 to 85 ± 33%, P = 0.4). Additionally, there was a marked gravitational influence on general PBF distribution so that clear dorsal dominance was observed both at rest and during adenosine infusion, but training status did not have an effect on this distribution. Regional blood flow heterogeneity was markedly lower in the high-perfusion dorsal areas, both at rest and during adenosine, in all subjects, but flow heterogeneity in dorsal area tended to further decrease in response to adenosine in athletes. In conclusion, reduced blood flow heterogeneity in response to adenosine in endurance athletes may be a reflection of capillary reserve, which is more extensively recruitable in athletes than in matched healthy control subjects.

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.

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

MRI images of pulmonary blood flow using arterial spin labeling (ASL) measure the delivery of magnetically tagged blood to an image plane during one systolic ejection period. However, the method potentially suffers from two problems, each of which may depend on the imaging plane location: 1) the inversion plane is thicker than the imaging plane, resulting in a gap that blood must cross to be detected in the image; and 2) ASL includes signal contributions from tagged blood in conduit vessels (arterial and venous). By using an in silico model of the pulmonary circulation we found the gap reduced the ASL signal to 64-74% of that in the absence of a gap in the sagittal plane and 53-84% in the coronal. The contribution of the conduit vessels varied markedly as a function of image plane ranging from ∼90% of the overall signal in image planes that encompass the central hilar vessels to <20% in peripheral image planes. A threshold cutoff removing voxels with intensities >35% of maximum reduced the conduit vessel contribution to the total ASL signal to ∼20% on average; however, planes with large contributions from conduit vessels underestimate acinar flow due to a high proportion of in-plane flow, making ASL measurements of perfusion impractical. In other image planes, perfusion dominated the resulting ASL images with good agreement between ASL and acinar flow. Similarly, heterogeneity of the ASL signal as measured by relative dispersion is a reliable measure of heterogeneity of the acinar flow distribution in the same image planes.

MRI of the lung (2/3). Why … when … how?

BACKGROUND

Among the modalities for lung imaging, proton magnetic resonance imaging (MRI) has been the latest to be introduced into clinical practice. Its value to replace X-ray and computed tomography (CT) when radiation exposure or iodinated contrast material is contra-indicated is well acknowledged: i.e. for paediatric patients and pregnant women or for scientific use. One of the reasons why MRI of the lung is still rarely used, except in a few centres, is the lack of consistent protocols customised to clinical needs.

METHODS

This article makes non-vendor-specific protocol suggestions for general use with state-of-the-art MRI scanners, based on the available literature and a consensus discussion within a panel of experts experienced in lung MRI.

RESULTS

Various sequences have been successfully tested within scientific or clinical environments. MRI of the lung with appropriate combinations of these sequences comprises morphological and functional imaging aspects in a single examination. It serves in difficult clinical problems encountered in daily routine, such as assessment of the mediastinum and chest wall, and even might challenge molecular imaging techniques in the near future.

CONCLUSION

This article helps new users to implement appropriate protocols on their own MRI platforms. Main Messages • MRI of the lung can be readily performed on state-of-the-art 1.5-T MRI scanners. • Protocol suggestions based on the available literature facilitate its use for routine • MRI offers solutions for complicated thoracic masses with atelectasis and chest wall invasion. • MRI is an option for paediatrics and science when CT is contra-indicated.

Rapid intravenous infusion of 20 mL/kg saline alters the distribution of perfusion in healthy supine humans

Rapid intravenous saline infusion, a model meant to replicate the initial changes leading to pulmonary interstitial edema, increases pulmonary arterial pressure in humans. We hypothesized that this would alter lung perfusion distribution. Six healthy subjects (29 ± 6 years) underwent magnetic resonance imaging to quantify perfusion using arterial spin labeling. Regional proton density was measured using a fast-gradient echo sequence, allowing blood delivered to the slice to be normalized for density and quantified in mL/min/g. Contributions from flow in large conduit vessels were minimized using a flow cutoff value (blood delivered > 35% maximum in mL/min/cm(3)) in order to obtain an estimate of blood delivered to the capillary bed (perfusion). Images were acquired supine at baseline, after infusion of 20 mL/kg saline, and after a short upright recovery period for a single sagittal slice in the right lung during breath-holds at functional residual capacity. Thoracic fluid content measured by impedance cardiography was elevated post-infusion by up to 13% (p<0.0001). Forced expiratory volume in 1s was reduced by 5.1% post-20 mL/kg (p=0.007). Infusion increased perfusion in nondependent lung by up to 16% (6.4 ± 1.6 mL/min/g baseline, 7.3 ± 1.8 post, 7.4 ± 1.7 recovery, p=0.03). Including conduit vessels, blood delivered in dependent lung was unchanged post-infusion; however, was increased at recovery (9.4 ± 2.7 mL/min/g baseline, 9.7 ± 2.0 post, 11.3 ± 2.2 recovery, p=0.01). After accounting for changes in conduit vessels, there were no significant changes in perfusion in dependent lung following infusion (7.8 ± 1.9 mL/min/g baseline, 7.9 ± 2.0 post, 8.5 ± 2.1 recovery, p=0.36). There were no significant changes in lung density. These data suggest that saline infusion increased perfusion to nondependent lung, consistent with an increase in intravascular pressures. Dependent lung may have been "protected" from increases in perfusion following infusion due to gravitational compression of the pulmonary vasculature.

Pulmonary capillary recruitment in response to hypoxia in healthy humans: a possible role for hypoxic pulmonary venoconstriction?

We examined mechanisms by which hypoxia may elicit pulmonary capillary recruitment in humans. On separate occasions, twenty-five healthy adults underwent exposure to intravenous saline infusion (30 ml/kg ∼ 15 min) or 17-h normobaric hypoxia ( [FIO2 = 12.5%). Cardiac output (Q) and pulmonary capillary blood volume (Vc) were measured before and after saline infusion and hypoxic-exposure by a rebreathing method. Pulmonary artery systolic pressure (sPpa) and left ventricular (LV) diastolic function were assessed before and after hypoxic-exposure via echocardiography. Saline infusion increased Q and Vc (P < 0.05) with no change in Vc/Q (P = 0.97). Hypoxic-exposure increased Vc (P < 0.01) despite no change in Q (P = 0.25), increased sPpa (P < 0.01), and impaired LV relaxation. Multiple regression suggested that ∼ 37% of the hypoxia-mediated increase in Vc was attributable to alterations in Q, sPpa and LV diastolic function. In conclusion, hypoxia-induced pulmonary capillary recruitment in humans is only partly accounted for by changes in Q, sPpa and LV diastolic function. We speculate that hypoxic pulmonary venoconstriction may play a role in such recruitment.

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

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

RATIONALE

Measurement of pulmonary blood flow and lung proton density using MR technique offers high spatial resolution images which can be quantified and the ability to perform repeated measurements under several different physiological conditions. In human studies, PET, SPECT, and CT are commonly used as the alternative techniques. However, these techniques involve exposure to ionizing radiation, and thus are not suitable for repeated measurements in human subjects.

Lung perfusion measured using magnetic resonance imaging: New tools for physiological insights into the pulmonary circulation

Since the lung receives the entire cardiac output, sophisticated imaging techniques are not required in order to measure total organ perfusion. However, for many years studying lung function has required physiologists to consider the lung as a single entity: in imaging terms as a single voxel. Since imaging, and in particular functional imaging, allows the acquisition of spatial information important for studying lung function, these techniques provide considerable promise and are of great interest for pulmonary physiologists. In particular, despite the challenges of low proton density and short T2* in the lung, noncontrast MRI techniques to measure pulmonary perfusion have several advantages including high reliability and the ability to make repeated measurements under a number of physiologic conditions. This brief review focuses on the application of a particular arterial spin labeling (ASL) technique, ASL-FAIRER (flow sensitive inversion recovery with an extra radiofrequency pulse), to answer physiologic questions related to pulmonary function in health and disease. The associated measurement of regional proton density to correct for gravitational-based lung deformation (the "Slinky" effect (Slinky is a registered trademark of Pauf-Slinky incorporated)) and issues related to absolute quantification are also discussed.

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

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

Pulmonary perfusion heterogeneity is increased by sustained, heavy exercise in humans

Exercise presents a considerable stress to the pulmonary system and ventilation-perfusion (Va/Q) heterogeneity increases with exercise, affecting the efficiency of gas exchange. In particular, prolonged heavy exercise and maximal exercise are known to increase Va/Q heterogeneity and these changes persist into recovery. We hypothesized that the spatial heterogeneity of pulmonary perfusion would be similarly elevated after prolonged exercise. To test this, athletic subjects (n = 6, Vo(2max) = 61 ml. kg(-1).min(-1)) with exercising Va/Q heterogeneity previously characterized by the multiple inert gas elimination technique (MIGET), performed 45 min of cycle exercise at approximately 70% Vo(2max). MRI arterial spin labeling measures of pulmonary perfusion were acquired pre- and postexercise (at 20, 40, 60 min post) to quantify the spatial distribution in isogravitational (coronal) and gravitationally dependent (sagittal) planes. Regional proton density measurements allowed perfusion to be normalized for density and quantified in milliliters per minute per gram. Mean lung density did not change significantly in either plane after exercise (P = 0.19). Density-normalized perfusion increased in the sagittal plane postexercise (P =or <0.01) but heterogeneity did not (all P >or= 0.18), likely because of perfusion redistribution and vascular recruitment. Density-normalized perfusion was unchanged in the coronal plane postexercise (P = 0.66), however, perfusion heterogeneity was significantly increased as measured by the relative dispersion [RD, pre 0.62(0.07), post 0.82(0.21), P < 0.0001] and geometric standard deviation [GSD, pre 1.74(0.14), post 2.30(0.56), P < 0.005]. These changes in heterogeneity were related to the exercise-induced changes of the log standard deviation of the ventilation distribution, an MIGET index of Va/Q heterogeneity (RD R(2) = 0.68, P < 0.05, GSD, R(2) = 0.55, P = 0.09). These data are consistent with but not proof of interstitial pulmonary edema as the mechanism underlying exercise-induced increases in both spatial perfusion heterogeneity and Va/Q heterogeneity.