To elucidate the characteristics of lymphocyte subsets in bronchopulmonary dysplasia (BPD) diagnosis following Jensen's criterion to understand the spectrum of lymphocytes in different degrees of BPD. This single-center retrospective cohort study included 120 neonates admitted to the neonatal intensive care unit between 1 July 2014 and 30 June 2021, who had undergone peripheral blood lymphocyte subpopulation detection. Thirty-one neonates were included in the control group, whereas 33 infants with BPD were included in the case group. In addition, we selected 56 infants with a gestational age (GA) <37 weeks without BPD who were receiving oxygen therapy. Among the three groups, the B cell and NK cell frequencies were significantly higher and the frequencies of T cells and CD4+ cells were significantly lower in the BPD group. In newborns without BPD, the distribution of T lymphocyte subsets was similar at different GAs. Comparing different degrees of BPD, the patients in the grades 2-3 BPD group had significantly lower percentages of T lymphocytes and CD4+ T cells than those in the other groups. Remarkably, the frequencies of NK cells were significantly higher in patients with grades 2-3 BPD, and the Treg cells slightly increased with BPD severity, although the differences were not significant. Healthy neonates had similar ratios of lymphocyte subsets among different GAs; although as the GAs increased, the percentage of lymphocytes increased slightly. Severe BPD was associated with lower CD4+ T cells and higher NK cells. However, whether such changes were the cause or the consequence of BPD has not been determined.