To further understand the effect of gene mutations on surrounding genes, we performed eQLT analysis for gene mutations with P < 0.01 and found that 24 deleterious mutations affected the expression of surrounding genes, namely TMEM17, CTD-2521M24.6, DYRK4, KLHDC4, AC007362.3, CNIH2, CTSF, SFI1, NGDN, IFNAR1, PLD4, F8A1, PRELID1, TMX4, RP11-510M2.2, STRCP1, ICAM3, SNRNP70, WWC3, DNM1P51, BIVM, LIFR, CASP1P2, and SMG8 (Table 4), these genes were the candidate genes in the pathogenesis of BPD. Here, WWC3 is linked to bronchopulmonary dysplasia.