Identification of Potential Key Genes in the Pathogenesis of Chronic Obstructive Pulmonary Disease Through Bioinformatics Analysis

Abstract

Chronic obstructive pulmonary disease (COPD) is a common respiratory disease with high morbidity and mortality. The etiology of COPD is complex, and the pathogenesis mechanisms remain unclear. In this study, we used rat and human COPD gene expression data from our laboratory and the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs) between individuals with COPD and healthy individuals. Then, protein–protein interaction (PPI) networks were constructed, and hub genes were identified. Cytoscape was used to construct the co-expressed network and competitive endogenous RNA (ceRNA) networks. A total of 198 DEGs were identified, and a PPI network with 144 nodes and 355 edges was constructed. Twelve hub genes were identified by the cytoHubba plugin in Cytoscape. Of these genes, CCR3, CCL2, COL4A2, VWF, IL1RN, IL2RA, and CCL13 were related to inflammation or immunity, or tissue-specific expression in lung tissue, and their messenger RNA (mRNA) levels were validated by qRT-PCR. COL4A2, VWF, and IL1RN were further verified by the GEO dataset GSE76925, and the ceRNA network was constructed with Cytoscape. These three genes were consistent with COPD rat model data compared with control data, and their dysregulation direction was reversed when the COPD rat model was treated with effective-component compatibility of Bufei Yishen formula III. This bioinformatics analysis strategy may be useful for elucidating novel mechanisms underlying COPD. We pinpointed three key genes that may play a role in COPD pathogenesis and therapy, which deserved to be further studied.