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). This evidence concerns the gene CCL2 and bronchopulmonary dysplasia.