Abstract

Mitochondrial ATP-sensitive potassium channel (mito) is a common end effector of many protective stimuli in myocardial ischemia-reperfusion injury (MIRI). However, the specific molecular mechanism underlying its myocardial protective effect is not well elucidated. We characterized an anoxia/reoxygenation (A/R) model using freshly isolated adult rat cardiomyocytes. Mito status was interfered with its specific opener diazoxide (DZ) or blocker 5-hydroxydecanote (5-HD). Digital gene expression (DGE) and bioinformatic analysis were deployed. Three energy metabolism related genes (MT-ND6, Idh2, and Acadl) were upregulated when mito opened. In addition, as many as 20 differentially expressed genes (DEGs) were significantly enriched in five energy homeostasis correlated pathways (PPAR, TCA cycle, fatty acid metabolism, and peroxisome). These findings indicated that mito opening in MIRI resulted in energy mobilization, which was confirmed by measuring ATP content in cardiomyocytes. These causal outcomes could be a molecular mechanism of myocardial protection of mito and suggested that the mito opening plays a physiologic role in triggering cardiomyocytes’ energy homeostasis during MIRI. Strategies of modulating energy expenditure during myocardial ischemia-reperfusion may be promising approaches to reduce MIRI.

1. Introduction

Myocardial infarction has been a leading cause of death worldwide. The prognosis of acute myocardial infarction has been dramatically improved due to the advances of both catheterization techniques and reperfusion therapy by coronary mechanical and pharmacological intervention methods. However, strategies to limit myocardial ischemia-reperfusion injury (MIRI), thus reducing infarct size, have not been well applied in clinical settings.

Although myocardial ischemia-reperfusion (IR) induces lethal injury in the heart, after some artificial interventions, the cardiomyocytes and the heart tissue therein have powerful endogenous mechanisms to protect themselves from oxidative stress, energy deficiency, protein aggregation, and organelle malfunction, thereby minimizing MIRI [1]. For example, Murry et al. in 1990 [2, 3] first proposed that ischemic preconditioning (IPC) may protect the heart by reducing myocardial energy demand during myocardial ischemia and decreasing cell death by preserving ATP content and/or reducing catabolite accumulation. Following the discovery of the mitoKATP channel locating at the inner mitochondrial membrane in 1991 [4], Garlid et al. and Liu et al. [5, 6] demonstrated it as a trigger of IPC. Pharmacological intervention mimicking the IPC has currently been considered as a promising modality for the treatment of MIRI. Similar myocardial protection can be produced by drugs such as diazoxide (DZ) that open mitoKATP [5, 7]. Conversely, mitoKATP blockers (5-hydroxydecanote (5-HD) or glibenclamide) cancelled the effect of preconditioning and pharmacological cardioprotection [5, 6, 8]. It is also demonstrated that the pharmacological inhibition of the mitoKATP in early reperfusion abolished the infarct-limiting effects of IPost [911].

We have reported that mitoKATP opening was cardioprotective in MIRI [1214], but our understanding of its specific mechanism remained quite preliminary. To date, the main proposed mechanisms of cardioprotection by mitoKATP were various: swelling of mitochondria increased fatty acid oxidation (FAO), mitochondrial respiration, and ATP production [15]; inhibition of ATP hydrolysis during ischemia [16, 17] preserved ATP and decreased Ca2+ uptake in the cardiomyocytes. However, other endogenous mechanisms of cardioprotection of mitoKATP activation during IR remain to be elucidated.

Most of the in vitro studies used neonatal cardiac cells or immortal cardiac cell lines such as H9c2, which is physiologically different from adult cardiomyocytes [18]. For example, it is reported that neonatal cardiomyocytes were more resistant to hypoxia in comparison to adult ones [19, 20]. So, it may limit the extrapolation of the research results. We developed an A/R model using adult cardiomyocytes freshly isolated from rat to mimic the IR microenvironment in vivo; after all, MIRI is present almost exclusively in the adult population.

Compared with microarray and PCR-based technologies, digital gene expression (DGE) platform can provide adequate sequence coverage and quantitative accuracy to capture subtle changes resulting from mitoKATP opening. In this study, a molecular and bioinformatic pipeline permitted comprehensive analysis of the myocardial mRNA expression. Next-generation sequencing technology was employed and the impact of mitoKATP on the myocardial transcriptome signature of MIRI was explored to crystallize cardioprotective effects of mitoKATP.

2. Materials and Methods

2.1. Experimental Animals

Male Sprague-Dawley rats (250–300 g, 16–20 weeks) were provided by the Third Military Medical University (Chongqing, China) and maintained in specific pathogen free (SPF) animal facility in Zunyi Medical College under standardized conditions with 12 h light/dark cycles and free access to rat chow and water. All experimental procedures were performed according to the “Guide for the Care and Use of Laboratory Animals” in China (no. 14924, 2001) and approved by the Experimental Animal Care and Use Committee of Zunyi Medical College.

2.2. Isolation of Adult Cardiomyocytes

Rats were anesthetized with sodium pentobarbital (60 mg/kg, combined with 250 U/kg heparin, peritoneal injection). When rats had been successfully anesthetized, the chest cavity was opened and the heart excised rapidly. Ventricular cardiomyocytes were obtained by enzymatic digestion as previously described [24], with some necessary modification. Briefly, hearts were retrogradely perfused with 0.1% type 2 collagenase (Sigma, USA) at constant pressure (9 mL/min/g) on a Langendorff apparatus (Alcott Biotech, China); then the ventricle was scissored out and digested by type 2 collagenase solution in a beaker with manually shaking. Cells were filtered through a piece of gauze and washed 5 times to get rid of collagenase. Cardiomyocytes from one heart were evenly titrated into four 60 mm laminin-precovered Petri dishes. Three mL serum free modified M199 medium (Hyclone, USA, with 2 mM carnitine, 2 mM glutamine, 5 mM taurine, 5 mM creatine, and 0.8 mM EGTA) was added. After 3 hours’ incubation, the medium was replaced to get rid of noncardiomyocytes. Cell quality was confirmed with trypan blue exclusion test.

2.3. Anoxia/Reoxygenation and DZ Postconditioning in Adult Rat Cardiomyocytes

For each test, the 4 Petri dishes were placed in normoxia incubator for 20 hours before being randomly distributed to 4 groups: Control (Con), anoxia/reoxygenation (A/R), diazoxide (DZ), and DZ + blocker 5-hydroxydecanote (5-HD) (DZ5HD). Cardiomyocytes of Con were continuously cultured in normoxia environment for 105 min; A/R: under anoxia for 45 min and then reoxygenated for 60 min; DZ: anoxia for 45 minutes, reoxygenated with 50 μM DZ (Sigma, USA) for 5 min, and then reoxygenated without DZ for another 55 min; DZ5HD: anoxia for 40 min, anoxia with 100 μM 5-HD (Sigma, USA) for 5 min and then reoxygenated with DZ for 5 min and reoxygenated without DZ for another 55 min (Figure 1). Oxygen deprivation and reoxygenation were achieved by series of changes of the medium and incubators. Normoxia was set in a normoxia incubator (O2/CO2 incubator containing a humidified atmosphere of 5% CO2 and 95% air at 37°C). Anoxia was achieved in an anoxic incubator (O2/CO2 incubator containing a humidified atmosphere of of 5% CO2, 1% O2, and 94% N2 at 37°C) (Figure 1).

2.4. Intracellular Free Calcium ([Ca2+]i) Test

At the end of reoxygenation, was detected as previously reported [25]. Briefly, M199 was removed; cardiomyocytes of 4 groups were washed twice with PBS and loaded with Fluo-3 AM (Biotium, USA) at a final concentration of 10 μM and incubated for 30 min at 37°C in O2/CO2 incubator. The solution containing the Ca2+ probe was removed and cells were washed twice again with PBS. The average fluorescence intensity of concentration in labeled cells was detected under a laser scanning confocal microscope (TCS SP2 AOBS, Leica, Germany). The wavelength of excitation was set at 488 nm and the emission wavelength was 525 nm for Fluo-3 fluorescence reading. More than 20 cells from each group were randomly chosen for data analysis; their outlines were circled out and the fluorescence density of Fluo-3 was calculated with Leica confocal software (Leica, Germany).

2.5. Cell Viability Detection

Adult cardiomyocytes’ viability was detected with Cell Counting Kit-8 (CCK-8, Beyotime, China) in accordance with the manufacturer’s instructions. The same amount of cells was seeded into 24-well plates. At the end point of reoxygenation, 30 μL WST-8 solution was added into M199 to form a 3% WST-8 final concentration. Cells were incubated for 1 h before the mixture’s OD value was detected at 450 nm wavelength. The replicate size was 6 for each group.

2.6. RNA Extraction

At the end of reoxygenation (see Section 2.3), cardiomyocyte samples (3 replicates for 4 groups) were homogenized in TRIzol reagent (Invitrogen, USA) and vortexed with chloroform. The mixture was prepared at room temperature for 2 min and then centrifuged at 12000 ×g at 4°C. The aqueous phase was mixed with 100% ethanol and then filtered with a Qiagen RNeasy column. Subsequent steps for extraction of total RNAs were carried out as the Qiagen RNeasy kit (Qiagen, Germany) instructions described.

2.7. Tag Library Construction

The tag-seq library was constructed in accordance with the manufacturer’s workflow as previously described [26]. Briefly, 6 mg extracted total RNA was used for mRNA capture with magnetic Oligo (dT) beads. Then cDNA was synthesized and the bead-bound cDNA was subsequently digested with NlaIII. Fragments attached to Oligo (dT) beads were washed away. GEX NlaIII adapter was ligated to the free 5′ end of the digested bead-bound cDNA fragments. Individual cDNA libraries were PCR amplified and purified on a 6% acrylamide gel. Attached DNA fragments were used to create a sequencing flow cell with millions of clusters, which contained about 1000 copies of the templates. Templates were sequenced by the Illumina HiSeq 2500 equipment using the four-color DNA sequencing-by-synthesis (SBS) technology. Each lane generated millions of raw reads.

2.8. Data Processing and Statistical Analysis

To obtain high quality and reliable data, raw reads were filtered to remove potentially erroneous reads. Briefly, the 3′ adaptor sequences were trimmed, low-quality tags containing N were abandoned, and small tags and only 1 copy tag were removed before obtaining the clean reads. After filtering, all reads were annotated to Rat Genome V3.4 Assembly (http://rgd.mcw.edu/). All the clean reads were mapped to the reference database; the unambiguous tags were annotated. Copy number of the clean tags of each gene was normalized with the RPKM (reads per kilobase of exon per million mapped reads) method [27] to get the final gene expression.

2.9. Identification of DEGs

According to the method by Tarazona et al. [28], the NOISeq-real algorithm was employed to determine the value (corresponding to the value in differential gene expression detection) and screen genes [29, 30] differentially expressed between Con and A/R, A/R and DZ, and DZ and DZ5HD. In the present study, we considered a gene differentially expressed if the value was more than 0.8.

2.10. Gene Annotation with Gene Ontology and KEGG Pathway

GO (http://www.geneontology.org) provides a dynamic, controlled vocabulary. It comprises 3 independent ontologies: Biological Process, Molecular Function, and Cellular Component, each of which contains hundreds of terms. These terms reflect our understanding of the gene function.

KEGG Pathway database is for systematical analysis of gene functions, linking genomic information with higher order functional information. It provides an indication of the main biochemical and signal transduction pathways that DEGs are involved in.

Finally, the DEGs were enriched with GO (into ontologies and terms) and KEGG Pathway database.

2.11. RT-qPCR Analysis

Twenty-five DEGs were randomly selected for real-time quantitative PCR (RT-qPCR). The total RNA used for sequencing was reused to validate DGE sequencing. 500 ng RNA was reverse-transcribed into cDNA using a cDNA synthesis kit (Takara, Japan) in a final volume of 10 μL according to the manufacturer’s protocol. RT-qPCR was performed with the CFX Connect Real-Time system (Bio-Rad, USA) using a SYBR green PrimeScript RT kit (Perfect Real Time, Takara, Japan) based on the manufacturer’s instructions. The PCR conditions included predenaturing at 95°C for 30 s followed by 40 cycles of denaturation at 95°C for 10 s and combined annealing/extension at 58°C for 30 s. All the mRNA expression levels were calculated based on the comparative quantification method (). The β-actin gene was used as an internal control. The primer sequences were listed in Table 4.

2.12. ATP Quantitation in Cardiomyocytes

At the end of reoxygenation, the cardiomyocytes were scraped off and centrifuged at 1000 ×g for 5 min; the supernatant (M199 medium) was abandoned. 1 mL precooled 0.4 M HClO4 was added into the pellet and followed by ultrasonication and centrifugation at 10000 ×g for 20 min. The supernatant was collected and its pH was adjusted to 6.0 to 7.0 with 0.7 mL 1 M K2HPO4 before centrifugation again at 10000 ×g for another 20 min. All the above-mentioned procedures were conducted at 4°C. The supernatant was filtered through 0.22 μm membrane before high performance liquid chromatography (HPLC) analysis. The chromatographic conditions were as follows: work station: LC 20A (Shimadzu, Japan); column: WondaSil C18-WR (150 mm × 4.6 mm, μm; GL Sciences, Japan); column temperature: 25°C; mobile phase: buffered phosphate at pH 7.0; flow rate: 1 mL/min; detection wavelength: 254 nm; sample size: 10 μL. The ATP peaks of samples were determined in reference to the ATP standards (Sigma, USA). The amount of ATP was determined based on the standard curve and regression equation from ATP standard’s concentration and peak area. Protein content was measured by using the same sample. ATP level of each sample is normalized to protein content.

2.13. Statistical Analysis

The quantitative data were expressed as mean ± SD. For experiments of cardiomyocytes of the four groups, one-way analysis of variance (ANOVA) was performed; LSD or Dunnett’s T3 method was used to make multiple comparisons. A value of less than 0.05 was considered to be statistically significant. All data analyses were carried out using SPSS v.19.0 (IBM, USA).

3. Results

3.1. Isolated Adult Rat Cardiomyocytes

A high percentage (70–80%) of rod-shaped adult cardiomyocytes with clear striations and sharp outlines without visible vesicles were obtained with our method (Figure 2(a)).

3.2. [Ca2+]i and Cell Viability Detection

We used Fluo-3 AM to examine Ca2+ mobilizations in cardiomyocytes. In Con group, the level of was the lowest. Compared with Con, increased significantly in A/R (). After the applying of DZ, fluorescence decreased dramatically () compared with A/R while there was an apparent increase () in DZ5HD compared with DZ (Figures 2(b) and 2(c)). It indicated that DZ strongly inhibited the levels in adult rat cardiomyocytes.

Cardiomyocytes in A/R group possessed lower level of cell viability () when compared with Con. DZ group contained higher level of cell viability () when compared with A/R group, while DZ5HD group showed lower level of cell viability () when compared with DZ group (Figure 2(d)).

3.3. Quality Evaluation of DGE Reads

A summary of the DGE reads and their mapping to the rat genome database is presented in Supplementary Table 1  available online at  http://dx.doi.org/10.1155/2014/756576. For each group, more than 4.4 million clean reads were sequenced. Low-quality reads accounted for no more than 1.6% and modified Q30 bases rate no less than 97% in all the 12 libraries (Supplementary Table 1). Besides, perfect matched reads accounted for 60% and unique matched reads occupied more than 70% of all reads mapping to rat genome (Table 1), which revealed that the sample preparation and the sequencings were in perfect condition.

3.4. Sequencing Saturation Analysis

Samples with replicates of sequencing, sequencing saturation analysis can be performed to test whether the detected genes’ percent increased with total reads number. As shown in Supplementary Figure 1, for 3 replicates of 4 groups, when the total tag number came to 3 million, the genes number started to level out. When the total tag number reached 4 million, gene number inclined to stabilization. It suggested that no more distinct genes would be identified when the total clean reads reached a certain number. For all of the 12 libraries, there were more than 4.4 million clean reads (Supplementary Table 1), which indicated that the deep sequencing results were comprehensive and saturated.

3.5. DEGs between Groups

All genes annotated to the rat genome (Supplementary Excel 1) were analyzed for an evidence of differential expression. A detailed description of DEGs between two groups was presented in Supplementary Excel 2 (Con versus A/R), Excel 3 (A/R versus DZ), and Excel 4 (DZ versus DZ5HD). Those genes were to some extent differentially expressed; they were considered significant with a value more than 0.8. A list of the top 10 DEGs between two groups was shown in Table 2. In these genes, Mt-nd6, Acadl, and Idh2 are energy metabolism correlated and their expression status is listed in Table 3.

3.6. RT-qPCR Analysis

To confirm the DEGs revealed by the Illumina sequencing, 25 genes were randomly selected (Table 4) and assayed by SYBR green based RT-qPCR (Figure 3). Except Cycs, Idh3B, Mgst3, and Pdk4, 21 out of the 25 genes were expressed well in accordance with the results from Illumina sequencing (Table 4).

3.7. GO Enrichment Analysis

Ontology and term enrichment of DEGs in GO is listed in Figure 4. GO enrichment showed many of the DEGs from Con versus A/R (Figure 4(a)), A/R versus DZ (Figure 4(b)), and DZ versus DZ5HD (Figure 4(c)) participating in the Biological Process ontology. Histogram presentation of Gene Ontology functional classification and DEGs’ enrichment is shown in Figures 4(d)4(f) ((d) Con versus A/R; (e) A/R versus DZ; (f) DZ versus DZ5HD). The significantly enriched (corrected ) terms for Con versus A/R (Supplementary Table 2), A/R versus DZ (Supplementary Table 3), and DZ versus DZ5HD (Supplementary Table 4) also listed.

3.8. Pathway Analysis

KEGG Pathway provides an indication of the main biochemical and signal transduction pathways that DEGs are involved in. Pathway enrichment for Con versus A/R, A/R versus DZ, and DZ versus DZ5HD were displayed in Supplementary Tables 5–7. For Con versus A/R, A/R versus DZ, and DZ versus DZ5HD, there were 40, 48, and 37 pathways highly enriched (, ), respectively.

In all the pathways, Metabolic Process was the DEGs most enriched one. It is not difficult to notice that many energy metabolism correlated pathways, such as fatty acid metabolism pathway, TCA cycle, proteasome, PPAR signaling pathway, and peroxisome pathway were highly enriched (Table 5).

3.9. ATP Detection

At the end of reoxygenation, the ATP levels of the A/R groups were much lower than Con (); ATP concentrations of DZ groups were much higher than A/R groups (). When 5-HD was applied, the beneficial effect of DZ was abolished in DZ5HD when compared with DZ () (Figure 5).

4. Discussion

MIRI is always one of the leading causes of morbidity. In vitro experimental A/R model is a powerful tool to mimic in vivo ischemia-reperfusion injury. We developed an A/R model using freshly isolated adult rat cardiomyocytes, which are more relevant to the in vivo IR conditions.

In present study, DGE and bioinformatics technologies were employed to analyse molecular change after A/R and mitoKATP opening. Our results demonstrated not only the robustness of next-generation sequencing in exploring the molecular change resulting from mitoKATP opening but also the potential of the combine of next-generation sequencing and KEGG Pathway analysis to provide clues into target finding of molecular mechanisms underlying the myocardium protective effect of mitoKATP.

MitoKATP opening or closing in cultured adult rat cardiomyocytes significantly resulted in gene expression change. Many of the genes were energy related. Metabolic Process was the DEGs most enriched GO ontology and energy metabolism correlated pathways were highly enriched too. We could not help doubting that mitoKATP might have interfered with the energy metabolism and we confirmed that by directly measuring ATP content of four groups at the end of reoxygenation.

Three energy metabolism correlated genes, Mt-nd6, Idh2, and Acadl, were all upregulated (A/R versus DZ). Mt-nd6 encodes NADH-quinone oxidoreductase (complex I) subunit 6 in mammal. In the respiratory chain, complex I is responsible for the oxidation of NADH and contributes to the formation of the proton gradient which drives ATP synthesis and passes electrons to ubiquinone [31]. Ischemia-reperfusion injury was characterized by decreased complex I respiration [32]. In this study, expression of Mt-nd6 decreased after A/R treatment, while it was upregulated tremendously in DZ compared with A/R. Complex I is extremely susceptible to oxidative damage and subsequently produces more ROS [33], leading to extensive mitochondrial dysfunction and the depletion of ATP. MitoKATP opening by DZ increased Mt-nd6 expression, which might have contributed to ATP synthesis and resulted in its myocardial protection.

Idh2 encodes isocitrate dehydrogenase in mitochondria. In present study, expression of Idh2 varyed: Con versus A/R downregulated; A/R versus DZ upregulated; DZ versus DZ5HD downregulated. Isocitrate dehydrogenase is the rate-limiting enzyme of TCA cycle, which catalyzes the oxidative decarboxylation of isocitrate to 2-oxoglutarate. Isocitrate dehydrogenase plays a role in intermediary metabolism and energy production. It had been reported that isocitrate dehydrogenase activity increased at the ischemia region when heart underwent ischemia [34, 35]; the authors deemed the increase that came from the increased need of energy.

Long Chain Acyl-CoA Dehydrogenase (Acadl) encodes long chain acyl-CoA dehydrogenase, which catalyzes the α- and β-dehydrogenation of acyl-CoA esters in fatty acid metabolism. It is the first rate-limiting enzyme in fatty acid β-oxidation reaction [36]. In present study, Acadl was downregulated after A/R, upregulated when mitoKATP opened, and downregulated again when mitoKATP was blocked by 5-HD. In physiological state, 60–70% of the total energy the heart needs comes from fatty acid β-oxidation [37]. In the ischemic condition, FAO is more indispensable. Ito et al. [38] demonstrated that high levels of fatty acids in the perfusate were capable of enhancing postischemic energy production and increasing contractile function. That study provided evidence that, in heart with limited oxidative capacity, increasing exogenous energy substrate supply and boosting FAO generated more ATP and quickly normalized energy production. From what is mentioned above, mitoKATP opening may alleviate the energy depletion when adult cardiomyocytes underwent A/R injury by boosting the fatty acid β-oxidation.

6 energy correlated pathways, Peroxisome pathway, PPAR signaling pathway, citrate cycle (TCA cycle) pathway, fatty acid metabolism pathway, and proteasome pathway were DEGs significantly enriched ().

TCA cycle and fatty acid metabolism directly generate energy. 7 DEGs from Con versus A/R were enriched in TCA cycle pathway. They were all downregulated after A/R injury. When mitoKATP was open, all of them were upregulated. It is obvious that A/R suppressed TCA cycle. This could be one of the reasons why A/R decreased ATP content. We could see that DZ saved TCA cycle. 8 DEGs from Con versus A/R and 14 DEGs from A/R versus DZ, including Acadl, were enriched in fatty acid metabolism pathway. It seemed that A/R suppressed these two pathways and mitoKATP reinforced them.

Peroxisome proliferator-activated receptors (PPARs), especially PPAR-, are sensitive to fatty acids and their derivatives. They are also ligand-activated transcription factors regulating cardiac FAO and energy homeostasis [39, 40]. PPAR- is expressed highly in the heart and evidence had showed that PPAR- was involved in the regulation of numerous genes encoding FAO enzymes [41]. Overexpression of PPAR- and its target metabolic genes promoted FAO as a source of energy under conditions of acute IR [42, 43]. Besides its well-known action on cardiac energy metabolism and lipid homeostasis, emerging evidence indicated that administration of PPAR- synthetic ligands was myocardial protective in an IR setting, as manifested by improved postischemic recovery of contractile function and reduced infarct size in both in vivo and ex vivo models [42, 44]. Mice overexpressing PPAR- in heart displayed increased FAO rates, accumulated triacylglycerides, and decreased glucose metabolism, and they eventually developed cardiomyopathy [45, 46]. Not surprisingly, mice lacking PPAR- had elevated free fatty acid levels as a consequence of inadequate FAO, rendering them hypoglycemic as a result of their reliance on glucose [47]. In present study, although PPAR- gene did not change after mitoKATP opening, 17 DEGs from A/R versus DZ were enriched in peroxisome pathway () and 16 DEGs significantly enriched in PPAR signaling pathway (). 12 DEGs from DZ versus DZ5HD were enriched in () peroxisome pathway and 14 enriched in PPAR signaling pathway (). Nevertheless, to test the relative activity change of PPAR- after mitoKATP opening, further studies are needed.

The ubiquitin proteasome system (UPS) degrades targeted abnormal and most normal proteins in cells. Most degradation via the UPS is ATP-dependent. This process involves ubiquitin ligases E1, E2, and E3, which function in concert with chaperones to identify and ubiquitinate appropriate target proteins [4850]. Then the resulting polyubiquitinated proteins are transferred to the 26S proteasome, where they are degraded into peptides and ubiquitin. Proteasome pathway enriched 14 DEGs from Con versus A/R; 13 DEGs were downregulated in A/R group while 7 DEGs upregulated after mitoKATP opening. It is obvious that A/R induced the downregulation of UPS and mitoKATP opening reactivated it. Proteasome that functioned insufficiently had been observed most consistently in MIRI [51, 52]. Such studies supported the hypothesis that IR decreased proteasome activity by reducing ATP levels, as well as oxidative unfolding and damaging proteasome proteins [53]. To test this hypothesis, proteasome gain-of-function or loss-of-function studies in animal models of MIRI were carried out. However, the results showed a paradox: gain-of-function using transgenic mice with increased proteasome activity showed protection from MIRI [54], whereas loss-of-function studies using pharmacological means also revealed protection from MIRI [5557]. In present study, 13 DEGs were downregulated after A/R (Con versus A/R); this should be a feedback of ATP depletion resulting from A/R. 7 DGEs were upregulated in DZ group (A/R versus DZ); this could be a consequence of ATP recovery after mitoKATP opening.

Energy was so desperately needed in A/R environment that Metabolic Process was the most enriched GO ontology in Con versus A/R, A/R versus DZ, and DZ versus DZ5HD. In addition, energy metabolism related genes and pathways were significantly interfered with each other (Figure 6). UPS is protein related; TCA cycle pathway controls aerobic metabolism of glucose, PPAR-α, and Acadl effect on β-dehydrogenation of acyl-CoA esters in FAO. To sum up, mitoKATP regulated the metabolism of 3 main nutriments: glucose, fatty acid, and protein and kept a balance between energy production and consumption at the setting of A/R in adult cardiomyocytes. Strategies to increase energy supply in MIRI may be a good choice. Metabolism correlated genes and pathway nodes may be promising therapeutic targets. At the same time, we must confess that, to assure the effects of specific gene and signaling pathway mentioned above in MIRI, further gain- or/and loss-of-function studies will be needed.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Song Cao and Yun Liu contributed equally to this study.

Acknowledgment

This work was supported by grants from Special Scientific Research Fund for Public Welfare, Ministry of Health of China (Grant no. 200802173).

Supplementary Materials

Excel 1: Raw data of the expression of all the gene (4 groups × 3 replicates of the sequencing) .

Excel 2: Significant differently expressed genes from Control versus A/R.

Excel 3: Significant differently expressed genes from A/R versus DZ.

Excel 4: Significant differently expressed genes from DZ versus DZ5HD.

Supplementary Figure 1: Assessment of the DGE saturation degree.

Supplementary Table 1: The raw data of DGE profile of 12 sequencings.

Supplementary Table 2: DEGs significantly enriched GO terms from Control vs A/R (P<0.05).

Supplementary Table 3: DEGs significantly enriched GO terms from A/R vs DZ (P<0.05).

Supplementary Table 4: DEGs significantly enriched GO terms from DZ vs DZ5HD (P<0.05).

Supplementary Table 5: Significantly enriched pathways for DEGs from Con vs A/R (P<0.01).

Supplementary Table 6: Significantly enriched pathways for DEGs from A/R vs DZ (P<0.01).

Supplementary Table 7: Significantly enriched pathways for DEGs from DZ vs DZ5HD (P<0.01).

  1. Supplementary Materials