The Novel Non-coding Transcriptional Regulator Gm18840 Drives Cardiomyocyte Apoptosis in Myocardial Infarction Post Ischemia/Reperfusion

Background Ischemia/reperfusion-mediated myocardial infarction (MIRI) is a major pathological factor implicated in the progression of ischemic heart disease, but the key factors dysregulated during MIRI have not been fully elucidated, especially those essential non-coding factors required for cardiovascular development. Methods A murine MIRI model and RNA sequencing (RNA-seq) were used to identify key lncRNAs after myocardial infarction. qRT-PCR was used to validate expression in cardiac muscle tissues and myocardial cells. The role of Gm18840 in HL-1 cell growth was determined by flow cytometry experiments in vitro. Full-length Gm18840 was identified by using a rapid amplification of cDNA ends (RACE) assay. The subcellular distribution of Gm18840 was examined by nuclear/cytoplasmic RNA fractionation and qRT-PCR. RNA pulldown and RNA immunoprecipitation (RIP)-qPCR assays were performed to identify Gm18840-interacting proteins. Chromatin isolation by RNA purification (ChIRP)-seq (chromatin isolation by RNA purification) was used to identify the genome-wide binding of Gm18840 to chromatin. The regulatory activity of Gm18840 in transcriptional regulation was examined by a luciferase reporter assay and qRT-PCR. Results Gm18840 was upregulated after myocardial infarction in both in vivo and in vitro MIRI models. Gm18840 was 1,471 nt in length and localized in both the cytoplasm and the nucleus of HL-1 cells. Functional studies showed that the knockdown of Gm18840 promoted the apoptosis of HL-1 cells. Gm18840 directly interacts with histones, including H2B, highlighting a potential function in transcriptional regulation. Further ChIRP-seq and RNA-seq analyses showed that Gm18840 is directly bound to the cis-regulatory regions of genes involved in developmental processes, such as Junb, Rras2, and Bcl3. Conclusion Gm18840, a novel transcriptional regulator, promoted the apoptosis of myocardial cells via direct transcriptional regulation of essential genes and might serve as a novel therapeutic target for MIRI.


INTRODUCTION
Myocardial infarction (MI) is a major cause of morbidity and mortality (Hausenloy et al., 2013;Reed et al., 2017) and affects more than 50 million people worldwide (Lopez et al., 2006). MI usually results in cardiomyocytes suffering from hypoxia, which induces a severe inflammatory response, promotes myocyte apoptosis and cardiac fibrosis, and finally leads to worse cardiac function and even heart failure (Matyas et al., 2017;Zhang et al., 2017;Ong et al., 2018). Thus, identifying key dysregulated programs and factors during MI is of great importance for the development of potential novel therapeutic strategies.
Emerging evidence has demonstrated that long non-coding RNAs (lncRNAs) display diverse cellular functions and widely participate in both physiological and pathophysiological processes (Li and Sun, 2020;Sayad et al., 2020;Zhang et al., 2020). LncRNAs are a diverse class of RNA transcripts longer than 200 nucleotides that have limited protein coding capacity (Zhou et al., 2020). Most lncRNAs do not engage in translation but exert critical roles in multiple regulatory functions, including microRNA sponging, transcriptional regulation, translational regulation, and protein-protein interactions. Notably, the transcriptional regulatory roles of lncRNAs have emerged as major functions in many key developmental processes. These lncRNAs can function by directly binding to histone-associated complexes, transcription factors, coregulators, or RNA polymerase II within distinct cis-regulatory regions.
Recently, lines of evidence have indicated that lncRNAs also play essential roles in the occurrence and development of cardiovascular diseases (Viereck et al., 2016;Kenneweg et al., 2019;Lu et al., 2019;Yang et al., 2019). For example, lncRNA MIAT is involved in the regulation of MI through the regulation of the ERK/MAPK signaling pathway (Fan et al., 2019).  reported that lncRNA HOTAIR protects MI rats by sponging miR-519d-3p. Wang et al. (2015) found that lncRNA AFP regulates MI by targeting miR-188-3p. Moreover, lncRNAs SLC8A1-AS1 , ZFAS1 (Jiao et al., 2019), and GAS5  play regulatory roles in the occurrence and development of MI in different ways. However, lncRNAs acting as non-coding transcriptional regulators in MI are rarely reported.
In this study, we identified a new lncRNA, Gm18840, as a non-coding transcriptional regulator that was significantly upregulated in MI myocardial tissue compared with the normal myocardial tissue of mice. With a combination of RNA sequencing (RNA-seq), bioinformatic analysis, and several molecular and cellular biology experiments, we elucidated novel mechanisms of MI.

MI Mouse Model
Animal experiments were approved by the Animal Care and Use Committee of Affiliated Liutie Central Hospital & Clinical Medical College of Guangxi Medical University. C57BL/6 male mice (8 weeks old) were obtained from Shanghai SLAC Laboratory Animal Co., Ltd. (Shanghai, China), fed under specific pathogen-free conditions, and randomly treated with MI. To anesthetize the mice, 50 µg/g body weight ketamine and 60 µg/g body weight sodium pentobarbital were injected intraperitoneally into the mice (n = 7). Afterward, the mice were intubated and placed onto a rodent ventilator. A 5mm incision was cut between the fourth and fifth intercostal spaces of the chest cavity. The left-anterior descending coronary artery (LAD) was observed by a microscope, followed by ligation using a 7-0 Prolene suture. Mice in the sham group were subjected to the same procedure without suture but with a moving suture behind the LAD artery. Finally, CO 2 asphyxiation was used to anesthetize and euthanize the mice to collect the hearts.

Cell Culture
The mouse cardiomyocyte cell line HL-1 was purchased from the American Type Culture Collection (ATCC) and cultured in a Claycomb medium (51800C, Signal-Aldrich) with 10% FBS, 0.1 mM norepinephrine, 2 mM L-glutamine, and antibiotics (100 U/ml penicillin, 100 µg/ml streptomycin). Cells were maintained in a humidified incubator (Thermo Fisher Scientific) at 37 • C with 5% CO 2 .

Establishment of an Ischemia/Reperfusion Cell Model
HL-1 cells were grown to 80% confluence and maintained in a quiescent state for 12 h. To generate hypoxic conditions, cells were transferred to an incubation chamber (Billups-Rothenberg MIC-101) and flushed with a hypoxic gas mixture (95% N2, 5% CO 2 ) for 24 h. Subsequently, the cells were subjected to reoxygenation treatment at 37 • C for 2 h to induce the establishment of an ischemia/reperfusion (I/R) model.

RNA Isolation and Quantitative Reverse Transcription PCR (qRT-PCR)
Total RNA was extracted with a TRIzol reagent (Life Technologies) according to the instructions of the manufacturer. qRT-PCR was performed using an AceQ Universal SYBR qPCR Master Mix (Vazyme) according to the protocols of the manufacturer. The relative expression of genes was normalized to GAPDH or U6 and expressed as fold change compared with the internal standard via the 2 − CT method in this study. The primers are shown in Table 1.

Apoptosis Analysis
To determine myocardium apoptosis levels, cultured HL-1 cells were collected after transfection and rinsed in chilled PBS, followed by double staining with an Annexin V fluorescein isothiocyanate (FITC)/propidium iodide (PI) detection kit (Invitrogen) in the dark for 15 min. Then, the rates of apoptotic cells were detected by a flow cytometer (Beckman Coulter, Brea, CA, United States).

RNA Sequencing
RNA was extracted using a TRIzol reagent (Life Technologies) according to the instructions of the manufacturer. The RNAseq library was constructed using the VAHTS mRNA-seq V2 Library Prep Kit for Illumina R (NR601, Vazyme) according to the instructions of the manufacturer. Sequencing was performed using Illumina's NovaSeq platform according to the instructions of the manufacturer. We have uploaded the RNA-seq data to a public dataset, SAR (PRJNA669226).

RNA Pulldown Assay
RNA pulldown assays were performed using the Pierce TM Magnetic RNA-Protein Pull-Down Kit (#20164, Thermo Fisher Scientific, United States) according to the instructions of the manufacturer. Full-length Gm18840 RNA was obtained using in vitro transcription with the T7 High Yield RNA Transcription Kit (TR101, Vazyme), after which the RNA was biotinylated using the Pierce TM RNA 3 End Desthiobiotinylation Kit (# 20163, Thermo Fisher Scientific, United States). Cell extracts were incubated with RNA for 20 min, followed by the addition of nucleic acid-compatible streptavidin magnetic beads (Thermo Fisher Scientific) for further incubation. After washing five times, Gm18840-associated proteins, which were retrieved from beads, were subjected to SDS-PAGE and silver staining. Subsequently, the protein bands were excised and identified by liquid chromatography-mass spectrometry (LC-MS/MS).

RNA Immunoprecipitation Assay
The RNA Immunoprecipitation (RIP) assay was conducted with the Magna RIP TM RNA-Binding Protein Immunoprecipitation Kit (Millipore, United States) according to the instructions of the manufacturer. Briefly, cells were lysed in a lysis buffer containing a protease inhibitor cocktail and an RNase inhibitor. Then, cell extracts were incubated with magnetic beads conjugated with control IgG or anti-Tim50 antibody (Proteintech Group, China). Immunoprecipitated RNAs were purified and quantified by qRT-PCR.

Luciferase Reporter Analysis
The sequences of Gm18840 and Junb, Rras2, and Bcl3 3 -UTRs were amplified using qRT-PCR and inserted into the pGL3 luciferase vector (Promega). HL-1 cells were transfected with a firefly luciferase vector (100 ng) and a Renilla luciferase expression vector (10 ng). Cells were transiently cotransfected with luciferase-reporter plasmids (Promega) with Lipofectamine 2,000 according to the instructions of the manufacturer. After cotransfection with the indicated plasmids for 48 h, the Dual-Luciferase Reporter Assay System (Promega) was used. Renilla luciferase (pRL-TK) acted as an internal control to normalize transfection efficiencies.

Bioinformatic Analysis
RNA sequencing libraries were constructed according to the instructions of the manufacturer for the VAHTS Total RNAseq (H/M/R) Library Prep Kit for Illumina (Vazyme, China). The libraries were sequenced on the Illumina Nova-seq 6000 platform. The raw sequenced reads were mapped to the mouse genome (mm10) using the STAR software (Dobin et al., 2013). Genes were annotated using the Ensemble database. The expression of the genes was quantified using HTseq (Anders et al., 2015), and the differentially expressed genes (DEGs) were identified using DEseq2 (Love et al., 2014). Genes with FDR < 0.01 and | Log 2 (fold change)| > 1 were considered DEGs. Gene ontology (GO) analysis was performed using the DAVID database (Huang da et al., 2009). For the ChIRP-seq analysis, raw reads were mapped to the mouse genome (mm10) using the bowtie2 software (Langmead and Salzberg, 2012). Peaks were called with the MACS suite (Zhang et al., 2008). Motif enrichment analysis was performed using the MEME suite (Bailey et al., 2015).

Statistical Analysis
Statistical analyses were performed using SPSS v.23.0 or Prism GraphPad 8.0. For most of the experiments, independent sample t-tests were used to calculate the p values. The experiments, with representative data shown in the figures, were repeated three times independently with similar results obtained. Statistical details and methods used are indicated in the figure legends, text, or methods.

Identification of lncRNAs That Are Differentially Expressed Upon Myocardial Ischemia-Reperfusion Injury
To identify key protein coding and non-coding genes related to myocardial ischemia-reperfusion injury, we performed transcriptomic analysis in myocardial tissues of mice with or without ischemia/reperfusion injury (IRI). As shown in Figure 1A, the cardiac tissues were dissected 7 days post-IRI treatment, and RNA-seq analysis was performed in the heart tissue of mice with or without myocardial ischemia-reperfusion injury (MIRI). The general quality of the RNA-seq experiments is shown in Supplementary Figure 1. A total of 1,936 upregulated protein-coding genes, 2,731 downregulated proteincoding genes, 356 upregulated non-coding genes, and 303 downregulated non-coding genes were identified (FDR < 0.01 and | Log 2 (Fold change)| > 1) (Supplementary Tables 1, 2). Differentially expressed protein-coding genes and lncRNAs were identified. A principal component analysis (PCA) plot showed the clustering of the RNA expression profiles of cardiac tissues in mice with or without IRI ( Figure 1B). In addition, differentially expressed genes in cardiac tissues with or without IRI are shown in the heatmap ( Figure 1C). Figure 1E illustrates the differentially expressed lncRNAs in cardiac tissues upon IRI treatment. GO functional enrichment analysis of these differentially expressed genes in cardiac tissues with or without IRI treatment showed that genes related to development processes, oxidation reduction processes, and anatomical structure development were significantly enriched ( Figure 1D and Supplementary Table 3), supporting the observed ischemia/reperfusion injury in heart tissues. Four lncRNAs, Gm18840, Meg3, Dnm3os, and Pvt1, were shown to be the most differentially expressed in cardiac tissues with or without IRI. We further validated the four top differentially expressed lncRNAs, GM18840, Meg3, Dnm3os, and Pvt1, using qRT-PCR. As shown in Figure 1F, the relative gene expression levels of GM18840, Meg3, Dnm3os, Pvt1, and, especially, Gm18840 were significantly upregulated in cardiac tissues with IRI treatment.

Validation of Identified lncRNAs Using a Simulated IRI Model in HL-1 Cells
To validate the function of GM18840, Meg3, Dnm3os, and Pvt1 in MIRI, we carried out RNA-seq analysis in an in vitro model of simulated ischemia. We investigated the DEGs in HL-1 cells under normoxia or hypoxia. Many genes showed an overlap between the mouse model and the cell line model (Supplementary Figure 2). The heatmap of DEGs between normoxia and hypoxia in HL-1 cells is shown in Figure 2A and the differentially expressed lncRNAs in the heart tissue of the mice with/without MIRI are shown in Supplementary Table 4. GO analysis was performed based on the above DEGs. Similar to those observed in the mouse model, genes related to the cellular responses to hypoxia, angiogenesis, negative regulation of cell proliferation, and positive regulation of cell migration were enriched (Figure 2B). In addition, we detected the relative expression levels of Gm18840, Meg3, Dnm3os, and Pvt1 in the murine cardiac muscle cell line HL-1 upon hypoxia treatment by qRT-PCR. The results revealed that these four differentially expressed lncRNAs, especially Gm18840, were also significantly upregulated in HL-1 cells upon hypoxia treatment ( Figure 2C). We further explored the function of Gm18840 in the regulation of myocardial ischemia-reperfusion injury. We found that the knockdown of Gm18840 significantly decreased hypoxiaresultant apoptosis in HL-1 cells according to flow cytometry experiments (Figures 3A,B). The morphology results showed that the cell density and the degree of damage of the Si-Gm18840 group were significantly improved compared with those of the hypoxia group. The morphology of HL-1 cells under normoxia, hypoxia with negative control (NC) siRNA, or hypoxia with siRNA targeting Gm18840 was consistent with the flow cytometry results (Figure 3C). Furthermore, Gm18840 is a novel non-coding RNA that has never been reported; hence, we chose Gm18840 as a target in the following experiment. The common deregulated genes in the different RNA-seq analyses are shown in Supplementary Figure 2.

Identification and Characterization of Gm18840
From the UCSC Genome Browser, 1 we found that Gm18840 is located on mouse chromosome 6F2 and composed of four exons  (Figure 4A). To identify the full-length Gm18840 transcript in HL-1 cells, a RACE assay was performed, and the PCR products of 5 -and 3 -RACE are shown in Figure 4B. Based on the Sanger sequencing results, the fulllength Gm18840 transcript was found to be 1,471 nt, with an additional 1,119 bp at the 5 end of exon 1 and an additional 320 bp at the 3 end of exon 4 (Figures 4A,C). However, the ORF finder from NCBI 2 failed to predict a protein from the Gm18840 sequence, as determined by RACE ( Figure 4D). The CPAT database 3 also showed that Gm18840 had limited protein-coding potential (coding probability = 0.141872378288) ( Figure 4E).

Interaction of Gm18840 With Histone H2B in HL-1 Cells
In the next step, RNA pulldown assays were performed coupled with LC-MS/MS to investigate the interaction of Gm18840 in HL-1 cells (Supplementary Table 5). As shown in Figure 5A, the Gm18840 lncRNA was hybridized with biotinylated DNA probes, and the interacting protein of Gm18840 was captured. Silver staining results for the proteins pulled down by Gm18840 indicated that Gm18840 can specifically bind to certain proteins ( Figure 5B). Moreover, seven proteins that could directly interact with Gm18840 with high confidence were revealed by the LC-MS/MS results ( Figure 5C). Among them, histone H2B was 2 https://www.ncbi.nlm.nih.gov/orffinder/ 3 http://lilab.research.bcm.edu/cpat/ the most reliable. Furthermore, RIP-qPCR analysis suggested that there was direct binding between the H2B protein and the Gm18840 lncRNA in HL-1 cells ( Figure 5D). Notably, the identified targets were core histone proteins, suggesting the possible transcriptional regulatory function of Gm18840. In addition, we determined the distribution of Gm18840 lncRNA in the cytoplasm and the nucleus. We found that Gm18840 was mainly expressed in the nucleus, further suggesting its potential function in transcriptional regulation ( Figure 5E).

Gm18840 Is Involved in Transcriptional Regulation Through Direct Binding to Regulatory Regions in Cardiac Muscle Cells
To research the direct targets of transcriptional regulation by Gm18840, ChIRP-seq experiments were conducted using biotinylated antisense DNA probe sets. A genome-browser screen showed the binding of Gm18840 in the promoters of representative genes among the regulatory regions of Atp5f1, Wdr77, Prdx2, Junb, Hook2, Ly6i, Ly6a, and Bcl3 ( Figure 6A). In addition, we identified 3,518 Gm18840 targets that directly bind to regions on chromatin (Supplementary Table 6). Of these regions, 13.9% were promoters, whereas others were regions on the gene (Figure 6B). Furthermore, motif enrichment analysis of these regions was performed, and the motifs EWSR1-FLI1, E2F6, ZNF263, ETS, SP1, and MEIS2 were significantly enriched (Figure 6C), implying a coregulation of Gm18840 with these factors. As shown in Figure 6D, the results of  GO analysis of genes associated with Gm18840-bound peaks revealed that Gm18840 plays an important role in the response to hypoxia, positive regulation of cell migration, positive regulation of autophagy, positive regulation of apoptotic processes, and positive regulation of angiogenesis, indicating an essential role for Gm18840 in MIRI. Luciferase reporter assays performed in HL-1 cells showed that Gm18840 directly regulated the transcriptional activity of its targets, such as Junb, Rras2, and Bcl3 ( Figure 6E and Supplementary Table 7), indicating an important function of Gm18840 in gene regulation.
To further understand the effects of Gm18840 on transcriptional regulation, we performed RNA-seq experiments in HL-1 cells under normoxia and HL-1 cells with Gm18840 knockdown. The PCA results showed that the transcriptomes of HL-1 cells under normoxia and HL-1 cells with Gm18840 knockdown under hypoxia were similar ( Figure 7A). In addition, as shown in Figure 7B and Supplementary Table 8, a heatmap displaying the expression of genes differentially expressed in normoxia compared to hypoxia indicated that the knockdown of Gm18840 under hypoxia partially restored gene expression to that under normoxia. With the knockdown of Gm18840 under hypoxia, the relative expression of Junb, Rras2, and Bcl3 increased, indicating that Gm18840 influenced the activity of the promoters of Junb and Bcl3 (Supplementary Figure 3). Gm18840 directly regulated genes that were defined by ChIRPseq binding and differentially expressed upon the Gm18840 knockdown. As shown in Figure 7C, the enriched GO items of Gm18840 directly regulated genes similar to those in hypoxia vs normoxia as identified in Figures 1D, 2B, 6D. These results suggested that Gm18840 partially restored the transcriptional pattern under hypoxia to that under normoxia.

DISCUSSION
MI is one of the major causes of hospitalization and death worldwide (Zidan et al., 2016;Chen et al., 2018). Therefore, studying molecular mechanisms and discovering impactful therapeutic targets are important for the treatment of cardiac ischemia-related disease. LncRNAs belong to the largest category of RNAs and are closely related to various human diseases (Liao et al., 2011;Huo et al., 2017). They are key regulators of various biological processes (Bhan et al., 2017;Kondo et al., 2017).
Recently, various studies have emphasized that multiple lncRNAs play important roles in the occurrence and development of cardiac ischemia-related disease (Jiang and Ning, 2015;Fa et al., 2020;Huang, 2020;Pierce and Feinberg, 2020). Chen et al. (2018) demonstrated that the expression of lncRNA UCA1 is correlated with I/R-induced cardiac injury in rats. Zhu and Zhao (2019) showed that the expression of lncRNA AK139128 in MIRI tissues and myocardial cells is upregulated, and MIR can be alleviated with AK139128 interference by inhibiting myocardial cell apoptosis. Hence, I/R injury might be improved by modulating the expression levels of lncRNAs. In this study, we first identified differentially expressed lncRNAs by performing RNA-seq analysis in the heart tissue of the mice with or without MIRI. Then, validation by qRT-PCR in cardiac tissues and RNA-seq and qRT-PCR detection in I/R HL-1 cells revealed the involvement of Gm18840, which was selected for further investigation in MI.
Apoptosis is a cell death mechanism induced by various stress conditions, such as oxidative stress and endoplasmic reticulum stress. It has been widely demonstrated that cardiac apoptosis is dramatically enhanced following myocardial ischemic injury (Frankenreiter et al., 2017;Tanimoto et al., 2017;Yao et al., 2017). Many studies have shown that lncRNAs, such as H19 , MALAT1 (Guo X. et al., 2019;Sun et al., 2019), and MEG3 (Wu et al., 2018), play multiple functions during myocardial ischemic injury. No previous studies have demonstrated that Gm18840 is related to MIRI. We reported that the regulatory role of lncRNA Gm18840 in cardiomyocyte apoptosis after hypoxia/reperfusion injury is meaningful. We found that Gm18840 could inhibit apoptosis in I/R HL-1 cells. Interestingly, we did not observe a significant increase in Gm18840 in the RNA-Seq results based on HL-1 cells but verified that Gm18840 was significantly upregulated in vivo and in vitro via qRT-PCR. The reason for this may be that there are more cell types in tissues, and the effect is more complicated; considering that the sample size in the RNA-seq experiment was relatively small, false-positive results were prone to occur. The results above provide an impetus to further study the function and mechanism of Gm18840 in MIRI.
The structure and localization of lncRNAs in cells are closely related to their functions and mechanisms (Chen, 2016;MacDonald and Mann, 2020). Gm18840 showed limited proteincoding potential. Nucleocytoplasmic separation showed that Gm18840 was mainly located in the nucleus. This may suggest that Gm18840 is involved in transcriptional regulation in I/R HL-1 cells. At the transcriptional level, lncRNA bridges DNA and protein by binding to chromatin and serving as a scaffold for the modification of protein complexes. Such a mechanism can bridge promoters to enhancers or enhancer-like non-coding genes by regulating chromatin looping, as well as conferring specificity upon histone modifying complexes by directing them to specific loci (Dykes and Emanueli, 2017).
Recent studies on the function of lncRNAs have revealed the direct histone binding and transcriptional regulation functions of lncRNAs. For example, MALAT1, which is an lncRNA, directly interacts with DBC1 and P53 and regulates the acetylation of the P53 protein . Here, we demonstrated that there was direct binding between the H2B protein and Gm18840 lncRNAs in HL-1 cells. In addition, the results showed that Gm18840 was mainly expressed in the nucleus, indicating that Gm18840 directly interacts with the nucleosome-related protein histone H2B and is involved in transcriptional regulation in HL-1 cells. Posttranslational regulation of H2B has been reported to be essential for myocardial development (Robson et al., 2019). Transcriptional regulation plays a vital role in almost all biological processes in myocardial development and disease. A large amount of evidence has shown that lncRNAs interact with histones to play a role in the development of disease (Jain et al., 2016;Sun et al., 2016). To further research the direct targets of transcriptional regulation by Gm18840, ChIRP-seq and luciferase reporter assays were conducted. Our ChIRP-seq and Luciferase reporter assays were performed in HL-1 cells. The representative regions bound by Gm18840, including Junb, Rras2, and Bcl3, were cloned into the pGL3 plasmid. The luciferase constructs were cotransfected with the pcDNA3.1 vector or pcDNA3.1-Gm18840 plasmids. The firefly luciferase activity was normalized to Renilla luciferase activity. **p < 0.01 and ***p < 0.001. RNA-seq results suggested that Gm18840, as a transcriptional modulator, directly binds to the chromatin regulation regions of many essential genes, such as Atp5f1, Junb, Ly6a, and Bcl3. Gm18840 regulated the promoter/enhancer activity of these genes, including Junb, Rras2, and Bcl3. Motif enrichment analysis showed that Gm18840 bound to the region enriched for the EWSR1-FLI1, E2F6, and ZNF263 motifs, indicating that Gm18840 might regulate the transcriptional activity of these transcription factors and contribute to myocardial ischemiareperfusion injury. Nevertheless, the effects of the Gm18840-H2B interaction must be investigated in the future. Gm18840 might interact with H2B and directly bind to the transcriptional regulatory regions of the genes essential for the promotion of apoptosis, such as Junb, Rras2, and Bcl3. Junb is a gene belonging to the AP-1 transcriptional complex that is suppressed and has been reported to suppress cell apoptosis in heart failure (Yan et al., 2018). Sabatasso et al. (2018) confirmed that JunB is an early, sensitive marker for myocardial ischemia. Schmidt et al. (2007) indicated that JunB is a target gene of hypoxia-induced signaling via NF-kappaB. Bcl3 is a proto-oncogene that belongs to the IκB family, and it has been reported that BCL3 can exert an anti-apoptotic effect. Lee et al. (2015) revealed that Bcl3 was upregulated under hypoxic conditions, resulting in the prevention of apoptosis induced by ischemic stress both in vitro and in vivo. We found that low expression of Gm18840 upregulated the expression of Junb and Bcl3, indicating that Gm18840 directly inhibits cell apoptosis in I/R HL-1 cells.
In conclusion, we profiled MIRI-related protein-coding and non-coding RNAs and identified an essential regulator, Gm18840, for MIRI. We demonstrated that Gm18840 was upregulated in MIRI tissues and hypoxic HL-1 cells. These results suggest that Gm18840 might promote myocardial cell apoptosis in MIRI and regulate cardiomyocyte apoptosis after hypoxia/reperfusion injury. However, the molecular mechanism of Gm18840 is still unclear. The role of the Gm18840-mediated regulatory network and its association with other apoptosis factors should be explored in the future.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: SRA Database (PRJNA669226).

ETHICS STATEMENT
The animal study was reviewed and approved by Animal Care and Use Committee of Affiliated Liutie Central Hospital & Clinical Medical College of Guangxi Medical University.

AUTHOR CONTRIBUTIONS
CL, XK, and SX designed the research. CL and XK carried out the analyses. YH and MD performed the molecular cell experiments. JZ, JC, and RY performed the visualization. SX and XK wrote the manuscript. All authors reviewed and approved the final manuscript.