Bioinformatics analysis of pathways of renal infiltrating macrophages in different renal disease models.

BACKGROUND: Recent studies have suggested that macrophages are significantly involved in different renal diseases. However, the role of these renal infiltrating macrophages has not been entirely uncovered. To further clarify the underlying mechanism and identify therapeutic targets, a bioinformatic analysis based on transcriptome profiles was performed.
METHODS: Three transcription profiling datasets, GSE27045, GSE51466 and GSE75808, were obtained from the Gene Expression Omnibus (GEO). Differentially expressed genes (DEGs) were assessed by Gene Ontology (GO) functional annotation, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and gene set enrichment analysis (GSEA).
RESULTS: The classic signaling pathways and metabolic pathways of macrophages infiltrating the kidney in different pathophysiological processes, including lupus nephritis (LN), renal crystal formation and renal ischemia-reperfusion injury (IRI), were analysed. Furthermore, the common classical pathways significantly altered in the three renal disorders were the oxidative phosphorylation, VEGF signaling and JAK/STAT signaling pathways, while the renin-angiotensin system was uniquely altered in LN, the glycolysis and gluconeogenesis pathways were uniquely altered in models of renal crystal formation, and the calcium signaling pathway was specific to renal IRI.
CONCLUSIONS: Via bioinformatics analysis, this study revealed the transcriptional features of macrophages in murine LN, renal crystal formation and IRI models, which may serve as promising targets for mechanistic research and the clinical treatment of multiple renal diseases. 2021 Translational Andrology and Urology. All rights reserved.

Introduction

Increasing evidence has shown that macrophages participate in kidney injury, repair and fibrosis in multiple experimental and human renal diseases (1). In previous research, macrophages have been categorized according to different functions: classically activated macrophages (proinflammatory M1 macrophages) and alternatively activated macrophages (profibrotic/wound healing M2 macrophages) (2). Polarization of macrophages into the M1 or M2 state can be observed in vivo. For example, Lee et al. (3) showed phenotypic switching of macrophages from the M1 to the M2 state in the repair phase of acute kidney injury (AKI), and transferring M2 macrophages induces the repair process (4). Additionally, macrophage polarization has been observed in a rhabdomyolysis-induced AKI mouse model (5). However, recent studies suggest that there is a spectrum of macrophage activation states (6). Clements et al. (7) revealed that the expression of M2-specific markers in renal CD11b+/Ly6Cintermediate and CD11b+/Ly6Clow macrophage populations did not enable clear clustering of traditional M1 and M2 macrophages by ex vivo standards. Therefore, the characteristics of macrophage activation states in renal tissue should be defined based on other factors (8), such as the signaling pathways activated in the macrophages.

Systemic lupus erythematosus (SLE) is a systemic autoimmune disorder characterized by autoantibodies against self-nucleic acids and leads to inflammation and injury (9). Lupus nephritis (LN) is a major type of end-stage organ damage in SLE and is a main risk factor for mortality (10). Recent studies have reported that infiltration of macrophages is related to poor prognosis in LN (11). During active LN in nephritic NZB/W mice, renal macrophages located throughout the interstitium secrete proinflammatory cytokines and chemokines, which can also contribute to tissue repair and degradation (12). Interestingly, the phenotype of these renal macrophages is different from that of classic inflammatory macrophages (12). Therefore, immune-regulating therapies are promising for LN treatment. A further understanding of the renal macrophages that modulate LN tissue damage is needed.

The high prevalence and recurrence rate of kidney stones have been recognized as problems worldwide, and kidney stones may lead to chronic kidney disease and end-stage renal disease (13). The pathology of kidney stone formation is diverse (14). Inflammatory processes, oxidative stress, mineral and lipid metabolism, and insulin resistance are involved in calcium oxalate (CaOx) crystal development (15). Taguchi et al. (16) found that macrophages suppress renal crystal formation by using hyperoxaluric CSF-1-deficient mice that lack M2 macrophages. The researchers then sorted renal macrophages by flow cytometry and revealed that macrophages affect renal crystal formation during hyperoxaluria by gene array profiling. Understanding the role of macrophages in renal crystal formation can aid the treatment of renal crystal deposition.

Renal ischemia-reperfusion injury (IRI) is a leading cause of AKI and is associated with poor prognosis. A recent report revealed that macrophages play a vital role in the renal injury, repair, and fibrosis associated with AKI, but the molecular mechanism remains unclear (17). Clements et al. (7) utilized a renal bilateral IRI murine model to investigate the function of macrophages in the pathological process by sorting macrophages from the kidney for microarray analysis. Sun et al. (18) also conducted an analysis of mRNA expression profiles to define differentially expressed genes (DEGs), biological processes and pathways associated with IRI. Since understanding how macrophages are altered may provide insight into promising treatments for AKI that modulate the immune response, we aimed to characterize the molecular function of macrophages after renal IRI.

The aims of this study were to explore the role of infiltrating macrophages in different kidney diseases by analysing signaling pathways, including metabolism-related pathways, assess the similarities and differences in macrophage function between different renal diseases, and provide therapeutic targets.

We present the following article in accordance with the STROBE reporting checklist (available at https://dx.doi.org/10.21037/tau-21-761).

Methods

Data preparation

To reveal the role of macrophages in different renal diseases, we employed expression profiling array datasets of renal macrophages. After searching the Gene Expression Omnibus (GEO) database using the keywords ‘renal’ and ‘macrophage’, we found three GEO datasets (GSE27045, GSE51466 and GSE75808) that met our criteria. The transcription profiles were obtained from the GEO (http://www.ncbi.nlm.nih.gov/geo), a public functional genomics data repository. According to the annotation information in GEO, the probes were converted into the corresponding gene symbol.

In the GSE27045 dataset, F4/80hi macrophages were isolated from NZB/W mice in the early stage of lupus, during lupus nephritis and after induction of remission using flow cytometry. RNA was extracted and processed for hybridization on the Affymetrix GeneChip Mouse Genome 430 2.0 Array. This dataset contained six “young mac” samples (GSM667532, GSM667533, GSM667534, GSM667535, GSM667536 and GSM667537), seven “sick mac” samples (GSM667538, GSM667539, GSM667540, GSM667541, GSM667542, GSM667543, and GSM667544) and four “rem mac” samples (GSM667545, GSM667546, GSM667547, and GSM667548).

In the GSE51466 dataset, renal CD11b+/CD11c+ macrophages from wild-type stone model mice and CSF-1-deficient (op/op) stone model mice were sorted and investigated with the Agilent-028005 SurePrint G3 Mouse GE 8x60K Microarray. GSE51466 contained three “op/op” samples (GSM1246087, GSM1246088 and GSM1246089) and three “wild type” samples (GSM1246084, GSM1246085, and GSM1246086).

In the GSE75808 dataset, renal CD11b+/Ly6C+ monocyte/macrophage populations were sorted by flow cytometry at 4 h, 1 day and 9 days after AKI. Samples were processed for the Affymetrix Mouse Genome 430 2.0 Array. GSE75808 contained nine “sham” samples (GSM1968231, GSM1968232, GSM1968233, GSM1968234, GSM1968235, GSM1968236, GSM1968237, GSM1968238 and GSM1968239), six “4 h IRI” samples (GSM1968225, GSM1968226, GSM1968227, GSM1968228, GSM1968229 and GSM1968230), three “24 h IRI” samples (GSM1968216, GSM1968217 and GSM1968218) and six “9-day IRI” samples (GSM1968219, GSM1968220, GSM1968221, GSM1968222, GSM1968223 and GSM1968224).

The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

Identification of DEGs

DEGs were identified by using GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r). GEO2R is an interactive web tool that helps users identify DEGs across experimental conditions by comparing two or more datasets in a GEO series. Probe sets without corresponding gene symbols were removed, while genes with more than one probe set were averaged. |LogFC (fold change)| >1 and P value <0.05 were considered statistically significant. Principal component analysis (PCA) was performed and volcano pots were generated with online tools (https://www.omicstudio.cn/tool).

GO and KEGG analyses and GSEA

To identify functional changes in the DEGs of macrophages, biological analyses were performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID; http://david.ncifcrf.gov) online database. Analyses were performed via the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A Benjamini-Hochberg corrected P value <0.05 was considered significant. The GO and KEGG analyses were performed under R version 3.6.1. Gene set enrichment analysis (GSEA) was performed via online tools (https://www.omicstudio.cn/tool). A P value <0.05 was considered statistically significant.

Statistical analysis

We used GEO2R to identify DEGs, and the thresholds were set as |logFC| >1 and P<0.05. P<0.05 was considered statistically significant for the GO and KEGG analyses and GSEA.

Results

The role of macrophages in murine systemic lupus erythematosus nephritis

GSE27045 had three groups: a group of renal isolated macrophages taken at an early stage of lupus (the young mac group), a group of renal isolated macrophages taken during lupus nephritis (the sick mac group) and a group of renal isolated macrophages taken after induction of remission (the rem mac group). Six young mac samples, seven sick mac samples and four rem mac samples were assessed in the current study. PCA was used to compare the transcription profiles between the young mac group, sick mac group and rem group. The plot showed a clear distinction among the three groups (). Volcano plots were generated to visualize the distribution of expressed genes (P<0.05 and |logFC| >1, ) between the young mac group and sick mac group and between the sick mac group and rem group. Compared to the young mac group, the sick mac group showed altered expression of functional genes, which were annotated based on GO biological process terms (). The results indicated that the upregulated genes were enriched in immune system process, the innate immune response, the inflammatory response, the adaptive immune response, the defense response to bacterium, phagocytosis recognition, the cellular response to interleukin-1, phagocytosis, engulfment, the immune response and the B cell receptor signaling pathway. Altered signaling pathways were also demonstrated in the KEGG pathway analysis (), including pathways related to tuberculosis, cytokine-cytokine receptor interactions, measles, Staphylococcus aureus infection, phagosomes, the haematopoietic cell lineage, cell adhesion molecules (CAMs), transcriptional misregulation in cancer, the intestinal immune network for IgA production and the chemokine signaling pathway. As we wanted to focus on classical signaling pathways and metabolic processes, GSEA was applied to reveal the crucial signaling pathways and metabolism-related pathways of macrophages in SLE nephritis (); the results included the terms Toll-like receptor signaling pathway [Normalized Enrichment Score (NES) =2.2851; P=0.0019], JAK-STAT signaling pathway (NES =1.8629; P=0.0139), VEGF signaling pathway (NES =1.7747; P=0.0274), Hedgehog signaling pathway (NES =−1.6290; P=0.0288), MTOR signaling pathway (NES =1.7115; P=0.0381), systemic lupus erythematosus (NES =1.6549; P=0.0477), ERBB signaling pathway (NES =1.6447; P=0.0438), oxidative phosphorylation (NES =2.5644; P=0.0003), sphingolipid metabolism (NES =2.1394; P=0.0086) and glycerolipid metabolism (NES =2.1033; P=0.0092).

Figure 1

The characteristic and pathways of macrophages involved in murine systemic lupus erythematosus nephritis. (A) Principal component analysis (PCA) plot of the separation between the group of samples taken in the early-stage lupus (the young mac group), the group of samples taken during lupus nephritis (the sick mac group) and the group of samples taken after induction of remission (the rem mac group); (B) volcano plot of differentially expressed genes (DEGs) between the young mac group and sick mac group and between the sick mac group and rem group (P<0.05 and |logFC| >1); (C) the enriched biological process terms of the DEGs as identified via Gene Ontology (GO) functional annotation; (D) the enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway terms of the DEGs; (E) gene set enrichment analysis (GSEA) was performed to identify significantly enriched pathways.

The role of macrophages in murine renal crystal formation

GSE51466 had two groups, including a group of op/op (CSF-1-deficient) renal macrophage that received glyoxylate (GOX) injection (op/op group) for 6 days and a group of wild-type renal macrophage that received GOX injection for 6 days (wild-type group). Three op/op samples and three wild-type samples were assessed in the present study. The PCA results showed a clear distinction between the samples (). As shown in the volcano plots in , upregulated and downregulated DEGs (P<0.05 and |logFC| >1) were identified in the dataset. To obtain deeper insights into the biological function of the DEGs, we analysed the GO terms. The DEGs were most significantly enriched in biological process terms such as immune response, inflammatory response, chemokine-mediated signaling pathway, positive regulation of smooth muscle cell proliferation, chemotaxis, cellular response to lipopolysaccharide, response to molecule of bacterial origin, cell chemotaxis, positive regulation of NF-kappaB import into nucleus, and cellular response to interferon-gamma (). KEGG pathway analysis was performed to investigate the involved signaling pathways, and the results included pathway related to the cytokine-cytokine receptor interaction, rheumatoid arthritis, graft-versus-host disease, the NOD-like receptor signaling pathway, the TNF signaling pathway, influenza A, antigen processing and presentation, Staphylococcus aureus infection, complement and coagulation cascades, and type I diabetes mellitus (). Furthermore, by using GSEA (), we found enrichment of the ERBB signaling pathway (NES =1.8573; P=0.0001), JAK-STAT signaling pathway (NES =1.5519; P=0.0002), MAPK signaling pathway (NES =1.5043; P=0.0001), MTOR signaling pathway (NES =1.9725; P=0.0001), PPAR signaling pathway (NES =1.5758; P=0.0025), Toll-like receptor signaling pathway (NES =1.9680; P=0.0001), and VEGF signaling pathway (NES =1.7056; P=0.0002), which are all well-known pathways related to the suppressive role of macrophages in renal crystal formation. Additionally, other enriched pathways suggested the involvement of macrophages in metabolism-related mechanisms, such as oxidative phosphorylation (NES =1.9923; P=0.0001), purine metabolism (NES =1.5497; P=0.0001), fatty acid metabolism (NES =1.7238; P=0.0023), fructose and mannose metabolism (NES =1.8738; P=0.0012) and pyrimidine metabolism (NES =1.8466; P=0.0001).

Figure 2

The characteristic and pathways of macrophages involved in murine renal crystal formation. (A) Principal component analysis (PCA) plot of the separation between the group of wild-type renal macrophages treated with GOX injection for 6 days (the wild-type group) and the group of op/op (CSF-1-deficient) renal macrophage treated with glyoxylate (GOX) injection for 6 days (the op/op group); (B) volcano plot of the differentially expressed genes (DEGs) between the wild-type group and the op/op group (P<0.05 and |logFC| >1); (C) the enriched biological process terms of the DEGs as identified via Gene Ontology (GO) functional annotation; (D) the enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway terms of the DEGs; (E) gene set enrichment analysis (GSEA) was performed to identify significantly enriched pathways.

The role of macrophages after murine renal ischemia-reperfusion injury

GSE75808 had four groups, including a group that received sham surgery (sham group), a group with samples taken 4 h postperfusion after IRI (4 h group), a group with samples taken 24 h postperfusion after IRI (1-day group) and a group of samples taken 9 days postperfusion after IRI (9-day group). Nine sham samples, six 4-h samples, three 1-day samples and six 9-day samples were assessed in the present study. The PCA plot showed separation between the sham, 4 h after IRI, 1-day and 9-day samples (). These results were also displayed visually in the form of a volcano plot (P<0.05 and |logFC| >1, ). The BP terms () enriched by the DEGs were mainly associated with RNA splicing, protein transport, mRNA processing, 2-oxoglutarate metabolic processes, transport, ubiquinone biosynthetic processes, the response to stress, oxidation–reduction processes, the tricarboxylic acid cycle, and isocitrate metabolic processes. Moreover, KEGG pathway analysis () illustrated significant enrichment of pathways related to the biosynthesis of antibiotics, carbon metabolism, the citrate cycle (TCA cycle), glutathione metabolism, metabolic pathways, the TNF signaling pathway, the spliceosome, endocrine- and other factor-regulated calcium reabsorption, chemokine signaling pathways and protein processing in the endoplasmic reticulum. The significantly enriched KEGG pathways in 1-day vs. sham samples identified in the GSEA are displayed in and mainly included the ERBB signaling pathway (NES =2.0081; P=0.0053), MAPK signaling pathway (NES =2.0732; P=0.0001), PPAR signaling pathway (NES =2.0377; P=0.0071), TGFβ signaling pathway (NES =2.0060; P=0.0050), Toll-like receptor signaling pathway (NES =2.1165; P=0.0036), WNT signaling pathway (NES =1.7735; P=0.0069), pathways related to arginine and proline metabolism (NES =2.4682; P=0.0023), pathways related to the biosynthesis of unsaturated fatty acids (NES =2.2817; P=0.0043), pathways related to glycerolipid metabolism (NES =2.2273; P=0.0047) and pathways related to purine metabolism (NES =1.7802; P=0.0070).

Figure 3

The characteristic and pathways of macrophages involved in murine renal ischemia-reperfusion injury (IRI). (A) Principal component analysis (PCA) plot of the separation between the sham surgery (sham) group, the group of samples taken 4 h postperfusion after IRI (the 4 h group), the group of samples taken 24 h postperfusion after IRI (the 1-day group) and the group of samples taken 9 days postperfusion after IRI (the 9-day group); (B) volcano plot of the differentially expressed genes (DEGs) between the 1-day group and sham group (P<0.05 and |logFC| >1); (C) the enriched biological process terms of the DEGs as identified via Gene Ontology (GO) functional annotation; (D) the enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway terms of the DEGs; (E) gene set enrichment analysis (GSEA) was performed to identify significantly enriched pathways; (F) the Venn diagram represents common and unique pathways that were significantly altered in different datasets.

Furthermore, a Venn diagram () of the common and unique pathways that were significantly altered in different diseases was generated. Among the 47 pathways identified from all the datasets, the oxidative phosphorylation pathway, VEGF signaling pathway and JAK/STAT signaling pathway were common to all datasets. Of the two pathways unique to the GSE27045 dataset, the renin-angiotensin system pathway is a classical pathway and was thus considered worthy of further study. Of the 34 pathways specific to the GSE51466 dataset, the glycolysis and gluconeogenesis pathways are well known. Of the 13 pathways specific to the GSE75808 dataset, the calcium signaling pathway is a classical pathway.

Discussion

Increasing evidence shows that multiple signals are involved in macrophage activation in various diseases (2), such as signals from the JAK/STAT, PI3K/AKT, NF-κB, JNK, Wnt/β-catenin and Notch signaling pathways (19-23). Interestingly, strategies inducing metabolic reprogramming have great potential for regulating macrophages and in the development of novel approaches. Metabolic pathways that modulate the activation of macrophages include the glycolysis, TCA cycle, and fatty acid metabolism pathways, while metabolites such as succinate, α-ketoglutarate, itaconate, amino acids, adenosine and lactate also drive macrophage activation (24-26). In this regard, therapeutic strategies targeting signaling pathways and metabolic processes involved in the recruitment and activation of macrophages deserve attention for the treatment of kidney disease in the future. In the current study, we investigated the role of infiltrating macrophages in different renal diseases (including LN, renal crystal formation and IRI) by analysing classic signaling pathways and metabolic pathways. Three GEO datasets (GSE27045, GSE51466 and GSE75808) were included in our study.

In GSE27045, Bethunaickan et al. (12) showed that the traditional inflammatory macrophage population found in most acute renal inflammatory models was decreased in NZB/W mice and other murine SLE models. In our study, macrophages in LN were found to be mainly involved in pathways such as the ERBB signaling pathway, Hedgehog signaling pathway, JAK-STAT signaling pathway, MTOR signaling pathway, Toll-like receptor signaling pathway and VEGF signaling pathway. Furthermore, glycerolipid metabolism, oxidative phosphorylation and sphingolipid metabolism were found to be related to macrophages in LN. Moreover, compared with the pathways enriched in macrophages in renal crystal formation and IRI, the renin-angiotensin system pathway was found to be uniquely enriched in macrophages in LN. The renin-angiotensin system may influence immune responses and particularly affect macrophages through different signaling pathways. Hammer et al. (27) identified that the Ang-[1-7]/Mas axis is an important intracellular signaling pathway for macrophage function in vascular system and the central nervous system inflammation. The renin-angiotensin system impairs macrophage lipid metabolism in age-related macular degeneration in mouse models (28). The renin-angiotensin system also plays a role in LN. Yue et al. (29) revealed that blocking the renin-angiotensin system early improves renal outcomes of SLE patients with antiphospholipid-associated nephropathy. In a murine model, Matsuki-Muramoto et al. (30) verified that bortezomib treatment prevents glomerulosclerosis related to LN by suppressing the immune and renin-angiotensin systems. Thus, identifying enrichment of these key pathways, particularly the enrichment of the renin-angiotensin system, may provide new targets for modulating macrophages to further control the LN process.

In the gene array profiling of GSE51466, Taguchi et al. (16) compared crystal formation and macrophage polarization in renal tissue between wild-type and CSF-1-deficient mice after inducing hyperoxaluria. The results indicated a suppressive function of macrophages in renal crystal formation. Moreover, we demonstrated that the differentially altered mRNAs between macrophages taken from the kidneys of wild-type or CSF-1-deficient mice were mainly involved in the ERBB signaling pathway, JAK-STAT signaling pathway, MAPK signaling pathway, MTOR signaling pathway, PPAR signaling pathway, Toll-like receptor signaling pathway, and VEGF signaling pathway. In addition, the oxidative phosphorylation, purine metabolism, fatty acid metabolism, fructose and mannose metabolism and pyrimidine metabolism pathways were identified as enriched metabolic pathways that may participate in the suppressive role of macrophages in renal crystal formation. Furthermore, compared with the pathways associated with macrophages in LN and IRI, glycolysis and gluconeogenesis were found to be activated in macrophages during renal crystal formation. According to stimuli from the microenvironment, macrophages can switch from aerobic oxidative phosphorylation to high rates of anaerobic glycolysis (31), and proinflammatory macrophages utilize glycolysis and the pentose phosphate pathway to meet their ATP requirements (32). The role of glycolysis in macrophages during renal crystal formation still needs further study.

With the GSE75808 dataset, Sun et al. (18) revealed that the proportion of macrophages peaked one day after IRI and was significantly decreased at Day 9 after IRI. Thus, we chose the macrophage profile data from one day after IRI to compare with those of the sham group. Studies have also shown that Atp5a1, Atp5o, Cdc42, Cox4i, Nhp2 and Rac2 are critical genes involved in oxidation-reduction, migration, apoptosis, and macrophage proliferation and polarization (18). We focused on specific signaling pathways to explore novel targets that can influence AKI. Compared to the sham group, the group with macrophages taken from the kidney 1 day after IRI showed significant enrichment of pathways including the ERBB signaling pathway, MAPK signaling pathway, PPAR signaling pathway, TGFβ signaling pathway, Toll-like receptor signaling pathway and WNT signaling pathway. Regarding metabolic pathways, the arginine and proline metabolism pathways, pathways related to the biosynthesis of unsaturated fatty acids, the glycerolipid metabolism pathway and the purine metabolism pathway were identified. Interestingly, compared with the pathways enriched in macrophages in LN and renal crystal formation, the calcium signaling pathway was found to be specifically enriched in macrophages in IRI. Calcium has been found to be central to various physiological processes, such as immune system activation and maintenance (33). Intracellular Ca2+ fluctuations act on many cellular functions in monocytes/macrophages, including chemotaxis and phagocytosis (34). Ca2+ also plays crucial roles in actin rearrangement and the activation of NADPH oxidase and protein kinase C in phagocytic cells like macrophages (35). In a mouse model of AKI, Wasnik et al. (36) demonstrated that blocking intracellular calcium improved survival. In primary tubular epithelial cells, TRPC6-mediated Ca2+ influx was found by Hou et al. (37) to suppress cytoprotective autophagy triggered by oxidative stress by downregulating PI3K/AKT and ERK signaling, increasing the autophagy level, and restraining the apoptosis mediated by mitochondria. These studies and results suggest that understanding the abovementioned signaling pathways, especially the calcium signaling pathway, in macrophages during IRI may provide new targets for treatment.

In summary, our study identified transcriptional features of macrophages in murine LN, renal crystal formation and IRI models. We focused on the classic signaling pathways and metabolic pathways of macrophages infiltrating the kidney during these pathophysiological processes. These results help uncover the comprehensive role of macrophages in multiple renal diseases at the transcript level, which could provide vital information for reversing different renal diseases. However, there are some limitations of this study. First, the identified signaling pathways require verification for further applications. Second, the abilities of macrophages in various renal diseases need to be further validated.

Conclusions

In this investigation, our data based on transcriptional profiles of macrophages provide a deeper understanding of the classic signaling pathways and metabolic pathways activated in renal infiltrating macrophages in LN, renal crystal formation and IRI models. We also explored the similarities and differences in macrophage function among these different renal diseases. Oxidative phosphorylation, the VEGF signaling pathway, and the JAK/STAT signaling pathway were found to be activated in macrophages in all the three renal diseases; the renin-angiotensin system was found to be activated in macrophages in LN; glycolysis and gluconeogenesis were found to be activated in macrophages during renal crystal formation; and the calcium signaling pathway was found to be activated in macrophages in renal IRI. These specific pathways deserve further study and may represent promising therapeutic targets.

The article’s supplementary files as