Characterizing pulmonary blood flow distribution measured using arterial spin labeling

Quantitative non-contrast perfusion MRI in the body using arterial spin labeling

Arterial spin labeling (ASL) is a non-invasive magnetic resonance imaging (MRI) method that enables the assessment and the quantification of perfusion without the need for an exogenous contrast agent. ASL was originally developed in the early 1990s to measure cerebral blood flow. The utility of ASL has since then broadened to encompass various organ systems, offering insights into physiological and pathological states. In this review article, we present a synopsis of ASL for quantitative non-contrast perfusion MRI, as a contribution to the special issue titled "Quantitative MRI-how to make it work in the body?" The article begins with an introduction to ASL principles, followed by different labeling strategies, such as pulsed, continuous, pseudo-continuous, and velocity-selective approaches, and their role in perfusion quantification. We proceed to address the technical challenges associated with ASL in the body and outline some of the innovative approaches devised to surmount these issues. Subsequently, we summarize potential clinical applications, challenges, and state-of-the-art ASL methods to quantify perfusion in some of the highly perfused organs in the thorax (lungs), abdomen (kidneys, liver, pancreas), and pelvis (placenta) of the human body. The article concludes by discussing future directions for successful translation of quantitative ASL in body imaging.

Assessing the pulmonary vascular responsiveness to oxygen with proton MRI

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

MR Perfusion Imaging of the Lung

Lung perfusion assessment is critical for diagnosing and monitoring a variety of respiratory conditions. MRI perfusion provides a radiation-free technique, making it an ideal choice for longitudinal imaging in younger populations. This review focuses on the techniques and applications of MRI perfusion, including contrast-enhanced (CE) MRI and non-CE methods such as arterial spin labeling (ASL), fourier decomposition (FD), and hyperpolarized 129-Xenon (129-Xe) MRI. ASL leverages endogenous water protons as tracers for a non-invasive measure of lung perfusion, while FD offers simultaneous measurements of lung perfusion and ventilation, enabling the generation of ventilation/perfusion mapsHyperpolarized 129-Xe MRI emerges as a novel tool for assessing regional gas exchange in the lungs. Despite the promise of MRI perfusion techniques, challenges persist, including competition with other imaging techniques and the need for additional validation and standardization. In conditions such as cystic fibrosis and lung cancer, MRI has displayed encouraging results, whereas in diseases like chronic obstructive pulmonary disease, further validation remains necessary. In conclusion, while MRI perfusion techniques hold immense potential for a comprehensive, non-invasive assessment of lung function and perfusion, their broader clinical adoption hinges on technological advancements, collaborative research, and rigorous validation.

On the application of pseudo-continuous arterial spin labeled MRI for pulmonary perfusion imaging

PURPOSE

To evaluate different approaches for the effective assessment of pulmonary perfusion with a pseudo-continuous arterial spin labeled (pCASL) MRI.

MATERIALS AND METHODS

Four different approaches were evaluated: 1) Cardiac-triggered inferior vena cava (IVC) labeling; 2) IVC labeling with cardiac-triggered acquisition; 3) Right pulmonary artery (RPA) labeling with cardiac-triggered acquisition; and 4) Cardiac-triggered RPA labeling with background suppression (BGS). Each approach was evaluated in 5 healthy volunteers (n = 20) using coefficient of variation (COV) across averages. Approach 4 was also compared against a flow alternating inversion recovery (FAIR).

RESULTS

The IVC labeling (Approach 1) achieved perfusion-weighted images of both lungs, although this approach was more sensitive to variations in heart rate. Cardiac-triggered acquisitions using IVC (Approach 2) and RPA (Approach 3) labeling improved signal consistencies, but were incompatible with BGS. The cardiac-triggered RPA labeling with BGS (Approach 4) achieved a COV of 0.34 ± 0.03 (p < 0.05 compared to IVC labeling approaches) and resulted in perfusion value of 434 ± 64 mL/100 g/min, which was comparable to 451 ± 181 mL/100 g/min measured by FAIR (p = 0.82).

DISCUSSION

Pulmonary perfusion imaging using pCASL-MRI is highly sensitive to cardiac phase, and requires approaches to minimize flow-induced signal variations. Cardiac-triggered RPA labeling with BGS achieves reduced COV and enables robust pulmonary perfusion imaging.

Update on state-of-the-art for arterial spin labeling (ASL) human perfusion imaging outside of the brain

This review article provides an overview of developments for arterial spin labeling (ASL) perfusion imaging in the body (i.e., outside of the brain). It is part of a series of review/recommendation papers from the International Society for Magnetic Resonance in Medicine (ISMRM) Perfusion Study Group. In this review, we focus on specific challenges and developments tailored for ASL in a variety of body locations. After presenting common challenges, organ-specific reviews of challenges and developments are presented, including kidneys, lungs, heart (myocardium), placenta, eye (retina), liver, pancreas, and muscle, which are regions that have seen the most developments outside of the brain. Summaries and recommendations of acquisition parameters (when appropriate) are provided for each organ. We then explore the possibilities for wider adoption of body ASL based on large standardization efforts, as well as the potential opportunities based on recent advances in high/low-field systems and machine-learning. This review seeks to provide an overview of the current state-of-the-art of ASL for applications in the body, highlighting ongoing challenges and solutions that aim to enable more widespread use of the technique in clinical practice.

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.

Velocity-selective arterial spin labeling perfusion MRI: A review of the state of the art and recommendations for clinical implementation

This review article provides an overview of the current status of velocity-selective arterial spin labeling (VSASL) perfusion MRI and is part of a wider effort arising from the International Society for Magnetic Resonance in Medicine (ISMRM) Perfusion Study Group. Since publication of the 2015 consensus paper on arterial spin labeling (ASL) for cerebral perfusion imaging, important advancements have been made in the field. The ASL community has, therefore, decided to provide an extended perspective on various aspects of technical development and application. Because VSASL has the potential to become a principal ASL method because of its unique advantages over traditional approaches, an in-depth discussion was warranted. VSASL labels blood based on its velocity and creates a magnetic bolus immediately proximal to the microvasculature within the imaging volume. VSASL is, therefore, insensitive to transit delay effects, in contrast to spatially selective pulsed and (pseudo-) continuous ASL approaches. Recent technical developments have improved the robustness and the labeling efficiency of VSASL, making it a potentially more favorable ASL approach in a wide range of applications where transit delay effects are of concern. In this review article, we (1) describe the concepts and theoretical basis of VSASL; (2) describe different variants of VSASL and their implementation; (3) provide recommended parameters and practices for clinical adoption; (4) describe challenges in developing and implementing VSASL; and (5) describe its current applications. As VSASL continues to undergo rapid development, the focus of this review is to summarize the fundamental concepts of VSASL, describe existing VSASL techniques and applications, and provide recommendations to help the clinical community adopt VSASL.

A novel nonlinear analysis of blood flow dynamics applied to the human lung

The spatial/temporal dynamics of blood flow in the human lung can be measured noninvasively with magnetic resonance imaging (MRI) using arterial spin labeling (ASL). We report a novel data analysis method using nonlinear prediction to identify dynamic interactions between blood flow units (image voxels), potentially providing a probe of underlying vascular control mechanisms. The approach first estimates the linear relationship (predictability) of one voxel time series with another using correlation analysis, and after removing the linear component estimates the nonlinear relationship with a numerical mutual information approach. Dimensionless global metrics for linear prediction (F_L) and nonlinear prediction (F_NL) represent the average amplitude of fluctuations in one voxel estimated by another voxel, as a percentage of the global average voxel flow. A proof-of-principle test of this approach analyzed experimental data from a study of high-altitude pulmonary edema (HAPE), providing two groups exhibiting known differences in vascular reactivity. Subjects were mountaineers divided into HAPE-susceptible (S, n=4) and HAPE-resistant (R, n=5) groups based on prior history at high altitude. Dynamic ASL measurements in the lung in normoxia (N, FIO₂=0.21) and hypoxia (H, FIO₂=0.13±0.01) were compared. The nonlinear prediction metric F_NL decreased with hypoxia (7.4±1.3(N) vs. 6.3±0.7(H), P=0.03) and was significantly different between groups (7.4±1.2 (R) vs. 6.2±14.1 (S), P=0.03). This proof-of-principle test demonstrates that this nonlinear analysis approach applied to ASL data is sensitive to physiological effects even in small subject cohorts, and potentially can be used in a wide range of studies in health and disease in the lung and other organs.

Pulmonary functional imaging (PFI): A historical review and perspective

PFI Pulmonary Functional Imaging (PFI) refers to visualization and measurement of ventilation, perfusion, gas flow and exchange as well as biomechanics. In this review, we will highlight the historical development of PFI, describing recent advances and listing the various techniques for PFI offered per modality. Challenges PFI is facing and requirements for PFI from a clinical point of view will be pointed out. Hereby the review is meant as an introduction to PFI.

Vaping disrupts ventilation-perfusion matching in asymptomatic users

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

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.

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.

Imaging Pulmonary Blood Flow Using Pseudocontinuous Arterial Spin Labeling (PCASL) With Balanced Steady-State Free-Precession (bSSFP) Readout at 1.5T

BACKGROUND

Quantitative assessment of pulmonary blood flow and visualization of its temporal and spatial distribution without contrast media is of clinical significance.

PURPOSE

To assess the potential of electrocardiogram (ECG)-triggered pseudocontinuous arterial spin labeling (PCASL) imaging with balanced steady-state free-precession (bSSFP) readout to measure lung perfusion under free-breathing (FB) conditions and to study temporal and spatial characteristics of pulmonary blood flow.

STUDY TYPE

Prospective, observational.

SUBJECTS

Fourteen volunteers; three patients with pulmonary embolism.

FIELD STRENGTH/SEQUENCES

1.5T, PCASL-bSSFP.

ASSESSMENT

The pulmonary trunk was labeled during systole. The following examinations were performed: 1) FB and timed breath-hold (TBH) examinations with a postlabeling delay (PLD) of 1000 msec, and 2) TBH examinations with multiple PLDs (100-1500 msec). Scan-rescan measurements were performed in four volunteers and one patient. Images were registered and the perfusion was evaluated in large vessels, small vessels, and parenchyma. Mean structural similarity indices (MSSIM) was computed and time-to-peak (TTP) of parenchymal perfusion in multiple PLDs was evaluated. Image quality reading was performed with three independent blinded readers.

STATISTICAL TESTS

Wilcoxon test to compare MSSIM, perfusion, and Likert scores. Spearman's correlation to correlate TTP and cardiac cycle duration. The repeatability coefficient (RC) and within-subject coefficient of variation (wCV) for scan-rescan measurements. Intraclass correlation coefficient (ICC) for interreader agreement.

RESULTS

Image registration resulted in a significant (P < 0.05) increase of MSSIM. FB perfusion values were 6% higher than TBH (3.28 ± 1.09 vs. 3.10 ± 0.99 mL/min/mL). TTP was highly correlated with individuals' cardiac cycle duration (Spearman = 0.89, P < 0.001). RC and wCV were better for TBH than FB (0.13-0.19 vs. 0.47-1.54 mL/min/mL; 6-7 vs. 19-60%). Image quality was rated very good, with ICCs 0.71-0.89.

DATA CONCLUSION

ECG-triggered PCASL-bSSFP imaging of the lung at 1.5T can provide very good image quality and quantitative perfusion maps even under FB. The course of labeled blood through the lung shows a strong dependence on the individuals' cardiac cycle duration.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 2 J. MAGN. RESON. IMAGING 2020;52:1767-1782.

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.

In silico modeling of oxygen-enhanced MRI of specific ventilation

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

Non-contrast quantitative pulmonary perfusion using flow alternating inversion recovery at 3T: A preliminary study

PURPOSE

To demonstrate the initial feasibility of non-contrast quantitative pulmonary perfusion imaging at 3T using flow alternating inversion recovery (FAIR), and to evaluate the intra-session and inter-session reliability of FAIR measurements at 3T.

MATERIALS AND METHODS

Nine healthy volunteers were imaged using our own implementation of FAIR pulse sequence at 3T. Quantitative FAIR perfusion, both with and without larger pulmonary vessels, was correlated with global phase contrast (PC) measured blood flow in the right pulmonary artery (RPA). The same volunteers were also imaged with SPECT perfusion using technetium-99m-macroaggregated albumin and relative dispersion (RD) was assessed between FAIR and SPECT perfusion. Four additional healthy volunteers were evaluated for FAIR repeatability, using intra-class correlation coefficient (ICC) and Bland-Altman analysis. p<0.05 was considered statistically significant.

RESULTS

FAIR perfusion across all subjects was 858±605mL/100g/min (with vessels) and 629±294mL/100g/min (without vessels) and correlated significantly with the PC measured blood flow in the RPA (r=0.62, p<0.01 with vessels; r=0.73, p<0.001 without vessels). The median RD of FAIR perfusion across all subjects was 0.73 (with vessels) and 0.49 (without vessels), compared against 0.23 with SPECT perfusion. The intra/inter-session ICC of FAIR perfusion with vessels was 0.95/0.59 and improved to 0.96/0.72, when vessels were removed.

CONCLUSIONS

Non-contrast quantitative pulmonary perfusion imaging using FAIR is feasible at 3T. This may serve as a reliable method to assess regional lung perfusion at 3T to characterize and monitor treatment response in chronic lung disease without the concerns of repeated exposure to ionizing radiation or the accumulation of exogenous contrast agent.

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.

Pulmonary perfusion quantification with flow-sensitive inversion recovery (FAIR) UTE MRI in small animal imaging

Blood perfusion in lung parenchyma is an important property for assessing lung function. In small animals, its quantitation is limited even with radioactive isotopes or dynamic contrast-enhanced MRI techniques. In this study, the feasibility flow-sensitive alternating inversion recovery (FAIR) for the quantification of blood flow in lung parenchyma in free breathing rats at 7 T has been investigated. In order to obtain sufficient signal from the short T₂ * lung parenchyma, a 2D ultra-short echo time (UTE) Look-Locker read-out has been implemented. Acquisitions were segmented to maintain acquisition time within an acceptable range. A method to perform retrospective respiratory gating (DC-SG) has been applied to investigate the impact of respiratory movement. Reproducibilities within and between sessions were estimated, and the ability of FAIR-UTE to identify the decrease of lung perfusion under hyperoxic conditions was tested. The implemented technique allowed for the visualization of lung parenchyma with excellent SNR and no respiratory artifact even in ungated acquisitions. Lung parenchyma perfusion was obtained as 32.54 ± 2.26 mL/g/min in the left lung, and 34.09 ± 2.75 mL/g/min in the right lung. Application of retrospective gating significantly but minimally changes the perfusion values, implying that respiratory gating may not be necessary with this center-our acquisition method. A decrease of 10% in lung perfusion was found between normoxic and hyperoxic conditions, proving the feasibility of the FAIR-UTE approach to quantify lung perfusion changes.

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.

Assessment of cerebral blood perfusion reserve with acetazolamide using 3D spiral ASL MRI: Preliminary experience in pediatric patients

PURPOSE

To demonstrate the clinical feasibility of a new non-Cartesian cylindrically-distributed spiral 3D pseudo-continuous arterial spin labeling (pCASL) magnetic resonance imaging (MRI) pulse sequence in pediatric patients in quantifying cerebral blood flow (CBF) response to an acetazolamide (ACZ) vasodilator challenge.

MATERIALS AND METHODS

MRI exams were performed on two 3 Tesla Philips Ingenia systems using 32 channel head coil arrays. After local institutional review board approval, the 3D spiral-based pCASL technique was added to a standard brain MRI exam and evaluated in 13 pediatric patients (average age: 11.7±6.4years, range: 1.4-22.2years). All patients were administered ACZ for clinically indicated reasons. Quantitative whole-brain CBF measurements were computed pre- and post-ACZ to assess cerebrovascular reserve.

RESULTS

3D spiral pCASL data were successfully reconstructed in all 13 cases. In 11 patients, CBF increased 2.8% to 93.2% after administration of ACZ. In the two remaining patients, CBF decreased by 2.4 to 6.0% after ACZ. The group average change in CBF due to ACZ was approximately 25.0% and individual changes were statistically significant (p<0.01) in all patients using a paired t-test analysis. CBF perfusion data were diagnostically useful in supporting conventional MR angiography and clinical findings.

CONCLUSION

3D cylindrically-distributed spiral pCASL MRI provides a robust approach to assess cerebral blood flow and reserve in pediatric patients.

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.

Dynamic contrast enhanced magnetic resonance perfusion imaging in high-risk smokers and smoking-related COPD: correlations with pulmonary function tests and quantitative computed tomography

The study aimed to prospectively evaluate correlations between dynamic contrast-enhanced (DCE) MR perfusion imaging, pulmonary function tests (PFT) and volume quantitative CT in smokers with or without chronic obstructive pulmonary disease (COPD) and to determine the value of DCE-MR perfusion imaging and CT volumetric imaging on the assessment of smokers. According to the ATS/ERS guidelines, 51 male smokers were categorized into five groups: At risk for COPD (n = 8), mild COPD (n = 9), moderate COPD (n = 12), severe COPD (n = 10), and very severe COPD (n = 12). Maximum slope of increase (MSI), positive enhancement integral (PEI), etc. were obtained from MR perfusion data. The signal intensity ratio (RSI) of the PDs and normal lung was calculated (RSI = SIPD/SInormal). Total lung volume (TLV), total emphysema volume (TEV) and emphysema index (EI) were obtained from volumetric CT data. For "at risk for COPD," the positive rate of PDs on MR perfusion images was higher than that of abnormal changes on non-enhanced CT images (p < 0.05). Moderate-to-strong positive correlations were found between all the PFT parameters and SIPD, or RSI (r range 0.445∼0.683, p ≤ 0.001). TEV and EI were negatively correlated better with FEV1/FVC than other PFT parameters (r range -0.48 --0.63, p < 0.001). There were significant differences in RSI and SIPD between "at risk for COPD" and "very severe COPD," and between "mild COPD" and "very severe COPD". Thus, MR perfusion imaging may be a good approach to identify early evidence of COPD and may have potential to assist in classification of COPD.

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.

Lung transplantation in adults and children: putting lung function into perspective

The number of lung transplants performed globally continues to increase year after year. Despite this growing experience, long-term outcomes following lung transplantation continue to fall far short of that described in other solid-organ transplant settings. Chronic lung allograft dysfunction (CLAD) remains common and is the end result of exposure to a multitude of potentially injurious insults that include alloreactivity and infection among others. Central to any description of the clinical performance of the transplanted lung is an assessment of its physiology by pulmonary function testing. Spirometry and the evaluation of forced expiratory volume in 1 s and forced vital capacity, remain core indices that are measured as part of routine clinical follow-up. Spirometry, while reproducible in detecting lung allograft dysfunction, lacks specificity in differentiating the different complications of lung transplantation such as rejection, infection and bronchiolitis obliterans. However, interpretation of spirometry is central to defining the different 'chronic rejection' phenotypes. It is becoming apparent that the maximal lung function achieved following transplantation, as measured by spirometry, is influenced by a number of donor and recipient factors as well as the type of surgery performed (single vs double vs lobar lung transplant). In this review, we discuss the wide range of variables that need to be considered when interpreting lung function testing in lung transplant recipients. Finally, we review a number of novel measurements of pulmonary function that may in the future serve as better biomarkers to detect and diagnose the cause of the failing lung allograft.

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.

Affine transformation registers small scale lung deformation

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

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.

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.

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 perfusion imaging using MRI: clinical application

BACKGROUND

Lung perfusion is one of the key components of oxygenation. It is hampered in pulmonary arterial diseases and secondary due to parenchymal diseases.

METHODS

Assessment is frequently required during the workup of a patient for either of these disease categories.

RESULTS

This review provides insight into imaging techniques, qualitative and quantitative evaluation, and focuses on clinical application of MR perfusion.

CONCLUSION

The two major techniques, non-contrast-enhanced (arterial spin labeling) and contrast-enhanced perfusion techniques, are discussed.

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.

Pulmonary functional imaging: qualitative comparison of Fourier decomposition MR imaging with SPECT/CT in porcine lung

PURPOSE

To compare unenhanced lung ventilation-weighted (VW) and perfusion-weighted (QW) imaging based on Fourier decomposition (FD) magnetic resonance (MR) imaging with the clinical reference standard single photon emission computed tomography (SPECT)/computed tomography (CT) in an animal experiment.

MATERIALS AND METHODS

The study was approved by the local animal care committee. Lung ventilation and perfusion was assessed in seven anesthetized pigs by using a 1.5-T MR imager and SPECT/CT. For time-resolved FD MR imaging, sets of lung images were acquired by using an untriggered two-dimensional balanced steady-state free precession sequence (repetition time, 1.9 msec; echo time, 0.8 msec; acquisition time per image, 118 msec; acquisition rate, 3.33 images per second; flip angle, 75°; section thickness, 12 mm; matrix, 128 × 128). Breathing displacement was corrected with nonrigid image registration. Parenchymal signal intensity was analyzed pixelwise with FD to separate periodic changes of proton density induced by respiration and periodic changes of blood flow. Spectral lines representing respiratory and cardiac frequencies were integrated to calculate VW and QW images. Ventilation and perfusion SPECT was performed after inhalation of dispersed technetium 99m ((99m)Tc) and injection of (99m)Tc-labeled macroaggregated albumin. FD MR imaging and SPECT data were independently analyzed by two physicians in consensus. A regional statistical analysis of homogeneity and pathologic signal changes was performed.

RESULTS

Images acquired in healthy animals by using FD MR imaging and SPECT showed a homogeneous distribution of VW and QW imaging and pulmonary ventilation and perfusion, respectively. The gravitation-dependent signal distribution of ventilation and perfusion in all animals was similarly observed at FD MR imaging and SPECT. Incidental ventilation and perfusion defects were identically visualized by using both modalities.

CONCLUSION

This animal experiment demonstrated qualitative agreement in the assessment of regional lung ventilation and perfusion between contrast media-free and radiation-free FD MR imaging and conventional SPECT/CT.

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.

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

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

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.