Differentially expressed miRNAs in sepsis-induced acute kidney injury target oxidative stress and mitochondrial dysfunction pathways

Abstract

Objective

To identify specific miRNAs involved in sepsis-induced AKI and to explore their targeting pathways.

Methods

The expression profiles of miRNAs in serum from patients with sepsis-induced AKI (n = 6), sepsis-non AKI (n = 6), and healthy volunteers (n = 3) were investigated by microarray assay and validated by quantitative PCR (qPCR). The targets of the differentially expressed miRNAs were predicted by Target Scan, mirbase and Miranda. Then the significant functions and involvement in signaling pathways of gene ontology (GO) and KEGG pathways were analyzed. Furthermore, eight miRNAs were randomly selected out of the differentially expressed miRNAs for further testing by qPCR.

Results

qPCR analysis confirmed that the expressions levels of hsa-miR-23a-3p, hsa-miR-4456, hsa-miR-142-5p, hsa-miR-22-3p and hsa-miR-191-5p were significantly lower in patients with sepsis compared with the healthy volunteers, while hsa-miR-4270, hsa-miR-4321, hsa-miR-3165 were higher in the sepsis patients. Statistically, miR-4321; miR-4270 were significantly upregulated in the sepsis-induced AKI compared with sepsis-non AKI, while only miR-4321 significantly overexpressed in the sepsis groups compared with control groups. GO analysis showed that biological processes regulated by the predicted target genes included diverse terms. They were related to kidney development, regulation of nitrogen compound metabolic process, regulation of cellular metabolic process, cellular response to oxidative stress, mitochondrial outer membrane permeabilization, etc. Pathway analysis showed that several significant pathways of the predicted target genes related to oxidative stress. miR-4321 was involved in regulating AKT1, mTOR and NOX5 expression while miR-4270 was involved in regulating PPARGC1A, AKT3, NOX5, PIK3C3, WNT1 expression. Function and pathway analysis highlighted the possible involvement of miRNA-deregulated mRNAs in oxidative stress and mitochondrial dysfunction.

Conclusion

This study might help to improve understanding of the relationship between serum miRNAs and sepsis-induced AKI, and laid an important foundation for further identification of the potential mechanisms of sepsis-induced AKI and oxidative stress and mitochondrial dysfunction.

About the Authors:

Qin-Min Ge

Affiliation: Department of Emergency, Xinhua Hospital, Shanghai Jiaotong University School of Medicine, Shanghai, China

Chun-Mei Huang

Affiliation: Department of Emergency, Xinhua Hospital, Shanghai Jiaotong University School of Medicine, Shanghai, China

Xiang-Yang Zhu

Affiliation: Division of Nephrology and Hypertension, Department of Internal Medicine, Mayo Clinic College of Medicine, Rochester, MN, United States of America

Fan Bian

Affiliation: Department of Nephrology, Xinhua Hospital, Shanghai Jiaotong University School of Medicine, Shanghai, China

Shu-Ming Pan

* E-mail: [email protected]

Affiliation: Department of Emergency, Xinhua Hospital, Shanghai Jiaotong University School of Medicine, Shanghai, China

Introduction

Despite advances in the treatment of critically ill patients, the development of acute kidney injury (AKI) shows a high mortality rate [1]. The mortality rate was 34% in patients with AKI versus 7% in patients without AKI [2]. The cause of AKI in critically ill patients is often multifactorial, while sepsis is the most common cause which is up to 50% [3, 4]. And sepsis-induced AKI is diagnosed in up to 47.9% of ICU patients, it is associated with increased progression to chronic kidney disease (CKD) [1]. Therefore sepsis-induced AKI should be intervened as early as possible. But the mechanisms underlying this event are not fully understood.

There were distinct pathophysiological changes in sepsis-induced AKI compared to non sepsis-induced AKI reflected by special biomarkers, including interleukin-18 (IL-18) whose concentrations were higher in urine of sepsis induced AKI patients [5] Various serum and urinary biomarkers have been found to be up-regulated in kidney injury, such as liver fatty-acid binding protein, cystatin C, kidney injury molecule-1 and neutrophil gelatinase-associated lipocalin [6]. However, they are not specific enough to reflect the pathophysiological changes during early stage of sepsis and AKI.

MicroRNAs (miRNA) are endogenous non-coding RNA molecules which regulate gene expression by degrading mRNA or inhibiting translation. Recent studies demonstrated circulating miRNAs had potential value in diagnosis and prognosis of sepsis [7, 8]. Changes in miRNA expression have been implicated in disease state and reflected the pathophysiological changes during early stage of sepsis when compared with other biomarkers [7, 9]. Circulating miRNAs detected in body fluids are novel biomarkers of renal diseases, including AKI [10, 11]. However, little is known about miRNA expression in septic AKI patients let alone its role in mitochondrial oxidative stress and endothelial dysfunction.

In this study, we profiled serum miRNAs expression levels in sepsis-induced AKI and sepsis-non AKI affected by Gram negative (G-) bacteria, identified the differentially expressed miRNAs in sepsis-induced AKI. In addition, we investigated the miRNA expression profiles to confirm the hypothesis that miRNAs alteration in serum might be involved in the mitochondrial oxidative stress and endothelial dysfunction during sepsis-induced AKI.

Materials and methods

Patients and blood samples collection

This study was approved by Shanghai Jiaotong University Xinhua Hospital Ethics Committee (XHEC-D-2017-006) and was carried out in accordance with the Declaration of Helsinki. Written informed consents were obtained from all participate in this study. Patients with sepsis were consecutively enrolled from the Xin-Hua Hospital affiliated with the Shanghai Jiaotong University School of Medicine, China, from May 2014 to January 2015. Patients were divided into two groups: the septic-AKI group (meeting the diagnosis of both sepsis and AKI, n = 35) and sepsis-non AKI group (with sepsis but without AKI; n = 30), AKI was defined with KDIGO criteria and sepsis with consensus criteria. Healthy subjects were recruited from the staff of hospital mentioned above. The patients with pre-existing chronic renal failure or chronic use of RRT were excluded from the study. The baseline characteristics and underlying diseases of patients with sepsis-induced AKI, sepsis-non AKI and healthy controls were shown in Table 1. All patients were evaluated with the peak C-reactive protein (CRP), procalcitonin (PCT) and the level of leukocytosis to determine the severity of inflammation. Additionally, the physiological and clinical data were collected, including those that compose the APACHEII (acute physiology and chronic health evaluation II) and SOFA (Sequential Organ Failure Assessment) score. To measure the levels of microRNAs and inflammatory cytokines, serum samples for the 1st day of diagnosis with SIRS or AKI were measured. The bacteria identified by blood culture in the patients with G- were Escherichia coli. There was no statistical significance in terms of age or gender among the groups.

[Figure omitted. See PDF.]

Table 1. Baseline characteristics.

https://doi.org/10.1371/journal.pone.0173292.t001

Plasma cytokines assay

Venous blood samples were taken at the time of admission for assessment of serum creatinine concentrations. Venous blood samples were obtained from each patient within 24 h from the first day of diagnosis with sepsis or AKI. The serum was separated, divided in aliquots, and immediately frozen at−80°C until the time of assay. The levels of cytokines including IL-6, IL-8, IL-2sR, IL-10, tumor necrosis factor (TNF)-α and IL-1β were quantified by using an IMMULITE 1000 Automated Immunoassay Analyzer in accordance with the manufacturer’s instructions.

Microarray analysis of mirna expression profile

15 blood samples, including 6 patients with G- sepsis-induced AKI, 6 patients with G- sepsis-non AKI and 3 healthy controls, randomly selected from the participants in Table 1, After blood collection, we separated serum, isolated total RNA from for miRNA and mRNA analysis. To detect the global miRNA expression profile of sepsis-induced AKI and sepsis-non AKI, we used the TRIzol Reagent (Invitrogen, USA) and a miRNeasy mini kit (Qiagen, Shanghai) to isolate total RNA according to the manufacturers’ instructions. The samples were labeled using the miRCURY™ Hy3™/Hy5™ Power labeling kit (Exiqon, Vedbaek, Denmark) and then hybridized on the miRCURY™ LNA Array (v.18.0) that contained 3100 capture probes. Following the washing steps the slides, the Axon GenePix 4000B microarray scanner (Axon Instruments, Foster City, CA) was used to scan the images. The data was analyzed in GenePix Pro 6.0 software (Axon) for grid alignment and data extraction.

After normalization, significant differentially expressed miRNAs were identified through Volcano Plot filtering between the two experimental groups. Finally, hierarchical clustering was performed to show distinguishable miRNA expression profiling among the samples. The differentially expressed miRNAs were identified by a fold change of ≥1.5.

Quantitative RT-PCR for miRNA expression

Total RNA was isolated from serum using TRIzol reagent (Invitrogen Life Technologies, USA). The RNA quality and quantity were measured using a NanoDrop ND-1000 spectrophotometer and then quantified using RT-PCR (16°C for 30 min, 42°C for 40 min, and 85°C for 5 min) and a Gene Amp PCR System 9700 (Applied Biosystems). Each RT reaction mixture contained 5 μl Master Mix, 0.5 μl each primer (10 μM), and 2 μl of the corresponding template cDNA, then was incubated at 95°C for 10 min, followed by 40 PCR cycles (95°C for 10 s and then 60°C for 60 s) on a ViiA 7 Real-time PCR System (Applied Biosystems). hsa-miR-423-5p level was used as the endogenous control. The bi-directional primer sequences are presented in S1 Table. The relative expression level for each miRNA was presented by using the comparative cycle threshold (CT) method (2−ΔΔCT).

Gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) pathway analysis of the deregulated miRNAs

The miRNA target gene database mirbase, miranda and targetscan were used to predict target genes of the differentially expressed miRNAs. Then the final targets were integrated from these three public databases. Using this identified target gene set, we analyzed gene ontology (GO) and KEGG pathway of the significant functions and its involvement in signaling pathways. From the Biological Process, Cellular Component and Molecular Function domains, GO analyze the main function of the predicted target genes and uncover the miRNA-gene regulatory network respectively. Pathway analysis showed that a significant enrichment in several pathways of the predicted target genes. The significance of GO term enrichment and the pathways analysis was defined by p value and FDR (a p-value<0.05 is recommended).

Statistical analysis

Data were presented as mean ± SD. The differences in levels of miRNAs among groups were analyzed using two-tailed t-tests or one-way ANOVA. A P value <0.05 was considered statistically significant.

Results

In the present study, 65 patients were diagnosed with sepsis, 53.8% patients developed septic-AKI. The main sources of sepsis were as follows: pulmonary (n = 20; 30.7%), urogenital (n = 15; 23.08%), gastrointestinal (n = 20; 30.7%), and others (n = 10; 15.4%). The patients’ Clinical characteristics and laboratory data were presented in Tables 1 and 2. The mean age of total 65 patients was 66.7±9.9 years and 40 patients (61.5%) were female. The characteristics of the patients were well balanced between the two study groups. Higher level of serum Creatinine, BUN, CRP and PCT were found among the patients with sepsis who developed AKI. The APACHEII score, SOFA score were significantly elevated in patients with septic-AKI compared with those of patients with sepsis but without AKI. However, there were no significant difference in the levels of Bilirubin, Lactate, and 24h urine output as well as 28d in-hospital mortality between the two groups.

[Figure omitted. See PDF.]

Table 2. Clinical characteristics of the septic-acute kidney injury (AKI) and the sepsis non-AKI groups.

https://doi.org/10.1371/journal.pone.0173292.t002

Comparison of the level of cytokines in the septic AKI group with those in the sepsis-non AKI and healthy control groups

To further investigate whether inflammation status was different between septic AKI patients and sepsis non-AKI patients, the levels of pro-inflammatory and anti-inflammatory cytokines as well as IL-2sR were studied. As mentioned, the level of cytokines was also compared to exclude the effect of decreased renal excretory function between the septic-AKI and healthy controls groups. In patients with septic AKI, the levels of TNF-α, IL-1β, IL-6, IL-2sR and IL-10 were significantly increased, while the levels of IL-8 was significantly decreased compared with sepsis patients without AKI. There was no statistical significance on the levels of IL-1β, IL-6, IL-2sR and IL-10 in sepsis-non AKI group compared with the control (Table 3).

[Figure omitted. See PDF.]

Table 3. Changes of cytokines in plasma in different groups (mean ± SD).

https://doi.org/10.1371/journal.pone.0173292.t003

Differential expression of serum miRNAs from the different groups

To identify the most significant candidates, miRNAs with at least 1.5 fold expression changes were selected (S2 Table). A total of 40 miRNAs were differentially expressed between sepsis-induced AKI group and sepsis-non AKI groups (FDR < 0.05). In the sepsis-induced AKI group, 13 of the 40 miRNAs, under the criteria, were up-regulated and 27 miRNAs were down-regulated compared with the sepsis-non AKI group (Fig 1, Tables 4 and 5) Then we profiled the miRNAs among these two groups and the control. A total of 37 miRNAs were differentially expressed among the sepsis-induced AKI group, sepsis-non AKI group and the control (FDR < 0.05) (Fig 2). Among the 37 miRNAs, there were similar expression trends of 8 miRNAs in both the sepsis-induced AKI and sepsis-non AKI group compared with the control. These 8 miRNAs differentially expressed between the sepsis-induced AKI group and the sepsis-non AKI group were hsa-miR-142-5p, hsa-miR-4270, hsa-miR-22-3p, hsa-miR-191-5p, hsa-miR-3165, hsa-miR-4321, hsa-miR-23a-3p and hsa-miR-4456. Three miRNAs (hsa-miR-4270, hsa-miR-4321, hsa-miR-3165) were up-regulated in the sepsis-induced AKI group; and 5 miRNAs (hsa-miR-142-5p, hsa-miR-22-3p, hsa-miR-191-5p, hsa-miR-23a-3p, hsa-miR-4456) were down-regulated.

[Figure omitted. See PDF.]

Fig 1. Differentially expressed miRNAs between two groups.

Septic AKI (n = 6); Sepsis-non AKI (n = 6). Red color indicated high relative expression and green color denoted low relative expression. miRNA with expression fold change >1.5 and with FDR <0.05 was considered statistically significant.

https://doi.org/10.1371/journal.pone.0173292.g001

[Figure omitted. See PDF.]

Fig 2. Differentially expressed miRNAs among the three groups.

Control (n = 3); Septic AKI (n = 6); Sepsis-non AKI (n = 6). Red color indicated high relative expression and green color denoted low relative expression. miRNA with expression fold change >1.5 and with FDR <0.05 was considered statistically significant. MiRNAs marked with arrow were randomly selected for further confirmation by RT-qPCR.

https://doi.org/10.1371/journal.pone.0173292.g002

[Figure omitted. See PDF.]

Table 4. Differentially expressed (fold change ≥ 1.5) miRNAs in sepsis-induced AKI and sepsis-non AKI.

https://doi.org/10.1371/journal.pone.0173292.t004

[Figure omitted. See PDF.]

Table 5. Differentially expressed (fold change ≥ 1.5) miRNAs in sepsis-induced AKI and sepsis-non AKI.

https://doi.org/10.1371/journal.pone.0173292.t005

A scatter plot was generated to assess the quality of the miRNA data after filtering, and a volcano plot was generated to visualize the differential expression between two different conditions, as shown in Fig 3.

[Figure omitted. See PDF.]

Fig 3.

Volcano plot (A) and scatter plot (B) of the miRNA microarray analysis. (A) Volcano plots are useful tool for visualizing differential expression patterns between two different conditions. The vertical lines correspond to 1.5-fold up and down, respectively, and the horizontal line represents a P value of 0.05. So the red point in the plot represents the differentially expressed miRNAs with statistical significance. (B) The scatter plot is a useful visualization for assessing the variation (or reproducibility) between chips. The axes of the scatter plot are the normalized signal values of the samples (the ratio scale).

https://doi.org/10.1371/journal.pone.0173292.g003

Validation of microarray results by RT-qPCR

Based on the predicted target genes of these miRNAs, we selected eight miRNAs (hsa-miR-4270, hsa-miR-4321, hsa-miR-3165, hsa-miR-23a-3p, hsa-miR-4456, hsa-miR-142-5p, hsa-miR-22-3p, hsa-miR-191-5p) randomly from the differentially expressed miRNAs for further validation.

Relative expression of each miRNA was normalized to hsa-miR-423-5p level. Among the eight selected miRNAs, five miRNAs including hsa-miR-23a-3p, hsa-miR-4456, hsa-miR-142-5p, hsa-miR-22-3p and hsa-miR-191-5p were significantly down-regulated in both the sepsis-induced AKI and sepsis-non AKI groups compared with the healthy controls (P < 0.05); three miRNAs including hsa-miR-4270, hsa-miR-4321, hsa-miR-3165 were significantly upregulated in the sepsis-induced AKI group compared with the healthy controls (P < 0.05). There was no statistical significance between sepsis-induced AKI and sepsis-non AKI groups on down-regulated miRNAs expression (P >0.05). As to the three up-regulated miRNAs, the expression of miR-4270 was up-regulated only in AKI group (P < 0.05) with no statistical significance in sepsis-non AKI group compared with the control. The expression of miR-3165 was increased in sepsis groups no matter whether AKI existed or not (P < 0.05). While statistical significance among the sepsis-induced AKI, sepsis-non AKI and the control groups was shown only on miR-4321 expression (P <0.05) (Fig 4). The miRNAs whose expression levels altered were related to at least 77 mitochondrial oxidative stress and dysfunction-associated target genes based on target predictions (Table 6).

[Figure omitted. See PDF.]

Fig 4. Confirmation of miRNA level by RT-qPCR.

Relative expression of each miRNA was normalized to hsa-miR-423-5p level. miR-142-5p, miR-22-3p, miR-191-5p, miR-23a-3p and miR-4456 were significantly decreased in both the sepsis-induced AKI and sepsis-non AKI groups compared with the healthy controls. miR-4270, miR-4321 and miR-3165 were obviously increased in the sepsis-induced AKI, and miR-3165 also was obviously increased in sepsis-non AKI group compared with the healthy controls. While statistical significance between the sepsis-induced AKI and sepsis-non AKI groups was shown only on miR-4321 expression. Note: a single asterisk indicates significant difference between the healthy controls (n = 3) and the sepsis-induced AKI group (n = 6), as well as between the healthy controls and the sepsis-non AKI group (n = 6). A pound sign indicates a statistical difference between the sepsis-induced AKI group and the sepsis-non AKI group. P < 0.05 was considered statistically significant.

https://doi.org/10.1371/journal.pone.0173292.g004

[Figure omitted. See PDF.]

Table 6. The differentially expressed miRNAs (sepsis-induced AKI and sepsis-non AKI >1.5) and their predicted target genes.

https://doi.org/10.1371/journal.pone.0173292.t006

GO and KEGG pathway analysis of the deregulated miRNAs

We predicted the potential target genes of these 8 differentially expressed miRNAs among the three groups to further explore the miRNAs function by using Target Scan, mirbase and Miranda. Then the significant functions and involvement in signaling pathways of GO and KEGG pathways were analyzed. GO analysis showed that regulation of cellular metabolic process, regulation of gene expression and regulation of nitrogen compound metabolic process represented the significantly enriched ones. (Fig 5A, S3 Table) Other biological processes regulated by the predicted target genes included kidney development, cellular response to oxidative stress, mitochondrial outer membrane permeabilization and etc. Pathway analysis showed that Wnt signaling pathway and FoxO signaling pathway were abundant among the significantly enriched ones. Other significant pathways of the predicted target genes included HIF-1 signaling pathway, PI3K-Akt signaling pathway, mTOR signaling pathway, TNF signaling pathway, TGF-beta signaling pathway (Fig 5B, S4 Table). They were found to be involved in response to oxidative stress and mitochondrial dysfunction (Fig 6). From the functional analysis of these DE miRNAs using KEGG pathway enrichment analysis, we inferred that circulating miRNAs might participate in septic AKI by influencing immune and inflammatory response indirectly. All target genes of the eight differentially expressed miRNAs in the three databases are shown in S5 Table. In addition, miRNA-mRNA gene network analysis also integrated all the miRNAs in the microarray and GOs by outlining the interactions of miRNA and GO-related genes (Fig 7).

[Figure omitted. See PDF.]

Fig 5. Bioinformatics GO and pathway analysis based on miRNA target genes.

Enrichment score is equal to -log10 (P-value) that represents the significant level of GOs and pathways. (A) Significant GOs. (B) Significant signaling pathways. KEGG: Kyoto Encyclopedia of Genes and Genomes; GO: Gene ontology. BH-corrected P < 0.05 was considered statistically significant.

https://doi.org/10.1371/journal.pone.0173292.g005

[Figure omitted. See PDF.]

Fig 6.

The pathways enriched in the FoxO signaling pathway (A), PI3K-Akt signaling pathway(B), Wnt signaling pathway (C) and mTOR signaling pathway.

https://doi.org/10.1371/journal.pone.0173292.g006

[Figure omitted. See PDF.]

Fig 7. miRNA-mRNA gene network analysis.

All the miRNAs in the microarray and GOs were integrated by outlining the interactions of miRNA and GO-related genes

https://doi.org/10.1371/journal.pone.0173292.g007

Discussion

Many recent clinical and experimental studies have shown that sepsis-induced AKI is characterized by endothelial injury with hemodynamic dysfunction [12, 13]. Oxidative stress, a state of reactive oxygen species (ROS) overproduction plays a key role in the pathogenesis of endothelial dysfunction. Release of inflammatory mediators, oxidative stress together with mitochondrial respiratory chain dysfunction can contribute to AKI [14, 15]. Sepsis is the most common cause of acute kidney injury in clinical practice. Inflammation and cytokines play key roles in sepsis pathophysiology [16]. Sepsis-induced AKI was considered as an integrated response, including transcriptional events, ROS signaling, mitochondrial activity and metabolic orientation such as apoptosis [17]. Despite the evidence that mitochondrial oxidative stress and dysfunction are implicated in the development of sepsis induced kidney injury, the precise pathophysiological mechanisms remain unclear.

Emerging evidence demonstrated that miRNAs are intensively engaged in mitochondrial oxidative stress and dysfunction [18]. Present studies demonstrated circulating miRNAs had potential value in diagnosis and prognosis of sepsis. Changes in miRNA expression have been implicated in disease states and reflected the pathophysiological changes during early stage of either sepsis or AKI when compared with other biomarkers. Some serum miRNAs, including miR-21, miR-494, miR-21, miR-155, miR-210, miR-101-1, miR-127-3p, miR-126, miR-26b, miR-29a, miR 146a, miR-27a, miR-93 and miR-10a were identified as the biomarkers for the diagnosis of AKI [19–21], but miRNAs in relation to septic AKI patients have not been previously reported yet.

In this study, we observed that hsa-miR-508-5p was the most significantly up-regulated, whereas, ebv-miR-BART1-3p was the most significantly down-regulated miRNAs in the sepsis-induced AKI group than healthy controls respectively; hsa-miR-508-5p was the most increased and ebv-miR-BART1-3p was the most decreased miRNAs in the sepsis-non AKI group than with the healthy controls respectively; hsa-miR-5094 was the most up-regulated, whereas, hsa-miR-326 was the most down-regulated miRNAs in the sepsis-induced AKI group than the sepsis-non AKI group. Many deregulated miRNAs in our study, for example, miR-222-3p and miR-221-3p have been reported to be involved in specific signaling pathways. miR-222-3p and miR-221-3p were found to promote intracellular ROS accumulation in Human Aortic Endothelial Cells (HAECs) by targeting PGC-1α [22]. However, the biological functions of most differentially expressed serum miRNAs, such as miR-4321 is still largely unknown. In this study, we found that a different subset of miRNAs was dysregulated in sepsis-induced AKI compared to reports of sepsis-non AKI patients. We identified 37 miRNAs differentially expressed among the three groups. While in our work, the majority of the deregulated miRNAs exhibited similar expression trend in sepsis groups compared with the healthy controls, which demonstrated that both sepsis-induced kidney injury and sepsis-non AKI seem to share, at least in part similar regulatory mechanisms. To identify serum miRNAs that might be correlated with resistance or susceptibility to sepsis-induced kidney injury, we focused on the miRNAs differentially expressed between the sepsis-induced kidney injury group and the sepsis-non AKI group. In the present study, we found that 49 miRNAs were differentially expressed among the three groups, of which 13 miRNAs were increased, whereas 36 miRNAs were decreased in the sepsis-induced kidney injury group compared with the sepsis-non AKI group. According to the KEGG pathway and GO analysis regarding target genes of these deregulated miRNAs, we found significantly enriched functions and signaling pathways were associated with oxidative stress and mitochondrial dysfunction. GO analysis showed that regulation of cellular metabolic process, regulation of gene expression and regulation of nitrogen compound metabolic process represented the significantly enriched ones. Other biological processes regulated by the predicted target genes, including kidney development, cellular response to oxidative stress, mitochondrial outer membrane permeabilization, etc. Pathway analysis showed that Wnt signaling pathway and FoxO signaling pathway were abundant among the significantly enriched ones. Other significant pathways of the predicted target genes related to oxidative stress, including HIF-1 signaling pathway, PI3K-Akt signaling pathway, mTOR signaling pathway, TNF signaling pathway, TGF-beta signaling pathway, etc. Function and pathway analysis highlighted the possible involvement of miRNA-deregulated mRNAs in mitochondrial oxidative stress and dysfunction.

To explore the effects of miRNAs on mitochondrial oxidative stress and dysfunction related target gene, we compared the results of miRNA and mRNA. As shown above, an individual miRNA may target multiple mRNAs. In our study, we found that the variation of the 8 miRNA expression were associated with 13 genes which were involved in the activation of mitochondrial oxidative stress and dysfunction response (PPARGC1A (PGC-1α), SIRT1, OXSR1, FOXO1), inflammation (IGF2BP3, IL8, IL6, PDCD4 and TNFAIP6). Notably, 5 miRNAs were involved in PGC-1α regulation; 4 miRNAs were involved in SIRT1 regulation. PGC-1α is a key regulator of mitochondrial biogenesis, the levels of PGC-1α are deregulated when there is damage to the mitochondria during sepsis [23, 24]. SIRT1, upstream of PGC-1α, is also an important modulator of energy metabolism and contributes to the mitochondrial biogenesis program. The SIRT1/PGC-1α axis was involved in mitochondrial biogenesis [25]. A few miRNAs have been found to be associated with mitochondrial oxidative stress and dysfunction in various organs. Overexpression of miR-23a was also found to damage mitochondria in cultured cardiomyocytes, ascribed to PGC-1α downregulation [26]. Overexpression of miR-142 led to a decrease in SIRT1 expression in the brain, and SIRT1 is a direct miR-142-5p target [27]. Our results showed that miR-142-5p, miR-191-5p, miR-3165, miR-4270 and miR-23a-3p were involved in regulating PGC-1α expression; miR-23a-3p, miR-22-3p, miR-4456 were involved in regulating SIRT1 expression. It was worth noting that difference between the sepsis-induced AKI and sepsis-non AKI groups was shown only on miR-4321; miR-4270 was only significantly overexpressed in the sepsis-induced AKI group compared with the sepsis and control groups. Our results indicated that miR-4321 was involved in regulating AKT1, mTOR and NOX5 expression while miR-4270 was involved in regulating PGC-1α, AKT3, NOX5, PIK3C3, WNT1 expression. Acute loss of renal function reduces leukocyte recruitment into inflamed tissue by reducing phosphorylation of Akt and p38 MAPK and transmigration [28]. mTOR is a key kinase downstream of PI3K/AKT, which regulates cell proliferation, growth, survival and angiogenesis. In our sepsis-induced AKI, many of the altered miRNAs were down-regulated which might decrease the expression of PGC-1α directly or indirectly, and weaken the mitochondrial oxidative stress and dysfunction.

Recent evidence from laboratory and clinical studies indicated that AKI in sepsis had a prominent inflammatory component in both the initiation and the extension phases of the kidney injury [29, 30]. The onset of AKI was preceded by early and remarkable inflammatory response and oxidative stress [31]. Inflammatory cytokines, including TNF, interleukin-1β (IL-1β) and IL-6, can induce oxidative stress, stimulate PGC-1α mRNA expression [22]. For a more comprehensive cytokine profile in sepsis, we measured the levels of pro-inflammatory and anti-inflammatory cytokines as well as IL-2sR. Our data showed that pro-inflammatory cytokines, such as TNF-α and IL-8 were significantly increased in the sepsis-non AKI compared with the healthy controls. With the aggravation of the disease, the levels of pro-inflammatory cytokine IL-8 was significantly decreased in the septic AKI compared with sepsis-non AKI group, but still higher than the normal control group; and the levels of pro-inflammatory cytokines TNF-α, IL-1β and IL-6 gradually increased, significantly higher than the normal control and sepsis-non AKI groups; anti-inflammatory cytokines IL-10 was significantly increased in the septic AKI compared with sepsis-non AKI and control group. IL-2sR was significantly increased in patients with septic AKI compared with those of either sepsis-non AKI or healthy controls patients. CD25 is an IL-2 receptor, sCD25 or IL-10 may be useful as a novel biomarker for the development of septic AKI [32]. Several large cohorts of critically ill patients demonstrated that IL-6 could be a robust predictor of AKI [33]. In keeping with our findings, Murugan et al [34] demonstrated that septic AKI patients had higher plasma levels of IL-6 and TNF-α than non-AKI patients. Our study extends the clinical observations by demonstrating the temporal relation between early, exaggerated inflammatory and oxidative stress response and AKI, thus supporting the cause-effect interplay between AKI and abnormal immune response.

To clarify the mechanism of septic AKI occurrence and evaluate the most significant candidates, miRNAs and their possible target genes, which were in the intersection of protein binding and FoxO signaling pathway and Wnt signaling pathway, might be the focus of the further research.

To summarize, we identified specific miRNAs among the sepsis-induced kidney injury group, the sepsis-non AKI group and the healthy controls. The different expression of miRNA might potentially be used to discriminate septic AKI from sepsis-non AKI and they were correlated with the regulation of mitochondrial oxidative stress and dysfunction, including PGC-1α, SIRT1, mTOR, OXSR1 and NOX5. Further study is needed to confirm the relationship between dysregulated miRNAs and their targets, together with the clinical significance of our observations. The interactions between miRNAs and their target genes are complex. The modulation of the oxidative stress and mitochondrial dysfunction by miRNAs is complex as well and requires further investigation.

Supporting information

[Figure omitted. See PDF.]

S1 Table. The bi-directional primer sequences for differentially expressed miRNAs.

https://doi.org/10.1371/journal.pone.0173292.s001

(XLSX)

S2 Table. Differentially expressed circulating miRNAs in septic AKI and sepsis-non AKI group.

https://doi.org/10.1371/journal.pone.0173292.s002

(XLSX)

S3 Table. Functional classification of the target genes by GO analysis.

https://doi.org/10.1371/journal.pone.0173292.s003

(XLS)

S4 Table. Pathway Enrichment analysis.

https://doi.org/10.1371/journal.pone.0173292.s004

(XLS)

S5 Table. All target genes of the eight differentially expressed miRNAs in the three databases.

https://doi.org/10.1371/journal.pone.0173292.s005

(XLSX)

Acknowledgments

The microarray experiments were performed by KangChen Bio-tech, Shanghai, China. This study was supported by National Natural Science Foundation of China (No.81000354), Shanghai Nature Science Foundation (14ZR1427000) and 2013–2014 National clinical key specialty construction project. This study greatly acknowledges the support from the Emergency nursing team of Xinhua Hospital, Shanghai Jiaotong University School of Medicine (Shanghai, China) for assistance with sample collection. The microarray experiments were performed by KangChen Bio-tech, Shanghai, China.

Author Contributions

1. Conceptualization: QMG SMP.

2. Data curation: CMH.

3. Formal analysis: QMG FB.

4. Funding acquisition: QMG SMP.

5. Investigation: QMG CMH FB.

6. Methodology: QMG XYZ.

7. Project administration: QMG SMP.

8. Resources: QMG CMH FB.

9. Software: FB.

10. Supervision: XYZ SMP.

11. Validation: QMG CMH.

12. Visualization: XYZ FB.

13. Writing – original draft: QMG CMH.

14. Writing – review & editing: FB SMP.

Citation: Ge Q-M, Huang C-M, Zhu X-Y, Bian F, Pan S-M (2017) Differentially expressed miRNAs in sepsis-induced acute kidney injury target oxidative stress and mitochondrial dysfunction pathways. PLoS ONE 12(3): e0173292. https://doi.org/10.1371/journal.pone.0173292

References

1. Keir I, Kellum JA. Acute kidney injury in severe sepsis: pathophysiology, diagnosis, and treatment recommendations. Journal of veterinary emergency and critical care (San Antonio, Tex: 2001). 2015;25(2):200–9. Epub 2015/04/08.

2. Chertow GM, Burdick E, Honour M, Bonventre JV, Bates DW. Acute kidney injury, mortality, length of stay, and costs in hospitalized patients. Journal of the American Society of Nephrology: JASN. 2005;16(11):3365–70. Epub 2005/09/24. pmid:16177006

3. Uchino S, Kellum JA, Bellomo R, Doig GS, Morimatsu H, Morgera S, et al. Acute renal failure in critically ill patients: a multinational, multicenter study. Jama. 2005;294(7):813–8. Epub 2005/08/18. pmid:16106006

4. Bagshaw SM, Uchino S, Bellomo R, Morimatsu H, Morgera S, Schetz M, et al. Septic acute kidney injury in critically ill patients: clinical characteristics and outcomes. Clinical journal of the American Society of Nephrology: CJASN. 2007;2(3):431–9. Epub 2007/08/19. pmid:17699448

5. Bagshaw SM, Langenberg C, Haase M, Wan L, May CN, Bellomo R. Urinary biomarkers in septic acute kidney injury. Intensive care medicine. 2007;33(7):1285–96. Epub 2007/05/10. pmid:17487471

6. Wan L, Bagshaw SM, Langenberg C, Saotome T, May C, Bellomo R. Pathophysiology of septic acute kidney injury: what do we really know? Critical care medicine. 2008;36(4 Suppl):S198–203. Epub 2008/04/11.

7. Essandoh K, Fan GC. Role of extracellular and intracellular microRNAs in sepsis. Biochimica et biophysica acta. 2014;1842(11):2155–62. Epub 2014/08/03. PubMed Central PMCID: PMCPMC4188718. pmid:25086335

8. Reid G, Kirschner MB, van Zandwijk N. Circulating microRNAs: Association with disease and potential use as biomarkers. Critical reviews in oncology/hematology. 2011;80(2):193–208. Epub 2010/12/15. pmid:21145252

9. Ardekani AM, Naeini MM. The Role of MicroRNAs in Human Diseases. Avicenna journal of medical biotechnology. 2010;2(4):161–79. Epub 2010/10/01. PubMed Central PMCID: PMCPMC3558168. pmid:23407304

10. Aguado-Fraile E, Ramos E, Conde E, Rodriguez M, Liano F, Garcia-Bermejo ML. MicroRNAs in the kidney: novel biomarkers of acute kidney injury. Nefrologia: publicacion oficial de la Sociedad Espanola Nefrologia. 2013;33(6):826–34. Epub 2013/10/26.

11. Bellinger MA, Bean JS, Rader MA, Heinz-Taheny KM, Nunes JS, Haas JV, et al. Concordant changes of plasma and kidney microRNA in the early stages of acute kidney injury: time course in a mouse model of bilateral renal ischemia-reperfusion. PloS one. 2014;9(4):e93297. Epub 2014/04/04. PubMed Central PMCID: PMCPMC3973691. pmid:24695114

12. Basile DP, Yoder MC. Renal endothelial dysfunction in acute kidney ischemia reperfusion injury. Cardiovascular & hematological disorders drug targets. 2014;14(1):3–14. Epub 2014/08/05. PubMed Central PMCID: PMCPMC4215733.

13. Pinto CF, Watanabe M, da Fonseca CD, Ogata CI, Vattimo Mde F. [The sepsis as cause of acute kidney injury: an experimental model]. Revista da Escola de Enfermagem da U S P. 2012;46 Spec No:86–90. Epub 2013/03/05. pmid:23250263

14. Tran M, Parikh SM. Mitochondrial biogenesis in the acutely injured kidney. Nephron Clinical practice. 2014;127(1–4):42–5. Epub 2014/10/25. pmid:25343819

15. Parikh SM, Yang Y, He L, Tang C, Zhan M, Dong Z. Mitochondrial function and disturbances in the septic kidney. Seminars in nephrology. 2015;35(1):108–19. Epub 2015/03/22. PubMed Central PMCID: PMCPMC4465453. pmid:25795504

16. Aziz M, Jacob A, Yang WL, Matsuda A, Wang P. Current trends in inflammatory and immunomodulatory mediators in sepsis. Journal of leukocyte biology. 2013;93(3):329–42. Epub 2012/11/09. PubMed Central PMCID: PMCPMC3579020. pmid:23136259

17. Quoilin C, Mouithys-Mickalad A, Lecart S, Fontaine-Aupart MP, Hoebeke M. Evidence of oxidative stress and mitochondrial respiratory chain dysfunction in an in vitro model of sepsis-induced kidney injury. Biochimica et biophysica acta. 2014;1837(10):1790–800. Epub 2014/07/16. pmid:25019585

18. Christian P, Su Q. MicroRNA regulation of mitochondrial and ER stress signaling pathways: implications for lipoprotein metabolism in metabolic syndrome. American journal of physiology Endocrinology and metabolism. 2014;307(9):E729–37. Epub 2014/09/04. pmid:25184990

19. Lorenzen JM, Kielstein JT, Hafer C, Gupta SK, Kumpers P, Faulhaber-Walter R, et al. Circulating miR-210 predicts survival in critically ill patients with acute kidney injury. Clinical journal of the American Society of Nephrology: CJASN. 2011;6(7):1540–6. Epub 2011/06/28. pmid:21700819

20. Saikumar J, Hoffmann D, Kim TM, Gonzalez VR, Zhang Q, Goering PL, et al. Expression, circulation, and excretion profile of microRNA-21, -155, and -18a following acute kidney injury. Toxicological sciences: an official journal of the Society of Toxicology. 2012;129(2):256–67. Epub 2012/06/19. PubMed Central PMCID: PMCPMC3499041.

21. Lan YF, Chen HH, Lai PF, Cheng CF, Huang YT, Lee YC, et al. MicroRNA-494 reduces ATF3 expression and promotes AKI. Journal of the American Society of Nephrology: JASN. 2012;23(12):2012–23. Epub 2012/11/20. PubMed Central PMCID: PMCPMC3507368. pmid:23160513

22. Xue Y, Wei Z, Ding H, Wang Q, Zhou Z, Zheng S, et al. MicroRNA-19b/221/222 induces endothelial cell dysfunction via suppression of PGC-1alpha in the progression of atherosclerosis. Atherosclerosis. 2015;241(2):671–81. Epub 2015/06/29. pmid:26117405

23. Haden DW, Suliman HB, Carraway MS, Welty-Wolf KE, Ali AS, Shitara H, et al. Mitochondrial biogenesis restores oxidative metabolism during Staphylococcus aureus sepsis. American journal of respiratory and critical care medicine. 2007;176(8):768–77. Epub 2007/06/30. PubMed Central PMCID: PMCPmc2020830. pmid:17600279

24. Sweeney TE, Suliman HB, Hollingsworth JW, Piantadosi CA. Differential regulation of the PGC family of genes in a mouse model of Staphylococcus aureus sepsis. PloS one. 2010;5(7):e11606. Epub 2010/07/27. PubMed Central PMCID: PMCPMC2905396. pmid:20657826

25. Fourcade S, Ferrer I, Pujol A. Oxidative stress, mitochondrial and proteostasis malfunction in adrenoleukodystrophy: A paradigm for axonal degeneration. Free radical biology & medicine. 2015;88(Pt A):18–29. Epub 2015/06/16.

26. Sun LY, Wang N, Ban T, Sun YH, Han Y, Sun LL, et al. MicroRNA-23a mediates mitochondrial compromise in estrogen deficiency-induced concentric remodeling via targeting PGC-1alpha. Journal of molecular and cellular cardiology. 2014;75:1–11. Epub 2014/07/02. pmid:24984145

27. Chaudhuri AD, Yelamanchili SV, Marcondes MC, Fox HS. Up-regulation of microRNA-142 in simian immunodeficiency virus encephalitis leads to repression of sirtuin1. FASEB journal: official publication of the Federation of American Societies for Experimental Biology. 2013;27(9):3720–9. Epub 2013/06/12. PubMed Central PMCID: PMCPMC3752547.

28. Rossaint J, Spelten O, Kassens N, Mueller H, Van Aken HK, Singbartl K, et al. Acute loss of renal function attenuates slow leukocyte rolling and transmigration by interfering with intracellular signaling. Kidney international. 2011;80(5):493–503. Epub 2011/05/13. PubMed Central PMCID: PMCPMC3156340. pmid:21562471

29. Holly MK, Dear JW, Hu X, Schechter AN, Gladwin MT, Hewitt SM, et al. Biomarker and drug-target discovery using proteomics in a new rat model of sepsis-induced acute renal failure. Kidney international. 2006;70(3):496–506. Epub 2006/06/09. PubMed Central PMCID: PMCPMC2276658. pmid:16760904

30. Wen X, Murugan R, Peng Z, Kellum JA. Pathophysiology of acute kidney injury: a new perspective. Contributions to nephrology. 2010;165:39–45. Epub 2010/04/30. pmid:20427954

31. Benes J, Chvojka J, Sykora R, Radej J, Krouzecky A, Novak I, et al. Searching for mechanisms that matter in early septic acute kidney injury: an experimental study. Critical care (London, England). 2011;15(5):R256. Epub 2011/10/28. PubMed Central PMCID: PMCPMC3334807.

32. Cho E, Lee JH, Lim HJ, Oh SW, Jo SK, Cho WY, et al. Soluble CD25 is increased in patients with sepsis-induced acute kidney injury. Nephrology (Carlton, Vic). 2014;19(6):318–24. Epub 2014/03/22.

33. Chawla LS, Seneff MG, Nelson DR, Williams M, Levy H, Kimmel PL, et al. Elevated plasma concentrations of IL-6 and elevated APACHE II score predict acute kidney injury in patients with severe sepsis. Clinical journal of the American Society of Nephrology: CJASN. 2007;2(1):22–30. Epub 2007/08/21. pmid:17699383

34. Murugan R, Karajala-Subramanyam V, Lee M, Yende S, Kong L, Carter M, et al. Acute kidney injury in non-severe pneumonia is associated with an increased immune response and lower survival. Kidney international. 2010;77(6):527–35. Epub 2009/12/25. PubMed Central PMCID: PMCPMC2871010. pmid:20032961