Chest MRI and CT Predictors of 10-Year All-Cause Mortality in COPD

Abstract Pulmonary imaging measurements using magnetic resonance imaging (MRI) and computed tomography (CT) have the potential to deepen our understanding of chronic obstructive pulmonary disease (COPD) by measuring airway and parenchymal pathologic information that cannot be provided by spirometry. Currently, MRI and CT measurements are not included in mortality risk predictions, diagnosis, or COPD staging. We evaluated baseline pulmonary function, MRI and CT measurements alongside imaging texture-features to predict 10-year all-cause mortality in ex-smokers with (n = 93; 31 females; 70 ± 9years) and without (n = 69; 29 females, 69 ± 9years) COPD. CT airway and vessel measurements, helium-3 (3He) MRI ventilation defect percent (VDP) and apparent diffusion coefficients (ADC) were quantified. MRI and CT texture-features were extracted using PyRadiomics (version2.2.0). Associations between 10-year all-cause mortality and all clinical and imaging measurements were evaluated using multivariable regression model odds-ratios. Machine-learning predictive models for 10-year all-cause mortality were evaluated using area-under-receiver-operator-characteristic-curve (AUC), sensitivity and specificity analyses. DLCO (%pred) (HR = 0.955, 95%CI: 0.934-0.976, p < 0.001), MRI ADC (HR = 1.843, 95%CI: 1.260-2.871, p < 0.001), and CT informational-measure-of-correlation (HR = 3.546, 95% CI: 1.660-7.573, p = 0.001) were the strongest predictors of 10-year mortality. A machine-learning model trained on clinical, imaging, and imaging textures was the best predictive model (AUC = 0.82, sensitivity = 83%, specificity = 84%) and outperformed the solely clinical model (AUC = 0.76, sensitivity = 77%, specificity = 79%). In ex-smokers, regardless of COPD status, addition of CT and MR imaging texture measurements to clinical models provided unique prognostic information of mortality risk that can allow for better clinical management.Clinical Trial Registration: www.clinicaltrials.gov NCT02279329


Introduction
Chronic obstructive pulmonary disease (COPD) is among the leading causes of global mortality [1].Currently, mortality risk is predicted using a variety of clinical-based models [2][3][4], the most common being the BODE index [5] (Body-mass-index, Obstruction measured via forced-expirator y-volume-in-one-second (FEV 1 ), Dyspnea-score and Exercise limitation measured via six-minute-walk-distance (6MWD)).A meta-analysis comparing various prognostic clinical models showed that compared to BODE, models based on the Age-Dyspnea-score-FEV 1 (ADO) index [4] may be stronger, but these were not statistically significantly different at predicting three-year survival [6].Decades after the predictions of the Fletcher-Peto model [7], the spirometry measurement of FEV 1 remains the clinical measurement of global lung function that helps diagnose and stratify COPD severity [8][9][10].
The small airways are considered the major site of airflow limitation in COPD [11] and spirometry measured at the mouth is not sensitive to small airway measurements.However, currently chest imaging measurements are not included in clinically-accepted mortality risk assessments, diagnosis, prognosis, nor staging of COPD.Chest x-ray computed tomography (CT) provides quantitative measurements of airway wall thinning [12], luminal narrowing and obliteration [13], and parenchymal measurements of gas-trapping and terminal airspace enlargement or emphysema [14].
Recent studies stemming from the COPD Genetic Epidemiology study (COPDGene) cohort showed that in ever-smokers with emphysema, emphysema progression over 5-years was associated with all-cause mortality [15].Emphysema also worsened more quickly in COPD patients with preexisting emphysema who continued to smoke [16].In addition, decreased number of distal vascular branches (i.e.pruning) on CT was associated with an increased mortality risk in both COPD and healthy adults [17].CT measurements of ex vivo lung tissue cores with micro-CT also suggested that COPD may initiate in the small airways and that small airway abnormalities precede the development of emphysema and airflow obstruction [18][19][20][21].
Qualitative visual CT scoring by radiologists shows a stronger association with both pulmonary function and mortality than standard quantitative CT measurements [22,23].This could, in part, be explained by visual estimates of emphysema describing both decreased tissue density and complexity of emphysema distribution (and the predominant emphysema type), which may only be identified by a highly-trained chest radiologist [24].Importantly, CT images consist of embedded electron-density textural features, which can be exploited using texture analysis tools to map voxel intensity (attenuation) and spatial relationships [25] that are not easily identified by expert observers.In this regard, CT texture analysis and machine-learning tandems have been shown to predict COPD severity [26], progression [27], and showed stronger association with lung function compared to conventional densitometry measures [28,29].Imaging textures have also been shown to differentiate emphysema types [30,31] and provided improved radiological finding assessments as a second-reader [32].
Pulmonary functional magnetic resonance imaging (MRI) using hyperpolarized helium ( 3 He) and xenon ( 129 Xe) gases provides a way to measure pulmonary microstructure, ventilation, perfusion, and gas-exchange [33][34][35].In patients with COPD, such MRI measurements of ventilation and parenchyma microstructure [36,37] are predictive of acute exacerbations [38], airway narrowing and remodeling [39], and symptoms and severity [40,41]; these measurements also correlate with longitudinal changes in quality-of-life [42].Importantly, MRI measurements are sensitive to COPD disease-related changes in patients in whom CT and pulmonary function test results have not changed [43,44].In addition, MRI ventilation texture features have been shown to predict longitudinal lung-function decline in ex-smokers with and without COPD [45].However, to the best of our knowledge, ventilation and diffusion-weighted MRI measurements and radiomics-based CT/MRI textures have not been investigated for the prediction of 10-year all-cause mortality in ex-smokers with and without COPD.
Given all of this previous evidence, we hypothesized that, regardless of COPD status, incorporating chest MRI and CT measurements and image texture analysis in combination with machine-learning would provide unique prognostic information for mortality risk assessments in ex-smokers.Hence, here we evaluated MRI and CT measurements and employed imaging texture-analysis to predict all-cause mortality in ex-smokers with and without COPD after 10-years.

Study participants
All participants provided written informed consent to a study protocol approved by a local research ethics board (Institutional Ethics Board #00000984) in compliance with the Health Canada approved and registered protocol (clinicaltrials.govNCT02279329).All participants were recruited from a tertiary-care academic center and by advertisement in London, Canada between 2009 and 2012 as a convenience sample.These participants were followed for 10-years.Inclusion criteria were age of 50-85 years and a history of cigarette smoking >10 pack-years at baseline visit.Exclusion criteria included current smokers, claustrophobia and any contraindications for MRI or CT.Death dates were obtained from the health electronic record and survival time was calculated from the participant's baseline visit date.Some longitudinal results from this study were previously reported [46][47][48].In contrast to previous evaluations, this study quantified CT and MR imaging texture features at baseline and their association with 10-year all-cause mortality.

Pulmonary function tests and questionnaires
Spirometry, plethysmography and the diffusing capacity of the lungs for carbon monoxide (DL CO ) were measured according to the American Thoracic Society/European Respiratory Society standardization document [49] using a whole-body plethysmography system (MedGraphics Corporation, St Paul, MN, USA) and attached gas analyzer [46].COPD was defined as post-bronchodilator spirometry according to the Global Initiative for Chronic Obstructive Lung Disease (GOLD) criteria [10].Abnormal DL CO was defined as DL CO <75% pred as previously reported [50].The 6MWD [51] test and St. George's Respiratory Questionnaire (SGRQ) [52] were administered under the supervision of study personnel.

CT Acquisition and analysis
Thoracic CT was acquired using a 64-slice Lightspeed VCT scanner (GE Healthcare, Milwaukee, WI, USA) (64 × 0.625 mm, 120 kVp, 100 effective mA, tube rotation time= 500 ms, pitch= 1.25, reconstructed using a standard convolution kernel to 1.25 mm slice thickness, slices= 200-250), as previously described [53].Images were acquired in the supine position under breath-hold after inhalation of a 1 L bag of N 2 (from FRC lung volume) in order to match the lung volume for MRI.Total effective dose was estimated as 1.8 mSv using the ImPACT CT patient dosimetry calculator (Health Protection Agency [UK] NRBP-SR250).CT data were quantitatively evaluated using VIDAvision2.2software (VIDA Diagnostics Inc., Coralville, IA, USA) to quantify lung density using the relative area of lung less than −950 Hounsfield units (RA 950 ), total lung volume (TLV) and create a binary lung mask.CT vessel measurements were automatically generated using Chest Imaging Platform (Brigham and Women's Hospital, Boston MA) [54].

MRI Acquisition and analysis
Anatomic proton ( 1 H) and hyperpolarized 3 He MR images were acquired using a whole-body 3.0 Tesla Discovery MR750 system (GE Healthcare, Milwaukee, WI, USA), a whole-body radiofrequency coil and a fast gradient recalled echo (FGRE) sequence with a partial echo implementation, with acquisition parameters as previously described [55].Hyperpolarized 3 He MRI was acquired using a linear bird-cage transmit/receive chest coil (RAPID Biomedical GmbH, Wuerzburg, Germany).A commercial system (HeliSpin™, Polarean Inc, Durham, NC, USA) was used to polarize 3 He gas to 30-40% and doses (5 mL/kg body weight) diluted with N 2 were administered in 1 L Tedlar® bags (from FRC lung volume).Hyperpolarized 3 He MRI diffusion-weighted imaging was performed using a 2D multi-slice fast gradient-echo method, as previously described [55], during breath-hold for acquisition of two interleaved images with and without additional diffusion sensitization with b = 1.6 sec/cm 2 (maximum gradient amplitude [G] = 1.94 G/cm, rise and fall-time = 0.5 ms, gradient duration = 0.46 ms, diffusion time = 1.46 ms) [48].
MR images were evaluated for the measurement of ventilation defect percent (VDP) using semi-automated custom-built software, as previously described [36].Briefly, the anatomic 1 H and functional 3 He images were first co-registered to segment and remove the large airways (trachea), then a k-means clustering approach was used to generate ventilation clusters, with the lowest cluster representing ventilation defects as previously described [36].Ventilation abnormalities were quantified as the ventilation defect volume (VDV) and VDP was calculated as VDV normalized to the MRI-measured volume of the thoracic cavity [36].Diffusion-weighted images were automatically processed to generate apparent diffusion coefficient (ADC) values and images, after the removal of the large airways, as previously described [46].Abnormal 3 He ADC was defined as ADC >0.25 cm 2 /s, as previously reported [56,57].

Texture feature extraction, selection and machinelearning
CT images were first pre-processed by extracting the lungs and removing the large airways using the binary lung mask generated from the segmented CT images.Signal normalization was applied to the MR images, while a threshold between −1000 Hounsfield Units (HU) and 0 HU was applied to the CT images.The binary mask was applied to the segmented CT images in order to create regions of interest (ROI) for feature extraction.Similarly, the ROI from MR images was generated using the binary lung mask created by co-registering the 1 H and 3 He MRI acquisitions, as previously described [45].Next, 110 unique, unfiltered texture features were extracted in a voxel-by-voxel manner from CT and MR images using the open-source PyRadiomics platform (version 2.2.0) [58].For the best compromise between differentiation and resolution, a fixed bin number (FBN) discretization approach [59] was utilized to extract features between CT and MR modalities and ensure that textures are assessed against similar contrasts within the modality ROI, as previously described [60].Histogram and shape, first-, and higher-order texture features from run-length, gap-length, size-zone, neighborhood-dependence and co-occurrence matrices were computed.In addition, 376 wavelet-based texture features were extracted using four high-and low-pass filter combinations applied to the original image in x-and y-directions for wavelet decomposition [61].
To maximize the model generalizability and avoid overfitting, a combination of principal component analysis (PCA) and Boruta analysis [62] were implemented for feature selection.Feature selection step is primarily required for removing any redundant features and/or misleading data for improved modeling accuracy and also used for dimensionality reduction of data, enabling more efficient computations.Feature selection included the generation of nine components from PCA, which explained >94% of the variance in the data.Components were generated using principal component scores for each participant using a Varimax rotation method with Kaiser-normalization that converged after 38 iterations.
All variables in the models, including the texture features and emergent components generated for every participant, were subjected to Boruta analysis for ranking.The Boruta algorithm generated shadow features for comparisons and used a two-step correction for multiple testing, with an optimizable random forest classifier for iterations (number of trees in the forest = 150, maximum iterations = 200, maximum tree depth = 10 [branches], percentage of shadow feature threshold = 95%, alpha-level = 0.05).
Once all the features and parameters were selected, four machine-learning models were generated using: 1) clinical measurements, 2) imaging measurements, 3) image texture measurements, and 4) a combination of all available measurements.Five-fold cross-validation was performed to avoid overfitting or selection bias of the machine-learning models during the training step.During cross-validation, all participant data (n = 162) were randomly and evenly divided into five groups (n = 32/33) and for each of the "folds", one group was withheld for testing and remaining groups were used for training iteratively.Each fold utilized a different combination of testing and training groups in order to avoid data contamination and insure that no information was carried over from training to testing steps.Single (Naïve Bayes [63], Support Vector Machines [SVM] [64], Decision Trees [65], K-Nearest Neighbors [KNN] [65]) and ensemble (Bagged Trees [66], subspace Discriminant [67], subspace K-Nearest-Neighbors [67], and Random Under-Sampling Boosted Trees [RUSBoosted] [68]) machine-learning classifiers were implemented for predicting 10-year all-cause mortality in ex-smokers.Data were standardized and hyperparameter search was performed using MATLAB2021a (Classification Learner App) for every machine-learning model.Classification performance was evaluated using the mean of the respective 5-fold cross-validation area under the receiver-operator characteristic curve (AUC), as well as sensitivity and specificity calculated from model's confusion matrix.

Statistical analysis
Statistical analysis was performed using SPSS Statistics v28.0 (IBM Statistics, Armonk, New York, USA).Predictors of 10-year all-cause mortality were evaluated using binary logistic regression to generate odds ratios (OR).Shapiro-Wilk tests were used to determine the normality of the data.The pvalue significance was determined using the Mann-Whitney U-test for non-parametric data, and a post-hoc analysis using Holm-Bonferroni correction was applied for multiple comparison tests for the selected texture features.Statistical significance was considered using a 5% Type-I error threshold (p < 0.05).

Participant demographics and mortality
A CONSORT diagram provided in Figure 1 shows that 266 ex-smokers were enrolled and 99 were excluded from analysis due to enrollment in a sub-study using oscillatory positive expiratory pressure device (n = 33), due to cancelation or not completing all protocol tests (n = 61), and due to poor image quality (n = 5).In addition, five participants were excluded because 10-year follow-up was outside the Dec 2009-2022 window.As shown in Figure 2, the last participant enrolled on December 12 th 2012 and the mortality data window closed on December 13 th , 2022.
Demographic, clinical and imaging data for survivor and deceased ex-smoker subgroups are provided in Table 1 and summarized by COPD status in Table S1 (online supplement).In total, 162 ex-smokers were evaluated, including 93 ex-smokers with spirometry evidence of COPD and 69 ex-smokers without COPD.As shown in Table 1, 52/162 (32%) ex-smokers died within the 10-year window.Just over half of these (20/52, 53%) had CT evidence of emphysema, defined by the published threshold, RA 950 >6.8% [69].There were no significant differences between survivors and deceased participants for sex, pack-years, total lung capacity, and the CT measurements of emphysema (lowest attenuating cluster) and the pulmonary vasculature (total blood volume, volume of blood in small and mid-sized vessels).25) 20 (53) .03bmi: body mass index; fev 1 : forced expiratory volume in 1 s; % pred : percent of predicted value; fvC: forced vital capacity; rv: residual volume; tlC: total lung capacity; Dl CO : diffusing capacity of lung for carbon-monoxide; 6mWD: six minute walk distance; sGrQ: st.George's respiratory questionnaire; GOlD: Global initiative for Chronic Obstructive lung Disease; tbv: total blood volume; bv5: vessel volume for vessels less than 5 mm 2 ; bv5-10: vessel volume for vessels between 5-10 mm 2 ; bv10: vessel volume for vessels greater than 10 mm 2 ; laC: lowest attenuating cluster; ra 950 : relative area of lung less than -950 hounsfield units; aDC: apparent diffusion coefficient; vDP: ventilation defect percent.p=uncorrected values showing significant differences between survivor and deceased ex-smokers.

Texture Extraction and selection
Table 2 shows the six MRI and seven CT texture features which were the highest performing machine-learning predictors of 10-year all-cause mortality.Texture feature predictors were grouped into fine and coarse textures based on their mathematical definitions and are described in Table S2.As shown in Table S3, most of these texture features were significantly correlated with clinically-relevant measures (BMI, FEV 1 , 6MWD, and SGRQ).Figures 3 and 4 provide context for the CT and MRI textures respectively.Figure 3 shows coronal chest CT slices and in the inserts, representative ROI and example textures.In the survivor with COPD (top panel: S13), there was negligible CT evidence of emphysema (RA 950 =10.8%) and the CT texture was visually homogeneous (e.g.GLCM-Imc2 = 0.77) similar to a survivor without COPD (top panel: S11, GLCM-Imc2 = 0.76).In contrast, for the deceased ex-smoker with COPD and quantitative evidence of emphysema (RA 950 =24.9%) (bottom panel: S35, GLCM-Imc2 = 0.84) and deceased ex-smoker with no evidence of COPD (bottom panel: S50, GLCM-Imc2 = 0.82), CT textures appeared patchy and heterogeneously coarse.

Predicting 10-year all-cause mortality
As shown in Table 3, machine-learning prediction models for 10-year all-cause mortality were generated using: 1) clinical measurements, 2) imaging measurements, 3) image texture measurements, and 4) a combination of clinical and imaging measurements.The best performing clinical model (77% accuracy) was based on age, pack years, BMI, FEV 1 % pred , TLC, DL CO % pred , 6MWD and SGRQ score.The best performing imaging-based model (77% accuracy), included MRI VDP and ADC as well as CT HU15 th percentile, TBV, BV 5 /TBV, BV 10 , BV 10 /TBV, and %RA 950 .The predictive model based exclusively on imaging texture features outperformed both clinical and imaging-based models (80% accuracy).Finally, the combined clinical-imaging model had the overall best performance (83% accuracy) which included: DL CO , MRI ADC, as well as MRI and CT texture features.
Table 4 summarizes the best performing measurements in descending order, across the key statistical tests and models generated.The common top performing measurements were DL CO , MRI-ADC, MRI-Wavelet-LL-FO-Skewness and CT-GLCM-Imc2.

Discussion
A recent investigation in patients with COPD [70] showed that unsupervised learning of chest CT measurements improved predictions of progression, exacerbation, and mortality risk.Another recent study [71] showed that machine learning models of clinical and CT imaging measurements outperformed the current standard (BODE and ADO indices) for predicting all-cause mortality.
We applied machine-learning algorithms to generate models for predicting all-cause mortality using clinical, CT, MRI and imaging textures.We observed that a combined model consisting of DL CO , MRI ADC and image textures, outperformed all other models.Surprisingly, none of the individual components of the BODE index were included in the "best performing" model after the feature selection step.Ensemble machine-learning models outperformed single machine-learning models, suggestive of complex, non-linear relationships between the individual imaging textures and mortality.
terminal airspace enlargement due to air-trapping, emphysema or both [43,46].This result also agrees with the COPDGene study results, whereby emphysema progression over 5-years was associated with mortality in ever-smokers with trace emphysema [15].This result also agrees with the finding that the MRI ADC measurement is highly sensitive to terminal airspace enlargement as previously shown [43] and by comparison with histological measurements of the mean linear intercept [72].
We observed that one of the "fine" CT textures (GLCM-Imc2) had the strongest independent association with 10-year all-cause mortality, which after adjusting for confounders resulted in a three-fold increased mortality risk (OR = 3.55, p = 0.001).Such "fine" CT textures may be intuitively considered as reflecting tissue attenuation heterogeneity.To provide context, previous work demonstrated that the spatial arrangement of low attenuating voxels, or the size and arrangement of emphysematous clusters can differentiate patients with similar COPD disease severity [73].However, unlike a more commonly used CT measurement of emphysema, RA 950 , CT texture features reflect complex spatial heterogeneity, which may be argued is more sensitive to emphysema [31,74] and perhaps similar to what is often visually interpreted as emphysema by expert chest CT radiologists [24].In contrast throughout this study, all MRI texture features selected were "coarse" and could be considered as reflecting the compactness and asymmetry of the ventilation distribution in the lungs.
Underscoring the power of MRI and CT measures as predictors of all-cause mortality, forest plots revealed significant associations with mortality for all selected MRI/CT textures, CT vascular measurements and MRI ADC.These results are consistent with those of recent studies reporting that CT emphysema [15,70,71] and vascular [17,54,75] measurements are associated with mortality and disease progression in COPD.Our findings contribute to the growing body of evidence about the utility of imaging for the management of patients at risk for, and with, a diagnosis of COPD.
In agreement with the ROC analysis and forest plots, Kaplan-Meier curves revealed that ex-smokers with abnormal DL CO , MRI ADC and a specific MRI texture were at greatest risk of mortality.Of note, across all models, the measurements which remained significant and highly associated with 10-year mortality were DL CO , MRI ADC, MRI Wavelet-LL-FO-Skewness and CT GLCM-Imc2.These consistent observations are indicative of the added-value of imaging measurements and textures for mortality risk assessments in ex-smokers with and without COPD.
We acknowledge a number of study limitations.First, direct comparison of our results with the BODE index should be undertaken with caution because we did not  3 he aDC, hu15 th percentile, tbv, bv 5 /tbv, bv 10 , bv 10 /tbv, and % ra 950 .§features selected for training included: mri wavelet-low-low-firs t-order-skewness, mri wavelet-low-low-kurtosis, Ct gray level co-occurrence matrix (GlCm)-informational measure of correlation (imc2), Ct wavelet-high-high-gray-level dependence matrix (GlDm)-dependence variance, Ct wavelet-high-high-run variance, Ct wavelet-low-low-firs t-order-skewness, Ct wavelet-low-low-GlCm-imc1. □Combined model included: Dl CO , aDC, mri Wavelet-ll-fO-skewness, Ct ll-GlCm-imc1, Ct shape-sphericity, Ct Wavelet-hh-GlDm-Dependence variance and mri Wavelet-ll-fO-Kurtosis. acquire the modified Medical Research Council dyspnea-score component of the BODE index.Second, our study included a relatively small sample size compared to other CT studies [76,77].Fortunately, we employed methods to avoid overfitting so the machine-learning methods we used may be generalized with larger datasets and multicenter data in the future.We also acknowledge that lung imaging measurements are influenced by lung volume [78], and in this study, data were captured at FRC + 1L which for many ex-smokers is within 90% of TLC [78,79].Finally, the availability of functional MRI for chest imaging is still limited to research sites.Thus, while such quantitative MR measurements are very helpful, CT measurements are certainly more readily available and generalizable to most clinical centers.Inclusion of at risk ex-smokers may help identify individuals who are at high-risk of death and can potentially be included in clinical trials [80].It is well established that asymptomatic ex-smokers with mildly abnormal DL CO are at-risk of developing COPD within four years [50].The abnormal MRI ADC values reported here have been shown to reflect mild air-trapping and/or subclinical emphysema [43].Unfortunately, we did not make small airways disease measurements, so we are unable to comment on whether small airway abnormalities have already initiated in these participants with normal spirometry [18][19][20][21].What we do know is that 14/69 (20%) of ex-smokers and 38/93 (41%) of COPD participants were deceased after 10 years and that similar MRI and CT measurements and textures indicative of emphysema helped explain risk.The results presented here expand on previous CT findings in COPD [15,17,54,70,71,75] and extend the application of such imaging measurements to ex-smokers without, but at risk of COPD.

Conclusion
In ex-smokers, regardless of COPD status, DL CO , CT and MR imaging measurements and textures resulted in high accuracy models for predicting mortality risk.Texture measurements provide a way to reveal MRI and CT lung Figure 6.logistic regression models for associations between all-cause mortality and clinical, imaging and textural measurements.all-cause mortality assessment was conducted in 162 ex-smokers, of whom 52 deceased across the longitudinal study duration (10-years).bolded values indicate categories where 95% Ci did not include 1.0 (P < 0.05).*all odds ratios were adjusted for age, bmi, sex, and pack-years.GlCm = gray level co-occurrence matrix; 6mWD = six minute walk distance; aDC = apparent diffusion coefficients; Dl CO =diffusing capacity of the lung for carbon monoxide; GlDm = gray level dependence matrix; hh = high-high pass filter; ll = low-low pass filter; svr = surface volume ratio; Dv = dependence variance; rv = run variance; fev 1 =forced expiratory volume in 1 second; fvC = forced vital capacity; sGrQ = st.George's respiratory questionnaire; laC = lowest attenuating cluster; ra 950 =relative area of lung less than -950 hounsfield units; all texture feature abbreviations and descriptions can be found in supplementary table s2.
pathologies that are not visible to the human eye and may help predict 10-year all-cause mortality in ex-smokers.

Figure 1 .
Figure 1.COnsOrt flow diagram.Of the 266 participants enrolled in the tinCan study, 33 were enrolled in a sub-study, 61 either canceled or did not complete all required tests during visit 1, and five had Ct or mri artifacts which precluded analysis.Of the 167 participants who completed visit 1, five were not yet within their 10-year follow-up timeframe.at follow-up, there were 52 deceased participants, of whom 14 were ex-smokers and 38 were ex-smokers with COPD.

Figure 2 .
Figure 2. Participant enrollment and follow-up timeframe.arrows showing the timeline of participant enrollment and follow-up period for mortality data collection within 10-years from the initial visit in all ex-smoker participants.Of the 162 ex-smokers analyzed, the first participant completed the baseline visit on December 1 st , 2009 and the last participant completed their baseline visit on December 12 th , 2012.

Figure 5 .
Figure 5. receiver-operator characteristic curves of texture features and clinical variables.top Panel: logistic regression analysis of individual clinical variables at predicting 10-year all-cause mortality in ex-smoker participants.Dl CO had the best auC=.736.middle Panel: logistic regression analysis of standard imaging measurements at predicting 10-year all-cause mortality in ex-smoker participants. 3he aDC had the best auC=.738.bottom Panel: logistic regression analysis of imaging texture features at predicting 10-year all-cause mortality in ex-smoker participants.Ct Wavelet-ll-GlCm-imc1 had best auC=.787.

Table 1 .
Participant demographics, pulmonary function and imaging measurements in survivors and deceased ex-smokers.

Table 2 .
mri and Ct texture features in survivors and deceased ex-smokers.