BACKGROUND: Chronic obstructive pulmonary disease (COPD) is one of the most prevalent chronic respiratory diseases and the fourth cause of mortality globally. Neutrophilic inflammation has a vital role in the occurrence and progression of COPD. This study aimed to identify the novel hub genes involved in neutrophilic inflammation in COPD through bioinformatic prediction and experimental validation.
METHODS: Both the single-cell RNA sequencing (scRNA-seq) dataset (GSE173896) and the RNA sequencing (RNA-seq) dataset (GSE57148) were downloaded from the Gene Expression Omnibus (GEO) database. The Seurat package was used for quality control, dimensions reduction, and cell identification of scRNA-seq. The irGSEA package was used for scoring individual cells. The Monocle2 package was used for the trajectory analysis of neutrophils. The CIBERSORT algorithm was used for analysis of immune cell infiltration in the lungs of COPD patients and controls in RNA-seq dataset, and weighted gene co-expression network analysis (WGCNA) correlated gene modules with neutrophil infiltration. The Mendelian randomization (MR) analysis explored the causal relationship between feature DEGs and COPD. The protein-protein interaction (PPI) network of novel hub genes was constructed, and real-time quantitative polymerase chain reaction (qRT-PCR) was used to validate novel hub genes in clinical specimens.
RESULTS: In scRNA-seq, the gene sets upregulated in COPD samples were related to the neutrophilic inflammatory response and TNF-α activation of the NF-κB signaling pathway. In RNA-seq, immune infiltration analysis showed neutrophils were upregulated in COPD lung tissue. We combined data from differential and modular genes and identified 51 differential genes associated with neutrophilic inflammation. Using MR analysis, 6 genes were explored to be causally associated with COPD. Meanwhile, 11 hub genes were identified by PPI network analysis, and all of them were upregulated. qRT-PCR experiments validated 9 out of 11 genes in peripheral blood leukocytes of COPD patients. Furthermore, 5 genes negatively correlated with lung function in COPD patients. Finally, a network of transcription factors for NAMPT and PTGS2 was constructed.
CONCLUSION: This study identified nine novel hub genes related to the neutrophilic inflammation in COPD, and two genes were risk factors of COPD, which may serve as potential biomarkers for the clinical severity of COPD.
