CT whole lung radiomic nomogram: a potential biomarker for lung function evaluation and identification of COPD

CT whole lung radiomic nomogram: a potential biomarker for lung function evaluation and identification of COPD

PDF

Tao-Hu Zhou¹^,²^,^†, Xiu-Xiu Zhou¹^,^†, Jiong Ni³, Yan-Qing Ma⁴, Fang-Yi Xu⁵, Bing Fan⁶, Yu Guan¹, Xin-Ang Jiang¹, Xiao-Qing Lin¹^,⁷, Jie Li¹^,⁷, Yi Xia¹, Xiang Wang¹, Yun Wang¹, Wen-Jun Huang¹^,⁸, Wen-Ting Tu¹, Peng Dong², Zhao-Bin Li⁹, Shi-Yuan Liu¹, Li Fan¹^,^*

Military Medical Research | 2025, 12(1) : 36 - 47

Less

Military Medical Research | 2025, 12(1): 36-47

• RESEARCH •

CT whole lung radiomic nomogram: a potential biomarker for lung function evaluation and identification of COPD

Full

Affiliations

¹Department of Radiology, the Second Affiliated Hospital of Naval Medical University, Shanghai 200003, China

²School of Medical Imaging, Shandong Second Medical University, Weifang 261053, Shandong, China

³Department of Radiology, School of Medicine, Tongji Hospital, Tongji University, Shanghai 200065, China

⁴Department of Radiology, Zhejiang Province People’s Hospital, Affiliated People’s Hospital of Hangzhou Medical College, Hangzhou 310014, China

⁵Department of Radiology, Sir Run Run Shaw Hospital, Zhejiang 310018, China

⁶Jiangxi Provincial People’s Hospital, the First Affiliated Hospital of Nanchang Medical College, Nanchang 330006, China

⁷College of Health Sciences and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China

⁸Department of Radiology, the Second People’s Hospital of Deyang, Deyang 618000, Sichuan, China

⁹Department of Radiation Oncology, Shanghai Jiao Tong University Affiliated Sixth People’s Hospital, Shanghai 200233, China

Published: 2025-01-10 doi: 10.1186/s40779-024-00516-9

Outline

Abstract

Less

Background: Computed tomography (CT) plays a great role in characterizing and quantifying changes in lung structure and function of chronic obstructive pulmonary disease (COPD). This study aimed to explore the performance of CT-based whole lung radiomic in discriminating COPD patients and non-COPD patients. Methods: This retrospective study was performed on 2785 patients who underwent pulmonary function examination in 5 hospitals and were divided into non-COPD group and COPD group. The radiomic features of the whole lung volume were extracted. Least absolute shrinkage and selection operator (LASSO) logistic regression was applied for feature selection and radiomic signature construction. A radiomic nomogram was established by combining the radiomic score and clinical factors. Receiver operating characteristic (ROC) curve analysis and decision curve analysis (DCA) were used to evaluate the predictive performance of the radiomic nomogram in the training, internal validation, and independent external validation cohorts. Results: Eighteen radiomic features were collected from the whole lung volume to construct a radiomic model. The area under the curve (AUC) of the radiomic model in the training, internal, and independent external validation cohorts were 0.888 [95% confidence interval (CI) 0.869–0.906], 0.874 (95%CI 0.844–0.904), and 0.846 (95%CI 0.822–0.870), respectively. All were higher than the clinical model (AUC were 0.732, 0.714, and 0.777, respectively, P<0.001). DCA demonstrated that the nomogram constructed by combining radiomic score, age, sex, height, and smoking status was superior to the clinical factor model. Conclusions: The intuitive nomogram constructed by CT-based whole-lung radiomic has shown good performance and high accuracy in identifying COPD in this multicenter study.

Key words

Chronic obstructive pulmonary disease (COPD) / Computed tomography (CT) / Radiomic

Cite this Article

Tao-Hu Zhou, Xiu-Xiu Zhou, Jiong Ni, Yan-Qing Ma, Fang-Yi Xu, Bing Fan, Yu Guan, Xin-Ang Jiang, Xiao-Qing Lin, Jie Li, Yi Xia, Xiang Wang, Yun Wang, Wen-Jun Huang, Wen-Ting Tu, Peng Dong, Zhao-Bin Li, Shi-Yuan Liu, Li Fan. CT whole lung radiomic nomogram: a potential biomarker for lung function evaluation and identification of COPD[J]. Military Medical Research, 2025 , 12 (1) : 36 -47 . DOI: 10.1186/s40779-024-00516-9

Full Text

Less

Background

Less

Chronic obstructive pulmonary disease (COPD) is a chronic inflammatory disorder with high heterogeneity and characterized by continuous airflow limitations. The current gold standard for diagnosing and evaluating COPD is the pulmonary function test (PFT)[1], which yields the ratio of forced expiratory volume in 1 s to forced vital capacity (FEV1/FVC) and the percentage of predicted FEV1 (FEV1% predicted). In China, the incidence of COPD in people ≥40 years old is 13.7%, however, the awareness rate of COPD is very low, less than 1%[2]. Many people have been underdiagnosed because the PFT is not widely used for screening in China. Based on the survey of PFT performance in people aged 40 years and above in China, the PFT rate in Chinese residents aged ≥40 years was 6.7% (95%CI 5.2%–8.2%) in 2019–2020, the overall PFT rate was still at a low level[3]. In contrast, the popularity of chest computed tomography (CT) is high, especially with large-scale lung cancer screening. Meanwhile, the 2023 Global Initiative for Chronic Obstructive Lung Disease report emphasized the importance of CT in evaluating patients with stable COPD, highlighting the role of imaging[4]. This evidence-based suggestion was made due to the limitations of PFT, which, despite being the gold standard, cannot be used for focal evaluation and does not show the lungs. Additionally, due to the high heterogeneity of COPD, focal and visual evaluations are considered to play an important role in guiding clinical decisions. As the most common and powerful imaging technique, CT has great potential in COPD with the rapid development of CT-based artificial intelligence (AI).

Furthermore, CT has a high anatomic resolution and can be used to evaluate the changes in lung parenchyma, small airway, and pulmonary blood vessels that occur with lung function decline and aging. However, it is not commonly used to simultaneously assess COPD abnormalities. In addition, the subjective evaluation of the lung parenchyma and small airway lesions was influenced by the experience of the radiologist. Especially as the number of patients with COPD continues to increase, visual assessment of lung lesions by radiologists is becoming more expensive and laborious. Therefore, identifying a method for quantitatively evaluating the whole lung is critical to obtaining a comprehensive evaluation of COPD.

Radiomic is a relatively novel approach that can rapidly collect quantitative high-throughput features from medical images (e.g., CT), such as complex patterns that are not easily recognized or quantified by the naked eye[5], showing great potential in clinical decision-making. At present, radiomic research on COPD is rapidly expanding, and several studies have indicated that such a method may have particular advantages in patients with COPD[6,7]. The radiomic features are extracted from segmented lesions. However, radiologists have found that manually segmenting diffuse and heterogeneous lung lesions, such as emphysema, interstitial lung disease, and coronavirus disease 2019[8,9], are difficult and timeconsuming due to unclear boundaries or low contrast on CT imaging. Additionally, different radiologists may segment the lesions differently when evaluating diffuse lung diseases[10]. Therefore, it is very important to develop a method to automatically segment diffuse lesions. It has been reported that the application of an AI-based system to detect diseases can reduce the workload of radiologists and maintain the accuracy of diagnoses[11]. Because COPD is a diffuse chronic lung disease, automatic segmentation of the whole lung region would help to comprehensively quantify lung abnormalities and aid clinical treatment decision-making. A recent editorial suggests that automated detection of COPD based on chest CT findings using radiomic or deep learning techniques has great potential to reduce the current underdiagnosis of COPD, particularly in high-risk cohorts[12]. The purpose of this study was to explore the performance of CT-based automatic segmentation of whole-lung radiomic in differentiating COPD from non-COPD and to assess the value of CT-based radiomic in lung function evaluation.

Materials and methods

Less

Patients

A total of 2941 patients who were admitted and underwent PFT at 5 centers, including the Second Affiliated Hospital of Naval Medical University, Tongji Hospital, School of Medicine, Tongji University, Zhejiang Province People’s Hospital, Sir Run Run Shaw Hospital, and the First Affiliated Hospital of Nanchang Medical College, were retrospectively recruited from February 2013 to December 2022 in the CSD-COPD cohort. The inclusion criteria were as follows: 1) chest CT and PFT both performed in the same hospital; 2) less than two weeks between PFT and chest CT; and 3) complete thin-slice (<2 mm) chest CT images. The exclusion criteria were as follows: 1) other comorbid thoracic diseases (e.g., pneumonia, pulmonary atelectasis, lung nodules larger than 6 mm or masses, asthma, and pleural effusion); 2) malignant tumors; and 3) spine implants or substantial image artifacts. Finally, 2785 patients were included in this study. Among them, 1714 patients from the Tongji Hospital, School of Medicine, Tongji University, Zhejiang Province People’s Hospital, Sir Run Run Shaw Hospital, and the First Affiliated Hospital of Nanchang Medical College were randomly assigned to the training cohort (n=1200) and the internal validation cohort (n=514) in a ratio of 7:3. Patients from the Second Affiliated Hospital of Naval Medical University were assigned to an independent external validation cohort (n=1071). Figure 1 shows the workflow for patient inclusion and exclusion. The basic clinical information of the patients, including age, sex, weight, height, body mass index, and smoking status, was collected through the electronic medical records system.

Pulmonary function parameters (FEV1, FVC) were measured with PFT apparatus (CHEST Multifunction Spirometer HI-801, Japan; Ganshorn Medizin Electronic GmbH; Carefusion GmbH, Hoechberg, Germany; Masterscreen PFT Pro, Carefusion, Netherlands), as well as CT acquisition parameters in Additional file 1: Table S1. The diagnostic criteria for PFT in COPD are as follows: FEV1/FVC <0.7 with an increase of FEV1 <200 ml after the use of a bronchodilator. In contrast, this study included patients with FEV1/FVC ≥0.7 and the FEV1% predicted ≥80% after bronchodilation as the non-COPD group. Participants in the training, internal validation, and independent external validation cohorts were divided into COPD and non-COPD groups according to these criteria.

This study was approved by the institutional review boards at 5 centers, and informed consent was waived due to the retrospective nature of this study (ChiCTR2300069929).

Whole-lung CT image segmentation and CT image preprocessing

Using a deep-learning model of open access U-net (R231) (https://github.com/JoHof/lungmask) for the automatic segmentations, which has been trained using different large-scale datasets covering a wide range of visual variability, the reliability of this method has been proved[13]. First, the right and left lungs were automatically segmented. Then, we merged the right and left lungs into a combined region of interest (ROI) (Fig. 2).

Since manual segmentation is often regarded as the ground truth, we assessed the consistency between manual and fully automatic segmentation in 20 randomly selected patients across the cohort. CT images of 20 patients were exported to ITK-SNAP software (version 3.8.0, www.itksnap.org) for manual segmentation. The consistency between manual and fully automatic segmentation was assessed using the Dice index, an objective measure that quantifies the spatial overlap between two contours. The remaining cases were then automatically segmented.

Before extracting the radiomic features, the images were preprocessed, which consisted of three steps. First, we used linear interpolation to resample the images to 1 mm×1 mm×1 mm. Second, we used gray-level discretization to convert continuous images into discrete integer values. Finally, log and wavelet image filters were used to eliminate the mixed noise in the process of image digitization and obtain low-frequency or high-frequency features.

Radiomic feature extraction and selection

A total of 1218 lung radiomic features were extracted from each volume of interest using the open-source package PyRadiomics (version 3.0.1, https://pyradiomics.readthedocs.io/en/latest/), including first-order, gray level co-occurrence matrix, gray level run length matrix, gray level size zone matrix, gray level dependence matrix, and shape features. The radiomic features extracted by this software are in accordance with the image biomarker standardization initiative. The Z score method was used to normalize the features and eliminate the difference in numerical scale.

The following three steps were used to select the best radiomic features. First, redundant features whose correlation coefficient with other features is greater than 0.90 were removed. Second, the maximal redundancy-minimal relevance algorithm was used to eliminate the redundant and irrelevant features. Minimal redundancy maximal relevance has been proven to be an effective and reliable feature selection method for radiomic, which can consider both the importance of features and the correlation between features to find the optimal feature subset[14,15]. Finally, the least absolute shrinkage and selection operator (LASSO) regression algorithm and penalty parameter adjustment were used for tenfold cross-validation. The optimal feature dataset with the smallest cross-validation binomial deviation was selected, and the non-zero coefficients were defined as the weight of the selected feature, representing the correlation between the feature and COPD. LASSO is a widely used embedded method for radiomic feature selection in high-dimensional data[16]. Finally, the Radscore of each patient was calculated by a linear combination of the selected feature and coefficient vectors, and the radiomic model was constructed.

Model construction, radiomic nomogram, and performance evaluation

Three models were constructed, including the clinical model, radiomic model, and combined model. Univariate logistic regression analysis was used to obtain statistically significant risk variables, and then multivariate analysis was performed to establish clinical and combined models. A radiomic nomogram was generated to visualize the combined model, graphically evaluate variable importance, and calculate prediction accuracy. The DeLong test was used to compare the area under the curves (AUCs) of the clinical model, radiomic model, and combined model. The calibration curves (Hosmer-Lemeshow test) were performed to evaluate the calibration of the nomogram. Decision curve analysis (DCA) was applied to evaluate the clinical practicability of the nomogram.

Statistical analysis

IBM SPSS Statistics (version 26.0; IBM Corp., New York, USA) and R software (version 4.2.2; http://www.Rproject.org) were employed for statistical analysis. Measurement variables are expressed as the mean±standard deviation (SD). Normally distributed continuous variables were compared using the Student’s unpaired t-test and non-normally distributed data were compared using the Mann-Whitney U test. Categorical variables were compared by the χ² test between groups. Independent predictors were identified from the clinical variables by multivariate logistic regression. P<0.05 indicated statistical significance. LASSO regression was conducted using the “glmnet” package. Additionally, the “rms” package was employed for drawing calibration plots and conducting multivariate logistic regression. The package of receiver operating characteristic (ROC) was utilized for drawing the ROC curves of the radiomic signatures, while the “rmda” package was utilized for DCA.

Results

Less

Clinical characteristics

In total, 2785 patients (male 1715, female 1070; non-COPD group 1377, COPD group 1408) with an average age of (65.4±11.4) years old were included. Table 1 displays the basic demographics of all patients studied. In the training and internal validation cohorts, the distribution of the patients in the 4 independent centers was as follows: 1529 patients from Tongji Hospital, School of Medicine, Tongji University, 73 patients from Zhejiang Province People’s Hospital, 42 patients from Sir Run Run Shaw Hospital, and 70 patients from the First Affiliated Hospital of Nanchang Medical College. The training cohort included 491 non-COPD patients and 709 COPD patients, and the internal validation cohort included 218 non-COPD patients and 296 COPD patients. A total of 1071 patients from the Second Affiliated Hospital of Naval Medical University were assigned to the independent external validation cohort, consisting of 668 patients without COPD and 403 patients with COPD. Significant differences were observed in age, sex, height, weight, body mass index, and smoking status between the non-COPD and COPD groups (P<0.05) in the training, internal and independent external cohorts. However, in the internal validation cohort, the difference between current and former smokers was not significant (Table 1).

Consistency assessment between manual and fully automatic segmentation

The segmentations were assessed using the Dice index, an objective measure that quantifies the spatial overlap between two contours. The mean Dice coefficient between manual and automatic segmentation was (0.97±0.06) (Additional file 1: Fig. S1).

Feature screening and radiomic signatures establishment

A total of 1218 radiomic features were normalized with the Z score method. After Pearson’s correlation analysis, 935 radiomic features (absolute value of Pearson correlation coefficients >0.9) were eliminated. Therefore, a total of 283 features were retained. Finally, 18 radiomic features with non-zero coefficients were selected by LASSO regression (Fig. 3). By linearly combining those features after weighting by their corresponding coefficients, we constructed the radiomic signature. The Radscore calculation formula is provided in the Additional file 1.

Performance comparison of radiomic model, clinical model, and combined model

The boxed scatter plots for the Radscore are shown in the Additional file 1: Fig. S2. As revealed by the Wilcoxon test, the Radscore exhibited significant differences between the COPD group and the non-COPD group (P<0.001). In addition, according to univariate and multivariate regression analysis, the Radscore was independently associated with COPD. Moreover, age, sex, height, and smoking status were identified as independent predictors of COPD by multivariate regression and included in the construction of the clinical model (Table 2). Last, the Radscore was integrated with these independent predictive factors to construct the combined model. The combined model calculation formula is described in Additional file 1.

Figure 4 and Table 3 showed the performances of the radiomic, clinical, and combined models. The constructed radiomic model contains 18 screened features and had a good degree of differentiation, with AUCs of 0.888 (95%CI 0.869–0.906), 0.874 (95%CI 0.844–0.904), and 0.846 (95%CI 0.822–0.870) in the training, internal and external validation cohorts, respectively. According to the DeLong test, there was a significant difference in the AUCs between the combined model and the clinical model (P<0.001 in the three cohorts). The DeLong test also showed that the AUC of the combined model and the radiomic model in the training cohort was significantly different [AUC=0.893 (95%CI 0.875–0.911) vs. AUC=0.888 (95%CI 0.869–0.906); P=0.02] and in the external validation cohort [AUC=0.853 (95%CI 0.830–0.877) vs. AUC=0.846 (95%CI 0.822–0.870); P=0.04], but there was no significant difference in the internal validation cohort [AUC=0.873 (95%CI 0.843–0.903) vs. AUC=0.874 (95%CI 0.844–0.904); P=0.71].

Development and performance of the nomogram

The visualization of the nomogram and the combination of radiomic and common clinical factors are helpful for doctors to conduct health education consultations for patients. The combined model was converted into a nomogram, and the total score obtained from the nomogram was used to predict the risk of COPD (Fig. 5a). The Hosmer-Lemeshow test showed that the calibration curves of the combined model for predicting COPD in the training, internal and external validation cohorts matched the actual data very well (P=0.972, 0.149 and 0.06, respectively) (Fig. 5b). According to DCA (Fig. 5c), the combined model showed a greater benefit than the clinical model in predicting COPD risk in the training cohort when the probability threshold in the clinical decision of the patient or physician was greater than 0.1. The nomogram showed the highest clinical net benefit across all threshold probability ranges in the training cohort, suggesting that the nomogram is a reliable tool for clinically predicting COPD. An example of the nomogram in use is shown in Fig. 6. Similar to the points scoring system, we assigned points for each predictor of COPD and then equated these predictors with the risk of COPD. We can read the top score scale upward from the predictors to determine the points score associated with patient age, height, smoking status, sex, and the Radscore. Once a score has been assigned to each predictor, an overall score is calculated. Then, the total score is converted to the probability of COPD by reading the associated probability of COPD from the total point scale.

Discussion

Less

COPD is a heterogeneous disease that causes a series of abnormalities, including small airway remodeling, lung vessel remodeling, and the formation of emphysema. The ability to comprehensively evaluate the disease is very important. Although PFT is the current clinical gold standard, CT plays an important role in the management of COPD due to its advantages of focal, accurate, and visual evaluations. In our large multicenter cohort of participants, we innovatively proposed the construction of a CT-based whole lung radiomic nomogram to identify COPD. The AUCs of the model were 0.893, 0.873, and 0.853 in the training, the internal validation, and the independent external validation cohorts, respectively. The subsequently constructed nomogram is intuitive, which can improve the value of CT in evaluating lung function and help to detect more underdiagnosis of COPD in clinical routine work.

The incidence and disease burden of COPD is high in China, the overall pulmonary function detection rate is still at a low level, and many people have been underdiagnosed. In contrast, the popularity of chest CT is very high, especially with the large-scale chest CT screening for lung cancer. Moreover, more and more community health service centers will be equipped with CT. Therefore, the most important clinical scenario is for the large-scale lung cancer screening population that usually does not perform PFT, and many underdiagnosed COPD can be found through our model prediction, which can help enhance the detection and early intervention of COPD, reduce the socioeconomic burden and improve the patient’s life quality. A recent study revealed that CT-based radiomic features extracted only from inspiratory CT scans outperformed existing advanced methods in detecting COPD on both standard- and low-dose CT scans. The model was constructed with the standard-dose CT radiomic feature[17].

Radiomic has great potential in obtaining useful medical information and enhancing the accuracy of clinical differential diagnosis. A previous study has identified the value of lung radiomic features based on CT imaging and clinical manifestations in the assessment of COPD[7]. CT finding of COPD patterns might be obscure and diffuse, making it difficult to accurately delineate abnormal areas. Li et al.[18] randomly selected 42 non-overlapping ROIs from 11 axial CT sections of every patient to extract radiomic features, with an AUC of 0.97. However, the approach they used could not comprehensively evaluate the disease in the whole lung. In contrast, the automatic segmentation of the whole lung into the whole ROI allows a comprehensive evaluation of the lung, with an AUC of 0.893 in this study. Automatic segmentation can improve efficiency and reduce inter- and intra-observer differences. However, the AUC of COPD identified in our study was lower than that of Li et al.[18], which may be related to the fact that they applied machine learning technologies for further radiomic feature selection based on LASSO regression. Nam et al.[19] trained and validated a deep learning method to predict the prognosis of COPD patients based on chest radiography, with an AUC of 0.76. Notably, a significant proportion of patients with COPD had normal chest X-rays based on subjective evaluations. Compared to chest X-ray images, chest CT scans are more sensitive to changes in COPD. Therefore, chest CT was used to evaluate COPD in this study. The greatest characteristic of radiomic in this study was that the whole lung was combined into one ROI to extract the radiomic features. Because COPD is a diffuse and heterogeneous disease involving the pulmonary parenchyma, small airway, and lung blood vessels, focal ROIs cannot fully represent the pathological changes induced by the disease. Moreover, we used the deep convolutional neural network extension based on the U-net architecture for lung segmentation[13,20].

According to our results, 283 potential radiomic features were selected on CT images, of which the LASSO regression model ultimately identified 18 predictors for constructing radiomic signatures. The radiomic features we screened were divided into four types (first-order, morphologic features, texture features, and wavelet features), which were significantly different between the non-COPD and COPD groups. These features essentially reflect information from the distribution of pixel intensity and texture morphology that radiologists cannot detect manually[21]. Morphologic features describe the size, volume, and shape of the volume of interest, while first-order features mainly reflect the internal texture of the lesions. Textural features, including the gray level co-occurrence matrix and gray level dependence matrix, describe the spatial relationship between each pixel and its neighbors. Wavelet features mainly reflect the time-frequency domain within the lesion[22]. Among the selected radiomic features, Wavelet LLL gldm LowGrayLevelEmphasis and wavelet LHL glcm ClusterShade have the highest significance and robustness in identifying COPD. They represent the intensity and textural features of lesions in high-intensity CT voxels. To some extent, radiomic is a quantitative method. Conventional quantitative CT evaluation has been applied in COPD diagnosis, severity evaluation, prognosis, and many other aspects. Cho et al.[23] reported the performance of an integrated model of quantitative features, such as emphysema, airway remodeling, pulmonary vascular diseases, and air trapping, extracted by fully automated inhouse software (AVIEW) with a radiomic approach as a predictor of survival in COPD patients, and they found that their integrated model outperformed a model constructed using only a single quantitative parameter. In our research, we established a whole-lung radiomic signature describing airways, blood vessels, and emphysema, similar to Cho et al.[23]. We found the radiomic model outperformed the clinical model, which is similar to a recent study by Amudala Puchakayala et al.[17].

In our model, age, sex, height, and smoking status were selected as independent risk factors for identifying COPD. Smoking is one of the most common clinical risk factors, but a significant proportion of patients with COPD have never smoked[24]. Therefore, the model used smoking status as a surrogate for total smoking exposure in our study. Additionally, as in previous studies, age, sex, and height were independent predictors of COPD[2,25,26].

The AUC of the combined model with clinical and radiomic features was 0.893, 0.873, and 0.853 respectively, which was superior to the clinical model in the three cohorts, and slightly better than the radiomic model in the training and external validation cohorts. These findings are similar to the recent study published in Radiology[17], which indicated radiomic alone is a potent tool for identifying COPD. However, to provide a predicting tool for the probability of COPD occurrence at the individual level, we still constructed a nomogram based on the combined model. Compared with traditional methods, the nomogram can predict more quickly, conveniently, and accurately, especially it can help the physician to consult patients for health education.

Additionally, the diagnostic accuracy of the different models was assessed using an external validation cohort. The combined model showed slightly better than the radiomic model in the external validation cohort (P=0.04), demonstrating that the model has good accuracy for this population. The results indicated the model has good predictive performance for new, unfitted data, as well as its high prediction capability and robustness. These findings are consistent with others utilizing radiomic to predict COPD survival[23,27], spirometry-based evaluation of emphysema and severity[28], COPD exacerbations[29], COPD stage classification[30,31], and analysis of COPD and resting heart rate[32].

Some limitations should be noted in this study. First, selection bias was inevitable due to its retrospective nature. The number of patients in the 5 centers was imbalanced, but the performance of the nomogram was good, confirming the universal applicability of our model. Second, only the CT radiomic features were evaluated, not common CT quantitative and qualitative parameters that are valuable in evaluating COPD. In future studies, we will extract common quantitative and qualitative parameters from paired inspiratory and expiratory CT images into our prediction model. Third, regarding the clinical variables, we only considered those that are most common and easy to acquire from CT scans. To more objectively evaluate the performance of the clinical model, future studies should incorporate more clinical variables, including symptoms. Fourth, the non-COPD patients in this study were defined as FEV1/FVC ≥0.7 and an FEV1% predicted ≥80% after bronchodilation, the preserved ratios of impaired spirometry (PRISm, FEV1/FVC ≥0.7 and an FEV1% predicted <80%) were not included. PRISm is considered to be a high-risk factor for COPD, so it is very important to distinguish it from COPD. Our team will perform a tri-classification study.

Conclusions

Less

In conclusion, whole-lung CT radiomic can be used as a good biomarker for the identification of COPD, not only for the lung cancer screening population but also for all the patients who performed chest CT examinations. With the gradual development of AI technology, the quantitative and intuitive nomogram based on whole-lung CT radiomic may have wide clinical application value and more research should be performed toward automatic detection of COPD.

Abbreviations

Less

AI: Artificial intelligence; AUC: Area under the curve; COPD: Chronic obstructive pulmonary disease; CT: Computed tomography; DCA: Decision curve analysis; FEV1/FVC: Ratio of forced expiratory volume in 1 s to forced vital capacity; FEV1%: predicted Percentage of predicted FEV1; LASSO: Least absolute shrinkage and selection operator; PFT: Pulmonary function test; ROI: Region of interest

Supplementary information

Less

The online version contains supplementary material available at https://doi.org/10.1186/s40779-024-00516-9.

Additional file 1: Table S1 CT protocols of the five centers. Fig. S1 Typical lung region segmentation results from the original chest HRCT images segmented fully automatically and manually in the transverse, coronal, and sagittal planes. Fig. S2 Boxplots show the whole lung CT radiomic signatures in COPD group were much higher than the non-COPD group in both the training (left) and test cohort (right). The calculation formula for the Radscore. The calculation formula for the combined model.

Acknowledgements

Less

Not applicable.

Authors’ contributions

Less

LF conceived, designed, and supervised the study. THZ and XXZ performed data analysis and drafted the manuscript. THZ conducted a major clinical experiment. JN, YQM, FYX, BF, YG, XAJ, XQL, JL, YX, XW, YW, WJH, and WTT contributed to the imaging and clinical data collection. PD, ZBL, and SYL supervised the literature review and data quality control. LF and THZ revised the manuscript. All the authors have read and approved the final manuscript.

Availability of data and materials

Less

All data generated or analyzed during this study are included in this article and its additional files.

Declarations
Ethics approval and consent to participate

Less

Ethics was approved and the informed consent was waived.

Consent for publication

Less

No individual participant data is reported that would require consent to publish from the participant (or legal parent or guardian for children).

Competing interests

Less

The authors declare that they have no competing interests.

Funding

Less

National Key Research and Development Program of China(2022YFC2010002)
National Key Research and Development Program of China(2022YFC2010000)
National Key Research and Development Program of China(2022YFC2010005)
National Natural Science Foundation of China(82171926)
National Natural Science Foundation of China(81930049)
National Natural Science Foundation of China(82202140)
Medical Imaging Database Construction Program of National Health Commission(YXFSC2022JJSJ002)
Clinical Innovative Project of Shanghai Changzheng Hospital(2020YLCYJ-Y24)
Program of Science and Technology Commission of Shanghai Municipality(21DZ2202600)
Shanghai Sailing Program(20YF1449000)

References

Less

Agustí

, Celli

, Criner

, Halpin

, Anzueto

, Barnes

, et al. Global initiative for chronic obstructive lung disease 2023 report: GOLD executive summary. Eur Respir J. 2023;61(4):2300239.

Wang

, Xu

, Yang

, Xu

, Zhang

, Bai

, et al. Prevalence and risk factors of chronic obstructive pulmonary disease in China (the China pulmonary health [CPH] study): a national cross-sectional study. Lancet. 2018;391(10131):1706-17.

Tong

, Cong

, Fang

, Fan

, Wang

, Zhao

, et al. Performance of pulmonary function test in people aged 40 years and above in China, 2019–2020. Zhonghua Liu Xing Bing Xue Za Zhi. 2023;44(5):727-34.

GOLD Global initiative for chronic obstructive lung disease—global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: 2023 report 2023. Available from: https://goldcopd.org/wp-content/uploads/2023/03/GOLD-2023-ver-1.3-17Feb2023_WMV.pdf

Mayerhoefer

, Materka

, Langs

, Haggstrom

, Szczypinski

, Gibbs

, et al. Introduction to radiomics. J Nucl Med. 2020;61(4):488-95.

, Tan

, Bourbeau

, Hogg

, Kirby

. Impact of image preprocessing methods on computed tomography radiomics features in chronic obstructive pulmonary disease. Phys Med Biol. 2021.

Yang

, Yang

, Kang

, Liang

, Wang

, Li

, et al. The value of radiomic features in chronic obstructive pulmonary disease assessment: a prospective study. Clin Radiol. 2022;77(6):e466-72.

Brown

, Kim

, Abtin

, Strange

, Galperin-Aizenberg

, Pais

, et al. Emphysema lung lobe volume reduction: effects on the ipsilateral and contralateral lobes. Eur Radiol. 2012;22(7):1547-55.

Choi

, Qi

, Yoon

, Park

, Lee

, Kim

, et al. Extension of coronavirus disease 2019 on chest CT and implications for chest radiographic interpretation. Radiol Cardiothorac Imaging. 2020;2(2):e200107.

10.

Cunliffe

, Al-Hallaq

, Labby

, Pelizzari

, Straus

, Sensakovic

, et al. Lung texture in serial thoracic CT scans: assessment of change introduced by image registration. Med Phys. 2012;39(8):4679-90.

11.

Exarchos

, Kostikas

. Artificial intelligence in COPD: possible applications and future prospects. Respirology. 2021;26(7):641-2.

12.

Vliegenthart

. Toward automated detection of chronic obstructive pulmonary disease in CT lung cancer screening. Radiology. 2023;307(5):e231350.

13.

Hofmanninger

, Prayer

, Pan

, Röhrich

, Prosch

, Langs

. Automatic lung segmentation in routine imaging is primarily a data diversity problem, not a methodology problem. Eur Radiol Exp. 2020;4(1):50.

14.

De Jay

, Papillon-Cavanagh

, Olsen

, El-Hachem

, Bontempi

, Haibe-Kains

. mRMRe: an R package for parallelized mRMR ensemble feature selection. Bioinformatics. 2013;29(18):2365-8.

15.

Zhang

, Zhang

, Cheng

, Li

, Teng

, Zhang

, et al. Artificial intelligence-driven radiomics study in cancer: the role of feature engineering and modeling. Mil Med Res. 2023;10(1):22.

16.

Remeseiro

, Bolon-Canedo

. A review of feature selection methods in medical applications. Comput Biol Med. 2019;112:103375.

17.

Amudala Puchakayala

, Sthanam

, Nakhmani

, Chaudhary

MFA

, Kizhakke Puliyakote

, Reinhardt

, et al. Radiomics for improved detection of chronic obstructive pulmonary disease in low-dose and standard-dose chest CT scans. Radiology. 2023;307(5):e222998.

18.

, Liu

, Zhang

, Yang

, Li

, Gao

, et al. A novel CT-based radiomics features analysis for identification and severity staging of COPD. Acad Radiol. 2022;29(5):663-73.

19.

Nam

, Kang

, Lee

, Kim

, Rhee

, Goo

, et al. Deep learning prediction of survival in patients with chronic obstructive pulmonary disease using chest radiographs. Radiology. 2022;305(1):199-208.

20.

Huang

, Huang

, Qiao

, Zhong

, Lu

, Wang

. Deep learning-based lung sound analysis for intelligent stethoscope. Mil Med Res. 2023;10(1):44.

21.

Song

, Zhu

, Mao

, Li

, Han

, Du

, et al. Clinical, conventional CT and radiomic feature-based machine learning models for predicting ALK rearrangement status in lung adenocarcinoma patients. Front Oncol. 2020;10:369.

22.

Park

, Park

, Lee

. Radiomics and deep learning: hepatic applications. Korean J Radiol. 2020;21(4):387-401.

23.

Cho

, Seo

, Lee

, Kim

, Yun

, Hwang

, et al. Radiomics approach for survival prediction in chronic obstructive pulmonary disease. Eur Radiol. 2021;31(10):7316-24.

24.

Do-Umehara

, Chen

, Zhang

, Misharin

, Abdala-Valencia

, Casalino-Matsuda

, et al. Epithelial cell-specific loss of function of Miz1 causes a spontaneous COPD-like phenotype and up-regulates Ace2 expression in mice. Sci Adv. 2020;6(33):7238.

25.

Miniati

, Bottai

, Pavlickova

, Monti

. Body height as risk factor for emphysema in COPD. Sci Rep. 2016;6:36896.

26.

Nacul

, Soljak

, Meade

. Model for estimating the population prevalence of chronic obstructive pulmonary disease: cross sectional data from the health survey for England. Popul Health Metr. 2007;5:8.

27.

Yun

, Cho

, Lee

, Hwang

, Lee

, Oh

, et al. Deep radiomics-based survival prediction in patients with chronic obstructive pulmonary disease. Sci Rep. 2021;11(1):15144.

28.

Occhipinti

, Paoletti

, Bartholmai

, Rajagopalan

, Karwoski

, Nardi

, et al. Spirometric assessment of emphysema presence and severity as measured by quantitative CT and CT-based radiomics in COPD. Respir Res. 2019;20(1):101.

29.

Liang

, Xu

, Wang

, Chen

, Tang

, Chen

, et al. Development of a radiomics model for predicting COPD exacerbations based on complementary visual information. Am J Respir Crit Care Med. 2021;203(9):A2296.

30.

Yang

, Li

, Guo

, Liu

, Li

, Yang

, et al. Early COPD risk decision for adults aged from 40 to 79 years based on lung radiomics features. Front Med. 2022;9:845286.

31.

Yang

, Li

, Guo

, Zeng

, Wang

, Chen

, et al. Lung radiomics features for characterizing and classifying COPD stage based on feature combination strategy and multi-layer perceptron classifier. Math Biosci Eng. 2022;19(8):7826-55.

32.

Yang

, Li

, Kang

, Guo

, Yang

, Li

, et al. A novel lung radiomics feature for characterizing resting heart rate and COPD stage evolution based on radiomics feature combination strategy. Math Biosci Eng. 2022;19(4):4145-65.

Appendix

Less

Year 2025 volume 12 Issue 1

PDF

Cite this Article

BibTeX

Article Info

doi: 10.1186/s40779-024-00516-9

Online Date：2025-11-04
Published：2025-01-10

Article Data

Affiliations

History

Funding

National Key Research and Development Program of China(2022YFC2010002)

National Key Research and Development Program of China(2022YFC2010000)

National Key Research and Development Program of China(2022YFC2010005)

National Natural Science Foundation of China(82171926)

National Natural Science Foundation of China(81930049)

National Natural Science Foundation of China(82202140)

Medical Imaging Database Construction Program of National Health Commission(YXFSC2022JJSJ002)

Clinical Innovative Project of Shanghai Changzheng Hospital(2020YLCYJ-Y24)

Program of Science and Technology Commission of Shanghai Municipality(21DZ2202600)

Shanghai Sailing Program(20YF1449000)

Affiliations