Comprehensive Analysis of N6-Methyladenosine RNA Methylation Regulators Expression Identify Distinct Molecular Subtypes of Myocardial Infarction

Background: Myocardial infarction (MI) is one of the leading threats to human health. N6-methyladenosine (m6A) modification, as a pivotal regulator of messenger RNA stability, protein expression, and cellular processes, exhibits important roles in the development of cardiac remodeling and cardiomyocyte contractile function. Methods: The expression levels of m6A regulators were analyzed using the GSE5406 database. We analyzed genome-wide association study data and single-cell sequencing data to confirm the functional importance of m6A regulators in MI. Three molecular subtypes with different clinical characteristics were established to tailor treatment strategies for patients with MI. We applied pathway analysis and differentially expressed gene (DEG) analysis to study the changes in gene expression and identified four common DEGs. Furthermore, we constructed the protein–protein interaction network and confirmed several hub genes in three clusters of MI. To lucubrate the potential functions, we performed a ClueGO analysis of these hub networks. Results: In this study, we identified that the levels of FTO, YTHDF3, ZC3H13, and WTAP were dramatically differently expressed in MI tissues compared with controls. Bioinformatics analysis showed that DEGs in MI were significantly related to modulating calcium signaling and chemokine signaling, and m6A regulators were related to regulating glucose measurement and elevated blood glucose levels. Furthermore, genome-wide association study data analysis showed that WTAP single-nucleotide polymorphism was significantly related to the progression of MI. In addition, single-cell sequencing found that WTAP is widely expressed in the heart tissues. Moreover, we conducted consensus clustering for MI in view of the dysregulated m6A regulators’ expression in MI. According to the expression levels, we found MI patients could be clustered into three subtypes. Pathway analysis showed the DEGs among different clusters in MI were assigned to HIF-1, IL-17, MAPK, PI3K-Akt signaling pathways, etc. The module analysis detected several genes, including BAG2, BAG3, MMP2, etc. We also found that MI-related network was significantly related to positive and negative regulation of angiogenesis and response to heat. The hub networks in MI clusters were significantly related to antigen processing and ubiquitin-mediated proteolysis, RNA splicing, and stability, indicating that these processes may contribute to the development of MI. Conclusion: Collectively, our study could provide more information for understanding the roles of m6A in MI, which may provide a novel insight into identifying biomarkers for MI treatment and diagnosis.


INTRODUCTION
Myocardial infarction (MI) is referred to as a heart attack event term in which myocardial cells die because of imbalances between oxygen supply and demand. The general definition of MI consists of five subtypes, among which type 1 and type 2 MI are the widely occurred sorts and also exhibit a tight relationship to clinicians (Gao et al., 2016). Type 1 MI is induced by acute atherosclerotic thrombosis, whereas type 2 MI results from the imbalanced supply and demand of myocardial oxygen without acute atherosclerotic thrombosis (Lu et al., 2015;Sandoval and Jaffe, 2019). The most common cause of MI is blood flow to a part of the heart is reduced or paused, leading to myocardial necrosis. Also, the result of a blood coagulum in the coronary artery is responsible for supplying the heart muscle area. The risk of MI usually increases with age, particularly at the age of more than 65 years (Lu et al., 2015).
Previous research has shown that epigenetic regulation plays importantly in the regulation of cardiovascular repair functions. Typically, as the ubiquitous and abundant transcription modification of messenger RNA (mRNA) and long non-coding RNA in the genome of eukaryotes, N6-methyladenosine (m6A) primarily occurs in the 3 -untranslated regions and nearby the stop codons of mRNAs (Meyer et al., 2012;Wang et al., 2014;Li et al., 2017;Liu et al., 2017;Zhang et al., 2017). Acting as a reversible modification, m6A is methylated by m6A methyltransferases (writers), demethylated by m6A demethylases (erasers), and recognized by m6A binding proteins (readers), which participates in various biological processes. Dysregulated m6A is gradually thought to be the etiopathogenesis of certain diseases, such as carcinomas and cardiovascular disease (CVD). Growing evidence indicated that the continuous dynamic modulation of m6A exerts an impact on specific genes' expression and some diseases' physiological and pathological processes, including MI, ischemic heart failure (HF), myocardial hypertrophy, and cardiomyogenesis (Mathiyalagan et al., 2019). For instance, WTAP promotes MI via modulating m6A modification of ATF4 mRNA . Dorn et al. (2019) showed that N6 adenosine methylation regulated by METTL3 is crucial for hypertrophy's pathological process in vivo and in vitro. Loss of METTL3 enhanced heart regeneration and repair after myocardial injury . Nonetheless, the roles of m6A regulators in MI remained largely unclear.
So far, the function of m6A in physiological and biological processes has been largely studied. Nevertheless, the research toward m6A is limited under pathological conditions in MI. Here, we systematically investigated the role of m6A epigenetic regulation in MI. Furthermore, we applied unsupervised consensus clustering on gene expression in MI based on m6A expression. In addition, we identified the potential roles of differentially expressed genes (DEGs) in different clusters based on m6A expression by performing Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and protein-protein interaction (PPI) network analysis. As far as we know, our literature, for the first time, comprehensively analyzed m6A's roles in MI, which may provide novel clues to identify biomarkers for MI treatment and diagnosis.

Data Collection
GSE5406 (Hannenhalli et al., 2006) included gene expression profile of human left ventricle tissue from 16 non-failing myocardium samples, 108 ischemic myocardium samples, and 86 idiopathic myocardium samples. Genome-wide association study (GWAS) data were downloaded from the GeneAtlas database. 1 The single-cell data were from The Single Cell Type Atlas.

Enrichment Analysis
We carried out functional enrichment analysis, including Gene Ontology (GO) and KEGG, using the DAVID system. 2 P-value < 0.05 meant significant difference.

Construction of Protein-Protein Interaction Network
The STRING database was utilized to construct the PPI network. The PPI network was displayed by Gephi software. We chose the hub genes based on the extent of genes' connectivity.

Unsupervised Consensus Clustering Analysis
We carried out unsupervised consensus analysis utilizing the ConsensusClusterPlus (Wilkerson and Hayes, 2010) R package. Taken briefly, the consistent matrix plots were shown in the light of the k-value. Also, empirical cumulative distribution function plots exhibited the uniform distributions for each k. What is more, the item tracking plot demonstrated the consistent clustering of items (columns) at each k (rows) to determine the clustering stability. The cluster-consensus plot illustrated the cluster-consensus value at different k-values. High cluster-consensus value meant clustering with low stability. Itemconsensus plot is the mean consensus value deriving from an item and a consensus cluster's members. An item indicated several item-consensus values with different ks.

Differentially Expressed Gene Analysis
The DEGs were determined by the Limma R package between different clusters in AMI, observing the cutoff value of | log2 (fold change) | ≥ 1 and P-value < 0.05. Ggplot R package was applied to draw the volcano plot, and pheatmap R package was used to plot the heatmap of DEGs.

Potential Hub Gene Identification
Each protein node's degree was evaluated by the CytoHubba plugin in Cytoscape software 3 (Chin et al., 2014;Ma et al., 2021). Nodes with a higher degree of connectivity are normally more important for keeping the entire network's stability. The 10 proteins with closet connection nodes here were the potential hub genes.

N6-Methyladenosine Regulators Expression Was Dysregulated in Myocardial Infarction
For assessing m6A regulators' biological function in MI's occurrence and development, we explored 15 m6A regulators' expression profiles systematically between MI and non-failing controls using GSE5406. Figure 1A illustrates m6A modulators' expression levels in MI and non-failing controls separately. Among them, we observed that FTO, YTHDF3, and ZC3H13 expression levels were largely lower in MI tissues than in non-failing controls (p < 0.01); however, WTAP expression levels were markedly higher in MI tissues relative to nonfailing controls (p < 0.05). In addition, there was no obvious difference between the non-failing controls and MI tissues with regard to RBM15, YTHDC1, HHRNPA2B1, IGF2BP2, METTL3, RBMX, YTHDF1, YTHDF2, HNRNPC, IGF2BP3, RBM15B, and YTHDC2 expression levels. To evaluate whether m6A modulators functioned crucially in the development of MI, we analyzed the correlation among them in non-failing controls and MI samples separately. The results demonstrated that the correlation among m6A regulators has a significant change between control and MI groups, implying that m6A modulators functioned crucially in the development of MI ( Figure 1B). The dysregulation of FTO, YTHDF3, ZC3H13, and WTAP in MI is presented in Figure 1C.

Functional Annotation of N6-Methyladenosine Regulators in Myocardial Infarction
To evaluate the potential functions of m6A regulators in MI, we performed bioinformatics analysis. As present in Figure 2A, the mountain map shows that the DEGs in MI were significantly enriched in regulation of locomotion, cation transport, response to organonitrogen compound, signaling receptor activity, cellular amide metabolic process, enzyme regulator activity, transporter activity, chemical homeostasis, oxoacid metabolic process, cytokine-cytokine receptor interaction, chemokine signaling, calcium signaling, phospholipase D signaling, thermogenesis, cAMP signaling, oxytocin signaling, lysosome, cGMP-PKG signaling, and adrenergic signaling in cardiomyocytes. Of note, Gene Set Enrichment Analysis further demonstrated that DEGs in MI were related to transport activity and calcium signaling ( Figure 2B). Moreover, by analyzing the relationship among the top 5 enriched signalings, we observed these pathways have strong crosstalk ( Figure 2C). Calcium signaling and chemokine signaling play a key role in MI ( Figure 2D).
We also predicted that the potential signaling related to m6A regulators and observed FTO and IGFBP2 are mainly related to glucose measurement, elevated blood glucose level, and chronic schizophrenia. HNRNPA2B1 and HNRNPC are mainly related to Astler-Coller B1 rectal carcinoma. IGF2BP3 was related to soft tissue neoplasms and peritoneal neoplasms ( Figure 2E).

Confirmation of the Functional Importance of N6-Methyladenosine Regulators in Myocardial Infarction
Next, we analyzed the GWAS data and single-cell sequencing data to confirm the functional importance of m6A regulators in MI. As present in Figure 3A, the Q-Q plot shows that the GWAS data can identify significantly related single-nucleotide polymorphism (SNP) sites. Through precise positioning of the GWAS data, causal SNPs are mainly distributed in the enriched area ( Figure 3B). The SNP annotation of the locus found that the m6A regulatory gene WTAP was located in the causal region on chromosome 6, suggesting it may be an important pathogenic gene of MI ( Figure 3C). Furthermore, we analyzed the single-cell sequencing of the heart tissues and found that WTAP is widely expressed in the heart tissues, whose expression in endothelial cells is the highest (Figure 3D). By combining the analysis mentioned earlier that WTAP was significantly lower in the MI group, this suggests that maintaining the normal expression of WTAP may be a prevention or improvement means for the treatment of MI.

Consensus Clustering Analysis for
Myocardial Infarction Based on the Expression of N6-Methyladenosine Regulators m6A RNA methylation exerted an important effect on modulating the mRNA stability (Wang et al., 2014), alternative splicing, and RNA structure (Liu et al., 2017), thus affecting the RNA expression of the human genome. Thus, we performed consensus clustering for MI in view of the dysregulated m6A modulators in MI.
We conducted the consensus clustering utilizing the Consensus Cluster Plus R package. The cumulative distribution function presented the lowest range ability at consensus index 0.2-0.6 with k = 3 ( Figure 4A). The delta area scores of 2.5 were the highest at k = 3 ( Figure 4B). Totally, 499 MI patients were clustered into three subtypes, including cluster1 (n = 125), cluster2 (n = 53), and cluster 3 (n = 16) on the basis of m6A regulators expression levels ( Figure 4C).

Identifying the Differentially Expressed Genes Between Different Clusters in Myocardial Infarction
The DEGs among cluster 1, cluster 2, and cluster 3 MI samples were identified based on m6A regulators expression in depth. Four hundred thirty-seven DEGs were identified in total, including 237 DEGs with upregulation and 200 DEGs with downregulation between MI and non-failing controls (Figures 5A,B). The top 10 overexpressed and suppressed genes between MI and non-failing controls are listed in Table 1.
Totally, 456 DEGs, comprising 420 DEGs with upregulation and 36 DEGs with downregulation, were screened after comparing cluster 1 with cluster 2 samples (Figures 5C,D). The top 10 overexpressed and suppressed between cluster 1 and cluster 2  are listed in Table 2. One thousand three hundred ninety-five DEGs in total, consisting of 884 DEGs with upregulation and 511 DEGs with downregulation, were obtained after comparing cluster 1 with cluster 3 samples (Figures 5E,F). The top 10 overexpressed and suppressed genes between cluster 1 and cluster 3 are listed in Table 3. One thousand six hundred seventy-eight DEGs were acquired, including 886 DEGs with upregulation and 792 DEGs with downregulation after comparing cluster 2 with cluster 3 samples (Figures 5G,H). The top 10 overexpressed and suppressed genes between cluster 2 and cluster 3 are listed in Table 4.

Bioinformatics Analysis of the Differentially Expressed Genes Between Different Clusters in Myocardial Infarction
We further performed KEGG analysis of DEGs between different clusters to evaluate their potential functions in MI. As present in Figure 4, the overexpressed genes in cluster 1 compared with cluster 2 were significantly related to fatty acid biosynthesis, RNA transport, ubiquitinmediated proteolysis, protein processing ( Figure 6A). After comparison with cluster 2, the downregulated genes in cluster 1 were significantly related to toxoplasmosis, cell adhesion molecules, antigen processing and presentation, and alcoholism ( Figure 6B). In comparison with cluster 3, the overexpressed genes in cluster 1 were significantly related to antigen processing and presentation, autophagy, and signaling pathways, including AGE-RAGE, AMPK, PI3K-Akt, TGF-beta, and cGMP-PKG signaling ( Figure 6C).  After comparison with cluster 3, the overexpressed genes in cluster 2 were significantly related to antigen processing and presentation and signaling pathways, including AGE-RAGE, AMPK, FoxO, HIF-1, PI3K-Akt, and TGF-beta signaling ( Figure 6E). The downregulated genes in cluster 1 compared with cluster 3 and cluster 2 compared with cluster 3 exhibited a significant relation to oxidative phosphorylation, spliceosome, and multiple neurodegenerative diseases, such as Alzheimer's disease, Huntington's disease, and Parkinson's disease (Figures 6D-F).

Screening of Key Genes in Different Clusters of Myocardial Infarction
We also investigated changes in gene expression among different clusters in MI based on the expression change of the GSE5406 dataset. As present in Figure 6, we observed that approximately 45.7% of DEGs were differently expressed in more than two different clusters ( Figure 7A). Venn diagram showed that there were four common DEGs among all clusters in MI, including BAG2, CD74, GOLGA8N, and PDK4 ( Figure 7B). The results showed that BAG2 and PDK4 have a similar expression pattern in MI, which were suppressed in all clusters of MI compared with non-failing controls, suppressed in cluster 2 in comparison with cluster 1 samples, and suppressed in cluster 3 after comparison with cluster 2 and cluster 1 samples (Figures 7C,F). Meanwhile, we found that CD74 was enhanced in clusters 1 and 2 after comparison with non-failing controls, enhanced in cluster 2 in comparison with cluster 1 samples, but suppressed in cluster 3 relative to cluster 1 and 2 samples ( Figure 7D). Finally, the results showed that GOLGA8N was enhanced in clusters 1 and 3 in comparison with non-failing controls but suppressed in cluster 2 in comparison with cluster 1 samples and enhanced in cluster 3 in comparison with cluster 1 and 2 samples ( Figure 7E).

Screening of Hub Networks Between Different Clusters of Myocardial Infarction
STRING tools were utilized to predict the PPI of the DEGs. Then, we used Cytoscape MCODE to further screen the core PPI network, and the results revealed four hub networks related to different clusters of MI. In total, 37 nodes and 259 edges were included in the PPI network consisting of DEGs in MI compared with normal samples (Figure 8A), 16 nodes and 120 edges were included in the PPI network consisting with DEGs in cluster 1 MI compared with cluster 2 MI samples (Figure 8B), 55 nodes and 739 edges were included in the PPI network consisting with DEGs in cluster 1 MI compared with cluster 3 MI samples (Figure 8C), and 63 nodes and 982 edges were included in the PPI network consisting with DEGs in cluster 2 MI compared with cluster 3 MI samples (Figure 8D). The module analysis filtered out several genes, including BAG2, BAG3, MMP2, FGF1, HSPB1, CXL12, PJA2, UBE2M, RNF14, RNF16, KLHL2, LTN1, RBM17, LSM2, TRIM37, and UBA3.

Functional Annotation of Hub Networks Between Different Clusters of Myocardial Infarction
To understand the potential functions involved in the hub networks and the connection between hub genes and biological functions, we performed a ClueGO analysis of these hub networks in MI. As present in Figure 8, the hub PPI network consisting of DEGs in MI compared with normal samples was significantly related to positive and negative modulation of angiogenesis and positive modulation of ATPase activity, response to heat, and extracellular matrix organization ( Figure 9A). Of note, we observed multiple HSP family genes, such as HSPH1, HSPA2, HSPA6, HSP90AA1, and HSPD1, were related to MI progression by modulating response to heat. The hub PPI network consisting of DEGs in cluster 1 MI compared with cluster 2 MI samples was significantly related to antigen processing and ubiquitin-mediated proteolysis ( Figure 9B). Also, the hub network consisting of DEGs in cluster 1 MI compared with cluster 3 MI samples was involved in regulating multiple biological processes related to protein ubiquitination and RNA splicing ( Figure 9C). Finally, the hub network consisting of DEGs in cluster 2 MI compared with cluster 3 MI samples was involved in regulating antigen processing and RNA splicing and stabilization ( Figure 9D). Interestingly, we found that C3 genes were significantly related to RNA splicing and stability, consistent with the crucial roles of m6A in RNA splicing and stability. These results indicated that the abnormal regulation of RNAs splicing and stability might contribute to the development of MI.

DISCUSSION
MI is still the leading threat to human health. m6A is considered one of the most common and abundant RNA methylation modifications in eukaryotes (Wang et al., 2014;Mathiyalagan et al., 2019;Wen et al., 2019). m6A is regulated by m6A methyltransferases and demethylases and controls the fate of target mRNA by affecting splicing, translation, and decay. Recently, some researches suggested that m6A modification exhibits importantly in the development of cardiac remodeling and cardiomyocyte contractile function (Dorn et al., 2019;Huang et al., 2019;Han et al., 2021;Wang et al., 2021). For example, ablation of METTL3 weakened MI-caused myocardial fibrosis by impeding the activation of cardiac fibroblasts (Dorn et al., 2019). WTAP promotes MI via modulating ATF4 . ALKBH5 mediated the modulation of heart regeneration via demethylating YTHDF1 mRNA (Han et al., 2021). However, there was lacking the comprehensive analysis of the roles of m6A in MI. In this study, we observed that the levels of FTO, YTHDF3, ZC3H13, and WTAP were dramatically differently expressed in MI tissues compared with non-failing controls, further demonstrating the crucial roles of m6A regulators in MI, which were consistent with previous reports. For example, FTO gene polymorphisms were related to HDL cholesterol concentration and high risk of CVD (Franczak et al., 2018) and could predict the incidence of CVD (Aijala et al., 2015). We also predicted the potential signaling related to m6A regulators and observed that FTO and IGFBP2 are mainly related to glucose measurement, elevated blood glucose level, and chronic schizophrenia. Of note, multiple previous reports had indicated the regulation among FTO and glucose metabolism (Khan et al., 2018;Mizuno, 2018). For example, GWASs showed that FTO mutation is related to impaired fasting glucose (Khan et al., 2018). In the liver, FTO regulated glucose and lipid metabolism, which was modulated by metabolic signals (Mizuno, 2018). HNRNPA2B1 and HNRNPC are mainly related to Astler-Coller B1 rectal carcinoma. IGF2BP3 was related to soft tissue neoplasms and peritoneal neoplasms. Next, we analyzed GWAS data and single-cell sequencing data to confirm the functional importance of m6A regulators in MI. The SNP annotation of the locus found that the m6A regulatory gene WTAP was located in the causal region on chromosome 6, suggesting it may be an important pathogenic gene of MI. Furthermore, we analyzed the single-cell sequencing of the heart tissues and found that WTAP is widely expressed in the heart tissues, whose expression in endothelial cells is the highest. By combining the analysis mentioned earlier that WTAP was significantly lower in the MI group, this suggests that maintaining the normal expression of WTAP may be a prevention or improvement means for the treatment of MI. Very interestingly, our findings were consistent with a recent report. Wang et al. (2021) reported that WTAP promotes MI via modulating ATF4. This study, for the first time, reveals novel clues to understand the mechanism of MI based on m6A expression.
The general definition of MI usually consists of five subtypes. In this study, we conducted consensus clustering for MI in view of the dysregulated m6A RNA methylation regulators expression in MI. According to m6A regulators' expression levels, we clustered 499 MI patients into three subtypes. Our findings demonstrated that FTO, YTHDF1, IGF2BP3, RBM15B, YTHDF3, IGF2BP2, and RBM15 were mostly upregulated in cluster 3. WTAP, HNRNPC, YTHDF2, YTHDC2, HNRNPA2B1, METTL3, ZC3H13, and YTHDC1 were most significantly upregulated in cluster 1. Also, WTAP and HNRNPC were significantly downregulated in cluster 3. Furthermore, our results showed that the DEGs in MI were significantly enriched in multiple signalings, such as the chemokine signaling pathway, calcium signaling pathway, and chemokine signaling. Interestingly, Gene Set Enrichment Analysis further demonstrated that DEGs in MI were related to transport activity and calcium signaling. Moreover, by analyzing the relationship among the top 5 enriched signaling, calcium signaling and chemokine signaling play a key role in MI. Calcium (Ca 2+ ), as a second messenger, played a key role in regulating cell proliferation, apoptosis, and survival (Sukumaran et al., 2021). Previous studies had demonstrated that calcium signaling served as a key regulator of MI (Vassalle and Lin, 2004;Li et al., 2021). Under physiological conditions, calcium signaling can effectively modulate the activity of vasodilation/contraction in vascular smooth muscle cells . However, under pathological conditions, Ca2 + overload induces apoptosis in cardiomyocytes, thus leading to MI (Vassalle and Lin, 2004;Jiao et al., 2019).
Over the past decades, multiple chemokines and their receptors were found to regulate MI (Altara et al., 2016;Wang et al., 2019;Liang et al., 2021). For example, CXCR3 and its ligands were reported to be valid biomarkers for HF (Altara et al., 2016). CXCL16 was reported to modulate the inflammatory responses in MI . CXCR4 blockade was found to induce tissue repair after MI by modulating immune-regulatory function (Wang et al., 2019).
Also, bioinformatics analysis indicated that the DEGs between these clusters in MI were related to modulating multiple signalings, such as PI3K-Akt, MAPK, and cGMP-PKG signalings. It is well-known that the PI3K-Akt signaling pathway is a survival-associated signal transduction pathway, protecting the myocardium from ischemic damage (Okumura et al., 2004;Quan et al., 2014;Cheng et al., 2018). Additionally, previous researches have shown that this signaling pathway is motivated when the estrogen receptor β is chronically upregulated. ERK1/2, JNK, and p38 were three subfamilies of MAPK (Zhang and Liu, 2002). The MAPK signaling pathway was reported to promote NF-κB, triggering additional inflammatory cytokines and causing extra damage in myocardial tissue (Verma et al., 2019). The ERK signaling pathway is currently the widely studied one, and it exhibited an association with a variety of biological processes' regulation, containing cell survival, growth and death, and inflammationassociated immune responses (Ren et al., 2019). cGMP-PKG signaling was considered to be a therapeutic target for myocardial ischemia-reperfusion injury (Ren et al., 2019). As previously described, some studies have shown that cGMP-PKG signaling exerted an effect on the function of the myocardial endoplasmic reticulum and decreased the level of endoplasmic reticulum stress under stress (Schlossmann and Desch, 2011). Moreover, the downregulated genes in cluster 1 compared with cluster 3 and cluster 2 compared with cluster 3 exhibited a significant relation to multiple neurodegenerative diseases. Of note, previous studies revealed abnormal regulation of m6A is related to neurodegenerative disease. For example, Han et al. (2020) reported that the m6A methylation enhanced the development of Alzheimer's disease. In addition, Du et al. (2021) revealed that the expressions of m6A regulators correlate with neurodegenerative pathways. These reports, together with our findings, revealed the potential crucial roles of m6A in MI.
In our study, we also investigated changes in gene expression among different clusters in MI based on the expression change of the GSE5406 dataset. As present in Figure 7, we observed that approximately 45.7% of DEGs in MI compared with non-failing samples were differently expressed in more than two different clusters. Venn diagram showed four common DEGs among all clusters in MI, including BAG2, CD74, GOLGA8N, and PDK4. Also, those genes were probably regarded as potential indicators for the prediction and diagnosis of MI. For example, BAG3 belongs to the BAG protein family acting as a chaperone molecular via physical interaction with Hsp70, HSPBs. It was shown that mutated BAG3 caused DCM, leading to systolic dysfunction, HF, and myofibrillar myopathy (Sturner and Behl, 2017;Diofano et al., 2020). BAG2 was also a member of the human BAG protein family and is expressed in brown adipose, lung, heart, and other tissues. It is shown that BAG2 presented a highly similar sequence and domain to BAG3, indicating BAG2 could compensate for the loss function upon the absence of BAG3 (Sturner and Behl, 2017). PDH kinases, including PDK1-4 mediated the phosphorylation and inactivation of PDH (Vary and Randle, 1985). It proves that the selectively overexpressed PDK4 in the heart gave rise to an obvious alteration in energy metabolism, including the increase in the utilization of fatty acids and the decrease in carbohydrate consumption with the increased pyruvate concentration. These results suggested that overexpressed PDK4 is sufficient to dramatically bring about the change of substrate utilization in the heart (Zhao et al., 2008). CD74 was reported to mediate the effect of macrophage migration inhibitory factor on inflammation and cell proliferation. Migration inhibitory factor-CD74 agonism supplied a potential treatment for acute myocardial ischemia by enhancing AMPK activation (Miller et al., 2008). These results indicated that targeting these genes may provide novel therapy strategies for MI.
Utilizing the STRING database, we constructed the PPI network. The module analysis filtered out several genes, including BAG2, BAG3, MMP2, FGF1, HSPB1, CXL12, PJA2, UBE2M, RNF14, RNF16, KLHL2, LTN1, RBM17, LSM2, TRIM37, UBA3, etc. Some of these genes in acute MI (AMI) have been reported, indicating that our conclusions were consistent with previous integrated bioinformatics analysis data. For instance, Wei Gong et al. (2018) revealed that trimetazidine could prevent cardiac rupture in AMI mice by inhibiting matrix metalloproteinase-2 (Mmp2) and Mmp9 expression levels, indicating that the MMP family was probably related to cardiac remodeling after AMI. Fibroblast growth factor 1 (FGF1), also named acidic FGF, took part in several physiological processes, e.g., development, morphogenesis, and others (Engel et al., 2006;Beenken and Mohammadi, 2009). Overexpression of FGF1 or intramyocardial infusion of FGF1 mimic could exert cardioprotective effects on MI mice models (Htun et al., 1998). Heat-shock proteins (HSPs) are largely expressed in various cells of the cardiovascular system, including endothelial cells, cardiomyocytes, monocytes, and platelets. Kraemer et al. (2019) for the first time, showed the upregulation and phosphorylation of HSP27 (HSPB1) of human platelets during MI on a cellular level ex vivo, showing a typical intracellular translocation pattern. Therefore, HSP27 (HSPB1) phenotype in platelets was a measurable stress response indicator in MI and other acute ischemic events [38,39]. After MI, the increased expression of CXCL12 was observed in the infarct area and served as a homing signal for endothelial progenitor cells (Jiao et al., 2019). To understand the potential functions involved in the hub networks and the connection between hub genes and biological functions, we performed ClueGO analysis of these hub networks in MIs. We found that MI-related network was significantly related to positive and negative regulation of angiogenesis and response to heat.
The hub networks in MI clusters were significantly related to antigen processing and ubiquitin-mediated proteolysis, RNA splicing, and stability. These results indicated that the abnormal regulation of RNAs splicing and stability might contribute to the development of MI.

CONCLUSION
In summary, we observed levels of FTO, YTHDF3, ZC3H13, and WTAP were dramatically differently expressed in MI tissues compared with non-failing controls, further demonstrating the crucial roles of m6A regulators in MI. Bioinformatics analysis showed that m6A regulators were related to regulating glucose measurement and elevated blood glucose levels. Furthermore, GWAS data analysis showed that WTAP SNP was significantly related to the progression of MI. In addition, the single-cell sequencing of the heart tissues found that WTAP is widely expressed in the heart tissues. By combining the analysis mentioned earlier that WTAP was significantly lower in the MI group, this suggests that maintaining the normal expression of WTAP may be a prevention or improvement means for the treatment of MI. Furthermore, three subtypes with different clinical characteristics were constructed based on m6A regulators' expression profiles, indicating that they were potential therapeutic strategies for MI patients. This study could provide more information for us to understand the roles of m6A in MI, which may provide novel clues to identify biomarkers for clinical treatment and diagnosis.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
BnH and L-HS conceived and designed the project and are responsible for the overall content. YC, XZ, CG, H-WN, FL, JX, and ZW analyzed and interpreted the data. XS, CG, BrH, and BnH prepared the manuscript. XS, YC, XZ, YL, JX, XW, and ZC contributed to revising the manuscript. All authors contributed to discussed the results, critically reviewed the manuscript, read and approved the final manuscript.