Transcriptome analysis demonstrated genes predicting BPD to be overrepresented in inflammatory pathways with increased disease severity characterized by the regulation of immune and defense response pathways and upstream regulator analysis confirmed TNF-α, interleukin (IL) -6, and interferon α as the highest activated cytokines in more severe disease. Here, TNF is linked to bronchopulmonary dysplasia.