The comprehensive imaging-based analysis of the lung: a forum for team science

Evaluating the role of sex-related structure-function differences on airway aerosol transport and deposition

Several experimental studies have found that females have higher particle deposition in the airways than males. This has implications for the delivery of aerosolized therapeutics and for understanding sex differences in respiratory system response to environmental exposures. This study evaluates several factors that potentially contribute to sex differences in particle deposition, using scale-specific structure-function models of one-dimensional (1-D) ventilation distribution, particle transport, and deposition. The impact of gravity, inhalation flow rate, and dead space are evaluated in 12 structure-based models (7 females; 5 males). Females were found to have significantly higher total, bronchial, and alveolar deposition than males across a particle size range from 0.01 to 10 μm. Results suggest that higher deposition fraction in females is due to higher alveolar deposition for smaller particle sizes and higher bronchial deposition for larger particles. Females had higher alveolar deposition in the lower lobes and slightly lower particle concentration in the left upper lobe. Males were found to be more sensitive to changes due to gravity, showing greater reduction in bronchial deposition fraction. Males were also more sensitive to change in inhalation flow rate and to scaling of dead space due to the larger male baseline airway size. Predictions of sex differences in particle deposition-that are consistent with the literature-suggest that sex-based characteristics of lung and airway size interacting with particle size gives rise to differences in regional deposition.NEW & NOTEWORTHY Sex differences in airway tract particle deposition are analyzed using computational models that account for scale-specific structure and function. We show that sex-related differences in lung and airway size can explain experimental observations of increased deposition fraction in females, with females tending toward enhanced fine particle deposition in the alveolar airways and enhanced bronchial deposition for larger particles.

An in silico approach to understanding the interaction between cardiovascular and pulmonary lymphatic dysfunction

The lung is extremely sensitive to interstitial fluid balance, yet the role of pulmonary lymphatics in lung fluid homeostasis and its interaction with cardiovascular pressures is poorly understood. In health, there is a fine balance between fluid extravasated from the pulmonary capillaries into the interstitium and the return of fluid to the circulation via the lymphatic vessels. This balance is maintained by an extremely interdependent system governed by pressures in the fluids (air and blood) and tissue (interstitium), lung motion during breathing, and the permeability of the tissues. Chronic elevation in left atrial pressure (LAP) due to left heart disease increases the capillary blood pressure. The consequent fluid accumulation in the delicate lung tissue increases its weight, decreases its compliance, and impairs gas exchange. This interdependent system is difficult, if not impossible, to study experimentally. Computational modeling provides a unique perspective to analyze fluid movement in the cardiopulmonary vasculature in health and disease. We have developed an initial in silico model of pulmonary lymphatic function using an anatomically structured model to represent ventilation and perfusion and underlying biophysical laws governing fluid transfer at the interstitium. This novel model was tested against increased LAP and noncardiogenic effects (increased permeability). The model returned physiologically reasonable values for all applications, predicting pulmonary edema when LAP reached 25 mmHg and with increased permeability.NEW & NOTEWORTHY This model presents a novel approach to understanding the interaction between cardiac dysfunction and pulmonary lymphatic function, using anatomically structured models and biophysical equations to estimate regional variation in fluid transport from blood to interstitial and lymphatic flux. This fluid transport model brings together advanced models of ventilation, perfusion, and lung mechanics to produce a detailed model of fluid transport in health and various altered pathological conditions.

Simulating Multi-Scale Pulmonary Vascular Function by Coupling Computational Fluid Dynamics With an Anatomic Network Model

The function of the pulmonary circulation is truly multi-scale, with blood transported through vessels from centimeter to micron scale. There are scale-dependent mechanisms that govern the flow in the pulmonary vascular system. However, very few computational models of pulmonary hemodynamics capture the physics of pulmonary perfusion across the spatial scales of functional importance in the lung. Here we present a multi-scale model that incorporates the 3-dimensional (3D) complexities of pulmonary blood flow in the major vessels, coupled to an anatomically-based vascular network model incorporating the multiple contributing factors to capillary perfusion, including gravity. Using the model we demonstrate how we can predict the impact of vascular remodeling and occlusion on both macro-scale functional drivers (flow distribution between lungs, and wall shear stress) and micro-scale contributors to gas exchange. The model predicts interactions between 3D and 1D models that lead to a redistribution of blood between postures, both on a macro- and a micro-scale. This allows us to estimate the effect of posture on left and right pulmonary artery wall shear stress, with predictions varying by 0.75-1.35 dyne/cm² between postures.

Quantitative CT imaging and advanced visualization methods: potential application in novel coronavirus disease 2019 (COVID-19) pneumonia

Increasingly, quantitative lung computed tomography (qCT)-derived metrics are providing novel insights into chronic inflammatory lung diseases, including chronic obstructive pulmonary disease, asthma, interstitial lung disease, and more. Metrics related to parenchymal, airway, and vascular anatomy together with various measures associated with lung function including regional parenchymal mechanics, air trapping associated with functional small airways disease, and dual-energy derived measures of perfused blood volume are offering the ability to characterize disease phenotypes associated with the chronic inflammatory pulmonary diseases. With the emergence of COVID-19, together with its widely varying degrees of severity, its rapid progression in some cases, and the potential for lengthy post-COVID-19 morbidity, there is a new role in applying well-established qCT-based metrics. Based on the utility of qCT tools in other lung diseases, previously validated supervised classical machine learning methods, and emerging unsupervised machine learning and deep-learning approaches, we are now able to provide desperately needed insight into the acute and the chronic phases of this inflammatory lung disease. The potential areas in which qCT imaging can be beneficial include improved accuracy of diagnosis, identification of clinically distinct phenotypes, improvement of disease prognosis, stratification of care, and early objective evaluation of intervention response. There is also a potential role for qCT in evaluating an increasing population of post-COVID-19 lung parenchymal changes such as fibrosis. In this work, we discuss the basis of various lung qCT methods, using case-examples to highlight their potential application as a tool for the exploration and characterization of COVID-19, and offer scanning protocols to serve as templates for imaging the lung such that these established qCT analyses have the best chance at yielding the much needed new insights.

Longitudinal assessment of interstitial lung disease in single lung transplant recipients with scleroderma

OBJECTIVE

To investigate the natural history of fibrotic lung disease in recipients of a single lung transplant for scleroderma-associated interstitial lung disease (ILD).

METHODS

Global ILD (including ground glass, nodular opacities and fibrosis) was categorized into severity quintiles on first and last post-transplant CT scans, and percent fibrosis by manual contouring was also determined, in nine single lung transplant recipients. Quantitative mean lung densities and volumes for the native and allograft lungs were also acquired.

RESULTS

In the native lung, global ILD severity quintile worsened in two cases and percent fibrosis worsened in four cases (range 5-28%). In the lung allograft, one case each developed mild, moderate and severe ILD; of these, new fibrotic ILD (involving <10% of lung) occurred in two cases and acute cellular rejection occurred in one. The average change in native lung density over time was +2.2 Hounsfield Units per year and lung volume +1.4 ml per year, whereas the allograft lung density changed by -5.5 Hounsfield Units per year and total volume +27 ml per year (P = 0.011 and P = 0.039 for native vs allograft density and volume comparisons, respectively).

CONCLUSIONS

While the course of ILD in the native and transplanted lungs varied in this series, these cases illustrate that disease progression is common in the native lung, suggesting that either the immune process continues to target autoantigens or ongoing fibrotic pathways are active in the native lung. Mild lung disease may occur in the allograft after several years due to either allograft rejection or recurrent mild ILD.

Anthropometry-based generation of personalized and population-specific human airway models

Understanding aerosol deposition in the human lung is of great significance in pulmonary toxicology and inhalation pharmacology. Adverse effects of inhaled environmental aerosols and pharmacological efficacy of inhaled therapeutics are dependent on aerosol properties as well as person-specific respiratory tract anatomy and physiology. Anatomical geometry and physiological function of human airways depend on age, gender, weight, fitness, health, and disease status. Tools for the generation of the population- and subject-specific virtual airway anatomical geometry based on anthropometric data and physiological vitals are invaluable in respiratory diagnostics, personalized pulmonary pharmacology, and model-based management of chronic respiratory diseases. Here we present a novel protocol and software framework for the generation of subject-specific airways based on anthropometric measurements of the subject's body, using the anatomical input, and the conventional spirometry, providing the functional (physiological) data. This model can be used for subject-specific simulations of respiration physiology, gas exchange, and aerosol inhalation and deposition.

Wave reflection in an anatomical model of the pulmonary circulation in local and global hypertensive disease

Pulmonary hypertension is a disease of the pulmonary vasculature which can occur for many different reasons, including pathological remodeling of the pulmonary vessels and occlusion of these vessels (amongst others). Pulmonary hypertension can lead to right heart failure and significantly reduces the quality of life of patients living with the condition. It is difficult to distinguish clinically between different classifications of pulmonary hypertension, and doing so accurately is critical for the management of an individual's condition. In addition, different presentations of the disease (e.g. occlusion versus remodeling) can put different strains on the right heart, despite patients having very similar elevations in pulmonary artery pressure. In this study we use an anatomically based model of the pulmonary circulation to predict pressure and flow wave transmission and reflection in two different kinds of pulmonary hypertension - primary pulmonary hypertension, and chronic thrombotic pulmonary embolism (CTEPH), to enable analysis of the impact of disease type on impedance spectra in the main pulmonary artery.

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.

Metrics of lung tissue heterogeneity depend on BMI but not age

Altered parenchymal microstructure and complexity have been observed in older age. How to distinguish between healthy, expected changes and early signs of pathology remains poorly understood. An objective quantitative analysis of computed tomography imaging was conducted to compare mean lung density, tissue density distributions, and tissue heterogeneity in 16 subjects, 8 aged >60 yr who were gender and body mass index matched with 8 subjects aged <30 yr. Subjects had never been smokers, with no prior respiratory disease, and no radiologically identified abnormalities on computed tomography. Volume-controlled breath hold imaging acquired at 80% vital capacity (end inspiration) and 55% vital capacity (end expiration) were used for analysis. Mean lung density was not different between the age groups at end inspiration ( P = 0.806) but was larger in the younger group at end expiration (0.26 ± 0.033 vs. 0.22 ± 0.026, P = 0.008), as is expected due to increased air trapping in the older population. However, gravitational gradients of tissue density did not differ with age; the only difference in distribution of tissue density between the two age groups was a lower density in the apices of the older group at end expiration. The heterogeneity of the lung tissue assessed using two metrics showed significant differences between end inspiration and end expiration, no dependence on age, and a significant relationship with body mass index at both lung volumes when heterogeneity was calculated using quadtree decomposition but only at end expiration when using a fractal dimension. NEW & NOTEWORTHY Changes to lung tissue heterogeneity can be a normal part of aging but can also be an early indicator of disease. We use novel techniques, which have previously not been used on thoracic computed tomography imaging, to quantify lung tissue heterogeneity in young and old healthy subjects. Our results show no dependence on age but a significant correlation with body mass index.

Morphological and functional properties of the conducting human airways investigated by in vivo computed tomography and in vitro MRI

The accurate representation of the human airway anatomy is crucial for understanding and modeling the structure-function relationship in both healthy and diseased lungs. The present knowledge in this area is based on morphometric studies of excised lung casts, partially complemented by in vivo studies in which computed tomography (CT) was used on a small number of subjects. In the present study, we analyzed CT scans of a cohort of healthy subjects and obtained comprehensive morphometric information down to the seventh generation of bronchial branching, including airway diameter, length, branching angle, and rotation angle. Although some of the geometric parameters (such as the child-to-parent branch diameter ratio) are found to be in line with accepted values, for others (such as the branch length-to-diameter ratio) our findings challenge the common assumptions. We also evaluated several metrics of self-similarity, including the fractal dimension of the airway tree. Additionally, we used phase-contrast magnetic resonance imaging (MRI) to obtain the volumetric flow field in the three-dimensional-printed airway model of one of the subjects during steady inhalation. This is used to relate structural and functional parameters and, in particular, to close the power-law relationship between branch flow rate and diameter. The diameter exponent is found to be significantly lower than in the usually assumed Poiseuille regime, which we attribute to the strong secondary (i.e., transverse) velocity component. The strength of the secondary velocity with respect to the axial component exceeds the levels found in idealized airway models and persists within the first seven generations. NEW & NOTEWORTHY We performed a comprehensive computed tomography-based study of the conductive airway morphology in normal human subjects, including branch diameter, length, and mutual angles. We found significant departure from classic homothetic relationships. We also carried out MRI measurements of the three-dimensional inspiratory flow in an anatomy-based model and directly assessed structure-function relationships that have so far been assumed. We found that strong secondary flows (i.e., transverse velocity components) persist through the first seven generations of bronchial branching.

Capturing complexity in pulmonary system modelling

Respiratory disease is a significant problem worldwide, and it is a problem with increasing prevalence. Pathology in the upper airways and lung is very difficult to diagnose and treat, as response to disease is often heterogeneous across patients. Computational models have long been used to help understand respiratory function, and these models have evolved alongside increases in the resolution of medical imaging and increased capability of functional imaging, advances in biological knowledge, mathematical techniques and computational power. The benefits of increasingly complex and realistic geometric and biophysical models of the respiratory system are that they are able to capture heterogeneity in patient response to disease and predict emergent function across spatial scales from the delicate alveolar structures to the whole organ level. However, with increasing complexity, models become harder to solve and in some cases harder to validate, which can reduce their impact clinically. Here, we review the evolution of complexity in computational models of the respiratory system, including successes in translation of models into the clinical arena. We also highlight major challenges in modelling the respiratory system, while making use of the evolving functional data that are available for model parameterisation and testing.

Pulmonary CT and MRI phenotypes that help explain chronic pulmonary obstruction disease pathophysiology and outcomes

Pulmonary x-ray computed tomographic (CT) and magnetic resonance imaging (MRI) research and development has been motivated, in part, by the quest to subphenotype common chronic lung diseases such as chronic obstructive pulmonary disease (COPD). For thoracic CT and MRI, the main COPD research tools, disease biomarkers are being validated that go beyond anatomy and structure to include pulmonary functional measurements such as regional ventilation, perfusion, and inflammation. In addition, there has also been a drive to improve spatial and contrast resolution while at the same time reducing or eliminating radiation exposure. Therefore, this review focuses on our evolving understanding of patient-relevant and clinically important COPD endpoints and how current and emerging MRI and CT tools and measurements may be exploited for their identification, quantification, and utilization. Since reviews of the imaging physics of pulmonary CT and MRI and reviews of other COPD imaging methods were previously published and well-summarized, we focus on the current clinical challenges in COPD and the potential of newly emerging MR and CT imaging measurements to address them. Here we summarize MRI and CT imaging methods and their clinical translation for generating reproducible and sensitive measurements of COPD related to pulmonary ventilation and perfusion as well as parenchyma morphology. The key clinical problems in COPD provide an important framework in which pulmonary imaging needs to rapidly move in order to address the staggering burden, costs, as well as the mortality and morbidity associated with COPD.

Quantifying normal geometric variation in human pulmonary lobar geometry from high resolution computed tomography

Previous studies of the ex vivo lung have suggested significant intersubject variability in lung lobe geometry. A quantitative description of normal lung lobe shape would therefore have value in improving the discrimination between normal population variability in shape and pathology. To quantify normal human lobe shape variability, a principal component analysis (PCA) was performed on high resolution computed tomography (HRCT) imaging of the lung at full inspiration. Volumetric imaging from 22 never-smoking subjects (10 female and 12 male) with normal lung function was included in the analysis. For each subject, an initial finite element mesh geometry was generated from a group of manually selected nodes that were placed at distinct anatomical locations on the lung surface. Each mesh used cubic shape functions to describe the surface curvilinearity, and the mesh was fitted to surface data for each lobe. A PCA was performed on the surface meshes for each lobe. Nine principal components (PCs) were sufficient to capture &gt;90% of the normal variation in each of the five lobes. The analysis shows that lobe size can explain between 20% and 50% of intersubject variability, depending on the lobe considered. Diaphragm shape was the next most significant intersubject difference. When the influence of lung size difference is removed, the angle of the fissures becomes the most significant shape difference, and the variability in relative lobe size becomes important. We also show how a lobe from an independent subject can be projected onto the study population's PCs, demonstrating potential for abnormalities in lobar geometry to be defined in a quantitative manner.

Quantifying tissue heterogeneity using quadtree decomposition

Volumetric computed tomography (CT) imaging provides a three-dimensional map of image intensities from which lung soft tissue density distribution can be estimated. The information gained from analyzing these images can prove valuable in diagnosis of conditions where lung tissue is damaged or has degenerated, and it is also necessary for modeling lung tissue mechanics. This paper presents a new technique for quantifying heterogeneity based on individual CT images, and investigates the heterogeneity of lung tissue in a group of healthy young subjects. It is intended that development of this technique leads to a standard model of classifying heterogeneity in lung tissue, while taking into account variables such as different imaging platforms and resolutions, and the position of the patient during imaging.

MDCT-based quantification of porcine pulmonary arterial morphometry and self-similarity of arterial branching geometry

The pig is frequently used as an experimental model for studies of the pulmonary circulation, yet the branching and dimensional geometry of the porcine pulmonary vasculature remains poorly defined. The purposes of this study are to improve the geometric definition of the porcine pulmonary arteries and to determine whether the arterial tree exhibits self-similarity in its branching geometry. Five animals were imaged using thin slice spiral computed tomography in the prone posture during airway inflation pressure at 25 cmH2O. The luminal diameter and distance from the inlet of the left and right pulmonary arteries were measured along the left and right main arterial pathway in each lung of each animal. A further six minor pathways were measured in a single animal. The similarity in the rate of reduction of diameter with distance of all minor pathways and the two main pathways, along with similarity in the number of branches arising along the pathways, supports self-similarity in the arterial tree. The rate of reduction in diameter with distance from the inlet was not significantly different among the five animals (P > 0.48) when normalized for main pulmonary artery diameter and total main artery pathlength, which supports intersubject similarity. Other metrics to quantify the tree geometry are strikingly similar to those from airways of other quadrupeds, with the exception of a significantly larger length to diameter ratio, which is more appropriate for the vascular tree. A simplifying self-similar model for the porcine pulmonary arteries is proposed to capture the important geometric features of the arterial tree.

Clinical feasibility of pulmonary perfusion analysis using dynamic computed tomography and a gamma residue function

PURPOSE

To create and determine the clinical feasibility of a model based on dynamic computed tomography (CT) and a bolus injection of iodine contrast medium for evaluation of pulmonary perfusion for healthy individuals and for patients with lung diseases.

MATERIALS AND METHODS

We analyzed pulmonary perfusion by means of dynamic 16-row multidetector CT scanning with a gamma residue function with adding a linear component (extended gamma function model) for 20 healthy individuals and in five patients.

RESULTS

Four types of the time-attenuation curve (TAC) were identified for the peripheral lung. Although the TACs of most pixels for the peripheral lung comprised a single peak or a single-peak followed by another increase, no peak was observed for a small proportion of pixels, which either increased linearly or resulted in a delayed peak for healthy subjects. The ratios of these linearly increasing or delayed peak types of lung fields increased for pathological lungs. The analytical results for pathological cases showed that changes in lung perfusion, difficult to identify using only CT imaging, could be detected.

CONCLUSIONS

The extended gamma function model adequately evaluated pulmonary perfusion not only for normal regions, but also for structurally abnormal regions. Regional changes in perfusion could be evaluated by use of our model, and we confirmed its clinical feasibility for pulmonary perfusion analysis.

Functional lung imaging during HFV in preterm rabbits

Although high frequency ventilation (HFV) is an effective mode of ventilation, there is limited information available in regard to lung dynamics during HFV. To improve the knowledge of lung function during HFV we have developed a novel lung imaging and analysis technique. The technique can determine complex lung motion information in vivo with a temporal resolution capable of observing HFV dynamics. Using high-speed synchrotron based phase contrast X-ray imaging and cross-correlation analysis, this method is capable of recording data in more than 60 independent regions across a preterm rabbit lung in excess of 300 frames per second (fps). This technique is utilised to determine regional intra-breath lung mechanics of preterm rabbit pups during HFV. Whilst ventilated at fixed pressures, each animal was ventilated at frequencies of 1, 3, 5 and 10 Hz. A 50% decrease in delivered tidal volume was measured at 10 Hz compared to 1 Hz, yet at the higher frequency a 500% increase in minute activity was measured. Additionally, HFV induced greater homogeneity of lung expansion activity suggesting this ventilation strategy potentially minimizes tissue damage and improves gas mixing. The development of this technique permits greater insight and further research into lung mechanics and may have implications for the improvement of ventilation strategies used to support severe pulmonary trauma and disease.

Blood flow redistribution and ventilation-perfusion mismatch during embolic pulmonary arterial occlusion

Acute pulmonary embolism causes redistribution of blood in the lung, which impairs ventilation/perfusion matching and gas exchange and can elevate pulmonary arterial pressure (PAP) by increasing pulmonary vascular resistance (PVR). An anatomically-based multi-scale model of the human pulmonary circulation was used to simulate pre- and post-occlusion flow, to study blood flow redistribution in the presence of an embolus, and to evaluate whether reduction in perfused vascular bed is sufficient to increase PAP to hypertensive levels, or whether other vasoconstrictive mechanisms are necessary. A model of oxygen transfer from air to blood was included to assess the impact of vascular occlusion on oxygen exchange. Emboli of 5, 7, and 10 mm radius were introduced to occlude increasing proportions of the vasculature. Blood flow redistribution was calculated after arterial occlusion, giving predictions of PAP, PVR, flow redistribution, and micro-circulatory flow dynamics. Because of the large flow reserve capacity (via both capillary recruitment and distension), approximately 55% of the vasculature was occluded before PAP reached clinically significant levels indicative of hypertension. In contrast, model predictions showed that even relatively low levels of occlusion could cause localized oxygen deficit. Flow preferentially redistributed to gravitationally non-dependent regions regardless of occlusion location, due to the greater potential for capillary recruitment in this region. Red blood cell transit times decreased below the minimum time for oxygen saturation (<0.25 s) and capillary pressures became high enough to initiate cell damage (which may result in edema) only after ~80% of the lung was occluded.

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.

An in-vivo computed tomography approach for quantifying porcine pulmonary arterial morphometry

The objective of this study was to develop an in vivo CT imaging-based approach for pulmonary arterial morphometry measurement, and to improve the geometrical basis for studies of the porcine vasculature. The luminal diameter and distance from the inlet of left and right pulmonary arteries, and pulmonary arteries within the lungs of two porcine subjects were measured at inflation pressure of 25 cmH(2)O. The results suggest that the porcine pulmonary arteries have geometric self-similarity, and that this approach will have utility for systematically quantifying pulmonary arterial vessel dimensions in vivo in a larger group of animals.

Pulmonary embolism: predicting disease severity

Pulmonary embolism (PE) is the most common cause of acute pulmonary hypertension, yet it is commonly undiagnosed, with risk of death if not recognized promptly and managed accordingly. Patients typically present with hypoxemia and hypomania, although the presentation varies greatly, being confounded by co-morbidities such as pre-existing cardio-respiratory disease. Previous studies have demonstrated variable patient outcomes in spite of similar extent and distribution of pulmonary vascular occlusion, but the path physiological determinants of outcome remain unclear. Computational models enable exact control over many of the compounding factors leading to functional outcomes and therefore provide a useful tool to understand and assess these mechanisms. We review the current state of pulmonary blood flow models. We present a pilot study within 10 patients presenting with acute PE, where patient-derived vascular occlusions are imposed onto an existing model of the pulmonary circulation enabling predictions of resultant haemodynamic after embolus occlusion. Results show that mechanical obstruction alone is not sufficient to cause pulmonary arterial hypertension, even when up to 65 per cent of lung tissue is occluded. Blood flow is found to preferentially redistribute to the gravitationally non-dependent regions. The presence of an additional downstream occlusion is found to significantly increase pressures.

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.

Regional deposition of particles in an image-based airway model: large-eddy simulation and left-right lung ventilation asymmetry

Regional deposition and ventilation of particles by generation, lobe and lung during steady inhalation in a computed tomography (CT) based human airway model are investigated numerically. The airway model consists of a seven-generation human airway tree, with oral cavity, pharynx and larynx. The turbulent flow in the upper respiratory tract is simulated by large-eddy simulation. The flow boundary conditions at the peripheral airways are derived from CT images at two lung volumes to produce physiologically-realistic regional ventilation. Particles with diameter equal to or greater than 2.5 microns are selected for study because smaller particles tend to penetrate to the more distal parts of the lung. The current generational particle deposition efficiencies agree well with existing measurement data. Generational deposition efficiencies exhibit similar dependence on particle Stokes number regardless of generation, whereas deposition and ventilation efficiencies vary by lobe and lung, depending on airway morphology and airflow ventilation. In particular, regardless of particle size, the left lung receives a greater proportion of the particle bolus as compared to the right lung in spite of greater flow ventilation to the right lung. This observation is supported by the left-right lung asymmetry of particle ventilation observed in medical imaging. It is found that the particle-laden turbulent laryngeal jet flow, coupled with the unique geometrical features of the airway, causes a disproportionate amount of particles to enter the left lung.

A mathematical model of pulmonary gas exchange under inflammatory stress

During a severe local or systemic inflammatory response, immune mediators target lung tissue. This process may lead to acute lung injury and impaired diffusion of gas molecules. Although several mathematical models of gas exchange have been described, none simulate acute lung injury following inflammatory stress. In view of recent laboratory and clinical progress in the understanding of the pathophysiology of acute lung injury, such a mathematical model would be useful. We first derived a partial differential equations model of gas exchange on a small physiological unit of the lung ( approximately 25 alveoli), which we refer to as a respiratory unit (RU). We next developed a simple model of the acute inflammatory response and implemented its effects within a RU, creating a single RU model. Linking multiple RUs with various ventilation/perfusion ratios and taking into account pulmonary venous blood remixing yielded our lung-scale model. Using the lung-scale model, we explored the predicted effects of inflammation on ventilation/perfusion distribution and the resulting pulmonary venous partial pressure oxygen level during systemic inflammatory stresses. This model represents a first step towards the development of anatomically faithful models of gas exchange and ventilation under a broad range of local and systemic inflammatory stimuli resulting in acute lung injury, such as infection and mechanical strain of lung tissue.

On intra- and intersubject variabilities of airflow in the human lungs

The effects of intra- and intersubject variabilities in airway geometry on airflow in the human lungs are investigated by large eddy simulation. The airway models of two human subjects consisting of extra- and intrathoracic airways are reconstructed from CT images. For intrasubject study, airflows at two inspiratory flow rates are simulated on the airway geometries of the same subject with four different levels of truncation. These airway models are the original complete geometry and three geometries obtained by truncating the original one at the subglottis, the supraglottis, and the laryngopharynx, respectively. A comparison of the airflows in the complete geometry model shows that the characteristics of the turbulent laryngeal jet in the trachea are similar regardless of Reynolds number in terms of mean velocities, turbulence statistics, coherent structures, and pressure distribution. The truncated airway models, however, do not produce the similar flow structures observed in the complete geometry. An improved inlet boundary condition is then proposed for the airway model truncated at the laryngopharynx to improve the accuracy of solution. The new boundary condition significantly improves the mean flow. The spectral analysis shows that turbulent characteristics are captured downstream away from the glottis. For intersubject study, although the overall flow characteristics are similar, two morphological factors are found to significantly affect the flows between subjects. These are the constriction ratio of the glottis with respect to the trachea and the curvature and shape of the airways.

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.

The lung physiome: merging imaging-based measures with predictive computational models

Global measurements of the lung provided by standard pulmonary function tests do not give insight into the regional basis of lung function and lung disease. Advances in imaging methodologies, computer technologies, and subject-specific simulations are creating new opportunities to study structure-function relationships in the lung through multidisciplinary research. The digital Human Lung Atlas is an image-based resource compiled from male and female subjects spanning several decades of age. The Atlas comprises both structural and functional measures, and includes computational models derived to match individual subjects for personalized prediction of function. The computational models in the Atlas form part of the Lung Physiome project, which is an international effort to develop integrative models of lung function at all levels of biological organization. The computational models provide mechanistic interpretation of imaging measures; the Atlas provides structural data on which to base model geometry, and functional data against which to test hypotheses. The example of simulating airflow on a subject-specific basis is considered. Methods for deriving multiscale models of the airway geometry for individual subjects in the Atlas are outlined, and methods for modeling turbulent flows in the airway are reviewed.

A SHAPE-NAVIGATED IMAGE DEFORMATION MODEL FOR 4D LUNG RESPIRATORY MOTION ESTIMATION

Intensity modulated radiation therapy (IMRT) for cancers in the lung remains challenging due to the complicated respiratory dynamics. We propose a shape-navigated dense image deformation model to estimate the patient-specific breathing motion using 4D respiratory correlated CT (RCCT) images. The idea is to use the shape change of the lungs, the major motion feature in the thorax image, as a surrogate to predict the corresponding dense image deformation from training.To build the statistical model, dense diffeomorphic deformations between images of all other time points to the image at end expiration are calculated, and the shapes of the lungs are automatically extracted. By correlating the shape variation with the temporally corresponding image deformation variation, a linear mapping function that maps a shape change to its corresponding image deformation is calculated from the training sample. Finally, given an extracted shape from the image at an arbitrary time point, its dense image deformation can be predicted from the pre-computed statistics.The method is carried out on two patients and evaluated in terms of the tumor and lung estimation accuracies. The result shows robustness of the model and suggests its potential for 4D lung radiation treatment planning.

Volumetric xenon-CT imaging of conventional and high-frequency oscillatory ventilation

RATIONALE AND OBJECTIVES

For mechanical ventilation of patients with pulmonary injuries, it has been proposed that high-frequency oscillatory ventilation (HFOV) offers advantages over conventional ventilation (CV); however, these advantages have been difficult to quantify. We used volumetric, dynamic imaging of Xenon (Xe) washout of the canine lung during both HFOV and CV to compare regional ventilation in the two modalities.

MATERIALS AND METHODS

Three anesthetized, mechanically ventilated animals were studied, each at three different ventilator settings. Imaging was performed on an experimental Toshiba 256-slice scanner at 80 kV, 250 mAs, and 0.5-second scans, yielding 12.8 cm of Z-axis coverage. Repeated images were acquired at increasing intervals between 1 and 10 seconds for 90 seconds during HFOV and using retrospective respiratory gating to end-expiration for 60 seconds during CV. Image series were analyzed to quantify regional specific ventilation (sV ) from the regional density washout time constants.

RESULTS

High-quality, high-resolution regional ventilation maps were obtained during both CV and HFOV. Overall ventilation decreased at smaller tidal volume, as expected. Regional sV was more uniform during HFOV compared to CV, but the underlying distribution of lung aeration was similar.

CONCLUSIONS

High-resolution volumetric ventilation maps of the lung may be obtained with the 256-slice multidetector computed tomographic scanner. There is a marked difference in the distribution of regional ventilation between CV and HFOV, with a significant gravitational ventilation gradient in CV that was not present during HFOV. This technique may be useful in exploring the mechanisms by which HFOV improves gas exchange.

A multivariate analysis of risk factors for the air-trapping asthmatic phenotype as measured by quantitative CT analysis

BACKGROUND

Patients with severe asthma have increased physiologically measured air trapping; however, a study using CT measures of air trapping has not been performed. This study was designed to address two hypotheses: (1) air trapping measured by multidetector CT (MDCT) quantitative methodology would be a predictor of a more severe asthma phenotype; and (2) historical, clinical, allergic, or inflammatory risk factors could be identified via multivariate analysis.

METHODS

MDCT scanning of a subset of Severe Asthma Research Program subjects was performed at functional residual capacity. Air trapping was defined as >or= 9.66% of the lung tissue < - 850 Hounsfield units (HU). Subjects classified as having air trapping were then compared to subjects without air trapping on clinical and demographic factors using both univariate and multivariate statistical analyses.

RESULTS

Subjects with air trapping were significantly more likely to have a history of asthma-related hospitalizations, ICU visits, and/or mechanical ventilation. Duration of asthma (odds ratio [OR], 1.42; 95% confidence interval [CI], 1.08 to 1.87), history of pneumonia (OR, 8.55; 95% CI, 2.07 to 35.26), high levels of airway neutrophils (OR, 8.67; 95% CI, 2.05 to 36.57), airflow obstruction (FEV(1)/FVC) [OR, 1.61; 95% CI, 1.21 to 2.14], and atopy (OR, 11.54; 95% CI, 1.97 to 67.70) were identified as independent risk factors associated with the air-trapping phenotype.

CONCLUSIONS

Quantitative CT-determined air trapping in asthmatic subjects identifies a group of individuals at high risk for severe disease. Several independent risk factors for the presence of this phenotype were identified: perhaps most interestingly, history of pneumonia, neutrophilic inflammation, and atopy.

Towards a virtual lung: multi-scale, multi-physics modelling of the pulmonary system

The essential function of the lung, gas exchange, is dependent on adequate matching of ventilation and perfusion, where air and blood are delivered through complex branching systems exposed to regionally varying transpulmonary and transmural pressures. Structure and function in the lung are intimately related, yet computational models in pulmonary physiology usually simplify or neglect structure. The geometries of the airway and vascular systems and their interaction with parenchymal tissue have an important bearing on regional distributions of air and blood, and therefore on whole lung gas exchange, but this has not yet been addressed by modelling studies. Models for gas exchange have typically incorporated considerable detail at the level of chemical reactions, with little thought for the influence of structure. To date, relatively little attention has been paid to modelling at the cellular or subcellular level in the lung, or to linking information from the protein structure/interaction and cellular levels to the operation of the whole lung. We review previous work in developing anatomically based models of the lung, airways, parenchyma and pulmonary vasculature, and some functional studies in which these models have been used. Models for gas exchange at several spatial scales are briefly reviewed, and the challenges and benefits from modelling cellular function in the lung are discussed.

Functional imaging: CT and MRI

Numerous imaging techniques permit evaluation of regional pulmonary function. Contrast-enhanced CT methods now allow assessment of vasculature and lung perfusion. Techniques using spirometric controlled multi-detector row CT allow for quantification of presence and distribution of parenchymal and airway pathology; xenon gas can be employed to assess regional ventilation of the lungs, and rapid bolus injections of iodinated contrast agent can provide a quantitative measure of regional parenchymal perfusion. Advances in MRI of the lung include gadolinium-enhanced perfusion imaging and hyperpolarized gas imaging, which allow functional assessment, including ventilation/perfusion, microscopic air space measurements, and gas flow and transport dynamics.

Two-dimensional extrapolation methods for texture analysis on CT scans

The application of texture analysis to medical images may require the calculation of texture descriptors on regions of interest (ROIs) that are not completely filled by the tissue under analysis. If a texture descriptor is calculated using such "deficient" ROIs, the accuracy and computational speed may be adversely affected. This study applied 198 texture descriptors from five texture classes (first-order statistical, second-order statistical, Fourier, fractal, and Laws' filtered) to lung parenchyma ROIs automatically extracted from the thoracic CT scans of ten patients. Statistically significant differences in the values of 138 of these texture descriptors were demonstrated when calculated on deficient ROIs. Three extrapolation methods (mean fill, tiled fill, and CLEAN deconvolution) then were applied to correct the deficient ROIs. Texture descriptor values were calculated and compared for the original, deficient, and corrected ROIs (based on the three extrapolation methods). Each extrapolation method induced statistically significant improvements in texture descriptor accuracy for some subset of texture descriptors. CLEAN deconvolution improved the greatest number of descriptors, demonstrated the best overall accuracy, and created ROIs that were visually most similar to the original ROIs.

Characteristics of the turbulent laryngeal jet and its effect on airflow in the human intra-thoracic airways

A computational fluid dynamics technique is applied to understand the relative importance of the upper and intra-thoracic airways and their role in determining central airflow patterns with particular attention paid to the importance of turbulence. The geometry of the human upper respiratory tract is derived from volumetric scans of a volunteer imaged via multidetector-row computed tomography. Geometry 1 consists of a mouthpiece, the mouth, the oropharynx, the larynx, and the intra-thoracic airways of up to six generations. Geometry 2 comprises only the intra-thoracic airways. The results show that a curved sheet-like turbulent laryngeal jet is observed only in geometry 1 with turbulence intensity in the trachea varying from 10% to 20%, whereas the turbulence in geometry 2 is negligible. The presence of turbulence is found to increase the maximum localised wall shear stress by three-folds. The proper orthogonal decomposition analysis reveals that the regions of high turbulence intensity are associated with Taylor-Görtler-like vortices. We conclude that turbulence induced by the laryngeal jet could significantly affect airway flow patterns as well as tracheal wall shear stress. Thus, airflow modeling, particularly subject specific evaluations, should consider upper as well as intra-thoracic airway geometry.

Characteristics of airflow in a CT-based ovine lung: a numerical study

The transient airflow in a rigid, asymmetric monopodial sheep (ovine) tracheobronchial tree of up to 13 generations was investigated numerically. The lung geometry was segmented and reconstructed from computed-tomographic (CT) images. The flow characteristics in the image-based sheep airway were compared with the flow patterns produced by a Weibel-based model at prime locations. Boundary conditions were prescribed 1) a velocity profile from experimental data at the inlet and 2) zero pressure at the bronchial outlets. A mesh convergence study was carried out to establish confidence in the model predictions, and gross left-right ventilation was validated against experimental xenon wash-in-washout data. Detailed flow characteristics were investigated at three points in the breathing cycle: 1) peak inhalation, 2) peak exhalation, and 3) transition. Simulation results revealed fundamental differences between airflow in monopodial and bipodial branching airways. Compared with idealized bipodial flow, the flow in the sheep airway was asymmetric and highly vortical, especially during exhalation and transition. The streak lines during the inhalation phase suggest that the left and right upper lobes are ventilated by airflow in the peripheral region of the trachea. This work may contribute to understanding the interplay between structure and function in the lung.

State of the Art. A structural and functional assessment of the lung via multidetector-row computed tomography: phenotyping chronic obstructive pulmonary disease

With advances in multidetector-row computed tomography (MDCT), it is now possible to image the lung in 10 s or less and accurately extract the lungs, lobes, and airway tree to the fifth- through seventh-generation bronchi and to regionally characterize lung density, texture, ventilation, and perfusion. These methods are now being used to phenotype the lung in health and disease and to gain insights into the etiology of pathologic processes. This article outlines the application of these methodologies with specific emphasis on chronic obstructive pulmonary disease. We demonstrate the use of our methods for assessing regional ventilation and perfusion and demonstrate early data that show, in a sheep model, a regionally intact hypoxic pulmonary vasoconstrictor (HPV) response with an apparent inhibition of HPV regionally in the presence of inflammation. We present the hypothesis that, in subjects with pulmonary emphysema, one major contributing factor leading to parenchymal destruction is the lack of a regional blunting of HPV when the regional hypoxia is related to regional inflammatory events (bronchiolitis or alveolar flooding). If maintaining adequate blood flow to inflamed lung regions is critical to the nondestructive resolution of inflammatory events, the pathologic condition whereby HPV is sustained in regions of inflammation would likely have its greatest effect in the lung apices where blood flow is already reduced in the upright body posture.

Computer analysis of computed tomography scans of the lung: a survey

Current computed tomography (CT) technology allows for near isotropic, submillimeter resolution acquisition of the complete chest in a single breath hold. These thin-slice chest scans have become indispensable in thoracic radiology, but have also substantially increased the data load for radiologists. Automating the analysis of such data is, therefore, a necessity and this has created a rapidly developing research area in medical imaging. This paper presents a review of the literature on computer analysis of the lungs in CT scans and addresses segmentation of various pulmonary structures, registration of chest scans, and applications aimed at detection, classification and quantification of chest abnormalities. In addition, research trends and challenges are identified and directions for future research are discussed.

Computational models of structure-function relationships in the pulmonary circulation and their validation

The pulmonary airway, arterial, venous and capillary networks are vast complex branching and converging systems that are mechanically coupled to the surrounding lung tissue. Early studies that examined vascular or airway geometry relied on measurements from casts, but medical imaging now enables measurement of the lung in vivo, at controlled lung volumes. The high-quality data that imaging provides have prompted development of increasingly sophisticated models of the geometry of the airway and pulmonary vascular trees. The accurate spatial relationships between airway, vessel and tissue in these imaging-derived models are necessary for computational analysis that aims to elucidate regional airway-vessel-tissue interactions. Predictions of blood flow through multiscale imaging-derived models of the pulmonary arteries and capillary bed reveal geometry-dependent patterns of perfusion in response to gravity and lung orientation that cannot be predicted with simplified, summary representations of the pulmonary transport trees. Validation of such predictions against measures from functional imaging holds significant potential for explaining and differentiating normal and disease-related heterogeneity in regional blood flow calculated using perfusion imaging.

Computational predictions of pulmonary blood flow gradients: gravity versus structure

A computational model of blood flow through the human pulmonary arterial tree has been developed to investigate the mechanisms contributing to regional pulmonary perfusion in the isolated network when the lung is in different orientations. The arterial geometric model was constructed using a combination of computed tomography and a volume-filling branching algorithm. Equations governing conservation of mass, momentum, and vessel distension, incorporating gravity, were solved to predict pressure, flow, and vessel radius. Analysis of results in the upright posture, with and without gravity, and in the inverted, prone, and supine postures reveals significant flow heterogeneity and a persistent decrease in flow in the cranial and caudal regions for all postures suggesting that vascular geometry makes a major contribution to regional flow with gravity having a lesser role. Results in the isolated arterial tree demonstrate that the vascular path lengths and therefore the positioning of the pulmonary trunk relative to the rest of the network play a significant role in the determination of flow.

Investigation of the relative effects of vascular branching structure and gravity on pulmonary arterial blood flow heterogeneity via an image-based computational model

RATIONALE AND OBJECTIVES

A computational model of blood flow through the human pulmonary arterial tree has been developed to investigate the relative influence of branching structure and gravity on blood flow distribution in the human lung.

MATERIALS AND METHODS

Geometric models of the largest arterial vessels and lobar boundaries were first derived using multidetector row x-ray computed tomography (MDCT) scans. Further accompanying arterial vessels were generated from the MDCT vessel endpoints into the lobar volumes using a volume-filling branching algorithm. Equations governing the conservation of mass and momentum were solved within the geometric model to calculate pressure, velocity, and vessel radius. Blood flow results in the anatomically based model, with and without gravity, and in a symmetric geometric model were compared to investigate their relative contributions to blood flow heterogeneity.

RESULTS

Results showed a persistent blood flow gradient and flow heterogeneity in the absence of gravitational forces in the anatomically based model. Comparison with flow results in the symmetric model revealed that the asymmetric vascular branching structure was largely responsible for producing this heterogeneity. Analysis of average results in varying slice thicknesses illustrated a clear flow gradient because of gravity in "lower resolution" data (thicker slices), but on examination of higher resolution data, a trend was less obvious.

CONCLUSIONS

Results suggest that although gravity does influence flow distribution, the influence of the tree branching structure is also a dominant factor. These results are consistent with high-resolution experimental studies that have demonstrated gravity to be only a minor determinant of blood flow distribution.

High-Resolution CT in Diagnosis of Diffuse Infiltrative Lung Disease

The chest radiograph remains the first imaging modality for the approach to diffuse infiltrative lung disease (DILD), but, 23 years after its introduction, high-resolution CT (HRCT) is still considered the best imaging tool for the evaluation of the pulmonary interstitium and to diagnose and assess DILD. The introduction of multidetector computed tomography (MDCT) has provided the thoracic radiologist with a powerful tool with which to image the lung. Moreover MDCT has enabled radiologists to understand better the functional information contained within CT images of DILD. By focusing on the HRCT signs, patterns, and distributions of abnormalities, and mentioning the clinical aspects and the new recent advances in pulmonary imaging, in this article we provide an overview of a practical approach to the interpretation of the DILD.