Aim Genetic variants contribute to the pathogenesis of bronchopulmonary dysplasia (BPD). The aim of this study is to evaluate the association of 45 SNPs with BPD susceptibility in a Turkish premature infant cohort. Methods Infants with gestational age <32 weeks were included. Patients were divided into BPD or no-BPD groups according to oxygen need at 28 days of life, and stratified according to the severity of BPD. We genotyped 45 SNPs, previously identified as BPD risk factors, in 192 infants. Results A total of eight SNPs were associated with BPD risk at allele level, two of which (rs4883955 on KLF12 and rs9953270 on CHST9) were also associated at the genotype level. Functional relationship maps suggested an interaction between five of these genes, converging on WNT5A, a member of the WNT pathway known to be implicated in BPD pathogenesis. Dysfunctional CHST9 and KLF12 variants may contribute to BPD pathogenesis through an interaction with WNT5A. Conclusions We suggest investigating the role of SNPs on different genes which are in relation with the Wnt pathway in BPD pathogenesis. We identified eight SNPs as risk factors for BPD in this study. In-silico functional maps show an interaction of the genes harboring these SNPs with the WNT pathway, supporting its role in BPD pathogenesis. Impact It is known that genetic factors may contribute to the development of BPD in preterm infants. Further studies are required to identify specific genes that play a role in the BPD pathway to evaluate them as a target for therapeutic interventions. Our study shows an association of BPD predisposition with certain polymorphisms on , , , and genes at allele level and polymorphisms on and genes at both allele and genotype level. In-silico functional mapping shows a functional relationship of these five genes with , suggesting that Wnt pathway disruption may play a role in BPD pathogenesis.