Bronchopulmonary dysplasia (BPD) is a disease of chronic lung dysfunction. We combined Weighted Gene Co-expression(WGCNA) and single-cell RNA analysis (scRNA-seq) to identify the biomarkers linked to BPD. Two data were obtained by searching the Gene Expression Omnibus (GEO) database with corresponding keywords, GSE190934 and GSE188944, the former for scRNA-seq and the latter for WGCNA. In the GSE190934 data, 10,762 differential genes showed up, with S100a9, S100a8, and Lyz2 being expressed primarily in pericytes. We have found that S100a8 was found to be significantly higher in the hyperoxia model mice. While Ccl5 was mainly located in CD4+ T cells, Ccl4 was mainly expressed in myeloid cells. In addition, the data has revealed that Emp2 is regarded as a marker of type I alveolar cells in BPD. Among the enriched signaling pathways by the GSE188944, several signaling communications were found to exhibit high expression patterns in BPD. The bioinformatics analysis identified some key genes and cellular communication components of BPD. These findings provide insights into the cellular and molecular underpinning of BPD pathological mechanisms.