Circular RNA Microarray Analyses in Hepatic Ischemia-Reperfusion Injury With Ischemic Preconditioning Prevention

Abstract

Ischemic preconditioning (IPC) represents an effective intervention to relieve hepatic ischemia-reperfusion injury (IRI). Systematic detection of circRNA expression revealing the protection effect of IPC still remains to be elucidated. Here, we applied a microarray to detect circRNA and mRNA expression in ischemic liver with and without IPC (n = 3 in each group). Compared with the sham group, there were 39 circRNAs and 432 mRNAs increased and 38 circRNAs and 254 mRNAs decreased (fold change ≥1.5, P < 0.05) in the group of hepatic IRI. As the result of IPC intervention, 43 circRNAs and 64 mRNAs were increased, and 7 circRNAs and 31 mRNAs were decreased in the IPC group when compared with IRI. We then identified circRNA_017753 as the most possible target that may closely relate to IPC protective signaling and predicted Jade1 as the target related to circRNA_017753. Three possible circRNA–miRNA–mRNA axes were constructed that may play a vital role in protective mechanisms in IPC. The study for the first time systematically detects the dysregulated circRNAs and mRNAs in response to hepatic IRI and IPC intervention. Our profile and bioinformatic analysis provide numerous novel clues to understanding the pathophysiologic mechanism of IPC protection against hepatic IRI.

Introduction

Hepatic ischemia-reperfusion injury (IRI) occurs in clinical circumstances, including hepatic resection, transplantation, liver trauma, or septic shock (1). The reperfusion aggravates hepatic injury after ischemia. Concerning mechanisms involve microcirculatory failure, inflammatory cytokine release, and reactive oxygen species accumulation (2, 3). Especially in the surgical procedures of hepatic resections and liver transplantation, IRI not only contributes to organ damage, but also reduces the long-term survival rates. Therefore, strategies to reduce hepatic IRI and improve patient outcomes is clinically important at any point.

Ischemic preconditioning (IPC) refers to an intrinsic procedure in the form of repeated short episodes of ischemia that increase the resistance of organs against IRI. This powerful intervention for protection against IRI was first reported in 1986 by Murry et al. (4). Though, today, the surgical techniques and clinical conditions improve rapidly concerning the intervention of IPC, the precise molecular mechanisms behind IRI and the protection effect of IPC still remain an important problem. Moreover, recent studies comparing different protective methods of IRI have demonstrated variable protective mechanisms in IPC and ischemic postconditioning (IPostC) (5, 6). The mechanism of comparative and collaborative research on IPC and IPostC are of great significance in the research on IRI.

Nowadays, high-throughput RNA sequencing and microarray of gene files have been widely applied among different species and diseases. Such transcriptome research can not only testify to gene functions and structures comprehensively, but also disclose the specific pathophysiologic mechanism underlying disease. However, a majority of the existing research remains limited to the coding RNA. On its rising slope, the quantity and quality of research on ncRNAs remains insufficient. Circular RNAs (circRNAs), the novel type of endogenous ncRNAs, has widely attracted scientists' attention lately (7). Compared with linear RNAs, the structure of covalently closed loops determines the stability of circRNAs by providing resistance to the exonuclease (8). With the development of this technique, an increasing amount of research has revealed the features of a majority of circRNAs in mammalian cells. They are endogenous, conserved, stable, and abundant, making them ideal therapeutic targets or potential biomarkers for the future (9–12).

A previous study applying microRNA microarrays has detected several miRNAs with significant expressive changes upon hepatic IPC following IRI (13). Moreover, a recent study conducted by Ye Z et al. investigated differentially expressed cirRNAs in liver IRI (14). Using a microarray, Zhang P et al. went further, comparing cirRNAs in hepatic IRI with or without IPostC and regarded mmu_circRNA_005186 as having a potential protective role in IPostC (15). However, systematic detection of the dysregulated circRNAs and mRNAs in response to IPC intervention still remains to be elucidated. Here, we suppose that the expressive alteration of circRNAs may closely relate to the pathophysiologic mechanism of hepatic IRI and may contribute to IPC-mediated protection against hepatic IRI. Comprehensive application of microarray, quantitative real-time PCR (qRT-PCR), and a progressive analysis were applied to reveal the expressive changes of circRNAs in the model of hepatic IRI with or without IPC and to investigate the protective circRNAs related to hepatic IPC intervention. Our data may offer a new understanding of the mechanical basis of IPC and may provide a possible research target for prospective studies.

Materials and Methods

Experimental Animals

C57BL/6 male mice were obtained for our study. Each mouse weighed from 18 to 22 g and was fostered in a temperature-constant room and nourished by a standard chow pellet diet. Preoperative fasting was performed 12 h before the surgical procedure without water prohibition. The experiments were performed in accordance with the Institutional Animal Committee guidelines. All animal protocols were approved by the Institutional Ethics Committee.

Model Construction and Experimental Groups

The hepatic IRI model was established by hepatic artery occlusion with clips as described previously (16). Briefly, after being anesthetized by pentobarbital (50 mg/kg body weight) and laparotomy, the mouse's portal vein and hepatic artery were isolated and temporarily clipped by atraumatic microvascular clip. After 1 h of ischemia to 70% of the liver, the vascular clip was wiped off to allow 6 h hepatic reperfusion. The ischemia extent was assessed by pale color of liver or pulselessness. Also, the recovery of the pale color and the pulse were considered to confirm valid hepatic reperfusion.

Three groups of the mouse model were constructed (n = 8 in each group): (1) sham group: sham model underwent portal vein and hepatic artery isolation with no occlusion; (2) IRI group: IRI model underwent IRI surgery without IPC; (3) IPC group: IPC was induced by 10 min of ischemia and 10 min of reperfusion before the intervention of IRI as reported previously (17). All the mice were sacrificed, and blood and organ samples were obtained at the end.

Histological Analysis

To achieve hepatic injury, we collected hepatic samples from the different groups of mice and embedded them in paraffin. Then, 5 μm thickness consecutive sections were stained separately with hematoxylin and eosin (H&E), immunohistochemistry of F4/80, and Ly6G + and evaluated by microscopy.

Analyses of Hepatic Injury and Cytokine Levels

Serum levels of alanine aminotransferase (ALT) and aspartate aminotransferase (AST) levels were tested using assay kits following the manufacturer's instructions (Nanjing Jiancheng). The levels of interleukin (IL)-6 and tumor necrosis factor (TNF)-α were measured using immunosorbent assay kits from the ENGTON Bio-engineering Limited Company (Shanghai, China) following the kit instructions.

RNA Extraction

Tissue RNA was extracted from each hepatic sample using Trizol (Invitrogen, Carlsbad, CA, USA) following the instructions. The concentration and purity of RNA were accessed by a NanoDrop 2000 spectrophotometer (Thermo Scientific, MA, USA). Genomic DNA (gDNA) residual and RNA integrity were measured by denaturing agarose gel electrophoresis.

Labeling and Hybridization

We selected three samples randomly from each group for microarray studies. RNA labeling and hybridization were done following the Agilent One-Color Microarray-Based Gene Expression Analysis construction. The microarray hybridization was applied following the manufacturer's instructions; five probes were used to increase the confidence of each transcript. To enrich circRNAs, total RNA of each sample was soaked with RNase R (Epicentre, Madison, WI, USA) to remove linear RNAs. Then, the Arraystar Labeling Kit (Arraystar, Rockville, USA) was used to transcribe circRNAs into fluorescent cRNA. Such cRNAs were then hybridized to the Arraystar circRNA array V₂.0 (8x15K). The Agilent Scanner was applied to detect the arrays following washing the slides. In the mRNA study, the whole mRNA microarray (Agilent Technology, CA, USA) was conventionally applied.

Microarray Information Analyses

Array data was processed by Agilent Feature Extraction software (version 11.0.1.1) and R software limma package. Low-intensity filtering was then applied after quantile normalization. The prominent different expressive circRNAs/mRNAs were sorted following the fold-change cutoff (fold change ≥ 1.5) or using volcano plot filtering. Distinctive circRNA/mRNA expression patterns were shown by hierarchical clustering.

Pathway Analyses

Gene Ontology (GO) analyses (http://www.geneontology.org) were conducted to set up gene annotation built upon different organisms. Gene functions include cellular components (CC), biological processes (BP), and molecular functions (MF). The -log10 (P-value) represented GO score, denoting the abundance in different genes. Also, pathway clusters were harvested using Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses relying on existing molecular networks from various gene data. The -log10 (P-value) indicated the KEGG score, denoting the significant relationship with the putative pathway.

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

The HiScript RT SuperMix for qPCR (+gDNA wiper) (Vazyme, Nanjing, China) was applied to reverse transcribe the total RNAs of each group following the protocol. QPCR using the SYBR Green kit (TaKaRa, Dalian, China) was applied to evaluate the expression levels of the circRNAs and mRNAs. Primers for mRNAs were routinely used. To amplify the circRNA junctions, we designed specific divergent primers based on the sequence obtained from the database “circBase” (http://www.circbase.org/). By using the 2-ΔΔCt method that normalizes against the expression of the β-actin gene, the circRNA and mRNA relative expressions were calculated. GenePharma (Shanghai, China) was applied to design, verify, and synthesize all primers. All sequences are as in Table 1.

Prediction of circRNA/miRNA Coactions

Arraystar's commercial miRNA coaction prediction software was applied to predict the circRNA/miRNA coaction. The database is based upon the existing online informatics tools miRanda (http://www.microrna.org/) and TargetScan (http://www.targetscan.org/). The Arraystar's miRNA scores of support vector regression (mirSVR) were used as a basis to rank miRNAs, and the miRNAs ranking the highest were considered for further analysis.

circRNA–miRNA–mRNA Pathway Prediction

Arraystar's software was applied to search circRNA_017753 MERs and select top potential target miRNAs based on match sequences. Then, top target miRNAs for the circRNA were selected according to the databases miRDB (http://mirdb.org/) and TargetScan. Based on the ceRNA theory, the direction of circRNA and mRNA changes in the same orientation. To reveal the beneficial mechanisms of the circRNA_017753 in hepatic IPC, we identified predicted downstream mRNA targets in the IRI group and elevated in the IPC group and constructed a Venn diagram. In the end, we selected the overlaying mRNAs and constructed the circRNA–miRNA–mRNA pathway.

Statistics

In the microarray data, a circRNA/mRNA fold change ≥ 1.5 was selected, and P < 0.05 was conventionally regarded as significant. GraphPad Prism 5.0 (GraphPad Software, CA, USA) was applied to analyze other data, and all data were indicated as the mean ± standard deviation. One-way analysis of variance followed by the Student–Newman–Keuls test were applied to compare data with normal distribution. The Kruskal–Wallis test and the Wilcoxon rank sum test with Bonferroni adjustments were used to compare data with nonnormal distributions. Means between two groups were compared by a two-tailed Student t-test. P < 0.05 was regarded as significant.

Results

Evaluation of Hepatic IRI and Hepatic Protection by IPC

The blood index of hepatic injury after hepatic IRI with or without IPC were tested to identify the beneficial effect of IPC in the study. As presented in our H&E staining, liver IRI induced a marked hepatocellular necrosis. In contrast, mice with IPC presented with only minor signs of ischemic congestion and necrosis. Our histochemical evaluation of F4/80 showed an intense inflammatory process occurring in ischemic parenchyma and alleviated levels in IPC. Also, the immunohistochemistry of Ly6G + showed a decreased level of infiltrating Ly6G + cells (a neutrophil plasma membrane biomarker) comparing the IPC and IRI groups (Figure 1A). The hepatic injury by IRI and the beneficial effect of IPC were further identified by serum examination. Serum ALT and AST were upregulated after hepatic IRI, and IPC remarkably decreased serum transaminase levels (Figure 1B). In addition, IPC ameliorated the TNF-α and IL-6 increase caused by IRI (Figure 1C). The results indicate that the current IPC intervention is effective enough to protect the liver against IRI. Hence, the model is suitable for comparing the circRNAs among IRI with or without IPC intervention.

Figure 1

Comparison of circRNA and mRNA in Hepatic IRI and IPC

To compare circRNA expression in the (1) sham, (2) IRI, and (3) IRI + IPC groups, we applied microarray analysis for circRNA and mRNA expression profiles of three samples in each. About 10,535 circRNA and 10,654 mRNA targets in the liver were spotted by microarray. Among them, 77 circRNAs and 686 mRNAs in the IRI group altered more than 1.5-fold compared with those in the sham group (P < 0.05). To be precise, 39 circRNAs and 432 mRNAs were increased, and 38 circRNAs and 254 mRNAs were decreased compared with the sham group (Figures 2A, 3A). Moreover, as a result of IPC intervention, a total of 50 circRNA and 95 mRNA alterations with significance (fold change ≥ 1.5, P < 0.05) were detected, in which 43 circRNAs and 64 mRNAs increased, 7 circRNAs and 31 mRNAs decreased in the IPC group when compared with IRI (Figures 2B, 3B). Based on fold change, we summarize the top 10 altered circRNAs in Table 2 and the top 10 altered mRNAs in Table 3. We then classified the significant altered circRNAs into various groups. Comparing the altered circRNAs in the IRI and sham groups, there were 77 exonic, 14 sense overlapping, 3 intronic, 5 antisense, and 1 intergenic and 74 exonic, 12 sense overlapping, 4 intronic, 6 antisense, and 4 intergenic when comparing the IPC and IRI groups (Figure 2C). Also, by using hierarchical clustering, we were able to cluster the differentially altered circRNAs and mRNAs of comparison groups and visualize the expression profiles. As indicated by the data, circRNAs and mRNAs in the IRI group were dramatically different from those of the sham group although IPC intervention significantly altered the IRI-induced expressions (Figures 2D, 3C).

Figure 2

Figure 3

circRNA and mRNA qRT-PCR Verification

qRT-PCR was further applied to confirm the microarray results. Seven circRNAs and five mRNAs were selected considering their fold change, P-value, and raw data. As qRT-PCR results show in Figure 4, the expression of circRNA_007095, Gclc, and Krt18 increased, and the circRNA_017753, circRNA_010415, circRNA_000895, circRNA_001946, and Hacl1 decreased enormously after hepatic IRI compared with the sham group. Otherwise, compared with IRI, IPC intervention remarkably elevated the expression of circRNA_027197, circRNA_010498, and Stim1, whereas Phc3 were reduced. The verified results are accordant with the data, indicating the dependability of our microarray profile.

Figure 4

Analysis of Variously Expressed mRNAs

Concerning that biologically similar genes share the same patterns of change, we decided to excavate the variations of coding gene expression among the IRI and IPC groups to reveal the behavior of these variously expressed circRNAs. GO enrichment analysis was used first for variously expressed mRNAs. In response to hepatic IRI, upregulated mRNAs were most related with response to organic substance and cellular response to chemical stimulus, and the downregulated mRNAs were most related to small molecule metabolic and organic acid metabolic processes (Figure 5A). Compared with the IRI, IPC mostly elevated the mRNAs related with response to hormone and regulation of signal transduction; meanwhile the downregulated mRNAs of IPC were most involved in regulation of cellular metabolic processes and organic cyclic compound metabolic processes (Figure 5B).

Figure 5

Applying the pathway analysis of KEGG, we found that fluid shear stress and atherosclerosis and protein processing in the endoplasmic reticulum were involved in the upregulated mRNAs, whereas the pathways of phenylalanine metabolism and peroxisome were most involved in the decreased mRNAs in the IRI group (Figure 6A). Notably, pathways related to the inflammatory process and apoptosis were the most involved in the IRI process as listed, such as pathways of PI3K-AKT, NF-Kappa B, IL-17 signaling, and apoptosis. Moreover, compared with the IRI group, IPC can significantly upregulate mRNAs involved in the pathways of ovarian steroidogenesis and lysine degradation and downregulate the mRNAs related with drug metabolism-cytochrome p450 and chemical carcinogenesis pathways (Figure 6B). Our computational analysis results may contribute dramatically to the pathogenesis of hepatic IRI and IPC intervention.

Figure 6

Identification of circRNAs Related to Hepatic Protection by IPC

Concerning the protection effect of IPC intervention against IRI, we compared the circRNA alteration and direction of alterations of the three groups to reveal the potential relevance between circRNA changes and IPC protection. As a result, we sorted circRNAs that showed up with opposite alteration directions between the comparison groups (sham vs. IRI and IRI vs. IPC). Following this method, we selected three circRNAs that were upregulated in IRI but downregulated in the IPC group and 12 circRNAs that showed a completely opposite direction (Figure 7A). Concerning the type of circRNAs and data sources of circRNAs, we then selected only one circRNA, circRNA_017753, that was considered significatively related to the IPC protective effect. The gnomic locus of circRNA_017753 is on chromosome 17, and it is spliced from Mapk14. The circRNA_017753 expression level was then confirmed by qRT-PCR; IPC intervention can significantly restore the downregulation caused by hepatic IRI. Our data indicate that circRNA_017753 may play a protective role in hepatic IRI that arouses our desire to learn the potential role of circRNA_017753. In addition, we also compared our data with the microarray data of hepatic IPostC intervention in IRI (GSE117524) (15). Though circRNA_010498 does not show a comparative fold change to circRNA_017753 in our data, the intersection of research (fold change ≥ 1.5, P-value <0.05) indicates its potential significant protective effect in both IPC and IPostC intervention and is worth further study (Supplementary Figure 1).

Figure 7

Prediction of miRNA and CircRNA–miRNA–mRNA Pathway for CircRNA_017753

It is well-known that circRNA regulates miRNAs by interacting with miRNA response elements (MREs). Such way of interaction can competitively suppress miRNAs' activity and is, thus, called miRNA sponges. To figure out the function of circRNA_017753, we applied Arraystar's prediction software, which combined TargetScan and miRanda databases. We predict and list the five highest ranking target miRNAs of circRNA_017753 and annotation of their circRNA/miRNA interactions (Table 4, Figure 8A).

Figure 8

To further investigate the potential mechanisms of circRNA_017753 in IPC, we construct the possible circRNA_017753–miRNA–mRNA pathways. Applying Arraystar's miRNA prediction software, 119 genes were selected related to the 5 miRNAs mentioned (Supplementary Table 1). Because of the widely accepted ceRNA concept that circRNAs may positively regulate mRNAs, we paid more attention to the mRNAs downregulated in IRI but significantly inhibited by IPC. By establishing the intersections of 65 mRNAs significantly increased by IPC and 119 predicted mRNA targets connected to circRNA_017753, we identified only one overlapping mRNA, Jade1, and predicted three following circRNA–miRNA–mRNA signals: circRNA_017753–miR-218-5p–Jade1, circRNA_017753–miR-7002-3p–Jade1, and circRNA_017753–miR-7008-3p–Jade1. These regulatory ceRNA signaling pathways may play important roles in the mechanisms of IPC protection and deserve further study.

Discussion

Hepatic IRI contributes significantly to organ damage in the surgical procedures of hepatic resections and liver transplantation. Such injury possesses a high mortality rate caused by an intense inflammatory process occurring in the ischemic liver. As presented in our research, not only did serum TNF-α and IL-6 increase in the IRI group, the evaluation of F4/80 and infiltrating of Ly6G + staining demonstrated an intense inflammatory process in IRI. IPC has been most investigated in the past decades as a life-saving intervention. Though widely applied clinically, the exact molecular mechanisms behind the protection effect of IPC still remain largely unclear. A majority of the existing research mainly focuses on protein-coding RNA while research on ncRNAs (long ncRNAs and circRNAs) remain insufficient. It was not until the mechanism of powerful miRNA sponges revealed circRNAs as a potential target for treatment and diagnosis scientists' attention was widely attracted. In our research, by utilizing circRNA microarray analysis, we first report circRNA alteration profiles induced by hepatic IRI and the protection intervention of IPC systematically. More importantly, through bioinformatic comparison of circRNA alteration and the direction of alteration, we identified one possible circRNA, circRNA_017753, related with hepatic protection by IPC intervention. By comparing the data with the previous microarray data of hepatic IPostC, we identified one possible circRNA, circRNA_010498, that may have a potential protective effect in both IPC and IPostC. Also, by integrated application of TargetScan and miRanda databases, we predicted three circRNA–miRNA–mRNA pathways that may take effect in the mechanisms of IPC protection against hepatic IRI.

In our research, the microarray profiles identified 77 circRNAs that were significantly altered (39 up- and 38 downregulated) after hepatic IRI, whereas a total of 50 circRNAs altered significantly due to the IPC intervention (43 up- and 7 downregulated). As previously reported, the majority of circRNAs arise from exons; our profiles identified that about 2/3 altered circRNAs were exonic (18). Though none of the circRNAs in our top dysregulated list have been reported before as the reason for nascent circRNA functional study, it is encouraging to see some mRNAs regarded as alternative transcripts of these circRNAs reported to take effect in different physiopathologic mechanisms related to the IRI process. For example, Bach1 is the alternative transcript of circRNA_29990 and circRNA_29992 (2 of the top 5 upregulated circRNAs in hepatic IRI), and its protein inhibits the transcription of HO-1 and related genes involved in the oxidative stress response by binding to Maf genome recognition elements (19). Bach1 deficiency may increase resistance to ischemic stresses by elevating HO-1 expression (20, 21). Several studies also reveal its important role in ischemic or oxidative damage (22–24). More importantly, Bach1 is identified to repress Wnt/β-Catenin signaling and angiogenesis in peripheral ischemic injury in a recent study (25). Magi2 is the transcript of circRNA_38159 (one of the top five downregulated circRNAs listed in IPC); it encodes scaffolding proteins binding to PTEN and is identified an important element in the ischemic injury of the central nervous system. As it is well-known that circRNAs may regulate its linear counterparts (26), the circRNAs listed may play crucial roles in hepatic IRI and IPC by regulating the transcription of the parent genes.

The data on dysregulated mRNAs are also inspiring. In the list of top 10 dysregulated mRNAs in hepatic IRI, G0s2 (G0/G1 switch gene 2) decreased the most among all downregulated mRNAs and is widely recognized as a direct activator of oxidative phosphorylation at the early phase of hypoxia (27, 28). A recent study just identified that its overexpression can alleviate ATP decrease in myocardial cells and increase their hypoxic resistance during ischemia (29), suggesting its crucial role in hepatic IRI. Also, we can identify some new clues on the protection mechanisms of IPC in our profile of dysregulated mRNAs between IPC and I/R groups. For instance, Ddit4 (DNA-damage inducible transcript 4), widely recognized as an autophagy regulator by negatively regulating mTORC1, is found downregulated in the IPC group (30). Our result is consistent with the previous research that Ddit4 is a novel protection molecule that prevents ischemic injury in hepatocytes (31). Also, the protective effect of Ddit4 is found to be vital in cardiac and cerebral IRI (32–34). All these consistencies prove the reliability of our profile and provide us the credibility of other listed mRNAs in the mechanisms of IPC protection against hepatic IRI, thus making it worth further study.

We can also profit from GO and KEGG pathway analysis of altered mRNAs, which may reveal the key processes of hepatic IRI pathogenesis. The data emphasize the crucial roles of the inflammatory process and apoptosis in the IRI process, such as pathways of PI3K-AKT, NF-Kappa B, IL-17 signaling, and apoptosis. Moreover, the pathways of the metabolic process and regulation of signaling are prominent when comparing the IPC and IRI groups, indicating that related genes may play vital roles in IPC protection mechanisms. Our observation was partially consistent with a previous study revealing the relationship with metabolic process and the IPC protective mechanism (35).

It is worth mentioning that, in our data, only one circRNA, circRNA_017753, was selected to have a relationship with the protective mechanisms of IPC with the most possibility. Our data on the microarray and qRT-PCR both confirm its decrease in the IRI group and upregulation with IPC intervention. We speculate that the decreased level of circRNA_017753 may reflect hepatic dysfunction during IRI, and IPC may alleviate the hepatic dysfunction by the circRNA adjustment. Thus, circRNA_017753 may play a crucial role in the protective mechanism of hepatic IRI and may be a possible therapeutic target of hepatic injury.

To make the analysis more comprehensive, we also compared our data with previous data on hepatic IRI and IPostC prevention (14, 15). Because of different experimental environments and model construction methods (such as reperfusion 4 h after ischemia 1 h by Zhang P et al.), though reported circRNAs, such as circRNA_005186, exhibit the same alteration trend in our study in both IRI and prevention, the fold change of these circRNAs does not meet our screening criteria (fold change ≥ 1.5, P-value <0.05). However, it is worth noting that, by constructing the intersection of our data and GSE117524, circRNA_010498 shows sufficient significance as a protective factor in both IPC and IPostC intervention.

It is widely recognized today that circRNAs regulate gene expression by serving as miRNA sponges. Based on sequence comparing and bioinformatic methods, we predicted the five potential candidate miRNAs for circRNA_017753 with the most possibility. Then, by applying prediction software that combines TargetScan and miRanda databases, we are able to construct three circRNA–miRNA–mRNA regulatory axes that may have a protective effect in IPC intervention. Even though three circRNA–miRNA–mRNA axes have not been previously reported, the predicted molecule of Jade1 was found to have a huge possibility to play a part in the hepatic IPC protection mechanism. As a component of the HBO1 complex, Jade1 was identified as a key regulator of apoptosis (36). In previous study, by promoting acetylation of histones, it functioned as a key regulator of cycle progression and redifferentiation in renal tubule regeneration after IRI (37). Further studies on Jade1 and related circRNA–miRNA–mRNA axes are in progress in our laboratory.

The current study has some limitations that should be admitted. First, our data was obtained in an animal model and may not fully represent human pathological pathways. However, there is a high degree of similarity between human and mouse research in circRNAs. As it is demonstrated in previous research that most circRNAs present to be conserved between mouse and human (38), most of the circRNAs detected in our microarray data are conserved and of great interest for human research. Second, our research focuses on integral hepatic parenchyma injury and IPC prevention, thus less consideration was put into the design of circRNA profiling in different hepatic zones. In the development of spatial transcriptomics and microarray methods, specific profiling for circRNA expression in different hepatic zones would be of great interest and worth further study.

Conclusion

Our study for the first time delineates the expression data of dysregulated circRNAs and mRNAs in response to hepatic IRI and IPC intervention. Our profile and bioinformatic analyses provide numerous novel clues on the pathophysiologic mechanism of IPC protection. Several potential circRNA–miRNA–mRNA axes we predicted may offer promising targets for hepatic ischemic prevention and treatment.

References

• 1.

  SakaiNVan SweringenHLSchusterRBlanchardJBurnsJMTevarADet al. Receptor activator of nuclear factor-κB ligand (RANKL) protects against hepatic ischemia/reperfusion injury in mice. Hepatology. (2012) 55:888–97. 10.1002/hep.24756

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2.

  WalshKBToledoAHRivera-ChavezFALopez-NeblinaFToledo-PereyraLH. Inflammatory mediators of liver ischemia-reperfusion injury. Exp Clin Transplant. (2009) 7:78–93. Available online at: https://apps.webofknowledge.com/full_record.do?product=UA&search_mode=GeneralSearch&qid=1&SID=5FlwRUV3HJhawCTnGz3&page=1&doc=4

  • Pubmed Abstract
  • Google Scholar
  • View reference in article

• 3.

  PerryBCSoltysDToledoAHToledo-PereyraLH. Tumor necrosis factor-α in liver ischemia/reperfusion injury. J Invest Surg. (2011) 24:178–88. 10.3109/08941939.2011.568594

  • CrossRef
  • Google Scholar
  • View reference in article

• 4.

  MurryCEJenningsRBReimerKA. Preconditioning with ischemia: a delay of lethal cell injury in ischemic myocardium. Circulation. (1986) 74:1124–36. 10.1161/01.cir.74.5.1124

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 5.

  ShawkyHMYounanSMRashedLAShoukryH. Effect of recombinant erythropoietin on ischemia-reperfusion-induced apoptosis in rat liver. J Physiol Biochem. (2012) 68:19–28. 10.1007/s13105-011-0114-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 6.

  KnudsenARKannerupASGrønbækHDutoitSHNyengaardJRFunch-JensenPet al. Quantitative histological assessment of hepatic ischemia-reperfusion injuries following ischemic pre- and post-conditioning in the rat liver. J Surg Res. (2013) 180:e11–20. 10.1016/j.jss.2012.03.036

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 7.

  QuSYangXLiXWangJGaoYShangRet al. Circular RNA: a new star of noncoding RNAs. Cancer Lett. (2015) 365:141–8. 10.1016/j.canlet.2015.06.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8.

  KristensenLSAndersenMSStagstedLVWEbbesenKKHansenTBKjemsJ. The biogenesis, biology and characterization of circular RNAs. Nat Rev Genet. (2019) 20:675–91. 10.1038/s41576-019-0158-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9.

  SalzmanJGawadCWangPLLacayoNBrownPO. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PLoS ONE. (2012) 7:e30733. 10.1371/journal.pone.0030733

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10.

  JeckWRSorrentinoJAWangKSlevinMKBurdCELiuJet al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA. (2013) 19:141–57. 10.1261/rna.035667.112

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11.

  MemczakSJensMElefsiniotiATortiFKruegerJRybakAet al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. (2013) 495:333–8. 10.1038/nature11928

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12.

  SalzmanJChenREOlsenMNWangPLBrownPO. Cell-type specific features of circular RNA expression. PLoS Genet. (2013) 9:e1003777. 10.1371/journal.pgen.1003777

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13.

  XuCFYuCHLiYM. Regulation of hepatic microRNA expression in response to ischemic preconditioning following ischemia/reperfusion injury in mice. Omics. (2009) 13:513–20. 10.1089/omi.2009.0035

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14.

  YeZKongQHanJDengJWuMDengH. Circular RNAs are differentially expressed in liver ischemia/reperfusion injury model. J Cell Biochem. (2018) 119:7397–405. 10.1002/jcb.27047

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15.

  ZhangPMingYYeQNiuY. Comprehensive circRNA expression profile during ischemic postconditioning attenuating hepatic ischemia/reperfusion injury. Sci Rep. (2019) 9:264. 10.1038/s41598-018-36443-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16.

  KurodaSTashiroHIgarashiYTanimotoYNambuJOshitaAet al. Rho inhibitor prevents ischemia-reperfusion injury in rat steatotic liver. J Hepatol. (2012) 56:146–52. 10.1016/j.jhep.2011.04.029

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17.

  NakayamaHYamamotoYKumeMYamagamiKYamamotoHKimotoSet al. Pharmacologic stimulation of adenosine A2 receptor supplants ischemic preconditioning in providing ischemic tolerance in rat livers. Surgery. (1999) 126:945–54. 10.1016/s0039-6060(99)70037-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18.

  Ashwal-FlussRMeyerMPamudurtiNRIvanovABartokOHananMet al. circRNA biogenesis competes with pre-mRNA splicing. Mol Cell. (2014) 56:55–66. 10.1016/j.molcel.2014.08.019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 19.

  OgawaKSunJTaketaniSNakajimaONishitaniCSassaSet al. Heme mediates derepression of Maf recognition element through direct binding to transcription repressor Bach1. Embo J. (2001) 20:2835–43. 10.1093/emboj/20.11.2835

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20.

  OmuraSSuzukiHToyofukuMOzonoRKohnoNIgarashiK. Effects of genetic ablation of bach1 upon smooth muscle cell proliferation and atherosclerosis after cuff injury. Genes Cells. (2005) 10:277–85. 10.1111/j.1365-2443.2005.00832.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21.

  YanoYOzonoROishiYKambeMYoshizumiMIshidaTet al. Genetic ablation of the transcription repressor Bach1 leads to myocardial protection against ischemia/reperfusion in mice. Genes Cells. (2006) 11:791–803. 10.1111/j.1365-2443.2006.00979.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22.

  TanimotoTHattoriNSenooTFuronakaMIshikawaNFujitakaKet al. Genetic ablation of the Bach1 gene reduces hyperoxic lung injury in mice: role of IL-6. Free Radic Biol Med. (2009) 46:1119–26. 10.1016/j.freeradbiomed.2009.01.017

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23.

  WarnatzHJSchmidtDMankeTPicciniISultanMBorodinaTet al. The BTB and CNC homology 1 (BACH1) target genes are involved in the oxidative stress response and in control of the cell cycle. J Biol Chem. (2011) 286:23521–32. 10.1074/jbc.M111.220178

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24.

  KondoKIshigakiYGaoJYamadaTImaiJSawadaSet al. Bach1 deficiency protects pancreatic β-cells from oxidative stress injury. Am J Physiol Endocrinol Metab. (2013) 305:E641–648. 10.1152/ajpendo.00120.2013

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25.

  JiangLYinMWeiXLiuJWangXNiuCet al. Bach1 represses Wnt/β-catenin signaling and angiogenesis. Circ Res. (2015) 117:364–75. 10.1161/circresaha.115.306829

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26.

  LiZHuangCBaoCChenLLinMWangXet al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. (2015) 22:256–64. 10.1038/nsmb.2959

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27.

  KiokaHKatoHFujikawaMTsukamotoOSuzukiTImamuraHet al. Evaluation of intramitochondrial ATP levels identifies G0/G1 switch gene 2 as a positive regulator of oxidative phosphorylation. Proc Natl Acad Sci USA. (2014) 111:273–8. 10.1073/pnas.1318547111

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28.

  HayashiTAsanoYShintaniYAoyamaHKiokaHTsukamotoOet al. Higd1a is a positive regulator of cytochrome c oxidase. Proc Natl Acad Sci USA. (2015) 112:1553–8. 10.1073/pnas.1419767112

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29.

  KamikuboKKatoHKiokaHYamazakiSTsukamotoONishidaYet al. A molecular triage process mediated by RING finger protein 126 and BCL2-associated athanogene 6 regulates degradation of G(0)/G(1) switch gene 2. J Biol Chem. (2019) 294:14562–73. 10.1074/jbc.RA119.008544

  • CrossRef
  • Google Scholar
  • View reference in article

• 30.

  DennisMDMcGheeNKJeffersonLSKimballSR. Regulated in DNA damage and development 1 (REDD1) promotes cell survival during serum deprivation by sustaining repression of signaling through the mechanistic target of rapamycin in complex 1 (mTORC1). Cell Signal. (2013) 25:2709–16. 10.1016/j.cellsig.2013.08.038

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 31.

  ChoSSKimKMYangJHKimJYParkSJKimSJet al. Induction of REDD1 via AP-1 prevents oxidative stress-mediated injury in hepatocytes. Free Radic Biol Med. (2018) 124:221–31. 10.1016/j.freeradbiomed.2018.06.014

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32.

  ChenQZhouYRichardsAMWangP. Up-regulation of miRNA-221 inhibits hypoxia/reoxygenation-induced autophagy through the DDIT4/mTORC1 and Tp53inp1/p62 pathways. Biochem Biophys Res Commun. (2016) 474:168–74. 10.1016/j.bbrc.2016.04.090

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33.

  ParkKMTeohJPWangYBroskovaZBayoumiASTangYet al. Carvedilol-responsive microRNAs, miR-199a-3p and−214 protect cardiomyocytes from simulated ischemia-reperfusion injury. Am J Physiol Heart Circ Physiol. (2016) 311:H371–383. 10.1152/ajpheart.00807.2015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34.

  ParkJALeeCH. Time-course change of redd1 expressions in the hippocampal CA1 region following chronic cerebral hypoperfusion. Cell Mol Neurobiol. (2017) 37:563–9. 10.1007/s10571-016-0385-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35.

  ZhangXJChengXYanZZFangJWangXWangWet al. An ALOX12-12-HETE-GPR31 signaling axis is a key mediator of hepatic ischemia-reperfusion injury. Nat Med. (2018) 24:73–83. 10.1038/nm.4451

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36.

  ZhouMIFoyRLChitaliaVCZhaoJPanchenkoMVWangHet al. Jade-1, a candidate renal tumor suppressor that promotes apoptosis. Proc Natl Acad Sci USA. (2005) 102:11035–40. 10.1073/pnas.0500757102

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 37.

  HavasiAHaegeleJAGallJMBlackmonSIchimuraTBonegioRGet al. Histone acetyl transferase (HAT) HBO1 and JADE1 in epithelial cell regeneration. Am J Pathol. (2013) 182:152–62. 10.1016/j.ajpath.2012.09.017

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 38.

  Rybak-WolfAStottmeisterCGlaŽarPJensMPinoNGiustiSet al. Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed. Mol Cell. (2015) 58:870–85. 10.1016/j.molcel.2015.03.027

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article