Glucagon-Like Peptide 1 Receptor Agonist Improves Renal Tubular Damage in Mice with Diabetic Kidney Disease

Abstract

Purpose

This study aims to investigate the renal protective effect of glucagon-like peptide 1 receptor agonist (GLP-1RA) on improving renal tubular damage in diabetic kidney disease (DKD) and to explore the potential mechanism of GLP-1RA on renal tubular protection.

Methods

Long-acting GLP-1RA was used to treat DKD mice for 12 weeks. The label-free quantitative proteomic analysis of renal proteins was conducted to explore the differentially expressed proteins (DEPs) in the renal tissues of the control, DKD and GLP-1RA groups. The DEPs and markers of renal tubular injury were verified by qPCR in vivo and in vitro. The expression of glucagon-likepeptide-1 receptor (GLP-1R) in renal tubules was determined by immunofluorescence staining.

Results

GLP-1RA treatment significantly improved the tubular damages in kidney tissues of DKD mice and mTEC cells stimulated by high glucose (HG). Proteomics analysis revealed that 30 proteins in kidney tissue were differentially expressed among three groups. Seminal vesicle secretory protein 6 (SVS6) was the most differentially expressed protein in kidney tissues among three groups of mice. The expression changes of Svs6 mRNA in vitro and in vivo detected by qPCR were consistent with the results of proteomic analysis. Furthermore, reduction of Svs6 expression by SVS6 siRNA could attenuate HG-stimulated tubular injury in mTEC cells. Immunofluorescence staining also found that GLP-1R was widely expressed in renal tubules in vitro and in vivo.

Conclusion

GLP-1RA significantly improved renal tubular damage in DKD mice. SVS6 may be a potential therapeutic target for GLP-1RA in the treatment of DKD.

Keywords:

• diabetic kidney disease
• glucagon-like peptide-1 receptor agonist
• renal tubular injury
• seminal vesicle secretory protein 6

Introduction

Diabetes Mellitus (DM) is recognized as a metabolic disorder characterized by hyperglycemia, which is caused by absolute or relative deficiency of insulin, and can affect people at different life stages.¹ Diabetic kidney disease (DKD) is the most common microvascular complication of DM, and the main cause of chronic kidney disease and end-stage renal disease worldwide.^2,3 The underlying pathogenesis of DKD is complex and multifactorial, including glycolipid metabolism disorder, oxidative stress, inflammatory response, autophagy and genetic factors.⁴ The clinical diagnostic criteria of DKD are defined as persistent proteinuria and decreased glomerular filtration rate (eGFR).^5,6 These alterations are associated with the production and composition of the extracellular matrix, leading to thickening of the glomerular basement membrane and expansion of the mesangial matrix, ultimately developing into glomerular hypertrophy and glomerulosclerosis.^4,6

Previous studies have mainly focused on glomerular lesions in the progression of DKD. However, increasing evidence has changed our traditional understanding of DKD, that is, in addition to glomerular damage, tubular interstitial damage is also closely related to the development of DKD.^7–9 Under the stimulation of high glucose, a series of pathological changes such as renal tubular and interstitial thickening, renal tubular atrophy, interstitial inflammation and fibrosis may occur in the proximal renal tubules, which accelerate renal failure.^10,11 A series of clinical and experimental studies have concluded that renal proximal tubules play a major role in the development of early initial albuminuria in DKD.^7,12 Therefore, in the process of inhibiting the occurrence and development of DKD, diabetic tubulopathy must be concerned as much as diabetic glomerulopathy.

Glucagon-like peptide 1 (GLP-1) is mainly secreted by L cells in the intestine, which stimulates insulin secretion and inhibits glucagon secretion in response to food intake.¹³ Other effects of GLP-1 include inhibiting food intake, retarding gastric emptying, limiting weight gain, lowering blood pressure, and reducing postprandial triglyceride and free fatty acid concentrations.¹⁴ Glucagon-like peptide-1 receptor agonist (GLP-1RA) is a novel type of anti-diabetic drug with the potential to prevent and treat DKD.¹³ GLP-1RAs may play an indirect renal protective effect by improving the traditional risk factors for DKD, such as blood glucose control, weight loss and blood pressure reduction.¹⁵ Moreover, it has been reported that glucagon-like peptide-1 receptor (GLP-1R) is widely expressed in the glomeruli and proximal tubules.^16,17 GLP-1RAs may also have a direct effect on DKD, including proximal tubular natriuretic stimulation, cAMP/PKA signaling regulation, and renin angiotensin system inhibition.¹⁵ Exenatide-loaded microsphere, a long-acting GLP-1RA, has a good regulatory effect on blood glucose and body weight by prolonging the stimulation of GLP-1R.^18,19 Accumulating data have indicated that Exendin-4 exerts a beneficial effect on DKD by acting on GLP-1R.^20,21 Recent clinical trials have demonstrated that Exenatide-loaded microspheres reduced albuminuria by 26% in patients with type 2 diabetes mellitus (T2DM) compared with non-GLP-1RA group, independent of their blood glucose-lowering effect.²² Most of previous studies have mainly investigated the glomerular damage in DM models and the mitigation effect of GLP-1RA on the glomerular injury, often neglecting its impact on renal tubular injury.^23–25 There are very few literatures suggesting that GLP-1RA also has a therapeutic effect on renal tubular injury. For example, the administration of Exendin-4 significantly alleviated renal tubular damage in DM mice induced by streptozotocin (STZ), which was shown to reduce the infiltration of macrophages and inhibit the production of chemokines in renal tubular cells.²¹

However, so far, the research on the therapeutic effect of GLP-1RA on renal tubular damage of DKD is still very limited. At the same time, the molecular mechanism of the protective effect of GLP-1RA on the renal tubules of DKD is still inconclusive. Therefore, we evaluated the therapeutic effect of GLP-1RA on renal tubular injury in DKD mice, and investigated the potential effects of GLP-1RA on the proteomics of the kidney in DKD mice, providing a favorable clue for finding possible therapeutic targets for DKD.

Materials and Methods

Animal

Eight-week-old male C57BL/6J mice (n=30) were obtained from Model Animal Research Center of Nanjing University and housed in Laboratory Animal Center at Tongji University Hubei Campus. All mice were maintained under stable conditions in the specific pathogen-free (SPF) environment (temperature: 21 ± 4^◦C, humidity: 55 ± 20%, 12/12 light-dark cycle). After 1 week of adaptation, mice were randomly divided into 3 groups: control group (n=10), DKD group (n=10) and GLP-1RA treated group (n=10, GLP-1RA group).

Mice in control group were fed with a normal chow diet containing 10% Kcal fat (Puluteng Biotech, China). Mice in DKD and GLP-1RA group were fed with a high fat diet (HFD) containing 42% Kcal fat (Puluteng Biotech, China). After 4 weeks, all HFD-fed mice were intraperitoneally injected with STZ (100mg/kg; Sigma Aldrich, USA) freshly diluted in sodium citrate buffer (0.1M, pH 4.5; Solarbio, China); mice of control mice were intraperitoneally injected with an equal amount of sodium citrate buffer. One week after STZ injection, fasting blood glucose level ≥11.1mmol/L was defined as successful models of T2DM. Mice with T2DM were fed with HFD for another 6 weeks to construct DKD models according to the previous study.²⁶ Then, mice in GLP-1RA group were intraperitoneally injected Exenatide-loaded microspheres once a week (10mg/kg/w; Eli Lilly, USA) for 12 weeks; mice in control group and DKD group were intraperitoneally injected with an equal amount of saline for 12 weeks (Figure 1A). During this period, mice in DKD group and GLP-1RA group continued to be fed with HFD, and fasting blood-glucose and body weight were monitored weekly. Immediately after the above 12 weeks, the three groups of mice were placed in metabolic cages and urine was collected for 24 hours. Next, the blood and kidney samples were collected and stored at − 80 °C. All animal experiments were approved by the Animal Ethics Committee of Tongji Hospital, Tongji University School of Medicine (2021-DW-(003)). This study was conducted in accordance with the Guide for the Care and Use of Laboratory Animals (NIH Publication, 8th Edition, 2011), and Animal Research: Reporting of In Vivo Experiments Guidelines. All international, national, and/or institutional guidelines applicable to the care and use of animals were followed.

Figure 1 GLP-1RA reduces body weight, lowers blood glucose, and improves kidney function. (A) Schematic presentation of the experimental design and protocol. (B) The gross appearance of the kidney among three groups. (C) Body weight, kidney weight and blood glucose of mice among three groups (n = 10). (D) Comparison of body weight and blood glucose of mice among three groups in different periods (n = 10). (E) Scr and BUN of mice among three groups (n = 10). Data are expressed as mean ± SD, *P<0.05, **P<0.01 vs control group; ^#P<0.05, ^##P<0.01 vs DKD group.

Abbreviations: GLP-1RA, glucagon-like peptide-1 receptor agonist; DKD, diabetic kidney disease; Scr, serum creatinine; BUN, blood uric nitrogen.

Cell Culture

Murine kidney proximal tubular epithelial (mTEC) cells, a gift from Department of Nephrology, Tongji Hospital were cultured in Dulbecco’s modified Eagle’s medium (DMEM) (Gibco, USA) with 10% FBS (Gibco, USA) and 1% penicillin/streptomycin (KeyGEN, China). Cells were incubated in an incubator with 5% CO2 at 37°C, and passaged every 3–5 days. The culture medium was replaced every 2 days. When mTEC cells were seeded at ~60% confluence in 6-well plates, they were harvested overnight and subsequently exposed to 5.5mM glucose (NG group), 30 mM glucose (HG group) and 30 mM glucose with 10nM Exenatide (GLP-1RA group) for an additional 48h.²⁰

For cell transfection, specific SVS6 siRNA, and negative control siRNA (Sangon Biotech, China) were transfected into mTEC cells at 20μM concentration in 6-well plates for 48h by using Lipofectamine™ 3000 reagent (Invitrogen, USA) according to the manufacturer’s protocol. mTEC cells were divided into the following groups: 1) control group (5.5mM glucose + negative control siRNA); 2) HG group (30 mM glucose + negative control siRNA); 3) GLP-1RA group (30 mM glucose with 10nM Exenatide + negative control siRNA); 4) HG+si-SVS6 (30 mM glucose+ specific SVS6 siRNA).

Biochemical Measurements

After fasting for at least 12h, blood samples were collected from all mice. Serum creatinine (Scr) and blood uric nitrogen (BUN) were measured at the central clinical laboratory of Tongji Hospital using fully automatic biochemical analyzer (PUZS-300X, Perlong, China). The levels of urine KIM-1 were measured using an enzyme-linked immunosorbent assay (ELISA) kit (CUSABIO, USA) according to the manufacturer’s instructions.

Histological Analysis

Kidney tissues were embedded in paraffin and sectioned for staining. Hematoxylin and Eosin (H&E), Periodic Acid-Silver Methenamine (PASM) and periodic acid-Schiff (PAS) staining were used to observe the morphological changes of kidney. To examine the effect of GLP-1RA on glomerular cross-sectional area, mesangial area and tubular injury area, PASM and PAS-stained kidney sections were analyzed respectively. Six mice were randomly selected from each group, and 20 glomeruli were randomly selected from each mouse for mesangial expansion assessment. Twenty cortico-medullary areas in each section were randomly assessed for tubulointerstitial lesions on a scale of 0 to 3. 0, normal; 1, less than 25% of the tubular area is injured; 2, 26–50% of the tubular area is injured; and 3, >50% of the tubular area is injured.²⁷ The evaluation of the mesangial area percentage (%) and the analysis of the area percentage (%) of the damaged renal tubular area were performed using Image-Pro Plus 6.0 software (Media Cybernetics, USA). The kidney sections were also stained with Sirius Red to assess the level of renal fibrosis. A semi-quantitative analysis of the collagen area on the Sirius Red stained sections was performed by Image-Pro Plus 6.0 software (Media Cybernetics, USA).

Immunofluorescence

Immunofluorescence staining technique was performed as previously described.²⁸ Kidney sections, ovary sections and mTEC cells were stained with GLP-1R antibody (NBP1-97308, Novus Biologicals, USA) overnight at 4°C, followed by incubation with Alexa-488 goat anti-rabbit (#111-545-003, Jackson ImmunoResearch, USA). Then, the sections were counterstained with DAPI-staining-solution (Boster, USA), and the images were captured using inverted optical microscope (Nikon, Japan).

Label-Free Quantification (LFQ)-Based Proteomics Analysis

Protein Extraction and Trypsin Digestion

Kidney tissues of three mice in each group were lysed to extract protein, and bicinchoninic acid (BCA) protein assay kit (Beyotime, China) was used for protein quantification according to manufacturer’s instructions. About 200μg protein of each sample was trypsinized using filter aided proteome preparation (FASP) method.²⁹ The peptide was desalted on C18 cartridge, concentrated by vacuum centrifugation, re-dissolved with 40μL, 0.1% formic acid solution, and quantified (OD280).

LC-MS/MS Analysis

The enzymatic hydrolysates were analyzed by LC-MS/MS analysis, as previously reported.³⁰ Each sample was separated by high-performance liquid chromatography (HPLC) in an Easy nLC system (Proxeon Biosystems). The chromatographic column was balanced with 95% liquid A (0.1% formic acid). The samples were automatically transferred to the loading column (Thermo Scientific Acclaim PepMap100, 100μm*2cm, nanoViper C18), separated by the analytical column (Thermo scientific EASY column, 10cm, ID75μm, 3μm, C18-A2) at a constant flow rate of 300 nL/min, and eluted with liquid B (0.1% formic acid in 84% acetonitrile) in different gradients. The eluted peptides were introduced into Q-Exactive mass spectrometer (MS) instrument (Thermo Fisher Scientific, USA). The instrument parameters were set as follows: the MS1 resolution ratio, 70,000 at 200 m/z; automatic gain control target (AGC), 1e6; maximum injection time (IT), 50ms; dynamic exclusion, 60s. The mass-to-charge ratio of polypeptides and fragments was collected as follows: 20 pieces of atlas were collected after each full scan (MS2 scan), MS2 activation type, high energy collisional dissociation (HCD); Isolation window, 2 m/z; The MS2 resolution ratio, 17,500 at 200 m/z; Normalized collision energy, 30eV; Underfill, 0.1%.

Bioinformatics Analysis

The original data collected by LC-MS/MS were RAW files, and MaxQuant software (version 1.5.3.17) was used for database identification and data quantitative analysis.³¹ The mass tolerance for MS1 was 10 ppm. The cut-off of false discovery rate (FDR) was 1% for both proteins and peptides. The quantitative information of the target protein dataset was normalized (normalized to the interval of (−1,1)) in order to assess the homogeneity of the designated groups. The data were analyzed with at least three valid values in each group. The protein expression level increased by 1 fold or decreased by 1 fold, indicating the difference among groups was significant (P<0.05). To generate a hierarchical clustering heatmap, the Complex heatmap R package (R Version 3.4) was used to classify the two dimensions of the sample and protein expression simultaneously (distance algorithm: Euclid; linkage method: average linkage). Venn diagrams were generated using a more flexible tool named InteractiVenn, which is a web-based tool for the analysis of sets. Gene ontology (GO, https://www.blast2go.com/) was used to determine the function of proteins. GO annotations were classified into the following three categories: biological process (BP), molecular function (MF), and cell population points (CC). Pathway analysis of proteins was conducted through the Kyoto Encyclopedia of Genes and Genomes (KEGG, https://www.kegg.jp/). Fisher’s exact test was used to compare the distribution of GO classification and KEGG pathway in the target protein dataset and the total protein dataset. Enrichment analysis was performed on GO annotations and KEGG pathway annotations in the target protein dataset.

Quantitative Real-Time PCR (qPCR)

Total RNA was extracted from kidney tissues and mTEC cells using a TRIzol reagent (Invitrogen, USA). RNA quality was measured using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, USA). Approximately 500 ng RNA was converted into 10ul cDNA using a reverse transcription kit (Takara RR037A, Japan), according to the manufacturer’s protocol. QPCR reaction systems were prepared with SYBR® Premix Ex Taq^TM (TAKARA RR820A, Japan) using a QuantStudio7Flex System (ABI, USA). The cDNA was diluted 10 times with RNase-free ddH2O.

The reaction mixture was configured according to the following components: TB Green^TM Premix Ex Taq II 5μL; PCR Forward Primer (10μM) 0.2μL; PCR Reverse Primer (10μM) 0.2μL; ROX Dye II 0.2μL; cDNA 2uL; ddH2O 2.4uL. The amplification conditions were set as follows: 95°C for 30 sec, followed by 40 cycles of 95°C for 5 sec, 60°C for 30 sec, and 72°C for 30 sec. Primer sequences were listed in Table 1. Expression levels of mRNA were quantified by the 2^−ΔΔCt method. Gapdh was used as an internal control for mRNAs.

Table 1 Gene-Specific Primers for qPCR Analysis

Gene Name	Forward (5ʹ-3ʹ)	Reverse (5ʹ-3ʹ)
Gapdh	AGGTCGGTGTGAACGGATTTG	TGTAGACCATGTAGTTGAGGTCA
Fabp-1	CCAGAAAGGGAAGGACATCA	GTCTCCAGTTCGCACTCCTC
Kim-1	ATCGTGGAATCACAACGACA	TGCCCCTTTAAGTTGTACCG
Ngal	TGGAAGAACCAAGGAGCTGT	GGTGGGGACAGAGAAGATGA
Sglt2	GTAGCACGCTCTTCACCATG	CACCCAGAGCCTTCCAACTA
Col1a1	ATGTGCCACTCTGACTGGAA	TCCATCGGTCATGCTCTCTC
Col3a1	CCCAACCCAGAGATCCCATT	GGTCACCATTTCTCCCAGGA
Fn1	TACCCTTCCACACCCCAATC	TGCCAGGAAGCTGAATACCA
Svs6	GGTGACAGAAACAGCTGCAA	AACCTCCTCATGGACAATGC
Cpt2	TCCTCGATCAAGATGGGAAC	GATCCTTCATCGGGAAGTCA
Apoa1	TATGTGGATGCGGTCAAAGA	ACGGTTGAACCCAGAGTGTC
Gstt1	TCCCTGTTTTCCTTGGTGAG	ATGCATCAGCTCTGTGATGG
Mup20	CTGCTGTGTTTGGGACTGAC	CTTTTGTCAGTGGCCAGCAT
Nampt	AGATACTGTGGCGGGAATTG	GCTGCTGGAACAGAATAGCC
Svs6 siRNA	AUUGUCCAUGAGGAGGUUUAUTT	AUAAACCUCCUCAUGGACAAUTT
Negative control siRNA	UUCUCCGAACGUGUCACGUTT	ACGUGACACGUUCGGAGAATT

Abbreviations: GAPDH, glyceraldehyde-3-phosphate dehydrogenase; FABP-1, fatty acid binding protein 1; KIM-1, kidney injury molecule-1; NGAL, neutrophil gelatinase-associated lipocalin; SGLT2, sodium-glucose co-transporter 2; COL1A1, collagen type I alpha 1 chain; COL3A1, collagen type III alpha 1 chain; FN1, fibronectin 1; SVS6, seminal vesicle secretory protein 6; CPT2, carnitine palmitoyl transferase 2; APOA1, apolipoprotein A1; GSTT1, glutathione S-transferase theta 1; MUP20, major urinary protein 20; NAMPT, nicotinamide phosphoribosyl transferase; SVS6, seminal vesicle secretory protein 6.

Statistical Analysis

All statistical analyses were performed with SPSS software version 22.0 (IBM, USA). All data were presented as mean ± standard deviation (SD). One-way analysis of variance (ANOVA) followed by LSD’s post hoc test or Dunnett’s T3 post hoc multiple comparison test was performed to compare the significant differences among three groups. P-values < 0.05 was set as statistical significance.

Results

GLP-1RA Reduces Body Weight, Lowers Blood Glucose, and Improves Kidney Function

The gross appearance of kidney is shown in Figure 1B. The body weight, the kidney weight and blood glucose were all significantly increased in DKD mice compared with controls, while they were all dramatically reduced after GLP-1RA treatment (Figure 1C and D). In addition, Scr and BUN were obviously elevated in DKD group compared to the control group, whereas they were significantly decreased after GLP-1RA treatment (Figure 1E).

GLP-1RA Improves Renal Tubular Injury in DKD Mice and mTEC Cells Stimulated by High Glucose

Glomerular mesangial damage and renal tubular injury were observed under the microscope by H&E, PASM and PAS staining. Glomerular mesangial expansion, disordered arrangement of renal tubular epithelial cells and vacuolar degeneration were significant in DKD mice compared with the control mice. In contrast, GLP-1RA treatment significantly inhibited mesangial expansion and improved tubular injury in DKD mice (Figure 2A–F).

Figure 2 GLP-1RA improves renal tubular injury in DKD mice. (A–C) H&E staining (A), PASM staining (B) and PAS staining (C) of kidney sections among three groups (×200) (n=6). (D–F) Relative quantification of the mesangial expansion (D), renal tubular injury area (%) (E) and renal tubular injury score (F) in the kidney among three groups. (G) The mRNA expression of Fabp1, Kim-1, Ngal and Sglt2 in kidney tissues among three groups (n=10). (H) Urine KIM-1 levels among three groups (n=10). (I) The mRNA expression of Fabp1, Kim-1, Ngal and Sglt2 in mTEC cells among three groups (n=10). (J) Sirius Red staining of kidney sections among three groups (×200) (n=3). (K) Collagen fiber area (%) among three groups (n=6). (L) The mRNA expression of Col1a1, Col3a1 and Fn1 in kidney tissues among three groups (n=10). Data are expressed as mean ± SD, *P<0.05, **P<0.01 vs control group; ^#P<0.05, ^##P<0.01 vs DKD group.

Abbreviations: GLP-1RA, glucagon-like peptide-1 receptor agonist; DKD, diabetic kidney disease; H&E, hematoxylin and eosin; PAS, periodic acid-Schiff; FABP1, fatty acid binding protein 1; KIM-1, kidney injury molecule-1; NGAL, neutrophil gelatinase-associated lipocalin; SGLT2, sodium-glucose co-transporter 2; COL1A1, collagen type I alpha 1 chain; COL3A1, collagen type III alpha 1 chain; FN1, fibronectin 1.

To further explore the therapeutic effect of GLP-1RA on renal tubular damage, we also detected the mRNA expression of renal tubular damage markers in DKD mice, such as fatty acid binding protein 1 (FABP1), kidney injury molecule-1 (KIM-1), neutrophil gelatinase-associated lipocalin (NGAL) and sodium-glucose co-transporter 2 (SGLT2). QPCR showed that the mRNA levels of Fabp1, Kim-1, Ngal and Sglt2 in renal tissues of DKD mice were significantly higher than those of control mice, and GLP-1RA treatment dramatically reduced the expression of Fabp1, Kim-1 and Sglt2 in DKD mice (Figure 2G). Meanwhile, the changes of KIM-1 levels in the urine of the three groups of mice were consistent with the gene expression changes of Kim-1 in the kidney tissues of mice among the three groups (Figure 2H). In addition, the mRNA levels of Fabp1, Kim-1, Ngal and Sglt2 in mTEC cells stimulated by high glucose (HG) were significantly elevated compared with the control group, whereas they were significantly reduced after GLP-1RA treatment (Figure 2I). Our results indicated that GLP-1RA can significantly improve diabetic renal tubular damage in vivo and in vitro.

Furthermore, the analysis of the results from Sirius red staining indicated that the renal collagen deposition was dramatically increased in DKD group compared with controls, whereas it was significantly decreased after GLP-1RA treatment (Figure 2J and K). The mRNA expressions of collagen type I alpha 1 chain (Col1a1), collagen type III alpha 1 chain (Col3a1) and fibronectin 1 (Fn1) were significantly up-regulated in renal tissues of the DKD group compared to the control group, whereas they were significantly down-regulated after GLP-1RA treatment, indicating that GLP-1RA treatment significantly ameliorates renal collagen deposition in DKD mice (Figure 2L).

Identification of Significant Differentially Expressed Proteins (DEPs) in Kidney of Mice Among Three Groups by LFQ-Based Proteomic Analysis

To explore the therapeutic mechanism of GLP-1RA on DKD mice, LFQ proteomics analysis was performed to detect the differentially expressed proteins (DEPs) in the kidneys of mice from the control, DKD and GLP-1RA groups. Protein identification indicated that a total of 3270 common proteins were observed in the kidneys of three groups (Figure 3A). To investigate the potential target proteins of GLP-1RA in the treatment of DKD, fold change (FC) > 1 (up-regulation> 1 or down-regulation < 1) and P value <0.05 were applied to define DEPs among the control, DKD and GLP-1RA groups. A heatmap was presented to show that 30 proteins significantly up-regulated in the kidneys of DKD group compared with the control group, while these proteins were significantly down-regulated in the GLP-RA group (Figure 3B). Seminal vesicle secretory protein 6 (SVS6) was the most significant differentially expressed protein among the three groups. In the results of proteomics analysis, the protein expression levels of SVS6, carnitine palmitoyl transferase 2 (CPT2), apolipoprotein A1 (APOA1), glutathione S-transferase theta 1 (GSTT1), major urinary protein 20 (MUP20), and nicotinamide phosphoribosyl transferase (NAMPT) in the kidney of mice among three groups are shown in Figure 3C. To further analyze the functions of these 30 DEPs, we performed GO annotations analysis (Figure S1A and B and Supplementary table 1). GO annotations were classified into three categories: BP, CC and MF. In the BP category, the most abundant DEPs were contained in cellular process, metabolic process and biological regulation. In the CC category, DEPs were mainly involved in the cell, cell part and organelle. In the MF category, DEPs displayed enrichment in binding and catalytic activity. In order to further understand the pathway information of these 30 DEPs, KEGG enrichment analysis was performed. The results of KEGG indicated that the majority of the DEPs were involved in chemical carcinogenesis, fatty acid metabolism and glutathione metabolism.

Figure 3 Identification of significant differentially expressed proteins (DEPs) in kidney of mice among three groups by LFQ-based proteomic analysis. (A) Venn diagram for overlapping proteins of the kidney among three groups (n=3). (B) Heatmap for the comparison of the DEPs among three groups (n=3). Red represents significantly up-regulated proteins, blue represents significantly down-regulated proteins, and gray represents no protein quantitative information. (C) The protein expression levels of SVS6, CPT2, APOA1, CSTT1, MUP20, and NAMPT in the kidney of mice among three groups were measured by proteomics analysis (n=3). Data are expressed as mean ± SD, *P<0.05, **P<0.01 vs control group; ^#P<0.05, ^##P<0.01 vs DKD group.

Abbreviations: GLP-1RA, glucagon-like peptide-1 receptor agonist; DKD, diabetic kidney disease; DEPs, differentially expressed proteins; LFQ, label-free quantification; SVS6, seminal vesicle secretory protein 6; CPT2, carnitine palmitoyl transferase 2; APOA1, apolipoprotein A1; GSTT1, glutathione S-transferase theta 1; MUP20, major urinary protein 20; NAMPT, nicotinamide phosphoribosyl transferase.

Validation of DEPs Among Three Groups of Renal Tissues and mTEC Cells by qPCR

QPCR analysis was used to verify the mRNA expression level of DEPs analyzed by the proteomics in vivo and in vitro. The results of qPCR showed that mRNA expressions of Svs6, Cpt2, Apoa1, Gstt1, Mup20, and Nampt were significantly increased in the kidneys of DKD group compared with the control group, whereas elevated expressions of Svs6, Cpt2, Apoa1, Gstt1 and Mup20 in renal tissues of the DKD mice were all dramatically decreased after GLP-1RA treatment. There was no significantly difference in the expression of Nampt between the DKD group and GLP-1RA group (Figure 4A). The qPCR verification results of the above DEPs in vivo were similar to the proteomics results. We also verified the mRNA expression levels of the above DEPs in three groups of mTEC cells, including the control group, HG group and GLP-1RA group. However, only the expression changes of Svs6 in the three groups of mTEC cells was consistent with the proteomics results (Figure 4B). As there is no commercial antibody of SVS6, the protein expression of SVS6 has not been verified by Western blot.

Figure 4 Validation of DEPs by qPCR in vivo and in vitro. (A) The mRNA expression of Svs6, Cpt2, Apoa1, Cstt1, Mup20, and Nampt in the kidney of mice among three groups was measured by qPCR (n = 10). (B) The mRNA expression of Svs6, Cpt2, Apoa1, Cstt1, Mup20, and Nampt in mTEC cells among three groups was measured by qPCR (n = 10). (C) The mRNA expression of Svs6 in mTEC cells after transfection with si-SVS6 by qPCR (n=10). (D) The mRNA expression of Kim-1, Ngal and Sglt2 in mTEC cells among NG group, HG group, GLP-1RA group and HG+si-SVS6 group (n=10). Data are expressed as mean ± SD, *P<0.05, **P<0.01 vs control group; ^#P<0.05, ^##P<0.01 vs DKD group or HG group.

Abbreviations: GLP-1RA, glucagon-like peptide-1 receptor agonist; DKD, diabetic kidney disease; DEPs, differentially expressed proteins; SVS6, seminal vesicle secretory protein 6; CPT2, carnitine palmitoyl transferase 2; APOA1, apolipoprotein A1; GSTT1, glutathione S-transferase theta 1; MUP20, major urinary protein 20; NAMPT, nicotinamide phosphoribosyl transferase.

In addition, mRNA expressions of Kim-1 and Sglt2 in the HG+si-SVS6 group were significantly lower than those in the HG group, while mRNA expression of Nagl in the HG+Si-SVS6 group showed a downward trend compared with the HG group (no statistical difference). Besides, there was no significant difference in the above renal tubular injury markers between the HG+ Si-SVS6 group and the GLP-1RA group (Figure 4C and D). Therefore, we speculated that reduction of SVS6 expression could alleviate HG-induced tubular injury in mTEC cells.

GLP-1R is Widely Expressed in Both Glomeruli and Renal Tubules of Mice

To explore whether GLP-1RA can directly acts on renal tubules through GLP-1R, we determined the expression of GLP-1R in the kidney of mice among three groups by immunofluorescence staining (Figure 5A). We found that GLP-1R was widespread expressed in the kidney of mice (including glomeruli and renal tubules) (Figure 5A). To further investigated the GLP-1R expression in renal tubules, we detected the expression of GLP-1R in mTEC cells (Figure 5B). Our data confirmed the widespread expression of GLP-1R in murine renal tubules in vitro (Figure 5B). The expression of GLP-1R in the ovary of mice was detected as a negative control according to a network bioinformatics tool (https://www.proteinatlas.org/) (Figure 5C).

Figure 5 GLP-1R is widely expressed in renal tubules of mice and mTEC cells. (A) Immunofluorescence staining of GLP-1R (green) was performed on the kidney sections of mice in the three groups (n = 3) (400×). (B) Immunofluorescence staining of GLP-1R (green) was performed on the mTEC cells (n = 3) (200×). (C) Immunofluorescence staining of GLP-1R (green) was performed on the ovary sections of mice as a negative control (n = 3) (400×).

Abbreviations: GLP-1RA, glucagon-like peptide-1 receptor agonist; DKD, diabetic kidney disease; GLP-1R, glucagon-like peptide-1 receptor.

Discussion

DKD complicates the course of numerous diabetic patients, imposes a huge economic burden on the family and society, and even increases the mortality of diabetic patients.^3,32 Although studies have shown that GLP-1RA may be beneficial to DKD through direct or indirect effects on the kidney, the specific mechanism remains unclear.¹⁵ Therefore, exploring the potential renal protection mechanism of GLP-1RA is of great significance for the effective treatment of DKD. In this study, a DKD mouse model induced by HFD combined with STZ was used. Our results showed that compared with the controls, the blood glucose, body weight, kidney weight, Scr and BUN were significantly increased in DKD group. In addition, mice in the DKD group showed significant mesangial expansion compared to the control group. According to the criteria for animal models of DKD developed by the Animal Models of Diabetic Complication Consortium (AMDCC),³³ we successfully established a mouse model of DKD. Our results also confirmed that GLP-1RA treatment significantly reduced blood glucose, body weight, kidney weight, Scr and BUN in DKD mice. Meanwhile, GLP-1RA treatment also dramatically improved the mesangial dilatation of DKD mice.

We also observed a series of structural changes in renal tubules of DKD mice, such as the disordered arrangement of renal tubular epithelial cells, vacuolar degeneration and tubular interstitial fibrosis. In addition, the mRNA levels of Fabp1, Kim-1, Ngal and Sglt2, which are key markers of renal tubular injury, in the kidneys of DKD group were significantly higher than those of the control group. Moreover, gene expressions of the above four renal tubular injury markers of mTEC cells were also significantly increased by high glucose stimulation. Therefore, we suggest that in addition to glomerular lesions, renal tubular damage also plays an important role in the pathogenesis of DKD. Several previous studies have suggested the importance of renal tubular injury in the pathogenesis of DKD. Gilbert et al¹⁰ emphasized the non-glomerular mechanism, especially tubular cell damage, in the pathogenesis of DKD, and indicated that renal tubular lesion could more accurately predict the changes in DKD renal function than glomerular damage. Moreover, Tang et al³⁴ proposed the concept of diabetic renal tubulopathy (DT), and suggested that tubulointerstitial injury in DKD was more serious than glomerular lesions, which possibly due to the multiple chemokines and injury signal factors secreted by the proximal tubular epithelial cells (PTECs), eventually lead to progressive interstitial inflammation and fibrosis. Therefore, renal tubular injury may become a new target for the therapy of DKD.⁶

Previous in vitro studies have demonstrated that hyperglycemia exerts both proinflammatory (via interleukin 6 and monocyte chemoattractant protein-1) and profibrotic effects (via the transforming growth factor-beta) on PTECs,^35,36 which are partly mediated via mitogen-activated protein kinase (MAPK) and protein kinase C (PKC) signaling.³⁵ Besides, other studies have revealed showed that HG induced renal tubular cell injury may be due to reduced autophagy.³⁷ Autophagy maintains the homeostasis and cell survival by degrading damaged organelles, so reduction or lack of autophagy may ultimately lead to cellular damage.³⁷ Meanwhile, overexpression of Klotho could dramatically enhance autophagy and AMP-activated protein kinase (AMPK) and extracellular signal-regulated kinase 1/2 (ERK1/2) activities in kidney tissues of T2DM mice induced by HFD/STZ and HG-induced PTECs.³⁷ In addition, bone morphogenetic protein 7 (BMP7) has been reported to inhibit advanced glycation end products (AGEs)-induced oxidative stress and multiple inflammatory signaling pathways, including p38 and p44/42 MAPK in PTECs.³⁸ However, so far, the molecular mechanism of renal tubular damage of DKD has not been fully determined.

Recently, GLP-1RA, an effective hypoglycemic agent, has been reported to have renal protection function in patients with T2DM.¹⁵ Our research showed that GLP-1RA treatment significantly improved renal tubular injury and tubulointerstitial fibrosis in DKD mice. Moreover, qPCR demonstrated that GLP-1RA could reduce the expression of renal tubule damage markers (Fabp1, Kim-1 and Sglt2) in the kidneys of DKD mice. In addition, our immunofluorescence staining confirmed the presence of GLP-1R in mouse glomeruli and renal tubules. Therefore, we speculated that GLP-1RA might improve renal function, renal tubular injury and fibrosis by directly acting on GLP-1R in renal tubules. A small amount of previous literature has revealed the effect of GLP-1RA on diabetic renal tubular damage. Serap et al²¹ observed that kidney damage mainly occurred in the renal tubule area of STZ-induced diabetic BALB/c mice rather than in the glomerulus, and found that Exendin-4 treatment improved renal tubular damage by reducing the production of reactive oxygen species and inflammation. Yijie et al²⁰ reported that Exendin-4 ameliorates the HG-induced FN1 and COL1A1 expression in PTECs by inhibiting the expression of miR-192, a p53-regulated microRNA that plays a role in renal fibrosis. Additionally, recombinant human GLP-1 has been shown to attenuate AGEs-induced tubulointerstitial injury by inhibiting phosphorylation of MAPK and nuclear factor-kappa B (NF-κB) in human proximal tubular cells (HK-2 cells).³⁹ Furthermore, liraglutide has been reported to mitigate renal fibrosis induced by unilateral ureteral obstruction and attenuate epithelial-mesenchymal transition (EMT) in recombinant transforming growth factor-beta 1 (TGF-β1)-treated renal tubular epithelial cells (NRK-52E).⁴⁰ However, the molecular mechanism of GLP-1RA in improving renal tubular damage of DKD remains controversial.

To investigate the mechanism of renal tubule injury in DKD and to elucidate the potential mechanism of GLP-1RA in improving renal tubular damage, we used LFQ-based proteomics to clarify the changes in kidney proteins among the control group, DKD group and GLP-1RA group. In this study, we identified a total of 30 DEPs of kidney among three groups. Then, we further validated several DEPs using qPCR, including SVS6, CPT2, APOA1, GSTT1, MUP20 and NAMPT. Among them, SVS6 was the most differentially expressed protein in the kidneys of the three groups of mice. Furthermore, the expression changes of Svs6 mRNA in mTEC cells among the control group, HG group and GLP-1RA group were consistent with LFQ-based proteomics results. When gene expression of Svs6 was inhibited by SVS6 siRNA in HG-stimulated mTEC cells, renal tubular damage markers were significantly decreased compared to the HG group. We speculated that SVS6 might be a candidate pathogenic protein for the development of DKD and a potential target for GLP-1RA in the treatment of DKD. Previous studies have found that the seminal vesicles of rodents can secrete 6 kinds of proteins, which are denoted as SVS1 to SVS6 in order of descending size.^41,42 It has been found in the past that the role of these proteins is to form mating plugs after mating.⁴¹ However, the expression and the role of SVS6 in the kidney have not been explored yet. We proposed for the first time that SVS6 was highly expressed in the kidneys of DKD mice and mTEC cells under HG conditions, whereas GLP-1RA reduced SVS6 expression and prevented tubular injury in vivo and in vitro. Therefore, SVS6 might be involved in the pathogenesis of DKD and in the treatment of DKD by GLP-1RA. However, further experiments are still needed to confirm the role of SVS6 in the pathogenesis of DKD.

Conclusions

In summary, our study found that there were abnormal renal function and obvious renal tubular damage in DKD mice, and GLP-1RA could significantly improve renal function and reduce renal tubular injury. The results of proteomics analysis indicated that the expression of SVS6 protein in the kidney of DKD group was significantly increased compared with the controls, and it was significantly decreased after GLP-1RA treatment. Changes in mRNA expression of SVS6 have been validated in vivo and in vitro. We also found that renal tubular damage markers were elevated in mTEC cells stimulated by HG, whereas tubular injury markers were obviously decreased after SVS6 expression inhibition by SVS6 siRNA in mTEC cells induced by HG. Our results suggest that SVS6 might be involved in the pathogenesis of DKD and in the treatment of DKD by GLP-1RA. Furthermore, we also found GLP-R was widely expressed in both glomeruli and renal tubules. Therefore, we speculate that GLP-1RA may directly bind to GLP-1R in the kidney to regulate the expression of SVS6, and ultimately improve renal tubular damage. SVS6 may be a potential therapeutic target for patients with DKD.

Abbreviations

GLP-1RA, glucagon-like peptide 1 receptor agonist; DKD, diabetic kidney disease; DEPs, the differentially expressed proteins; mTEC, murine kidney proximal tubular epithelial cells; GLP-R, glucagon-like peptide-1 receptor; eGFR, glomerular filtration rate; T2DM, type 2 diabetes mellitus; STZ, streptozotocin; SPF, the specific pathogen-free; HFD, high-fat diet; Scr, Serum creatinine; BUN, blood uric nitrogen; H&E, Hematoxylin and Eosin; PASM, Periodic Acid-Silver Methenamine; PAS, periodic acid-Schiff; LFQ, label-free quantification; GO, gene ontology; BP, biological process; MF, molecular function; CC, cell population points; KEGG, the Kyoto Encyclopedia of Genes and Genomes; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; FABP-1, fatty acid binding protein 1; KIM-1, kidney injury molecule-1; NGAL, neutrophil gelatinase-associated lipocalin; SGLT2, sodium-glucose co-transporter 2; COL1A1, collagen type I alpha 1 chain; COL3A1, collagen type III alpha 1 chain; FN1, fibronectin 1; SVS6, seminal vesicle secretory protein 6; CPT2, carnitine palmitoyl transferase 2; APOA1, apolipoprotein A1; GSTT1, glutathione S-transferase theta 1; MUP20, major urinary protein 20; NAMPT, nicotinamide phosphoribosyl transferase.

Ethics Approval

All animal experiments were approved by the Animal Ethics Committee of Tongji Hospital, Tongji University School of Medicine (2021-DW-(003)). This study was conducted in accordance with the Guide for the Care and Use of Laboratory Animals (NIH Publication, 8th Edition, 2011), and Animal Research: Reporting of In Vivo Experiments guidelines. All applicable international, national, and/or institutional guidelines for the care and use of animals were followed.

Disclosure

The authors declare that they have no conflicts of interest for this work.

Additional information

Funding

This study was sponsored by the National Natural Science Foundation of China (Grant No. 81974105, 81670805, 81800705), Shanghai Pujiang Program (Grant No. 2019PJD050), Novo Nordisk China Diabetes Young Scientific Talent Research Funding, Psychosomatic Medicine Project of Key Developing Disciplines of Shanghai Municipal Health Commission (2019ZB0202), the Science and Technology Innovation Program of Shanghai (Grant No. 16411954700), Scientific research project of Shanghai Municipal Health Commission (Grant No. 20194Y0460), Cultivation project for National Natural Science Foundation of Shanghai Tongji Hospital, and Doctoral Fund of Northern Jiangsu Hospital (Grant No. BSQDJ0142).