Multivariate modeling revealed FGF2, sIL-2 Rα, MCP-1, MIP1a, and TNF-α as predictors of BPD when considering GA as a strong influencing factor (VIP). Here, CCL3 is linked to bronchopulmonary dysplasia.