Temporal Exploration of COPD Phenotypes: Insights from the COPDGene and SPIROMICS Cohorts.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) exhibits considerable progression heterogeneity. We hypothesized that elastic principal graph analysis (EPGA) would identify distinct clinical phenotypes and their longitudinal relationships. METHODS: Cross-sectional data from 8,972 tobacco-exposed COPDGene participants, with and without COPD, were used to train a model with EPGA, using thirty clinical, physiologic and CT features. Principal component analysis (PCA) was used to reduce data dimensionality to six principal components. An elastic principal tree was fitted to the reduced space. 4,585 participants from COPDGene Phase 2 were used to test longitudinal trajectories. 2,652 participants from SPIROMICS tested external reproducibility. RESULTS: Our analysis used cross-sectional data to create an elastic principal tree, where the concept of time is represented by distance on the tree. Six clinically distinct tree segments were identified that differed by lung function, symptoms, and CT features: 1) Subclinical (SC); 2) Parenchymal Abnormality (PA); 3) Chronic Bronchitis (CB); 4) Emphysema Male (EM); 5) Emphysema Female (EF); and 6) Severe Airways (SA) disease. Cross-sectional SPIROMICS data confirmed similar groupings. 5-year data from COPDGene mapped longitudinal changes onto the tree. 29% of patients changed segment during follow-up; longitudinal trajectories confirmed a net flow of patients along the tree, from SC towards Emphysema, although alternative trajectories were noted, through airway disease predominant phenotypes, CB and SA. CONCLUSION: This novel analytic methodology provides an approach to defining longitudinal phenotypic trajectories using cross sectional data. These insights are clinically relevant and could facilitate precision therapy and future trials to modify disease progression.

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.

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.

Relationship between diffusion capacity and small airway abnormality in COPDGene.

Impaired single breath carbon monoxide diffusing capacity (DLCO) is associated with emphysema. Small airways disease (SAD) may be a precursor lesion to emphysema, but the relationship between SAD and DLCO is undescribed. We hypothesized that in mild COPD, functional SAD (fSAD) defined by computed tomography (CT) and Parametric Response Mapping methodology would correlate with impaired DLCO. Using data from ever-smokers in the COPDGene cohort, we established that fSAD correlated significantly with lower DLCO among both non-obstructed and GOLD 1-2 subjects. The relationship between DLCO with CT-defined emphysema was present in all GOLD stages, but most prominent in severe disease. TRIAL REGISTRATION: NCT00608764. Registry: COPDGene. Registered 06 February 2008, retrospectively registered.

Topologic Parametric Response Mapping Identifies Tissue Subtypes Associated with Emphysema Progression.

RATIONALE AND OBJECTIVES: Small airways disease (SAD) and emphysema are significant components of chronic obstructive pulmonary disease (COPD), a heterogenous disease where predicting progression is difficult. SAD, a principal cause of airflow obstruction in mild COPD, has been identified as a precursor to emphysema. Parametric Response Mapping (PRM) of chest computed tomography (CT) can help distinguish SAD from emphysema. Specifically, topologic PRM can define local patterns of both diseases to characterize how and in whom COPD progresses. We aimed to determine if distribution of CT-based PRM of functional SAD (fSAD) is associated with emphysema progression. MATERIALS AND METHODS: We analyzed paired inspiratory-expiratory chest CT scans at baseline and 5-year follow up in 1495 COPDGene subjects using topological analyses of PRM classifications. By spatially aligning temporal scans, we mapped local emphysema at year five to baseline lobar PRM-derived topological readouts. K-means clustering was applied to all observations. Subjects were subtyped based on predominant PRM cluster assignments and assessed using non-parametric statistical tests to determine differences in PRM values, pulmonary function metrics, and clinical measures. RESULTS: We identified distinct lobar imaging patterns and classified subjects into three radiologic subtypes: emphysema-dominant (ED), fSAD-dominant (FD), and fSAD-transition (FT: transition from healthy lung to fSAD). Relative to year five emphysema, FT showed rapid local emphysema progression (-57.5% ± 1.1) compared to FD (-49.9% ± 0.5) and ED (-33.1% ± 0.4). FT consisted primarily of at-risk subjects (roughly 60%) with normal spirometry. CONCLUSION: The FT subtype of COPD may allow earlier identification of individuals without spirometrically-defined COPD at-risk for developing emphysema.

Fully automated pipeline for measurement of the thoracic aorta using joint segmentation and localization neural network.

Purpose: Diagnosis and surveillance of thoracic aortic aneurysm (TAA) involves measuring the aortic diameter at various locations along the length of the aorta, often using computed tomography angiography (CTA). Currently, measurements are performed by human raters using specialized software for three-dimensional analysis, a time-consuming process, requiring 15 to 45 min of focused effort. Thus, we aimed to develop a convolutional neural network (CNN)-based algorithm for fully automated and accurate aortic measurements. Approach: Using 212 CTA scans, we trained a CNN to perform segmentation and localization of key landmarks jointly. Segmentation mask and landmarks are subsequently used to obtain the centerline and cross-sectional diameters of the aorta. Subsequently, a cubic spline is fit to the aortic boundary at the sinuses of Valsalva to avoid errors related inclusions of coronary artery origins. Performance was evaluated on a test set of 60 scans with automated measurements compared against expert manual raters. Result: Compared to training separate networks for each task, joint training yielded higher accuracy for segmentation, especially at the boundary (p<0.001), but a marginally worse (0.2 to 0.5 mm) accuracy for landmark localization (p<0.001). Mean absolute error between human and automated was ≤1 mm at six of nine standard clinical measurement locations. However, higher errors were noted in the aortic root and arch regions, ranging between 1.4 and 2.2 mm, although agreement of manual raters was also lower in these regions. Conclusion: Fully automated aortic diameter measurements in TAA are feasible using a CNN-based algorithm. Automated measurements demonstrated low errors that are comparable in magnitude to those with manual raters; however, measurement error was highest in the aortic root and arch.

Comparison of hyperpolarized 3He-MRI, CT based parametric response mapping, and mucus scores in asthmatics.

Purpose: The purpose of this study was to anatomically correlate ventilation defects with regions of air trapping by whole lung, lung lobe, and airway segment in the context of airway mucus plugging in asthma. Methods: A total of 34 asthmatics [13M:21F, 13 mild/moderate, median age (range) of 49.5 (36.8-53.3) years and 21 severe, 56.1 (47.1-62.6) years] and 4 healthy subjects [1M:3F, 38.5 (26.6-52.2) years] underwent HP 3He MRI and CT imaging. HP 3He MRI was assessed for ventilation defects using a semi-automated k-means clustering algorithm. Inspiratory and expiratory CTs were analyzed using parametric response mapping (PRM) to quantify markers of emphysema and functional small airways disease (fSAD). Segmental and lobar lung masks were obtained from CT and registered to HP 3He MRI in order to localize ventilation defect percent (VDP), at the lobar and segmental level, to regions of fSAD and mucus plugging. Spearman's correlation was utilized to compare biomarkers on a global and lobar level, and a multivariate analysis was conducted to predict segmental fSAD given segmental VDP (sVDP) and mucus score as variables in order to further understand the functional relationships between regional measures of obstruction. Results: On a global level, fSAD was correlated with whole lung VDP (r = 0.65, p < 0.001), mucus score (r = 0.55, p < 0.01), and moderately correlated (-0.60 ≤ r ≤ -0.56, p < 0.001) to percent predicted (%p) FEV1, FEF25-75 and FEV1/FVC, and more weakly correlated to FVC%p (-0.38 ≤ r ≤ -0.35, p < 0.001) as expected from previous work. On a regional level, lobar VDP, mucus scores, and fSAD were also moderately correlated (r from 0.45-0.66, p < 0.01). For segmental colocalization, the model of best fit was a piecewise quadratic model, which suggests that sVDP may be increasing due to local airway obstruction that does not manifest as fSAD until more extensive disease is present. sVDP was more sensitive to the presence of a mucus plugs overall, but the prediction of fSAD using multivariate regression showed an interaction in the presence of a mucus plugs when sVDP was between 4% and 10% (p < 0.001). Conclusion: This multi-modality study in asthma confirmed that areas of ventilation defects are spatially correlated with air trapping at the level of the airway segment and suggests VDP and fSAD are sensitive to specific sources of airway obstruction in asthma, including mucus plugs.

Lung cancer lesion detection in histopathology images using graph-based sparse PCA network.

Early detection of lung cancer is critical for improvement of patient survival. To address the clinical need for efficacious treatments, genetically engineered mouse models (GEMM) have become integral in identifying and evaluating the molecular underpinnings of this complex disease that may be exploited as therapeutic targets. Assessment of GEMM tumor burden on histopathological sections performed by manual inspection is both time consuming and prone to subjective bias. Therefore, an interplay of needs and challenges exists for computer-aided diagnostic tools, for accurate and efficient analysis of these histopathology images. In this paper, we propose a simple machine learning approach called the graph-based sparse principal component analysis (GS-PCA) network, for automated detection of cancerous lesions on histological lung slides stained by hematoxylin and eosin (H&E). Our method comprises four steps: 1) cascaded graph-based sparse PCA, 2) PCA binary hashing, 3) block-wise histograms, and 4) support vector machine (SVM) classification. In our proposed architecture, graph-based sparse PCA is employed to learn the filter banks of the multiple stages of a convolutional network. This is followed by PCA hashing and block histograms for indexing and pooling. The meaningful features extracted from this GS-PCA are then fed to an SVM classifier. We evaluate the performance of the proposed algorithm on H&E slides obtained from an inducible K-rasG12D lung cancer mouse model using precision/recall rates, Fß-score, Tanimoto coefficient, and area under the curve (AUC) of the receiver operator characteristic (ROC) and show that our algorithm is efficient and provides improved detection accuracy compared to existing algorithms.

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.

Validation of a robust method for quantification of three-dimensional growth of the thoracic aorta using deformable image registration.

PURPOSE: Accurate assessment of thoracic aortic aneurysm (TAA) growth is important for appropriate clinical management. Maximal aortic diameter is the primary metric that is used to assess growth, but it suffers from substantial measurement variability. A recently proposed technique, termed vascular deformation mapping (VDM), is able to quantify three-dimensional aortic growth using clinical computed tomography angiography (CTA) data using an approach based on deformable image registration (DIR). However, the accuracy and robustness of VDM remains undefined given the lack of ground truth from clinical CTA data, and, furthermore, the performance of VDM relative to standard manual diameter measurements is unknown. METHODS: To evaluate the performance of the VDM pipeline for quantifying aortic growth, we developed a novel and systematic evaluation process to generate 76 unique synthetic CTA growth phantoms (based on 10 unique cases) with variable degrees and locations of aortic wall deformation. Aortic deformation was quantified using two metrics: area ratio (AR), defined as the ratio of surface area in triangular mesh elements and the magnitude of deformation in the normal direction (DiN) relative to the aortic surface. Using these phantoms, we further investigated the effects on VDM's measurement accuracy resulting from factors that influence the quality of clinical CTA data such as respiratory translations, slice thickness, and image noise. Lastly, we compare the measurement error of VDM TAA growth assessments against two expert raters performing standard diameter measurements of synthetic phantom images. RESULTS: Across our population of 76 synthetic growth phantoms, the median absolute error was 0.063 (IQR: 0.073-0.054) for AR and 0.181 mm (interquartile range [IQR]: 0.214-0.143 mm) for DiN. Median relative error was 1.4% for AR and 3.3 % $3.3\%$ for DiN at the highest tested noise level (contrast-to-noise ratio [CNR] = 2.66). Error in VDM output increased with slice thickness, with the highest median relative error of 1.5% for AR and 4.1% for DiN at a slice thickness of 2.0 mm. Respiratory motion of the aorta resulted in maximal absolute error of 3% AR and 0.6 mm in DiN, but bulk translations in aortic position had a very small effect on measured AR and DiN values (relative errors < 1 % $< 1\%$ ). VDM-derived measurements of magnitude and location of maximal diameter change demonstrated significantly high accuracy and lower variability compared to two expert manual raters ( p < 0.03 $p<0.03$ across all comparisons). CONCLUSIONS: VDM yields an accurate, three-dimensional assessment of aortic growth in TAA patients and is robust to factors such as image noise, respiration-induced translations, and differences in patient position. Further, VDM significantly outperformed two expert manual raters in assessing the magnitude and location of aortic growth despite optimized experimental measurement conditions. These results support validation of the VDM technique for accurate quantification of aortic growth in patients and highlight several important advantages over diameter measurements.

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.

Comparison of CT Lung Density Measurements between Standard Full-Dose and Reduced-Dose Protocols.

PURPOSE: To evaluate the reproducibility and predicted clinical outcomes of CT-based quantitative lung density measurements using standard fixed-dose (FD) and reduced-dose (RD) scans. MATERIALS AND METHODS: In this retrospective analysis of prospectively acquired data, 1205 participants (mean age, 65 years ± 9 [standard deviation]; 618 men) enrolled in the COPDGene study who underwent FD and RD CT image acquisition protocols between November 2014 and July 2017 were included. Of these, the RD scans of 640 participants were also reconstructed using iterative reconstruction (IR). Median filtering was applied to the RD scans (RD-MF) to investigate an alternative noise reduction strategy. CT attenuation at the 15th percentile of the lung CT histogram (Perc15) was computed for all image types (FD, RD, RD-MF, and RD-IR). Reproducibility coefficients were calculated to quantify the measurement differences between FD and RD scans. The ability of Perc15 to predict chronic obstructive pulmonary disease (COPD) diagnosis and exacerbation frequency was investigated using receiver operating characteristic analysis. RESULTS: The Perc15 reproducibility coefficients with and without volume adjustment were as follows: RD, 29.43 HU ± 0.62 versus 32.81 HU ± 1.70; RD-MF, 7.42 HU ± 0.42 versus 19.40 HU ± 2.65; and RD-IR, 7.10 HU ± 0.52 versus 22.46 HU ± 3.91. Receiver operating characteristic curve analysis indicated that Perc15 on volume-adjusted FD and RD scans were both predictive for COPD diagnosis (area under the receiver operating characteristic curve [AUC]: FD, 0.724 ± 0.045; RD, 0.739 ± 0.045) and for having one or more exacerbation per year (AUCs: FD, 0.593 ± 0.068; RD, 0.589 ± 0.066). Similar trends were observed when volume adjustment was not applied. CONCLUSION: A combination of volume adjustment and noise reduction filtering improved the reproducibility of lung density measurements computed using serial FD and RD CT scans.Supplemental material is available for this article.© RSNA, 2021.

Improved detection of air trapping on expiratory computed tomography using deep learning.

BACKGROUND: Radiologic evidence of air trapping (AT) on expiratory computed tomography (CT) scans is associated with early pulmonary dysfunction in patients with cystic fibrosis (CF). However, standard techniques for quantitative assessment of AT are highly variable, resulting in limited efficacy for monitoring disease progression. OBJECTIVE: To investigate the effectiveness of a convolutional neural network (CNN) model for quantifying and monitoring AT, and to compare it with other quantitative AT measures obtained from threshold-based techniques. MATERIALS AND METHODS: Paired volumetric whole lung inspiratory and expiratory CT scans were obtained at four time points (0, 3, 12 and 24 months) on 36 subjects with mild CF lung disease. A densely connected CNN (DN) was trained using AT segmentation maps generated from a personalized threshold-based method (PTM). Quantitative AT (QAT) values, presented as the relative volume of AT over the lungs, from the DN approach were compared to QAT values from the PTM method. Radiographic assessment, spirometric measures, and clinical scores were correlated to the DN QAT values using a linear mixed effects model. RESULTS: QAT values from the DN were found to increase from 8.65% ± 1.38% to 21.38% ± 1.82%, respectively, over a two-year period. Comparison of CNN model results to intensity-based measures demonstrated a systematic drop in the Dice coefficient over time (decreased from 0.86 ± 0.03 to 0.45 ± 0.04). The trends observed in DN QAT values were consistent with clinical scores for AT, bronchiectasis, and mucus plugging. In addition, the DN approach was found to be less susceptible to variations in expiratory deflation levels than the threshold-based approach. CONCLUSION: The CNN model effectively delineated AT on expiratory CT scans, which provides an automated and objective approach for assessing and monitoring AT in CF patients.

Quantitative Emphysema on Low-Dose CT Imaging of the Chest and Risk of Lung Cancer and Airflow Obstruction: An Analysis of the National Lung Screening Trial.

BACKGROUND: Lung cancer risk prediction models do not routinely incorporate imaging metrics available on low-dose CT (LDCT) imaging of the chest ordered for lung cancer screening. RESEARCH QUESTION: What is the association between quantitative emphysema measured on LDCT imaging and lung cancer incidence and mortality, all-cause mortality, and airflow obstruction in individuals who currently or formerly smoked and are undergoing lung cancer screening? STUDY DESIGN AND METHODS: In 7,262 participants in the CT arm of the National Lung Screening Trial, percent low attenuation area (%LAA) was defined as the percentage of lung volume with voxels less than -950 Hounsfield units on the baseline examination. Multivariable Cox proportional hazards models, adjusting for competing risks where appropriate, were built to test for association between %LAA and lung cancer incidence, lung cancer mortality, and all-cause mortality with censoring at 6 years. In addition, multivariable logistic regression models were built to test the cross-sectional association between %LAA and airflow obstruction on spirometry, which was available in 2,700 participants. RESULTS: The median %LAA was 0.8% (interquartile range, 0.2%-2.7%). Every 1% increase in %LAA was independently associated with higher hazards of lung cancer incidence (hazard ratio [HR], 1.02; 95% CI, 1.01-1.03; P = .004), lung cancer mortality (HR, 1.02; 95% CI, 1.00-1.05; P = .045), and all-cause mortality (HR, 1.01; 95% CI, 1.00-1.03; P = .042). Among participants with spirometry, 892 had airflow obstruction. The likelihood of airflow obstruction increased with every 1% increase in %LAA (odds ratio, 1.07; 95% CI, 1.06-1.09; P < .001). A %LAA cutoff of 1% had the best discriminative accuracy for airflow obstruction in participants aged > 65 years. INTERPRETATION: Quantitative emphysema measured on LDCT imaging of the chest can be leveraged to improve lung cancer risk prediction and help diagnose COPD in individuals who currently or formerly smoked and are undergoing lung cancer screening.

Progression of Emphysema and Small Airways Disease in Cigarette Smokers.

BACKGROUND: Little is known about factors associated with emphysema progression in cigarette smokers. We evaluated factors associated with change in emphysema and forced expiratory volume in 1 second (FEV1) in participants with and without chronic obstructive pulmonary disease (COPD). METHODS: This retrospective study included individuals participating in the COPD Genetic Epidemiology study who completed the 5-year follow-up, including inspiratory and expiratory computed tomography (CT) and spirometry. All paired CT scans were analyzed using micro-mapping, which classifies individual voxels as emphysema or functional small airway disease (fSAD). Presence and progression of emphysema and FEV1 were determined based on comparison to nonsmoker values. Logistic regression analyses were used to identify clinical parameters associated with disease progression. RESULTS: A total of 3088 participants were included with a mean ± SD age of 60.7±8.9 years, including 72 nonsmokers. In all Global initiative for chronic Obstructive Lung Disease (GOLD) stages, the presence of emphysema at baseline was associated with emphysema progression (odds ratio [OR]: GOLD 0: 4.32; preserved ratio-impaired spirometry [PRISm]; 5.73; GOLD 1: 5.16; GOLD 2: 5.69; GOLD 3/4: 5.55; all p ≤0.01). If there was no emphysema at baseline, the amount of fSAD at baseline was associated with emphysema progression (OR for 1% increase: GOLD 0: 1.06; PRISm: 1.20; GOLD 1: 1.7; GOLD 3/4: 1.08; all p ≤ 0.03).In 1735 participants without spirometric COPD, progression in emphysema occurred in 105 (6.1%) participants and only 21 (1.2%) had progression in both emphysema and FEV1. CONCLUSIONS: The presence of emphysema is an important predictor of emphysema progression. In patients without emphysema, fSAD is associated with the development of emphysema. In participants without spirometric COPD, emphysema progression occurred independently of FEV1 decline.

Reprint of: Voxel-Wise Longitudinal Parametric Response Mapping Analysis of Chest Computed Tomography in Smokers.

RATIONALE AND OBJECTIVES: Chronic obstructive pulmonary disease is a heterogeneous disease characterized by small airway abnormality and emphysema. We hypothesized that a voxel-wise computed tomography analytic approach would identify patterns of disease progression in smokers. MATERIALS AND METHODS: We analyzed 725 smokers in spirometric GOLD stages 0-4 with two chest CTs 5 years apart. Baseline inspiration, follow-up inspiration and follow-up expiration images were spatially registered to baseline expiration so that each voxel had correspondences across all time points and respiratory phases. Voxel-wise Parametric Response Mapping (PRM) was then generated for the baseline and follow-up scans. PRM classifies lung as normal, functional small airway disease (PRMfSAD), and emphysema (PRMEMPH). RESULTS: Subjects with low baseline PRMfSAD and PRMEMPH predominantly had an increase in PRMfSAD on follow-up; those with higher baseline PRMfSAD and PRMEMPH mostly had increases in PRMEMPH. For GOLD 0 participants (n = 419), mean 5-year increases in PRMfSAD and PRMEMPH were 0.3% for both; for GOLD 1-4 participants (n = 306), they were 0.6% and 1.6%, respectively. Eighty GOLD 0 subjects (19.1%) had overall radiologic progression (30.0% to PRMfSAD, 52.5% to PRMEMPH, and 17.5% to both); 153 GOLD 1-4 subjects (50.0%) experienced progression (17.6% to PRMfSAD, 48.4% to PRMEMPH, and 34.0% to both). In a multivariable model, both baseline PRMfSAD and PRMEMPH were associated with development of PRMEMPH on follow-up, although this relationship was diminished at higher levels of baseline PRMEMPH. CONCLUSION: A voxel-wise longitudinal PRM analytic approach can identify patterns of disease progression in smokers with and without chronic obstructive pulmonary disease.

Voxel-Wise Longitudinal Parametric Response Mapping Analysis of Chest Computed Tomography in Smokers.

RATIONALE AND OBJECTIVES: Chronic obstructive pulmonary disease is a heterogeneous disease characterized by small airway abnormality and emphysema. We hypothesized that a voxel-wise computed tomography analytic approach would identify patterns of disease progression in smokers. MATERIALS AND METHODS: We analyzed 725 smokers in spirometric GOLD stages 0-4 with two chest CTs 5 years apart. Baseline inspiration, follow-up inspiration and follow-up expiration images were spatially registered to baseline expiration so that each voxel had correspondences across all time points and respiratory phases. Voxel-wise Parametric Response Mapping (PRM) was then generated for the baseline and follow-up scans. PRM classifies lung as normal, functional small airway disease (PRMfSAD), and emphysema (PRMEMPH). RESULTS: Subjects with low baseline PRMfSAD and PRMEMPH predominantly had an increase in PRMfSAD on follow-up; those with higher baseline PRMfSAD and PRMEMPH mostly had increases in PRMEMPH. For GOLD 0 participants (n = 419), mean 5-year increases in PRMfSAD and PRMEMPH were 0.3% for both; for GOLD 1-4 participants (n = 306), they were 0.6% and 1.6%, respectively. Eighty GOLD 0 subjects (19.1%) had overall radiologic progression (30.0% to PRMfSAD, 52.5% to PRMEMPH, and 17.5% to both); 153 GOLD 1-4 subjects (50.0%) experienced progression (17.6% to PRMfSAD, 48.4% to PRMEMPH, and 34.0% to both). In a multivariable model, both baseline PRMfSAD and PRMEMPH were associated with development of PRMEMPH on follow-up, although this relationship was diminished at higher levels of baseline PRMEMPH. CONCLUSION: A voxel-wise longitudinal PRM analytic approach can identify patterns of disease progression in smokers with and without chronic obstructive pulmonary disease.

A dynamic model-based approach to motion and deformation tracking of prosthetic valves from biplane x-ray images.

PURPOSE: Transcatheter aortic valve replacement (TAVR) is a minimally invasive procedure in which a prosthetic heart valve is placed and expanded within a defective aortic valve. The device placement is commonly performed using two-dimensional (2D) fluoroscopic imaging. Within this work, we propose a novel technique to track the motion and deformation of the prosthetic valve in three dimensions based on biplane fluoroscopic image sequences. METHODS: The tracking approach uses a parameterized point cloud model of the valve stent which can undergo rigid three-dimensional (3D) transformation and different modes of expansion. Rigid elements of the model are individually rotated and translated in three dimensions to approximate the motions of the stent. Tracking is performed using an iterative 2D-3D registration procedure which estimates the model parameters by minimizing the mean-squared image values at the positions of the forward-projected model points. Additionally, an initialization technique is proposed, which locates clusters of salient features to determine the initial position and orientation of the model. RESULTS: The proposed algorithms were evaluated based on simulations using a digital 4D CT phantom as well as experimentally acquired images of a prosthetic valve inside a chest phantom with anatomical background features. The target registration error was 0.12 ± 0.04 mm in the simulations and 0.64 ± 0.09 mm in the experimental data. CONCLUSIONS: The proposed algorithm could be used to generate 3D visualization of the prosthetic valve from two projections. In combination with soft-tissue sensitive-imaging techniques like transesophageal echocardiography, this technique could enable 3D image guidance during TAVR procedures.

Single-view geometric calibration for C-arm inverse geometry CT.

Accurate and artifact-free reconstruction of tomographic images requires precise knowledge of the imaging system geometry. A projection matrix-based calibration method to enable C-arm inverse geometry CT (IGCT) is proposed. The method is evaluated for scanning-beam digital x-ray (SBDX), a C-arm mounted inverse geometry fluoroscopic technology. A helical configuration of fiducials is imaged at each gantry angle in a rotational acquisition. For each gantry angle, digital tomosynthesis is performed at multiple planes and a composite image analogous to a cone-beam projection is generated from the plane stack. The geometry of the C-arm, source array, and detector array is determined at each angle by constructing a parameterized three-dimensional-to-two-dimensional projection matrix that minimizes the sum-of-squared deviations between measured and projected fiducial coordinates. Simulations were used to evaluate calibration performance with translations and rotations of the source and detector. The relative root-mean-square error in a reconstruction of a numerical thorax phantom was 0.4% using the calibration method versus 7.7% without calibration. In phantom studies, reconstruction of SBDX projections using the proposed method eliminated artifacts present in noncalibrated reconstructions. The proposed IGCT calibration method reduces image artifacts when uncertainties exist in system geometry.

Correcting Nonpathological Variation in Longitudinal Parametric Response Maps of CT Scans in COPD Subjects: SPIROMICS.

Small airways disease (SAD) is one of the leading causes of airflow limitations in patients diagnosed with chronic obstructive pulmonary disease (COPD). Parametric response mapping (PRM) of computed tomography (CT) scans allows for the quantification of this previously invisible COPD component. Although PRM is being investigated as a diagnostic tool for COPD, variability in the longitudinal measurements of SAD by PRM has been reported. Here, we show a method for correcting longitudinal PRM data because of non-pathological variations in serial CT scans. In this study, serial whole-lung high-resolution CT scans over a 30-day interval were obtained from 90 subjects with and without COPD accrued as part of SPIROMICS. It was assumed in all subjects that the COPD did not progress between examinations. CT scans were acquired at inspiration and expiration, spatially aligned to a single geometric frame, and analyzed using PRM. By modeling variability in longitudinal CT scans, our method could identify, at the voxel-level, shifts in PRM classification over the 30-day interval. In the absence of any correction, PRM generated serial percent volumes of functional SAD with differences as high as 15%. Applying the correction strategy significantly mitigated this effect with differences ~1%. At the voxel-level, significant differences were found between baseline PRM classifications and the follow-up map computed with and without correction (P <. 01 over GOLD). This strategy of accounting for nonpathological sources of variability in longitudinal PRM may improve the quantification of COPD phenotypes transitioning with disease progression.