Bioinformatics integration reveals key genes associated with mitophagy in myocardial ischemia-reperfusion injury

Abstract

Background Myocardial ischemia is a prevalent cardiovascular disorder associated with significant morbidity and mortality. While prompt restoration of blood flow is essential for improving patient outcomes, the subsequent reperfusion process can result in myocardial ischemia–reperfusion injury (MIRI). Mitophagy, a specialized autophagic mechanism, has consistently been implicated in various cardiovascular disorders. However, the specific connection between ischemia–reperfusion and mitophagy remains elusive. This study aims to elucidate and validate central mitophagy-related genes associated with MIRI through comprehensive bioinformatics analysis. Methods We acquired the microarray expression profile dataset (GSE108940) from the Gene Expression Omnibus (GEO) and identified differentially expressed genes (DEGs) using GEO2R. Subsequently, these DEGs were cross-referenced with the mitophagy database, and differential nucleotide sequence analysis was performed through enrichment analysis. Protein–protein interaction (PPI) network analysis was employed to identify hub genes, followed by clustering of these hub genes using cytoHubba and MCODE within Cytoscape software. Gene set enrichment analysis (GSEA) was conducted on central genes. Additionally, Western blotting, immunofluorescence, and quantitative polymerase chain reaction (qPCR) analyses were conducted to validate the expression patterns of pivotal genes in MIRI rat model and H9C2 cardiomyocytes. Results A total of 2719 DEGs and 61 mitophagy-DEGs were identified, followed by enrichment analyses and the construction of a PPI network. HSP90AA1, RPS27A, EEF2, EIF4A1, EIF2S1, HIF-1α, and BNIP3 emerged as the seven hub genes identified by cytoHubba and MCODE of Cytoscape software. Functional clustering analysis of HIF-1α and BNIP3 yielded a score of 9.647, as determined by Cytoscape (MCODE). In our MIRI rat model, Western blot and immunofluorescence analyses confirmed a significant elevation in the expression of HIF-1α and BNIP3, accompanied by a notable increase in the ratio of LC3II to LC3I. Subsequently, qPCR confirmed a significant upregulation of HIF-1α, BNIP3, and LC3 mRNA in the MIRI group. Activation of the HIF-1α/BNIP3 pathway mediates the regulation of the degree of Mitophagy, thereby effectively reducing apoptosis in rat H9C2 cardiomyocytes. Conclusions This study has identified seven central genes among mitophagy-related DEGs that may play a pivotal role in MIRI, suggesting a correlation between the HIF-1α/BNIP3 pathway of mitophagy and the pathogenesis of MIRI. The findings highlight the potential importance of mitophagy in MIRI and provide valuable insights into underlying mechanisms and potential therapeutic targets for further exploration in future studies.