Identification of ceRNA (lncRNA-miRNA-mRNA) Regulatory Network in Myocardial Fibrosis After Acute Myocardial Infarction.

PURPOSE: Myocardial fibrosis (MF) after acute myocardial infarction (AMI) ultimately results in heart failure, which is a serious threat to human life. This study aimed to excavate critical biomarkers associated with MF after AMI.
MATERIALS AND METHODS: RNA-sequencing was performed to obtain differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs) and lncRNAs (DElncRNAs) in AMI and MF after AMI.
RESULTS: Abundant DEmRNAs, DEmiRNAs and DElncRNAs were identified in AMI and MF after AMI. The ceRNA network, which contained 9 lncRNA-miRNA pairs and 9 miRNA-mRNA pairs, was acquired. In AMI, all candidate markers generally exhibited the same pattern as that in our RNA-seq results; while in MF after AMI, except for CENPB, JAK2 and hsa-miR-197-3p, the expression of the others in the qRT-PCR results exhibited the same pattern as that in our RNA-seq results.
CONCLUSION: We speculated that LINC00664/hsa-miR-197-3p/JAK2 and GAS6-AS1/SNHG22/hsa-miR-135a-5p/CENPB/BCL9L interaction pairs may serve as potential biomarkers in MF after AMI.

Introduction

Acute myocardial infarction (AMI) is a life-threatening cardiovascular disease with high mortality and morbidity globally, featured by myocardial necrosis caused by acute or persistent ischemia and hypoxia of cardiac muscle.1 Myocardial fibrosis (MF) is a main pathological changes appear in the heart after AMI, characterized by increasing cardiac fibroblasts activity, increasing the risk of heart failure.2 Thus, preventing and reversing cardiac fibrosis may have potential therapeutic effects on improving heart function and even reducing mortality risk.

With the rapid development of RNA-sequencing technology, the roles of non-coding RNAs (ncRNAs) in cardiovascular disease has been extensively studied. Chistiakov et al have highlighted the important role of miRNA in cardiac function under physiological conditions and cardiac dysfunction under pathological conditions.3 Wang et al concluded the regulatory functions or diagnostic potential of miRNAs, siRNAs, lncRNAs and circRNAs in cardiovascular diseases.4 Similarly, a review indicated that ncRNAs, including lncRNAs, miRNAs and circRNAs, are emerging key regulators of gene expression of several physiological and pathological processes in AMI, as well as lncRNA/circRNA-miRNA-mediated interaction (ie, competitive endogenous RNA (ceRNA) mechanism).5 Thum reviewed the role of ncRNAs in MF and suggested that ncRNAs can serve as powerful therapeutic targets to treat MF.6 PLAUR was found to play a vital role in transducing inter-cellular signals from endothelial cells and fibroblast cells to intra-cellular pathways of myocardial cells, leading to gene expression involved in cellular response to hypoxia.7 Zhuo et al indicated that lncRNA SNHG8 may serve as diagnostic or prognostic biomarkers of AMI.8 It was reported that miR-223 enhanced cell proliferation, migration, and differentiation in cardiac fibrosis, thus mediated cardiac fibrosis after MI partially via the involvement of RASA1.9 Based on the above-mentioned findings, therefore, we performed this current work.

In current study, differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs) and lncRNAs (DElncRNAs) in AMI and MF after AMI were acquired, respectively. Then, shared DEmRNAs/DEmiRNAs/DElncRNAs in AMI and MF after AMI were further determined. Finally, a ceRNA (DElncRNA-DEmiRNA-DEmRNA) regulatory network was constructed to explore the potential molecular mechanisms of MF after AMI.

Materials and Methods

Subjects and Samples

The cohort subjected to RNA-Seq comprised 3 patients with AMI, 3 patients with MF after AMI and 3 healthy individuals. The inclusion criteria of AMI patients were as follows: (1) time of chest pain or distress >30min within 24h, 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. 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 admission, at discharge, and 6 months after MI. Those patients diagnosed with MF (1 year after MI) were enrolled. All samples were collected after obtaining written informed consent from every participant. This study was approved by the ethics committee of Shijiazhuang People’s Hospital and performed in accordance with the Declaration of Helsinki. Peripheral whole blood (2.5 mL) drawn from each subject was collected in PAXgene® RNA blood tubes. Specifically, the blood sample from patients with AMI was collected within 12 h after onset of AMI.

RNA Isolation and Sequencing

According to the manufacturer’s protocol, RNA was extracted with PAXgene blood RNA kit. With Agilent 2100 Bioanalyzer (Agilent RNA 6000 Nano Kit), the concentration, integrity and RNA integrity number (RIN) values of RNA were assessed. After total RNA DNase I treatment, 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, DNBSEQ platform (PE100 strategy) was used to perform RNA sequencing for mRNA and lncRNA. 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%. HISAT2 was applied to align the clean reads with the human reference genome (GRCh38). Expression of mRNAs and lncRNAs was normalized and outputted with StringTie. The Rfam was used for annotation analysis on measured small RNA. The mature miRNA and miRNA precursor sequences were downloaded from miRBase. The expression of miRNA was quantified with miRDeep2. Finally, a total of 19,951 mRNAs, 1313 miRNAs and 16,054 lncRNAs were obtained.

Identification of DEmRNAs, DEmiRNAs and DElncRNAs in AMI and MF After AMI

DEGseq2 () was applied to identify DEmRNAs, DEmiRNAs and DElncRNAs in AMI and MF after AMI with |log2FC| ≥ 1 and p-value < 0.05. Hierarchical clustering analysis of DEmRNAs, DEmiRNAs and DElncRNAs was performed with R () package “pheatmap”. Then, shared DEmRNAs in AMI and MF after AMI were further identified with Venny 2.1.0 (). CPDB () was used to perform GO and KEGG enrichment analysis for up-regulated and down-regulated shared DEmRNAs in AMI and MF after AMI. A value of p < 0.05 was considered to be represented statistically significant.

CeRNA (DElncRNA-DEmiRNA-DEmRNA) Regulatory Network

Six bioinformatic algorithms (PITA, RNA22, miRmap, microT, miRanda and PicTar) were used to predict the targeted DEmRNAs of DEmiRNAs. Then, the DEmiRNA-DEmRNA pairs recorded by ≥4 algorithms were remained for network construction. DIANA-LncBase v2.0 was applied to predict DElncRNA-DEmiRNA interaction pairs. Based on ceRNA theory, the ceRNA (DElncRNA-DEmiRNA-DEmRNA) regulatory network was constructed with Cytoscape, by using a combination of lncRNA-miRNA pairs and miRNA-mRNA pairs.

Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) Validation

Thirteen blood samples were obtained from 5 patients with AMI (3 for RNA sequence), 3 patients with MF after AMI (2 for RNA sequence) and 5 healthy controls (3 for RNA sequence). Total RNA was isolated with the Trizol reagent (Invitrogen, USA) following manufacturer’s protocol. The qRT-PCR reactions were performed in ABI 7300 Real-time PCR Detection System with SuperReal PreMix Plus (SYBR GREEN) (Invitrogen, USA). Relative gene expression was analyzed by 2−ΔΔCT method. The human GAPDH and ACTB were used as endogenous controls for mRNA and lncRNA expression in analysis. The human U6 was used as endogenous controls for miRNA in analysis. The PCR primers were displayed in .

Validation in the GEO Database

To further validate our findings, GSE66360, including miRNA and mRNA expression profiles and consisting of 49 patients experiencing acute myocardial infarction and 50 healthy cohorts, was downloaded from GEO database. Then, the expression levels of selected mRNAs and lncRNAs were validated with GSE66360.

Results

Clinical Characteristics of Enrolled Individuals for RNA Sequencing and qRT-PCR Validation

A total of 14 individuals (5 patients with AMI, 4 patients with MF after AMI, and 5 healthy controls) were enrolled in RNA sequencing and qRT-PCR validation. The clinical information of these individuals was presented in .

Compared to normal controls, 1990 DEmRNAs (1047 up- and 943 down-regulated) were identified in AMI, 1477 DEmRNAs (290 up- and 1187 down-regulated) were identified in MF after AMI, respectively (Figure 1A and B). Of these, PFKFB2 and SLC4A10, TVP23C-CDRT4 and UBAP1L was the most up- and down-regulated DEmRNA in AMI and MF after AMI, respectively (Table 1). A total of 509 shared DEmRNAs (84 up- and 425 down-regulated) in AMI and MF after AMI were acquired ().

Figure 1

The heatmap of the DEmRNAs in AMI (A) and MF after AMI (B), DEmiRNAs in AMI (C) and MF after AMI (D), and DElncRNAs in AMI (E) and MF after AMI (F). Row and column represented DEmRNAs/DEmiRNAs/DElncRNAs and samples, respectively. The color scale represented the expression levels.

Table 1

Top 10 Up- and Down-Regulated DEmRNAs

Symbol	log2fc	p-value	FDR	Regulation
AMI vs Nor
PFKFB2	3.696883	3.30E-34	4.33E-30	Up
CD177	4.177989	9.71E-26	6.37E-22	Up
IL18R1	3.83739	1.83E-23	8.01E-20	Up
BMX	3.691837	1.25E-21	4.09E-18	Up
DAAM2	4.74454	7.00E-20	1.84E-16	Up
IRAK3	2.869636	8.59E-20	1.88E-16	Up
CNTNAP3	3.313797	1.12E-16	2.11E-13	Up
TLR4	2.40295	1.42E-16	2.33E-13	Up
OLAH	7.821776	2.68E-16	3.91E-13	Up
TPST1	3.137414	9.84E-16	1.29E-12	Up
SLC4A10	−5.00668	2.25E-15	2.46E-12	Down
CTSW	−2.68749	2.82E-14	1.85E-11	Down
SMPD3	−4.00009	2.40E-13	1.26E-10	Down
IL32	−2.28019	5.48E-11	1.44E-08	Down
IL2RB	−2.48908	1.30E-10	2.99E-08	Down
GZMM	−3.20545	1.86E-10	4.00E-08	Down
CD8B	−2.21323	2.39E-10	4.97E-08	Down
TICAM2	−4.03595	3.22E-10	6.40E-08	Down
NKG7	−2.50365	4.97E-10	9.07E-08	Down
RORC	−3.29499	2.05E-09	2.89E-07	Down
MF after AMI vs Nor
TVP23C-CDRT4	2.957244	6.58E-08	0.000144	Up
EVI2A	1.813807	2.76E-07	0.000377	Up
LIN7A	1.568656	3.64E-07	0.000443	Up
TFEC	1.686724	5.73E-07	0.000608	Up
NEIL3	1.649259	3.84E-06	0.002807	Up
ANKRD22	2.064817	7.93E-06	0.00436	Up
AP1S2	1.501505	9.73E-06	0.005077	Up
TGFBR1	1.322111	1.09E-05	0.005171	Up
TMEM170B	1.317456	1.14E-05	0.005196	Up
CMC1	1.806918	1.44E-05	0.006087	Up
UBAP1L	−4.19409	2.78E-22	3.05E-18	Down
SARM1	−3.21831	4.97E-21	2.72E-17	Down
CDRT4	−4.16786	6.15E-12	2.25E-08	Down
RNF182	−3.10996	5.71E-08	0.000144	Down
PLEKHB1	−2.85527	1.26E-07	0.000231	Down
JPT2	−2.36819	1.81E-07	0.000283	Down
RAP1GAP	−2.87297	6.10E-07	0.000608	Down
KCNAB2	−1.83602	1.05E-06	0.000963	Down
TMEM120B	−1.89758	1.90E-06	0.001605	Down
MFNG	−1.54803	3.63E-06	0.002807	Down

Abbreviations: DEmRNAs, differentially expressed mRNAs; FC, fold change, FDR, false discovery rate; AMI, acute myocardial infarction; MF, myocardial fibrosis.

Continue.

Significantly enriched GO terms and KEGG pathways of up-regulated (A–D) and down-regulated (E–H) shared DEmRNAs in AMI and MF after AMI. (A and E) BP, biological process; (B and F) CC, cellular component; (C and G) MF, molecular function; (D and H) KEGG pathways. The x-axis shows counts of DEmRNAs enriched in GO terms or KEGG pathways and the y-axis shows GO terms or KEGG pathways. The color scale represented -lg p-value.

Compared to normal controls, 88 DEmiRNAs (22 up- and 66 down-regulated) were identified in AMI, 36 DEmiRNAs (16 up- and 20 down-regulated) were identified in MF after AMI, respectively (Figure 1C and D). Of these, hsa-miR-223-3p and hsa-miR-6802-3p, hsa-miR-1299 and hsa-miR-197-3p was the most up- and down-regulated DEmiRNA in AMI and MF after AMI, respectively (Table 2). A total of 22 shared DEmiRNAs (6 up- and 16 down-regulated) in AMI and MF after AMI were acquired ().

Table 2

Top 10 Up- and Down-Regulated DEmiRNAs

ID	log2FC	p-value	FDR	Regulation
AMI vs Nor
hsa-miR-223-3p	1.028602	8.25E-07	0.000209	Up
hsa-miR-100-5p	2.550814	4.54E-05	0.002554	Up
hsa-miR-143-3p	1.598589	7.45E-05	0.003542	Up
hsa-miR-340-5p	1.514196	0.000274	0.009894	Up
hsa-miR-199b-5p	1.881356	0.000719	0.021392	Up
hsa-let-7f-2-3p	1.287981	0.001282	0.025949	Up
hsa-miR-1179	1.971241	0.001811	0.028291	Up
hsa-miR-135a-5p	1.143364	0.003014	0.040129	Up
hsa-miR-133b	3.440573	0.004444	1	Up
hsa-miR-3688-5p	1.554589	0.004881	0.056126	Up
hsa-miR-6802-3p	−3.33479	1.10E-07	5.58E-05	Down
hsa-miR-1306-5p	−3.3865	1.42E-06	0.000239	Down
hsa-miR-6511b-3p	−2.81604	2.83E-06	0.000358	Down
hsa-miR-1229-3p	−3.23234	6.25E-06	0.000633	Down
hsa-miR-6741-3p	−3.2005	3.14E-05	0.002321	Down
hsa-miR-4433a-5p	−4.61984	4.08E-05	0.002554	Down
hsa-miR-4685-3p	−3.15279	7.70E-05	0.003542	Down
hsa-miR-625-3p	−2.84779	0.000478	0.016132	Down
hsa-miR-1908-5p	−3.22467	0.000712	0.021392	Down
hsa-miR-197-3p	−2.47058	0.00089	0.024591	Down
MF after AMI vs Nor
hsa-miR-1299	3.003464	0.000208	0.083425	Up
hsa-miR-3074-5p	2.0904	0.000384	0.115634	Up
hsa-miR-362-5p	1.169846	0.001426	0.255323	Up
hsa-miR-1179	1.729352	0.01264	0.813026	Up
hsa-miR-371b-5p	3.075288	0.013425	0.813026	Up
hsa-miR-3690	1.676461	0.018571	0.851651	Up
hsa-miR-135a-5p	1.445663	0.02038	0.877059	Up
hsa-miR-132-3p	1.31628	0.022345	0.889135	Up
hsa-miR-643	1.936435	0.022874	0.889135	Up
hsa-miR-152-3p	1.08017	0.029264	0.992927	Up
hsa-miR-197-3p	−2.72354	0.000121	0.083425	Down
hsa-miR-6802-3p	−2.79884	0.000201	0.083425	Down
hsa-miR-1229-3p	−2.83905	0.001058	0.254904	Down
hsa-miR-4433a-5p	−3.31172	0.001736	0.255323	Down
hsa-miR-4433b-5p	−2.60415	0.001845	0.255323	Down
hsa-miR-6741-3p	−2.52569	0.001951	0.255323	Down
hsa-miR-339-5p	−2.63087	0.002119	0.255323	Down
hsa-miR-144-3p	−3.65992	0.003463	0.37935	Down
hsa-miR-1306-5p	−2.75723	0.004393	0.441152	Down
hsa-miR-1249-3p	−2.6606	0.008338	0.717627	Down

Abbreviations: DEmiRNAs, differentially expressed miRNAs; FC, fold change, FDR, false discovery rate; AMI, acute myocardial infarction; MF, myocardial fibrosis.

Compared to normal controls, 473 DElncRNAs (315 up- and 158 down-regulated) were identified in AMI, 502 DElncRNAs (212 up- and 290 down-regulated) were identified in MF after AMI, respectively (Figure 1E and F). Of these, DLEU1 and FP236383.3, SBF2-AS1 and CTBP1-DT was the most up- and down-regulated DElncRNA in AMI and MF after AMI, respectively (Table 3). A total of 90 shared DElncRNAs (39 up- and 51 down-regulated) in AMI and MF after AMI were acquired ().

Table 3

Top 10 Up- and Down-Regulated DElncRNAs

Symbol	log2FC	p-value	FDR	Regulation
AMI vs Nor
DLEU1	1.257401	0.000113	0.002616	Up
MRVI1-AS1	1.529697	0.000261	0.004902	Up
ZFAS1	1.458439	1.36E-06	9.85E-05	Up
LINC00174	1.268871	7.73E-07	6.34E-05	Up
AL049647.1	1.951559	0.001781	0.021287	Up
LINC02363	1.287814	0.005704	0.047997	Up
AC005037.1	1.577588	3.59E-05	0.001244	Up
CYYR1-AS1	4.441044	5.40E-05	0.001645	Up
LINC00862	1.537288	0.000675	0.010503	Up
AL645568.1	1.224525	5.82E-05	0.001724	Up
FP236383.3	−2.51449	7.76E-13	3.31E-10	Down
AC005944.1	−2.54495	1.02E-10	2.56E-08	Down
AC010619.2	−3.17969	6.82E-10	1.16E-07	Down
SNHG22	−2.07379	3.09E-09	4.26E-07	Down
CTBP1-DT	−3.42961	9.28E-09	1.24E-06	Down
AC109460.3	−2.05074	3.69E-07	3.43E-05	Down
AP003068.1	−2.42275	2.23E-06	0.000144	Down
AC010536.1	−2.6355	3.15E-06	0.000184	Down
LINC00235	−2.28106	4.91E-06	0.000269	Down
IL21R-AS1	−1.9457	6.18E-06	0.000314	Down
MF after AMI vs Nor
SBF2-AS1	1.571598	1.47E-05	0.006153	Up
AC244453.2	1.162292	4.07E-05	0.010654	Up
AL121985.1	1.040738	4.03E-05	0.010654	Up
LINC02432	1.75043	6.48E-05	0.013566	Up
AP005329.1	1.916811	0.000153	0.020131	Up
AC007336.1	1.093402	0.000154	0.020131	Up
AC243772.2	1.070516	0.000254	0.026612	Up
AL138820.1	1.419275	0.000287	0.027303	Up
AC006511.5	1.113412	0.000345	0.031444	Up
PRAL	1.110538	0.000427	0.032653	Up
CTBP1-DT	−3.15908	7.53E-08	0.000125	Down
FP236383.3	−1.94351	1.19E-07	0.000125	Down
LEF1-AS1	−2.25841	4.41E-07	0.000308	Down
AC011558.1	−1.93535	1.38E-05	0.006153	Down
DTX2P1-UPK3BP1-PMS2P11	−1.49489	3.61E-05	0.010654	Down
AC136475.3	−1.90278	4.69E-05	0.010902	Down
AC006449.3	−1.86753	8.12E-05	0.015451	Down
AC010619.2	−1.98165	0.000129	0.020131	Down
GAS6-AS1	−1.597	0.000141	0.020131	Down
TNRC6C-AS1	−1.205	0.000136	0.020131	Down

Abbreviations: DElncRNAs, differentially expressed lncRNAs; FC, fold change, FDR, false discovery rate; AMI, acute myocardial infarction; MF, myocardial fibrosis.

Functional Annotation of Shared DEmRNAs in AMI and MF After AMI

For up-regulated shared DEmRNAs group, the result indicated that these genes were enriched in response to stress (p = 9.94E-06), cell activation (p = 1.25E-04), enzyme regulator activity (p = 2.19E-04), Pantothenate and CoA biosynthesis (p = 3.53E-03) and Leishmaniasis (p = 4.09E-03) (Figure 2A–D). These, down-regulated genes were enriched in lymphocyte costimulation (p = 2.45E-05), leukocyte activation (p = 1.28E-04), protein binding (p = 3.69E-05), Th17 cell differentiation (p = 1.85E-03), NF-kappa B signaling pathway (p = 5.95E-03) (Figure 2E–H).

A total of 509 shared DEmRNAs, 22 shared DEmiRNAs and 90 shared DElncRNAs in AMI and MF after AMI were used to construct the ceRNA regulatory network. By overlapping DElncRNA-DEmiRNA interaction pairs and DEmiRNA-DEmRNA interaction pairs, a ceRNA network was acquired. The ceRNA network contained 9 lncRNA-miRNA pairs and 9 miRNA-mRNA pairs, including 9 lncRNAs, 3 miRNAs, and 8 mRNAs (Figure 3).

Figure 3

CeRNA (DElncRNA-DEmiRNA-DEmRNA) regulatory network. The inverted triangles, rectangles and ellipses indicate DElncRNAs, DEmiRNAs and DEmRNAs, respectively. Red and green color represents up-regulation and down-regulation, respectively.

QRT-PCR Validation and Validation in GEO Database

Three DEmRNAs (BCL9L, CENPB and JAK2), three DEmiRNAs (hsa-miR-135a-5p, hsa-miR-144-3p and hsa-miR-197-3p) and two DElncRNAs (IL21R-AS1 and LEF1-AS1) were selected randomly to use for qRT-PCR validation. Based on our RNA-seq results, JAK2 and hsa-miR-135a-5p were up-regulated while BCL9L, CENPB, hsa-miR-144-3p, hsa-miR-197-3p, IL21R-AS1 and LEF1-AS1 were down-regulated in AMI and MF after AMI. In AMI, all these candidate markers generally exhibited the same pattern as that in our RNA-seq results; while in MF after AMI, except for CENPB, JAK2 and hsa-miR-197-3p, the expression of the others in the qRT-PCR results exhibited the same pattern as that in our RNA-seq results (Figure 4). As shown in Figure 5, the expression of BCL9L, CENPB, LEF1-AS1, LINC00664 and SNHG22 were found to exhibit the same pattern as that in our RNA-seq results.

Figure 4

The qRT-PCR results of the DEmRNAs, DEmiRNAs and DElncRNAs in AMI (A) and MF after AMI (B). The x-axis represents the DEmRNAs/DEmiRNAs/DElncRNAs and the y-axis represents log2 (fold change). *Indicates p < 0.05, **Indicates that p < 0.01.

Figure 5

Validation of mRNAs and lncRNAs in GEO. The x-axes represent normal control and AMI groups. The y-axes represent gene expression levels. (A) BCL9L, (B) CENPB, (C) LEF1-AS1, (D) LINC00664, (E) SNHG22.

Discussion

MF is a common hallmark in various heart diseases and predominant cause of morbidity and mortality in the world. The potential biomarkers for cardiovascular disease have been explored widely. It has been reported that ADMA is an independent predictor of cardiovascular morbidity and mortality.10,11 Zinellu et al suggested that plasma ADMA was involved in carotid narrowing after carotid endarterectomy intervention.12 Bivona et al indicated that serum Gal 3 levels rise immediately after AMI and decrease significantly within 5 days after the acute event.13 In addition, it has been concluded that H-FABP is not a reliable marker for AMI diagnosis.14 Despite available drug treatment strategy for MF, there is still an urgent need for development of new biomarkers to prevent and reverse MF after AMI.

Many lncRNAs have been reported to be involved in the development of various diseases. In this study, abundant lncRNAs were identified to be associated with AMI and MF after AMI, such as IL21R-AS1 and LEF1-AS1. Cai et al indicated that IL21R-AS1 was a candidate biomarker for coronary artery disease.15 Overexpressed LEF1-AS1 was detected in coronary artery atherosclerosis, which regulated vascular smooth muscle cell proliferation and migration.16 In agreement with previous studies, we observed IL21R-AS1 and LEF1-AS1 were significantly dyregulated and validated in the qRT-PCR as well in this analysis, which reminds us to focus on the roles of IL21R-AS1 and LEF1-AS1 in AMI and MF after AMI.

High LINC00664 expression was associated with shorter survival in endometrial and clear cell carcinomas patients, whereas, it is converse for high-grade serous carcinomas patients.17 Zhang et al determined an eight-lncRNA signature with high prognostic value for colorectal cancer, including LINC00664.18 Chen et al indicated LINC00664 was a prognostic marker for ovarian carcinoma.19 Chouvarine et al suggested that hypoxia alone could drive cardiac inflammation and identified miR-197-3p as a hypoxia-regulated miRNA.20 McManus et al examined associations of miRNA levels with 6 cardiometabolic (CM) traits and determined that miR-197-3p was significantly associated with body mass index, triglycerides, and systolic and diastolic blood pressures.21 The role of JAK/STAT signaling pathway in the heart has been extensively studied. Negoro et al demonstrated that JAK/STAT signaling pathway was activated in the rat model of AMI.22 Mascareno et al suggested that activation of the JAK/STAT signaling pathway contributed to the pathogenesis of myocardial ischemia.23 Chen et al indicated that activation of JAK2/STAT3 pathway can prevent apoptosis after AMI.24 Zhang et al revealed inhibition of JAK2/STAT3 pathway can attenuate myocardial hypertrophy induced by transverse aortic constriction.25 Chen et al implied that JAK/STAT pathway made contribution to left-atrial-fibrogenesis.26 Eid et al reported that activation of JAK2/STAT3 signaling contributed to reduce MI-induced left ventricle injury.27 In this study, JAK2 was a target of hsa-miR-197-3p, and LINC00664 may function as a ceRNA to capture hsa-miR-197-3p.

Altered expression of GAS6-AS1 has been correlated with poor outcome in various tumors. Han et al indicated that decreased GAS6-AS1 was related to poor prognosis in non-small cell lung cancer patients.28 Zhang et al demonstrated that high GAS6-AS1 expression was associated with advanced TNM stage in gastric cancer.29 Ai et al suggested that up-regulated GAS6-AS1 was negatively correlated with the overall survival of hepatocellular carcinoma patients.30 Similarly, SNHG22, located on chromosome 18q21.1, has been associated with multiple cancers as well. Zhang et al indicted that high expressed SNHG22 predicted poor prognosis in patients with epithelial ovarian carcinoma.31 Fang et al found that SNHG22 exerted carcinogenic effects in triple-negative breast cancer.32 Gao et al suggested that SNHG22 overexpression was correlated with worse overall survival in papillary thyroid cancer patients.33 However, to our best knowledge, the role of GAS6-AS1 and SNHG22 in the heart has not been reported yet.

Recent studies indicated that miR-135a-5p may play important roles in a variety of tumors, including gliomas, lung cancer and head and neck squamous cell carcinoma.34–36 In addition, Wu et al suggested that miR-135a inhibited cardiac fibrosis and it may be a therapeutic target for the amelioration of cardiac fibrosis.37 Xu et al demonstrated that miR-135a-5p negatively regulated hypoxia-induced cardiomyocyte injury in H9c2 cells.38 CENPB, a cell cycle-related gene, was reported to be differentially regulated in the heart tissue in serum-response factor null embryonic stem cells.39 BCL9L is a nuclear co-activator in the Wnt/β-catenin pathway, which was shown to be associated with various cancers. BCL9L knockdown was reported to inhibit pancreatic cancer progression and liver metastasis.40 Huge et al indicated that high expression of BCL9L was associated with poor overall survival in hepatocellular carcinoma patients.41 In addition, Cantù et al suggested a causative link between mutations in BCL9L and human congenital heart disease.42 Both CENPB and BCL9L were targeted by hsa-miR-135a-5p in this analysis. We speculated that GAS6-AS1 and SNHG22 may participate in development of MF after AMI via acting as ceRNAs of CENPB and BCL9L by sponging hsa-miR-135a-5p.

Conclusions

In conclusion, we speculated that LINC00664/hsa-miR-197-3p/JAK2 and GAS6-AS1/SNHG22/hsa-miR-135a-5p/CENPB/BCL9L interaction pairs may serve as potential biomarkers in MF after AMI. Due to small sample size included in this study, further research with a larger cohort is essential to be performed to confirm the conclusion. In addition, the specific mechanism of the lncRNA-miRNA-mRNA axes in MF after AMI need to be validated in vivo and in vitro.