Abstract: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.

Genomic modulators of gene expression in human neutrophils

Identification of Expression Quantitative Trait Loci (eqtls) in Human Peripheral Blood Mononuclear Cells (pbmcs) and Shared with Liver and Brain

Functional annotation of genetic associations by transcriptome-wide association analysis provides insights into neutrophil development regulation

Gene Expression during the Generation and Activation of Mouse Neutrophils: Implication of Novel Functional and Regulatory Pathways

Single-cell analysis reveals the range of transcriptional states of circulating human neutrophils

Genetic regulation of B cell receptor signaling pathway: Insights from expression quantitative trait locus analysis using a mixed model

Single-cell RNA-sequencing of peripheral blood mononuclear cells reveals widespread, context-specific gene expression regulation upon pathogenic exposure

Single-cell transcriptome profiling reveals neutrophil heterogeneity in homeostasis and infection

Genetic variants associated with cell-type-specific intra-individual gene expression variability reveal new mechanisms of genome regulation

Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells

The impact of cell type and context-dependent regulatory variants on human immune traits

Impact of Genetic Polymorphisms on Human Immune Cell Gene Expression

Identification of novel biomarkers related to neutrophilic inflammation in COPD

Shared and distinct molecular effects of regulatory genetic variants provide insight into mechanisms of distal enhancer-promoter communication

Single-cell expression quantitative trait loci (eQTL) analysis of SLE-risk loci in lupus patient monocytes

Exploring key genes associated with neutrophil function and neutrophil extracellular traps in heart failure: a comprehensive analysis of single-cell and bulk sequencing data

Whole-genome Sequencing in Diverse Subjects Identifies Genetic Correlates of Leukocyte Traits: the NHLBI TOPMed Program

Identification of Target Genes at Juvenile Idiopathic Arthritis GWAS Loci in Human Neutrophils.

Complexity of the neutrophil transcriptome in early and severe rheumatoid arthritis. A role for microRNAs?

Variance quantitative trait loci reveal gene-gene interactions which alter blood traits

Distant regulatory effects of genetic variation in multiple human tissues