Bioinformatic analysis reveals the relationship between macrophage infiltration and Cybb downregulation in hyperoxia-induced bronchopulmonary dysplasia

Abstract Bronchopulmonary dysplasia (BPD) is the most common sequela of prematurity and is characterized by alveolar simplification and lung angiogenesis failure. The aim of this study was to explore the immune signatures of BPD. Differentially expressed gene analysis and immune infiltration analysis were conducted to identify key immune cell types and related genes by using the mRNA-seq dataset GSE25286. The expression patterns of key genes were validated in the scRNA-seq dataset GSE209664 and in experiments. The cell–cell crosstalk of key immune cells was explored with CellChat. We found that differentially expressed genes between BPD mice and controls were mostly enriched in leukocyte migration and M1 macrophages were highly enriched in BPD lungs. Hub genes (Cybb, Papss2, F7 and Fpr2) were validated at the single-cell level, among which the downregulation of Cybb was most closely related to macrophage infiltration. The reduced mRNA and protein levels of Cybb were further validated in animal experiments. Colocalization analysis of Cybb and macrophage markers demonstrated a significant decrease of Cybb in M1 macrophages. Cell–cell crosstalk found that alveolar epithelial cells interacted actively with macrophages through MIF-(CD74 + CD44) signalling. In conclusion, M1 macrophages played important roles in promoting BPD-like lung injury, which was correlated with a specific reduction of Cybb in macrophages and the potential activation of MIF signalling.