CT-Based Local Distribution Metric Improves Characterization of COPD

Local heterogeneity of normal lung parenchyma and small airways disease are associated with COPD severity and progression

BACKGROUND

Small airways disease (SAD) is a major cause of airflow obstruction in COPD patients and has been identified as a precursor to emphysema. Although the amount of SAD in the lungs can be quantified using our Parametric Response Mapping (PRM) approach, the full breadth of this readout as a measure of emphysema and COPD progression has yet to be explored. We evaluated topological features of PRM-derived normal parenchyma and SAD as surrogates of emphysema and predictors of spirometric decline.

METHODS

PRM metrics of normal lung (PRMNorm) and functional SAD (PRMfSAD) were generated from CT scans collected as part of the COPDGene study (n = 8956). Volume density (V) and Euler-Poincaré Characteristic (χ) image maps, measures of the extent and coalescence of pocket formations (i.e., topologies), respectively, were determined for both PRMNorm and PRMfSAD. Association with COPD severity, emphysema, and spirometric measures were assessed via multivariable regression models. Readouts were evaluated as inputs for predicting FEV1 decline using a machine learning model.

RESULTS

Multivariable cross-sectional analysis of COPD subjects showed that V and χ measures for PRMfSAD and PRMNorm were independently associated with the amount of emphysema. Readouts χfSAD (β of 0.106, p < 0.001) and VfSAD (β of 0.065, p = 0.004) were also independently associated with FEV1% predicted. The machine learning model using PRM topologies as inputs predicted FEV1 decline over five years with an AUC of 0.69.

CONCLUSIONS

We demonstrated that V and χ of fSAD and Norm have independent value when associated with lung function and emphysema. In addition, we demonstrated that these readouts are predictive of spirometric decline when used as inputs in a ML model. Our topological PRM approach using PRMfSAD and PRMNorm may show promise as an early indicator of emphysema onset and COPD progression.

Quantitative CT of Normal Lung Parenchyma and Small Airways Disease Topologies are Associated With COPD Severity and Progression

Objectives

Small airways disease (SAD) is a major cause of airflow obstruction in COPD patients, and has been identified as a precursor to emphysema. Although the amount of SAD in the lungs can be quantified using our Parametric Response Mapping (PRM) approach, the full breadth of this readout as a measure of emphysema and COPD progression has yet to be explored. We evaluated topological features of PRM-derived normal parenchyma and SAD as surrogates of emphysema and predictors of spirometric decline.

Materials and Methods

PRM metrics of normal lung (PRMNorm) and functional SAD (PRMfSAD) were generated from CT scans collected as part of the COPDGene study (n=8956). Volume density (V) and Euler-Poincaré Characteristic (χ) image maps, measures of the extent and coalescence of pocket formations (i.e., topologies), respectively, were determined for both PRMNorm and PRMfSAD. Association with COPD severity, emphysema, and spirometric measures were assessed via multivariable regression models. Readouts were evaluated as inputs for predicting FEV1 decline using a machine learning model.

Results

Multivariable cross-sectional analysis of COPD subjects showed that V and χ measures for PRMfSAD and PRMNorm were independently associated with the amount of emphysema. Readouts χfSAD (β of 0.106, p<0.001) and VfSAD (β of 0.065, p=0.004) were also independently associated with FEV1% predicted. The machine learning model using PRM topologies as inputs predicted FEV1 decline over five years with an AUC of 0.69.

Conclusions

We demonstrated that V and χ of fSAD and Norm have independent value when associated with lung function and emphysema. In addition, we demonstrated that these readouts are predictive of spirometric decline when used as inputs in a ML model. Our topological PRM approach using PRMfSAD and PRMNorm may show promise as an early indicator of emphysema onset and COPD progression.

Quantitative CT Characteristics of Cluster Phenotypes in the Severe Asthma Research Program Cohorts

Background Clustering key clinical characteristics of participants in the Severe Asthma Research Program (SARP), a large, multicenter prospective observational study of patients with asthma and healthy controls, has led to the identification of novel asthma phenotypes. Purpose To determine whether quantitative CT (qCT) could help distinguish between clinical asthma phenotypes. Materials and Methods A retrospective cross-sectional analysis was conducted with the use of qCT images (maximal bronchodilation at total lung capacity [TLC], or inspiration, and functional residual capacity [FRC], or expiration) from the cluster phenotypes of SARP participants (cluster 1: minimal disease; cluster 2: mild, reversible; cluster 3: obese asthma; cluster 4: severe, reversible; cluster 5: severe, irreversible) enrolled between September 2001 and December 2015. Airway morphometry was performed along standard paths (RB1, RB4, RB10, LB1, and LB10). Corresponding voxels from TLC and FRC images were mapped with use of deformable image registration to characterize disease probability maps (DPMs) of functional small airway disease (fSAD), voxel-level volume changes (Jacobian), and isotropy (anisotropic deformation index [ADI]). The association between cluster assignment and qCT measures was evaluated using linear mixed models. Results A total of 455 participants were evaluated with cluster assignments and CT (mean age ± SD, 42.1 years ± 14.7; 270 women). Airway morphometry had limited ability to help discern between clusters. DPM fSAD was highest in cluster 5 (cluster 1 in SARP III: 19.0% ± 20.6; cluster 2: 18.9% ± 13.3; cluster 3: 24.9% ± 13.1; cluster 4: 24.1% ± 8.4; cluster 5: 38.8% ± 14.4; P < .001). Lower whole-lung Jacobian and ADI values were associated with greater cluster severity. Compared to cluster 1, cluster 5 lung expansion was 31% smaller (Jacobian in SARP III cohort: 2.31 ± 0.6 vs 1.61 ± 0.3, respectively, P < .001) and 34% more isotropic (ADI in SARP III cohort: 0.40 ± 0.1 vs 0.61 ± 0.2, P < .001). Within-lung Jacobian and ADI SDs decreased as severity worsened (Jacobian SD in SARP III cohort: 0.90 ± 0.4 for cluster 1; 0.79 ± 0.3 for cluster 2; 0.62 ± 0.2 for cluster 3; 0.63 ± 0.2 for cluster 4; and 0.41 ± 0.2 for cluster 5; P < .001). Conclusion Quantitative CT assessments of the degree and intraindividual regional variability of lung expansion distinguished between well-established clinical phenotypes among participants with asthma from the Severe Asthma Research Program study. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Verschakelen in this issue.

Quantitative CT Correlates with Local Inflammation in Lung of Patients with Subtypes of Chronic Lung Allograft Dysfunction

Chronic rejection of lung allografts has two major subtypes, bronchiolitis obliterans syndrome (BOS) and restrictive allograft syndrome (RAS), which present radiologically either as air trapping with small airways disease or with persistent pleuroparenchymal opacities. Parametric response mapping (PRM), a computed tomography (CT) methodology, has been demonstrated as an objective readout of BOS and RAS and bears prognostic importance, but has yet to be correlated to biological measures. Using a topological technique, we evaluate the distribution and arrangement of PRM-derived classifications of pulmonary abnormalities from lung transplant recipients undergoing redo-transplantation for end-stage BOS (N = 6) or RAS (N = 6). Topological metrics were determined from each PRM classification and compared to structural and biological markers determined from microCT and histopathology of lung core samples. Whole-lung measurements of PRM-defined functional small airways disease (fSAD), which serves as a readout of BOS, were significantly elevated in BOS versus RAS patients (p = 0.01). At the core-level, PRM-defined parenchymal disease, a potential readout of RAS, was found to correlate to neutrophil and collagen I levels (p < 0.05). We demonstrate the relationship of structural and biological markers to the CT-based distribution and arrangement of PRM-derived readouts of BOS and RAS.

Bullous Parametric Response Map for Functional Localization of COPD

Advanced bronchoscopic lung volume reduction treatment (BLVR) is now a routine care option for treating patients with severe emphysema. Patterns of low attenuation clusters indicating emphysema and functional small airway disease (fSAD) on paired CT, which may provide additional insights to the target selection of the segmental or subsegmental lobe of the treatments, require further investigation. The low attenuation clusters (LACS) were segmented to identify the scalar and spatial distribution of the lung destructions, in terms of 10 fractions scales of low attenuation density (LAD) located in upper lobes and lower lobes. The LACs of functional small airway disease (fSAD) were delineated by applying the technique of parametric response map (PRM) on the co-registered CT image data. Both emphysematous LACs of inspiratory CT and fSAD LACs on expiratory CT were used to derive the coefficients of the predictive model for estimating the airflow limitation. The voxel-wise severity is then predicted using the regional LACs on the co-registered CT to indicate the functional localization, namely, the bullous parametric response map (BPRM). A total of 100 subjects, 88 patients with mild to very severe COPD and 12 control participants with normal lung functions (FEV₁/FVC % > 70%), were evaluated. Pearson’s correlations between FEV₁/FVC% and LAV%_HU-950 of severe emphysema are − 0.55 comparing to − 0.67 and − 0.62 of LAV%_HU-856 of air-trapping and LAV%_fSAD respectively. Pearson’s correlation between FEV₁/FVC% and FEV₁/FVC% predicted by the proposed model using LAD% of HU-950 and fSAD on BPRM is 0.82 (p < 0.01). The result of the Bullous Parametric Response Map (BPRM) is capable of identifying the less functional area of the lung, where the BLVR treatment is aimed at removing from a hyperinflated area of emphysematous regions.

Feasibility of function-guided lung treatment planning with parametric response mapping

Purpose

Recent advancements in functional lung imaging have been developed to improve clinicians’ knowledge of patient pulmonary condition prior to treatment. Ultimately, it may be possible to employ these functional imaging modalities to tailor radiation treatment plans to optimize patient outcome and mitigate pulmonary complications. Parametric response mapping (PRM) is a computed tomography (CT)–based functional lung imaging method that utilizes a voxel‐wise image analysis technique to classify lung abnormality phenotypes, and has previously been shown to be effective at assessing lung complication risk in diagnostic applications. The purpose of this work was to demonstrate the implementation of PRM guidance in radiotherapy treatment planning.

Methods and materials

A retrospective study was performed with 18 lung cancer patients to test the incorporation of PRM into a radiotherapy planning workflow. Paired inspiration/expiration pretreatment CT scans were acquired and PRM analysis was utilized to classify each voxel as normal, parenchymal disease, small airway disease, and emphysema. Density maps were generated for each PRM classification to contour high density regions of pulmonary abnormalities. Conventional volumetric‐modulated arc therapy and PRM‐guided treatment plans were designed for each patient.

Results

PRM guidance was successfully implemented into the treatment planning process. The inclusion of PRM priorities resulted in statistically significant (p < 0.05) improvements to the V20Gy within the PRM avoidance contours. On average, reductions of 5.4% in the V20Gy(%) were found. The PRM‐guided treatment plans did not significantly increase the dose to the organs at risk or result in insufficient planning target volume coverage, but did increase plan complexity.

Conclusions

PRM guidance was successfully implemented into a treatment planning workflow and shown to be effective for dose redistribution within the lung. This work has provided a framework for the potential clinical implementation of PRM‐guided treatment planning.

Vascular Deformation Mapping for CT Surveillance of Thoracic Aortic Aneurysm Growth

Background Aortic diameter measurements in patients with a thoracic aortic aneurysm (TAA) show wide variation. There is no technique to quantify aortic growth in a three-dimensional (3D) manner. Purpose To validate a CT-based technique for quantification of 3D growth based on deformable registration in patients with TAA. Materials and Methods Patients with ascending and descending TAA with two or more CT angiography studies between 2006 and 2020 were retrospectively identified. The 3D aortic growth was quantified using vascular deformation mapping (VDM), a technique that uses deformable registration to warp a mesh constructed from baseline aortic anatomy. Growth assessments between VDM and clinical CT diameter measurements were compared. Aortic growth was quantified as the ratio of change in surface area at each mesh element (area ratio). Manual segmentations were performed by independent raters to assess interrater reproducibility. Registration error was assessed using manually placed landmarks. Agreement between VDM and clinical diameter measurements was assessed using Pearson correlation and Cohen κ coefficients. Results A total of 38 patients (68 surveillance intervals) were evaluated (mean age, 69 years ± 9 [standard deviation]; 21 women), with TAA involving the ascending aorta (n = 26), descending aorta (n = 10), or both (n = 2). VDM was technically successful in 35 of 38 (92%) patients and 58 of 68 intervals (85%). Median registration error was 0.77 mm (interquartile range, 0.54-1.10 mm). Interrater agreement was high for aortic segmentation (Dice similarity coefficient = 0.97 ± 0.02) and VDM-derived area ratio (bias = 0.0, limits of agreement: -0.03 to 0.03). There was strong agreement (r = 0.85, P < .001) between peak area ratio values and diameter change. VDM detected growth in 14 of 58 (24%) intervals. VDM revealed growth outside the maximally dilated segment in six of 14 (36%) growth intervals, none of which were detected with diameter measurements. Conclusion Vascular deformation mapping provided reliable and comprehensive quantitative assessment of three-dimensional aortic growth and growth patterns in patients with thoracic aortic aneurysms undergoing CT surveillance. Published under a CC BY 4.0 license Online supplemental material is available for this article. See also the editorial by Wieben in this issue.

Computed tomography-identified phenotypes of small airway obstructions in chronic obstructive pulmonary disease

ABSTRACT

Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease characteristic of small airway inflammation, obstruction, and emphysema. It is well known that spirometry alone cannot differentiate each separate component. Computed tomography (CT) is widely used to determine the extent of emphysema and small airway involvement in COPD. Compared with the pulmonary function test, small airway CT phenotypes can accurately reflect disease severity in patients with COPD, which is conducive to improving the prognosis of this disease. CT measurement of central airway morphology has been applied in clinical, epidemiologic, and genetic investigations as an inference of the presence and severity of small airway disease. This review will focus on presenting the current knowledge and methodologies in chest CT that aid in identifying discrete COPD phenotypes.

Why We Should Target Small Airways Disease in Our Management of Chronic Obstructive Pulmonary Disease

For more than 50 years, small airways disease has been considered a key feature of chronic obstructive pulmonary disease (COPD) and a major cause of airway obstruction. Both preventable and treatable, small airways disease has important clinical consequences if left unchecked. Small airways disease is associated with poor spirometry results, increased lung hyperinflation, and poor health status, making the small airways an important treatment target in COPD. The early detection of small airways disease remains the key barrier; if detected early, treatments designed to target small airways may help reduce symptoms and allow patients to maintain their activities. Studies are needed to evaluate the possible role of new drugs and novel drug formulations, inhalers, and inhalation devices for treating small airways disease. These developments will help to improve our management of small airways disease in patients with COPD.

Radiation-induced Hounsfield unit change correlates with dynamic CT perfusion better than 4DCT-based ventilation measures in a novel-swine model

To analyze radiation induced changes in Hounsfield units and determine their correlation with changes in perfusion and ventilation. Additionally, to compare the post-RT changes in human subjects to those measured in a swine model used to quantify perfusion changes and validate their use as a preclinical model. A cohort of 5 Wisconsin Miniature Swine (WMS) were studied. Additionally, 19 human subjects were recruited as part of an IRB approved clinical trial studying functional avoidance radiation therapy for lung cancer and were treated with SBRT. Imaging (a contrast enhanced dynamic perfusion CT in the swine and 4DCT in the humans) was performed prior to and post-RT. Jacobian elasticity maps were calculated on all 4DCT images. Contours were created from the isodose lines to discretize analysis into 10 Gy dose bins. B-spline deformable image registration allowed for voxel-by-voxel comparative analysis in these contours between timepoints. The WMS underwent a research course of 60 Gy in 5 fractions delivered locally to a target in the lung using an MRI-LINAC system. In the WMS subjects, the dose-bin contours were copied onto the contralateral lung, which received < 5 Gy for comparison. Changes in HU and changes in Jacobian were analyzed in these contours. Statistically significant (p < 0.05) changes in the mean HU value post-RT compared to pre-RT were observed in both the human and WMS groups at all timepoints analyzed. The HU increased linearly with dose for both groups. Strong linear correlation was observed between the changes seen in the swine and humans (Pearson coefficient > 0.97, p < 0.05) at all timepoints. Changes seen in the swine closely modeled the changes seen in the humans at 12 months post RT (slope = 0.95). Jacobian analysis showed between 30 and 60% of voxels were damaged post-RT. Perfusion analysis in the swine showed a statistically significant (p < 0.05) reduction in contrast inside the vasculature 3 months post-RT compared to pre-RT. The increases in contrast outside the vasculature was strongly correlated (Pearson Correlation 0.88) with the reduction in HU inside the vasculature but were not correlated with the changes in Jacobians. Radiation induces changes in pulmonary anatomy at 3 months post-RT, with a strong linear correlation with dose. The change in HU seen in the non-vessel lung parenchyma suggests this metric is a potential biomarker for change in perfusion. Finally, this work suggests that the WMS swine model is a promising pre-clinical model for analyzing radiation-induced changes in humans and poses several benefits over conventional swine models.

Unraveling the interplay of image formation, data representation and learning in CT-based COPD phenotyping automation: The need for a meta-strategy

PURPOSE

In the literature on automated phenotyping of chronic obstructive pulmonary disease (COPD), there is a multitude of isolated classical machine learning and deep learning techniques, mostly investigating individual phenotypes, with small study cohorts and heterogeneous meta-parameters, e.g., different scan protocols or segmented regions. The objective is to compare the impact of different experimental setups, i.e., varying meta-parameters related to image formation and data representation, with the impact of the learning technique for subtyping automation for a variety of phenotypes. The identified associations of these parameters with automation performance and their interactions might be a first step towards a determination of optimal meta-parameters, i.e., a meta-strategy.

METHODS

A clinical cohort of 981 patients (53.8 ± 15.1 years, 554 male) was examined. The inspiratory CT images were analyzed to automate the diagnosis of 13 COPD phenotypes given by two radiologists. A benchmark feature set that integrates many quantitative criteria was extracted from the lung and trained a variety of learning algorithms on the first 654 patients (two thirds) and the respective algorithm retrospectively assessed the remaining 327 patients (one third). The automation performance was evaluated by the area under the receiver operating characteristic curve (AUC). 1717 experiments were conducted with varying meta-parameters such as reconstruction kernel, segmented regions and input dimensionality, i.e., number of extracted features. The association of the meta-parameters with the automation performance was analyzed by multivariable general linear model decomposition of the automation performance in the contributions of meta-parameters and the learning technique.

RESULTS

The automation performance varied strongly for varying meta-parameters. For emphysema-predominant phenotypes, an AUC of 93%-95% could be achieved for the best meta-configuration. The airways-predominant phenotypes led to a lower performance of 65%-85%, while smooth kernel configurations on average were unexpectedly superior to those with sharp kernels. The performance impact of meta-parameters, even that of often neglected ones like the missing-data imputation, was in general larger than that of the learning technique. Advanced learning techniques like 3D deep learning or automated machine learning yielded inferior automation performance for non-optimal meta-configurations in comparison to simple techniques with suitable meta-configurations. The best automation performance was achieved by a combination of modern learning techniques and a suitable meta-configuration.

CONCLUSIONS

Our results indicate that for COPD phenotype automation, study design parameters such as reconstruction kernel and the model input dimensionality should be adapted to the learning technique and may be more important than the technique itself. To achieve optimal automation and prediction results, the interaction between input those meta-parameters and the learning technique should be considered. This might be particularly relevant for the development of specific scan protocols for novel learning algorithms, and towards an understanding of good study design for automated phenotyping.

Latent traits of lung tissue patterns in former smokers derived by dual channel deep learning in computed tomography images

Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease and the traditional variables extracted from computed tomography (CT) images may not be sufficient to describe all the topological features of lung tissues in COPD patients. We employed an unsupervised three-dimensional (3D) convolutional autoencoder (CAE)-feature constructor (FC) deep learning network to learn from CT data and derive tissue pattern-clusters jointly. We then applied exploratory factor analysis (EFA) to discover the unobserved latent traits (factors) among pattern-clusters. CT images at total lung capacity (TLC) and residual volume (RV) of 541 former smokers and 59 healthy non-smokers from the cohort of the SubPopulations and Intermediate Outcome Measures in the COPD Study (SPIROMICS) were analyzed. TLC and RV images were registered to calculate the Jacobian (determinant) values for all the voxels in TLC images. 3D Regions of interest (ROIs) with two data channels of CT intensity and Jacobian value were randomly extracted from training images and were fed to the 3D CAE-FC model. 80 pattern-clusters and 7 factors were identified. Factor scores computed for individual subjects were able to predict spirometry-measured pulmonary functions. Two factors which correlated with various emphysema subtypes, parametric response mapping (PRM) metrics, airway variants, and airway tree to lung volume ratio were discriminants of patients across all severity stages. Our findings suggest the potential of developing factor-based surrogate markers for new COPD phenotypes.

The Emerging Role of Quantification of Imaging for Assessing the Severity and Disease Activity of Emphysema, Airway Disease, and Interstitial Lung Disease

There has been an explosion of use for quantitative image analysis in the setting of lung disease due to advances in acquisition protocols and postprocessing technology, including machine and deep learning. Despite the plethora of published papers, it is important to understand which approach has clinical validation and can be used in clinical practice. This paper provides an introduction to quantitative image analysis techniques being used in the investigation of lung disease and focusses on the techniques that have a reasonable clinical validation for being used in clinical trials and patient care.

Chronic lung diseases: prospects for regeneration and repair

COPD and idiopathic pulmonary fibrosis (IPF) together represent a considerable unmet medical need, and advances in their treatment lag well behind those of other chronic conditions. Both diseases involve maladaptive repair mechanisms leading to progressive and irreversible damage. However, our understanding of the complex underlying disease mechanisms is incomplete; with current diagnostic approaches, COPD and IPF are often discovered at an advanced stage and existing definitions of COPD and IPF can be misleading. To halt or reverse disease progression and achieve lung regeneration, there is a need for earlier identification and treatment of these diseases. A precision medicine approach to treatment is also important, involving the recognition of disease subtypes, or endotypes, according to underlying disease mechanisms, rather than the current “one-size-fits-all” approach. This review is based on discussions at a meeting involving 38 leading global experts in chronic lung disease mechanisms, and describes advances in the understanding of the pathology and molecular mechanisms of COPD and IPF to identify potential targets for reversing disease degeneration and promoting tissue repair and lung regeneration. We also discuss limitations of existing disease measures, technical advances in understanding disease pathology, and novel methods for targeted drug delivery.

Treatment outcomes with COPD and IPF are suboptimal. Better understanding of the diseases, such as targetable repair mechanisms, may generate novel therapies, and earlier diagnosis and treatment is needed to stop or even reverse disease progression.https://bit.ly/2Ga8J1g

Gene expression profiling of bronchial brushes is associated with the level of emphysema measured by computed tomography-based parametric response mapping

Parametric response mapping (PRM) is a computed tomography (CT)-based method to phenotype patients with chronic obstructive pulmonary disease (COPD). It is capable of differentiating emphysema-related air trapping with nonemphysematous air trapping (small airway disease), which helps to identify the extent and localization of the disease. Most studies evaluating the gene expression in smokers and COPD patients related this to spirometric measurements, but none have investigated the relationship with CT-based measurements of lung structure. The current study aimed to examine gene expression profiles of brushed bronchial epithelial cells in association with the PRM-defined CT-based measurements of emphysema (PRM^Emph) and small airway disease (PRM^fSAD). Using the Top Institute Pharma (TIP) study cohort (COPD = 12 and asymptomatic smokers = 32), we identified a gene expression signature of bronchial brushings, which was associated with PRM^Emph in the lungs. One hundred thirty-three genes were identified to be associated with PRM^Emph. Among the most significantly associated genes, CXCL11 is a potent chemokine involved with CD8⁺ T cell activation during inflammation in COPD, indicating that it may play an essential role in the development of emphysema. The PRM^Emph signature was then replicated in two independent data sets. Pathway analysis showed that the PRM^Emph signature is associated with proinflammatory and notch signaling pathways. Together these findings indicate that airway epithelium may play a role in the development of emphysema and/or may act as a biomarker for the presence of emphysema. In contrast, its role in relation to functional small airways disease is less clear.

Predictive Modelling of Lung Function using Emphysematous Density Distribution

Target lung tissue selection remains a challenging task to perform for treating severe emphysema with lung volume reduction (LVR). In order to target the treatment candidate, the percentage of low attenuation volume (LAV%) representing the proportion of emphysema volume to whole lung volume is measured using computed tomography (CT) images. Although LAV% have shown to have a correlation with lung function in patients with chronic obstructive pulmonary disease (COPD), similar measurements of LAV% in whole lung or lobes may have large variations in lung function due to emphysema heterogeneity. The functional information of regional emphysema destruction is required for supporting the choice of optimal target. The purpose of this study is to develop an emphysema heterogeneity descriptor for the three-dimensional emphysematous bullae according to the size variations of emphysematous density (ED) and their spatial distribution. The second purpose is to derive a predictive model of airflow limitation based on the regional emphysema heterogeneity. Deriving the bullous representation and grouping them into four scales in the upper and lower lobes, a predictive model is computed using the linear model fitting to estimate the severity of lung function. A total of 99 subjects, 87 patients with mild to very severe COPD (Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage I~IV) and 12 control participants with normal lung functions (forced expiratory volume in one second (FEV₁)/forced vital capacity (FVC) > 0.7) were evaluated. The final model was trained with stratified cross-validation on randomly selected 75% of the dataset (n = 76) and tested on the remaining dataset (n = 23). The dispersed cases of LAV% inconsistent with their lung function outcome were evaluated, and the correlation study suggests that comparing to LAV of larger bullae, the widely spread smaller bullae with equivalent LAV has a larger impact on lung function. The testing dataset has the correlation of r = -0.76 (p < 0.01) between the whole lung LAV% and FEV₁/FVC, whereas using two ED % of scales and location-dependent variables to predict the emphysema-associated FEV₁/FVC, the results shows their correlation of 0.82 (p < 0.001) with clinical FEV₁/FVC.

Influence of Inspiratory/Expiratory CT Registration on Quantitative Air Trapping

RATIONALE AND OBJECTIVES

The aim of this study was to assess variability in quantitative air trapping (QAT) measurements derived from spatially aligned expiration CT scans.

MATERIALS AND METHODS

Sixty-four paired CT examinations, from 16 school-age cystic fibrosis subjects examined at four separate time intervals, were used in this study. For each pair, visually inspected lobe segmentation maps were generated and expiration CT data were registered to the inspiration CT frame. Measurements of QAT, the percentage of voxels on the expiration CT scan below a set threshold were calculated for each lobe and whole-lung from the registered expiration CT and compared to the true values from the unregistered data.

RESULTS

A mathematical model, which simulates the effect of variable regions of lung deformation on QAT values calculated from aligned to those from unaligned data, showed the potential for large bias. Assessment of experimental QAT measurements using Bland-Altman plots corroborated the model simulations, demonstrating biases greater than 5% when QAT was approximately 40% of lung volume. These biases were removed when calculating QAT from aligned expiration CT data using the determinant of the Jacobian matrix. We found, by Dice coefficient analysis, good agreement between aligned expiration and inspiration segmentation maps for the whole-lung and all but one lobe (Dice coefficient > 0.9), with only the lingula generating a value below 0.9 (mean and standard deviation of 0.85 ± 0.06).

CONCLUSION

The subtle and predictable variability in corrected QAT observed in this study suggests that image registration is reliable in preserving the accuracy of the quantitative metrics.

Structural and Functional Features on Quantitative Chest Computed Tomography in the Korean Asian versus the White American Healthy Non-Smokers

OBJECTIVE

Considering the different prevalence rates of diseases such as asthma and chronic obstructive pulmonary disease in Asians relative to other races, Koreans may have unique airway structure and lung function. This study aimed to investigate unique features of airway structure and lung function based on quantitative computed tomography (QCT)-imaging metrics in the Korean Asian population (Koreans) as compared with the White American population (Whites).

MATERIALS AND METHODS

QCT data of healthy non-smokers (223 Koreans vs. 70 Whites) were collected, including QCT structural variables of wall thickness (WT) and hydraulic diameter (Dh) and functional variables of air volume, total air volume change in the lung (ΔVair), percent emphysema-like lung (Emph%), and percent functional small airway disease-like lung (fSAD%). Mann-Whitney U tests were performed to compare the two groups.

RESULTS

As compared with Whites, Koreans had smaller volume at inspiration, ΔVair between inspiration and expiration (p < 0.001), and Emph% at inspiration (p < 0.001). Especially, Korean females had a decrease of ΔVair in the lower lobes (p < 0.001), associated with fSAD% at the lower lobes (p < 0.05). In addition, Koreans had smaller Dh and WT of the trachea (both, p < 0.05), correlated with the forced expiratory volume in 1 second (R = 0.49, 0.39; all p < 0.001) and forced vital capacity (R = 0.55, 0.45; all p < 0.001).

CONCLUSION

Koreans had unique features of airway structure and lung function as compared with Whites, and the difference was clearer in female individuals. Discriminating structural and functional features between Koreans and Whites enables exploration of inter-racial differences of pulmonary disease in terms of severity, distribution, and phenotype.

Pulmonary quantitative CT imaging in focal and diffuse disease: current research and clinical applications

The frenetic development of imaging technology-both hardware and software-provides exceptional potential for investigation of the lung. In the last two decades, CT was exploited for detailed characterization of pulmonary structures and description of respiratory disease. The introduction of volumetric acquisition allowed increasingly sophisticated analysis of CT data by means of computerized algorithm, namely quantitative CT (QCT). Hundreds of thousands of CTs have been analysed for characterization of focal and diffuse disease of the lung. Several QCT metrics were developed and tested against clinical, functional and prognostic descriptors. Computer-aided detection of nodules, textural analysis of focal lesions, densitometric analysis and airway segmentation in obstructive pulmonary disease and textural analysis in interstitial lung disease are the major chapters of this discipline. The validation of QCT metrics for specific clinical and investigational needs prompted the translation of such metrics from research field to patient care. The present review summarizes the state of the art of QCT in both focal and diffuse lung disease, including a dedicated discussion about application of QCT metrics as parameters for clinical care and outcomes in clinical trials.

Computed Tomography Image Matching in Chronic Obstructive Pulmonary Disease

Chronic obstructive pulmonary disease (COPD), characterized by progressive airflow obstruction due to the combined effects of emphysema and small airways disease, is associated with high morbidity and mortality. The complex link between emphysema and airways disease is associated with significant heterogeneity in clinical presentation. Spirometry is the current gold standard for diagnosis and stratification of the severity of airflow obstruction in COPD. Although spirometry is simple to use, it does not enable the separation of emphysema from airways disease. Computed tomography (CT), on the other hand, provides the anatomic localization of disease and has been increasingly used to phenotype COPD. The majority of current CT measures are extracted from a single-volume CT scan and although useful to characterize emphysema and airways disease, they do not link structural and functional abnormalities. Alternatively, CT image matching combines information from both inspiratory and expiratory CT scans, thus enabling determination of functional changes such as regional ventilation and mechanical properties of the lung. In this review, we discuss recent applications of CT image matching that provide clinically meaningful information beyond spirometry and single-volume CT scan measures.