Ferroptosis-Related Genes in Bronchoalveolar Lavage Fluid Serves as Prognostic Biomarkers for Idiopathic Pulmonary Fibrosis

Li, Meng; Wang, Ke; Zhang, Yanpeng; Fan, Meng; Li, Anqi; Zhou, Jiejun; Yang, Tian; Shi, Puyu; Li, Dan; Zhang, Guangjian; Chen, Mingwei; Ren, Hui

Background: Idiopathic pulmonary fibrosis (IPF) is a chronic progressive disease with unknown etiology and unfavorable prognosis. Ferroptosis is a form of regulated cell death with an iron-dependent way that is involved in the development of various diseases. Whereas the prognostic value of ferroptosis-related genes (FRGs) in IPF remains uncertain and needs to be further elucidated. Methods: The FerrDb database and the previous studies were screened to explore the FRGs. The data of patients with IPF were obtained from the GSE70866 dataset. Wilcoxon's test and univariate Cox regression analysis were applied to identify the FRGs that are differentially expressed between normal and patients with IPF and associated with prognosis. Next, a multigene signature was constructed by the least absolute shrinkage and selection operator (LASSO)-penalized Cox model in the training cohort and evaluated by using calibration and receiver operating characteristic (ROC) curves. Then, 30% of the dataset samples were randomly selected for internal validation. Finally, the potential function and pathways that might be affected by the risk score-related differently expressed genes (DEGs) were further explored. Results: A total of 183 FRGs were identified by the FerrDb database and the previous studies, and 19 of them were differentially expressed in bronchoalveolar lavage fluid (BALF) between IPF and healthy controls and associated with prognosis (p < 0.05). There were five FRGs (aconitase 1 [ACO1], neuroblastoma RAS viral (v-ras) oncogene homolog [NRAS], Ectonucleotide pyrophosphatase/phosphodiesterase 2 [ENPP2], Mucin 1 [MUC1], and ZFP36 ring finger protein [ZFP36]) identified as risk signatures and stratified patients with IPF into the two risk groups. The overall survival rate in patients with high risk was significantly lower than that in patients with low risk (p < 0.001). The calibration and ROC curve analysis confirmed the predictive capacity of this signature, and the results were further verified in the validation group. Risk score-related DEGs were found enriched in ECM-receptor interaction and focal adhesion pathways. Conclusion: The five FRGs in BALF can be used for prognostic prediction in IPF, which may contribute to improving the management strategies of IPF.