- Open Access
SGLT2 inhibitors improve kidney function and morphology by regulating renal metabolic reprogramming in mice with diabetic kidney disease
Journal of Translational Medicine volume 20, Article number: 420 (2022)
Diabetic kidney disease (DKD) is the leading cause of end-stage renal disease (ESRD) worldwide. SGLT2 inhibitors are clinically effective in halting DKD progression. However, the underlying mechanisms remain unclear. The serum and kidneys of mice with DKD were analyzed using liquid chromatography with tandem mass spectrometry (LC–MS/MS)-based metabolomic and proteomic analyses. Three groups were established: placebo-treated littermate db/m mice, placebo-treated db/db mice and EMPA-treated db/db mice. Empagliflozin (EMPA) and placebo (10 mg/kg/d) were administered for 12 weeks. EMPA treatment decreased Cys-C and urinary albumin excretion compared with placebo by 78.60% and 57.12%, respectively (p < 0.001 in all cases). Renal glomerular area, interstitial fibrosis and glomerulosclerosis were decreased by 16.47%, 68.50% and 62.82%, respectively (p < 0.05 in all cases). Multi-omic analysis revealed that EMPA treatment altered the protein and metabolic profiles in the db/db group, including 32 renal proteins, 51 serum proteins, 94 renal metabolites and 37 serum metabolites. Five EMPA-related metabolic pathways were identified by integrating proteomic and metabolomic analyses, which are involved in renal purine metabolism; pyrimidine metabolism; tryptophan metabolism; nicotinate and nicotinamide metabolism, and glycine, serine and threonine metabolism in serum. In conclusion, this study demonstrated metabolic reprogramming in mice with DKD. EMPA treatment improved kidney function and morphology by regulating metabolic reprogramming, including regulation of renal reductive stress, alleviation of mitochondrial dysfunction and reduction in renal oxidative stress reaction.
Diabetes mellitus (DM) is a common chronic metabolic disease characterized by persistent hyperglycemia caused by absolute or relative insulin deficiency and insulin resistance. The estimated number of patients with DM is projected to increase to 700 million by 2045 worldwide . In addition, DM is the leading cause of chronic kidney disease (CKD) and end-stage renal disease (ESRD). Approximately 40% of patients with DM gradually progress to CKD and even develop ESRD . Unfortunately, the treatment of DKD remains an unresolved challenge worldwide. To date, strict glycaemic and blood pressure control, especially with the use of angiotensin-converting enzyme inhibitors and angiotensin II receptor blockers, has been shown to delay the progression of DKD [3, 4]. However, the benefit of these interventions to DKD outcomes is limited . Therefore, developing new treatment strategies to delay the progression of DKD is necessary.
The pathogenesis of DKD is multifactorial and still unclear. Accumulating evidence elucidates that various cellular stress responses triggered by metabolic abnormalities contribute to the pathophysiological mechanism of DKD . The upregulation of sodium-glucose cotransporter 2 (SGLT2), when appearing in diabetes, leads to high transmembrane transport in the proximal tubule and generates a dramatic energy requirement in the renal cortex . Meanwhile, metabolic reprogramming characterized by tricarboxylic acid cycle (TCA cycle) inhibition and glycolysis enhancement occurs in the renal cortex . This phenomenon of mitochondrial oxidative phosphorylation activity mismatched with increased glycolytic flux results in the accumulation of intermediates from glucose metabolism and triggers multiple pathogenic signaling pathways, thereby inducing oxidative stress, reductive stress, and various cellular stress responses mediated by energy insufficiency [9, 10]. Moreover, insulin resistance shifts the pattern of renal biofuel utilization from glucose-mediated aerobic metabolism to fatty acid β-oxidation and amino acid metabolism . Reversing these metabolic abnormalities has demonstrated encouraging results in improving the outcomes of DKD. For example, inhibited KIM-1-mediated fatty acid uptake in renal tubular epithelial cells can block a series of kidney injury-related events triggered by abnormal lipid metabolism in DKD ; similarly, the reversal of abnormal glycolysis and lipid metabolism in DKD was protective against the pathophysiology of DKD .
SGLT2 inhibitors originally developed for treating type 2 Diabetes mellitus (T2DM) are clinically very effective in halting the progression of DKD [12,13,14]. Recent studies have demonstrated that SGLT2 inhibitors exert both direct or indirect protective effects on the cardiovascular and renal systems in T2DM, which are independent of their glucose control effects but related to their effects on blood pressure control, glomerular hemodynamic amelioration, RAAS regulation, and anti-inflammation and in reducing glucotoxicity, lipotoxicity, and uric acid. Furthermore, SGLT2 inhibitors downregulate hepcidin and promote erythropoiesis [15, 16], thereby alleviating renal hypoxia by improving the circulating oxygen supply . Interestingly, a recent study has demonstrated that dapagliflozin prevents the high glucose-induced metabolic transition shift from lipid oxidation to glycolysis in renal tubular cells by inhibiting HIF-1α . In contrast, empagliflozin inhibit the over-activated nutrient-sensing signaling pathway by increasing endogenous ketone body production, thus improving energy utilization and reducing renal oxidative stress . This suggests that the SGLT2 inhibitor effectively corrects metabolic disorders in DKD. However, the mechanisms behind the renoprotective effect of SGLT2 inhibitors are unclear.
This study was aimed to investigate the metabolic effects of SGLT2 inhibitors in DKD via liquid chromatography with tandem mass spectrometry (LC–MS/MS)-based metabolomic and proteomic analyzes of serum and the kidneys, and to provide a deeper understanding of the renoprotective mechanism of SGLT2 inhibitors.
Research design and methods
Diabetic kidney disease model
Six-week-old male db/db mice (C57BLKS/J-leprdb/leprdb, n = 30) and matched littermate db/m mice (C57BLKS/J-leprdb/leprm, n = 15, as normal controls) were purchased from the Nanjing University Experimental Animal Center. After 2 weeks of adaptive feeding, db/db mice were randomly divided into the db/db (model) and empagliflozin (EMPA; EMPA-treated) groups (n = 15). Mice in the EMPA group were administered 10-mg/kg EMPA per day via oral gavage for 12 weeks. Similarly, mice in the model and control groups were administered 10-mL/kg sterile water was dosed to the mice in the model and control group in the same way. All mice were housed in a controlled SPF-grade environment (temperature, 22 ± 2 ℃; humidity, 55 ± 10%; with a 12-h light/dark cycle) without dietary and drinking water restriction. At 3 weeks and every 4 weeks thereafter, whole blood was collected from the tail veins of mice after overnight fasting, and their fasting blood glucose levels were measured using a portable glucometer. All mice were weighed once a week beginning at the third week. In addition, 24-h urine samples were obtained from each mouse using metabolic cages after 12 weeks of treatment and stored at − 80 ℃ after they were snap-frozen in liquid nitrogen. Subsequently, the mice were anesthetized using 2,2,2-tribromoethanol via intraperitoneal injection, and whole blood samples were obtained via cardiac puncture. The samples were incubated at room temperature for 1 h in procoagulant tubes and centrifugated at 3000 r/min for 15 min to collect the upper serum. For kidney specimens, one side of the kidney was snap-frozen in liquid nitrogen after weighting, whereas the other side was fixed in 4% paraformaldehyde overnight and subsequently embedded in paraffin. All samples were stored at − 80 ℃ before performing LC–MS/MS. All experiments in this study were conducted according to the requirements of the National Law for Laboratory Animal Experimentation and were approved by the Commission on Experimental Animal Ethics of Jinan University (No. 202069-04). A detailed flowchart of animal experiments is shown in Fig. 1A.
Biochemical analysis and histopathological examination
The fasting blood glucose level was detected using a portable glucometer (BAYER, Germany) with blood glucose test strips (BAYER, Germany). The serum and urine samples were thawed at 4 ℃ before detection. Serum cystosin C (Cys-C) levels were measured with a mouse Cys-C ELISA kit (E-EL-M3024, Elabscience, Wuhan, China). Urine creatinine (UCR) and urine microalbuminuria albumin levels were measured using a UCR enzyme-linked immunoassay kit (MM-44289M1, Elabscience, Wuhan, China) and MAU/ALB enzyme-linked immunoassay kit (MM-0705M1, Elabscience, Wuhan, China), respectively.
For renal histopathological analysis, paraffin-embedded renal tissues were cut into 4-um sections and subjected to haematoxylin–eosin (H&E), Masson and PAS staining. Subsequently, 10 randomly selected non-overlapping areas (25 × 25 μm) on HE-stained images (200×) were used to evaluate the extent of glomerular damage by calculating the average glomerular perimeter and area using the ImageJ software. The same method was used to select 10 fields of view on Masson-stained images and determine the percentage of the fibrotic area (blue-stained area) in each observed field for assessing the degree of glomerular and interstitial fibrosis. Renal histopathological analysis was performed by two independent investigators using a blinded method.
Proteomic analysis was conducted as reported previously . Briefly, an appropriate amount of renal tissue solubilised with 300 uL of RIPA lysis buffer was homogenised before sonication. Thereafter, the homogenates were centrifuged at 4 ℃ and 12,000 rpm for 10 min to collect the supernatant. The supernatant was precipitated using 80% pre-cold acetone and washed thrice, and the precipitates were redissolved in protein resolve buffer (HEPESSDC1). Furthermore, serum samples were incubated in a water bath and used to extract proteins using the BIO-RAD de-high abundance protein kit. Thereafter, we determined protein concentration using a BCA kit. Subsequently, 100-ug total protein from each sample was diluted 5 times with 50-mM NH4HCO3 and reduced and alkylated with DTT (5 mM) and IAA (10 mM), respectively. Thereafter, the sample was mixed with 0.5-ug/uL trypsin (v:v, 50:1) and incubated overnight at 37 ℃ to obtain peptide solution for further LC–MS/MS analysis, and the reaction was quenched by adding 1-uL TFA.
LC–MS/MS for proteomic analysis
We used a nano-UPLC system (EASYnLC1200) in combination with a mass spectrometer (QExactive HFX) equipped with a nano-litre ion source for proteomic data acquisition. A 100-um ID × 15 cm reversed-phase column (ReprosilPur 120 C18AQ, 1.9 um, Dr. Math) was used for chromatographic separation. Mobile phase A contained 0.1% formic acid aqueous solution with 2% acetonitrile, whereas mobile phase B contained 0.1% formic acid aqueous solution with 80% acetonitrile, and the flow rate was 300 nL/min. MS was performed based on data-dependent acquisition (DDA) in the positive ion mode. The parameters were set as follows: full scan range, 350–1600 m/z; resolution, 120 K; AGC target, 3E6; maximum ion implantation time, 50 ms. The top 20 highest-intensity ions were fragmented via HCD with the normalised crash energy (NCE) of 27%, and the quadrupole isolation window was 1.2 m/z. The dynamic exclusion time was set to 45 s, and single-charged peaks and peaks with charges > 6 were excluded from DDA.
Subsequently, raw MS data were processed using the Proteome Discoverer (PD) software (Version 188.8.131.525) and the plug-in Sequest HT search engine. The MS data were matched according to the species-level UniProt FASTA databases (uniport-Mus + musculus-10090-2020-10.fasta). Carbaminomethyl (C) was considered a fixed modification, whereas oxidation (M) and acetyl (protein N term) were considered variable modifications. Trypsin was used as a protease. Peptide tolerance was set to 10 ppm, and MS/MS tolerance was set to 0.02 Da. Only two missed cleavages were acceptable. The PSM and peptide levels were set as the false discovery rate (FDR) of 0.01. The quantification of proteins and normalisation of the total peptide amount were performed based on the unique and razor peptides. Other parameters were maintained as default. The detailed process of metabolomic analysis is described below.
For metabolite extraction, 25 mg of each kidney sample was taken and mixed with 500 uL of the extract solution (methanol:acetonitrile:water = 2:2:1, with isotopically-labelled internal standard mixture). The samples were homogenised at 35 Hz for 4 min and subsequently sonicated in an ice-water bath for 5 min. These steps were repeated twice or thrice. The serum samples were thawed at room temperature, and 50 uL of each sample was mixed with 200 μL of the extract solution in an EP tube, vortexed for 30 s and sonicated in an ice-water bath for 5 min. Thereafter, the processed kidney and serum samples were incubated at − 40 ℃ for 1 h and centrifuged at 12,000 rpm for 15 min at 4 ℃. The supernatant was collected in a sample bottle for further detection. In addition, an equal volume of all samples was stored for subsequent quality control (QC) of LC–MS/MS analysis.
A UHPLC system (Vanquish, Thermo Fisher Scientific) with UPLC BEH Amide (2.1 mm × 100 mm, 1.7 μm) combined with a Q Exactive HFX mass spectrometer (Orbitrap MS, Thermo) was used for LC–MS/MS analysis. Mobile phase A comprised 25-mmol/L ammonium acetate and 25 ammonia hydroxide in water (pH = 9.75); mobile phase B was acetonitrile. The temperature of the auto-sampler was maintained at 4 ℃, and the injection volume was 3 μL. All data, including MS and MS/MS data, were obtained using the QE HFX mass spectrometer under the control of an acquisition software (Xcalibu, Thermo). The detailed parameters of ESI source conditions were as follows: sheath gas flow rate, 30 Arb; aux gas flow rate, 25 Arb; capillary temperature, 350 ℃, full MS resolution, 60,000; MS/MS resolution, 7500; collision energy, 10/30/60 in the NCE mode; spray voltage, 3.6 kV (positive) or − 3.2 kV (negative).
The ProteoWizard software was used to convert the format of the raw data to mzXML. The converted data were processed using an in-house R package based on XCMS for peak detection, extraction, alignment and integration. Subsequently, metabolites were annotated by matching with a secondary MS database based on the self-built BiotreeDB (V2.1) with the cutoff for annotation set to 0.3. Missing values in raw data were filled using half of the minimum value, and all data were normalised using the total ion intensity (TIC) of every single sample for the convenience of further analysis.
Normalization and imputation of proteomic and metabolomic analyses
For the raw proteomic and metabolomic data, the relative intensity of each protein was normalised by log2 conversion, whereas the relative intensity of each metabolite was normalised based on the total ion current (TIC) of each metabolite. A protein or metabolite was discarded if it was detected in < 50% samples. Otherwise, the protein or metabolite was retained, and the missing value was replaced with half of the minimum positive value of each variable.
Multivariate statistical analysis was performed, in which unsupervised principal component analysis (PCA) was conducted to observe the metabolic profile of all subjects, identify abnormal outliers and assess the stability of QC samples. In addition, supervised partial least square-based discriminant analysis (PLS-DA) was used for characterising metabolic feature separation among different groups and screening biomarkers. Q2 and R2 were the parameters for evaluating the predictability and interpretability and testing the fit of the PLS-DA model via leave-one-out cross-validation. The variable importance in projection (VIP) score was calculated based on PLS-DA to quantify the contribution of metabolites in the model and was used to screen for differentially expressed metabolites. Metabolites with VIP scores > 1, P-value < 0.05 (Student’s t‐test) and |fold change|> 1.5 were considered significantly differentially expressed metabolites (DEMs). Similarly, proteins with |fold change|> 1.5 and adjusted P-value < 0.05 were considered differentially expressed proteins (DEPs). The Short Time-series Expression Miner (STEM) program was used to identify target proteins regulated by EMPA. A Venn diagram, volcano plots and a hierarchically clustered heat map were generated using the R packages ‘ggplot’ and ‘Complexheatmap’ to visualise the feature of identified DEPs and DEMs. Gene Ontology (GO) enrichment analysis and KEGG enrichment analysis were used to identify the potential pathways of differentially expressed proteins and metabolites and EMPA-altered proteins and metabolites. To quantify the overall expression of each enriched metabolic pathway, the average intensity of all metabolites in a pathway with at least three metabolites was assessed and used to calculate the fold change of each pathway, which indicated the activity of the pathway . The identified metabolites were classified according to the HMDB database, and the percentage of each metabolite category was calculated. Thereafter, the Fisher’s exact test or chi-square test was used to evaluate changes in the percentage of each metabolic category before and after EMPA intervention.
Biochemical indicators were analyzed using SPSS Statistics version 25.0 (IBM, USA), GraphPad Prism 9 and R software (R Foundation for Statistical Computing, Vienna, Austria, v3.5.3). The Student’s t-test was used for comparing data between two groups. Normally distributed data were expressed as mean ± SD. Comparison among multiple groups and pairwise comparison were performed using one‐way ANOVA and LSD test, respectively. The relationship between the serum and kidney in both metabolomic and proteomic analyses was assessed using Pearson correlation analysis. A P-value of < 0.05 was considered statistically significant.
Effects of EMPA on biochemical indicators and pathological characteristics in db/db mice
The flowchart of the animal experiment is shown in Fig. 1A. Mice in the db/db and EMPA groups had higher body weight than mice in the db/m group, and no difference was observed in body weight between the db/db and EMPA groups (Fig. 1B). EMPA treatment improved fasting blood sugar (FBG), blood urea nitrogen (BUN), Cys-C and serum creatinine levels and urine albumin-to-creatinine ratio (P < 0.05) (Fig. 1C–G).
Subsequently, renal histopathological examination (Fig. 1H) revealed that db/db mice had conspicuous pathological alterations, such as glomerular dilatation, mesangial cell hypercellularity, mesangial matrix hyperplasia and tubular basement membrane thickening, which were corrected after EMPA treatment (Fig. 1H, a–f). As shown in Fig. 1I and J, EMPA treatment normalized the average glomerular perimeters (P = 0.070) and glomerular area (P = 0.045) in db/db mice, suggesting that EMPA mitigates glomerular hypertrophy in DKD. In addition, Masson staining showed a significant increase in renal fibrosis in db/db mice (Fig. 1H, g–l). EMPA intervention reduced the level of glomerular and interstitial fibrosis by 62.8% and 68.5%, respectively (Fig. 1K, L). These results indicate that EMPA can reduce urinary protein excretion, improve renal function and reverse the renal pathological alterations resulting from DKD to some extent, which was consistent with a previous study .
Effects of EMPA on the protein expression profile of db/db mice
A total of 6009 renal proteins were detected in 6 db/m, 6 db/db and 5 EMPA samples using the DIA-based non-labeled quantitative method (Additional file 1: Fig. S1A, Additional file 2: Table S1). Of the 6009 renal proteins, 5010 (83.4%) were supported by at least 2 peptides, with an average number of 10 peptides (Additional file 1: Fig. S1B). In addition, 84% of proteins (5046) were detected in 17 kidney samples (Additional file 1: Fig. S1C). Similarly, 1481 serum proteins were detected in 6 db/m, 5 db/db and 6 EMPA samples (Additional file 1: Fig. S1D, Additional file 2: Table S2); of which, 74% of proteins were supported by > 2 peptides, and 762 proteins were co-expressed in 17 serum samples (Additional file 1: Fig. S1E). Therefore, the quantification results of the proteomic and metabolomic analyses of the kidney and serum are highly credible.
In renal proteomic analysis, 501 DEPs were identified in the db/db group (Additional file 2: Table S3), and metabolism-related proteins (e.g. ugt1a9, ugt1a2, gsta2, acsf2, npl, acsf2, cbr1, acsm3, amacr and acy3) dominated the top altered proteins, revealing dramatic changes in metabolic processes in the kidney of db/db mice (Fig. 2A). To identify target proteins regulated by EMPA, we first identified 91 DEPs in the Empa group (Additional file 2: Table S3), and 576 trend proteins that varied between groups (Fig. 2B); of which, 32 proteins were significantly altered in the kidneys of db/db mice and the effects were considerably reversed after EMPA treatment(Fig. 2B and C, Additional file 2: Table S4). Among these proteins, high levels of fibrinogen proteins (Fgg and Fga), cathepsin proteins (ctss and ctsa), H2-Ab1 and prcp were considered risk factors for diabetic nephropathy and the predictors of urinary protein progression [22,23,24,25,26]. Pathway and GO enrichment analyses showed that altered proteins in the db/db group were significantly enriched in metabolic pathways, the PPAR signaling pathway and the renin–angiotensin system (Fig. 2D). However, EMPA therapy remarkably altered the metabolic process and immune response in the kidney of db/db mice (Fig. 2E).
Simultaneously, we detected 1481 proteins in serum samples, including 235 and 83 DEPs in the db/db and EMPA groups, respectively (Fig. 3A, Additional file 2: Table S5). Among the intersection of altered serum proteins in the db/db and EMPA groups, 33 downregulated and 18 upregulated proteins were identified in the db/db group, whose expression was significantly reversed after EMPA therapy (Fig. 3B, Additional file 2: Table S6). Among these proteins, decreased levels of cavin2, rap1b and prdx2 were reported to be associated with insulin resistance , renal tubular mitochondrial dysfunction  and oxidative damage , respectively, whereas evaluated levels of Agt and s100a9 showed significant association with diabetic nephropathy  and diabetic retinopathy , respectively. In the db/db group, 20 proteins, including known diabetes-associated proteins such as Itih2, Apoc1, Timp3, B2M, Adh1, Bhmt, Apoh, C4bpa, Ctsd, Ctsb, Lgals3 and Cyb5a, were significantly altered in both serum and kidney samples, whereas only three identical proteins (Igk-V19-17, Eif5 and Ighg1) were detected in the EMPA group (Fig. 3A). This finding suggested that the drug effects of EMPA on the renal and circulating protein profiles were mutually independent. Furthermore, correlation analysis showed no significant correlation between circulating and renal tissue proteins in the db/db and EMPA groups (Fig. 3C). Pathway analysis revealed that altered serum proteins in the db/db group were significantly enriched in metabolic pathways, complement and coagulation cascades and leukocyte transendothelial migration (Fig. 3D). EMPA treatment substantially changed the amino acid metabolism and inflammatory response in the circulatory system of db/db mice (Fig. 3E). These results indicated that EMPA exhibited strong effects on metabolic processes in db/db mice.
To systematically characterize the metabolic profile in DKD and the potential metabolic alterations after EMPA treatment, we performed non-targeted metabolomic analysis of the serum and kidney tissue in three experimental groups. We identified 897 (positive ion mode: 666 and negative ion mode: 231) and 678 (positive ion mode: 457, negative ion mode: 221) metabolites from the kidney and serum samples, respectively (Additional file 2: Tables S7 and S8).
In both positive and negative ion modes, extracted ion chromatograms (EIC) for the internal standard (IS) in quality control (QC) samples showed great consistency in retention time and response intensity (Additional file 1: Fig. S2A and B), indicating good stability for the acquisition of MS data in this study. In addition, Pearson correlation analysis revealed a high correlation among all QC samples (Additional file 1: Fig. S2C and D). Furthermore, three ISs were monitored, and their retention time and mass-to-charge ratio had satisfactory RSD (RSD median ≤ 15%) (Additional file 1: Fig. S2E and F). These results indicated that the sample quality, experimental method and system stability were reliable and suitable for metabolomic analysis of the kidney and serum.
Overall metabolic alterations in db/db mice
PLS-DA models showed noticeable differences in the distribution between db/db and db/m groups in the metabolomic analysis of both kidney (Q2 = 0.902, R2 = 0.999) and serum (R2 = 0.999, Q2 = 0.898) (Fig. 4A , B). According to the screening criteria of VIP > 1, P < 0.05 and |fold change|> 1.5, 389 and 174 DEMs were identified between the db/m and db/db groups in the kidney and serum, respectively. Of these DEMs, 264 renal and 174 serum metabolites were finally matched by comparing the KEGG and HDMB databases (Fig. 4C , D, Additional file 2: Tables S9 and S10).
To assess the aberrant alterations in metabolic pathways in DKD, we submitted 264 and 174 DEMs to Metaboanalyst 5.0 for KEGG pathway enrichment analysis. DEMs in the kidney were significantly enriched in 26 pathways (P < 0.05 and impact > 0.1; Fig. 4E, Additional file 2: Table S10), involving five types of metabolic abnormalities, including energy metabolism (glucose, lipid and amino acid metabolism), cofactor and vitamin metabolism, nucleotide metabolism and metabolism of other amino acids. The results indicated a reconfiguration of overall metabolism in db/db mice. Subsequently, we calculated the fold change of enriched pathways with at least three metabolites to assess their activity in the db/db and db/m groups (Fig. 4F). We found that the metabolic pathways of various amino acids, including Try, alanine, aspartic acid, glutamic acid, cysteine, methionine, arginine, proline, histidine, beta-alanine, glycine, serine, threonine and tyrosine metabolism; taurine and hypotaurine metabolism; lysine degradation and arginine biosynthesis, were altered in the kidney and serum, which highlighted the role of protein–energy wasting in DKD pathogenesis. In addition, nucleotide metabolism (e.g. purine and pyrimidine metabolism) and cofactor and vitamin metabolism (e.g. nicotinate and nicotinamide metabolism) were also altered in db/db mice.
Subsequently, we established a 3D PLS-DA model to examine the metabolic profiles of groups and observed a significant separation trend among them (Figs. 2H and 4G), which indicated that the construction of the DKD model was successful, and db/db mice had evident metabolic disorders. However, EMPA treatment reversed the changes resulting from DKD.
EMPA ameliorated kidney metabolism
We constructed a PLS-DA model to compare the metabolic characteristics between the db/db and EMPA groups in the kidney (Q2 = 0.54, R2 = 0.848) and identified 194 DEMs with the screening conditions mentioned earlier (Additional file 1: Fig. S3A and 3B, Additional file 2: Table S11). Of the 194 DEMs, 94 metabolites, co-expressed with 264 DEMs, were identified as EMPA-altered metabolites (EAMs) (Additional file 1: Fig. S3B); of which, the expression of 41 upregulated and 53 downregulated metabolites in the db/db group were significantly reversed after EMPA treatment (Fig. 5A, Additional file 2: Table S12). In addition, the expression of renal damage-related metabolites (e.g. l-phenylalanine, l-arginine and d-serine) was downregulated after EMPA treatment, whereas the expression of renal protection-related metabolites (e.g. adenosine, adenine and acetylcysteine) was increased in the kidney (Fig. 5A) [32,33,34,35,36].
Metabolism-related proteins (e.g. ugt1a9, ugt1a2, gsta2, acsf2, npl, acsf2, cbr1, acsm3, amacr and acy3) dominated the top altered proteins, revealing dramatic changes in metabolic processes in the kidney of db/db mice.
Subsequently, we classified 264 DEMs and 94 EAMs in the kidneys according to the HMDB database and used the chi-square test or Fisher’s exact test to estimate changes in the proportion of metabolite categories before and after EMPA treatment (Fig. 5B). We found that elevated proportion of the ‘nucleoside, nucleotides and derivatives’ (from 15.2% to 22.3%, P = 0.045) and reduced proportion of the ‘lipids and lipid-like molecules’ (from 24.7% to 11.7%, P = 0.003) categories were statistically significant.
Furthermore, hierarchical cluster analysis was performed on 94 renal EAMs (Fig. 5C). In the ‘nucleoside, nucleotides and derivatives’ category, most downregulated metabolites were upregulated after EMPA treatment, when compared to the db/db and db/m groups. KEGG enrichment analysis revealed that these nucleotide metabolites were significantly enriched in pathways related to purine and pyrimidine metabolism (Fig. 5D). In the ‘lipids and lipid-like molecules’ category, the expression of three fatty acids and three phospholipids (including two lysoPE) was elevated in db/db mice but decreased after EMPA treatment. In addition, the expression of two types of acylcarnitines and three fatty acyl classes (especially acetyl-CoA, which is the hub of energy metabolism in the body) was increased after EMPA intervention.
The ‘amino acids and derivatives’ category accounted for the largest proportion both before and after treatment. Among these amino metabolites, the proportion of dipeptides was the highest. The most noticeable change was observed in the expression of valyl-tyrosine (Val-Tyr, from downregulation to upregulation, FC [db/db versus db/m] = 0.24 and FC [EMPA versus db/db] = 6.94). In addition, there were 10 types of protein amino acids and 14 types of derivatives of α-amino acids in EMPA regulation. Five N-acyl-alpha amino acids, except for N-acetyl-l-alanine, were upregulated after EMPA treatment. It is noteworthy that N-acetylcysteine is a direct precursor to glutathione. We also observed that the expression of two known uremia toxins, argininosuccinic acid and impic dimethylarginine, was reduced after EMPA treatment (Additional file 2: Table S12).
We further performed KEGG pathway enrichment analysis on 94 EAMs in the kidney and assessed the activity of enriched pathways (Fig. 5D, E). The overall activity of purine metabolism (log2[FC] = 2.278, P = 0.009), nicotinate and nicotinamide metabolism (log2[FC] = 2.081, P = 0.0004), pyrimidine metabolism (log2[FC] = − 0.395, P = 0.0118) and tryptophan metabolism (log2[FC] = − 1.151, P = 0.006) was altered after EMPA treatment, as compared with the model group.
EMPA ameliorated serum metabolism
We analyzed alterations in serum metabolism in db/db mice treated with EMPA. The results of PLS-DA in serum were consistent with metabolic alterations observed in the kidney (R2 = 0.994, Q2 = 0.483) (Additional file 1: Fig. S3C). A total of 174 serum metabolites were identified as DEMs between the db/db and EMPA groups using the same screening criteria mentioned earlier. Of the 174 metabolites, 37 (DEMs in the db/m and db/db groups) were identified as EAMs in serum (Additional file 1: Fig. S3D). As demonstrated by the volcano plot, the expression of 13 EAMs in serum was increased and that of 14 EAMs was decreased after EMPA treatment (Fig. 6A). In addition, only 8 EAMs were significantly altered in both serum and kidney samples (Fig. 6B). The correlation analysis revealed no significant relationship between serum and renal metabolites in the db/db and EMPA groups, suggesting the mutually independent effects of EMPA in renal and serum metabolic profiles (Fig. 6C).
Similarly, we classified 174 DEMs and 37 EAMs in serum according to the HMDB database and performed the Fisher’s exact test to distinguish the percentage change in metabolite classification. The results revealed that the proportion of the ‘amino acids and derivatives’ category was increased (from 27.6 to 45.9%, P = 0.051) after EMPA treatment (Fig. 6D) and was highest before and after treatment. In comparison, the proportion of other metabolite categories was not significantly different before and after treatment. A heat map demonstrating hierarchical clustering of 37 TRMs showed that the expression of all metabolites in the ‘amino acids and derivatives’ category was increased after EMPA treatment (Fig. 6E). The proportion of the ‘amino acids and derivatives’ category was ranked in order of ‘derivatives of α amino acids’ (35%, n = 6), ‘dipeptide’ (24%, n = 4), ‘protein amino acids’ (18%, n = 3), ‘N-acyl-alpha amino’ (12%, n = 2) and ‘non-protein amino (12%, n = 2)’.
Subsequently, KEGG pathway enrichment analysis was performed and the activity state of pathways associated with 37 serum EAMs were predicted. We observed that all significantly altered pathways in serum were related to amino metabolism, such as glycine, serine and threonine metabolism (P = 0.017, impact = 0.207); arginine and proline metabolism (P = 0.029, impact = 0.188); aminoacyl-tRNA biosynthesis (P = 0.0005, impact = 0.1667); beta-alanine metabolism (P = 0.011, impact = 0.153) and histidine metabolism (P = 0.017, impact = 0.123). In these significantly enriched pathways, the overall levels of metabolites in ‘glycine, serine and threonine metabolism’ (log2[FC] = 0.765) and ‘aminoacyl-tRNA biosynthesis’ (log2[FC] = 0.688) pathways showed a trend opposite to that of the db/db group after EMPA intervention (Figs. 4E and 6G).
Potential mechanisms of EMPA in protection against DKD
As described above, we eventually identified four metabolic pathways associated with EMPA therapy in the kidney, including purine, nicotinate, nicotinamide, pyrimidine, tyrosine and β-alanine metabolism. These metabolism pathways are associated with the tricarboxylic acid (TCA) cycle through acetyl-CoA (Fig. 7). The ‘glycine, serine and threonine metabolism’ pathway was the most significantly altered pathway that responded to EMPA treatment in serum. It was related to the conformational transition between l-serine and d-serine. Altogether, we speculate that these metabolic pathways may play a critical role in the progression of DKD, and EMPA may exert beneficial effects by regulating these pathways.
In this work, we observed that Empa treatment effectively alleviates the renal pathological changes caused by T2DM; simultaneously, we characterized the specific feature of metabolic reprogramming that occurs in kidney and serum in the mice with DKD via multi-omics analysis and identified five metabolic pathways related to Empa treatment, which involve in renal purine metabolism (local purinergic signaling); pyrimidine metabolism (mitochondrial function); tryptophan metabolism (reductive stress); nicotinate and nicotinamide metabolism (reductive stress) and serum glycine, serine and threonine metabolism (nephrotoxicity). Our results suggest that EMPA treatment improves renal function and morphology by regulating renal metabolic reprogramming in mice with DKD.
In physiological conditions, the filtered amino acids are almost reabsorbed by tubular epithelial cells. Furthermore, the kidney is a crucial organ that participates in the biosynthesis and catabolism of multiple amino acids and their derivatives . A previous study found that patients with DKD were more predisposed to develop amino acid metabolic derangements than patients with DM or healthy people . In addition, a recent study confirmed that dapagliflozin intervention could reverse the abnormal expression of multiple amino acid transporters and amino acid degrading enzymes in the renal cortex of diabetic mice . In this study, amino acids and their derivatives had the highest proportion of metabolites before and after EMPA treatment, whether in the serum or kidney. Among them, N-Acetylcysteine is the precursor to glutathione in vivo, which has a potent effect of scavenging oxygen free radicals and antioxidants by providing sulfhydryl to increase the synthesis of endogenous glutathione and thus reducing the oxidative damage caused by mitochondrial dysfunction . A recent study confirmed that Canagliflozin could reverse kidney oxidation damage in isoproterenol-induced mice by elevating the concentrations of endogenous glutathione . In addition, the expression of most dipeptides in the kidney was considerably different after EMPA treatment. Particularly, Val-Tyr (downregulation to upregulation, FC = from 0.24 to 6.94), known as Ang-(3–4), is the shortest peptide segment of local angiotensin (Ang) I and Ang II derivatives, which can bind to AT2R to activate Ca2+-ATPase in renal tubular epithelial cells via the cAMP–PKA pathway and antagonist Ang II-mediated Na+ reabsorption and affects the downstream signal for RAAS activity to promote urinary sodium excretion and vasorelaxant effect, thereby declining blood pressure [42,43,44]. In this study, Val-Tyr expression increased remarkably after EMPA intervention. According to previous studies , EMPA can promote sodium excretion and decrease blood pressure. Moreover, in a previous study, EMPA treatment inhibited RAAS overactivation in Apo−/− mice but did not influence the expression of Ang II . Therefore, we speculate that EMPA serves as an antagonist of RAAS by improving the metabolism of some dipeptides or stimulating their synthesis.
Tryptophan and its metabolites are precursors to various microbial biosynthetic products and metabolites of the host . The following three pathways are currently associated with tryptophan catabolism: (1) serotonin (5-hydroxytryptamine) metabolism, (2) the indole pathway and (3) the kynurenine (KYN) pathway . Approximately 90% of 5-HT in the body is metabolized by gut bacteria, with the most typical bacteria being Clostridium sporogenes belonging to the phylum Firmicutes . In addition, Clostridium sporogenes can also utilize tryptophan to produce indoleacetic acid (IAA) and indole-3-propionic acid (IPA) via the indole metabolic pathway . A recent study demonstrated that the proportion of gut bacteria belonging to Firmicutes in db/db mice declined after dapagliflozin intervention , which partly explains the downregulation of intermediates in both serotonin (5-hydroxy-l-tryptophan) and indole (indole-3-acetaldehyde) pathways after EMPA treatment in this study. High expression of 5-HT and its end product metabolite (5-hydroxyindoleacetic acid) in the serum is strongly associated with the pathogenesis of renal dysfunction in DM , whereas 5-HTR antagonists exert renoprotective effects in DKD , suggesting that EMPA plays a renoprotective role by regulating the metabolism of gut bacteria. Most of the tryptophan is degraded via the KYN pathway in the cycle.
Tryptophan, which acts on a series of enzymes, eventually produces acetyl-CoA through the KYN pathway, which completely decomposes to CO2 and water after entering the TCA cycle (Fig. 7). Some intermediate metabolites are produced with this process, especially quinolinic acid, which mediates most of the de novo synthesis of NAD+. NAD+ in the kidney is majorly derived from de novo synthesis of NAD+ mediated by tryptophan metabolism [51, 52]. However, a previous study found that P5P phosphatase (a cofactor for KYNase and a key enzyme regulating the 3-hydroxykynurenine metabolism downstream) is inhibited owing to chronic inflammation in DKD, which blocks tryptophan metabolism to produce NAD+ through the KYN pathway . This finding is consistent with the accumulation of KYN and its metabolites (3-hydroxykynurenine [3-HKYN] and kynurenic acid [KYNA]) observed in chronic kidney disease such as DKD . The accumulation of KYN and its metabolites mediate and enhance oxidative stress, immune activation and inflammatory reaction to exacerbate renal damage . In this study, the expression of acetyl CoA (end product of the KYN pathway) and nicotinic acid mononucleotide (NAMN) (downstream metabolites of quinolinic acid) was upregulated in the kidney of db/db mice after EMPA intervention, indicating that EMPA may reverse renal dysfunction by promoting KYN metabolism to ameliorate NAD+ levels.
NAD+/NADH imbalance in DKD manifests as NAD+ decline and NADH overload, which aggravates oxidative stress mediated by ROC in the kidneys and causes mitochondrial dysfunction, energy metabolic disorders and inactivation of various NAD+-dependent enzymes . In a study, impaired NAD+ synthesis in mice with CKD was associated with decreased expression of key enzymes, including quinolinate phosphoribosyl transferase (QPRT) in the NAD+ de novo synthesis pathway and nicotinamide nucleotide adenylyl transferase 1 (NMNAT1) and NMNAT 3 in the salvage synthesis pathway . In this study, the expression of NAMN (precursor to the NAD+ de novo synthesis pathway) and niacinamide (NAM, precursor to the salvage synthesis pathway) was increased after EMPA treatment, indicating that these key enzymes may act as targets for EMPA to amend energy metabolism. In addition, most NAM is mainly produced by NAD+-consuming enzymes, including the superfamily members of sirtuins (SIRT), poly (ADP-ribose) polymerase (PARP) and CD38 . However, in a study, canagliflozin treatment reversed the reduction in SIRT1 expression in db/db mice and improved local NAD+ metabolism in the kidney . We believe this study will provide novel insights into the renoprotection mechanisms of SGLT2 inhibitors.
As shown in Fig. 7, we discovered a dominant accumulation of metabolites after purine degradation in db/db mice, such as xanthine, hypoxanthine and uric acid. However, metabolites involved in ATP degradation-mediated purine nucleotide conversion pathway, including adenosine diphosphate (ADP), adenosine monophosphate (AMP), adenosine (ADO) and adenine (AD), were downregulated after EMPA treatment. EMPA blocks the catabolism of hypoxanthine but improves the recycling of hypoxanthine, which is consistent with the amelioration of ischemia–reperfusion renal injury in mice after intervention with XOR inhibitors . Previous studies have confirmed that the uric acid-decreasing effect of SGLT2 inhibitors is associated with its promoting effects on uric acid excretion and reabsorption inhibition [58, 59]. Based on the data of this study, we speculate that EMPA may reduce uric acid by blocking purine degradation (similar to XOR inhibitors).
Inhibition of the uric acid generation pathway promotes the increase of ATP degradation-mediated purine nucleotide conversion pathway. In the kidneys, the local paracrine release of purine nucleotides is regulated by connexins . Abnormal expression of connexins and the subsequent interruption of intercellular communication in DKD result in the up-regulated expression and activity of hemichannels, which mediate the imbalances for the release of ATP and ADO and may have contributed in part to the pathology of DKD [61,62,63]. ATP and its degradation metabolite ADO are effective extracellular signaling molecules for purinergic signaling in the kidney, which activate P2 and P1 receptors, respectively, to perform contradictory physiological functions, such as the pro-inflammatory response for ATP and the anti-inflammatory response for ADO . Downregulation of ADO and its receptor A1R is involved in ultrafiltration in the early stages of DKD . A2AR (an adenosine receptor) agonists can attenuate proteinuria and reduce the number of pro-inflammatory cytokines in DKD . Moreover, ADO can be produced by NAD+ metabolism, and enhancing NAD+ metabolism-mediated ADO production has been shown to prevent ischemia-induced acute kidney injury . Researchers have demonstrated that blocking CX43 and its mediated local ATP release from renal tubular epithelial cells can ameliorate renal fibrosis . Interestingly, a recent study observed that EMPA treatment repaired ventricular myocytes' gap junctional intercellular communication and attenuated ventricular fibrosis in mice with metabolic syndrome . Unfortunately, no studies have explored the effects of SGLT2 inhibitors on renal CX43. However, this is undoubtedly an interesting research direction, as the modulation of the local purinergic system by SGLT2 inhibitors may be optimized by using nanotechnology to embed drugs with Cx43 hemichannel blocking effects into nanomaterials with specific targeting to the renal tubules.
A previous study described kidney-specific metabolic reprogramming associated with mitochondrial dysfunctional in db/db mice, manifesting as a compensatory increase in glycolysis and fatty acid metabolism, which was a response to diminished production of ATP induced by dysfunction of the mitochondrial electron transport chain and uncoupling of oxidative phosphorylation . Dihydroorotate dehydrogenase (DHODH) is the first rate-limiting enzyme of pyrimidine de novo synthesis, catalyzing the dihydroorotate oxidized to orotate and further producing various downstream pyrimidine nucleotides . DHODH delivers electrons to ubiquinone during this process and provides reductive ubiquinone for compounds I and III of the respiratory chain, thus coupling the pyrimidine metabolism with mitochondrial phosphorylation . It has been demonstrated that the knockdown of intracellular DHODH partially inhibits the activity of respiratory chain complex III and increases mitochondrial ROS production . However, DHODH inhibitors decline the level of pyrimidine nucleotide but increase the level of upstream metabolites of dihydroorotate . The results of the abovementioned studies are consistent with those observed in this study (the kidney of db/db mice). Furthermore, we found that dihydroorotate levels were significantly increased, whereas the levels of its downstream metabolites were increased, including orotate (despite P = 0.052, it had an increasing trend, with VIP = 1.28 and FC = 1.45), uridine monophosphate (UMP), uracil dinucleotide (UDP) and cytidine monophosphate (CMP). However, the metabolites uridine triphosphate (UTP) and cytidine triphosphate (CTP) involved in the pyrimidine salvage pathway were not detected. Therefore, we speculate that EMPA improves mitochondrial dysfunction and alleviates the metabolic reprogramming of the kidneys in DKD by promoting DHODH-mediated de novo synthesis of pyrimidine to increase mitochondrial electron transport.
Furthermore, in this study, metabolic alterations were observed less in the serum than in the kidney after EMPA treatment. The potential interfered pathway of serum metabolites was ‘glycine, serine and threonine metabolism’, specifically involving the biological transformation of d-serine and l-serine. d-Serine, the most abundant d-amino acid in mammals, is produced from l-serine mediated by serine racemase. The kidney mostly regulates d-serine levels by excreting it through urine or degrading it via d-amino acid oxidase (DAAO, highly expressed in the kidney). However, the intermediate metabolites during d-serine oxidative decomposition, mainly H2O2, have strong renal toxicity. Previous studies have demonstrated that d-serine concentration is positively correlated with GFR and can serve as a clinical diagnostic biomarker for CKD . However, Tomonori Kimura et al. discovered that the level of d-serine did not increase with a decline in GFR in a proportion of patients with CKD . In this study, we found that d-serine levels increased in serum but were reduced in the kidney in db/db mice. Moreover, the ratio of l-serine to d-serine did not change before and after EMPA treatment (Fig. 7), which is consistent with the results of the aforementioned study , collectively indicating that a compensation mechanism in the kidney maintains a certain ratio of l-serine to d-serine. We also found that EMPA significantly reduced the level of local d-serine in the kidney, indicating that it is associated with improving renal function or alterations in l-serine metabolism.
Our study has some limitations. First, the untargeted metabolomics we used is a relative quantification method of the metabolites and needs to be followed up with target validation. Second, the sample size of our study is small, and more animal or clinical samples can be included in the follow-up and combined with some new histological techniques (e.g., spatial omics techniques) for the in-depth study.
In conclusion, this study demonstrated conspicuous metabolic reprogramming in mice with DKD. EMPA treatment improved kidney function and morphology by regulating metabolic reprogramming, including regulation of renal reductive stress, alleviation of mitochondrial dysfunction and reduction in renal oxidative stress reaction. Therefore, this study provides an essential reference for understanding the mechanism of EMPA in renoprotection.
Availability of data and materials
The data supporting the findings of this study are available in Additional files and from the corresponding author.
Saeedi P, Petersohn I, Salpea P, et al. Global and regional diabetes prevalence estimates for 2019 and projections for 2030 and 2045: results from the International Diabetes Federation Diabetes Atlas, 9(Th) Edition. Diabetes Res Clin Pract. 2019;157:107843.
Gheith O, Farouk N, Nampoory N, et al. Diabetic kidney disease: world wide difference of prevalence and risk factors. J Nephropharmacol. 2016;5(1):49–56.
Oellgaard J, Gaede P, Rossing P, et al. Intensified multifactorial intervention in type 2 diabetics with microalbuminuria leads to long-term renal benefits. Kidney Int. 2017;91(4):982–8.
Afkarian M, Zelnick LR, Hall YN, et al. Clinical manifestations of kidney disease among Us adults with diabetes, 1988–2014. JAMA. 2016;316(6):602–10.
Parving HH, de Zeeuw D, Cooper ME, et al. Ace gene polymorphism and losartan treatment in type 2 diabetic patients with nephropathy. J Am Soc Nephrol. 2008;19(4):771–9.
Sakashita M, Tanaka T, Inagi R. Metabolic changes and oxidative stress in diabetic kidney disease. Antioxidants. 2021;10(7):1143.
Hesp AC, Schaub JA, Prasad PV, et al. The role of renal hypoxia in the pathogenesis of diabetic kidney disease: a promising target for newer renoprotective agents including Sglt2 inhibitors? Kidney Int. 2020;98(3):579–89.
Hasegawa S, Tanaka T, Saito T, et al. The oral hypoxia-inducible factor prolyl hydroxylase inhibitor enarodustat counteracts alterations in renal energy metabolism in the early stages of diabetic kidney disease. Kidney Int. 2020;97(5):934–50.
Ahmad AA, Draves SO, Rosca M. Mitochondria in diabetic kidney disease. Cells. 2021;10(11):2945.
Schaffer SW, Jong CJ, Mozaffari M. Role of oxidative stress in diabetes-mediated vascular dysfunction: unifying hypothesis of diabetes revisited. Vascul Pharmacol. 2012;57(5–6):139–49.
Mori Y, Ajay AK, Chang JH, et al. Kim-1 mediates fatty acid uptake by renal tubular cells to promote progressive diabetic kidney disease. Cell Metab. 2021;33(5):1042–61.
Zinman B, Wanner C, Lachin JM, et al. Empagliflozin, cardiovascular outcomes, and mortality in type 2 diabetes. N Engl J Med. 2015;373(22):2117–28.
Mayer GJ, Wanner C, Weir MR, et al. Analysis from the Empa-Reg Outcome((R)) trial indicates empagliflozin may assist in preventing the progression of chronic kidney disease in patients with type 2 diabetes irrespective of medications that alter intrarenal hemodynamics. Kidney Int. 2019;96(2):489–504.
Wanner C, Inzucchi SE, Lachin JM, et al. Empagliflozin and progression of kidney disease in type 2 diabetes. N Engl J Med. 2016;375(4):323–34.
Karamzad N, Eftekhari A, Ashrafi-Asgarabad A, et al. Serum hepcidin, the hepcidin/ferritin ratio and the risk of type 2 diabetes: a systematic review and meta-analysis. Curr Med Chem. 2021;28(6):1224–33.
Ghanim H, Abuaysheh S, Hejna J, et al. Dapagliflozin suppresses hepcidin and increases erythropoiesis. J Clin Endocrinol Metab. 2020;105(4):dgaa057.
Cai T, Ke Q, Fang Y, et al. Sodium-glucose cotransporter 2 inhibition suppresses Hif-1alpha-mediated metabolic switch from lipid oxidation to glycolysis in kidney tubule cells of diabetic mice. Cell Death Dis. 2020;11(5):390.
Tomita I, Kume S, Sugahara S, et al. Sglt2 inhibition mediates protection from diabetic kidney disease by promoting ketone body-induced Mtorc1 inhibition. Cell Metab. 2020;32(3):404–19.
Zhang Z, Xie H, Zuo W, et al. Lysine 2-hydroxyisobutyrylation proteomics reveals protein modification alteration in the actin cytoskeleton pathway of oral squamous cell carcinoma. J Proteomics. 2021;249:104371.
Liang L, Rasmussen MH, Piening B, et al. Metabolic dynamics and prediction of gestational age and time to delivery in pregnant women. Cell. 2020;181(7):1680–92.
Korbut AI, Taskaeva IS, Bgatova NP, et al. sglt2 inhibitor empagliflozin and Dpp4 inhibitor linagliptin reactivate glomerular autophagy in Db/Db Mice, a model of type 2 diabetes. Int J Mol Sci. 2020;21(8):2987.
Klein RL, Hunter SJ, Jenkins AJ, et al. Fibrinogen is a marker for nephropathy and peripheral vascular disease in type 1 diabetes: studies of plasma fibrinogen and fibrinogen gene polymorphism in the Dcct/Edic cohort. Diabetes Care. 2003;26(5):1439–48.
Shlipak MG, Matsushita K, Arnlov J, et al. Cystatin C versus creatinine in determining risk based on kidney function. N Engl J Med. 2013;369(10):932–43.
Inker LA, Schmid CH, Tighiouart H, et al. Estimating glomerular filtration rate from serum creatinine and cystatin C. N Engl J Med. 2012;367(1):20–9.
Zhou Y, Luo Z, Liao C, et al. Mhc Class Ii in renal tubules plays an essential role in renal fibrosis. Cell Mol Immunol. 2021;18(11):2530–40.
Xu S, Lind L, Zhao L, et al. Plasma prolylcarboxypeptidase (Angiotensinase C) is increased in obesity and diabetes mellitus and related to cardiovascular dysfunction. Clin Chem. 2012;58(7):1110–5.
Li X, Jia Z, Shen Y, et al. Coordinate suppression of Sdpr and Fhl1 expression in tumors of the breast, kidney, and prostate. Cancer Sci. 2008;99(7):1326–33.
Regazzetti C, Dumas K, Lacas-Gervais S, et al. Hypoxia inhibits Cavin-1 and Cavin-2 expression and down-regulates caveolae in adipocytes. Endocrinology. 2015;156(3):789–801.
Almomen SM, Guan Q, Liang P, et al. Daily intake of grape powder prevents the progression of kidney disease in obese type 2 diabetic Zsf1 rats. Nutrients. 2017;9(4):345.
Tziastoudi M, Stefanidis I, Zintzaras E. The genetic map of diabetic nephropathy: evidence from a systematic review and meta-analysis of genetic association studies. Clin Kidney J. 2020;13(5):768–81.
Lim RR, Vaidya T, Gadde SG, et al. Correlation between systemic S100a8 and S100a9 levels and severity of diabetic retinopathy in patients with type 2 diabetes mellitus. Diabetes Metab Syndr. 2019;13(2):1581–9.
Monaghan MT, Bailey MA, Unwin RJ. Purinergic signalling in the kidney: in physiology and disease. Biochem Pharmacol. 2021;187:114389.
Rushworth GF, Megson IL. Existing and potential therapeutic uses for N-acetylcysteine: the need for conversion to intracellular glutathione for antioxidant benefits. Pharmacol Ther. 2014;141(2):150–9.
Barba C, Benoit B, Bres E, et al. A low aromatic amino-acid diet improves renal function and prevent kidney fibrosis in mice with chronic kidney disease. Sci Rep. 2021;11(1):19184.
Liu X, Xu X, Shang R, et al. Asymmetric Dimethylarginine (Adma) as an important risk factor for the increased cardiovascular diseases and heart failure in chronic kidney disease. Nitric Oxide. 2018;78:113–20.
Kimura T, Hesaka A, Isaka Y. D-Amino acids and kidney diseases. Clin Exp Nephrol. 2020;24(5):404–10.
Makrides V, Camargo SM, Verrey F. Transport of amino acids in the kidney. Compr Physiol. 2014;4(1):367–403.
Zhou C, Zhang Q, Lu L, et al. Metabolomic profiling of amino acids in human plasma distinguishes diabetic kidney disease from type 2 diabetes mellitus. Front Med. 2021;8:765873.
Kogot-Levin A, Hinden L, Riahi Y, et al. Proximal tubule Mtorc1 is a central player in the pathophysiology of diabetic nephropathy and its correction by Sglt2 inhibitors. Cell Rep. 2020;32(4):107954.
Fard JK, Hamzeiy H, Sattari M, et al. Triazole rizatriptan induces liver toxicity through lysosomal/mitochondrial dysfunction. Drug Res. 2016;66(9):470–8.
Hasan R, Lasker S, Hasan A, et al. Canagliflozin ameliorates renal oxidative stress and inflammation by stimulating Ampk-Akt-Enos pathway in the isoprenaline-induced oxidative stress model. Sci Rep. 2020;10(1):14659.
Axelband F, Dias J, Miranda F, et al. Angiotensin-(3–4) counteracts the angiotensin Ii inhibitory action on renal Ca2+-Atpase through a Camp/Pka pathway. Regul Pept. 2012;177(1–3):27–34.
Luzes R, Crisostomo T, Silva PA, et al. Angiotensin-(3–4) normalizes blood pressure, decreases Na(+) and energy intake, but preserves urinary Na(+) excretion in overweight hypertensive rats. Biochim Biophys Acta Mol Basis Dis. 2021;1867(3):166012.
Shaltout HA, Westwood BM, Averill DB, et al. Angiotensin metabolism in renal proximal tubules, urine, and serum of sheep: evidence for Ace2-dependent processing of angiotensin Ii. Am J Physiol Renal Physiol. 2007;292(1):F82-91.
Liu Y, Wu M, Xu B, et al. Empagliflozin alleviates atherosclerosis progression by inhibiting inflammation and sympathetic activity in a normoglycemic mouse model. J Inflamm Res. 2021;14:2277–87.
Agus A, Planchais J, Sokol H. Gut microbiota regulation of tryptophan metabolism in health and disease. Cell Host Microbe. 2018;23(6):716–24.
Kwiatkowska I, Hermanowicz JM, Mysliwiec M, et al. Oxidative storm induced by tryptophan metabolites: missing link between atherosclerosis and chronic kidney disease. Oxid Med Cell Longev. 2020;2020:6656033.
Lee DM, Battson ML, Jarrell DK, et al. Sglt2 inhibition via dapagliflozin improves generalized vascular dysfunction and alters the gut microbiota in type 2 diabetic mice. Cardiovasc Diabetol. 2018;17(1):62.
Saito J, Suzuki E, Tajima Y, et al. Increased plasma serotonin metabolite 5-hydroxyindole acetic acid concentrations are associated with impaired systolic and late diastolic forward flows during cardiac cycle and elevated resistive index at popliteal artery and renal insufficiency in type 2 diabetic patients with microalbuminuria. Endocr J. 2016;63(1):69–76.
Yang Y, Huang H, Xu Z, et al. Serotonin and its receptor as a new antioxidant therapeutic target for diabetic kidney disease. J Diabetes Res. 2017;2017:7680576.
Xu J, Kitada M, Koya D. Nad(+) homeostasis in diabetic kidney disease. Front Med. 2021;8:703076.
Liu L, Su X, Quinn WJ 3rd, et al. Quantitative analysis of nad synthesis-breakdown fluxes. Cell Metab. 2018;27(5):1067–80.
Oxenkrug GF. Increased plasma levels of xanthurenic and kynurenic acids in type 2 diabetes. Mol Neurobiol. 2015;52(2):805–10.
Liu X, Luo D, Huang S, et al. Impaired nicotinamide adenine dinucleotide biosynthesis in the kidney of chronic kidney disease. Front Physiol. 2021;12:723690.
Liu JR, Deng ZH, Zhu XJ, et al. Roles of nicotinamide N-methyltransferase in obesity and type 2 diabetes. Biomed Res Int. 2021;2021:9924314.
Umino H, Hasegawa K, Minakuchi H, et al. High basolateral glucose increases sodium-glucose cotransporter 2 and reduces sirtuin-1 in renal tubules through glucose transporter-2 detection. Sci Rep. 2018;8(1):6791.
Tani T, Okamoto K, Fujiwara M, et al. Metabolomics analysis elucidates unique influences on purine/pyrimidine metabolism by xanthine oxidoreductase inhibitors in a rat model of renal ischemia-reperfusion injury. Mol Med. 2019;25(1):40.
Cowie MR, Fisher M. Sglt2 inhibitors: mechanisms of cardiovascular benefit beyond glycaemic control. Nat Rev Cardiol. 2020;17(12):761–72.
Novikov A, Fu Y, Huang W, et al. Sglt2 inhibition and renal urate excretion: role of luminal glucose, Glut9, and Urat1. Am J Physiol Renal Physiol. 2019;316(1):F173–85.
Hills CE, Price GW, Squire PE. Mind the gap: connexins and cell-cell communication in the diabetic kidney. Diabetologia. 2015;58(2):233–41.
Peleli M, Carlstrom M. Adenosine signaling in diabetes mellitus and associated cardiovascular and renal complications. Mol Aspects Med. 2017;55:62–74.
Dwyer KM, Kishore BK, Robson SC. Conversion of extracellular Atp into adenosine: a master switch in renal health and disease. Nat Rev Nephrol. 2020;16(9):509–24.
Eftekhari A, Vahed SZ, Kavetskyy T, et al. Cell junction proteins: crossing the glomerular filtration barrier in diabetic nephropathy. Int J Biol Macromol. 2020;148:475–82.
Bell TD, Welch WJ. Regulation of renal arteriolar tone by adenosine: novel role for type 2 receptors. Kidney Int. 2009;75(8):769–71.
Persson P, Friederich-Persson M, Fasching A, et al. Adenosine A2 a receptor stimulation prevents proteinuria in diabetic rats by promoting an anti-inflammatory phenotype without affecting oxidative stress. Acta Physiol. 2015;214(3):311–8.
Ralto KM, Rhee EP, Parikh SM. Nad(+) homeostasis in renal health and disease. Nat Rev Nephrol. 2020;16(2):99–111.
Xu H, Wang M, Li Y, et al. Blocking connexin 43 and its promotion of atp release from renal tubular epithelial cells ameliorates renal fibrosis. Cell Death Dis. 2022;13(5):511.
Jhuo SJ, Liu IH, Tasi WC, et al. Characteristics of ventricular electrophysiological substrates in metabolic mice treated with empagliflozin. Int J Mol Sci. 2021;22(11):6105.
Sas KM, Kayampilly P, Byun J, et al. Tissue-specific metabolic reprogramming drives nutrient flux in diabetic complications. JCI Insight. 2016;1(15):e86976.
Zhang J, Teran G, Popa M, et al. Dhodh inhibition modulates glucose metabolism and circulating Gdf15, and improves metabolic balance. iScience. 2021;24(5):102494.
Khutornenko AA, Roudko VV, Chernyak BV, et al. Pyrimidine biosynthesis links mitochondrial respiration to the P53 pathway. Proc Natl Acad Sci USA. 2010;107(29):12828–33.
Fang J, Uchiumi T, Yagi M, et al. Dihydro-orotate dehydrogenase is physically associated with the respiratory complex and its loss leads to mitochondrial dysfunction. Biosci Rep. 2013;33(2):e00021.
Bennett CF, O’Malley KE, Perry EA, et al. Peroxisomal-derived ether phospholipids link nucleotides to respirasome assembly. Nat Chem Biol. 2021;17(6):703–10.
Kimura T, Hamase K, Miyoshi Y, et al. Chiral amino acid metabolomics for novel biomarker screening in the prognosis of chronic kidney disease. Sci Rep. 2016;6:26137.
Hesaka A, Sakai S, Hamase K, et al. D-Serine reflects kidney function and diseases. Sci Rep. 2019;9(1):5104.
This work was supported by the National Natural Science Foundation of China (Grant No. 8210032527), GuangDong Basic and Applied Basic Research Foundation (Grant No. 2020A1515111209), Young Innovative Talents Project of General Colleges and Universities in Guangdong Province (Grant No. 2018KQNCX010), Postdoctoral Fund the Seventh Affiliated Hospital, Sun Yat-sen University (Grant No. ZSQYRSFPD0030), and General program of China Postdoctoral Science Foundation (Grant No. 2021M703749).
Consent for publication
All the authors have read and approved the paper and declare no potential conflicts of interest in the paper. All the authors agree to publish this paper.
The authors have declared that no conflict of interest exists.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1.
Quality control of kidney and serum proteomics. Figure S2. Quality control of kidney and serum metabolomics. Figure S3. The identification of differential expressed metabolites between db/db and Empa group.
Additional file 2: Table S1.
Information of 6009 quantified proteins in renal proteomics. Table S2. Information of 1481 quantified proteins in serum proteomics. Table S3. Information of 501 and 91 DEPs identified in db/db and Empa groups in renal proteomics, respectively. Table S4. Information of 32 EAPs between Empa and db/db groups in kidney. Table S5. Information of 235 and 83 DEPs identified in db/db and Empa groups in serum proteomics, respectively. Table S6. Information of 51 EAPs between Empa and db/db groups in serums. Table S7. Information of 897 quantified metabolites in renal metabolomics. Table S8. Information of 678 quantified metabolites in serum metabolomics. Table S9. Information of 264 identified differentially expressed metabolites between db/db and db/m groups in kidney metabolomics. Table S10. Information of 174 identified differentially expressed metabolites between db/db and db/m groups in serum metabolomics. Table S11. Information of 94 Empa-altered metabolites between Empa and db/db groups in kidney. Table S12. Information of 37 Empa-altered metabolites between Empa and db/db groups in serum.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Lu, YP., Zhang, ZY., Wu, HW. et al. SGLT2 inhibitors improve kidney function and morphology by regulating renal metabolic reprogramming in mice with diabetic kidney disease. J Transl Med 20, 420 (2022). https://doi.org/10.1186/s12967-022-03629-8