Identification and Validation of Prognostic Markers for Lung Squamous Cell Carcinoma Associated with Chronic Obstructive Pulmonary Disease

Abstract

Background. Globally, the incidence and associated mortality of chronic obstructive pulmonary disease (COPD) and lung carcinoma are showing a worsening trend. There is increasing evidence that COPD is an independent risk factor for the occurrence and progression of lung carcinoma. This study aimed to identify and validate the gene signatures associated with COPD, which may serve as potential new biomarkers for the prediction of prognosis in patients with lung carcinoma. Methods. A total of 111 COPD patient samples and 40 control samples were obtained from the GSE76925 cohort, and a total of 4933 genes were included in the study. The weighted gene coexpression network analysis (WGCNA) was performed to identify the modular genes that were significantly associated with COPD. The KEGG pathway and GO functional enrichment analyses were also performed. The RNAseq and clinicopathological data of 490 lung squamous cell carcinoma patients were obtained from the TCGA database. Further, univariate Cox regression and Lasso analyses were performed to screen for marker genes and construct a survival analysis model. Finally, the Human Protein Atlas (HPA) database was used to assess the gene expression in normal and tumor tissues of the lungs. Results. A 6-gene signature (DVL1, MRPL4, NRTN, NSUN3, RPH3A, and SNX32) was identified based on the Cox proportional risk analysis to construct the prognostic RiskScore survival model associated with COPD. Kaplan–Meier survival analysis indicated that the model could significantly differentiate between the prognoses of patients with lung carcinoma, wherein higher RiskScore samples were associated with a worse prognosis. Additionally, the model had a good predictive performance and reliability, as indicated by a high AUC, and these were validated in both internal and external sets. The 6-gene signature had a good predictive ability across clinical signs and could be considered an independent factor of prognostic risk. Finally, the protein expressions of the six genes were analyzed based on the HPA database. The expressions of DVL1, MRPL4, and NSUN3 were relatively higher, while that of RPH3A was relatively lower in the tumor tissues. The expression of SNX32 was high in both the tumor and paracarcinoma tissues. Results of the analyses using TCGA and GSE31446 databases were consistent with the expressions reported in the HPA database. Conclusion. Novel COPD-associated gene markers for lung carcinoma were identified and validated in this study. The genes may be considered potential biomarkers to evaluate the prognostic risk of patients with lung carcinoma. Furthermore, some of these genes may have implications as new therapeutic targets and can be used to guide clinical applications.