OBJECTIVETo explore the mechanism underlying the therapeutic effect of Bufei Yishen Formula III combined with exercise rehabilitation (ECC-BYF III + ER) on chronic obstructive pulmonary disease (COPD) and further identify hub genes.MATERIALS AND METHODSGene Set Enrichment Analysis was used to identify the COPD-associated pathways and reversal pathways after ECC-BYF III + ER treatment. Protein-protein interaction network analysis and cytoHubba were used to identify the hub genes. These genes were verified using independent datasets, molecular docking and quantitative real-time polymerase chain reaction experiment.RESULTSUsing the high-throughput sequencing data of COPD rats from our laboratory, 49 significantly disturbed pathways were identified in COPD model compared with control group via gene set enrichment analysis (false discovery rate < 0.05). The 34 pathways were reversed after ECC-BYF III + ER treatment. In the 2306 genes of these 34 pathways, 121 of them were differentially expressed in COPD rats compared with control samples. A protein-protein interaction network comprising 111 nodes and 274 edges was created, and 34 candidate genes were identified. Finally, seven COPD hub genes (Il1b, Ccl2, Cxcl1, Apoe, Ccl7, Ccl12, and Ccl4) were well identified and verified in independent COPD rat data from our laboratory and the public dataset GSE178513. The area under the receiver operating characteristic curve values ranged from 0.86 to 1 and from 0.67 to 1, respectively. The reliability of the mentioned genes, which can bind to the active ingredients of ECC-BYF III through molecular docking, were further verified through qRT-PCR experiments.CONCLUSIONThirty-four COPD-related pathways and seven hub genes that may be regulated by ECC-BYF III + ER were identified and well verified. The findings of this study may provide insights into the treatment and mechanism underlying COPD.