The analysis of a time-course transcriptome profile by systems biology approaches reveals key molecular processes in acute kidney injury
Kobra Moradzadeh1, Yousof Gheisari2
1 Department of Genetics and Molecular Biology, Isfahan University of Medical Sciences, Isfahan, Iran
2 Department of Genetics and Molecular Biology, Isfahan University of Medical Sciences; Regenerative Medicine Research Center, Isfahan University of Medical Sciences, Isfahan, Iran
|Date of Submission||27-Aug-2018|
|Date of Decision||14-Sep-2018|
|Date of Acceptance||07-Oct-2018|
|Date of Web Publication||31-Jan-2019|
Dr. Yousof Gheisari
Department of Genetics and Molecular Biology, Isfahan University of Medical Sciences, Isfahan 8174673461
Source of Support: None, Conflict of Interest: None
Background: Acute kidney injury is a common debilitating disease with no curative treatment. The recent development of big biological data is expected to expand our understanding of the disorder if appropriately analyzed to generate translational knowledge. We have here re-analyzed a time-course microarray data on mRNA expression of rat kidneys exposed to ischemia-reperfusion to identify key underlying biological processes. Materials and Methods: The dataset was quality controlled by principal component analysis and hierarchical clustering. Using limma R package, differentially expressed (DE) genes were detected which were then clustered according to their expression trajectories. The biological processes related to each cluster were harvested using gene ontology enrichment analysis. In addition, the interaction map of proteins encoded by the DE genes was constructed, and the functions related to network central genes were determined. Furthermore, signaling pathways related to the DE genes were harvested using pathway enrichment analysis. Results: We found 8139 DE genes that drive critical processes such as the control of blood circulation, reactive species metabolism, mitochondrial respiration, apoptosis, cell proliferation, as well as inflammatory and immunological reactions. The role of less recognized pathways such as olfactory signaling in acute kidney injury is also proposed that remains to be investigated in future studies. Conclusion: Using systems biology top-down approach, we have suggested novel potential genes and pathways to be intervened toward kidney regeneration.
Keywords: Acute kidney injury, reperfusion injury, systems biology
|How to cite this article:|
Moradzadeh K, Gheisari Y. The analysis of a time-course transcriptome profile by systems biology approaches reveals key molecular processes in acute kidney injury. J Res Med Sci 2019;24:3
|How to cite this URL:|
Moradzadeh K, Gheisari Y. The analysis of a time-course transcriptome profile by systems biology approaches reveals key molecular processes in acute kidney injury. J Res Med Sci [serial online] 2019 [cited 2019 Apr 21];24:3. Available from: http://www.jmsjournal.net/text.asp?2019/24/1/3/251201
| Introduction|| |
Despite huge investigations on acute kidney injury, it is yet a common cause of death in hospital-admitted patients,,, and the existing therapeutic options are just supportive. The insufficiency of current therapies for this clinical condition is partly due to our limited understanding of the molecular mechanisms. The recent availability of omics data in line with the emerging hypothesis-free paradigm of systems biology paves the way for more comprehensive insights toward the molecular pathophysiology of the disorder. This holistic approach is expected to accelerate the development of novel efficient therapeutics.
Although a huge amount of biological big data has been generated in the last few years, time-course expression profiling studies encompass a minority of the datasets. Considering the dynamic nature of biological processes, time-series evaluations are more appreciated as they represent the course of events more precisely. Here, we selected a dataset deposited by Speir et al. in which the mRNA expression profile of rat kidneys exposed to ischemia-reperfusion (IR) injury is assessed in several time steps following the insult. Following quality control steps, we have re-analyze this dataset to identify key role player genes as well as critical signaling pathways and molecular processes related to acute kidney injury.
| Materials and Methods|| |
mRNA expression profiling deposited by Speir et al. with the accession number GSE58438 was retrieved from the Gene Expression Omnibus database., The series matrix file of the data has been normalized by the method of Sketch-Quantile. Dataset quality was assessed by principal component analysis using ggplot2 package and hierarchical clustering using pheatmap package of R software which are accessible through CRAN repository. The clustering method was “complete” with “Euclidean” distance measure. Differentially expressed (DE) genes were identified using limma package of R software, available in Bioconductor repository. In the analysis, IR rats harvested at 3, 24, and 120 h postischemia as well as normal controls were compared by the method of F-statistic. Benjamini–Hochberg false discovery rate (FDR) was then used for P value correction. The genes with adjusted (adj.) P < 0.05 were considered as DE. Venn diagram of DE genes was depicted using Eulerr package of CRAN repository.
For clustering of the DE genes, Short Time-Series Expression Miner version 1.3.11 was used. For this method, the matrix files of DE genes containing the expression values of three samples in each group were provided. FDR was used as the correction method in model profile section.
Gene ontology enrichment analysis
The ClueGO plugin version 2.3.5 of Cytoscape version 3.2.1 was applied to retrieve the gene ontology (GO) data from GO database. Bonferroni step-downmethod was chosen for P value correction and adj. P < 0.05 was considered significant. REViGO was used to find the parent terms of the retrieved GO terms and the similarity index was set at 0.4.
Protein-protein interaction network
Cytoscape CluePedia plugin version 1.3.5 was used for the construction of protein-protein interaction (PPI) network. The activation, inhibition, binding, and posttranslational modification interactions were retrieved from STRING database, and the confidence cutoff was set at 0.8. The network topology was analyzed using network analyzer tool.
Pathway enrichment analysis
Cytoscape ClueGO plugin was employed for pathway enrichment analysis. KEGG, and REACTOME databases were used to find related pathways. Bonferroni step-downmethod was applied for P value correction, and pathways with adj. P < 0.05 were selected. Ram and Wickham package was applied for drawing the scatter plot.
| Results|| |
In this study, we have re-analyzed the microarray dataset “GSE58438” which includes the mRNA profiles of rat kidneys at 3, 24, and 120 h following IR injury as well as normal controls. To assess the quality of the data, principal component analysis and hierarchical clustering were employed. Both these unsupervised methods revealed the segregation of the samples in four subsets matching with the time steps of the study [Figure 1]a. However, three samples that were not scattered according to the experiment group were omitted in downstream analyses to improve the quality of the dataset.
|Figure 1: The quality of the microarray dataset was satisfactory as most samples were segregated according to the experimental groups in unsupervised principal component analysis and hierarchical clustering. Three samples that did not follow this trend were omitted in downstream analyses (a). The numbers of DE genes in each IR time point are depicted in the Venn diagram (b)|
Click here to view
To identify DE genes, limma package of R software was employed which revealed 8139 genes to be variably expressed with FDR <0.05 [Supplementary Table 1 [Additional file 1]]; [Figure 1]b. As the animal model used in this study is generated with a transient episode of ischemia, we were interested to focus on the profiles of DE genes related to this condition and plotted their expression trajectories over time [Figure 2]; Actb, the activator of Siah2, increases following IR which is associated with a transient Siah2 overexpression. In addition, Hif1an (Fih1) declines and Ard1a increases shortly after injury. These alterations can lead to Egln2 (Phd1) inactivation and Hif1a stabilization. Therefore, the expression of these genes is according to their role for protecting the cells against hypoxia. In addition, Cxcl12 (Sdf1) shows a declining trend which is in line with our previous data indicating the downregulation of this gene both at mRNA and protein levels following kidney IR. Furthermore, Tgfb is observed to be overexpressed which is in agreement with previous studies [Figure 2].,
|Figure 2: The temporal expression profiles of some differentially expressed genes related to tissue hypoxia are shown|
Click here to view
Next, to provide a holistic view, the DE genes were clustered based on their temporal expression patterns, and the functions attributable to each group were identified by GO term enrichment analysis. Interestingly, in each cluster, one or two processes were highlighted, indicating that the genes with similar trajectories drive similar functions. The identified processes were in accordance with the molecular pathophysiology of the disorder and included reactive species metabolism, lipid peroxidation, mitochondrial metabolism, cell proliferation, regeneration, cilium organization, and immunological reactions [Figure 3].
|Figure 3: The DE genes were clustered based on their expression trajectories and gene ontology enrichment analysis was performed for each cluster. Parent gene ontology terms are demonstrated|
Click here to view
To assess the functional relationships between the proteins encoded by the DE genes, a PPI map was constructed, and its topology was analyzed. As expected for biological networks, the node degree in the constructed network followed a power low distribution with the form y = 1548.5x− 1.516, where y is the number of nodes and x is degree (correlation: 0.965). We have previously shown that central nodes in PPI are critical for cell functions., Therefore, based on degree, betweenness, and closeness centrality parameters, top 1% DE genes were selected [Figure 4]a and their ontologies were assessed [Figure 4]b. These central genes were mainly related to regulation of blood circulation, angiotensin signaling, reactive oxygen species metabolic processes, mitochondria organization, apoptosis, and cell proliferation. These processes are known to be of special importance in the pathogenesis of acute kidney injury.,,
|Figure 4: Based on the degree, betweenness centrality, and closeness centrality, top 1% genes in the PPI network were determined (a). Gene ontology enrichment analysis identified key biological processes, cellular components, and molecular functions related to these central genes (b)|
Click here to view
To further provide clues regarding to the biological functions in kidney ischemia, pathway enrichment analysis was performed with the DE genes and the signaling pathways with FDR <0.05 were harvested [Figure 5]. In agreement with previous studies,,,,, this analysis highlights the importance of processes such as mitochondrial respiration, apoptosis, cell cycle, DNA repair, and nuclear factor-kappa B activation in acute kidney injury. In addition, some pathways such as olfactory signaling are enriched that can be proposed as novel role players in this disorder. Although olfactory receptor family encompasses a large amount of genes, the employment of P value correction in the enrichment analyses precludes the chance that the enrichment of this pathway is a false-positive record.
|Figure 5: Pathway enrichment analysis was performed for all differentially expressed genes and signaling pathways with false discovery rate <0.05 are demonstrated|
Click here to view
| Discussion|| |
The 20th century was a time for great achievements in medicine such as the control of communicable diseases and the development of surgical procedures and diagnostic techniques. However, the treatment options for most complex disorders are yet supportive rather than curative. The recent emergence of systems biology has raised hopes to address this limitation through more comprehensive investigations that may lead to novel therapeutics. Although a huge amount of big data has been generated in the last years, they are not commonly appropriately analyzed to generate translational knowledge. Indeed, the speed of data generation has surpassed that of data analysis and so it has been advised to focus on the interpretation of available datasets rather than the generation of new ones., In this study, we have re-analyzed a time-course microarray mRNA profiling of acute kidney injury and generated a comprehensive outlook of molecular interactions, key processes, and signaling pathways to address the scarcity of holistic investigations of this disorder.
Although the current dataset is originally developed as a part of a study on the effect of valproic acid on kidney ischemia, we were here interested on the natural course of acute kidney ischemia and focused on the samples derived from rats exposed to IR but not receiving the drug. We have shown that the selected samples are acceptable for analysis according to unsupervised quality control measures. Another advantage is that the samples are harvested in different time points following the insult and so can represent the dynamic nature of the biological events. Indeed, single time-point expression profiling can be misleading due to complex trajectories and noise in genes expression. However, time-course microarray experiments are not common due to practical limitations. In addition, the analysis of such data is challenging and requires specific tools. The use of limma in the current work is based on a previous study indicating that this R package has an acceptable performance for the identification of DE genes in time-course microarray datasets and outperforms other similar algorithms.
Following the detection of DE genes by limma, we focused on the expression patterns of special genes related to tissue hypoxia among which Cxcl12 is interesting. This chemokine is positively regulated by Hif1 and is responsible to attract stem cells and leukocytes to injured tissues., Therefore, it is expected to be upregulated in response to ischemia, but it showed a declining trend in this dataset. We previously had a similar observation on Cxcl12 downregulation following kidney ischemia and using a stochastic Petri net mathematical model, proposing that this unexpected observation is due to the inhibitory effect of histone deacetylase on this gene which dominates Hif1-induced overexpression. This prediction is recently experimentally corroborated (unpublished data).
To explore the underlying mechanisms of this disorder, the DE genes were clustered based on their temporal expression patterns and the functions attributable to each cluster were detected. Furthermore, the interaction map of proteins encoded by the DE genes was constructed, and the ontology of the central network genes was assessed. Furthermore, enrichment analysis was performed to identify the signaling pathways in which the DE genes are involved. These investigations highlight the importance of essential processes such as free radical metabolism, cellular respiration, immunological responses, cell proliferation, and apoptosis. In addition, the involvement of olfactory receptors in acute kidney injury is here proposed which to the best of our knowledge has not been previously indicated. Interestingly, olfactory receptors are shown to be involved in normal renal function,, as well as polycystic kidney disease and nondiabetic chronic kidney disease. Future experimental studies are required to investigate if these receptors also contribute in the pathogenesis of acute kidney injury.
Financial support and sponsorship
This study was supported by Isfahan University of Medical Sciences (grant number 395115) and Iran National Science Foundation (grant number 9600608).
Conflicts of interest
There are no conflicts of interest.
| References|| |
Bellomo R, Kellum JA, Ronco C. Acute kidney injury. Lancet 2012;380:756-66.
Lameire NH, Bagga A, Cruz D, De Maeseneer J, Endre Z, Kellum JA, et al.
Acute kidney injury: An increasing global concern. Lancet 2013;382:170-9.
Li PK, Burdmann EA, Mehta RL; World Kidney Day Steering Committee 2013. Acute kidney injury: Global health alert. Kidney Int 2013;83:372-6.
Heydari B, Khalili H, Beigmohammadi MT, Abdollahi A, Karimzadeh I. Effects of atorvastatin on biomarkers of acute kidney injury in amikacin recipients: A pilot, randomized, placebo-controlled, clinical trial. J Res Med Sci 2017;22:39.
] [Full text]
Saberi K, Salehi M, Rahmanian M, Bakhshandeh AR, Massoumi GR. Appropriate blood component therapy can reduce postcardiac surgery acute kidney injury through packed cell transfusion reduction. J Res Med Sci 2017;22:80.
] [Full text]
He JC, Chuang PY, Ma'ayan A, Iyengar R. Systems biology of kidney diseases. Kidney Int 2012;81:22-39.
Speir RW, Stallings JD, Andrews JM, Gelnett MS, Brand TC, Salgar SK, et al.
Effects of valproic acid and dexamethasone administration on early bio-markers and gene expression profile in acute kidney ischemia-reperfusion injury in the rat. PLoS One 2015;10:e0126622.
Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res 2002;30:207-10.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al.
NCBI GEO: Archive for functional genomics data sets – Update. Nucleic Acids Res 2013;41:D991-5.
Wickham H. Ggplot2: Elegant Graphics for Data Analysis. New York: Springer-Verlag; 2009.
Team RC. R: A Language and Environment for Statistical Computing. Foundation for Statistical Computing, Vienna, Austria; 2016. Available from: https://cran.r-project.org/
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al.
Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 2015;43:e47.
Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al.
ClueGO: A cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 2009;25:1091-3.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al.
Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res 2003;13:2498-504.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al.
Gene ontology: Tool for the unification of biology. The gene ontology consortium. Nat Genet 2000;25:25-9.
Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One 2011;6:e21800.
Bindea G, Galon J, Mlecnik B. CluePedia cytoscape plugin: Pathway insights using integrated experimental and in silico
data. Bioinformatics 2013;29:661-3.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al.
The STRING database in 2017: Quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res 2017;45:D362-8.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res 2000;28:27-30.
Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res 2016;44:D457-62.
Fabregat A, Jupe S, Matthews L, Sidiropoulos K, Gillespie M, Garapati P, et al.
The reactome pathway knowledgebase. Nucleic Acids Res 2018;46:D649-55.
Heidary Z, Ghaisari J, Moein S, Naderi M, Gheisari Y. Stochastic petri net modeling of hypoxia pathway predicts a novel incoherent feed-forward loop controlling SDF-1 expression in acute kidney injury. IEEE Trans Nanobioscience 2016;15:19-26.
Gentle ME, Shi S, Daehn I, Zhang T, Qi H, Yu L, et al.
Epithelial cell TGFβ signaling induces acute tubular injury and interstitial inflammation. J Am Soc Nephrol 2013;24:787-99.
Spurgeon KR, Donohoe DL, Basile DP. Transforming growth factor-beta in acute renal failure: Receptor expression, effects on proliferation, cellularity, and vascularization after recovery from injury. Am J Physiol Renal Physiol 2005;288:F568-77.
Barabási AL, Oltvai ZN. Network biology: Understanding the cell's functional organization. Nat Rev Genet 2004;5:101-13.
Rabieian R, Abedi M, Gheisari Y. Central nodes in protein interaction networks drive critical functions in transforming growth factor beta-1 stimulated kidney cells. Cell J 2017;18:514-31.
Abedi M, Gheisari Y. Nodes with high centrality in protein interaction networks are responsible for driving signaling pathways in diabetic nephropathy. PeerJ 2015;3:e1284.
Makris K, Spanou L. Acute kidney injury: Definition, pathophysiology and clinical phenotypes. Clin Biochem Rev 2016;37:85-98.
Chatauret N, Badet L, Barrou B, Hauet T. Ischemia-reperfusion: From cell biology to acute kidney injury. Prog Urol 2014;24 Suppl 1:S4-12.
Kalogeris T, Baines CP, Krenz M, Korthuis RJ. Cell biology of ischemia/reperfusion injury. Int Rev Cell Mol Biol 2012;298:229-317.
Devarajan P. Update on mechanisms of ischemic acute kidney injury. J Am Soc Nephrol 2006;17:1503-20.
El Sabbahy M, Vaidya VS. Ischemic kidney injury and mechanisms of tissue repair. Wiley Interdiscip Rev Syst Biol Med 2011;3:606-18.
Zuk A, Bonventre JV. Acute kidney injury. Annu Rev Med 2016;67:293-307.
Bonavia A, Singbartl K. A review of the role of immune cells in acute kidney injury. Pediatr Nephrol 2018;33:1629-39.
Rung J, Brazma A. Reuse of public genome-wide gene expression data. Nat Rev Genet 2013;14:89-99.
Martens L, Vizcaíno JA. A golden age for working with public proteomics data. Trends Biochem Sci 2017;42:333-41.
Ladbury JE, Arold ST. Noise in cellular signaling pathways: Causes and effects. Trends Biochem Sci 2012;37:173-8.
Bar-Joseph Z, Gitter A, Simon I. Studying and modelling dynamic biological processes using time-series gene expression data. Nat Rev Genet 2012;13:552-64.
Moradzadeh K, Moein S, Nickaeen N, Gheisari Y. Analysis of time-course microarray data: Comparison of common tools. Genomics 2018. doi: 10.1016/j.ygeno.2018.03.021.
Cencioni C, Capogrossi MC, Napolitano M. The SDF-1/CXCR4 axis in stem cell preconditioning. Cardiovasc Res 2012;94:400-7.
Yellowley C. CXCL12/CXCR4 signaling and other recruitment and homing pathways in fracture repair. Bonekey Rep 2013;2:300.
Togel FE, Westenfelder C. Role of SDF-1 as a regulatory chemokine in renal regeneration after acute kidney injury. Kidney Int Suppl (2011) 2011;1:87-9.
Shepard BD, Pluznick JL. How does your kidney smell? Emerging roles for olfactory receptors in renal function. Pediatr Nephrol 2016;31:715-23.
Shepard BD, Cheval L, Peterlin Z, Firestein S, Koepsell H, Doucet A, et al.
Arenal olfactory receptor aids in kidney glucose handling. Sci Rep 2016;6:35215.
Pluznick JL, Rodriguez-Gil DJ, Hull M, Mistry K, Gattone V, Johnson CA, et al.
Renal cystic disease proteins play critical roles in the organization of the olfactory epithelium. PLoS One 2011;6:e19694.
Koseoglu S, Derin S, Huddam B, Sahan M. The effect of non-diabetic chronic renal failure on olfactory function. Eur Ann Otorhinolaryngol Head Neck Dis 2017;134:161-4.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5]