Pulmonary blood flow remains fractal down to the level of gas exchange

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.

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.

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.

Integrative Computational Models of Lung Structure-Function Interactions

Anatomically based integrative models of the lung and their interaction with other key components of the respiratory system provide unique capabilities for investigating both normal and abnormal lung function. There is substantial regional variability in both structure and function within the normal lung, yet it remains capable of relatively efficient gas exchange by providing close matching of air delivery (ventilation) and blood delivery (perfusion) to regions of gas exchange tissue from the scale of the whole organ to the smallest continuous gas exchange units. This is despite remarkably different mechanisms of air and blood delivery, different fluid properties, and unique scale-dependent anatomical structures through which the blood and air are transported. This inherent heterogeneity can be exacerbated in the presence of disease or when the body is under stress. Current computational power and data availability allow for the construction of sophisticated data-driven integrative models that can mimic respiratory system structure, function, and response to intervention. Computational models do not have the same technical and ethical issues that can limit experimental studies and biomedical imaging, and if they are solidly grounded in physiology and physics they facilitate investigation of the underlying interaction between mechanisms that determine respiratory function and dysfunction, and to estimate otherwise difficult-to-access measures. © 2021 American Physiological Society. Compr Physiol 11:1501-1530, 2021.

Lung lobe torsion in 15 dogs: Peripheral band sign on ultrasound

The diagnosis of lung lobe torsion in dogs is usually based on radiological, endoscopic, and CT features. Few ultrasonographic descriptions have been published. The purpose of this multicenter, retrospective, and prospective observational study was to investigate the presence of a hypoechoic area forming a pulmonary band or line at the periphery of the twisted lobe on ultrasonography and assess its significance by comparing it to CT and histological findings. Fifteen dogs with lung lobe torsion confirmed surgically or postmortem were included. All had received ultrasonography and CT examinations; 13 had additional histopathological examination performed. In 14 cases, thoracic ultrasonography revealed a peripheral hypoechoic band, overlying areas of scattered, hyperreflecting interfaces in the affected lobe. On CT, central emphysema was surrounded by a peripheral, soft tissue attenuation band, affecting the periphery in 14 cases. No band was observed in one case, in which the lobe was entirely consolidated. Histological examination yielded a comparable peripheral band, consisting of a thickened visceral pleura with or without hemorrhagic necrosis of the underlying pulmonary parenchyma. This peripheral band may be related to the specific fractal organization of airways and vessels, which plays an important role in lung perfusion and ventilation and makes the lung periphery more prone to ischemia. Our findings suggest that the presence of a peripheral hypoechoic band, associated with central emphysema in a noncollapsed lung lobe on ultrasonography, is suggestive of compromised blood supply and air flow, and lung lobe torsion should therefore be suspected.

Using the prone position could help to combat the development of fast hypoxia in some patients with COVID-19

The world is facing an explosive COVID‐19 pandemic. Some cases rapidly develop deteriorating lung function, which causes deep hypoxaemia and requires urgent treatment. Many centres have started treating patients in the prone position, and oxygenation has improved considerably in some cases. Questions have been raised regarding the mechanisms behind this. The mini review provides some insights into the role of supine and prone body positions and summarises the latest understanding of the responsible mechanisms. The scope for discussion is outside the neonatal period and entirely based on experimental and clinical experiences related to adults. The human respiratory system is a complex interplay of many different variables. Therefore, this mini review has prioritised previous and ongoing research to find explanations based on three scientific areas: gravity, lung structure and fractal geometry and vascular regulation. It concludes that gravity is one of the variables responsible for ventilation/perfusion matching but in concert with lung structure and fractal geometry, ventilation and regulation of lung vascular tone. Since ventilation distribution does not change between supine and prone positions, the higher expression of nitric oxide in dorsal lung vessels than in ventral vessels is likely to be the most important mechanism behind enhanced oxygenation in the prone position.

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.

Resilience of three-dimensional sinusoidal networks in liver tissue

Can three-dimensional, microvasculature networks still ensure blood supply if individual links fail? We address this question in the sinusoidal network, a plexus-like microvasculature network, which transports nutrient-rich blood to every hepatocyte in liver tissue, by building on recent advances in high-resolution imaging and digital reconstruction of adult mice liver tissue. We find that the topology of the three-dimensional sinusoidal network reflects its two design requirements of a space-filling network that connects all hepatocytes, while using shortest transport routes: sinusoidal networks are sub-graphs of the Delaunay graph of their set of branching points, and also contain the corresponding minimum spanning tree, both to good approximation. To overcome the spatial limitations of experimental samples and generate arbitrarily-sized networks, we developed a network generation algorithm that reproduces the statistical features of 0.3-mm-sized samples of sinusoidal networks, using multi-objective optimization for node degree and edge length distribution. Nematic order in these simulated networks implies anisotropic transport properties, characterized by an empirical linear relation between a nematic order parameter and the anisotropy of the permeability tensor. Under the assumption that all sinusoid tubes have a constant and equal flow resistance, we predict that the distribution of currents in the network is very inhomogeneous, with a small number of edges carrying a substantial part of the flow-a feature known for hierarchical networks, but unexpected for plexus-like networks. We quantify network resilience in terms of a permeability-at-risk, i.e., permeability as function of the fraction of removed edges. We find that sinusoidal networks are resilient to random removal of edges, but vulnerable to the removal of high-current edges. Our findings suggest the existence of a mechanism counteracting flow inhomogeneity to balance metabolic load on the liver.

Pulmonary Vascular Dynamics

The pulmonary circulation carries deoxygenated blood from the systemic veins through the pulmonary arteries to be oxygenated in the capillaries that line the walls of the pulmonary alveoli. The pulmonary circulation carries the cardiac output with a relatively low driving pressure, and so differs considerably in structure and function from the systemic circulation to maintain a low-resistance vascular system. The pulmonary circulation is often considered to be a quasi-static system in both experimental and computational studies of pulmonary perfusion and its matching to ventilation (air flow) for exchange. However, the system is highly dynamic, with cardiac output and regional perfusion changing with posture, exercise, and over time. Here we review this dynamic system, with a focus on understanding the physiology of pulmonary vascular dynamics across spatial and temporal scales, and the changes to these dynamics that are reflective of disease. © 2019 American Physiological Society. Compr Physiol 9:1081-1100, 2019.

Gravity outweighs the contribution of structure to passive ventilation-perfusion matching in the supine adult human lung

Gravity and matched airway/vascular tree geometries are both hypothesized to be key contributors to ventilation-perfusion (V̇/Q̇) matching in the lung, but their relative contributions are challenging to quantify experimentally. We used a structure-based model to conduct an analysis of the relative contributions of tissue deformation (the "Slinky" effect), other gravitational mechanisms (weight of blood and gravitational gradient in tissue elastic recoil), and matched airway and arterial tree geometry to V̇/Q̇ matching and therefore to total lung oxygen exchange. Our results showed that the heterogeneity in V̇ and Q̇ were lowest and the correlation between V̇ and Q̇ was highest when the only mechanism for V̇/Q̇ matching was either tissue deformation or matched geometry. Heterogeneity in V̇ and Q̇ was highest and their correlation was poorest when all mechanisms were active (that is, at baseline). Eliminating the contribution of matched geometry did not change the correlation between V̇ and Q̇ at baseline. Despite the much larger heterogeneities in V̇ and Q̇ at baseline, the contribution of in-common (to V̇ and Q̇) gravitational mechanisms provided sufficient compensatory V̇/Q̇ matching to minimize the impact on oxygen transfer. In summary, this model predicts that during supine normal breathing under gravitational loading, passive V̇/Q̇ matching is predominantly determined by shared gravitationally induced tissue deformation, compliance distribution, and the effect of the hydrostatic pressure gradient on vessel and capillary size and blood pressures. Contribution from the matching airway and arterial tree geometries in this model is minor under normal gravity in the supine adult human lung. NEW & NOTEWORTHY We use a computational model to systematically analyze contributors to ventilation-perfusion matching in the lung. The model predicts that the multiple effects of gravity are the predominant mechanism in providing passive ventilation-perfusion matching in the supine adult human lung under normal gravitational loads, while geometric matching of airway and arterial trees plays a minor role.

Technical Note: Contrast free angiography of the pulmonary vasculature in live mice using a laboratory x-ray source

Purpose:

In vivo imaging of the pulmonary vasculature in small animals is difficult yet highly desirable in order to allow study of the effects of a host of dynamic biological processes such as hypoxic pulmonary vasoconstriction. Here the authors present an approach for the quantification of changes in the vasculature.

Methods:

A contrast free angiography technique is validated in silico through the use of computer-generated images and in vivo through microcomputed tomography (μCT) of live mice conducted using a laboratory-based x-ray source. Subsequent image processing on μCT data allowed for the quantification of the caliber of pulmonary vasculature without the need for external contrast agents. These measures were validated by comparing with quantitative contrast microangiography in the same mice.

Results:

Quantification of arterial diameters from the method proposed in this study is validated against laboratory-based x-ray contrast microangiography. The authors find that there is a high degree of correlation (R = 0.91) between measures from microangiography and their contrast free method.

Conclusions:

A technique for quantification of murine pulmonary vasculature without the need for contrast is presented. As such, this technique could be applied for longitudinal studies of animals to study changes to vasculature without the risk of premature death in sensitive mouse models of disease. This approach may also be of value in the clinical setting.

Physiology for the pulmonary functional imager

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

Stochastic tracking of infection in a CF lung

Magnetic Resonance Imaging (MRI) and Computed Tomography (CT) scan are the two ubiquitous imaging sources that physicians use to diagnose patients with Cystic Fibrosis (CF) or any other Chronic Obstructive Pulmonary Disease (COPD). Unfortunately the cost constraints limit the frequent usage of these medical imaging procedures. In addition, even though both CT scan and MRI provide mesoscopic details of a lung, in order to obtain microscopic information a very high resolution is required. Neither MRI nor CT scans provide micro level information about the location of infection in a binary tree structure the binary tree structure of the human lung. In this paper we present an algorithm that enhances the current imaging results by providing estimated micro level information concerning the location of the infection. The estimate is based on a calculation of the distribution of possible mucus blockages consistent with available information using an offline Metropolis-Hastings algorithm in combination with a real-time interpolation scheme. When supplemented with growth rates for the pockets of mucus, the algorithm can also be used to estimate how lung functionality as manifested in spirometric tests will change in patients with CF or COPD.

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

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

Acute lung injury and pulmonary vascular permeability: use of transgenic models

Acute lung injury is a general term that describes injurious conditions that can range from mild interstitial edema to massive inflammatory tissue destruction. This review will cover theoretical considerations and quantitative and semi-quantitative methods for assessing edema formation and increased vascular permeability during lung injury. Pulmonary edema can be quantitated directly using gravimetric methods, or indirectly by descriptive microscopy, quantitative morphometric microscopy, altered lung mechanics, high-resolution computed tomography, magnetic resonance imaging, positron emission tomography, or x-ray films. Lung vascular permeability to fluid can be evaluated by measuring the filtration coefficient (Kf) and permeability to solutes evaluated from their blood to lung clearances. Albumin clearances can then be used to calculate specific permeability-surface area products (PS) and reflection coefficients (σ). These methods as applied to a wide variety of transgenic mice subjected to acute lung injury by hyperoxic exposure, sepsis, ischemia-reperfusion, acid aspiration, oleic acid infusion, repeated lung lavage, and bleomycin are reviewed. These commonly used animal models simulate features of the acute respiratory distress syndrome, and the preparation of genetically modified mice and their use for defining specific pathways in these disease models are outlined. Although the initiating events differ widely, many of the subsequent inflammatory processes causing lung injury and increased vascular permeability are surprisingly similar for many etiologies.

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.

Spatial distribution of ventilation and perfusion: mechanisms and regulation

With increasing spatial resolution of regional ventilation and perfusion, it has become more apparent that ventilation and blood flow are quite heterogeneous in the lung. A number of mechanisms contribute to this regional variability, including hydrostatic gradients, pleural pressure gradients, lung compressibility, and the geometry of the airway and vascular trees. Despite this marked heterogeneity in both ventilation and perfusion, efficient gas exchange is possible through the close regional matching of the two. Passive mechanisms, such as the shared effect of gravity and the matched branching of vascular and airway trees, create efficient gas exchange through the strong correlation between ventilation and perfusion. Active mechanisms that match local ventilation and perfusion play little if no role in the normal healthy lung but are important under pathologic conditions.

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.

High correlation of the response of upper and lower lobe small airway epithelium to smoking

The distribution of lung disease induced by inhaled cigarette smoke is complex, depending on many factors. With the knowledge that the small airway epithelium (SAE) is the earliest site of smoking-induced lung disease, and that the SAE gene expression is likely sensitive to inhaled cigarette smoke, we compared upper vs. lower lobe gene expression in the SAE within the same cigarette smokers to determine if the gene expression patterns were similar or different. Active smokers (n = 11) with early evidence of smoking-induced lung disease (normal spirometry but low diffusing capacity) underwent bronchoscopy and brushing of the upper and lower lobe SAE in order to compare upper vs lower lobe genome-wide and smoking-responsive gene expression by microarray. Cluster and principal component analysis demonstrated that, for each individual, the expression of the known SAE smoking-responsive genes were highly correlated in upper and lower lobe pairs, although, as expected, there were differences in the smoking-induced changes in gene expression from individual to individual. These observations support the concept that the heterogeneity observed among smokers in the anatomic distribution of smoking-induced disease are not secondary to the topographic differences in the effects of cigarette smoke on the airway epithelium.

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.

Design of vascular networks: a mathematical model approach

In this paper, methods for the modeling of a realistic vascular tree in 3D space and for hemodynamic calculating throughout the simulated tree structure have been developed. Vascular trees are generated based on the power law relationship. Variations in branching asymmetry and segment length of the vasculature are precisely controlled by the designed gaussian distributions. The resolution limit of current imaging techniques for vessel detectability is simulated by designed pruning technique. On the basis of the generated diameters and lengths, the space locations of the vessel segments are calculated by optimizing the out-of-plane angles of two daughter branches. The generated vascular tree not only follows the power law relationship, but also maximizes the filling volume of the tree structure in 3D space. From the hemodynamic calculation in the simulated vasculature, the processes for which structural changes affect hemodynamic distributions are studied in detail. And also, the fractal nature and resistance of the vascular trees are quantified and compared. The developed method provides some insight into the design of the vascular trees in biology and may be used as a reference for the study of vascular diseases.

Visualizing regional myocardial blood flow in the mouse

RATIONALE

The spatial distribution of blood flow in the hearts of genetically modified mice is a phenotype of interest because derangements in blood flow may precede detectable changes in organ function. However, quantifying the regional distribution of blood flow within organs of mice is challenging because of the small organ volume and the high resolution required to observe spatial differences in flow. Traditional microsphere methods in which the numbers of microspheres per region are indirectly estimated from radioactive counts or extracted fluorescence have been limited to larger organs for 2 reasons; to ensure statistical confidence in the measured flow per region and to be able to physically dissect the organ to acquire spatial information.

OBJECTIVE

To develop methods to quantify and statistically compare the spatial distribution of blood flow within organs of mice.

METHODS AND RESULTS

We developed and validated statistical methods to compare blood flow between regions and with the same regions over time using 15-µm fluorescent microspheres. We then tested this approach by injecting fluorescent microspheres into isolated perfused mice hearts, determining the spatial location of every microsphere in the hearts, and then visualizing regional flow patterns. We demonstrated application of these statistical and visualizing methods in a coronary artery ligation model in mice.

CONCLUSIONS

These new methods provide tools to investigate the spatial and temporal changes in blood flow within organs of mice at a much higher spatial resolution than currently available by other methods.

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.

Cystic fibrosis therapy: a community ecology perspective

Current therapy for cystic fibrosis (CF) focuses on minimizing the microbial community and the host's immune response through the aggressive use of airway clearance techniques, broad-spectrum antibiotics, and treatments that break down the pervasive endobronchial biofilm. Antibiotic selection is typically based on the susceptibility of individual microbial strains to specific antibiotics in vitro. Often this approach cannot accurately predict medical outcomes because of factors both technical and biological. Recent culture-independent assessments of the airway microbial and viral communities demonstrated that the CF airway infection is considerably more complex and dynamic than previously appreciated. Understanding the ecological and evolutionary pressures that shape these communities is critically important for the optimal use of current therapies (in both the choice of therapy and timing of administration) and the development of newer strategies. The climax-attack model (CAM) presented here, grounded in basic ecological principles, postulates the existence of two major functional communities. The attack community consists of transient viral and microbial populations that induce strong innate immune responses. The resultant intense immune response creates microenvironments that facilitate the establishment of a climax community that is slower-growing and inherently resistant to antibiotic therapy. Newer methodologies, including sequence-based metagenomic analysis, can track not only the taxonomic composition but also the metabolic capabilities of these changing viral and microbial communities over time. Collecting this information for CF airways will enable the mathematical modeling of microbial community dynamics during disease progression. The resultant understanding of airway communities and their effects on lung physiology will facilitate the optimization of CF therapies.

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.

A computational model of the topographic distribution of ventilation in healthy human lungs

The topographic distribution of ventilation in the lungs is determined by the interaction of several factors, including lung shape, airway tree geometry, posture, and tissue deformation. Inter-species differences in lung structure-function and technical difficulty in obtaining high resolution imaging of the upright human lung means that it is not straightforward to experimentally determine the contribution of each of these factors to ventilation distribution. We present a mathematical model for predicting the topological distribution of inhaled air in the upright healthy human lung, based on anatomically structured model geometries and biophysical equations for model function. Gravitational deformation of the lung tissue is predicted using a continuum model. Airflow is simulated in anatomically based conducting airways coupled to geometrically simplified terminal acinar units with varying volume-dependent compliances. The predicted ventilation distribution is hence governed by local tissue density and elastic recoil pressure, airway resistance and acinar compliance. Results suggest that there is significant spatial variation in intrinsic tissue properties in the lungs. The model confirms experimental evidence that in the healthy lungs tissue compliance has a far greater effect than airway resistance on the spatial distribution of ventilation, and hence a realistic description of tissue deformation is essential in models of ventilation.

Computational models of the pulmonary circulation: Insights and the move towards clinically directed studies

Biophysically-based computational models provide a tool for integrating and explaining experimental data, observations, and hypotheses. Computational models of the pulmonary circulation have evolved from minimal and efficient constructs that have been used to study individual mechanisms that contribute to lung perfusion, to sophisticated multi-scale and -physics structure-based models that predict integrated structure-function relationships within a heterogeneous organ. This review considers the utility of computational models in providing new insights into the function of the pulmonary circulation, and their application in clinically motivated studies. We review mathematical and computational models of the pulmonary circulation based on their application; we begin with models that seek to answer questions in basic science and physiology and progress to models that aim to have clinical application. In looking forward, we discuss the relative merits and clinical relevance of computational models: what important features are still lacking; and how these models may ultimately be applied to further increasing our understanding of the mechanisms occurring in disease of the pulmonary circulation.

The interdependent contributions of gravitational and structural features to perfusion distribution in a multiscale model of the pulmonary circulation

Recent experimental and imaging studies suggest that the influence of gravity on the measured distribution of blood flow in the lung is largely through deformation of the parenchymal tissue. To study the contribution of hydrostatic effects to regional perfusion in the presence of tissue deformation, we have developed an anatomically structured computational model of the pulmonary circulation (arteries, capillaries, veins), coupled to a continuum model of tissue deformation, and including scale-appropriate fluid dynamics for blood flow in each vessel type. The model demonstrates that both structural and the multiple effects of gravity on the pulmonary circulation make a distinct contribution to the distribution of blood. It shows that postural differences in perfusion gradients can be explained by the combined effect of tissue deformation and extra-acinar blood vessel resistance to flow in the dependent tissue. However, gravitational perfusion gradients persist when the effect of tissue deformation is eliminated, highlighting the importance of the hydrostatic effects of gravity on blood distribution in the pulmonary circulation. Coupling of large- and small-scale models reveals variation in microcirculatory driving pressures within isogravitational planes due to extra-acinar vessel resistance. Variation in driving pressures is due to heterogeneous large-vessel resistance as a consequence of geometric asymmetry in the vascular trees and is amplified by the complex balance of pressures, distension, and flow at the microcirculatory level.

High-resolution spatial measurements of ventilation-perfusion heterogeneity in rats

This study was designed to validate a high-resolution method to measure regional ventilation (VA) in small laboratory animals, and to compare regional Va and perfusion (Q) before and after methacholine-induced bronchoconstriction. A mixture of two different colors of 0.04-microm fluorescent microspheres (FMS) was aerosolized and administered to five anesthetized, mechanically ventilated rats. Those rats also received an intravenous injection of a mixture of two different colors of 15-microm FMS to measure regional blood flow (Q). Five additional rats were labeled with aerosol and intravenous FMS, injected with intravenous methacholine, and then relabeled with a second pair of aerosol and intravenous FMS colors. After death, the lungs were reinflated, frozen, and sequentially sliced in 16-microm intervals on an imaging cryomicrotome set to acquire signal for each of the FMS colors. The reconstructed lung images were sampled using randomly placed 3-mm radius spheres. Va within each sphere was estimated from the aerosol fluorescence signal, and Q was estimated from the number of 15-microm FMS within each sphere. Method error ranged from 6 to 8% for Q and 0.5 to 4.0% for Va. The mean coefficient of variation for Q was 17%, and for Va was 34%. The administration of methacholine altered the distribution of both VA and Q within lung regions, with a change in Va distribution nearly twice as large as that seen for Q. The methacholine-induced changes in Va were not associated with compensatory shifts in Q. Cryomicrotome images of FMS markers provide a high-resolution, anatomically specific means of measuring regional VA/Q responses in the rat.

Effects of posture on regional pulmonary blood flow in rats as measured by PET

Using small animal PET with (68)Ga-radiolabeled human albumin microspheres (Ga-68-microspheres), we investigated the effect of posture on regional pulmonary blood flow (PBF) in normal rats. This in vivo method is noninvasive and quantitative, and it allows for repeated longitudinal measurements. The purpose of the experiment was to quantify spatial differences in PBF in small animals in different postures. Two studies were performed in anesthetized, spontaneously breathing Wistar rats. Study 1 was designed to determine PBF in the prone and supine positions. Ga-68-microspheres were given to five prone and eight supine animals. We found that PBF increased in dorsal regions of supine animals (0.75) more than in prone animals (0.70; P = 0.037), according to a steeper vertical gradient of flow in supine than in prone animals. No differences in spatial heterogeneity were detected. Study 2 was designed to determine the effects of tissue distribution on PBF measurements. Because microspheres remained fixed in the lung, PET was performed on animals in the position in which they received Ga-68-microsphere injections and thereafter in the opposite posture. The distribution of PBF showed a preference for dorsal regions in both positions, but the distribution was dependent on the position during administration of the microspheres. We conclude that PET using Ga-68-microspheres can detect and quantify regional PBF in animals as small as the rat. PBF distributions differed between the prone and supine postures and were influenced by the distribution of lung tissue within the thorax.

Species-specific pulmonary arterial asymmetry determines species differences in regional pulmonary perfusion

The functional significance of differences in pulmonary vascular branching and diameter asymmetry between the human and quadruped lung has not previously been addressed. To evaluate the contribution of branching asymmetry to observable species differences in blood flow gradients, computed distributions of blood flow were compared in structure-based models of the human and ovine pulmonary arteries. The models were derived using a combination of computed tomography and a volume-filling algorithm. Pressure, flow, and deformed vessel diameter were calculated in both species models using equations representing conservation of mass and momentum, and a pressure-diameter relationship. The major difference between the human and ovine results was the presence of a large region of "zone 4" flow and higher mean flows in the central region of the ovine lung compared to that in the human. Heterogeneity in tissue perfusion and the contribution of gravity were similar in both species models; however, the gravitationally directed gradients of perfusion in the human and ovine models were different and each consistent with human and quadruped measurements, respectively. The results suggest that measured species differences in pulmonary perfusion gradients are largely determined by differences in branching asymmetry.

HEMORRHAGE PROGRESSIVELY DISTURBS INTERALVEOLAR PERFUSION IN THE LUNGS OF RATS

Acute hemorrhage is often followed by devastating lung injury. However, why blood loss should lead to lung injury is not known. One possibility is that hemorrhage rapidly disturbs the distribution of microvascular perfusion at the alveolar level, which may be a triggering event for subsequent injury. We showed previously that a 30% blood loss in rats caused significant maldistribution of interalveolar perfusion within 45 min (J Trauma 60:158, 2006). In this report, we describe results of further exploration of this phenomenon. We wanted to know if perfusion distribution was disturbed at 15 min, when vascular pressures were significantly reduced by the blood loss, compared with those at 45 min, when the pressures had returned substantially toward normal. We hemorrhaged rats by removing 30% of their blood volume. We quantified interalveolar perfusion distribution by statistically analyzing the trapping patterns of 4-microm-diameter fluorescent latex particles infused into the pulmonary circulation 15 (red particles) and 45 min (green particles) after blood removal. We used confocal fluorescence microscopy to digitally image the trapping patterns in sections of the air-dried lungs and used pattern analysis to quantify the patterns in tissue image volumes that ranged from 1,300 alveoli to less than 1 alveolus. LogDI, a measure of perfusion maldistribution, increased from 1.00 +/- 0.15 at 15 min after blood loss to 1.62 +/- 0.24 at 45 min (P < 0.001). These values were 0.86 +/- 0.22 (15 min) and 1.12 +/- 0.24 (45 min) in control rats (P = 0.03). Hemorrhage caused the green (45 min)-to-red (15 min) particle distance to decrease from 35.9 +/- 6.5 to 28.0 +/- 5.1 microm (P = 0.024) and the red-to-green particle distance to remain unchanged (30.2 +/- 5.7 microm [red]; 31.5 +/- 10.0 microm [green] [n.s.]). We conclude that hemorrhage caused a progressive increase in interalveolar perfusion maldistribution over 45 min that did not correspond to reduced arterial pressures or altered blood gases. Our particle distance measurements led us to further conclude that this maldistribution occurred in areas that were perfused at 15 min rather than in previously unperfused areas .

Quantification of lung injury using ventilation and perfusion distributions obtained from gamma scintigraphy

This paper explores the potential of isotope V/Q lung scans to quantify lung disease. Areas of restricted perfusion in subjects with a pulmonary embolus (PE) were identified in 3D reconstructions of V/Q images achieved using anatomical data from the Visible Human Project. From these, the extent of lung damage was quantified. Significant differences in the values of both LogSD V and LogSD Q (p > 0.05) obtained from plots of V and Q against Log(V/Q) were found between normal subjects and subjects with a PE, but no correlation was found between either of these parameters and the degree of lung damage in subjects with a PE (p > 0.05). Whilst V/Q values were log normally distributed, the V/Q distributions from the subjects with a PE failed to show the bimodal distribution predicted from theoretical considerations and MIGET measurements previously reported. There was a statistically significant difference in the mean and standard deviation values of the V/Q distributions between normal subject and subjects with a PE (p < 0.05) but not in the median values (p > 0.05). There was no correlation between the mean, median and standard deviation of the distributions from the subjects with a PE and the percentage of damage present (p > 0.05).

Spatial and temporal heterogeneity of regional pulmonary blood flow in piglets

Regional pulmonary blood flow (PBF) in adult animals varies over space and time, following a fractal pattern. We hypothesized that PBF would follow a fractal pattern in young animals. Five, two-week old piglets were sedated and mechanically ventilated. After stabilization, fluorescent microspheres were injected via the right atrium at baseline and then again at 5, 20, 20.5, 40 and 60 min. The lungs were subsequently excised, dried, inflated, and cored into 0.12-cm3 pieces (mean n=561+/-106 per animal) with the spatial coordinates recorded for each piece. Regional PBF was spatially and temporally heterogeneous with a spatial coefficient of variation of 43.3+/-7.9% and a temporal coefficient of variation of 14.3+/-0.4%. PBF followed a fractal pattern with a fractal dimension of 1.20+/-0.06 at 20 min, remaining stable throughout the experiment. PBF decreased with distance from the hilum but did not follow a lobar pattern. Temporal heterogeneity did not significantly increase with time but low flow regions demonstrated the greatest temporal variability throughout the study. Hence, PBF in young piglets was characterized both spatial and temporal heterogeneity.

Biologically variable ventilation improves gas exchange and respiratory mechanics in a model of severe bronchospasm*

OBJECTIVE

Mechanical ventilation can be lifesaving for status asthmaticus, but how best to accomplish mechanical ventilation is unclear. Biologically variable ventilation (mechanical ventilation that emulates healthy variation) and conventional control mode ventilation (monotonously regular) were compared in an animal model of bronchospasm to determine which approach yields better gas exchange and respiratory mechanics.

DESIGN

A randomized prospective trial of biologically variable ventilation vs. control mode ventilation in swine.

SETTING

University research laboratory.

SUBJECTS

Eighteen farm-raised pigs.

INTERVENTIONS

Methacholine was administered as a nebulized aerosol to initiate bronchospasm, defined as doubling of peak inspiratory pressure and respiratory system resistance, and then randomized (n = 9 each group) to either continue control mode ventilation or switch to biologically variable ventilation at the same minute ventilation. Over the next 4 hrs, hemodynamics, blood gases, respiratory mechanics, and carbon dioxide expirograms were recorded hourly. At end-experiment, tracheobronchial lavage was undertaken to determine interleukin-6 and -10 concentrations.

MEASUREMENTS AND MAIN RESULTS

Measurements of physiologic variables and inflammatory cytokines showed that biologically variable ventilation significantly improved gas exchange, with greater arterial oxygen tensions (p = .006; group x time interaction), lower arterial carbon dioxide tensions (p = .0003; group effect), lower peak inspiratory pressures (p = .0001; group x time), greater static compliance (p = .0001; group x time), greater dynamic compliance (p = .0001; group x time), and lower total respiratory system resistance (p = .028; group x time), compared with conventional ventilation. The appearance of inflammatory cytokines in bronchoalveolar lavage fluid (interleukin-6 and -10) was not affected by mode of ventilation.

CONCLUSIONS

In this experimental model, biologically variable ventilation was superior to control mode ventilation in terms of gas exchange and respiratory mechanics during severe bronchospasm.

Measurement of nonlinear pO2 decay in mouse lungs using 3He-MRI

Spatial and temporal variations in oxygen partial pressure (pO(2)) during breath-hold can be exploited to obtain important regional parameters of lung function. In the course of apnea, the oxygen concentration is known to decay exponentially. Therefore, the initial pO(2) (p(0)) can be used to represent local ventilation, and the oxygen depletion time constant can characterize perfusion. The protocol, based on a nonlinear model of pO(2) decay, was validated in six healthy mice. Parametric maps of p(0) and oxygen depletion time constant were obtained for pure (3)He and (3)He/air mixture. The mean measured values of p(0) were 77 +/- 9 mbar for the pure (3)He insufflation and 107 +/- 5 mbar for (3)He/air mixture, in agreement with the predefined p(0) values: 75 +/- 15 mbar and 123 +/- 15 mbar, respectively. The mean measured oxygen depletion time constants were 6.5 +/- 0.2 s for pure (3)He and 7.1 +/- 0.8 s for the (3)He/air mixture, in agreement with physiology.

Quantifying the genetic influence on mammalian vascular tree structure

The ubiquity of fractal vascular trees throughout the plant and animal kingdoms is postulated to be due to evolutionary advantages conferred through efficient distribution of nutrients to multicellular organisms. The implicit, and untested, assertion in this theory is that the geometry of vascular trees is heritable. Because vascular trees are constructed through the iterative use of signaling pathways modified by local factors at each step of the branching process, we sought to investigate how genetic and nongenetic influences are balanced to create vascular trees and the regional distribution of nutrients through them. We studied the spatial distribution of organ blood flow in armadillos because they have genetically identical littermates, allowing us to quantify the genetic influence. We determined that the regional distribution of blood flow is strongly correlated between littermates (r(2) = 0.56) and less correlated between unrelated animals (r(2) = 0.36). Using an ANOVA model, we estimate that 67% of the regional variability in organ blood flow is genetically controlled. We also used fractal analysis to characterize the distribution of organ blood flow and found shared patterns within the lungs and hearts of related animals, suggesting common control over the vascular development of these two organs. We conclude that the geometries of fractal vascular trees are heritable and could be selected through evolutionary pressures. Furthermore, considerable postgenetic modifications may allow vascular trees to adapt to local factors and provide a flexibility that would not be possible in a rigid system.

Alveolar oxygen partial pressure and oxygen depletion rate mapping in rats using 3He ventilation imaging

A hyperpolarized 3He ventilation imaging protocol was implemented to assess alveolar pO2 values and the oxygen depletion rate in rats. The imaging protocol, which is based on spiral k-space sampling, was designed to acquire a high signal-to-noise ratio (SNR) T1-weighted ventilation series of images in a single breath-hold. Simulations were performed to estimate the accuracy and dependence of the pO2 imaging protocol on the image SNR and the RF flip-angle determination. The imaging protocol was validated in vitro in phantoms and in vivo in rats. Imaging sessions were carried out for different inhaled O2 concentrations ranging from 20% to 40%. Parametric maps of alveolar pO2 and oxygen depletion rate were generated from the series of images. For each investigated animal, the differences in measured alveolar pO2 values are in agreement with the changes in inhaled O2 concentration. The oxygen depletion rates, ranging between 0.7 and 8.0 mbar s-1, are in close agreement with the published values for healthy rats.

Distribution of blood flow and ventilation in the lung: gravity is not the only factor

Current textbooks in anaesthesia describe how gravity affects the regional distribution of ventilation and blood flow in the lung, in terms of vertical gradients of pleural pressure and pulmonary vascular pressures. This concept fails to explain some of the clinical features of disturbed lung function. Evidence now suggests that gravity has a less important role in the variation of regional distribution than structural features of the airways and blood vessels. We review more recent studies that used a variety of methods: external radioactive counters, measurements using inhaled and injected particles, and computer tomography scans. These give a higher spatial resolution of regional blood flow and ventilation. The matching between ventilation and blood flow in these small units of lung is considered; the effects of microgravity, increased gravity, and different postures are reviewed, and the application of these findings to conditions such as acute lung injury is discussed. Down to the scale of the acinus, there is considerable heterogeneity in the distribution of both ventilation and blood flow. However, the matching of blood flow with ventilation is well maintained and may result from a common pattern of asymmetric branching of the airways and blood vessels. Disruption of this pattern may explain impaired gas exchange after acute lung injury and explain how the prone position improves gas exchange.

Microsphere maps of regional blood flow and regional ventilation

Systematically mapped samples cut from lungs previously labeled with intravascular and aerosol microspheres can be used to create high-resolution maps of regional perfusion and regional ventilation. With multiple radioactive or fluorescent microsphere labels available, this methodology can compare regional flow responses to different interventions without partial volume effects or registration errors that complicate interpretation of in vivo imaging measurements. Microsphere blood flow maps examined at different levels of spatial resolution have revealed that regional flow heterogeneity increases progressively down to an acinar level of scale. This pattern of scale-dependent heterogeneity is characteristic of a fractal distribution network, and it suggests that the anatomic configuration of the pulmonary vascular tree is the primary determinant of high-resolution regional flow heterogeneity. At approximately 2-cm(3) resolution, the large-scale gravitational gradients of blood flow per unit weight of alveolar tissue account for <5% of the overall flow heterogeneity. Furthermore, regional blood flow per gram of alveolar tissue remains relatively constant with different body positions, gravitational stresses, and exercise. Regional alveolar ventilation is accurately represented by the deposition of inhaled 1.0-microm fluorescent microsphere aerosols, at least down to the approximately 2-cm(3) level of scale. Analysis of these ventilation maps has revealed the same scale-dependent property of regional alveolar ventilation heterogeneity, with a strong correlation between ventilation and blood flow maintained at all levels of scale. The ventilation-perfusion (VA/Q) distributions obtained from microsphere flow maps of normal animals agree with simultaneously acquired multiple inert-gas elimination technique VA/Q distributions, but they underestimate gas-exchange impairment in diffuse lung injury.

Magnetic Resonance Imaging of Uneven Pulmonary Perfusion in Hypoxia in Humans

RATIONALE

Inhomogeneous hypoxic pulmonary vasoconstriction causing regional overperfusion and high capillary pressure is postulated for explaining how high pulmonary artery pressure leads to high-altitude pulmonary edema in susceptible (HAPE-S) individuals.

OBJECTIVE

Because different species of animals also show inhomogeneous hypoxic pulmonary vasoconstriction, we hypothesized that inhomogeneity of lung perfusion in general increases in hypoxia, but is more pronounced in HAPE-S. For best temporal and spatial resolution, regional pulmonary perfusion was assessed by dynamic contrast-enhanced magnetic resonance imaging.

METHODS

Dynamic contrast-enhanced magnetic resonance imaging and echocardiography were performed during normoxia and after 2 h of hypoxia (Fi(O2) = 0.12) in 11 HAPE-S individuals and 10 control subjects. As a measure for perfusion inhomogeneity, the coefficient of variation for two perfusion parameters (peak signal intensity, time-to-peak) was determined for the whole lung and isogravitational slices.

RESULTS

There were no differences in perfusion inhomogeneity between the groups in normoxia. In hypoxia, analysis of coefficients of variation indicated a greater inhomogeneity in all subjects, which was more pronounced in HAPE-S compared with control subjects. Discrimination between HAPE-S and control subjects was best in gravity-dependent lung areas. Pulmonary artery pressure during hypoxia increased from 22 +/- 3 to 53 +/- 9 mm Hg in HAPE-S and 24 +/- 4 to 33 +/- 6 mm Hg in control subjects (mean +/- SD; p < 0.001), respectively.

CONCLUSION

This study shows that hypoxic pulmonary vasoconstriction is inhomogeneous in hypoxia in humans, particularly in HAPE-S individuals where it is accompanied by a greater increase in pulmonary artery pressure compared with control subjects. These findings support the hypothesis of exaggerated and uneven hypoxic pulmonary vasoconstriction in HAPE-S individuals.

Regional pulmonary blood flow in dogs by 4D-X-ray CT

ECG-triggered computed tomography (CT) was used during passage of iodinated contrast to determine regional pulmonary blood flow (PBF) in anesthetized prone/supine dogs. PBF was evaluated as a function of height within the lung (supine and prone) as a function of various normalization methods: raw unit volume data (PBFraw) or PBF normalized to regional fraction air (PBFair), fractional non-air (PBFgm), or relative number of alveoli (PBFalv). The coefficient of variation of PBFraw, PBFair, PBFalv, and PBFgm ranged between 30 and 50% in both lungs and both body postures. The position of maximal flow along the height of the lung (MFP) was calculated for PBFraw, PBFair, PBFalv, and PBFgm. Only PBFgm showed a significantly different MFP height supine vs. prone (whole lung: 2.60 +/- 1.08 cm supine vs. 5.08 +/- 1.61 cm prone, P < 0.01). Mean slopes (ml/min/gm water content/cm) of PBFgm were steeper supine vs. prone in the right (RL) but not left lung (LL) (RL: -0.65 +/- 0.29 supine vs. -0.26 +/- 0.25 prone, P < 0.02; LL: -0.47 +/- 0.21 supine vs. -0.32 +/- 0.26 prone, P > 0.10). Mean slopes of PBFgm vs. vertical lung height were not different prone vs. supine above this vertical height of MFP (VMFP), but PBFgm slopes were steeper in the supine position below the VMFP in the RL. We conclude that PBFgm distribution was posture dependent in RL but not LL. Support of the heart may play a role. We demonstrate that normalization factors can lead to differing attributions of gravitational effects on PBF heterogeneity.

Comparing microsphere deposition and flow modeling in 3D vascular trees

Blood perfusion in organs has been shown to be heterogeneous in a number of cases. At the same time, a number of models of vascular structure and flow have been proposed that also generate heterogeneous perfusion. Although a relationship between local perfusion and vascular structure has to exist, no model has yet been validated as an accurate description of this relationship. A study of perfusion and three-dimensional (3D) arterial structure in individual rat kidneys is presented, which allows comparison between local measurements of perfusion and model-based predictions. High-resolution computed tomography is used to obtain images of both deposited microspheres and of an arterial cast in the same organ. Microsphere deposition is used as an estimate of local perfusion. A 3D cylindrical pipe model of the arterial tree is generated based on an image of the arterial cast. Results of a flow model are compared with local microsphere deposition. High correlation (r(2) > 0.94) was observed between measured and modeled flows through the vascular tree segments. However, the relative dispersion of the microsphere perfusion measurement was two- to threefold higher than perfusion heterogeneity calculated in the flow model. Also, there was no correlation in the residual deviations between the methods. This study illustrates the importance of comparing models of local perfusion with in vivo measurements of perfusion in the same biologically realistic vascular tree.

Gravity-dependent signal gradients on MR images of the lung in supine and prone positions: a comparison with isogravitational signal variability

PURPOSE

To investigate the tendency of proton MR signal intensity (SI) gradients to be steeper in the supine than in the prone body position, and to quantify the relation between gravity-related and isogravitational changes of SI on proton MR images of the lung.

MATERIALS AND METHODS

In eight healthy volunteers, MR images were obtained in the supine and prone positions using a multiple inversion recovery turbo spin-echo (TSE) sequence. The variation in SI along the gravity-dependent direction and within isogravitational planes was measured on a pixel-by-pixel basis. Ratios of slopes were calculated for comparisons among volunteers. Comparisons of ratios were made using Fisher's exact test. Isogravitational variability was compared with the mean SI, the signal-to-noise ratio (SNR), and the image noise.

RESULTS

The average ratios of slopes showed that the overall SI gradient was steeper in the supine than the prone position, with a substantial difference in the supine/prone ratios between inspiration (1.21) and expiration (1.72). In both the supine and prone positions, gravity-dependent gradients were steeper in expiration than in inspiration (P = 0.001). The SI variability along the gravitational direction was larger than the isogravitational variability. The isogravitational variability in turn was larger than the image noise but smaller than the mean SI of the MR images.

CONCLUSION

Gravity-dependent gradients in proton MR SI are steeper in the supine than in the prone position. The magnitudes of these gradients were larger than the isogravitational signal variability, showing that MRI is sensitive to gravitationally induced effects.

Dynamic ventilation imaging from four-dimensional computed tomography

A novel method for dynamic ventilation imaging of the full respiratory cycle from four-dimensional computed tomography (4D CT) acquired without added contrast is presented. Three cases with 4D CT images obtained with respiratory gated acquisition for radiotherapy treatment planning were selected. Each of the 4D CT data sets was acquired during resting tidal breathing. A deformable image registration algorithm mapped each (voxel) corresponding tissue element across the 4D CT data set. From local average CT values, the change in fraction of air per voxel (i.e. local ventilation) was calculated. A 4D ventilation image set was calculated using pairs formed with the maximum expiration image volume, first the exhalation then the inhalation phases representing a complete breath cycle. A preliminary validation using manually determined lung volumes was performed. The calculated total ventilation was compared to the change in contoured lung volumes between the CT pairs (measured volume). A linear regression resulted in a slope of 1.01 and a correlation coefficient of 0.984 for the ventilation images. The spatial distribution of ventilation was found to be case specific and a 30% difference in mass-specific ventilation between the lower and upper lung halves was found. These images may be useful in radiotherapy planning.