- Research article
- Open Access
Analyses of changes in myocardial long non-coding RNA and mRNA profiles after severe hemorrhagic shock and resuscitation via RNA sequencing in a rat model
BMC Molecular Biology volume 19, Article number: 11 (2018)
Ischemia–reperfusion injury has been proven to induce organ dysfunction and death, although the mechanism is not fully understood. Long non-coding RNAs (lncRNAs) have drawn wide attention with their important roles in the gene expression of some biological processes and diseases, including myocardial ischemia–reperfusion (I/R) injury. In this paper, a total of 26 Sprague–Dawley (SD) rats were randomized into two groups: sham and ischemia–reperfusion (I/R) injury. Hemorrhagic shock was induced by removing 45% of the estimated total blood volume followed by reinfusion of shed blood. High-throughput RNA sequencing was used to analyze differentially expressed (DE) lncRNAs and messenger RNAs (mRNAs) in the heart tissue 4 h after reperfusion. Myocardial function was also evaluated.
After resuscitation, the decline of myocardial function of shocked animals, expressed by cardiac output, ejection fraction, and myocardial performance index (MPI), was significant (p < 0.05). DE lncRNAs and mRNAs were identified by absolute value of fold change ≥ 2 and the false discovery rate ≤ 0.001. In rats from the I/R injury group, 851 lncRNAs and 1015 mRNAs were significantly up-regulated while 1533 lncRNAs and 1702 m RNAs were significantly down-regulated when compared to the sham group. Among the DE lncRNAs, we found 12 location-associated with some known apoptosis-related protein-coding genes which were up-regulated or down-regulated accordingly, including STAT3 and Il1r1. Real time PCR assays confirmed that the expression levels of five location-associated lncRNAs (NONRATT006032.2, NONRATT006033.2, NONRATT006034.2, NONRATT006035.2 and NONRATT029969.2) and their location-associated mRNAs (STAT3 and Il1r1) in the rats from the I/R injury group were all significantly up-regulated versus the sham group.
The DE lncRNAs (NONRATT006032.2, NONRATT006033.2, NONRATT006034.2 and NONRATT006035.2) could be compatible with their role in myocardial protection by stimulating their co-located gene (STAT3) after hemorrhagic shock and resuscitation. The final prognosis of I/R injury might be regulated by different genes, which is regarded as a complex network.
Severe hemorrhagic shock, which is a major cause of morbidity and mortality for patients suffering from trauma, visceral hemorrhage, or major surgery, can lead to myocardial dysfunction [1,2,3,4]. The immediate presence of adequate perfusion of the tissue (reperfusion) in patients of hemorrhagic shock is the fundamental goal of acute resuscitation, however it may result in an injury that potentially induces organ dysfunction and death, known as ischemia–reperfusion (I/R) injury . Suggested mechanisms of I/R injury include the generation of reactive oxygen species, intracellular Ca2+ overload, the increase of proinflammatory cytokines, and mitochondrial dysfunction, which subsequently results in apoptosis [6,7,8]. Although defining its complex pathophysiological process remains a challenge, ascertaining how to relieve I/R injury is of great importance. To that end, thorough knowledge of the genetic basis of I/R injury may provide therapeutic targets and tools.
Long non-coding RNAs (lncRNAs) are a family of RNA molecules (containing more than 200 nucleotides) which do not encode proteins, but functionally participate in many processes controlling gene expression and cell differentiation . In recent years, lncRNAs have drawn wide attention with their critical roles in gene expression of various biological processes and diseases, including Parkinson’s disease , liver disease , gynecological diseases  and tumors [13, 14]. They serve as biomarkers, prognostic assessments, and even therapeutic targets. Wang et al.  reported that an lncRNA named necrosis-related factor (NRF) reduced myocardial infarct size upon I/R injury in the animal model, with miR-873 as one of its targets suppressing the expression of RIPK1 (receptor-interacting serine/threonine-protein kinase 1)/RIPK3 (receptor-interacting serine/threonine-protein kinase 3) and decreasing death of cardiac cells. However, the operating function of lncRNAs may be a complex network involving different facets of cell biology. More lncRNAs and their effects in the mechanism of I/R injury remain to be identified.
In our previous study , we preliminarily confirmed the correlation between the phosphorylation of protein STAT3 and the myocardial protective mechanism after myocardial I/R injury. In this study, we identified changes in the expressions of mRNAs and lncRNAs in myocardia after severe hemorrhagic shock and resuscitation models compared with sham rats using RNA sequencing techniques, aiming to explore the upstream regulator of protein STAT3 for cardioprotective effects. Our findings predict the regulatory role of lncRNAs. It could lead to further research to explore therapeutic targets for cardioprotection in I/R injuries induced by hemorrhagic shock and resuscitation.
Model identification of I/R injury caused by hemorrhagic shock and resuscitation
Between the I/R injury and sham groups, there was no difference in the baseline characteristics (Table 1). In the hemorrhage and shock phases, mean arterial pressure (MAP) of the I/R injury group was significantly lower compared to the sham group (p < 0.05). In the resuscitation phase, MAP of the I/R injury group rapidly returned to normal, nearly matching the baseline (Fig. 1a). After resuscitation, the CO and EF levels of the I/R injury group had significantly decreased compared to the sham group and baseline value (p < 0.05). Meanwhile, the MPI of the I/R injury group showed a significant raise compared to the sham group since the resuscitation phase (p < 0.05, Fig. 1b).
DE lncRNAs and mRNAs
We examined all significantly different expression levels of the transcripts in the heart tissue of I/R injury rats as previously described . Differentially expressed lncRNAs and mRNAs were identified by the absolute value of fold change bigger than 2 and the FDR < 0.001. The total reads (mRNAs and lncRNAs) of the two samples were 163, 128, 480 and 84, 501, 062, respectively, while the clean bases of the two samples were 24.47G and 12.68G. The raw data in this study are available in the NCBI SRA database.
Our data revealed that 851 lncRNAs were up-regulated in rats from the I/R injury group, and 1533 lncRNAs were down-regulated when compared to the sham group (Fig. 2a). Meanwhile, 1015 mRNAs were significantly up-regulated and 1702 m RNAs were significantly down-regulated versus the sham group (Fig. 2b). The top 5 up-regulated and down-regulated lncRNAs and mRNAs are listed in Tables 2 and 3.
DE lncRNA target mRNAs
LncRNAs have been shown to function as regulators of adjacent protein-coding genes through RNA–protein interactions [18, 19]. In our study, lncRNAs located within 5 kb of known protein-coding genes, or those contained in the introns of genes, were considered to be location-associated lncRNAs. Among the DE lncRNAs, we found 12 of them location-associated with some apoptosis—related protein-coding genes (Table 4). For instance, four lncRNAs named NONRATT006032.2, NONRATT006033.2, NONRATT006034.2, and NONRATT006035.2 were location-associated with STAT3; NONRATT029969.2 was location-associated with Il1r1. STAT3 codes STAT3 proteins and our previous research showed they were involved in the mechanism of cardio-protection. Consistent with the lncRNAs, expressions of the location-associated mRNAs were up-regulated or down-regulated in I/R injury group. On these grounds we affirm that lncRNAs could participate in expressions of co-located protein-coding genes. The sequence read data of the research can be obtained from the NCBI database with the accession code: SRP164364.
Real-time quantitative PCR
To validate the accuracy of lncRNA and mRNA expression profiles determined by next generation RNA sequencing, several specific location-associated lncRNAs and mRNAs in the heart segments of intact rats were confirmed by real-time PCR assay. The expression levels of five location-associated lncRNAs in the I/R injury group, NONRATT006032.2, NONRATT006033.2, NONRATT006034.2, NONRATT006035.2 and NONRATT029969.2, were all significantly up-regulated versus the sham group (Fig. 3a–e). The expression levels of two location-associated mRNAs in the I/R injury group, STAT3 and Il1r1, were also significantly enhanced compared to the sham group (Fig. 3f, g).
Several studies have highlighted that I/R injury, ultimately leading to heart failure, is a major cause of morbidity and mortality [20, 21]. However, the molecular mechanism is not fully understood. In the present study, which is the first of its kind, we utilized next-generation RNA sequencing technology and bioinformatics analysis to investigate whole mRNA and lncRNA profiles of heart tissue in rat models suffering from I/R injury caused by severe hemorrhagic shock. The results of our research indicated differential expression of mRNAs and lncRNAs in the cardiomyocytes between the I/R injury group and sham group after hemorrhagic shock and resuscitation. In the I/R injury group 851 up-regulated and 1533 down-regulated lncRNAs were deemed to be involved in the I/R injury. Meanwhile, 1015 up-regulated transcriptions and 1702 down-regulated transcriptions in the I/R injury group were found to be DE mRNAs. To further validate the reliability of the sequencing results, expressions of several DE lncRNAs and mRNAs were verified by qRT-PCR.
A family of RNA, which had long been considered as a platform for protein production until recently, was reported to play a major role in pathophysiology of human diseases including complex heart diseases [7, 22, 23]. Meanwhile, their potential roles in I/R injuries were also proven [20, 21]. However, myocardial lncRNA expression and how it functions in hemorrhagic shock induced I/R injuries haven’t yet been fully explored. Liu et al.  investigated lncRNA expression in the infarct region of mice after myocardial ischemia by RNA-sequencing analysis. They reported 31,423 arrayed lncRNAs, of which 151 were aberrantly expressed (64 up-regulated and 87 down-regulated). In contrast, our research found a remarkably greater number of dysregulated lncRNAs. There are several possible reasons for the divergence. Firstly, different species of animal were used in the two studies, and the gene expression varies between species. Secondly, we analyzed DE profiles by high throughput RNA-sequencing techniques rather than gene microarray. The former could detect a great quantity of novel transcriptions more comprehensively . Thirdly, our research was based on I/R models induced by bleeding (removing almost half estimated total blood volume) followed by resuscitation of shed blood, and the diffusely ischemic heart tissue was used instead of the focal infarct sections. The Liu’s I/R models were performed via occlusion followed by reperfusion of the left anterior coronary artery, which led to myocardial infarction. Our animal models simulated more physiopathology processes found in severe hemorrhagic shock. The significant difference of the DE lncRNA numbers between the two studies indicated that aberrant expression of lncRNA may has tissue specificity. Finally, the two studies screened DE lncRNAs at different time points (2 h and 4 h after I/R injury respectively), in which durations of ischemia were also different (45 min and 60 min respectively). This discrepancy uncovered that the expression of lncRNAs involved in I/R injury might be correlated with the degree of injury. Moreover, how to select the optimal timing for treatment of I/R injury is also an interesting topic, which is absolutely worth further study.
The balance between cell survival and death is generally regarded as the prime determinant of prognosis for ischemic injury . As a main type of programmed cell death, apoptosis is proven to be crucial to an I/R injury [27, 28], which is regulated by several signaling pathways including the janus kinase (JAK)/signal transducer and activator of transcription (STAT) pathway. Evidence has demonstrated that the activation of STAT3 limits apoptosis in rat models of myocardial infarction , while others have highlighted the importance of JAK2-STAT3 activation in apoptosis of post-ischemic damage .
The JAK/STAT pathway is activated after the binding of a ligand to the receptor in the plasma membrane. The ligand phosphorylates the receptor by activating JAK kinases then binds cytosolic STAT proteins to it. The STAT proteins are phosphorylated and translocated into the nucleus to regulate the transcription of target DNA sequences [31, 32]. Actually about 10 years ago, Oshima  reported that cardiac-specific transgenic mice expressing constitutively active STAT3 exhibited resistance against I/R injury, compared with non-transgenic littermates. The data of Xie et al.  revealed that elevated expression of JAK2/STAT3 mRNA and an activated JAK2/STAT3 signaling pathway acted as a cardioprotective factor in I/R injury mouse models. Our group recently reported the existence of a cardioprotective process called remote ischemic post-conditioning (RIPostC) in rats : the RIPostC group demonstrated attenuated myocardial injuries, increased survival rate and enhanced phosphorylation of STAT3 4 h after reperfusion. We verified that RIPostC protected cardiomyocytes partially by phosphorylating STAT3 through the SAFE pathway. This work result was consistent with the hypothesis that lncRNAs function as an upstream contributor of the JAK/STAT signaling pathway to confer cardio-protection in an I/R injury.
In the present study, we identified increased levels of 4 relevant lncRNAs (NONRATT006032.2, NONRATT006033.2, NONRATT006034.2 and NONRATT006035.2) in myocardial cells after an I/R injury. Meanwhile, a significant increase in the abundance of STAT3, predicted to be lncRNA co-located mRNA, was found. It is reported that the positions of about 65% lncRNAs are within 10 kb of known protein-coding genes, which may be involved in cis-acting or trans-regulatory relationships with their neighbor protein-coding genes [18, 19]. This means lncRNAs may modulate the expression patterns of these location-associated genes. We concluded that the aforementioned 4 enhanced lncRNAs up-regulated the expression of location-associated mRNA STAT3. They might proceed to stimulate the expression of STAT3 protein which was verified to relate to the anti-apoptotic mechanism after I/R injury.
However, the I/R injury group in our research showed a significant decrease rather than improvement in myocardial function despite cardio-protection of STAT3 during ischemia and reperfusion. Not surprisingly, we expect that other of the DE lncRNAs involved in pathogenesis of I/R injury may conspire to cause the decline of myocardial function. For instance, we confirmed that NONRATT029969.2 co-located with and significantly up-regulated its adjacent protein-coding gene IL1R1 by both high-throughput sequencing and qPCR. Stimulation of the interleukin 1 receptor (IL-1R) leads to the recruitment of adaptor proteins such as myeloid differentiation factor 88 (MyD88), which triggers large-scale downstream signaling cascades and production of proinflammatory cytokines [35, 36]. IL-1R antagonist was reported to have protective effects against I/R injury by attenuating the inflammatory response, which was associated with decreased apoptosis [37, 38].
Our analyses may indicate the complex network relationships among lncRNAs, mRNAs and involved proteins after I/R injury caused by severe hemorrhagic shock. Moreover, certain lncRNAs are likely to modulate their mRNA targets, not only neighbor protein-coding genes, by interacting with miRNAs and inhibiting the ability of them. This suggested a novel perspective to screen and predict target mRNAs involved in physiopathology mechanisms of I/R injury in our future research.
There are several limitations in our study. Firstly, due to the limited funds and resources, we chose to validate only 5 lncRNAs and 2 mRNAs by PCR, which might be involved in our previous study results, remaining other DE mRNAs to be analyzed in our further research. Secondly, knockdown of the 4 lncRNAs (NONRATT006032.2, NONRATT006033.2, NONRATT006034.2 and NONRATT006035.2) was not performed. Whether the lncRNA knockdown would increase cell apoptosis and reverse cardioprotective effects after hemorrhagic shock remains unknown. Thirdly, in this study we tested only three samples from the I/R injury and sham groups for RNA sequencing. More samples will be needed to verify the results in future research. Finally, the number of time points for sample collecting was relatively limited. Further studies should detect the lncRNA expression levels at more time points to determine if the regulation effect of lncRNAs is associated with length of ischemia and reperfusion time.
Our research provides evidence of location-associated modulation roles of lncRNAs (NONRATT006032.2, NONRATT006033.2, NONRATT006034.2 and NONRATT006035.2) in I/R injury rat models caused by hemorrhagic shock by targeting their co-located gene (STAT3). We predict that they may confer myocardial protection by stimulating expressions of downstream proteins (STAT3 protein) and reducing apoptosis through activated JAK2/STAT3 pathway. However, the final prognosis of I/R injury could be regulated by different genes, which is regarded as a complex network.
Experiments were performed on 26 male Sprague–Dawley rats weighing between 350 and 450 g (the Experimental Animal Center of Sun Yat-sen University, Guangzhou, China). The rats were kept in a specific sterile laboratory, at a uniform temperature of 20–22 °C and on a light/dark cycle of 12 h. Our research was approved by the Institutional Animal Care and Use Committee, established in the Tang Wanchun Laboratories of Emergency Critical Care Medicine (Sun Yat-sen Memorial Hospital, Sun Yat-sen University, IACUC-2016R1604).
After fasting overnight, with the exception of free access to water, the rats were anesthetized by intraperitoneal injection of sodium pentobarbital (45 mg/kg) and then stretched supinely on a surgical board. To maintain anesthesia, an additional administration of sodium pentobarbital (10 mg/kg) was required. The trachea was orally intubated with a 14-G cannula (Abbocath-T; Abbott Hospital Products Division, North Chicago, IL) mounted on a blunt needle with a 145° angled tip. End-tidal carbon dioxide pressure (ETCO2) was measured by the capnometer module of a BeneView T5 patient monitor (Mindray, Shenzhen, China), and a conventional lead II electrocardiogram was continuously monitored. Spontaneous respiration of room air was maintained. The rectal temperature was kept at 37 °C ± 0.5 °C with infrared surface heating lamps.
The left femoral artery was surgically exposed, into which a polyethylene catheter (PE-50; Becton–Dickinson, Franklin Lakes, NJ) was advanced to the descending aorta to measure aortic pressure and blood gas. Another PE-50 catheter was advanced into the right atrium through the left external jugular vein to measure the right atrial pressure. Saline with 2.5 IU/ml of crystalline bovine heparin was flushed intermittently to all catheters.
After surgical preparation, the animals were randomized into two groups: I/R injury and sham group. In the I/R injury group the animals were heparinized with 100 U/kg of bovine heparin. An estimated 45% of total blood volume, according to the blood volume of each animal (EBV = 6.12 ml/100 g BW) , was withdrawn over an interval of 60 min. An infusion/withdrawal dual syringe pump (LongerPump LSP01-1A, Longer Corporation, China) was used for bleeding and recovery of shed blood. Blood from the left femoral artery was allowed to flow into a sterile 20 ml syringe. The reinfusion of shed blood was performed 40 min after the completion of bleeding, over the ensuing 40 min. In the sham group, the procedures of animal preparation rather than bleeding were performed. All the animals were then monitored for an additional 4 h. ETCO2, electrocardiography, temperature, and aortic pressures were recorded at intervals of 5-min during hemorrhage and 15-min thereafter. Myocardial functions were measured at baseline and hourly after the ending of reinfusion. After the 4 h period, the animals were euthanized by intraperitoneal injection of pentobarbital (150 mg/kg) for the collection of heart tissue. Three rats from both groups were randomly selected for RNA sequencing while others remained for quantitative real-time RT-PCR evaluation.
We used a PC-based data-acquisition system supported by WINDAQ software (DATAQ, Akron, OH) to make continuous records of ETCO2 levels, electrocardiogram output, and aortic and right atrial pressures. Myocardial function was non-invasively measured with an Ultrasonix ultrasound system (Model130-4311, Ultrasonix Medical Corpration, Canada) by two separate units, covering cardiac output (CO), ejection fraction (EF), and myocardial performance index (MPI). CO and EF served as indexes of myocardial contractile function and MPI as a marker of both systolic and diastolic function . Aortic blood pH, PCO2, PO2, hemoglobin, and lactate concentrations were measured in 200 μL aliquots of blood with a Stat Profile pHOx Plus L analyzer (Model RADIOMETER ABL80FLEX; Radiometer Medical ApS, Bronshoj, Denmark).
Tissue collection and RNA isolation
Three rats from each group were randomly selected to collect heart segments. Total RNA was extracted from the heart tissue by Trizol Reagent (TR118-500, MRC) as shown in the manufacturer’s instructions. Quantity and quality of RNA were assessed by a Nano-100 (AoSheng Biotech Co. Ltd.) while RNA integrity was confirmed by 2100 bioanalyzer (Agilent).
Library construction and sequencing
We prepared total RNAs (1 μg per sample) as input material, removed ribosomal RNAs and constructed sequencing libraries using a VAHTS Total RNA-seq (H/M/R) Library Prep Kit from Illumina® (NR603-02, Vazyme) according to the manufacturer’s recommendation. The first chain of cDNA was generated using random hexamer primers and Reverse Transcriptase. The second chain was synthesized using RNase H and DNA pol I. dNTPs with dTTP were replaced by dUTP and the reaction buffer. The 3′ ends of DNA fragments were adenylated, and an Illumina sequencing Adapter was ligated for hybridization. VAHTS DNA Clean Beads (N411-03, Vazyme) were used to purify the library fragments to select cDNA fragments of 350–400 bp. Then PCR amplification was performed to establish the complete sequence of the cDNA library. Library quality was assessed on a 2100 bioanalyzer (Agilent). The library was sequenced on a Nextseq500 platform (Illumina Company, USA), and 150 bp paired-end reads were generated. We analyzed the sequencing results and dynamically removed joint sequence fragments and low-quality segments from the 3′ end. FastQC software was used for quality control of the raw data.
Comparison with reference genome
Reference genome sequencing alignment was performed against the whole genome sequence using TopHat software. The reference genome was Rattus norvegicus (rn6) (ftp://ftp.ensembl.org/pub/release-81/fasta/rattus_norvegicus/dna/). The known transcript information in the genome location was downloaded from Ensembl (ftp://ftp.ensembl.org/pub/release-81/gtf/rattus_norvegicus/Rattus_norvegicus.Rnor_6.0.81.abinitio.gtf.gz).
There followed a list of required parameters: mismatch of 2 bases was allowed; during each read optimal match record up to 20 bp was allowed; given the variable shear, fragment length of 25 bp was used; in each fragment mismatch of less than 2 bp was allowed; insertions and deletions were limited to 3 bp; to completely guarantee comparison, alternative splicing position mismatch was limited to 0; the allowed length of intron was 50–50,000 bp.
Assembly of transcripts
The mapped reads of each sample were assembled by Cufflinks . Considering read coverage gaps in the regions of the reference gene,we compensated the incompletely-assembled transcripts by performing reference annotation based transcripts (RABT) assembly with the reference gene annotation. The features could be missing in the sequencing data, attributed to low coverage. To avoid that, faux-reads were generated from reference transcripts. These reads and the (aligned) sequenced reads were merged for assembly.
Quantitative real-time PCR
We performed quantitative analysis by quantitative real-time PCR to confirm the expression levels of several DE lncRNAs and their target mRNAs. Total RNA reverse transcription was performed using M-MLV Reverse Transcriptase (Promega, M1705) following the manufacturer’s instructions with either orligo (dT) primers or specific RT primers. cDNA was quantified in an Opticon real-time PCR machine (MJ Research, Waltham, MA, USA). The primers used are listed in Table 5. Each sample was run in triplicate in 20 μl reactions with 0.4 μM forward and reverse PCR primers and 10 μl of GoTaq® qPCR Master Mix (Promega, A6002)in a ABI StepOne Plus real-time PCR system. The PCR cycle parameters were set as follows: initial denaturation at 95 °C for 2 min followed by 40 cycles of 95 °C for 15 s and 60 °C for 30 s. Relative expression was determined by normalization to GAPDH using the 2−ΔΔCt method. The experiments were performed in triplicate.
Normally distributed data were presented as mean ± SD while non-normally distributed data were presented as a median (25th, 75th percentiles). Statistical analysis was performed using a t-test between two cohorts. RT-PCR results were performed with one-way analysis of variance (ANOVA) followed by Tukey’s multiple comparison test. Data analysis and subsequent quantitative normalization were performed using GeneSpring GX v12.1 software (Agilent Technologies). Differentially expressed lncRNAs and mRNAs were identified by absolute value of fold change ≥ 2 and FDR ≤ 0.001. Statistical significance was set at a p < 0.05.
long non-coding RNAs
myocardial performance index
false discovery rate
mean arterial pressure
Aimo A, Borrelli C, Giannoni A, Pastormerlo LE, Barison A, Mirizzi G, et al. Cardioprotection by remote ischemic conditioning: mechanisms and clinical evidences. World J Cardiol. 2015;7:621.
Merchant RM, Yang L, Becker LB, Berg RA, Nadkarni V, Nichol G, et al. Incidence of treated cardiac arrest in hospitalized patients in the United States. Crit Care Med. 2011;39:2401–6.
Robinson S, Swain AH, Hoyle SR, Larsen PD. Survival from out-of-hospital cardiac arrest in New Zealand following the 2005 resuscitation guideline changes. Resuscitation. 2010;81:1648.
Wutzler S, Maegele M, Wafaisade A, Wyen H, Marzi I, Lefering R. Risk stratification in trauma and haemorrhagic shock: scoring systems derived from the TraumaRegister DGU ®. Injury. 2014;45:S29–34.
Eltzschig HK, Eckle T. Ischemia and reperfusion—from mechanism to translation. Nat Med. 2011;17:1391–401.
Garcia-Dorado D, Rodriguez-Sinovas A, Ruiz-Meana M, Inserte J. Protection against myocardial ischemia–reperfusion injury in clinical practice. Rev Esp Cardiol (Engl Ed). 2014;67:394–404.
Yang Y, Zhao L, Ma J. Penehyclidine hydrochloride preconditioning provides cardiac protection in a rat model of myocardial ischemia/reperfusion injury via the mechanism of mitochondrial dynamics mechanism. Eur J Pharmacol. 2017;813:130–9.
Wang Z, Lin D, Zhang L, Liu W, Tan H, Ma J. Penehyclidine hydrochloride prevents anoxia/reoxygenation injury and induces H9c2 cardiomyocyte apoptosis via a mitochondrial pathway. Eur J Pharmacol. 2017;797:115–23.
Splinter E, de Wit E, Nora EP, Klous P, van de Werken HJ, Zhu Y, et al. The inactive X chromosome adopts a unique three-dimensional conformation that is dependent on Xist RNA. Genes Dev. 2011;25:1371–83.
Liu W, Zhang Q, Zhang J, Pan W, Zhao J, Xu Y. Long non-coding RNA MALAT1 contributes to cell apoptosis by sponging miR-124 in Parkinson disease. Cell Biosci. 2017;7:19.
Chen Y, Huang H, Xu C, Yu C, Li Y. Long non-coding RNA profiling in a non-alcoholic fatty liver disease rodent model: new insight into pathogenesis. Int J Mol Sci. 2017;18:21.
Wang WT, Sun YM, Huang W, He B, Zhao YN, Chen YQ. Genome-wide long non-coding RNA analysis identified circulating LncRNAs as novel non-invasive diagnostic biomarkers for gynecological disease. Sci Rep. 2016;6:23343.
Wang F, Xie C, Zhao W, Deng Z, Yang H, Fang Q. Long non-coding RNA CARLo-5 expression is associated with disease progression and predicts outcome in hepatocellular carcinoma patients. Clin Exp Med. 2017;17:33–43.
Ma B, Li M, Zhang L, Huang M, Lei JB, Fu GH, et al. Upregulation of long non-coding RNA TUG1 correlates with poor prognosis and disease status in osteosarcoma. Tumour Biol. 2016;37:4445–55.
Wang K, Liu F, Liu CY, An T, Zhang J, Zhou LY, et al. The long noncoding RNA NRF regulates programmed necrosis and myocardial injury during ischemia and reperfusion by targeting miR-873. Cell Death Differ. 2016;23:1394–405.
Huang J, Xu D, Guo Q, Ou B, Ling Q, Li J, et al. Remote ischemic post-conditioning improves myocardial dysfunction via the risk and safe pathways in a rat model of severe hemorrhagic shock. Shock. 2017;49(4):460–5.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Wang KC, Yang YW, Liu B, Sanyal A, Corces-Zimmerman R, Chen Y, et al. A long noncoding RNA maintains active chromatin to coordinate homeotic gene expression. Nature. 2011;472:120–4.
Engreitz JM, Haines JE, Perez EM, Munson G, Chen J, Kane M, et al. Local regulation of gene expression by lncRNA promoters, transcription and splicing. Nature. 2016;539:452–5.
Hausenloy DJ, Yellon DM. Myocardial ischemia–reperfusion injury: a neglected therapeutic target. J Clin Invest. 2013;123:92–100.
Li X, Dai Y, Yan S, Shi Y, Han B, Li J, et al. Down-regulation of lncRNA KCNQ1OT1 protects against myocardial ischemia/reperfusion injury following acute myocardial infarction. Biochem Biophys Res Commun. 2017;491:1026–33.
Congrains A, Kamide K, Katsuya T, Yasuda O, Oguro R, Yamamoto K, et al. CVD-associated non-coding RNA, ANRIL, modulates expression of atherogenic pathways in VSMC. Biochem Biophys Res Commun. 2012;419:612–6.
Schonrock N, Harvey RP, Mattick JS. Long noncoding RNAs in cardiac development and pathophysiology. Circ Res. 2012;111:1349–62.
Liu Y, Li G, Lu H, Li W, Li X, Liu H, et al. Expression profiling and ontology analysis of long noncoding RNAs in post-ischemic heart and their implied roles in ischemia/reperfusion injury. Gene. 2014;543:15–21.
LaCroix-Fralish ML, Austin JS, Zheng FY, Levitin DJ, Mogil JS. Patterns of pain: meta-analysis of microarray studies of pain. Pain. 2011;152:1888–98.
Koike M, Shibata M, Tadakoshi M, Gotoh K, Komatsu M, Waguri S, et al. Inhibition of autophagy prevents hippocampal pyramidal neuron death after hypoxic-ischemic injury. Am J Pathol. 2008;172:454–69.
Lu SF, Huang Y, Wang N, Shen WX, Fu SP, Li Q, et al. Cardioprotective effect of electroacupuncture pretreatment on myocardial ischemia/reperfusion injury via antiapoptotic signaling. Evid Based Complement Alternat Med. 2016;2016:4609784.
Badalzadeh R, Mokhtari B, Yavari R. Contribution of apoptosis in myocardial reperfusion injury and loss of cardioprotection in diabetes mellitus. J Physiol Sci. 2015;65:201–15.
Tian Y, Zhang W, Xia D, Modi P, Liang D, Wei M. Postconditioning inhibits myocardial apoptosis during prolonged reperfusion via a JAK2-STAT3-Bcl-2 pathway. J Biomed Sci. 2011;18:53.
Zhang S, Liu X, Goldstein S, Li Y, Ge J, He B, et al. Role of the JAK/STAT signaling pathway in the pathogenesis of acute myocardial infarction in rats and its effect on NF-kappaB expression. Mol Med Rep. 2013;7:93–8.
Fischer P, Hilfiker-Kleiner D. Survival pathways in hypertrophy and heart failure: the gp130-STAT3 axis. Basic Res Cardiol. 2007;102:279–97.
O’Shea JJ, Gadina M, Schreiber RD. Cytokine signaling in 2002: new surprises in the Jak/Stat pathway. Cell. 2002;109(Suppl):S121–31.
Oshima Y, Fujio Y, Nakanishi T, Itoh N, Yamamoto Y, Negoro S, et al. STAT3 mediates cardioprotection against ischemia/reperfusion injury through metallothionein induction in the heart. Cardiovasc Res. 2005;65:428–35.
Xie XJ, Fan DM, Xi K, Chen YW, Qi PW, Li QH, et al. Suppression of microRNA-135b-5p protects against myocardial ischemia/reperfusion injury by activating JAK2/STAT3 signaling pathway in mice during sevoflurane anesthesia. Biosci Rep. 2017. https://doi.org/10.1042/BSR20170186.
Akira S, Uematsu S, Takeuchi O. Pathogen recognition and innate immunity. Cell. 2006;124:783–801.
Na Y, Bai R, Zhao Z, Wei Y, Li D, Wang Y, et al. 36、IL1R1 gene polymorphisms are associated with knee osteoarthritis risk in the Chinese Han population. Oncotarget. 2017;8:4228–33.
Meimei H, Dejin M, Erzhen C, Minmin S, Songyao J, Jianfang L, et al. Inhibiting the toll-like receptor 4 Toll/interleukin-1 receptor domain protects against hepatic warm ischemia and reperfusion injury in mice. Crit Care Med. 2014;42:e123–31.
Grothusen C, Hagemann A, Attmann T, Braesen J, Broch O, Cremer J, et al. Impact of an interleukin-1 receptor antagonist and erythropoietin on experimental myocardial ischemia/reperfusion injury. Sci World J. 2012;2012:737585.
Steensrud T, Li J, Dai X, Manlhiot C, Kharbanda RK, Tropak M, et al. Pretreatment with the nitric oxide donor SNAP or nerve transection blocks humoral preconditioning by remote limb ischemia or intra-arterial adenosine. Am J Physiol Heart Circ Physiol. 2010;299:H1598–603.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.
LL and ZY wrote this manuscript, and also contributed in analyzing and interpreting the data regarding this manuscript. GZ, YZ, YW and JL performed experiments, and was a major contributor in writing the manuscript. RC and WT were responsible for the design. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
Our supporting raw data are contained in the manuscript and in Additional file 1: Table S1, Additional file 2: Table S2, Additional file 3: Table S3, Additional file 4: Table S4 and Additional file 5: Table S5. The sequence read data of the research can be obtained from the NCBI database with the accession code: SRP164364. The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
Our research was approved by the Institutional Animal Care and Use Committee, established in the Tang Wanchun Laboratories of Emergency Critical Care Medicine (Sun Yat-sen Memorial Hospital, Sun Yat-sen University, IACUC-2016R1604).
The design of this study and the collection, analysis, and interpretation of data were supported by research grants from a project of Leading Talents in Pearl River Talent Plan of Guangdong Province (No. 81000-42020004) and Guangzhou Science and Technology Plan (No. 201804010471).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The raw data of mean arterial pressure, myocardial function and QPCR validation of sham group and I/R group.
The different expression levels of all lncRNAs between the sham group and I/R group.
The raw data of differentially expressed lncRNAs.
The different expression levels of all mRNAs between the sham group and I/R group.
The raw data of differentially expressed mRNAs.
About this article
Cite this article
Lin, L., Yang, Z., Zheng, G. et al. Analyses of changes in myocardial long non-coding RNA and mRNA profiles after severe hemorrhagic shock and resuscitation via RNA sequencing in a rat model. BMC Molecular Biol 19, 11 (2018). https://doi.org/10.1186/s12867-018-0113-8
- Hemorrhagic shock
- Ischemia reperfusion injury
- Long non-coding RNA
- High-throughput RNA sequencing