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.

Loss of Airway Phylogenetic Diversity Is Associated with Clinical and Pathobiological Markers of Disease Development in Chronic Obstructive Pulmonary Disease.

Rationale: The airway microbiome has the potential to shape chronic obstructive pulmonary disease (COPD) pathogenesis, but its relationship to outcomes in milder disease is unestablished. Objectives: To identify sputum microbiome characteristics associated with markers of COPD in participants of the Subpopulations and Intermediate Outcome Measures of COPD Study (SPIROMICS). Methods: Sputum DNA from 877 participants was analyzed using 16S ribosomal RNA gene sequencing. Relationships between baseline airway microbiota composition and clinical, radiographic, and mucoinflammatory markers, including longitudinal lung function trajectory, were examined. Measurements and Main Results: Participant data represented predominantly milder disease (Global Initiative for Chronic Obstructive Lung Disease stage 0-2 obstruction in 732 of 877 participants). Phylogenetic diversity (i.e., range of different species within a sample) correlated positively with baseline lung function, decreased with higher Global Initiative for Chronic Obstructive Lung Disease stage, and correlated negatively with symptom burden, radiographic markers of airway disease, and total mucin concentrations (P < 0.001). In covariate-adjusted regression models, organisms robustly associated with better lung function included Alloprevotella, Oribacterium, and Veillonella species. Conversely, lower lung function, greater symptoms, and radiographic measures of small airway disease were associated with enrichment in members of Streptococcus, Actinobacillus, Actinomyces, and other genera. Baseline sputum microbiota features were also associated with lung function trajectory during SPIROMICS follow-up (stable/improved, decline, or rapid decline groups). The stable/improved group (slope of FEV1 regression ⩾66th percentile) had greater bacterial diversity at baseline associated with enrichment in Prevotella, Leptotrichia, and Neisseria species. In contrast, the rapid decline group (FEV1 slope ⩽33rd percentile) had significantly lower baseline diversity associated with enrichment in Streptococcus species. Conclusions: In SPIROMICS, baseline airway microbiota features demonstrate divergent associations with better or worse COPD-related outcomes.

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.

Multicenter evaluation of parametric response mapping as an indicator of bronchiolitis obliterans syndrome after hematopoietic stem cell transplantation.

Parametric response mapping (PRM) is a novel computed tomography (CT) technology that has shown potential for assessment of bronchiolitis obliterans syndrome (BOS) after hematopoietic stem cell transplantation (HCT). The primary aim of this study was to evaluate whether variations in image acquisition under real-world conditions affect the PRM measurements of clinically diagnosed BOS. CT scans were obtained retrospectively from 72 HCT recipients with BOS and graft-versus-host disease from Fred Hutchinson Cancer Research Center, Karolinska Institute, and the University of Michigan. Whole lung volumetric scans were performed at inspiration and expiration using site-specific acquisition and reconstruction protocols. PRM and pulmonary function measurements were assessed. Patients with moderately severe BOS at diagnosis (median forced expiratory volume at 1 second [FEV1] 53.5% predicted) had similar characteristics between sites. Variations in site-specific CT acquisition protocols had a negligible effect on the PRM-derived small airways disease (SAD), that is, BOS measurements. PRM-derived SAD was found to correlate with FEV1% predicted and FEV1/ forced vital capacity (R = -0.236, P = .046; and R = -0.689, P < .0001, respectively), which suggests that elevated levels in the PRM measurements are primarily affected by BOS airflow obstruction and not CT scan acquisition parameters. Based on these results, PRM may be applied broadly for post-HCT diagnosis and monitoring of BOS.

Noninvasive Imaging Biomarker Identifies Small Airway Damage in Severe Chronic Obstructive Pulmonary Disease.

Rationale: Evidence suggests damage to small airways is a key pathologic lesion in chronic obstructive pulmonary disease (COPD). Computed tomography densitometry has been demonstrated to identify emphysema, but no such studies have been performed linking an imaging metric to small airway abnormality.Objectives: To correlate ex vivo parametric response mapping (PRM) analysis to in vivo lung tissue measurements of patients with severe COPD treated by lung transplantation and control subjects.Methods: Resected lungs were inflated, frozen, and systematically sampled, generating 33 COPD (n = 11 subjects) and 22 control tissue samples (n = 3 subjects) for micro-computed tomography analysis of terminal bronchioles (TBs; last generation of conducting airways) and emphysema.Measurements and Main Results: PRM analysis was conducted to differentiate functional small airways disease (PRMfSAD) from emphysema (PRMEmph). In COPD lungs, TB numbers were reduced (P = 0.01); surviving TBs had increased wall area percentage (P < 0.001), decreased circularity (P < 0.001), reduced cross-sectional luminal area (P < 0.001), and greater airway obstruction (P = 0.008). COPD lungs had increased airspace size (P < 0.001) and decreased alveolar surface area (P < 0.001). Regression analyses demonstrated unique correlations between PRMfSAD and TBs, with decreased circularity (P < 0.001), decreased luminal area (P < 0.001), and complete obstruction (P = 0.008). PRMEmph correlated with increased airspace size (P < 0.001), decreased alveolar surface area (P = 0.003), and fewer alveolar attachments per TB (P = 0.01).Conclusions: PRMfSAD identifies areas of lung tissue with TB loss, luminal narrowing, and obstruction. This is the first confirmation that an imaging biomarker can identify terminal bronchial pathology in established COPD and provides a noninvasive imaging methodology to identify small airway damage in COPD.

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.

A Deep Learning-Based Radiomic Classifier for Usual Interstitial Pneumonia.

BACKGROUND: Because chest CT scan has largely supplanted surgical lung biopsy for diagnosing most cases of interstitial lung disease (ILD), tools to standardize CT scan interpretation are urgently needed. RESEARCH QUESTION: Does a deep learning (DL)-based classifier for usual interstitial pneumonia (UIP) derived using CT scan features accurately discriminate radiologist-determined visual UIP? STUDY DESIGN AND METHODS: A retrospective cohort study was performed. Chest CT scans acquired in individuals with and without ILD were drawn from a variety of public and private data sources. Using radiologist-determined visual UIP as ground truth, a convolutional neural network was used to learn discrete CT scan features of UIP, with outputs used to predict the likelihood of UIP using a linear support vector machine. Test performance characteristics were assessed in an independent performance cohort and multicenter ILD clinical cohort. Transplant-free survival was compared between UIP classification approaches using the Kaplan-Meier estimator and Cox proportional hazards regression. RESULTS: A total of 2,907 chest CT scans were included in the training (n = 1,934), validation (n = 408), and performance (n = 565) data sets. The prevalence of radiologist-determined visual UIP was 12.4% and 37.1% in the performance and ILD clinical cohorts, respectively. The DL-based UIP classifier predicted visual UIP in the performance cohort with sensitivity and specificity of 93% and 86%, respectively, and in the multicenter ILD clinical cohort with 81% and 77%, respectively. DL-based and visual UIP classification similarly discriminated survival, and outcomes were consistent among cases with positive DL-based UIP classification irrespective of visual classification. INTERPRETATION: A DL-based classifier for UIP demonstrated good test performance across a wide range of UIP prevalence and similarly discriminated survival when compared with radiologist-determined UIP. This automated tool could efficiently screen for UIP in patients undergoing chest CT scan and identify a high-risk phenotype among those with known ILD.

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.

Machine learning for screening of at-risk, mild and moderate COPD patients at risk of FEV1 decline: results from COPDGene and SPIROMICS.

Purpose: The purpose of this study was to train and validate machine learning models for predicting rapid decline of forced expiratory volume in 1 s (FEV1) in individuals with a smoking history at-risk-for chronic obstructive pulmonary disease (COPD), Global Initiative for Chronic Obstructive Lung Disease (GOLD 0), or with mild-to-moderate (GOLD 1-2) COPD. We trained multiple models to predict rapid FEV1 decline using demographic, clinical and radiologic biomarker data. Training and internal validation data were obtained from the COPDGene study and prediction models were validated against the SPIROMICS cohort. Methods: We used GOLD 0-2 participants (n = 3,821) from COPDGene (60.0 ± 8.8 years, 49.9% male) for variable selection and model training. Accelerated lung function decline was defined as a mean drop in FEV1% predicted of > 1.5%/year at 5-year follow-up. We built logistic regression models predicting accelerated decline based on 22 chest CT imaging biomarker, pulmonary function, symptom, and demographic features. Models were validated using n = 885 SPIROMICS subjects (63.6 ± 8.6 years, 47.8% male). Results: The most important variables for predicting FEV1 decline in GOLD 0 participants were bronchodilator responsiveness (BDR), post bronchodilator FEV1% predicted (FEV1.pp.post), and CT-derived expiratory lung volume; among GOLD 1 and 2 subjects, they were BDR, age, and PRMlower lobes fSAD. In the validation cohort, GOLD 0 and GOLD 1-2 full variable models had significant predictive performance with AUCs of 0.620 ± 0.081 (p = 0.041) and 0.640 ± 0.059 (p < 0.001). Subjects with higher model-derived risk scores had significantly greater odds of FEV1 decline than those with lower scores. Conclusion: Predicting FEV1 decline in at-risk patients remains challenging but a combination of clinical, physiologic and imaging variables provided the best performance across two COPD cohorts.

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.

A Framework for Evaluating the Technical Performance of Multiparameter Quantitative Imaging Biomarkers (mp-QIBs).

Multiparameter quantitative imaging incorporates anatomical, functional, and/or behavioral biomarkers to characterize tissue, detect disease, identify phenotypes, define longitudinal change, or predict outcome. Multiple imaging parameters are sometimes considered separately but ideally are evaluated collectively. Often, they are transformed as Likert interpretations, ignoring the correlations of quantitative properties that may result in better reproducibility or outcome prediction. In this paper we present three use cases of multiparameter quantitative imaging: i) multidimensional descriptor, ii) phenotype classification, and iii) risk prediction. A fourth application based on data-driven markers from radiomics is also presented. We describe the technical performance characteristics and their metrics common to all use cases, and provide a structure for the development, estimation, and testing of multiparameter quantitative imaging. This paper serves as an overview for a series of individual articles on the four applications, providing the statistical framework for multiparameter imaging applications in medicine.

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.