Lipopolysaccharide Alters the m6A Epitranscriptomic Tagging of RNAs in Cardiac Tissue

Abstract

N6-methyladenosine (m⁶A) modification plays important roles in the pathology of a variety of diseases. However, the roles of m⁶A modification in sepsis-induced myocardial dysfunction are not well defined. Rats were divided into control and lipopolysaccharide (LPS)-induced sepsis group. Global m⁶A levels of left ventricle tissue were measured by LC-MS/MS, and transcriptome-wide m⁶A modifications were profiled using epitranscriptomic microarrays (mRNAs and lncRNAs). Bioinformatics analysis was conducted to understand the functional implications of m⁶A modifications during sepsis. Methylated lncRNAs and mRNAs were measured by m⁶A single-base site qPCR. The global m⁶A levels in left ventricle tissue were significantly decreased in the LPS group. While 27 transcripts (23 mRNAs and four lncRNAs) were hypermethylated, 46 transcripts (39 mRNAs and 7 lncRNAs) were hypomethylated in the LPS group. The mRNA expression of writers and readers was significantly decreased in the LPS group. The m⁶A modification of Clec1b, Stk38l and Tnfrsf26 was associated with platelet activation and apoptotic pathways. Moreover, the decrease in m⁶A modification of lncRNA XR_346,771 may be related to cation import in cardiac tissue. Our data provide novel information regarding changes to m⁶A modifications in cardiac tissue during sepsis, and m⁶A modifications might be promising therapeutic targets.

Introduction

Sepsis is a life-threatening condition involving organ dysfunction that is caused by a dysregulated host response to infection (Cohen et al., 2015), It has been reported that approximately 50% of patients with septic shock are diagnosed with septic cardiomyopathy, which causes myocardial dysfunction, and these patients tend to have poor prognoses (Zaky et al., 2014). Myocardial dysfunction is an important factor that contributes to the high mortality rate of sepsis. Therefore, it is necessary to elucidate the pathogenesis of myocardial dysfunction in order to develop new treatment strategies to reduce the mortality rate of sepsis.

To date, a series of reversible posttranscriptional modifications located in RNAs have gained increasing attention. More than 140 chemical modifications have been discovered in RNAs. The modifications range from simple methylation or isomerisation, such as N6-methyladenosine (m⁶A), 5-methylcytosine (m5C), N1-methyladenosine (m¹A), pseudouridine (Ψ), 5-methyluridine (m⁵U),1-methylguanosine, and 7-methylguanosine (m¹G and m⁷G, respectively) and inosine (I), to complex multistep chemical modifications, such as N⁶-threonylcarbamoyladenosine and 5-methoxycarbonyl-methyl-2-thiouridine (mcm5s²U) (Pan, 2018). Among them, the m⁶A is one of the most abundant and influential modifications in eukaryotes. M⁶A modification is involved in regulating a variety of posttranscriptional events, including pre-miRNA processing and RNA stability, translation and alternative splicing (Chokkalla et al., 2019). M⁶A modifications occur via the m⁶A methyltransferases called “writers”; they are removed by the demethylases called “erasers” and are recognized by m6A-binding proteins called “readers”. The orchestrated interplay among writers, erasers, and readers drives the dynamics and outcomes of the m⁶A modification of RNAs. An increasing number of studies have reported that m⁶A modification plays important roles in biological and pathophysiological processes, such as tumorigenesis, embryonic stem cell differentiation, and viral infection (Kennedy et al., 2016). In mouse GC-1 SPG cells, meclofenamic acid inhibits spermatogonial proliferation by affecting CDKs expression through a m⁶A-dependent mRNA degradation pathway (Huang et al., 2019). In liver cancer, overexpression of METTL3 leads to the degradation of suppressor of cytokine signalling two and promotes tumor growth (Lin et al., 2019).

M⁶A modification is also associated with the pathophysiological processes of cardiac differentiation and a variety of cardiovascular diseases. Methyltransferase-like 3 (METTL3) is a major factor involved in abnormal m⁶A modification. Silencing METTL3 enhances the autophagic flux and represses apoptosis in hypoxia/reoxygenation-treated cardiomyocytes (Song et al., 2019). Overexpression METTL3 increased m⁶A levels in mRNAs isolated from transverse aortic constriction (TAC) mice hearts. TAC inducd pathological hypertrophic cellular growth was attenuated in hearts of METTL3-overexpressing mice, as evidenced by the crosssectional area of myocytes (Kmietczyk et al., 2019). In addition, aging caused METTL3/METTL14 down-regulation in aorta and atria in male animals while the differentiation-induced increased level of METTL16 (Arcidiacono et al., 2020). Prabhu M et al. found that downregulation of fat mass and obesity associated protein (FTO) expression in failing mammalian hearts and hypoxic cardiomyocytes increased the m⁶A modification of RNAs and reduced the contractile function of cardiomyocytes (Mathiyalagan et al., 2019). FTO also played a vital role in cardiac contractile function during homeostasis and remodeling. FTO overexpression attenuated the ischemia-induced elevation in m⁶A modification and significantly improved cardiac function of post-myocardial infraction (Mathiyalagan et al., 2019). Berulava et al. reported that differently expressed genes of m⁶A methylation are involved in heart failure development. Bioinformatics analysis has revealed that the differentially m⁶A are mainly involved in metabolism and cardiac signaling (Berulava et al., 2020). However, the regulatory role of the m⁶A modification in cardiac function is still unclear.

The contraction of the ventricle transports blood to the capillaries of the body or lungs, while the contraction of the atria send blood into the ventricle. In addition, the right ventricle mainly supplies blood to the lungs, and the blood in the right ventricle is venous blood. The left ventricle mainly supplies blood to the organs and tissues of the body, and the blood in the left ventricle is arterial blood. Therefore, the function of the left ventricle is very important to maintain organ perfusion. If the left ventricular function declines, it will leads to organ perfusion insufficiency and, eventually, organ failure. In this direction, we focused on the mechanism underlying left ventricle injury in sepsis to find any potential clinical significance for improving organ perfusion. The overall goal of this study was to explore a new layer of epigenetic alterations through the genome-wide screening of the altered m⁶A-tagged transcript profiles in lipopolysaccharide (LPS)-induced myocardial dysfunction. We successfully mapped m⁶A transcripts in left ventricle tissue in LPS-induced sepsis. Some potential roles of the m⁶A modification in the physiological and pathological mechanisms of sepsis were revealed.

Materials and Methods

Animals

Eight-week-old male Wistar rats weighing 250–350 g were used in this study. All the rats used in our study were obtained from Charles River Laboratories (Beijing, China). Water and standard laboratory food were freely available to the animals. The Experimental Animal Welfare Ethics Branch and the Biomedical Ethics Committee of Peking University approved our study protocol (LA 2020343). All the procedures for handling the animals were in accordance with the Guide for the Care and Use of Laboratory Animals (NIH Publication No. 85–23, revised 1996).

After acclimating for 1 week, 10 Wistar rats were randomly divided into two groups. The rats in the LPS group (n = 5) were intraperitoneally injected with LPS (10 mg/kg, 5 mg LPS dissolved in 1 ml 0.9% saline). The rats in the control group (n = 5) were intraperitoneally injected with 0.9% saline (2 ml/kg). 24 hours after LPS injection, mean blood pressure was measured by the tail-cuff method with a noninvasive blood pressure measurement system. The left ventricle tissue samples were immediately transferred to liquid nitrogen and stored at −80°C for preservation.

RNA Extraction and Quality Control

Total RNA was isolated from the left ventricle tissues of the LPS (n = 5) and control (n = 5) groups using TRIzol Reagent (Invitrogen, United States). The quantity and purity of the total RNA samples were measured by a NanoDrop ND-1000 (ThermoFisher, United States).

Quantification of Global m⁶A Levels

Liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based mRNA modification detection was performed according to the Aksomics standard protocol. Briefly, 5 μg of total RNA from the heart tissues of the LPS (n = 5) and control (n = 5) groups was used to isolate mRNA using the NEBNext Poly(A) mRNA Magnetic Isolation Module (NEB, E7490, United States). The purified mRNA was quantified using the Qubit RNA HS Assay kit (ThermoFisher, United States) and digested to single dephosphorylated nucleosides by an enzyme mixture. Pretreated nucleosides solution was deproteinized using Satorius 10,000-Da MWCO spin filter. LC-MS/MS analysis was performed on an Agilent 6460 QQQ mass spectrometer with an Agilent 1260 HPLC system using Multi reaction monitoring (MRM) detection mode (Agilent, United States. The nucleosides were quantified by using retention time and the nucleoside to base ion mass transitions of 268-136 (A) and 282-150 (m⁶A). Quantification was performed in comparison with the standard curve obtained from pure nucleoside standards running with the same batch of samples. The m⁶A level was calculated as the ratio of m⁶A to A based on the calibrated concentrations.

Quantitative Real-Time PCR

After extraction, the total RNA was reverse-transcribed into cDNA using SuperScriptTM III Reverse Transcriptase (Invitrogen) according to the manufacturer’s instructions. Amplification and detection were performed using 2× SYBR Green PCR Master Mix (Arraysta, United States) on a ViiA 7 Real-Time PCR System (Applied Biosystems). The sequences of the gene-specific primers are listed in Supplementary Table S1. GAPDH was used as the control housekeeping gene.

Western Blot Analysis

Protein extraction of left ventricle tissue and the protein concentration was determined by Bradford method (M&C Gene Technology Ltd.). The protein samples were then subjected to SDS-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to polyvinylidene difluoride membrane. The membranes were probed with antibodies of interest. The antibodies used were as follows: METTL3, METTL14, Wilms tumor 1-associated protein (WTAP), and YTH N6-methyladenosine RNA binding protein 1 (YTHDF1), YTHDF3 (all1:1,000; Cell Signal Technology), and β-actin (1:4,000; Abcam).

M⁶A mRNA&lncRNA Epitranscriptomic Microarray

Sample preparation and microarray hybridization were performed based on Arraystar’s standard protocols (Arraystar). Briefly, total RNA from the heart tissues of the LPS (n = 5) and control (n = 5) groups was immunoprecipitated with an anti-m⁶A antibody (Synaptic Systems, 202003). The modified RNAs immunoprecipitated by the magnetic beads were labelled as “IP”, and the unmodified RNAs in the supernatant were labelled as “Sup”. The “IP” and “Sup” RNAs were labeled with Cy5 and Cy3, respectively, as cRNAs in separate reactions using the Arraystar Super RNA Labeling Kit (ArrayStar). The cRNAs were combined and hybridized onto an Arraystar Rat mRNA&lncRNA Epitranscriptomic Microarray (4 × 44K, Arraystar) that contained 27,770 mRNA and 10,582 lncRNA degenerate probes. The hybridized arrays were scanned in two-color channels by an Agilent Scanner G2505C.

Microarray Data Analysis

Agilent Feature Extraction software (version 11.0.1.1, United States) was used to analyze the acquired array images. Each probe signal was evaluated and flagged as present, absent or marginal in at least 5 out of 10 samples. The raw intensities of the IP (immunoprecipitated, Cy5-labeled) and Sup (supernatant, Cy3-labeled) samples were normalized by the average of log2-scaled spike-in RNA intensities. The differentially m⁶A-methylated RNAs between the LPS and control groups were identified by filtering with thresholds of fold change >1.5 and statistical significance (p < 0.05).

Bioinformatics Analysis

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed using some package in the R environment for statistical computing and graphics. The differentially m⁶A-methylated mRNAs, as well as the downstream mRNAs predicted by the competing endogenous RNA (ceRNA) network, were classified into different GO terms and enriched in certain biological pathways. Protein-protein interaction (PPI) analysis was performed using the STRING database (https://string-db.org). The SRAMP database (http://www.cuilab.cn/sramp) was used to predict the number of m⁶A sites on the differentially m⁶A-methylated transcripts (Zhou et al., 2016).

M⁶A Single-Base Site qPCR

Before amplification and detection, the LPS (n = 5) and control (n = 5) group samples were treated with an Escherichia coli toxin and RNA endoribonuclease, MazF. MazF is reported to be sensitive to m⁶A modification within the ACA motif (Imanishi et al., 2017). The MazF-digested mRNA samples and the nondigested samples were subjected to reverse transcription using SuperScriptTM III Reverse Transcriptase (Invitrogen) for qPCR as described above. The targeted lncRNAs and mRNAs were predicted by SRAMP to identify the ACA motif and m⁶A position. The primer sequences specific for the methylated lncRNAs and mRNAs are listed in Supplementary Table S2. The relative expression levels were calculated using the 2^−△△Ct method, and the tested genes were calibrated with MazF as follows:

Competing Endogenous RNA Network Construction

The targeted lncRNAs were verified by m⁶A single-base site qPCR to construct the ceRNA network. Only three steps were used for ceRNA network construction. First, the potential target microRNAs of lncRNAs were predicted with Aksomics’s homemade miRNA target prediction software based on TargetScan and miRanda (Enright et al., 2003; Pasquinelli, 2012). Second, the target genes of the miRNAs involved in the lncRNA-miRNA interaction network were predicted with the database described above. Finally, the lncRNA-microRNA-mRNA interaction networks were constructed by Cytoscape v2.8.3. The targeted mRNAs were also analyzed by GO and KEGG to completely understand the ceRNA effects.

Statistical Analysis

The statistical significance of the GO and KEGG analysis of mRNAs enrichment were calculated by Fisher’s exact test p < 0.05 and -log10(p) transformed as the enrichment score as well as mRNAs predicted by the ceRNA network for GO and KEGG analysis. The significance of the differences in the expression and methylation levels between the LPS and control groups was evaluated with an unpaired two-sided t-test for LC-MS/MS, qRT-PCR, microarray analysis and m⁶A single-base site qPCR, and the recommended p-value threshold was less than 0.05. “*” indicates a significant difference compared with control group (p < 0.05).

Results

The Global Levels of m⁶A Modification Were Decreased

In about 24 h after giving intraperitoneal injection of LPS, we measured the MAP of rats. As previous paper, the rats with a MAP declined to 25–30% or below were chosen as the sepsis model (Nie et al., 2020). Furthermore, signs of shock such as lassitude, tachycardia and a sharp drop in body temperature were observed in all sepsis rats. In the control group, there was no significant change in MAP in rats injected with saline. Previous study showed that there is a positive relationship between MAP and heart function (Yang et al., 2018; Zhang et al., 2019). Therefore, we isolated heart tissue and chose the left ventricle for identification of m⁶A modification levels (Figure 1A). There was no measurement of cardiac function or sepsis physiology in these experiments.

FIGURE 1

To determine the difference in the m⁶A modification levels between the LPS and control groups, we used LC-MS/MS to measure these levels. We found that the global m⁶A modification level in the rat left ventricle tissues from the LPS group was significantly decreased compared with control group (Figure 1B).

The Expression of m⁶A Writers and Readers Was Downregulated

The m⁶A methylation modification is dynamically modulated by RNA methyltransferases, RNA-binding proteins and demethylases. We analyzed the mRNAs level of RNA methyltransferases (writers: METTL3, METTL14, and WTAP), RNA-binding proteins (readers: YTHDF 1 and YTHDF 3) and a demethylase (eraser: FTO). It was found that the expression of the METTL3/METTL14/WTAP/YTHDF 1, and YTHDF 3 in the LPS group was significantly downregulated compared with control group. However, the expression of FTO has no significant difference (Figure 1C). In addition, we measured the protein expression of METTL3/14, WTAP, and YTHDF1/3. We found that protein expressions of these proteins significantly downregulated in LPS group (Supplementary Figure S1). The expressions of METTL3/14, WTAP, and YTHDF1/3 were decreased by 42.2, 64, 47.6, 49.6 and 56.6%, respectively. The results of protein expression and mRNA expression are consistent.

M⁶A Modification Profiles of lncRNAs and mRNAs

Probes specific for 27,770 mRNAs and 10,582 lncRNAs were used to analyze the samples from the LPS and control groups using an m⁶A mRNA and lncRNA Arraystar epitranscriptomic microarray. The results showed that the m⁶A modification levels of 62 mRNAs and 11 lncRNAs were significantly altered in the LPS group compared with control group (fold change >1.5, p-value < 0.05) (Figures 2A,B). The altered lncRNAs and mRNAs were aligned for cluster analysis (Figures 2C,D). Based on their fold changes, we selected 11 lncRNAs with significantly altered m⁶A modification levels and the top 10 up- and downregulated mRNAs to present the results (Tables 1, 2).

FIGURE 2

The data that support the findings of this study are openly available in the GenBank databases under accession number GSE159309.

GO and KEGG Analyses of Differentially Methylated mRNAs

GO analysis of the mRNAs with increased and decreased m⁶A modification levels was performed. Among the mRNAs with decreased m⁶A modification levels, the “cytokine-mediated signaling pathway” in BP, “immune response” in BP and “cytokine activity” in MF were revealed to have the maximum enrichment scores for GO terms (Figures 3A,B; Tables 3, 4).

FIGURE 3

KEGG analysis of the mRNAs with increased and decreased m⁶A methylation modification levels were also conducted. The results showed that the mRNAs with increased m⁶A modification levels were enriched in only one pathway, namely, the “PPAR signaling pathway” (Figure 3C). The mRNAs with decreased m⁶A modification levels were enriched in 27 pathways, and the following main pathways reached higher enrichment scores: “IL-17 signaling pathway”, “TNF signaling pathway” and “Cytokine-cytokine receptor interaction” (Figure 3D).

qPCR Results of a Single-Base m⁶A Site

By using the SRAMP database, we predicted the mRNAs and lncRNAs with differential m⁶A modification in the order from high to low fold change of m⁶A modification. As a result, we identified four lncRNAs and six mRNAs with m⁶A ACA sites and with high confidence scores of qPCR validation of single-base m⁶A sites. The Detail information of lncRNAs and mRNAs detected by m⁶A single-base site qPCR was shown in Table 5. The results showed that among the lncRNAs, m⁶A modification levels were significantly increased in XR_346771. Among the mRNAs, m⁶A modification levels were significantly increased in C-type lectin domain family 1 member B (Clec1b) and tumor necrosis factor receptor superfamily member 26 (Tnfrsf26) and significantly decreased in serine/threonine kinase 38 like (STK38L) (Figure 4).

FIGURE 4

PPI Analysis of Targeted Genes

By employing the STRING database, we performed PPI analysis of Clec1b, Stk38l, and Tnfrsf26 (Figures 5A–C). KEGG analysis revealed that the Clec1b interaction-related proteins were mainly enriched in “platelet activation” and other pathways. The Stk38l interaction-related proteins were mainly enriched in the “hippo signaling pathway”, and the Tnfrsf26 interaction-related proteins were mainly enriched in “apoptosis” (Table 6).

FIGURE 5

ceRNA Analysis of lncRNA XR_346771

A ceRNA network of the lncRNA-microRNA-mRNA interaction network was constructed for XR_346,771 (Figure 6A). GO and KEGG analyses were performed on the predicted mRNAs (Figures 6B,C). The “inorganic cation import across plasma membrane”, “transmembrane transport”, and “import across plasma membrane” in BP were significantly enriched for GO terms (Figure 6B). FAM155B and PICALM are downstream proteins related to transmembrane transport that are regulated by XR_346771, and the corresponding sponged miRNAs are shown in Table 7.

FIGURE 6

Discussion

Our study suggested that the m⁶A modification level was decreased in septic heart tissues. Our genome-wide screening of m⁶A-tagged transcript profiles indicated that the m⁶A modification levels in 39 mRNAs and 7 lncRNAs were significantly decreased and that the m⁶A modification levels in 23 mRNAs and four lncRNAs were significantly increased.

Pathway analyses of mRNAs with altered m6A modification were conducted in our study. The results showed that the mRNAs with decreased m⁶A modification levels were enriched in many pathways related to immune and inflammatory responses. In recent years, the roles of m⁶A modification in autoimmune and inflammatory diseases have attracted substantial attention (Rubio et al., 2018). METTL3 is a key enzyme of m⁶A methylation modification and is involved in immune and inflammatory regulation. Feng et al. found that METTL3 inhibits inflammation by affecting the alternative splicing of MyD88 (Feng et al., 2018). METTL3 expression is significantly upregulated in blood samples from patients with rheumatoid arthritis. In addition, METTL3 attenuates LPS-induced macrophage-mediated inflammation through the NF-κB pathway (Wang et al., 2019a). The m⁶A “reader” YT521-B homology domain family 2 (YTHDF2) inhibits HCC cells and tumor vasculature by degrading the mRNAs of IL11 and serpin family E member 2 (Wang et al., 2019b). These findings indicate the important role of m⁶A modification in inflammatory and immune disorders. Sepsis is mainly caused by the host’s unbalanced response to infection. In sepsis, microbial infections or necrotic tissues release a large amount of harmful substances, leading to the activation of the systemic immune response and the excessive activation of immune cells. The excessive release of cytokines is destructive. There is currently no effective treatment for sepsis because in the past few decades, attempts to use anti-inflammatory treatments to limit the tissue damage caused by excessive inflammation have failed (Zeni et al., 1997). Our study suggests that alterations in m⁶A modification in sepsis are closely related to inflammatory and immune responses. Since we study cardiac tissue, we speculate that m⁶A modification plays an important role in sepsis-induced cardiac function dysfunction.

We selected four lncRNAs and six mRNAs for the qPCR validation of single-base m⁶A sites in our study. The results showed that among the lncRNAs, the m⁶A modification levels were significantly increased in XR_346771. Among the mRNAs, the m⁶A modification levels were significantly increased in Clec1b and Tnfrsf26 and significantly decreased in Stk38l. These results are consistent with those of high-throughput sequencing; therefore, we attempted to conduct an in-depth analysis of this lncRNA and the three mRNAs.

First, we performed PPI analysis of Clec1b, STK38L, and Tnfrsf26 to investigate the downstream possible regulatory genes and biological functions of these proteins. The Clec1b-interacting proteins are enriched in the pathway of “platelet activation”. C-type lectin-like receptor 2 (CLEC-2) is a protein encoded by the Clec1b gene. Platelets play a critical role in innate and adaptive immunity (Semple et al., 2011). Platelets can alleviate LPS-induced septic shock by regulating the ability of macrophages to engulf and kill bacteria (Xiang et al., 2013). In two mouse models of sepsis (intraperitoneal LPS injection and cecal ligation and puncture), platelet CLEC-2 reduced the severity of sepsis by controlling the migration of monocytes/macrophages to the infection site, the expression of inflammatory mediators, and damage to organs (Rayes et al., 2017). In our sepsis model, we found that the m⁶A modification level of Clec1b mRNA was significantly increased. Considering that the degree and pattern of m⁶A modification may affect the transport, splicing, storage, translation, stability, and decay of mRNAs, we speculated that the elevated level of Clec1b m⁶A modification affects the functions of Clec1b and associate with activation of platelets. The interacting proteins of Clec1b were glycoprotein VI (Gp6), lymphocyte cytosolic protein 2 (Lcp2), and spleen tyrosine kinase (Syk). These proteins might be used as target proteins for in-depth research.

Our study suggests that the m⁶A modification level of Tnfrsf26 mRNA is increased and of Stk38l mRNA is decreased. According to PPI network analysis, Stk38l-interacting proteins are enriched in the “Hippo signaling pathway” pathway. Hippo signaling pathway may promote cell differentiation and death and inhibit cell proliferation (Chang et al., 2019). Tnfrsf26-interacting proteins are mainly enriched in the pathway of apoptosis. Sepsis-induced cardiac dysfunction has been reported by many studies to be closely related to apoptosis. When LPS stimulates H9C2 cells (rat embryonic cardiomyoblasts), the overexpression of miR-146a suppresses apoptosis and helps attenuate myocardial depression (An et al., 2018). Chao et al. found that the overexpression of Tid1-S can enhance ER-α to activate p-PI3K/p-Akt, attenuating LPS-induced apoptosis in H9C2 cells (Chao et al., 2019). In recent years, as the understanding of m⁶A modification has gradually deepened, studies have found that m⁶A modification plays an important role in the process of apoptosis. In the human glioma cell line U251, lower level of m⁶A modification can reduce apoptosis and promote cell proliferation (Li et al., 2019). Moreover, METTL3 can regulate the expression of several important proteins that regulate the survival, apoptosis, invasion, and proliferation of lung cancer cells (Lin et al., 2016). We hypothesized that the m⁶A modification of Tnfrsf26 and Stk38l might be related to cardiomyocyte apoptosis.

LncRNAs are a class of noncoding RNAs that are more than 200 nucleotides in length. In recent years, m⁶A-modified lncRNAs have received extensive attention. M⁶A modification may control gene expression by regulating the translation efficiency and stability of lncRNAs (Coker et al., 2019). For example, the m⁶A-modified lncRNA MALAT1 is implicated in ischemia reperfusion injury-induced inflammation (Yang et al., 2020). LncRNAs contain miRNA-responsive elements, which function as ceRNAs, interact with miRNAs and indirectly regulate mRNAs. In the present study, we established and analyzed a potential ceRNA network of XR_346771 through bioinformatics analysis. We performed GO and KEGG analyses of the predicted mRNAs to further explore the function of XR_346771.

GO analysis showed that “inorganic cation import across plasma membrane”, “inorganic ion import across plasma membrane”, “transmembrane transport”, and “import across plasma membrane” were significantly enriched. Cardiac contraction requires Ca²⁺, so abnormal Ca²⁺ homeostasis may play a key role in the pathogenesis of common cardiovascular diseases. Studies have reported that m⁶A modification promotes neurite extension and neuronal differentiation by acting as a Ca²⁺ channel (Mukobata et al., 2002). In 6-OHDA-induced PC12 cells and the cerebral striatum of rats with Parkinson’s disease, decreased m⁶A modification induces the expression of N-methyl-d-aspartate receptor 1 and increases oxidative stress and Ca²⁺ influx, leading to apoptosis of dopaminergic neurons (Chen et al., 2019). These studies indicate that m⁶A modification is involved in Ca²⁺ transport. We speculated that the m⁶A modification of XR_346771 might associate with abnormal Ca²⁺ homeostasis in cardiac tissues.

It is interesting that even though the methyltransferase complex was downregulated upon sepsis induced myocardial dysfunction; there were still 27 transcripts that were significantly hypermethylated. Chokkalla et al. found that, compared with sham group, global m⁶A increased significantly at 12 and 24 h of reperfusion in transient middle cerebral artery occlusion in C57BL/6J mice. The FTO decreased significantly after stroke compared with sham. While 139 transcripts (122 mRNAs and 17 lncRNAs) were hypermethylated, 8 transcripts (5 mRNAs and 3 lncRNAs) were hypomethylated in the ischemic brain at 12 h reperfusion (Chokkalla et al., 2019). Wang et al. found that METTL3 was downregulated after traumatic brain injury in mice. In total, 922 m⁶A peaks were differentially expressed as determined by m⁶A-RIP-seq, with 370 upregulated and 552 downregulated methylated mRNA (Wang et al., 2019b). These studies are consistent with our research. Previous research indicated that the methylation modification of m⁶A has been proved to be reversible, as it modulated by RNA methyltransferases, RNA-binding proteins and demethylases. The m⁶A methyltransferase, also known as “Writers”, including METTL3/14, WTAP, KIAA1429 and RNA binding motifs protein 15/15B (RBM15/15B), E3 ubiquitin-protein ligase Hakai (HAKAI) and zinc finger CCCH-type containing 13 (ZC3H13). The RNA-binding proteins, also known as “readers”, including: YTHDF1/2/3, YTH domain-containing reader proteins 1/2 (YTHDC1/2), insulin-like growth factor 2 mRNA-binding proteins 1/2/3 (IGF2BP1/2/3), epithelium-specific ETS (ESE) transcription factors (ELF3), heterogeneous nuclear ribonucleoprotein A2B1 (hnRNPA2B1). The m⁶A demethylase, also known as “eraser”, including: FTO, α-ketoglutarate-dependent dioxygenase alkB homologue 5 (ALKHB5). In our research we tested mRNA level of some enzymes, but not all enzymes. Therefore, there may be other enzymes involved in m⁶A modification upon sepsis induced myocardial dysfunction. The m⁶A modification status of specific RNA was catalyzed by different enzymes. METTL3 positively regulates expression of MYD88, a critical upstream regulator of NF-κB signaling, by facilitating m⁶A methylation modification to MYD88-RNA in mesenchymal stem cells (Yu et al., 2020). However, reduced expression of FTO increased Nanog mRNA m⁶A methylation under TNF-α stimulation, decreased Nanog mRNA and protein expression, and significantly inhibited mesenchymal stem cells capacity for differentiation to sweat gland cells (Wang et al., 2020). In our research, these transcripts that have undergone m⁶A modification might have been modified by different enzymes.

In conclusion, our experiments revealed that the level of m⁶A modification is significantly decreased in septic cardiac tissue. Through the genome-wide profiling of m⁶A-tagged mRNAs and lncRNAs and subsequent bioinformatics analysis, we revealed some potential functions of transcripts with altered m⁶A modification.

References

• 1

  AnR.FengJ.XiC.XuJ.SunL. (2018). miR-146a Attenuates Sepsis-Induced Myocardial Dysfunction by Suppressing IRAK1 and TRAF6 via Targeting ErbB4 Expression. Oxid. Med. Cell Longev.2018, 1–9. 10.1155/2018/7163057

  • CrossRef
  • Google Scholar

• 2

  ArcidiaconoO. A.KrejčíJ.BártováE. (2020). The Distinct Function and Localization of METTL3/METTL14 and METTL16 Enzymes in Cardiomyocytes. Int. J. Mol. Sci.21, 8139. 10.3390/ijms21218139

  • CrossRef
  • Google Scholar

• 3

  BerulavaT.BuchholzE.ElerdashviliV.PenaT.IslamM. R.LbikD.et al (2020). Changes in m6A RNA Methylation Contribute to Heart Failure Progression by Modulating Translation. Eur. J. Heart Fail.22, 54–66. 10.1002/ejhf.1672

  • CrossRef
  • Google Scholar

• 4

  ChangY.FuX. R.CuiM.LiW. M.ZhangL.LiX.et al (2019). Activated Hippo Signal Pathway Inhibits Cell Proliferation and Promotes Apoptosis in NK/T Cell Lymphoma Cells. Cancer Med.8, 3892–3904. 10.1002/cam4.2174

  • CrossRef
  • Google Scholar

• 5

  ChaoC. N.LoJ. F.KhanF. B.DayC. H.LaiC. H.ChenC. H.et al (2019). Tid1‐S Attenuates LPS‐induced Cardiac Hypertrophy and Apoptosis through ER‐a Mediated Modulation of p‐PI3K/p‐Akt Signaling cascade. J. Cel Biochem.120, 16703–16710. 10.1002/jcb.28928

  • CrossRef
  • Google Scholar

• 6

  ChenX.YuC.GuoM.ZhengX.AliS.HuangH.et al (2019). Down-Regulation of m6A mRNA Methylation Is Involved in Dopaminergic Neuronal Death. ACS Chem. Neurosci.10, 2355–2363. 10.1021/acschemneuro.8b00657

  • CrossRef
  • Google Scholar

• 7

  ChokkallaA. K.MehtaS. L.KimT.ChelluboinaB.KimJ.VemugantiR. (2019). Transient Focal Ischemia Significantly Alters the M 6 A Epitranscriptomic Tagging of RNAs in the Brain. Stroke50, 2912–2921. 10.1161/strokeaha.119.026433

  • CrossRef
  • Google Scholar

• 8

  CohenJ.VincentJ.-L.AdhikariN. K. J.MachadoF. R.AngusD. C.CalandraT.et al (2015). Sepsis: a Roadmap for Future Research. Lancet Infect. Dis.15, 581–614. 10.1016/s1473-3099(15)70112-x

  • CrossRef
  • Google Scholar

• 9

  CokerH.WeiG.BrockdorffN. (2019). m6A Modification of Non-coding RNA and the Control of Mammalian Gene Expression. Biochim. Biophys. Acta Gene Regul. Mech.1862, 310–318. 10.1016/j.bbagrm.2018.12.002

  • CrossRef
  • Google Scholar

• 10

  EnrightA. J.JohnB.GaulU.TuschlT.SanderC.MarksD. S. (2003). MicroRNA Targets in Drosophila. Genome Biol.5, R1. 10.1186/gb-2003-5-1-r1

  • CrossRef
  • Google Scholar

• 11

  FengZ.LiQ.MengR.YiB.XuQ. (2018). METTL 3 Regulates Alternative Splicing of MyD88 upon the Lipopolysaccharide‐induced Inflammatory Response in Human Dental Pulp Cells. J. Cel Mol. Med.22, 2558–2568. 10.1111/jcmm.13491

  • CrossRef
  • Google Scholar

• 12

  HuangT.GuoJ.LvY.ZhengY.FengT.GaoQ.et al (2019). Meclofenamic Acid Represses Spermatogonial Proliferation through Modulating m6A RNA Modification. J. Anim. Sci. Biotechnol.10, 63. 10.1186/s40104-019-0361-6

  • CrossRef
  • Google Scholar

• 13

  ImanishiM.TsujiS.SudaA.FutakiS. (2017). Detection ofN6-Methyladenosine Based on the Methyl-Sensitivity of MazF RNA Endonuclease. Chem. Commun.53, 12930–12933. 10.1039/c7cc07699a

  • CrossRef
  • Google Scholar

• 14

  KennedyE. M.BogerdH. P.KornepatiA. V. R.KangD.GhoshalD.MarshallJ. B.et al (2016). Posttranscriptional M 6 A Editing of HIV-1 mRNAs Enhances Viral Gene Expression. Cell Host Microbe19, 675–685. 10.1016/j.chom.2016.04.002

  • CrossRef
  • Google Scholar

• 15

  KmietczykV.RiechertE.KalinskiL.BoileauE.MalovrhE.MaloneB.et al (2019). m6A-mRNA Methylation Regulates Cardiac Gene Expression and Cellular Growth. Life Sci. Alliance2, e201800233. 10.26508/lsa.201800233

  • CrossRef
  • Google Scholar

• 16

  LiF.ZhangC.ZhangG. (2019). m6A RNA Methylation Controls Proliferation of Human Glioma Cells by Influencing Cell Apoptosis. Cytogenet. Genome Res.159, 119–125. 10.1159/000499062

  • CrossRef
  • Google Scholar

• 17

  LinS.ChoeJ.DuP.TribouletR.GregoryR. I. (2016). The M 6 A Methyltransferase METTL3 Promotes Translation in Human Cancer Cells. Mol. Cel62, 335–345. 10.1016/j.molcel.2016.03.021

  • CrossRef
  • Google Scholar

• 18

  LinX.ChaiG.WuY.LiJ.ChenF.LiuJ.et al (2019). RNA m6A Methylation Regulates the Epithelial Mesenchymal Transition of Cancer Cells and Translation of Snail. Nat. Commun.10, 2065. 10.1038/s41467-019-09865-9

  • CrossRef
  • Google Scholar

• 19

  MathiyalaganP.AdamiakM.MayourianJ.SassiY.LiangY.AgarwalN.et al (2019). FTO-dependent N 6 -Methyladenosine Regulates Cardiac Function during Remodeling and Repair. Circulation139, 518–532. 10.1161/circulationaha.118.033794

  • CrossRef
  • Google Scholar

• 20

  MukobataS.HibinoT.SugiyamaA.UranoY.InatomiA.KanaiY.et al (2002). M6a Acts as a Nerve Growth Factor-Gated Ca2+ Channel in Neuronal Differentiation. Biochem. Biophys. Res. Commun.297, 722–728. 10.1016/s0006-291x(02)02284-2

  • CrossRef
  • Google Scholar

• 21

  NieM. W.HanY. C.ShenZ. J.XieH. Z. (2020). Identification of circRNA and mRNA Expression Profiles and Functional Networks of Vascular Tissue in Lipopolysaccharide‐induced Sepsis. J. Cel Mol. Med.24, 7915–7927. 10.1111/jcmm.15424

  • CrossRef
  • Google Scholar

• 22

  PanT. (2018). Modifications and Functional Genomics of Human Transfer RNA. Cell Res28, 395–404. 10.1038/s41422-018-0013-y

  • CrossRef
  • Google Scholar

• 23

  PasquinelliA. E. (2012). MicroRNAs and Their Targets: Recognition, Regulation and an Emerging Reciprocal Relationship. Nat. Rev. Genet.13, 271–282. 10.1038/nrg3162

  • CrossRef
  • Google Scholar

• 24

  RayesJ.LaxS.WichaiyoS.WatsonS. K.DiY.LombardS.et al (2017). The Podoplanin-CLEC-2 axis Inhibits Inflammation in Sepsis. Nat. Commun.8, 2239. 10.1038/s41467-017-02402-6

  • CrossRef
  • Google Scholar

• 25

  RubioR. M.DepledgeD. P.BiancoC.ThompsonL.MohrI. (2018). RNA m6A Modification Enzymes Shape Innate Responses to DNA by Regulating Interferon β. Genes Dev.32, 1472–1484. 10.1101/gad.319475.118

  • CrossRef
  • Google Scholar

• 26

  SempleJ. W.ItalianoJ. E.Jr.FreedmanJ. (2011). Platelets and the Immune Continuum. Nat. Rev. Immunol.11, 264–274. 10.1038/nri2956

  • CrossRef
  • Google Scholar

• 27

  SongH.FengX.ZhangH.LuoY.HuangJ.LinM.et al (2019). METTL3 and ALKBH5 Oppositely Regulate m6A Modification of TFEB mRNA, Which Dictates the Fate of Hypoxia/reoxygenation-Treated Cardiomyocytes. Autophagy15, 1419–1437. 10.1080/15548627.2019.1586246

  • CrossRef
  • Google Scholar

• 28

  WangJ.YanS.LuH.WangS.XuD. (2019a). METTL3 Attenuates LPS-Induced Inflammatory Response in Macrophages via NF-κB Signaling Pathway. Mediators Inflamm.2019, 1–8. 10.1155/2019/3120391

  • CrossRef
  • Google Scholar

• 29

  WangY.MaoJ.WangX.LinY.HouG.ZhuJ.et al (2019b). Genome-wide Screening of Altered m6A-Tagged Transcript Profiles in the hippocampus after Traumatic Brain Injury in Mice. Epigenomics11, 805–819. 10.2217/epi-2019-0002

  • CrossRef
  • Google Scholar

• 30

  WangY.WangR.YaoB.HuT.LiZ.LiuY.et al (2020). TNF-α Suppresses Sweat Gland Differentiation of MSCs by Reducing FTO-Mediated m6A-Demethylation of Nanog mRNA. Sci. China Life Sci.63, 80–91. 10.1007/s11427-019-9826-7

  • CrossRef
  • Google Scholar

• 31

  XiangB.ZhangG.GuoL.LiX.-A.MorrisA. J.DaughertyA.et al (2013). Platelets Protect from Septic Shock by Inhibiting Macrophage-dependent Inflammation via the Cyclooxygenase 1 Signalling Pathway. Nat. Commun.4, 2657. 10.1038/ncomms3657

  • CrossRef
  • Google Scholar

• 32

  YangC.FanZ.YangJ. (2020). m6A Modification of LncRNA MALAT1: A Novel Therapeutic Target for Myocardial Ischemia-Reperfusion Injury. Int. J. Cardiol.306, 162. 10.1016/j.ijcard.2019.11.140

  • CrossRef
  • Google Scholar

• 33

  YangN.ShiX.-L.ZhangB.-L.RongJ.ZhangT.-N.XuW.et al (2018). The Trend of β3-Adrenergic Receptor in the Development of Septic Myocardial Depression: A Lipopolysaccharide-Induced Rat Septic Shock Model. Cardiology139, 234–244. 10.1159/000487126

  • CrossRef
  • Google Scholar

• 34

  YuJ.ShenL.LiuY.MingH.ZhuX.ChuM.et al (2020). The m6A Methyltransferase METTL3 Cooperates with Demethylase ALKBH5 to Regulate Osteogenic Differentiation through NF-κB Signaling. Mol. Cel Biochem.463, 203–210. 10.1007/s11010-019-03641-5

  • CrossRef
  • Google Scholar

• 35

  ZakyA.DeemS.BendjelidK.TreggiariM. M. (2014). Characterization of Cardiac Dysfunction in Sepsis. Shock41, 12–24. 10.1097/shk.0000000000000065

  • CrossRef
  • Google Scholar

• 36

  ZeniF.FreemanB.NatansonC. (1997). Anti-inflammatory Therapies to Treat Sepsis and Septic Shock. Crit. Care Med.25, 1095–1100. 10.1097/00003246-199707000-00001

  • CrossRef
  • Google Scholar

• 37

  ZhangT.-N.YangN.GoodwinJ. E.MahrerK.LiD.XiaJ.et al (2019). Characterization of Circular RNA and microRNA Profiles in Septic Myocardial Depression: a Lipopolysaccharide-Induced Rat Septic Shock Model. Inflammation42, 1990–2002. 10.1007/s10753-019-01060-8

  • CrossRef
  • Google Scholar

• 38

  ZhouY.ZengP.LiY.-H.ZhangZ.CuiQ. (2016). SRAMP: Prediction of Mammalian N6-Methyladenosine (m6A) Sites Based on Sequence-Derived Features. Nucleic Acids Res.44, e91. 10.1093/nar/gkw104

  • CrossRef
  • Google Scholar