Identification of Hub Genes in the Pathogenesis of Bronchiolitis Obliterans via Bioinformatic Analysis and Experimental Verification
Zhongji Wu,Xiaowen Chen,Kangkang Zhang,Zhenwei Liu,Haidi Zhang,Zhaocong Zheng,Xiaodie Zhang,Yubiao Chen,Yinghui Peng,Hui Li,Kaiyin Huang,Sixiang Tang,Li Zhao,Dehui Chen
DOI: https://doi.org/10.2147/jir.s419845
IF: 4.5
2023-08-08
Journal of Inflammation Research
Abstract:Zhongji Wu, 1, &ast Xiaowen Chen, 1, &ast Kangkang Zhang, 1 Zhenwei Liu, 1 Haidi Zhang, 1 Zhaocong Zheng, 2 Xiaodie Zhang, 2 Yubiao Chen, 3 Yinghui Peng, 1 Hui Li, 1 Kaiyin Huang, 1 Sixiang Tang, 1 Li Zhao, 1 Dehui Chen 1 1 The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, 510000, People's Republic of China; 2 Guangzhou Medical University, Guangzhou, 510000, People's Republic of China; 3 State Key Laboratory of Respiratory Diseases, Guangzhou, 510000, People's Republic of China &astThese authors contributed equally to this work Correspondence: Dehui Chen, 151 Yanjiang West Road, Yuexiu District, Guangzhou City, Guangdong Province, 510000, People's Republic of China, Tel +86 13922732810, Email 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 development of distinct BO diagnostic markers and novel therapeutic avenues. Keywords: bronchiolitis obliterans, biomarkers, pathways, transcriptome, therapeutic targets Bronchiolitis obliterans (BO) is a chronic pulmonary disease characterized by irreversible narrowing of the small airways, which leads to reduced long-term survival after lung transplantation. 1,2 It is also a common complication in children diagnosed with severe adenovirus pneumonia. 3 While certain diagnosis and treatments were available for BO patients, significant challenges were still present for patients complete recovery, hence the poor prognosis for BO patients. Therefore, efforts to understand the underlying molecular mechanisms of BO may be crucial to provide a significant breakthrough in developing better treatments for patients with BO. With the current advancements in information technology, data from pooled biological samples can be easily accessed through bioinformatics approaches. Bioinformatics is an interdisciplinary field that combines life sciences and computer science to analyze genomic, transcriptomic, proteomic, metabolomic, and other sequencing data that may reveal the underlying biological significance of a particular gene or regulatory pathways. 4 Previous studies that rely on transcriptomic data analysis using bioinformatics approaches have also succeeded in yielding promising results, specifically by exploring BO molecular mechanisms. 5–7 However, identification and verification of distinct molecular markers are still lacking, primarily because of limitations such as small sample sizes, inadequate transcriptome data analysis, and the absence of verification processes. To address these limitations, multiple BO-related transcriptomic datasets were integrated from the Gene Expression Omnibus (GEO) database. We sought to explore the underlying molecular mechanisms of BO and subsequently identify potential biomarkers by employing strategies that incorporate bioinformatic analysis from multiple facets. Our study also utilized two distinct BO experimental mouse models that were commonly constructed to generate diseased mouse models. 8,9 Total RNA was then extracted, and expression levels for our target genes of interest were examined in hopes t -Abstract Truncated-
immunology