Integrated Bioinformatics Analysis and Verification of Gene Targets for Myocardial Ischemia-Reperfusion Injury

Wang, Jianru; Li, Xiaohui; Peng, Guangcao; Fan, Genhao; Zhang, Mengmeng; Chen, Jian

doi:https://doi.org/10.1155/2022/2056630

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Therapeutic Strategies for the Treatment of Cardio/Cerebrovascular-Related Diseases

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 2056630 | https://doi.org/10.1155/2022/2056630

Integrated Bioinformatics Analysis and Verification of Gene Targets for Myocardial Ischemia-Reperfusion Injury

Jianru Wang,¹Xiaohui Li,¹Guangcao Peng,¹Genhao Fan,¹Mengmeng Zhang,¹and Jian Chen^2,3

Academic Editor: Thanasekaran Jayakumar

Received16 Sept 2021

Revised11 Mar 2022

Accepted28 Mar 2022

Published15 Apr 2022

Abstract

Background. Myocardial ischemia-reperfusion injury (MIRI) has become a thorny and unsolved clinical problem. The pathological mechanisms of MIRI are intricate and unclear, so it is of great significance to explore potential hub genes and search for some natural products that exhibit potential therapeutic efficacy on MIRI via targeting the hub genes. Methods. First, the differential expression genes (DEGs) from GSE58486, GSE108940, and GSE115568 were screened and integrated via a robust rank aggregation algorithm. Then, the hub genes were identified and verified by the functional experiment of the MIRI mice. Finally, natural products with protective effects against MIRI were retrieved, and molecular docking simulations between hub genes and natural products were performed. Results. 230 integrated DEGs and 9 hub genes were identified. After verification, Emr1, Tyrobp, Itgb2, Fcgr2b, Cybb, and Fcer1g might be the most significant genes during MIRI. A total of 75 natural products were discovered. Most of them (especially araloside C, glycyrrhizic acid, ophiopogonin D, polyphyllin I, and punicalagin) showed good ability to bind the hub genes. Conclusions. Emr1, Tyrobp, Itgb2, Fcgr2b, Cybb, and Fcer1g might be critical in the pathological process of MIRI, and the natural products (araloside C, glycyrrhizic acid, ophiopogonin D, polyphyllin I, and punicalagin) targeting these hub genes exhibited potential therapeutic efficacy on MIRI. Our findings provided new insights to explore the mechanism and treatments for MIRI and revealed new therapeutic targets for natural products with protective properties against MIRI.

1. Introduction

Coronary heart disease (CHD) as a global public health problem is one of the main causes of decreased quality of life and death worldwide [1]. Acute myocardial infarction (AMI) is a kind of acute and critical manifestation of CHD with high morbidity and mortality. Early rapid restoration of coronary blood flow can reduce the size of myocardial infarction and prevent further cardiac injury. With the development of reperfusion strategies such as thrombolysis and percutaneous coronary intervention, the mortality rate of AMI has been greatly reduced [2]. However, myocardial reperfusion can further cause secondary myocardial injury, which can account for 50% of the final myocardial infarct area [3]. This phenomenon, termed myocardial ischemia-reperfusion injury (MIRI), has become an unsolved clinical problem and a major cause of cardiac insufficiency. Currently, there is no specific treatment for MIRI in pharmacopeia [4]. Therefore, clinical and basic studies in the prevention and treatment of MIRI are active areas in the medical field.

It is currently believed that the pathogenesis of MIRI is a multi-factorial complex process. The latest evidence indicates that epigenetic regulation including histone modification, DNA methylation, non-coding RNAs, and N6-methyladenosine (m6A) methylation plays a key role in MIRI and can be used as new therapeutic targets for MIRI [5]. Iron, as an essential mineral, is an important player in the physiological functions of many tissues and organs in the human body. Studies in recent years have shown that imbalance in iron metabolism homeostasis is intimately related to the pathological process of MIRI [6]. Ferroptosis, a newly discovered form of regulated cell death, results from the imbalance in iron metabolism [6]. The occurrence of ferroptosis plays an important role during MIRI and can serve as a potential therapy for MIRI [7]. Furthermore, there are numerous other reported pathological mechanisms of MIRI including immunoreaction, autophagy, cell-derived exosomes, dysfunctional mitochondria, different forms of cell death, inflammation, oxidative stress, intracellular calcium overload, and so on [4, 8–10]. Although the pathological mechanisms of MIRI have been extensively explored, its specific pathogenesis has not been fully elucidated. Therefore, plentiful work is still required to better understand the potential pathologic mechanism of MIRI at the molecular level, which may be imperative to develop more effective treatments and seek early diagnostic markers.

In recent years, with the help of bioinformatics and high-throughput techniques, a comprehensive understanding of the molecular mechanism and progression of various diseases has become possible [11]. Molecular docking, a very popular and useful tool, is used to research the interaction between a small-molecule ligand and a target receptor in the drug discovery arena [12]. In our previous study, the molecular docking technique was successfully used to predict the binding modes between 12 COVID-19 targets and the 20 compounds of Qingfei Paidu decoction that have been widely used to treat COVID-19 in China [13]. The analyses of microarray and RNA-seq gene expression data by bioinformatics have been extensively used not only to explore the diagnostic and prognostic biomarkers but also to probe crucial genes and biological processes in many diseases [14–17]. However, some reasons, including the limited sample quantities, the heterogeneity of experimental samples, the use of different detection platforms, and so on, may generate deviation in the results [18, 19]. The integration of the results of multiple gene expression data sets by a robust rank aggregation (RRA) algorithm based on a statistical model that naturally allows evaluating the significance of the results is a promising strategy to overcome these shortcomings [20]. RRA has been widely used in the comprehensive analysis of multiple data sets for both oncology and non-oncology diseases [18, 19, 21]. However, there are no reports on the use of the RRA algorithm in MIRI.

In this study, expression profile analyses and bioinformatics methods were combined to explore hub genes and their functions in MIRI. We subsequently performed the validation of the results via the functional experiment of the MIRI mice. Meanwhile, we retrieved the natural products in Chinese herbal medicine exerting protection for MIRI from the PubMed database. Finally, molecular docking simulations between natural products and hub genes were performed (Figure 1). To sum up, the ultimate goal of the research is to provide novel insights into new potential therapeutic targets and the pathogenesis of MIRI and develop more effective anti-MIRI drugs.

Figure 1

Flowchart of the study design. First, the differential expression genes (DEGs) from three arrays were screened for myocardial ischemia-reperfusion injury (MIRI) and integrated to explore the potential pathogenesis of MIRI using the RRA algorithm. Second, protein-protein interaction network and transcription factors analysis of the integrated DEGs were performed. Then, the most significant DEGs were identified and verified by the functional experiment of the MIRI mice. Finally, natural products exerting protection for MIRI were retrieved and molecular docking simulations were performed.

2. Materials and Methods

2.1. Affymetrix Microarray Data

We collected the expression profile data of GSE58486 (GPL18802), GSE108940 (GPL7202), and GSE115568 (GPL16570) from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). GSE58486 contains six MIRI samples and three sham surgery samples. GSE108940 contains six MIRI samples and six sham surgery samples. GSE115568 contains three MIRI samples and three sham surgery samples. All samples were collected from murine hearts.

2.2. Screening for Differential Expression Genes (DEGs) and Functional Enrichment Analyses

GEO2R online analysis (https://www.ncbi.nlm.nih.gov/geo/geo2r/) was firstly used to obtain the .txt files of the DEGs for the three microarray data sets. The -value <0.01 and the absolute fold-change (FC) >1.5 were set to determine DEGs. Then, multi-list enrichment analysis of the three gene lists was performed using Metascape (https://metascape.org) [22]. , a minimum count of 3, and the enrichment factor >1.5 were set as the thresholds for enrichment analyses.

2.3. Integration of DEGs in Three Data Sets

The TXT files of all gene lists ranked by logFC in three data sets were collected for RRA analysis (https://www.icesi.edu.co/CRAN/web/packages/RobustRankAggreg/) in R version 3.5.3 (https://cran.rproject.org/), and the -value <0.01 was set to define the significant integrated DEGs. The code used for the RRA analysis is presented in Supplementary Table 1. Gene ontology (GO) and pathway enrichment analyses of integrated DEGs were performed using Metascape. For pathway enrichment analyses, the KEGG pathway, Reactome gene sets, and WikiPathways were employed. The bubble plots were drawn with the ggplot2 package.

2.4. PPI Network and TFs Analysis

The interaction network of the integrated DEGs was performed based on the Search Tool for the Retrieval of Interacting Genes (STRING) [23] to mine their target genes, and interaction pairs with a combined score of >0.4 were considered significant. Then, Cytoscape (version 3.2.1) [24] was applied to visualize the molecular interactions. Next, MCODE analysis [25] was carried out to find densely connected modules in the PPI network. Finally, transcription factors (TFs) analysis of the MCODE modules were analyzed by the iRegulon plugin [26] of Cytoscape with the default criteria. The top five TFs of each module with the higher normalized enrichment scores (NES > 4) were listed.

2.5. Hub Protein Analysis

CytoHubba app [27] was used to identify the top 20 ranking genes using the degree, closeness, edge percolated component (EPC), maximum neighborhood component (MNC), and maximal clique centrality (MCC). The overlapping genes from five algorithms were taken as hub genes. At last, the functional enrichment analyses of hub genes were performed in Metascape, using GO biological processes, KEGG pathway, Reactome gene sets, and WikiPathways.

2.6. Establishment of MIRI Model in Mice

The animal experiments were approved by the Experimental Animal Welfare Ethics Review Committee of the First Affiliated Hospital of Henan University of CM (approval number, YFYDW2020004) and followed the National Institutes of Health Guide for the Care and Use of Laboratory Animals (NIH Publications No. 8023, revised 1978). Male C57BL/6 J mice (6–8 weeks) were purchased from Huaxing Experimental Animal Farm (Zhengzhou, China; qualified production number, SCXK-(Yu)-2019-0002). The mice were fed under a 12 h cycle of light/dark in IVC condition and had free access to food and water. Before modeling, the mice were randomly divided into two groups (n = 9/group), namely sham group and MIRI group. MIRI mouse model was established by ligating the left anterior descending coronary artery, as previously described [28]. After ischemia for 30 min, the ligation was released to allow reperfusion. Mice were sacrificed following reperfusion 24 h. Sham‐operated mice underwent the same procedure without ligation. Blood samples were collected for measurement of lactate dehydrogenase (LDH) activities using a kit (Nanjing Jiancheng Bioengineering Institute, Nanjing, China) following the manufacturer’s instructions. The infarct area was determined by staining the heart tissues with 2% triphenyl tetrazolium chloride (TTC; Beijing Solarbio Science & Technology Co. Ltd., Beijing, China) for 15 min at 37°C. The infarct size was evaluated according to the method reported in previous literature by Image-Pro Plus 6.0 software (Media Cybernetics, Silver Spring, MD, USA) [29, 30].

2.7. Verifying of Hub Genes by Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR)

Total RNA was extracted from each myocardial tissue using TRIzol (Invitrogen Corporation, CA, USA). cDNA was synthesized using NovoScript^® Plus1st Strand cDNA Synthesis SuperMix (Novoprotein Scientific Inc., Shanghai, China). RT-qPCR was performed using SYBR High-Sensitivity qPCR SuperMix (Novoprotein Scientific Inc., Shanghai, China). The primers were synthesized by Shanghai Sangon Biological Engineering Technology (Shanghai, China; Table 1). The mRNA expression level was normalized to that of β-actin in the same sample. The relative mRNA expression level was calculated using the 2^−ΔΔCT method.

2.8. Molecular Docking Simulation

The natural products in Chinese herbal medicine exerting protection for MIRI were retrieved from the PubMed database (https://pubmed.ncbi.nlm.nih.gov/). Then, the 2D chemical structures of natural products were obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/). The 3D structures of hub genes were obtained from the Protein Data Bank (PDB) database (https://www.rcsb.org/). Molecular docking simulation of each hub gene with natural products was performed using AutoDock Vina software.

2.9. Statistical Analysis

All data were presented as mean ± SEM and analyzed using IBM SPSS statistics 21.0 software. The Student’s t-test was utilized to compare data between two groups. A value < 0.05 was accepted as statistically significant.

3. Results

3.1. Identification and Functional Enrichment Analyses of DEGs

After screening with the threshold of |FC| >1.5 and , 1986, 1703, and 402 DEGs were identified in the GSE58486, GSE108940, and GSE115568, respectively (Figure 2(a)–2(c)). To facilitate the understanding of the functional mechanisms that are shared between, or selectively ascribed to, specific gene lists, a multi-gene-list meta-analysis was performed. Figure 2(d) shows that eight GO terms were common in all studies. Figure 2(e) shows that eight pathways were common in all studies.

(a)

(b)

(c)

(d)

(e)

3.2. Identification of Integrated DEGs in MIRI

Table 2 presents that 230 integrated DEGs were identified using the RRA (). The top 20 up- and downregulated hub genes were drawn into a heatmap (Figure 3(a)). The top five BP terms showed that integrated DEGs were mainly related to inflammatory response, cell chemotaxis, and so on (Figure 3(b)). The top five CC terms showed that integrated DEGs were mainly related to the extracellular matrix, external encapsulating structure, and so on (Figure 3(b)). The top five MF terms showed that integrated DEGs were mainly related to extracellular matrix structural constituent, cell adhesion molecule binding, and so on (Figure 3(b)). We enriched 35, 55, and 15 pathways from the KEGG, Reactome, and WikiPathways databases, respectively, setting a -value < 0.01, a minimum count of 3, and an enrichment factor >1.5 as screening thresholds. The top 10 enriched Reactome pathways closely associated with MIRI were mainly neutrophil degranulation, extracellular matrix organization, collagen degradation, and so on (Figure 3(c)). The top 10 enriched KEGG pathways closely associated with MIRI were mainly IL-17 signaling pathway, AGE-RAGE signaling pathway, ECM-receptor interaction, and so on (Figure 3(d)). The top 10 enriched WikiPathways closely associated with MIRI were mainly focal adhesion, PI3K-Akt-mTOR signaling pathway, inflammatory response pathway, and so on (Figure 3(e)).

(a)

(b)

(c)

(d)

(e)

Figure 3

Identification and functional enrichment analyses of integrated DEGs. (a) Heatmap of top 20 up- and downregulated DEGs of three microarrays screening by RRA. Red boxes represent upregulated genes and green boxes represent downregulated genes. The values in the boxes represent logFC values. (b) The top five enriched GO terms for integrated DEGs in biological process, molecular function, and cellular component categories. (c) The top ten enriched Reactome pathways for integrated DEGs. (d): The top ten enriched KEGG pathways for integrated DEGs. (e) The top ten enriched WikiPathways for integrated DEGs.

3.3. PPI Network and TFs Analysis

There were 193 nodes and 1,203 edges in the PPI network of the 230 integrated DEGs in MIRI (Figure 4). MCODE analysis showed that four meaningful functional modules were selected (score > 5). The top 5 TFs were predicted to target these MCODE modules, including Pura, Srf, Dbp, Stat4, and so on (Figure 5(a)–5(d)).

(a)

(b)

(c)

(d)

3.4. Hub Protein Analysis

The hub genes were identified by CytoHubba from 230 integrated DEGs. The top 20 ranking genes, which were selected based on closeness, degree, EPC, MCC, and MNC, are displayed in Figure 6(a)–6(e). There were 9 intersection genes (Tyrobp, Itgb2, Emr1, Fcer1g, Fcgr1, Fcgr2b, Cybb, Cd68, and Ptprc) among 5 algorithms (Figure 6(f)). These genes, also known as hub genes, were selected for further validation in the MIRI model. In addition, Metascape results showed that 9 genes were mainly related to neutrophil degranulation, inflammatory response, superoxide anion generation, leukocyte mediated cytotoxicity, and so on (Figure 6(g)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

3.5. The Successful Establishment of the MIRI Model in Mice

As illustrated in Figure 7(a)–7(b), the infarct areas of hearts in the MIRI group were significantly larger than those in the sham group. Moreover, the serum LDH activities were significantly increased in mice subjected to MIRI (Figure 7(c)). Taken together, the obtained data suggested that the MIRI mouse model was successfully established.

(a)

(b)

(c)

(d)

3.6. The mRNA Expression Levels of Six Hub Genes Were Upregulated in the MIRI Model

To verify the results of bioinformatics analysis, we detected the mRNA expression level of hub genes by RT-qPCR in the myocardial tissues. Compared with the sham group, the mRNA expression levels of Emr1, Tyrobp, Itgb2, Fcgr2b, Cybb, and Fcer1g were higher in the MIRI group (Figure 7(d)). Unexpectedly, CD68, Ptprc, and Fcgr1 were not markedly differentially expressed. However, the changing trend was consistent with the result of bioinformatics results (Figure 7(d)).

3.7. Molecular Docking Results

A total of 75 natural products were discovered, and their multiple mechanisms of action for treating MIRI were complex (Supplementary Table 2). To verify the effect of natural products on hub genes, molecular docking was performed. Unfortunately, the 3D structure of Emr1 was not obtained from the PDB database. The molecular docking results showed that there were different degrees of binding between natural products and hub genes, and most of the natural products showed a good ability to bind the hub genes, particularly araloside C, glycyrrhizic acid, ophiopogonin D, polyphyllin I, and punicalagin (Figure 8).

4. Discussion

MIRI, as an inevitable phenomenon, is a major challenge to the treatment of AMI and plays a crucial role in the damage, repair, and remodeling of myocardial tissue. MIRI not only can contribute to myocardial infarction but also can lead to myocardial stunning, the no-reflow phenomenon, reperfusion arrhythmia, lethal reperfusion injury, coronary capillary rupture, and hemorrhages [10, 31]. Despite remarkable progress in understanding the pathogenesis of MIRI, the optimal treatment strategies to effectively restrict MIRI remain elusive [32]. It has been very disappointing to translate effective strategies that have been proven in many preclinical studies into better clinical outcomes [32]. There is an undisputed unmet need between therapy and understanding the pathophysiological mechanism of MIRI. Only with a better understanding of the mechanisms of MIRI, more appropriate strategies to reduce myocardial damage can be developed. Therefore, researching the mechanisms underlying MIRI development is very important.

In this study, GSE58486, GSE108940, and GSE115568 data sets were firstly downloaded from the GEO database, and DEGs were determined. Then, the multi-gene-list meta-analyses were performed to find the common GO and KEGG terms of the DEGs. Moreover, the three data sets were analyzed using the RRA method, and integrated DEGs were found. Function annotation and pathway enrichment analyses of integrated DEGs were then conducted. The results of the above functional enrichment analyses showed that DEGs might be mainly involved in regulating the extracellular matrix, cell adhesion, cell chemotaxis, inflammatory response, AGE-RAGE signaling pathway, and PI3K-Akt-mTOR signaling pathway. MIRI induces a series of sterile inflammatory reactions, which can lead to further myocardial injury. Many risk factors, including chemokines, inflammatory cytokines, inflammatory cells, complement cascades, inflammatory pathways, and so on, play important roles in the inflammatory response during MIRI through multiple interacting pathways [33, 34]. Adhesion molecules are transmembrane glycoproteins with a variety of biological activities, which are produced by cells. They can mediate the interaction and binding between cells and cells and between cells and extracellular matrix. Adhesion molecules E-selectin and intercellular cell adhesion molecule-1 expressed on activated endothelial cells can regulate the leukocyte adhesion cascade and thus aggravate tissue inflammation in the pathological process of MIRI [35]. A study has shown that HMGB1/TLR4 signaling, as an inflammatory signaling pathway, promotes the release of inflammatory factors and aggravates MIRI by regulating the migration, adhesion, and aggregation of dendritic cells to the myocardium [36]. Substantial evidence implies the AGE-RAGE signaling pathway plays a crucial role in MIRI [37]. AGE-RAGE pathway can regulate ventricular arrhythmias, cardiomyocyte apoptosis, and contractile impairment following MIRI, which is a potential therapeutic strategy for ameliorating MIRI [37]. More recent evidence suggests that PI3K/Akt/mTOR signaling pathway protects cardiac injury from MIRI by regulating apoptosis and autophagy in cardiomyocytes [38, 39]. Thus, it can serve as a potential target in the setting of MIRI. According to the above studies, the results of the identified enrichment analyses in the current research play a recognized role in the pathogenesis of MIRI.

We also constructed a PPI network with integrated DEGs and identified the following nine hub genes: Tyrobp, Itgb2, Emr1, Fcer1g, Fcgr1, Fcgr2b, Cybb, Cd68, and Ptprc. Consistent with the bioinformatics results, MIRI mouse model results proved that Emr1, Tyrobp, Itgb2 (synonym: CD18), Fcgr2b, Cybb, and Fcer1g were significantly higher in the MIRI group. The changing trend of CD68, Ptprc, and Fcgr1 was consistent with the bioinformatics results. These results suggest that the results of our bioinformatics analysis may provide a reliable basis for the study of the mechanism of MIRI. The results of pathway and process enrichment analysis showed that hub genes were mainly related to microglia pathogen phagocytosis pathway, neutrophil degranulation, regulation of B-cell-mediated immunity, inflammatory response, response to lipoprotein particle, superoxide anion generation, leukocyte mediated cytotoxicity, and so on. Neutrophils, as the primary responders of MIRI, represent an important component in the protracted inflammatory response and severity related to MIRI [40]. Early research suggested that specific neutrophils degranulation within the ischemic myocardium occurred in MIRI [41]. In addition, superoxide anion generation that occurs during MIRI can further augment MIRI by enhancing the pro-oxidant activity of aconitase [42]. In summary, the enriched pathway and process of these hub genes were, in part, consistent with the pathological process of MIRI shown in previous studies.

It is noticeable that the role of several hub genes in MIRI has been reported in the literature. Palazzo et al. [43] found that CD18^-/- mice subjected to MIRI showed a marked reduction in neutrophils accumulation and myocardial necrosis. The results clearly demonstrate that deficiency of CD18 protects MIRI. Cybb expression levels were found to increase with MIRI [44]. Fortunately, myocardial ischemic postconditioning, melatonin, and RIPerC have protective effects on the myocardium by decreasing Cybb expression in the rat model of MIRI [45]. Fcer1g, adapter protein containing an immunoreceptor tyrosine-based activation motif, plays a pivotal role in the extension of MIRI at least partly through mediating collagen-induced platelet activation [46]. Tyrobp expressed in circulating immune cells is a cell membrane-associated protein, which has been reported to be involved in ischemia-reperfusion injury including the liver, lung, and kidney, but not the heart [47–49]. There are no reports on the involvement of Emr1 and Fcgr2b in MIRI. However, our findings suggest otherwise. Thus, Emr1, Tyrobp, and Fcgr2b are fascinating therapeutic targets, the inhibition of which may protect the myocardium from ischemia-reperfusion injury. Moreover, some important TFs, such as Pura, Srf, Dbp, and Stat4 were uncovered in the present study. Among these TFs, Srf has been reported to mediate MIRI by regulating its target genes (Tagln, Fos, NCX1, Slc8a1, and Egr1) [50]. To sum up, these TFs can further enhance our understanding of MIRI pathogenesis and have the potential to offer novel treatment strategies for MIRI.

Currently, there remains an unmet clinical need for the treatment methods and drugs for MIRI. Natural products especially the extracts of Chinese herbal medicine have been an important source of medicines and drug templates with proven success, which has increasingly drawn widespread attention [51, 52]. In the study, we searched the natural products in Chinese herbal medicine exerting protection for MIRI and discovered 75 natural products. To validate the therapeutic efficacy of natural products against MIRI via targeting the hub genes identified in this research, molecular docking simulations were performed. The results showed that most of the natural products showed a good binding ability to the hub genes, and the binding score values of araloside C, glycyrrhizic acid, ophiopogonin D, polyphyllin I, and punicalagin were the highest. We summarized the mechanism of action of natural products for the treatment of MIRI and found that anti-apoptosis, anti-inflammatory, anti-oxidative stress, regulating energy metabolism, reducing Ca²⁺ overload, and so on were principal mechanisms, especially araloside C, glycyrrhizic acid, ophiopogonin D, polyphyllin I, and punicalagin (Supplementary Table 2). Interestingly, the pathways and processes of the hub genes enrichment analyses were basically consistent with the mechanism of action of natural products. In summary, these natural products (especially araloside C, glycyrrhizic acid, ophiopogonin D, polyphyllin I, and punicalagin) provide a broad application prospect for developing more effective anti-MIRI drugs and will undoubtedly deserve in-depth investigation in the future.

The following are the novelties and multiple strengths of our current research. First, the RRA algorithm, a promising strategy used in the integrated analysis of multiple data sets, was first applied to explore DEGs in MIRI. Second, the functional experiment of the MIRI mice was accomplished to validate the obtained hub genes, which could avoid the bias from the results of pure bioinformatics analysis and thus yield more reliable results. Third, the molecular docking simulations of hub genes and natural products retrieved from the PubChem database were performed to analyze the potential binding effects, which could lead to better prevention and treatment for MIRI. However, our results had several limitations, which should be taken into consideration. On the one hand, our findings are based on currently available gene expression data from the microarray. Some genes, although playing potential roles in pathological processes of MIRI, were ignored in our study because they were not detected by the microarray. On the other hand, further molecular biology experiments are required to validate the function of hub genes in the progression of MIRI and the protective effects of natural products that showed a good ability to bind the hub genes against MIRI.

In conclusion, this study identified several key candidate genes, TFs, and biological pathways and searched for some natural products that exhibited potential therapeutic efficacy on MIRI via targeting the hub genes. Through our research is a preliminary investigation, the findings help us acquire a better understanding of pathogenesis, biomarkers discovery, and therapeutic targets for MIRI, reveal new therapeutic targets for natural products with protective properties against MIRI, and develop more effective anti-MIRI drugs.

Data Availability

GSE58486, GSE108940, and GSE115568 data sets used in this study were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo). The data used to support the findings of this study are included in the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest in the paper.

Authors’ Contributions

Jian Chen conceived and designed the study. Jianru Wang analyzed the data and wrote the manuscript. Xiaohui Li completed the molecular docking simulation. Genhao Fan and Mengmeng Zhang contributed significantly to the functional experiment of the MIRI mice, and Guangcao Peng provided important suggestions.

Acknowledgments

The authors would like to thank the authors of those data sets and the people who manage them. This work was supported by the Key R&D and Promotion Special Project of Henan Province (Scientific Problem Tackling) Project (No. 202102310492), the National Natural Science Foundation of China (Nos. 82004311 and 82104873), the Traditional Chinese Medicine Scientific Research Special Project of Henan Province (Nos. 2019JDZX2013 and 20–21ZY2187), the Shanghai Sailing Program (No. 20YF1445900), the Budgetary Projects of Shanghai University of Traditional Chinese Medicine (No. 2019LK046), the Special Clinical Research Project of Health Profession of Shanghai Municipal Commission of Health (20194Y0081), and Scientific Research Project of Traditional Chinese Medicine of Hongkou District Health (HKQ-ZYY-2020-19).

Supplementary Materials

Supplementary Table 1. The R code used for the RRA analyses. Supplementary Table 2. Summary of natural products in Chinese herbal medicine exerting protection for MIRI. (Supplementary Materials)

References

J. Schmidt-RioValle, M. Abu Ejheisheh, M. J. Membrive-Jiménez et al., “Quality of life after coronary artery bypass surgery: a systematic review and meta-analysis,” International Journal of Environmental Research and Public Health, vol. 17, no. 22, p. 8439, 2020.
View at: Publisher Site | Google Scholar
P. K. Thrombolysis, “In ST-elevation myocardial infarction is not dead,” EuroIntervention, vol. 16, no. 14, pp. 1129-1130, 2021.
View at: Google Scholar
D. M. Yellon and D. J. Hausenloy, “Myocardial reperfusion injury,” New England Journal of Medicine, vol. 357, no. 11, pp. 1121–1135, 2007.
View at: Publisher Site | Google Scholar
C. D. Sánchez-Hernández, L. A. Torres-Alarcón, A. González-Cortés, and A. N. Peón, “Ischemia/reperfusion injury: pathophysiology, current clinical management, and potential preventive approaches,” Mediators Inflamm, vol. 2020, Article ID 8405370, 2020.
View at: Google Scholar
K. Wang, Y. Li, T. Qiang, J. Chen, and X. Wang, “Role of epigenetic regulation in myocardial ischemia/reperfusion injury,” Pharmacological Research, vol. 170, Article ID 105743, 2021.
View at: Publisher Site | Google Scholar
J.-y. Li, S.-q. Liu, R.-q. Yao, Y.-p. Tian, and Y.-m. Yao, “A novel insight into the fate of cardiomyocytes in ischemia-reperfusion injury: from iron metabolism to ferroptosis,” Frontiers in Cell and Developmental Biology, vol. 9, Article ID 799499, 2021.
View at: Publisher Site | Google Scholar
X. Li, N. Ma, J. Xu et al., “Targeting ferroptosis: pathological mechanism and treatment of ischemia-reperfusion injury,” Oxidative Medicine and Cellular Longevity, vol. 2021, Article ID 1587922, 14 pages, 2021.
View at: Publisher Site | Google Scholar
C.-L. Chen, L. Zhang, Z. Jin, T. Kasumov, and Y.-R. Chen, “Mitochondrial redox regulation and myocardial ischemia-reperfusion injury,” American Journal of Physiology–Cell Physiology, vol. 322, no. 1, pp. C12–c23, 2022.
View at: Publisher Site | Google Scholar
M. Cui, Y. Han, J. Yang, G. Li, and C. Yang, “A narrative review of the research status of exosomes in cardiovascular disease,” Annals of Palliative Medicine, vol. 11, no. 1, pp. 363–377, 2022.
View at: Publisher Site | Google Scholar
G. Heusch, “Myocardial ischaemia-reperfusion injury and cardioprotection in perspective,” Nature Reviews Cardiology, vol. 17, no. 12, pp. 773–789, 2020.
View at: Publisher Site | Google Scholar
M. Yousef, A. Kumar, and B. Bakir-Gungor, “Application of biological domain knowledge based feature selection on gene expression data,” Entropy, vol. 23, 1 page, 2020.
View at: Publisher Site | Google Scholar
H. Yu, C. Li, X. Wang et al., “Techniques and strategies for potential protein target discovery and active pharmaceutical molecule screening in a pandemic,” Journal of Proteome Research, vol. 19, no. 11, pp. 4242–4258, 2020.
View at: Publisher Site | Google Scholar
J. Chen, Y.-k. Wang, Y. Gao et al., “Protection against COVID-19 injury by qingfei paidu decoction via anti-viral, anti-inflammatory activity and metabolic programming,” Biomedicine & Pharmacotherapy, vol. 129, Article ID 110281, 2020.
View at: Publisher Site | Google Scholar
M. H. Rahman, S. Peng, X. Hu et al., “A network-based bioinformatics approach to identify molecular biomarkers for type 2 diabetes that are linked to the progression of neurological diseases,” International Journal of Environmental Research and Public Health, vol. 17, no. 3, 2020.
View at: Publisher Site | Google Scholar
M. H. Rahman, H. K. Rana, S. Peng et al., “Bioinformatics and system biology approaches to identify pathophysiological impact of COVID-19 to the progression and severity of neurological diseases,” Computers in Biology and Medicine, vol. 138, Article ID 104859, 2021.
View at: Publisher Site | Google Scholar
M. H. Rahman, H. K. Rana, S. Peng et al., “Bioinformatics and machine learning methodologies to identify the effects of central nervous system disorders on glioblastoma progression,” Briefings in Bioinformatics, vol. 22, no. 5, 2021.
View at: Publisher Site | Google Scholar
H.-Q. Mu, Z.-Q. Liang, Q.-P. Xie et al., “Identification of potential crucial genes associated with the pathogenesis and prognosis of prostate cancer,” Biomarkers in Medicine, vol. 14, no. 5, pp. 353–369, 2020.
View at: Publisher Site | Google Scholar
J. Niu, T. Yan, W. Guo et al., “Identification of potential therapeutic targets and immune cell infiltration characteristics in osteosarcoma using bioinformatics strategy,” Frontiers in Oncology, vol. 10, p. 1628, 2020.
View at: Publisher Site | Google Scholar
Y. Tang, Y. Zhang, and X. Hu, “Identification of potential hub genes related to diagnosis and prognosis of hepatitis B virus-related hepatocellular carcinoma via integrated bioinformatics analysis,” BioMed Research International, vol. 2020, Article ID 4251761, 19 pages, 2020.
View at: Publisher Site | Google Scholar
R. Kolde, S. Laur, P. Adler, and J. Vilo, “Robust rank aggregation for gene list integration and meta-analysis,” Bioinformatics, vol. 28, no. 4, pp. 573–580, 2012.
View at: Publisher Site | Google Scholar
H. Zeng, Y. Fu, L. Shen, and S. Quan, “Integrated analysis of multiple microarrays based on raw data identified novel gene signatures in recurrent implantation failure,” Frontiers in Endocrinology, vol. 13, Article ID 785462, 2022.
View at: Publisher Site | Google Scholar
Y. Zhou, B. Zhou, L. Pache et al., “Metascape provides a biologist-oriented resource for the analysis of systems-level datasets,” Nature Communications, vol. 10, no. 1, p. 1523, 2019.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, D. Lyon et al., “STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–d613, 2019.
View at: Publisher Site | Google Scholar
M. Kohl, S. Wiese, and B. Warscheid, “Cytoscape: software for visualization and analysis of biological networks,” Methods in Molecular Biology, vol. 696, pp. 291–303, 2011.
View at: Publisher Site | Google Scholar
G. D. Bader and C. W. Hogue, “An automated method for finding molecular complexes in large protein interaction networks,” BMC Bioinformatics, vol. 4, no. 1, p. 2, 2003.
View at: Publisher Site | Google Scholar
R. s. Janky, A. Verfaillie, H. Imrichová et al., “iRegulon: from a gene list to a gene regulatory network using large motif and track collections,” PLoS Computational Biology, vol. 10, no. 7, Article ID e1003731, 2014.
View at: Publisher Site | Google Scholar
C. H. Chin, S. H. Chen, H. H. Wu, C. W. Ho, M. T. Ko, and C. Y. L. cytoHubba, “Identifying hub objects and sub-networks from complex interactome,” BMC Systems Biology, vol. 8, 2014.
View at: Publisher Site | Google Scholar
X. Fang, H. Wang, D. Han et al., “Ferroptosis as a target for protection against cardiomyopathy,” Proceedings of the National Academy of Sciences, vol. 116, no. 7, pp. 2672–2680, 2019.
View at: Publisher Site | Google Scholar
J. Zhang, Y. Zhou, Y. Sun et al., “Beneficial effects of oridonin on myocardial ischemia/reperfusion injury: insight gained by metabolomic approaches,” European Journal of Pharmacology, vol. 861, Article ID 172587, 2019.
View at: Publisher Site | Google Scholar
H. C. Fang, B. Q. Wu, Y. L. Hao et al., “KRT1 gene silencing ameliorates myocardial ischemia-reperfusion injury via the activation of the Notch signaling pathway in mouse models,” Journal of Cellular Physiology, vol. 234, no. 4, pp. 3634–3646, 2019.
View at: Publisher Site | Google Scholar
C.-F. Yang, “Clinical manifestations and basic mechanisms of myocardial ischemia/reperfusion injury,” Tzu Chi Medical Journal, vol. 30, no. 4, pp. 209–215, 2018.
View at: Publisher Site | Google Scholar
A. Rout, U. S. Tantry, M. Novakovic, A. Sukhi, and P. A. Gurbel, “Targeted pharmacotherapy for ischemia reperfusion injury in acute myocardial infarction,” Expert Opinion on Pharmacotherapy, vol. 21, no. 15, pp. 1851–1865, 2020.
View at: Publisher Site | Google Scholar
M. Algoet, S. Janssens, U. Himmelreich et al., “Myocardial ischemia-reperfusion injury and the influence of inflammation,” Trends in Cardiovascular Medicine, 2022.
View at: Google Scholar
N. J. Pluijmert, D. E. Atsma, and P. H. A. Quax, “Post-ischemic myocardial inflammatory response: a complex and dynamic process susceptible to immunomodulatory therapies,” Frontiers in Cardiovascular Medicine, vol. 8, Article ID 647785, 2021.
View at: Publisher Site | Google Scholar
L. Arnold, M. Weberbauer, M. Herkel et al., “Endothelial BMP4 promotes leukocyte rolling and adhesion and is elevated in patients after survived out-of-hospital cardiac arrest,” Inflammation, vol. 43, no. 6, pp. 2379–2391, 2020.
View at: Publisher Site | Google Scholar
J. Xue, H. Ge, Z. Lin et al., “The role of dendritic cells regulated by HMGB1/TLR4 signalling pathway in myocardial ischaemia reperfusion injury,” Journal of Cellular and Molecular Medicine, vol. 23, no. 4, pp. 2849–2862, 2019.
View at: Publisher Site | Google Scholar
T.-W. Lee, Y.-H. Kao, Y.-J. Chen, T.-F. Chao, and T.-I. Lee, “Therapeutic potential of vitamin D in AGE/RAGE-related cardiovascular diseases,” Cellular and Molecular Life Sciences, vol. 76, no. 20, pp. 4103–4115, 2019.
View at: Publisher Site | Google Scholar
L. Gong, X. Wang, J. Pan et al., “The co-treatment of rosuvastatin with dapagliflozin synergistically inhibited apoptosis via activating the PI3K/AKt/mTOR signaling pathway in myocardial ischemia/reperfusion injury rats,” Open Medicine, vol. 15, no. 1, pp. 47–57, 2021.
View at: Google Scholar
F. F. Ramli, A. Ali, and N. I. Ibrahim, “Molecular-signaling pathways of ginsenosides Rb in myocardial ischemia-reperfusion injury: a mini review,” International Journal of Medical Sciences, vol. 19, no. 1, pp. 65–73, 2022.
View at: Publisher Site | Google Scholar
Z. V. Schofield, T. M. Woodruff, R. Halai, M. C.-L. Wu, and M. A. Cooper, “Neutrophils-A key component of ischemia-reperfusion injury,” Shock, vol. 40, no. 6, pp. 463–470, 2013.
View at: Publisher Site | Google Scholar
P. Dhôte-Burger, A. Vuilleminot, T. Lecompte et al., “Neutrophil degranulation related to the reperfusion of ischemic human heart during cardiopulmonary bypass,” Journal of Cardiovascular Pharmacology, vol. 25, pp. S124–S129, 1995.
View at: Google Scholar
Y.-R. Chen and J. L. Zweier, “Cardiac mitochondria and reactive oxygen species generation,” Circulation Research, vol. 114, no. 3, pp. 524–537, 2014.
View at: Publisher Site | Google Scholar
A. J. Palazzo, S. P. Jones, W. G. Girod, D. C. Anderson, D. N. Granger, and D. J. Lefer, “Myocardial ischemia-reperfusion injury in CD18- and ICAM-1-deficient mice,” American Journal of Physiology - Heart and Circulatory Physiology, vol. 275, no. 6, pp. H2300–H2307, 1998.
View at: Publisher Site | Google Scholar
G. Aslan, A. Atessahin, and E. Sahna, “The inhibition of apoptosis through myocardial postconditioning by affecting Fas/FasIg signaling through miR139-3p and miR181a-1,” Journal of Cardiac Surgery, vol. 35, no. 3, pp. 564–570, 2020.
View at: Publisher Site | Google Scholar
K. Gul-Kahraman, M. Yilmaz-Bozoglan, and E. Sahna, “Physiological and pharmacological effects of melatonin on remote ischemic perconditioning after myocardial ischemia-reperfusion injury in rats: role of Cybb, Fas, NfκB, Irisin signaling pathway,” Journal of Pineal Research, vol. 67, no. 2, Article ID e12589, 2019.
View at: Google Scholar
N. Takaya, Y. Katoh, K. Iwabuchi et al., “Platelets activated by collagen through the immunoreceptor tyrosine-based activation motif in the Fc receptor γ-chain play a pivotal role in the development of myocardial ischemia-reperfusion injury,” Journal of Molecular and Cellular Cardiology, vol. 39, no. 6, pp. 856–864, 2005.
View at: Publisher Site | Google Scholar
J. H. Spahn, W. Li, A. C. Bribriesco et al., “DAP12 expression in lung macrophages mediates ischemia/reperfusion injury by promoting neutrophil extravasation,” The Journal of Immunology, vol. 194, no. 8, pp. 4039–4048, 2015.
View at: Publisher Site | Google Scholar
T. Nakao, Y. Ono, H. Dai et al., “DNAX activating protein of 12 kDa/triggering receptor expressed on myeloid cells 2 expression by mouse and human liver dendritic cells: functional implications and regulation of liver ischemia-reperfusion injury,” Hepatology, vol. 70, no. 2, pp. 696–710, 2019.
View at: Publisher Site | Google Scholar
D. N. Grigoryev, D. I. Cheranova, D. P. Heruth et al., “Meta-analysis of molecular response of kidney to ischemia reperfusion injury for the identification of new candidate genes,” BMC Nephrology, vol. 14, no. 1, p. 231, 2013.
View at: Publisher Site | Google Scholar
K. M. O’Shea, R. Ananthakrishnan, Q. Li et al., “The formin, DIAPH1, is a key modulator of myocardial ischemia/reperfusion injury,” EBioMedicine, vol. 26, pp. 165–174, 2017.
View at: Publisher Site | Google Scholar
A. A. Izzo, M. Teixeira, S. P. H. Alexander et al., “A practical guide for transparent reporting of research on natural products in the British Journal of Pharmacology: reproducibility of natural product research,” British Journal of Pharmacology, vol. 177, no. 10, pp. 2169–2178, 2020.
View at: Publisher Site | Google Scholar
C. Zhao, S. Li, J. Zhang et al., “Current state and future perspective of cardiovascular medicines derived from natural products,” Pharmacology & Therapeutics, vol. 216, Article ID 107698, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Jianru Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

501

Downloads

712

Citations