Identification of accurate biomarkers is still particularly urgent for improving the poor survival of chronic obstructive pulmonary disease (COPD) patients. In this investigation, we aimed to identity the potential biomarkers in COPD via bioinformatics and next generation sequencing (NGS) data analysis. In this investigation, the differentially expressed genes (DEGs) in COPD were identified using NGS dataset (GSE239897) from Gene Expression Omnibus (GEO) database. Subsequently, gene ontology (GO) and pathway enrichment analysis was conducted to evaluate the underlying molecular mechanisms involved in progression of COPD. Protein-protein interaction (PPI), modules, miRNA-hub gene regulatory network and TF-hub gene regulatory network analysis were performed to determine the hub genes, miRNAs and TFs. The receiver operating characteristic (ROC) analysis was performed to determine the diagnostic value of hub genes. A total of 956 overlapping DEGs (478 up regulated and 478 down regulated genes) were identified in the NGS dataset. DEGs were mainly associated with GO functional terms and pathways in cellular response to stimulus. response to stimulus, immune system and neutrophil degranulation. There were 10 hub genes (MYC, LMNA, VCAM1, MAPK6, DDX3X, SHMT2, PHGDH, S100A9, FKBP5 and RPS6KA2) identified by PPI, modules, miRNA-hub gene regulatory network and TF-hub gene regulatory network analysis. In conclusion, the DEGs, relative GO terms, pathways and hub genes identified in the present investigation might aid in understanding of the molecular mechanisms underlying COPD progression and provide potential molecular targets and biomarkers for COPD.