Transcriptional progressive patterns from mild to severe renal ischemia/reperfusion-induced kidney injury in mice

Abstract

The renal ischemia/reperfusion (I/R)-induced acute kidney injury incidence after nephron-sparing surgery for localized renal tumors is 20%, but the biological determinant process of postoperative acute kidney injury remains unclear. Using Gene Expression Omnibus database (GSE192883) and several bioinformatics analyses (discrete time points analysis, gene set enrichment analysis, dynamic network biomarker analysis, etc), combined with the establishment of the I/R model for verification, we identified three progressive patterns involving five core pathways confirmed using gene set enrichment analysis and six key genes (S100a10, Pcna, Abat, Kmo, Acadm, and Adhfe1) verified using quantitative polymerase chain reaction The dynamic network biomarker (DNB) subnetwork composite index value is the highest in the 22-min ischemia group, suggesting the transcriptome expression level fluctuated sharply in this group, which means 22-min ischemia is an critical warning point. This study illustrates the core molecular progressive patterns from mild to severe I/R kidney injury, laying the foundation for precautionary biomarkers and molecular intervention targets for exploration. In addition, the safe renal artery blocking time of nephron-sparing surgery that we currently accept may not be safe anymore.

Introduction

To date, nephron-sparing surgery (NSS) has been the standard treatment for localized renal tumors (Campbell et al., 2009), and offers equivalent oncologic outcomes (Scosyrev et al., 2014) with better preservation of renal function (Capitanio et al., 2015) compared with radical nephrectomy. Renal ischemia/reperfusion (I/R)-induced acute kidney injury (AKI) manifests as a rapid decline in renal function after NSS and is associated with high mortality and an increased risk of chronic kidney disease (Bedford et al., 2012; Birnie et al., 2014; Chawla et al., 2014).

Approximately 20% of patients develop AKI after NSS, which is mainly caused by an insufficient compensatory capacity of the contralateral kidney (Simmons et al., 2013; Bravi, 2020). However, the impact on renal function may be hidden by compensatory mechanisms in the contralateral kidney and evolve into unilateral kidney atrophy (Lane et al., 2011; Porpiglia et al., 2012). Therefore, it is important to study the biological progressive patterns from mild to severe I/R kidney injury to determine the occurrence of postoperative AKI.

Ischemic time has been associated with the degree of postoperative kidney injury in mice; a 22-min ischemia time leads to restorable kidney injury, and a 25/30-min ischemia time leads to non-restorable kidney injury (Wei and Dong, 2012; Liu et al., 2017). For humans, recommendations suggest that ischemia time should be controlled to 20–25 min (Thompson et al., 2010). Different ischemic times result in different degrees of kidney injury in different patients. Therefore, studying the molecular mechanism of renal I/R injury is the key to accurately evaluating the potential degree of kidney injury and to discovering potential molecular targets.

Previous studies have focused on the mechanisms of severe kidney injury and recovery after reperfusion, and current mainstream theories include cell death, inflammatory response, and fibrous repairing (Kumar et al., 2014; Venkatachalam et al., 2015; Cippa et al., 2018). For sequencing-related studies, Andrew et al. revealed molecular characterization of the transition from acute to chronic kidney injury following I/R by studying changes in molecular patterns, including cell death and proliferation (Krt8, Krt20, and Sox9), cell cycle and wound repair (Havcr1 and Lcn2), cell adhesion and inflammation (Timp2), adaptive immune responses (Ptprc, Cd3d, Cd74, and Cd48) and tubular function (Kap and Lrp2) (Liu et al., 2017). They clearly revealed the molecular progressive patterns from AKI to chronic kidney disease after renal I/R but did not focus on molecular changes in mild I/R kidney injury or the progressive patterns from mild to severe I/R kidney injury.

I/R injury is a cascade amplification process that involves progressive patterns from mild to severe I/R kidney injury. Unlike previous studies investigating a single molecule or several molecules, it seems that progressive patterns from mild to severe I/R kidney injury are more conducive for urologists to understand postoperative AKI and explore precautionary markers, even monitoring the progress of molecular biology in real time during NSS.

Herein, we used RNA sequencing to detect transcriptome expression in continuous ischemic time groups from mild to severe injury in a mouse I/R model as reported previously (Battistone et al., 2020). The observed transcriptional differences were reproduced using specific gene and protein expression analyses. We found gene sets with significant changes in mild and severe kidney injury and described and verified the progressive pattern from mild to severe injury through functional enrichment analysis, discrete time points analysis (Kumar and Matthisa, 2007), dynamic network biomarker (DNB) (Liu et al., 2014), and gene set enrichment analysis (GSEA) (Subramanian et al., 2005), laying the foundation for further research.

Materials and methods

Animals

All animal experiments were performed in accordance with international guidelines and approved by the Institutional Animal Care and Use Committee of the PLA General Hospital. Male C57BL/6 mice aged 6–8 weeks were purchased from GemPharmatech Co., Ltd. Mice were housed in a pathogen-free, constant temperature environment with a 12 h light/dark cycle and allowed to acclimatize for a week in the animal facility before the operation.

Ischemic renal injury model

Mice were subjected to bilateral renal I/R injury as previously described (Zhou et al., 2012). Briefly, mice were anesthetized with isoflurane on a thermostatic blanket while maintaining the temperature of the mice at approximately 37°C. Bilateral renal pedicles were dissociated through a middle abdominal incision, and the bilateral renal pedicles were clamped with atraumatic vascular clips for 18 or 30 min (n = 3). The incision was closed in two layers, and the mice were injected subcutaneously with 0.3 ml of warm saline on the back after surgery for volume supplementation. Sham-operated mice underwent the same procedure but without clamping of the renal pedicles.

Bioinformatics analysis

The edgeR package from Bioconductor was used to explore DEGs. DEGs were analyzed using the edgeR package (Robinson et al., 2010). DEGs were finally determined using parameters of false discovery rate adjusted p value <0.01, and log2 fold change >1.2 or < −1.2 unless specified. The Mfuzz (Kumar and M) package from Bioconductor was used for the discrete time points analysis. GO and KEGG analyses were performed to enrich the biological processes and molecular pathways of different gene sets (Ashburner et al., 2000; Yamada et al., 2006). GSEA was used to identify the activation or inhibition of the screened pathways (Subramanian et al., 2005). Enrichment p-values were adjusted using the Benjamini and Hochberg method, and p-adjusted values of less than 0.05 were determined to be significantly enriched.

DNB analysis was performed as previously reported (Ideker and Sharan, 2008; Jin et al., 2008; Liu et al., 2014). DNB is a dynamic analysis conducted through the variation of genes themselves and the correlation between genes, A gene group with larger variation and less correlation with other genes means that the gene group has outlier performance and becomes unstable, thus determining that the biological process is transforming to another stage. When another stage of the disease is reached, this phenomenon disappears. Briefly, we selected genes with significantly high deviations at each time point. Next, we built a distance matrix with Pearson correlation coefficients (PCC) for each pair of genes and set a threshold of 0.9 to group genes. The composite index (CI) was calculated from the average standard deviation (SD) of DNB molecules, the average absolute PCC value inside the DNB cluster (PCC₁), and the average absolute PCC value of the DNB cluster and others (PCC₀), as follows:

Quantitative RT-PCR assays

After 24 h of reperfusion, the kidneys of all mice were used for RNA extraction. TRIzol Reagent (Invitrogen, Carlsbad, CA, United States) was used in accordance with the manufacturer’s instructions, and reverse transcription (RT) was performed using the QuantiTect Reverse Transcription Kit. The cDNA samples were used for gene expression quantification via real-time PCR with a 2x Super SYBR Green qPCR Master Mix purchased from ESscience (ES-QP002, Shanghai, China), and RT-PCR was performed using an Applied Biosystems QuantStudio3 (Applied Biosystems). Each gene was run in three technical replicates and normalized to Ppia, and the fold-change relative to the sham group was calculated using the 2^−ΔΔCt formula. Primers for RT-PCR were designed using Primer-BLAST (NCBI).

Statistical analysis

Each analysis was performed in at least three independent experiments. All data are presented as the mean ± SEM (standard error of the mean). Differences between two groups were determined by two-tailed Student’s t-test using the Prism five software (GraphPad, La Jolla, CA, United States). Two-way analysis of variance (ANOVA) was used for analysis between the three groups. Statistical significance was set at p < 0.05.

Results

Differentially expressed genes and discrete time points analysis

To identify the differentially expressed genes at various times after renal I/R in mice, we performed RNA sequencing to explore the transcriptional progressive patterns from mild to severe kidney injury using different ischemia-time treatments. RNA sequencing data were deposited into the Gene Expression Omnibus database (GSE192883). Using differential gene analysis, we found that compared with the sham group, there were 1,156 differentially expressed genes (DEGs; 571 upregulated and 585 downregulated) in the 18-min ischemia group, 1,256 DEGs (635 upregulated and 621 downregulated) in the 22-min ischemia group, 1,325 DEGs (674 upregulated) in the 26-min ischemia group, and 651 DEGs (downregulated) and 2,119 DEGs (1,063 upregulated and 1,056 downregulated) in the 30-min ischemia group (Figure 1A). Among them, the top 10 DEGs with the greatest differences were more than 9-fold upregulated or downregulated (Table 1). Combined with principal component analysis (Figure 1B), it can be concluded that a large number of DEGs exist not only in severe renal I/R injury but also in mild renal I/R injury, indicating that many biological processes have been initiated in mild renal I/R injury and continue to change with the aggravation of the injury.

FIGURE 1

Among the top 50 DEGs in each group, the proportion of overlapping genes between the 18-min and 22-min ischemia groups was 67.1%; the proportion of overlapping genes among the 18-min, 22-min, and 26-min ischemia groups was 64.3%; and the proportion of overlapping genes among the four experimental groups was 42.9%. This suggests that DEGs gradually changed with the aggravation of injury but the most different genes changed unsustainably and discontinuously as time changes (Figure 1C). The heat map of overlapping DEGs among the four experimental groups shows that the expression levels of these DEGs changed regularly as the injury intensified (Figure 1D), supporting the results of the discrete time points analysis.

During the time-series analysis, all genes were automated and clustered into eight groups. Cluster1 indicated that the change in expression level was constantly accompanied by injury, which was a gene set related to I/R injury (Figure 1E). Cluster2 and Cluster3 showed gene sets that changed only in severe injury and represented the gene sets related to severe I/R injury (Figure 1F). Cluster4 to Cluster8 gene sets relate to the progression from mild to severe I/R kidney injury and showed a trend of continuous upregulation (Cluster4 and Cluster5), upregulation followed by downregulation (Cluster6 and Cluster7), and continuous downregulation (Cluster8). Moreover, the 18-min and 30-min ischemia groups were the two time points with more obvious changes, suggesting progressive patterns from mild to severe renal I/R injury (Figure 1G). Phenotypically, the injury initially occurred in the IR30 group but not the IR18, IR22 or IR26 group suggests that the transcriptome in the IR18 group has undergone changes, while they are not reflected in the pathological sections (Supplementary Figure S4). The results of the Cluster1–8 Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) functional enrichment analyses are presented in Table 2.

GSEA and qPCR verification

The core pathway of Cluster4 was identified as “Ribonucleoprotein complex biogenesis” through the overlapping pathways and pathway network analysis (Figure 2A) of the 18-min and 30-min ischemia groups; the GSEA result (Figure 2B) showed that the activation level of this pathway in the 30-min ischemia group was higher than that in the 18-min ischemia group. This suggests that the biogenesis function of the ribonucleoprotein complex was activated during the transition from mild to severe I/R kidney injury. The heat map of core enrichment genes, which contributed the most to the enrichment results, showed the same transitional pattern (Figure 2C). The core pathways of Cluster5 and Cluster8 were identified as “Cell-substrate junction organization” and “Organic acid catabolic process” through the same method as Cluster4 (Figures 2D,G); GSEA results (Figures 2E,H) showed that the activation level of “Cell-substrate junction organization” in the 30-min ischemia group was higher than that in the 18-min ischemia group and the inhibition level of “Organic acid catabolic process” in the 30-min ischemia group was higher than that in the 18-min ischemia group. This indicates that the adhesive function between the cell and substrate was activated and that the organic acid catabolic process was inhibited during the transition from mild to severe I/R kidney injury.

FIGURE 2

The heat map of the core enrichment genes showed the same transitional pattern (Figures 2F,I). The core pathways of Cluster6 and Cluster7 were identified as “Chromosome segregation” and “DNA repair” through the same method (Figures 3A,D), and the GSEA results (Figures 3B,E) showed that there were no significant differences between the activation level of these two pathways in the 18-min and the 30-min ischemia groups; however, in the leading edge subset, the pathway activation level of the 18-min group was significantly higher than that of the 30-min group. The heat map shows a similar transitional pattern (Figures 3C,F). The other representative pathways of the five clusters with significant differences, as well as the representative pathways of other clusters, are also shown (Figure 3G, Supplementary Figures S1, S2). Genes with a limit fold-change (LFC) greater than 0.7 in the IR18 group, LFC greater than one in the IR30 group or LFC greater than one in the IR18 group, and LFC greater than 0.7 in the IR30 group were screened out in the leading edge subset of each core pathway as key genes (Supplementary Figure S3), and quantitative polymerase chain reaction (qPCR) verification was conducted. The results showed that S100a10 was continuously upregulated in the IR18 and IR30 groups, which is consistent with the sequencing results. Pcna was upregulated in the IR18 group and downregulated in the IR30 group. Abat, Kmo, Acadm, and Adhfe1 were continuously downregulated (Figure 4B), which is consistent with the sequencing results for these molecules (Figure 4A). However, there were no significant changes between IR18, IR22 and IR26 groups.

FIGURE 3

FIGURE 4

DNB analysis

Specifically, we used a time-course (16–30 min) of mice for the I/R injury with a renal pedicle clip. Using our sequencing data, we obtained gene subnetworks (DNB) with fluctuations in different ischemic time groups according to DNB analysis methods. The CI value of the DNB subnetwork was the highest in the 22-min ischemia group, suggesting that the transcriptome expression level fluctuated sharply in this group and that 22 min of ischemia may be the key time point in the transition from mild to severe I/R kidney injury. When the ischemia time was prolonged and the injury was aggravated, the CI value of DNB began to decline, confirming that DNB only specifically recognizes the state of “pre-disease”, as previously reported (Liu et al., 2014) (Figure 5B). The molecular network diagram shows that except for the 22-min ischemia group, DNB in the other ischemic groups had strong connections with other genes, indicating that the internal network resonance of DNB in this group was strong, while its relationship with other genes was significantly weakened (Figure 5A). Detailed calculation results are shown in Table 3. However, the standard deviation and internal correlation of DNB in the 22-min ischemia group were not the highest, suggesting that the key contribution to the transition from mild to severe injury during renal I/R lies in the weakened relationship between DNB and other genes. Briefly, during the transition from mild to severe I/R kidney injury, renal transcriptome expression levels fluctuated strongly when the ischemia time reached 22 min, and the injury began to transition from mild to severe (Figure 5C). Therefore, 22 min of ischemia can be considered a warning point for severe irreversible I/R kidney injury.

FIGURE 5

Discussion

AKI can be caused by a variety of diseases, such as sepsis (Hoste et al., 2015), cardiorenal syndrome (McCullough et al., 2013), hepatorenal syndrome (Amin et al., 2019), and cardiac surgery-associated AKI(Hundemer et al., 2021). Different diseases have different characteristics in the process that leads to AKI. For urology, NSS-related I/R induced kidney injury is possessed of “controllability” and “monistic nature,” because the ischemia time in NSS is controllable, and the injury factor is relatively single compared with other diseases. Therefore, urologists should pay more attention to progressive patterns from mild to severe I/R kidney injury than to progressive patterns of kidney injury after I/R. This study focused on the progressive patterns from mild to severe I/R kidney injury after NSS, which provides a novel direction for monitoring and early warning of severe renal I/R injury, as well as the potential molecular targets for extending the time window of controllable kidney injury and striving for precious operation time for surgeons.

Previous studies on NSS-related AKI have suggested that ischemia time and renal tolerance to ischemia are decisive factors for the occurrence of postoperative AKI(Yossepowitch et al., 2006; Iida et al., 2008), but the internal molecular mechanism has not yet been explored. Most studies on molecular mechanisms have focused on biological changes in the kidney after reperfusion and the molecular or signaling pathways that play key roles in them. For example, KIM-1 is significantly upregulated in injured proximal renal tubular epithelial cells, and plays an important role in the clearance of necrotic tubular epithelial cell debris; it thus affects the repair process of renal injury (Yang et al., 2015). Many studies have been conducted on innate immunity-related molecules, such as IL-34, P2Y14, and PP2Acα, which play important roles in the process of renal injury and repair by regulating the function of macrophages and neutrophils (Baek et al., 2015; Battistone et al., 2020; Liang et al., 2021).

In contrast to the above studies, we propose a progressive pattern in which the core pathway best reflects the transition from mild to severe I/R kidney injury. Although these pathways may not play a decisive role in I/R injury and repair, they and the molecules involved may play a crucial role in determining the degree of I/R kidney injury. In other words, the core pathway involved in the pattern shown in our study is a key criterion for evaluating the degree of I/R kidney injury, and may also be an important target molecule that slows down the process from mild to severe I/R kidney injury. Pcna and Kmo have been reported to be closely associated with renal I/R injury and may be used as standards to evaluate the degree of injury and as targets for molecular intervention in the future (Li et al., 2017; Zheng X. et al., 2019).

However, the remaining four molecules have not yet been reported, among which S100a10 is associated with cell adhesion to the matrix and is also involved in the regulation of cell cycle and differentiation. Abat, Acadm, and Adhfe1 are metabolism-related molecules, among which Abat and Adhfe1 are involved in the metabolism of butyrate in cells. According to literature reports, butyrate metabolism is closely related to I/R kidney injury (Zheng Y. et al., 2019; Sun et al., 2022); therefore, Abat and Adhfe1 may also play an important role in I/R kidney injury, suggesting that they can be used as important indices to evaluate the potential degree of I/R kidney injury and important molecular targets for intervention. Meanwhile, for therapy, the downregulation of Abat and Adhfe1 leads to butyrate metabolic imbalances, which significantly enhanced renal dysfunction and histologic damage induced by renal IRI. Repairing metabolic imbalances through interventions can help increase butyrate which can cause a significant attenuation of neutrophil infiltration, which was reflected by the reduction of renal MPO activity, reduce apoptotic tubular cell death and improve caspase-3 activation. Therefore, Abat and Adhfe1 protein may be a potential therapeutic agent for preventing renal IRI (Zheng Y. et al., 2019).

Additionally, DNB analysis was introduced in this study to verify the existence of a “pre-disease” state through a new theory in renal I/R models, and it is proposed that before the occurrence of biological progressive patterns, the organism has already sent signals (DNB violent fluctuation). This provides a new idea for clinical warning, monitoring, and intervention for renal I/R injury. In addition, the safe renal artery blocking time of nephron-sparing surgery that we currently accept may not be safe anymore. However, our study also has some limitations. For example, the number of mice in each group was not sufficiently large, which may have caused bias in the analysis process. At the verification level, only RNA expression was verified, and further in-depth and comprehensive studies are required.

References

• 1

  AminA. A.AlabsawyE. I.JalanR.DavenportA. (2019). Epidemiology, pathophysiology, and management of hepatorenal syndrome. Semin. Nephrol.39 (1), 17–30. 10.1016/j.semnephrol.2018.10.002PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 2

  AshburnerM.BallC. A.BlakeJ. A.BotsteinD.ButlerH.CherryJ. M.et al (2000). Gene ontology: tool for the unification of biology. The gene Ontology consortium. Nat. Genet.25 (1), 25–29. 10.1038/75556PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 3

  BaekJ. H.ZengR.Weinmann-MenkeJ.ValeriusM. T.WadaY.AjayA. K.et al (2015). IL-34 mediates acute kidney injury and worsens subsequent chronic kidney disease. J. Clin. Invest.125 (8), 3198–3214. 10.1172/JCI81166PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 4

  BattistoneM. A.MendelsohnA. C.SpallanzaniR. G.AllegrettiA. S.LibermanR. N.SesmaJ.et al (2020). Proinflammatory P2Y14 receptor inhibition protects against ischemic acute kidney injury in mice. J. Clin. Invest.130 (7), 3734–3749. 10.1172/JCI134791PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 5

  BedfordM.FarmerC.LevinA.AliT.StevensP. (2012). Acute kidney injury and CKD: Chicken or egg?Am. J. Kidney Dis.59 (4), 485–491. 10.1053/j.ajkd.2011.09.010PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 6

  BirnieK.VerheydenV.PaganoD.BhabraM.TillingK.SterneJ. A.et al (2014). Predictive models for kidney disease: Improving global outcomes (KDIGO) defined acute kidney injury in UK cardiac surgery. Crit. Care18 (6), 606. 10.1186/s13054-014-0606-xPubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 7

  BraviC. A. (2020). Reply to Won Ho Kim, Hyun-Kyu Yoon, Chang Wook Jeong's Letter to Editor re: Carlo Bravi, Emily Vertosick, Nicole Benfante, et al. Impact of Acute Kidney Injury, Its Duration on Long-term Renal Function After Partial Nephrectomy. Eur. Urol.76, 398–403. Eur Urol 77(1), e16-e17. 10.1016/j.eururo.2019.09.022PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 8

  CampbellS. C.NovickA. C.BelldegrunA.BluteM. L.ChowG. K.DerweeshI. H.et al (2009). Guideline for management of the clinical T1 renal mass. J. Urol.182 (4), 1271–1279. 10.1016/j.juro.2009.07.004PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 9

  CapitanioU.TerroneC.AntonelliA.MinerviniA.VolpeA.FurlanM.et al (2015). Nephron-sparing techniques independently decrease the risk of cardiovascular events relative to radical nephrectomy in patients with a T1a-T1b renal mass and normal preoperative renal function. Eur. Urol.67 (4), 683–689. 10.1016/j.eururo.2014.09.027PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 10

  ChawlaL. S.EggersP. W.StarR. A.KimmelP. L. (2014). Acute kidney injury and chronic kidney disease as interconnected syndromes. N. Engl. J. Med.371 (1), 58–66. 10.1056/NEJMra1214243PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 11

  CippaP. E.SunB.LiuJ.ChenL.NaesensM.McMahonA. P.et al (2018). Transcriptional trajectories of human kidney injury progression. JCI Insight3 (22), 123151. 10.1172/jci.insight.123151PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 12

  HosteE. A.BagshawS. M.BellomoR.CelyC. M.ColmanR.CruzD. N.et al (2015). Epidemiology of acute kidney injury in critically ill patients: the multinational AKI-EPI study. Intensive Care Med.41 (8), 1411–1423. 10.1007/s00134-015-3934-7PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 13

  HundemerG. L.SrivastavaA.JacobK. A.KrishnasamudramN.AhmedS.BoergerE.et al (2021). Acute kidney injury in renal transplant recipients undergoing cardiac surgery. Nephrol. Dial. Transpl.36 (1), 185–196. 10.1093/ndt/gfaa063CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 14

  IdekerT.SharanR. (2008). Protein networks in disease. Genome Res.18 (4), 644–652. 10.1101/gr.071852.107PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 15

  IidaS.KondoT.AmanoH.NakazawaH.ItoF.HashimotoY.et al (2008). Minimal effect of cold ischemia time on progression to late-stage chronic kidney disease observed long term after partial nephrectomy. Urology72 (5), 1083–1088. ; discussion 1088-1089. 10.1016/j.urology.2008.06.074PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 16

  JinG.ZhouX.WangH.ZhaoH.CuiK.ZhangX. S.et al (2008). The knowledge-integrated network biomarkers discovery for major adverse cardiac events. J. Proteome Res.7 (9), 4013–4021. 10.1021/pr8002886PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 17

  KumarL.MatthisaE. F. (2007). Mfuzz: a software package for soft clustering of microarray data. Bioinformation2 (1), 5–7. 10.6026/97320630002005PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 18

  KumarS.LiuJ.McMahonA. P. (2014). Defining the acute kidney injury and repair transcriptome. Semin. Nephrol.34 (4), 404–417. 10.1016/j.semnephrol.2014.06.007PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 19

  LaneB. R.RussoP.UzzoR. G.HernandezA. V.BoorjianS. A.ThompsonR. H.et al (2011). Comparison of cold and warm ischemia during partial nephrectomy in 660 solitary kidneys reveals predominant role of nonmodifiable factors in determining ultimate renal function. J. Urol.185 (2), 421–427. 10.1016/j.juro.2010.09.131PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 20

  LiW. F.YangK.ZhuP.ZhaoH. Q.SongY. H.LiuK. C.et al (2017). Genistein ameliorates ischemia/reperfusion-induced renal injury in a SIRT1-dependent manner. Nutrients9 (4), E403. 10.3390/nu9040403PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 21

  LiangY.SunX.WangM.LuQ.GuM.ZhouL.et al (2021). PP2Acα promotes macrophage accumulation and activation to exacerbate tubular cell death and kidney fibrosis through activating Rap1 and TNFα production. Cell Death Differ.28 (9), 2728–2744. 10.1038/s41418-021-00780-5PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 22

  LiuJ.KumarS.DolzhenkoE.AlvaradoG. F.GuoJ.LuC.et al (2017). Molecular characterization of the transition from acute to chronic kidney injury following ischemia/reperfusion. JCI Insight2 (18), 94716. 10.1172/jci.insight.94716PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 23

  LiuR.WangX.AiharaK.ChenL. (2014). Early diagnosis of complex diseases by molecular biomarkers, network biomarkers, and dynamical network biomarkers. Med. Res. Rev.34 (3), 455–478. 10.1002/med.21293PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 24

  McCulloughP. A.KellumJ. A.HaaseM.MullerC.DammanK.MurrayP. T.et al (2013). Pathophysiology of the cardiorenal syndromes: Executive summary from the eleventh consensus conference of the acute dialysis quality initiative (ADQI). Contrib. Nephrol.182, 82–98. 10.1159/000349966PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 25

  PorpigliaF.FioriC.BertoloR.MorraI.RussoR.PiccoliG.et al (2012). Long-term functional evaluation of the treated kidney in a prospective series of patients who underwent laparoscopic partial nephrectomy for small renal tumors. Eur. Urol.62 (1), 130–135. 10.1016/j.eururo.2012.02.001PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 26

  RobinsonM. D.McCarthyD. J.SmythG. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics26 (1), 139–140. 10.1093/bioinformatics/btp616PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 27

  ScosyrevE.MessingE. M.SylvesterR.CampbellS.Van PoppelH. (2014). Renal function after nephron-sparing surgery versus radical nephrectomy: results from EORTC randomized trial 30904. Eur. Urol.65 (2), 372–377. 10.1016/j.eururo.2013.06.044PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 28

  SimmonsM. N.LieserG. C.FerganyA. F.KaoukJ.CampbellS. C. (2013). Association between warm ischemia time and renal parenchymal atrophy after partial nephrectomy. J. Urol.189 (5), 1638–1642. 10.1016/j.juro.2012.11.042PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 29

  SubramanianA.TamayoP.MoothaV. K.MukherjeeS.EbertB. L.GilletteM. A.et al (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A.102 (43), 15545–15550. 10.1073/pnas.0506580102PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 30

  SunY.ZhouC.ChenY.HeX.GaoF.XueD.et al (2022). Quantitative increase in short-chain fatty acids, especially butyrate protects kidney from ischemia/reperfusion injury. J. Investig. Med.70 (1), 29–35. 10.1136/jim-2020-001715PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 31

  ThompsonR. H.LaneB. R.LohseC. M.LeibovichB. C.FerganyA.FrankI.et al (2010). Every minute counts when the renal hilum is clamped during partial nephrectomy. Eur. Urol.58 (3), 340–345. 10.1016/j.eururo.2010.05.047PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 32

  VenkatachalamM. A.WeinbergJ. M.KrizW.BidaniA. K. (2015). Failed tubule recovery, AKI-CKD transition, and kidney disease progression. J. Am. Soc. Nephrol.26 (8), 1765–1776. 10.1681/ASN.2015010006PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 33

  WeiQ.DongZ. (2012). Mouse model of ischemic acute kidney injury: technical notes and tricks. Am. J. Physiol. Ren. Physiol.303 (11), F1487–F1494. 10.1152/ajprenal.00352.2012PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 34

  YamadaT.KanehisaM.GotoS. (2006). Extraction of phylogenetic network modules from the metabolic network. BMC Bioinforma.7, 130. 10.1186/1471-2105-7-130PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 35

  YangL.BrooksC. R.XiaoS.SabbisettiV.YeungM. Y.HsiaoL. L.et al (2015). KIM-1-mediated phagocytosis reduces acute injury to the kidney. J. Clin. Invest.125 (4), 1620–1636. 10.1172/JCI75417PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 36

  YossepowitchO.EggenerS. E.SerioA.HuangW. C.SnyderM. E.VickersA. J.et al (2006). Temporary renal ischemia during nephron sparing surgery is associated with short-term but not long-term impairment in renal function. J. Urol.176 (4 Pt 1), 1339–1343. ; discussion 1343. 10.1016/j.juro.2006.06.046PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 37

  ZhengX.ZhangA.BinnieM.McGuireK.WebsterS. P.HughesJ.et al (2019a). Kynurenine 3-monooxygenase is a critical regulator of renal ischemia-reperfusion injury. Exp. Mol. Med.51 (2), 1–14. 10.1038/s12276-019-0210-xPubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 38

  ZhengY.ZhangZ.ZhangN. (2019b). Protective effects of butyrate on renal ischemia-reperfusion injury in rats. Urol. Int.102 (3), 348–355. 10.1159/000497476PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 39

  ZhouD.LiY.LinL.ZhouL.IgarashiP.LiuY.et al (2012). Tubule-specific ablation of endogenous beta-catenin aggravates acute kidney injury in mice. Kidney Int.82 (5), 537–547. 10.1038/ki.2012.173PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar