Identification of Hub Genes in the Pathogenesis of Bronchiolitis Obliterans via Bioinformatic Analysis and Experimental Verification

Background: Bronchiolitis obliterans (BO) is a chronic disease that can arise as a complication of severe childhood pneumonia and can also impact the long-term survival of patients after lung transplantation. However, the precise molecular mechanism underlying BO remains unclear. We aimed to identify BO-associated hub genes and their molecular mechanisms. Methods: BO-associated transcriptome datasets (GSE52761, GSE137169, and GSE94557) were downloaded from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs). Additional bioinformatics analyses, such as Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Protein-Protein Interaction (PPI) analyses, were performed to determine functional roles and DEG-associated regulatory networks. Prediction of hub genes using the 12 algorithms available in the Cytohubba plugin of Cytoscape software was also performed. Verification was performed using the BO mouse model. Results: Our results revealed 57 DEGs associated with BO, of which 18 were down-regulated and 39 were up-regulated. The Cytohubba plugin data further narrowed down the 57 DEGs into 9 prominent hub genes (CCR2, CD1D, GM2A, TFEC, MPEG1, CTSS, GPNMB, BIRC2, and CTSZ). Genes such as CCR2, TFEC, MPEG1, CTSS, and CTSZ were dysregulated in 2,3-butanedione-induced BO mice, whereas TFEC, CTSS, and CTSZ were dysregulated in nitric acid-induced BO mouse models. Conclusion: Our study identified and validated four novel BO biomarkers, which may allow further investigation into the develop-ment of distinct BO diagnostic markers and novel therapeutic avenues.