Uncovering Potential lncRNAs and mRNAs in the Progression From Acute Myocardial Infarction to Myocardial Fibrosis to Heart Failure

Background: Morbidity and mortality of heart failure (HF) post-myocardial infarction (MI) remain elevated. The aim of this study was to find potential long non-coding RNAs (lncRNAs) and mRNAs in the progression from acute myocardial infarction (AMI) to myocardial fibrosis (MF) to HF. Methods: Firstly, blood samples from AMI, MF, and HF patients were used for RNA sequencing. Secondly, differentially expressed lncRNAs and mRNAs were obtained in MF vs. AMI and HF vs. MF, followed by functional analysis of shared differentially expressed mRNAs between two groups. Thirdly, interaction networks of lncRNA-nearby targeted mRNA and lncRNA-co-expressed mRNA were constructed in MF vs. AMI and HF vs. MF. Finally, expression validation and diagnostic capability analysis of selected lncRNAs and mRNAs were performed. Results: Several lncRNA-co-expressed/nearby targeted mRNA pairs including AC005392.3/AC007278.2-IL18R1, AL356356.1/AL137145.2-PFKFB3, and MKNK1-AS1/LINC01127-IL1R2 were identified. Several signaling pathways including TNF and cytokine–cytokine receptor interaction, fructose and mannose metabolism and HIF-1, hematopoietic cell lineage and fluid shear stress, and atherosclerosis and estrogen were selected. IL1R2, IRAK3, LRG1, and PLAC4 had a potential diagnostic value for both AMI and HF. Conclusion: Identified AC005392.3/AC007278.2-IL18R1, AL356356.1/AL137145.2-PFKFB3, and MKNK1-AS1/LINC01127-IL1R2 lncRNA-co-expressed/nearby targeted mRNA pairs may play crucial roles in the development of AMI, MF, and HF.


INTRODUCTION
Heart failure (HF), the terminal stage of a wide range of cardiovascular diseases, leads to the decompensation of the heart's ability to contract or relax (1). In most clinical cases, HF is caused by myocardial infarction (MI). Following MI, a lot of cardiac cells will die in response to ischemia (2). As the adult heart has a very limited capacity to regenerate after injury, the lost cardiac cells will be replaced by fibrotic scars. The process is accompanied by remodeling of the surrounding myocardium and eventually results in impaired cardiac function. HF is significantly associated with cardiac hypertrophy, valvular disease, dilated cardiomyopathy, hypertrophic cardiomyopathy, age, obesity, diabetes, and hypertension (3,4). Despite effective treatments targeted at improving HF, the 5-year mortality rate of the disease varies from 45 to 60% (5). Therefore, in-depth understanding of pathological mechanisms of HF is very important to develop effective therapy and assessment of patient prognosis. Long non-coding RNAs (lncRNAs) regulate mRNA expression at the transcriptional and post-transcriptional levels. It has been shown that lncRNAs are emerging as essential players in heart development, atherosclerosis, and HF (6). In addition, the major component of HF disease management is diagnosis and monitoring of disease progression. In this study, we tried to find potential differentially expressed lncRNAs and mRNAs in the blood of patients with acute myocardial infarction (AMI), myocardial fibrosis (MF), and HF. We supposed that these potential differentially expressed lncRNAs and interacted mRNAs may be involved in the disease progression from AMI, MF and HF.

Subjects
In this study, three patients with AMI, three patients with MF, and three patients with HF were enrolled. The inclusion criteria of AMI patients were as follows: (1) time of chest pain or distress > 30 min within 24 h, and the level of the myocardial enzyme of creatine kinase (CK)-MB and cardiac troponin T (cTnT) was higher than the normal range; (2) patients had their first episode; (3) patients received no medical or surgical treatment prior to admission; (4) patients had blood samples before admission, at discharge, and 6 months after MI; and (5) patients had complete clinical data, including gender, age, height, and weight. The exclusion criteria of AMI patients were as follows: (1) patients with myocarditis and other diseases caused by chest pain or distress; (2) patients with a history of renal failure, advanced liver disease, malignant tumors, and other inflammatory diseases such as psoriasis and rheumatoid arthritis; (3) recurrent patients; (4) patients with incomplete clinical data; and (5) patients with missing blood samples before hospitalization, at discharge, and 6 months after MI. Those patients diagnosed with MF (1 year after MI) were enrolled. Depending on the diagnostic criteria, those patients with HF caused by MI were also enrolled.

Blood Samples
The blood of AMI, MF, and HF patients was collected. Ethical approval was obtained from the ethics committee of the Shijiazhuang People's Hospital. In addition, informed written consent was obtained from all individuals.

RNA Isolation and Sequencing
The total RNA was extracted from the blood sample using the RNAliquid Reagent according to the manufacturer's protocols. Agilent 2100 was utilized to test the RNA integrity number (RIN) (7). RNA conditions to be met for sequencing were as follows: (1) the total amount (required for single database construction) is 5 µg; (2) the concentration was ≥200 ng/µl; and (3) the value of OD260/280 was between 1.8 and 2.2. After total RNA DNase I digestion, removal of rRNA, RNA interruption, synthesis of reverse transcriptional one and two strands, end repair, cDNA with an "A" at the 3 ′ end, connection of cDNA 5 ′ adapter, digestion of two strands of cDNA, PCR reaction, and product recovery and library quality inspection, the Illumina Hiseq x-ten platform (PE150 strategy) (8) was used to perform RNA sequencing for lncRNA. After RNA fragment selection, 3' connector connection, reverse primer annealing, 5' connector connection, the synthesis of a strand of cDNA, PCR amplification, library fragment selection, library quantification, and pooling annularization, BGISEQ-500 platform (SE50 strategy) (9) was used to perform RNA sequencing for mRNA. The Combat function for the SVA package was used to eliminate batch effect (10). The Fastx-Toolkit was used to trim 5' and 3' segments of reads to remove bases with mass <20 and delete reads with N >10%. For lncRNA analysis, the cleaned gene sequencing reads were aligned to the reference genome (GRCh38) for matching via HISAT2 (11) (https:// ccb.jhu.edu/software/hisat2/index.shtml). Stringtie (12) (http:// www.ccb.jhu.edu/software/stringtie/) was utilized to quantify the expression levels of lncRNA and mRNA. For mRNA analysis, the Rfam (13) was used for annotation analysis on measured small RNA. The mature mRNA and mRNA precursor sequences were downloaded from miRBase (14). The expression of mRNA was quantified with miRDeep2 (15). Finally, DEGSeq2 package (16) was used to compare the expression difference of lncRNA and mRNA between the two groups in the R environment. The value of fragments per kilobase of exon per million mapped reads FIGURE 1 | The volcano plot and hierarchical clustering diagram of lncRNAs (A) and mRNAs (B) in MF vs. AMI. In the volcano, the x-axis and y-axis present Log2 fold change and -log10 p-value, respectively. Red and blue color represent upregulation and downregulation, respectively. In the hierarchical clustering diagram, the diagram presents the result of a two-way hierarchical clustering of lncRNAs/mRNAs and sample. Each row represents a differentially expressed lncRNA/mRNA, and each column represents a sample. Differentially expressed lncRNAs/mRNAs clustering tree is shown on the top right corner. The color scale illustrates the relative level of differentially expressed genes: red, below the reference channel; green, higher than the reference.
Frontiers in Cardiovascular Medicine | www.frontiersin.org (FPKM) of each gene/transcript in the sample was calculated according to the comparison results of all samples with the reference genome. The value was regarded as the expression amount of gene/transcript in the sample.

Identification of Differentially Expressed lncRNAs and mRNAs in MF vs. AMI and HF vs. MF
Differential expression analysis of lncRNAs and mRNAs in MF vs. AMI and HF vs. MF was evaluated in the R-bioconductor package DESeq2. Firstly, in order to normalize for sequencing depth, the original read count performed normalization. Then, the hypothesis testing probability (p-value) was obtained by the statistical model. Thirdly, multiple hypothesis test correction was carried out to obtain the corrected p-value [false discovery rate (FDR)]. Those differentially expressed lncRNAs and mRNAs were identified under the criterion of p-value < 0.05 and |log2Fold change| ≥ 1.

Functional Analysis of Common Differentially Expressed mRNAs in MF vs. AMI and HF vs. MF
In order to study the biological function of common differentially expressed mRNAs in MF vs. AMI and HF vs. MF, functional analysis including Gene Ontology (GO) classification (17) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis (18) was performed through CPDB (http://cpdb.molgen.mpg.de/CPDB). Significantly enriched GO and KEGG terms were identified under the criterion of p-value < 0.05.

Nearby Target and Co-expression Analysis
Between lncRNA and mRNA in MF vs. AMI and HF vs. MF In order to identify nearby targeted differentially expressed mRNAs of differentially expressed lncRNAs in MF vs. AMI and HF vs. MF, the mRNAs transcribed within a 100-kb window up-or downstream of lncRNAs (19) were searched by using BEDTools (20). In addition, the Pearson correlation coefficient (cor) of lncRNAs and mRNA was calculated using the cor.test function in R. The correlation pairs with |cor| ≥ 0.98 and p-value < 0.05 were considered in that lncRNAs were significantly coexpressed with mRNAs. Interaction networks of lncRNA-nearby mRNA and lncRNA-co-expressed mRNA were visualized by using Cytoscape (21).

In vitro Validation of Key Differentially Expressed lncRNAs and mRNAs
In total, five patients with AMI (three patients from RNA sequencing), three patients with MF (two patients from RNA sequencing), and five patients with HF (three patients from RNA sequencing) were enrolled. The inclusion criteria of patients with AMI, MF, and HF were consistent with that for RNA sequencing. Ethical approval was obtained from the ethics Frontiers in Cardiovascular Medicine | www.frontiersin.org committee of the Shijiazhuang People's Hospital and informed written consent was obtained from all individuals. Total RNA of the blood samples was extracted according to the manufacturer's protocols. Then, qPCR was performed in an ABI 7300 qPCR system with SYBR R Green PCR Master Mix. All reactions were performed in triplicate. ACTB and GAPDH were used for the internal reference. Relative gene expression was analyzed by the fold change method. p < 0.05 is the significance level.

Electronic Validation and Diagnostic Analysis of Key Differentially Expressed lncRNAs and mRNAs
In order to further validate the expression of the identified key differentially expressed lncRNAs and mRNAs, the GSE21125 dataset (HF cases:normal controls = 18:9) and the GSE48060 dataset (AMI cases:normal controls = 31:21) were used for expression validation. The expression pattern was displayed in the box plots. In addition, we also performed the ROC analyses to assess the diagnostic value of lncRNAs and mRNAs in the above two datasets.

Clinical Information of Enrolled Individuals for RNA Sequencing and in vitro Validation
A total of 14 patients (5 AMI patients, 4 MF patients, and 5 HF patients) were enrolled in both RNA sequencing and in vitro validation. The clinical information of these individuals is presented in Table 1. The chisq.test of gender showed that there FIGURE 2 | The volcano plot and hierarchical clustering diagram of lncRNAs (A) and mRNAs (B) in HF vs. MF. In the volcano, the x-axis and y-axis present Log2 fold change and -log10 p-value, respectively. Red and blue color represent upregulation and downregulation, respectively. In the hierarchical clustering diagram, the diagram presents the result of a two-way hierarchical clustering of lncRNAs/mRNAs and sample. Each row represents a differentially expressed lncRNA/mRNA, and each column represents a sample. Differentially expressed lncRNAs/mRNAs clustering tree is shown on the top right corner. The color scale illustrates the relative level of differentially expressed genes: red, below the reference channel; green, higher than the reference. was no statistical difference of gender between male and female patients in three groups (p-value = 0.638).

Differentially Expressed lncRNAs and mRNAs in MF vs. AMI and HF vs. MF
A total of 248 differentially expressed lncRNAs (Supplementary Table 1) and 936 differentially expressed mRNAs (Supplementary Table 2 (Figure 3). In addition, several enriched signaling pathways in the KEGG analysis were identified, including TNF signaling pathway and cytokine-cytokine receptor interaction (involving IL18R1), fructose and mannose metabolism and HIF-1 signaling pathway (involving PFKFB3), hematopoietic cell lineage and fluid shear stress and atherosclerosis (involving IL1R2), and estrogen signaling pathway (involving FKBP5). All significantly enriched KEGG signaling pathways and involved mRNAs are listed in Table 6.

Nearby Target and Co-expressed Analysis Between lncRNAs and mRNAs in MF vs. AMI and HF vs. MF
A total of 209 lncRNA-nearby targeted mRNA pairs (including 150 lncRNA and 155 mRNAs) and 171 lncRNA-nearby targeted mRNA pairs (including 109 lncRNA and 139 mRNAs) were obtained in MF vs. AMI and HF vs. MF, respectively.    Interestingly, some common lncRNA-nearby targeted mRNA pairs between MF vs. AMI and HF vs. MF were found, such as AC007278.2-IL18R1, AL137145.2-PFKFB3, and LINC01127-IL1R2. In addition, a total of 1214 lncRNA-co-expressed mRNA pairs (including 80 lncRNA and 300 mRNAs) and 2069 lncRNA-co-expressed mRNA pairs (including 78 lncRNA and 288 mRNAs) were obtained in MF vs. AMI and HF vs. MF, respectively. Significantly, a total of 186 common lncRNA-co-expressed mRNA pairs were identified between MF vs. AMI and HF vs. MF. For instance, AC005392.3-IL18R1, AL356356.1-PFKFB3, and MKNK1-AS1-IL1R2 were found. The interaction network of common lncRNA-co-expressed mRNA pairs between MF vs. AMI and HF vs. MF is shown in Figure 4.

In vitro Validation
Five patients with AMI (three patients from RNA sequencing), three patients with MF (two patients from RNA sequencing), and five patients with HF (three patients from RNA sequencing) were included. Based on bioinformatic analysis and literature retrieval results, we selected one differentially expressed lncRNA (PLAC4) and eight differentially expressed mRNA (IL18R1, PFKFB3, IL1R2, IRAK3, FKBP5, LRG1, IKZF2, and RNASE1) for validation ( Figure 5). PLAC4 was upregulated in MF vs. AMI, while it was downregulated in HF vs. MF. Besides FKBP5, all above mRNAs were downregulated in MF vs. AMI, while they were upregulated in HF vs. MF. The validated result was basically in line with the RNA sequencing analysis. In addition, the expression validation of IL18R1, IL1R2, LRG1, and RNASE1 in additional HF and MF patients was also performed (Supplementary Figure 1). The method was similar to the description before. The expression trend of these genes was upregulated. Some genes were not expressed significantly. Small sample size may lead to the bias. Larger numbers of samples are further needed.

Electronic Validation and Diagnostic Analysis of Key Differentially Expressed lncRNAs and mRNAs
Firstly, four key differentially expressed mRNAs including FKBP5, IL1R2, IRAK3, and LRG1 were randomly selected for expression validation in the GSE21125 of HF dataset (Figure 6). FKBP5, IL1R2, IRAK3, and LRG1 were upregulated in HF patients compared with normal controls. This further suggested that upregulation of these mRNAs may play a critical role in the development of HF. In addition, ROC curve analysis was carried out to assess the diagnosis ability of FKBP5, IL1R2, IRAK3, LRG1, RNASE1, and PLAC4 for HF (Figure 7). AUC values of FKBP5, IL1R2, IRAK3, LRG1, RNASE1, and PLAC4 were all more than 0.6, which suggested that they have a potential diagnostic value for HF. Secondly, six key differentially expressed mRNAs including IKZF2, IL1R2, IL18R1, IRAK3, LRG1, and PFKFB3 were randomly selected for expression validation in the GSE48060 of the AMI dataset (Figure 8). IL1R2, IRAK3, LRG1, and PFKFB3 were significantly upregulated in AMI patients compared with normal controls. IKZF2 and IL18R1 were remarkably downregulated in AMI patients compared with normal controls. This indicated that these mRNAs may play a significant role in AMI. In addition, ROC curve analysis  was performed to assess the diagnosis ability of IKZF2, IL1R2, IL18R1, IRAK3, LRG1, PFKFB3, LINC01127, and PLAC4 for AMI (Figure 9). AUC values of IKZF2, IL1R2, IL18R1, IRAK3, LRG1, PFKFB3, LINC01127, and PLAC4 were all more than 0.6, which suggested that they have a potential diagnostic value for AMI. It is noted that IL1R2, IRAK3, LRG1, and PLAC4 could be considered as potential diagnostic markers for both AMI and HF.

DISCUSSION
Up to now, no reports about AC005392.3 have been recorded. The expression level of AC007278.2 is upregulated in peripheral blood mononuclear cells in multiple sclerosis (22). The IL-18 cytokine, coded by interleukin 18 receptor 1 (IL18R1), has been found to be involved in atherosclerosis and cardiovascular complications (23). Bochmann et al. found that IL18R1 is a high-confidence marker for epicardial cells (24). The single nucleotide polymorphism of the IL18R1 gene has been found in cardiovascular disease, such as coronary artery disease (25,26). In addition, the expression level of IL18R1 is significantly altered in systemic and arterial bed-dependent atherosclerosis plaques (27). IL18R1 is upregulated in patients with acute ST elevation MI (28,29). In this study, we first found the association between AC005392.3, AC007278.2, and IL18R1. This suggested that AC007278.2 and IL18R1 may be associated with atherosclerosis and cardiovascular complications, which may provide a novel field in understanding the disease progression of AMI, MF, and HF. Moreover, IL18R1 could be considered as a potential diagnostic biomarker for AMI patients. AL356356.1 is significantly related to the clinical traits of clear cell renal cell carcinoma patients and could serve as a predictor for the disease (30). The highest percentage of candidate fusion transcripts of AL137145.2-6-phosphofructo-2kinase/fructose-2, 6-biphosphatase 3 (PFKFB3) was found in non-small cell lung cancer (31). In cerebrovascular disease of moyamoya disease, AL137145.2 acts as ceRNAs to regulate C-X-C motif chemokine ligand 16 (CXCL16) through miR-107 (31). It is noted that PFKFB3-regulated glycolysis in endothelial cells also plays an important role in vascular sprouting (32,33). PFKFB3 also maintains endothelial cells in a quiescent state to reduce injury-and inflammation-induced pathological processes including atherosclerosis (34)(35)(36). Higher levels of PFKFB3 have been detected in patients with symptomatic atherosclerosis and incident MI (37,38). In the present study, we found that PFKFB3 was co-expressed with AL356356.1. Moreover, it is a nearby target mRNA of AL137145.2. Furthermore, PFKFB3 had a potential diagnostic value for AMI patients. Interaction pairs of AL356356.1/AL137145.2-PFKFB3 may be involved in maintaining the function of endothelial cells in the process of AMI, MF, and HF.
Dysregulated MKNK1-AS1 could predict adverse cardiovascular outcomes in patients with end-stage renal disease (39). It is reported that LINC01127 promotes the development of ovarian tumors (40). Upregulation of LINC01127 has been observed in clear cell renal cell carcinoma (41). Interleukin 1 receptor type 2 (IL1R2) is associated with the inflammatory response, which leads to the development of atherosclerotic plaques (42,43). Additionally, DNA methylation and single nucleotide polymorphism are respectively found in subclinical cardiovascular disease and coronary heart disease (44,45). It is pointed out that IL1R2 is upregulated in vasculitis and AMI patients (46,47). Herein, we found the interaction pairs of MKNK1-AS1/LINC01127-IL1R2 in AMI, MF, and HF, which indicated that these lncRNAs and mRNAs play a role in atherosclerosis. In addition, IL1R2 could be considered as a potential diagnostic biomarker for both AMI and HF patients. Besides interaction with co-expressed/nearby lncRNAs, IL18R1, PFKFB3, and IL1R2 were also significantly involved in some signaling pathways according to the KEGG analysis. For example, IL18R1 was enriched in the TNF signaling pathway and cytokine-cytokine receptor interaction; PFKFB3 was enriched in fructose and mannose metabolism and HIF-1 signaling pathway; IL1R2 was enriched in hematopoietic cell lineage and fluid shear stress and atherosclerosis. The TNF signaling pathway is linked to cardiac remodeling following MI (48). It has been found that TNF-α is involved in the pathogenesis and progression of MF (49). The clinical trial testing of anti-TNF-α drugs has been performed in patients with chronic HF (50,51). It has been demonstrated that cytokine-cytokine receptor interaction plays an important role in the occurrence and development of AMI and HF (47,52). Increased activity of fructose and mannose metabolism is found in acute myocardial ischemia injury (53). Parisi et al. found that HIF-1 had the diagnostic potential effect of ischemia, which can effectively reflect the process of AMI to a certain extent (54). Moreover, lack of HIF-1 causes angiogenesis disorders and MF, which leads to HF (55). Hematopoietic stem cells have a potential regenerative ability in the AMI xenotransplantation (56). In addition, hematopoietic cells can trans-differentiate into cardiac fibroblasts when ischemic injury occurs (57). Westenbrink et al. found that HF affected multiple hematopoietic lineages (58). Reduced shear stress and low wall shear stress gradient are considered to increase the possibility of atherosclerosis (59,60). Atherosclerosis is the most common cause of MI and HF (61,62). Thus, it can be seen that the abovementioned signaling pathways and enriched genes play an important role in cardiovascular disease.
Between MF vs. AMI and HF vs. MF, we also found several other common differentially expressed lncRNAs (such as PLAC4) and differentially expressed mRNAs including interleukin 1 receptor associated kinase 3 (IRAK3), leucine-rich alpha-2glycoprotein 1 (LRG1), ribonuclease A family member 1 (RNASE1), FKBP prolyl isomerase 5 (FKBP5), and IKAROS family zinc finger (IKZF2). PLAC4 is associated with congenital heart disease and hypertrophic cardiomyopathy (63,64). The expression level of IRAK3 is significantly increased 2 days post-MI (65). Upregulation of LRG1 has been found in incident MI (38). In mouse models of MI, LRG1 ablation leads to aggravate MF and heart dysfunction (66). In addition, LRG1 participates in vascular remodeling during HF (67). RNASE1 is upregulated in patients with a history of early MI (68). It is speculated that RNASE1 could be considered as a novel prognostic biomarker for the process of HF after AMI (69). It has been demonstrated that FKBP5 is upregulated and plays a crucial role in endothelial platelet aggregation in AMI patients (70,71). IKZF2 is overexpressed in the left ventricular myocardium (72). In this study, we found that IRAK3, LRG1, RNASE1, and FKBP5 were downregulated in MF vs. AMI, while they were upregulated in HF vs. MF. On the contrary, PLAC4 and IKZF2 were upregulated in MF vs. AMI, while they downregulated in HF vs. MF. Thus, it can be seen that these genes may be involved in the function of vascular endothelial cells.
In addition, FKBP5 was remarkably enriched in the estrogen signaling pathway. In the serum, estrogen receptor β activation reduces the infarct size and lowers the levels of myocardial enzymes (73). In addition, estrogen attenuates MF in HF mainly via estrogen receptor β (74). It is indicated that the estrogen signaling pathway (in which FKBP5 is involved) could reduce the infarct size. It is worth mentioning that PLAC4, IRAK3, and LRG1 had a potential diagnostic value for both AMI and HF in our study.
In summary, identified lncRNA-co-expressed/nearby targeted mRNAs pairs including AC005392.3/AC007278.2-IL18R1, AL356356.1/AL137145.2-PFKFB3, and MKNK1-AS1/LINC01127-IL1R2 and several signaling pathways (TNF signaling pathway, cytokine-cytokine receptor interaction, fructose and mannose metabolism, HIF-1 signaling pathway, hematopoietic cell lineage, fluid shear stress and atherosclerosis, and estrogen signaling pathway) may be involved in the progression from AMI to MF to HF. Significantly, IL1R2, IRAK3, LRG1, and PLAC4 had a potential diagnostic value for both AMI and HF. Our study may provide a theoretical basis for the diagnosis and mechanism of research for the development process of AMI-MF-HF patients, which can be valuable for patient prognosis.
However, there are limitations to our study. Firstly, the sample size in RNA sequencing and qPCR is small. Large numbers of blood samples are further needed for RNA sequencing and in vitro validation. Secondly, we did not validate the expression of identified lncRNAs and mRNAs and explore deeper pathologic mechanisms of the disease progression. More functional and biological validation experiments, such as immunostaining, ELISA, or Western blot, are further needed. In addition, further mice model experiments are needed to validate and investigate the detailed molecular mechanism of these lncRNAs and mRNAs.

DATA AVAILABILITY STATEMENT
We have deposited the RNA seq data in a public, communitysupported repository of GEO (GSE168281).

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by Shijiazhuang People's Hospital. The patients/participants provided their written informed consent to participate in this study.

AUTHOR CONTRIBUTIONS
EW, QC, and YY analyzed the data. LX, XZ, RW, and ZW interpreted the data. SW was a major contributor in writing the manuscript. XH designed the project. All authors read and approved the final manuscript.

FUNDING
This study was funded by Hebei Medical Science Research Project in 2021 "Related molecular mechanisms between lncRNA/mRNA/miRNA/circRNA and development progression from myocardial fibrosis to myocardial infarction to heart failure based on high flux transcriptome sequencing technology" (20210683).