Skip to main content

Single-cell transcriptomics reveals a mechanosensitive injury signaling pathway in early diabetic nephropathy

Abstract

Background

Diabetic nephropathy (DN) is the leading cause of end-stage renal disease, and histopathologic glomerular lesions are among the earliest structural alterations of DN. However, the signaling pathways that initiate these glomerular alterations are incompletely understood.

Methods

To delineate the cellular and molecular basis for DN initiation, we performed single-cell and bulk RNA sequencing of renal cells from type 2 diabetes mice (BTBR ob/ob) at the early stage of DN.

Results

Analysis of differentially expressed genes revealed glucose-independent responses in glomerular cell types. The gene regulatory network upstream of glomerular cell programs suggested the activation of mechanosensitive transcriptional pathway MRTF-SRF predominantly taking place in mesangial cells. Importantly, activation of MRTF-SRF transcriptional pathway was also identified in DN glomeruli in independent patient cohort datasets. Furthermore, ex vivo kidney perfusion suggested that the regulation of MRTF-SRF is a common mechanism in response to glomerular hyperfiltration.

Conclusions

Overall, our study presents a comprehensive single-cell transcriptomic landscape of early DN, highlighting mechanosensitive signaling pathways as novel targets of diabetic glomerulopathy.

Background

Diabetic nephropathy (DN) is one of the major complications in diabetic patients and is the leading cause of end-stage renal disease worldwide [1]. DN is a complex disease and its progression, particularly in type 2 diabetes, is confounded by multiple pathogenic factors, such as insulin resistance and obesity [1, 2]. Previous studies used single-cell or single-nucleus RNA sequencing (scRNA-seq or snRNA-seq) to investigate cellular changes or drug responses in diabetic kidneys. Single-nucleus transcriptomics on early human diabetic kidneys suggests increased potassium secretion in the distal nephron and pro-angiogenic signaling in diverse kidney cells [3]. A following study combining single-nucleus RNA and assay for transposase-accessible chromatin (ATAC) sequencing highlights the glucocorticoid signaling in the proximal tubule [4]. Stefansson et al. provide molecular programs associated with glomerular hyperfiltration in human DN [5], and Fu et al. show the single-cell transcriptomics of eNOS − / − mice with streptozotocin-induced diabetes, which is a type 1 diabetes model without obesity nor insulin-resistance [6]. Single-cell/nucleus transcriptomics on drug-treated murine diabetic kidneys suggest important alterations in proximal tubule cells [7, 8].

DN is characterized by distinct histopathological lesions dominating the renal glomerulus, such as glomerular basement membrane (GBM) thickening, mesangial expansion, and nodular glomerulosclerosis [2]. Early changes such as glomerular hyperfiltration have been proposed to be critical for the subsequent development of glomerulosclerosis [1]. Preservation of glomerular filtration rate is of great importance for renoprotection [2]. However, the signaling pathways that initiate these glomerular alterations are incompletely understood. Therefore, it is particularly warranted to understand the cellular changes and molecular mechanisms that initiate DN in the glomerulus.

In this study, we conducted a comprehensive analysis of cellular changes in type 2 diabetic glomeruli. The type 2 diabetes mouse model BTBR ob/ob exhibits insulin resistance, hyperglycemia, and obesity, as well as the rapid progression of diabetic glomerulopathy [9,10,11]. We generated single-cell transcriptomics of kidneys and bulk transcriptomics of purified glomerular cell types. We compared mouse DN with human DN and validated our findings in different patient cohort datasets at both mRNA and protein levels. Our study reveals the upregulation of the mechanosensitive signaling pathway MRTF-SRF in response to glomerular hyperfiltration, which is associated with glomerulopathy in early DN.

Methods

Animals

All animal experiments were conducted according to the National Institutes of Health Guide for the Care and Use of Laboratory Animals as well as the German law for the welfare of animals. Mice were housed in a specific pathogen-free facility with free access to chow and water and a 12-h/12-h day/night cycle. Breeding and genotyping were performed according to standard procedures.

BTBR ob/ob podocyte-reporter mice

Animal experiments were approved by the veterinary administration of the City of Hamburg under the license Ü003-2018. BTBR ob/ + (BTBR.Cg-Lepob/wt WiscJ) heterozygous animals (Jax No. 004824) and corresponding wild-type mice BTBR (Jax No 002282) were purchased from JAX (Bar Harbor, ME, USA). Podocyte-reporter mice (Gt(ROSA)26Sortm4(ACTB−tdTomato,−EGFP)Luo;Tg(NPHS2-cre)295Lbh) were also purchased from JAX (Jax No. 007576) and crossed for at least 7 generations with BTBR wild-type animals before BTBR.Cg-Lepob/wt WiscJ;Gt(ROSA)26Sortm4(ACTB−tdTomato,−EGFP)Luo;Tg(NPHS2-cre)295Lbh animals were generated. The resultant offspring was intercrossed to yield BTBR.Cg-Lepob/ob WiscJ;Gt(ROSA)26Sortm4(ACTB−tdTomato,−EGFP)Luo;Tg(NPHS2-cre)295Lbh and BTBR.Cg-Lepwt/wt WiscJ;Gt(ROSA)26Sortm4(ACTB−tdTomato,−EGFP)Luo;Tg(NPHS2-cre)295Lbh animals.

Mouse kidney normothermic machine perfusion

All procedures involving mouse kidney perfusion described in this manuscript were conducted according to German and Hamburg law and approved by the veterinary administration of the City of Hamburg under the license N002/2020. As previously described [12], male C57BL/6 kidneys were collected under anesthesia (i.p. injection of 100 μL/10 g body weight of a solution containing 10 mg/ml Ketamine and 1.6 mg/ml Xylazine). Briefly, a laparotomy was performed, the aorta was cannulated, ligatures placed and closed around all main vessels, and perfusion initiated only in the right kidney. The left unperfused kidney was then removed and used as control. After this, the right kidney was removed from the mouse during ongoing perfusion. Perfusion took place ex vivo at 100 mmHg at 37 °C with 100% GBSS solution without the addition of albumin or red blood cells (Sigma G9779) continuously gassed in a dialyzer with 100% carbogen (95% O2, 5% CO2). Perfusion was stopped after 60 min. Isolation of glomeruli was performed as described in Additional file 1: Supplementary methods (FACS-sorted glomerular cells).

Pig kidney normothermic machine perfusion

Pig kidneys were collected from a German slaughterhouse. Kidneys were received directly after slaughter and immediately flushed with 400–500 ml 1 × PBS with 2 ml Heparin 10.000 IE/ml. Kidney were stored at + 4 °C for around 2 h, after which normothermic perfusion was initiated. Perfusion took place at 100 mmHg at 37 °C with increased glomerular pressure and hyperfiltration as described above. Renal cortex tissue was collected either before perfusion (control) or 120 min after start of perfusion (perfused).

Single-cell/nucleus suspension and library preparation

Single-cell suspension

Glomeruli were enriched by depleting the PT segment from the kidneys using Percoll density gradient centrifugation described previously [13]. Briefly, the mouse was euthanized by cervical dislocation under anesthesia (4% isoflurane for 2 min). Both kidneys were harvested and the renal papilla of each kidney was removed. The renal tissue was minced into small pieces and transferred to a C tube (130–093-237, Miltenyi Biotec) filled with 5 ml of collagenase solution (1 mg/ml collagenase from Clostridium histolyticum [C9263, Sigma-Aldrich] and 0.25% bovine serum albumin [BSA] in DMEM/F-12 medium [11320074, Gibco]). The renal tissue was incubated at 37 °C for 5 min and homogenized using a gentleMACS™ Dissociator (Miltenyi Biotec) with the program “m. spleen 01_01”; this process was repeated once. Subsequently, the renal tissue was passed through a 250-μm sieve and separated by Percoll density gradient centrifugation using a solution that contained 45% Percoll® (GE17-0891–02, GE Healthcare) and 55% 2 × PBS-glucose at 17,500 rpm for 15 min at 4 °C. Red blood cells were removed with a Percoll density gradient. PTs that accumulated in the interphase of the gradient were discarded. The rest of the kidney tissue that accumulated in the upper phase of the gradient was taken, washed twice in ice-cold HBSS, and centrifuged at 400 × g for 4 min at 4 °C. A single-cell suspension was prepared using cold-active protease (CAP) as described in previous studies [14]. The pellet was suspended in 2 ml of CAP solution (10 mg/ml Bacillus licheniformis protease [P5380, Sigma-Aldrich], 5 mM CaCl2, 20 U/ml DNase [4716728001, Roche] in PBS) for 20 min on ice with repeated trituration steps for 20 s every 5 min. The digestion was neutralized by 15 ml PBS supplemented with 10% FBS and passed through a 40-μm Corning® cell strainer. Cells were centrifuged at 400 × g for 4 min at 4 °C and washed twice in 20 ml PBS supplemented with 0.5% BSA. Dead cells were removed with a dead cell removal kit (130–090-101, Miltenyi Biotec) according to the manufacturer’s protocol. The living cells were passed through a 30-μm cell strainer (04–004-2326, Sysmex), and the cell concentration was determined using a TC20™ automated cell counter (Bio-Rad). The cells were loaded onto a 10 × Genomics Chromium single-cell instrument. All steps were performed according to the standard protocol of the Chromium single-cell 3’ v2/v3 kits to generate high-quality cDNA libraries.

Single-nucleus suspension

A 1 × 3 mm piece of snap frozen renal cortex tissue was thawed and chopped with a razor blade in a petri dish on ice and homogenized using a Dounce homogenizer (D8938-1 SET, Sigma-Aldrich) in 200-μl ice-cold lysis solution and incubated on ice for 20 min with additional 3.8 ml of ice-cold lysis solution. Lysis solution was prepared with Nuclei PURE lysis buffer (NUC-201, Sigma-Aldrich), 1 mM dithiotreitanol (D9779, Sigma-Aldrich), and 0.1% Triton X-100 (NUC-201, Sigma-Aldrich) according to manufacturer’s protocol and a RNAse inhibitor mix (0.04 U/μl SUPERaseIN RNAse Inhibitor [AM 2696, Thermo Fisher]; 0.04 U/μl RNAsin Plus RNAse Inhibitor [N2615, Promega]) was added. The single nuclei suspension was filtered through a 30-μm filter (04–004-2326, Sysmex) and centrifuged at 500 g for 5 min at 4 °C. The pellet was resuspended and incubated for 2 min in 1 ml red blood cell lysing buffer hybri-max™ (R7757-100 ml, Sigma-Aldrich), filtered through a 5-μm filter (04–004-2323, Sysmex), and washed with 4 ml of ice-cold 0.01% BSA (AM2616, Thermo Fisher) in DPBS (59331C; Sigma) with 0.04 U/μl SUPERaseIN RNAse Inhibitor and 0.04 U/μl RNAsin Plus RNAse Inhibitor at 500 g for 5 min at 4 °C. The pellet was resuspended in 1% BSA in DPBS with 0.04 U/μl SUPERaseIN RNAse Inhibitor and 0.04 U/μl RNAsin Plus RNAse Inhibitor, and nuclei number was counted. The libraries were prepared with the Chromium NEXT GEM Single Cell 3’ Reagent kits v3.1 according to manufacturer’s protocol. The libraries were sequenced on an Illumina Novaseq6000 platform as symmetrically paired end runs (150 bases) with 200 million raw sequencing reads per sample.

Single-cell data analysis

Preprocessing and quality control (QC) of scRNA-seq data

10 × Genomics raw sequencing data were processed using CellRanger software (version 3.0.2, 10 × Genomics, Pleasanton, CA), and the 10 × Genomics mouse genome mm10 3.0.0 release was used as the reference genome (function cellranger count). The matrices of cells and the unique molecular identifier (UMI) count were obtained and further processed with the R package Seurat (version 3.1.1) [15]. For QC, we first filtered out genes detected in fewer than 3 cells and data for cells in which fewer than 200 genes had nonzero counts. To remove potential doublets, cells with more than 7000 expressed genes (nFeature) were excluded. We removed low-quality cells with more than 50% mitochondrial genes among all detected genes, as is recommended for kidney tubular cells [16]. After clustering and cell type identification, we performed curated doublet and high-mitochondrial gene cell removal (see below) based on known lineage-restricted markers.

Dimensionality reduction and clustering

The Seurat R package (version 4.0.2) was used to perform unsupervised clustering analysis on scRNA-seq data. In brief, gene counts for cells that passed QC were normalized by library size and log-transformed (function NormalizeData, normalization.method = “LogNormalize”, scale.factor = 10,000). Then, highly variable genes were detected (function FindVariableFeatures, selection.method = “vst”, nfeatures = 2000). To reduce batch effects, we applied the integration method implemented in Seurat version 3 (functions FindIntegrationAnchors and IntegrateData, dims = 1:30). The integrated matrix was then scaled with the ScaleData function (default parameters). PCA was performed on the scaled data (function RunPCA, npcs = 30) to reduce dimensionality.

The number of principal components used for each clustering round was dataset dependent and was determined on the basis of the elbow of a PCA scree plot. The selected principal components were then used to compute the KNN graph based on the Euclidean distance (function FindNeighbors). Cell clusters were subsequently generated using the function FindClusters. The resolution of the FindClusters function for each dataset was also determined by the exploration of the top marker genes of each cluster. Uniform Manifold Approximation and Projection (UMAP) was used to visualize the clustering results. The top DEGs in each cluster were found using the FindAllMarkers function (min.pct = 0.25, logfc.threshold = 0.25) with Wilcoxon rank-sum tests. The most highly expressed genes were then used to determine the cell type of each cluster.

Curated doublet and high-mitochondrial gene cell removal

After the cell type was determined for clusters, we performed additional curated doublet and high-mitochondrial-cell removal. Based on the literature and the exploration of our datasets, we created a lineage-restricted marker gene list for tubular (Cubn, Epcam) and nontubular (Pecam1, Pdgfrb, Nphs2, Ptprc) cell types. We removed the cells that expressed markers of the opposite lineage. For nontubular cells, we discarded cells with more than 20% mitochondrial genes among all detected genes. There were 86,508 cells before curated doublet removal, 71,831 cells after curated doublet removal, and 770,944 cells after high-mitochondrial gene cell removal.

Differential gene expression analysis

The Seurat FindMarker function was used to perform differential gene expression analysis for each cell type between the control and DN groups. The test method used in the FindMarker function was MAST (v1.12.0) [17].

Single-nucleus data analysis

Preprocessing and QC

10 × Genomics raw sequencing data were processed using CellRanger software (version 5.0.1, 10 × Genomics, Pleasanton, CA). Pig genome was built by function “cellranger mkref” using the Sscrofa v103 FASTA file and GTF file (https://www.ensembl.org) and then mapping was done and the count matrices were generated by function “cellranger count” with parameter “include-introns.” For QC, firstly, we applied soupX [18] with default parameters to remove ambient RNA contamination. Then, the nuclei were filtered out if the number of genes detected was less than 500 or greater than 6000 or the percentage of mitochondrial genes detected exceeded 7%. To further remove potential doublets, Scrublet [19] was applied to the data and cells were excluded if they were identified as doublets (default parameters).

Dimensionality reduction and clustering

An updated version of Seurat R package (version 4.0.2) was used to perform unsupervised clustering analysis on snRNA-seq data since the experiments were performed later than the mouse scRNA-seq analysis. Count data normalization, scaling, highly variable gene selection, and sample integration were the same as the mouse single-cell data process described above. PCA was performed on the scaled data (function RunPCA, npcs = 30) and the first 15 PCs were used for clustering. The KNN graph was calculated (function FindNeighbors) and then the clustering result is obtained (function FindClusters, resolution = 0.5). To visualize the clustering result, UMAP coordinates were calculated (function RunUMAP, dims = 1:15). Wilcoxon method was used to perform the differential gene analysis. For each cluster, the marker gene list was determined by log2FoldChange > 0.25 and adj. p value < 0.05 (FindAllMarkers). The most highly expressed genes were then used to determine the cell type of each cluster.

Bulk RNA-seq

Library preparation and RNA-seq

A small amount of RNA (2 ng) was used as input material, and libraries were prepared with a SMART-Seq Stranded Kit according to the user manual (Takara Bio USA, Mountain View, CA, USA). In brief, samples were fragmented at 85 °C for 6 min prior to first-strand synthesis. Illumina adaptors and indexes were added to single-stranded cDNA via 5 cycles of PCR. After library purification with AMPure beads and depletion of ribosomal cDNA with scZapR, final RNA-seq library amplification (13 cycles) was conducted, and the final RNA-seq library was purified with AMPure beads. The library samples were quantified using Quant-iT PicoGreen dsDNA Reagent (Invitrogen; Thermo Fisher Scientific, Waltham, MA, USA) on a ClarioStar microplate reader according to the manufacturer’s instructions (BMG LABTECH, Ortenberg, Germany). The quality, including fragment size, of the cDNA was assessed on an Agilent Technologies Bioanalyzer 2100 using an Agilent DNA 1000 kit according to the manufacturer’s instructions (Agilent Technologies, Palo Alto, CA, USA).

Pooled samples were quantified with a Qubit 1X dsDNA HS Assay Kit on a Qubit fluorometer (Thermo Fisher Scientific, Waltham, MA, USA). Single-read sequencing was performed on a NovaSeq 6000 device using an S2 Reagent kit (100 cycles) according to the manufacturer’s instructions (Illumina Inc., CA, USA).

Bulk RNA-seq data analysis

The quality of the bulk RNA-seq reads was assessed using FastQC (v0.11.5), and the reads were aligned to the mouse reference genome (mm10) with Bowtie2 (v2.3.3.1) [20] using RSEM (v1.3.0) [21] with the default parameters. The function rsem-calculate-expression was used to align the reads and quantify the gene and isoform abundance. The output of rsem-calculate-expression separately gives the read count and transcripts per million (TPM) value for each gene and isoform. Differential expression analysis was carried out using gene read counts with the DESeq2 package (v1.22.2) [22] to produce log2FC values and corresponding p values and adjusted p values. Principal component analysis (PCA) was performed using regularized log transformation of the count data, and the results were visualized using gplots (v3.0.1.1).

Integrated analysis

Correlation of bulk RNA-seq and scRNA-seq data

For podocytes, ECs, and mesangial cells, the single-cell data were correlated with bulk RNA-seq data. The cell type markers (top DEGs of each cell type identified by the Seurat FindAllMarkers function (min.pct = 0.25, logfc.threshold = 0.25) were used for correlation analysis. The normalized gene expression count matrices of cell type markers were extracted from the single-cell data and bulk RNA-seq data separately, and then Pearson correlation coefficients were calculated between every pair of single cells and replicates of bulk RNA-seq data.

CellChat

We applied Cellchat (v 1.4.0) [23] to infer cell–cell communications across all kidney cell types and the glomerular cell types. We used the mouse database curated in CellChat including the “Secreted Signaling,” “Cell–cell Contact,” and “ECM-Receptor.” The comparison analysis between diseased samples and control samples was performed according to the Cellchat tutorial (https://htmlpreview.github.io/?https://github.com/sqjin/CellChat/blob/master/tutorial/Comparison_analysis_of_multiple_datasets.html). The communication probabilities were compared to identify the unregulated and downregulated signaling ligand-receptor pairs between the two conditions.

SCENIC

The gene regulatory network was inferred using pySCENIC (v 0.11.2, a lightning-fast python implementation of the SCENIC pipeline) [24, 25]. Firstly, the GRNboost2 algorithm was used to infer gene regulatory network and generate co-expression modules (pyscenic grn). Next, the regulon prediction step was performed (pyscenic ctx, using the default parameters and mm10__refseq-r80__10kb_up_and_down_tss.mc9nr and mm10__refseq-r80__500bp_up_and_100bp_down_tss.mc9nr motif collections). Finally, the AUCell matrix was generated with a threshold of 0.01 (pyscenic aucell). Cell type-specific regulators were identified based on Z-score of the AUCell values for the cells of a given type.

Upstream analysis

We performed Qiagen’s ingenuity Pathway Analysis (IPA) tool to identify transcriptional regulators (significance: p value < 0.05). Positive and negative z-score values indicated the activation or inhibition of transcriptional regulators, respectively.

Enrichment analysis

GO biological processes, KEGG, WikiPathways, and Reactome Gene Sets were performed with Metascape [26].

Human data

ERCB-Kröner-Fresenius Biopsy Bank

Human renal indication biopsy specimens were collected in an international multicenter study and deposited in the ERCB-Kröner-Fresenius Biopsy Bank (ERCB-KFB [27, 28]; for participating centers, see Shved et al. [29]). Glomerular samples from patient with different renal diseases were analyzed for mRNA expression levels (GSE32591, GSE35489, GSE37463, GSE47185, GSE99340). The analysis included gene expression profiles from patients with DN (n = 14), hypertensive nephropathy (HT; n = 15), minimal change disease (MCD; n = 14), IgA nephropathy (IgA; n = 27), focal segmental glomerulosclerosis (FSGS; n = 23), membranous nephropathy (MGN; n = 21), lupus nephritis (SLE; n = 32), ANCA-associated glomerulonephritis (RPGN; n = 23), and controls (LDs; n = 41).

Early diabetic nephropathy in Pima Indians

Protocol kidney biopsy specimens were collected from Pima Indians with type 2 diabetes. The study subjects participated in an intervention trial (Renoprotection in Early Diabetic Nephropathy in Pima Indians, clinicaltrials.gov, NCT00340678) [30, 31]. Biopsies were obtained from patients after informed consent was obtained and with the approval of the local ethics committees. The analysis included gene expression profiles from DN patients (n = 68) and LD controls (n = 18). To account for ethical considerations, to ensure privacy protection, and to avoid identifying individual study participants in this vulnerable population of Pima people, the Institutional Review Board of the National Institute of Diabetes and Digestive and Kidney Diseases has stipulated that individual-level gene expression and genotype data from this study cannot be made publicly available.

Microarray analysis

Tissue processing and microdissection protocols were performed similarly on both ERCB and Pima biopsy material. Following renal biopsy, the tissue was transferred to an RNase inhibitor and microdissected into glomeruli and tubulointerstitial tissue. Total RNA was isolated, reverse-transcribed, and amplified. Fragmentation, hybridization, staining, and imaging were performed according to the Affymetrix Expression Analysis Technical Manual (Affymetrix, Santa Clara, CA, USA). CEL file normalization was performed with the robust multichip average method using RMAExpress (Version 1.0.5) and the human Entrez-Gene custom CDF annotation from Brain Array version 18 (PIMA) and 25 (ERCB) (http://brainarray.mbni.med.umich.edu/Brainarray/default.asp). The log-transformed ERCB dataset was corrected for batch effect using ComBat from the GenePattern pipeline (http://www.broadinstitute.org/cancer/software/genepattern/). To identify DEGs, the Significance Analysis of Microarrays (SAM) method [32] was applied using the SAM function in Multiple Experiment Viewer (TiGR MeV, Version 4.9). A q-value below 5% was considered to indicate statistical significance.

Immunofluorescent staining

Paraffin sections from human biopsies or from experimental mice (3 μm thick) were deparaffinized and rehydrated in water. Following antigen retrieval with DAKO pH6 for 30 min at 98 °C, unspecific binding was blocked for 30 min at RT in blocking buffer (0.05% TritonX-100, 5% normal horse serum (Vector) in PBS). Primary antibodies for human biopsies (MRTFA [Sigma HPA030782]; MRTFB [Invitrogen PA5-113519]; CD34 [Leica Biosystems NCL-L-END]; NPHS1 [Progene GP-N2]) and for murine sections (MRTFA [Sigma HPA030782]; MRTFB [Invitrogen PA5-113519]; Endomucin-AF546 [Santa Cruz sc-65495]; NPHS1 [Progene GP-N2]) were diluted in blocking buffer and incubated overnight at 4 °C. Stainings were visualized following incubation with appropriate fluorochrome coupled secondary donkey antibodies (Jackson Immunoresearch Laboratories) and Hoechst (Molecular Probes) and mounted with fluoromount (Thermo Fisher). Analyses were performed on a LSM800 airyscan 1 using the ZENBlue software (all ZEISS).

A piece of pig cortex (1 cm deep) was snap frozen in OCT. Seven-micrometer cryo-sections were fixed immediately in 4% PFA for 10 min at RT and washed out with PBS. The sections were incubated for 30 min with 5% BSA in PBS + 0.1% Triton for blocking and permeabilization and then incubated overnight at 4 °C with primary antibody MEIS1 (Invitrogen MA5-27,191). After washing out the primary antibody, sections were incubated for 1 h at RT with secondary antibody (Life Technologies) and DAPI (4′,6-diamidino-2-phenylindole) (Invitrogen). Slides were mounted and analyzed using the Leica TCS SP5 microscope.

Results

scRNA-seq of BTBR ob/ob mouse kidneys with early DN

We generated BTBR ob/ob podocyte-reporter mice (BTBR.Cg-Lepob/ob WiscJ;Gt (ROSA)26Sortm4(ACTB−tdTomato,−EGFP)Luo;Tg(NPHS2-cre)295Lbh; hereafter referred to as BTBR ob/ob or DN mice) to investigate the cellular and molecular changes that occur in DN. By 6 weeks of age, both male and female BTBR ob/ob mice exhibited obesity, hyperglycemia, and albuminuria (Additional file 1: Figure S1a). From 12 weeks of age, obvious glomerular hypertrophy was observed in histological analysis (Additional file 1: Figure S1b). Mesangial expansion with accumulation of collagen IV was detected in the DN glomeruli (Additional file 1: Figure S1c and d). Based on these observations, we defined 6 weeks of age as the onset stage of DN and 12 weeks of age as the early stage of DN.

Kidney tissue was sampled from a total of 16 mice, including male and female each of BTBR ob/ob (DN) and BTBR WT (control) mice at 6 and 12 weeks of age. The samples were partially depleted of the proximal tubule (PT) segment using a Percoll gradient [13] to enrich glomeruli. Single-cell suspensions were prepared using cold-active protease (CAP) to reduce dissociation-induced artifacts in kidney cells [33]. A total of 70,944 single cells were profiled after data pre-processing and quality control (Fig. 1b). Unsupervised clustering identified 21 clusters (Additional file 1: Figure S2a, b and c; see the “Methods” section). The cells in the 21 clusters were classified into 18 cell types and annotated on the basis of cell-specific marker genes reported in previous kidney and glomerulus single-cell data [14,15,16]. Defining marker genes for each cell type and unbiased marker genes for each cluster are demonstrated in a dot plot (Fig. 1c) and in a heatmap (Additional file 1: Figure S2d), respectively.

Fig. 1
figure 1

scRNA-seq of BTBR ob/ob mouse kidneys with early DN. a Experimental scheme. b UMAP plot of annotated cell types. c Dot plot of defining marker genes for each cell type. d Total number of significant DEGs (nDEGs) in each cell type. Top enriched pathways in glomerular e and PT cell types f at 6 and 12 weeks. UMAP, uniform manifold approximation and projection; Ctrl, control; DN, diabetic nephropathy; Podo, podocyte; EC, endothelial cell; Mesan, mesangial cell; Int, interstitial cell; PT, proximal tubule; dLOH, descending limb of loop of Henle; aLOH, ascending limb of loop of Henle; DCT, distal convoluted tubule; CNT, connecting tubule; PC, collecting duct principal cell; IC-A, A-type collecting duct intercalated cell; IC-B, B-type collecting duct intercalated cell; Trans, transition cell; Imm, immune cell; Mitotic, mitotic cell

We identified all major components in the kidneys. The endothelial cell (EC) subtypes expressed distinct marker genes, such as Ehd3 for EC_1 (glomerular ECs/gECs), Plvap for EC_2 (fenestrated ECs with diaphragm in veins and peritubular capillaries), and Fbln2 for EC_3 (arteriolar and arterial ECs) [34]. Mesangial cells defined by marker genes Gata3, Itga8, and Adamts5 are clearly distinguished from interstitial cells, which were composed of vascular smooth muscle cells marked by Acta2 and Tagln, and renin cells marked by Ren1 [35, 36]. The subtype PT_1 and PT_2 were composed of cells from S1/2 and S2/3 PT segments, respectively, as they showed a continuum of marker gene expression consistent with known S1 (Slc5a2, Slc5a12 and Gatm) to S3 (Slc13a3, Cyp4a10, and Slc5a8) PT segments.

Due to PT depletion, the percentage of PT cells was reduced to 22.4%, while glomerular cell types were enriched to 15.6% (Additional file 1: Figure S2e). The fractions of control and DN cells in each cell type were similar except for podocytes, which included more control (approximately 77%, 807 cells) than DN cells (approximately 23%, 235 cells) (Additional file 1: Figure S2f). All cell types were identified in individual kidney samples (Additional file 1: Figure S3a), and cells from individual kidney samples were also identified for each cell type (Additional file 1: Fig. S3b). The overall low expression of dissociation-induced stress genes in individual cell types (Additional file 1: Fig. S3c) indicates the advantages of our modified single-cell preparation protocol (Methods), which shortened the 37 °C digestion time and omitted the steps stressful to glomerular cell dissociation.

Cell type-specific alterations in DN kidneys suggest glucose-dependent and glucose-independent responses in diabetes

Next, we compared gene expression levels in DN versus control kidneys by cell type. In total, we identified 447 and 740 differentially expressed genes (DEGs) across all cell types at 6 and 12 weeks, respectively (Fig. 1d). Among all cell types, mesangial cells and PT cells had the highest numbers of DEGs (Fig. 1d), indicating that these two cell types were predominantly affected in the early stage of DN. Enrichment analysis identified significant pathways in glomerular cell types and PT cell types at 6 and 12 weeks (Fig. 1e and f). Glomerular cell types showed regulations of extracellular matrix molecules, cell adhesion, cell proliferation, etc., and cellular responses to mechanic stress, cytokines, and blood pressure. PT cell types showed multiple metabolic processes such as biological oxidations, and responses involved in oxidative stress, such as NRF2 and PPAR signaling pathways. These differences suggest glucose-dependent and glucose-independent responses in PT and glomerular cell types in early DN, respectively.

Shared features of human and experimental DN

Early changes such as glomerular hypertrophy, GBM thickening, and mesangial expansion occur in both human and animal DN [37]. Thus, we conducted intensive analyses on glomerular cell types.

We isolated glomeruli from BTBR ob/ob and BTBR WT kidneys and subsequently purified podocytes, mesangial cells, and gECs by fluorescence-activated cell sorting (FACS) (Additional file 1: Figure S4, Additional file 1: Supplementary methods). Bulk RNA-seq was performed for each purified glomerular cell type. The similarity and purity of individual glomerular cell type samples were controlled by principal component analysis (PCA) and the expression levels of specific marker genes (Additional file 1: Figure S5). The cell gene expression data from scRNA-seq were compared to the bulk RNA-seq data. Podocytes, EC_1/gECs, and mesangial cells showed positive Pearson correlation coefficients, suggesting a linear correlation of the two datasets (Fig. 2a). Overlapping the glomerular DEGs from single-cell and bulk RNA-seq datasets resulted in 194 unique genes (Fig. 2b), which were mapped to the European Renal cDNA Bank (ERCB) patient dataset. The mRNA expression of 106 genes was significantly regulated in microdissected glomeruli from DN patients (Fig. 2b and c). Importantly, most DEGs detected in human DN glomeruli were significantly changed in mesangial cells in both single-cell and bulk RNA-seq data (Fig. 2c).

Fig. 2
figure 2

Shared features of human and experimental DN. a Heatmap displaying the Pearson correlation coefficient between scRNA-seq and bulk RNA-seq data. b Overview of significant glomerular DEGs detected in both single-cell and bulk RNA-seq datasets, as well as in the ERCB patient dataset. cHeatmap showing the significant regulation of DEGs identified in microdissected glomeruli and in glomerular cells from DN patients and mice, respectively. Nonsignificant genes are shown in grey. ERCB, European Renal cDNA Bank

Mechanosensitive transcriptional regulators are activated in DN glomeruli

To identify the transcriptional regulations responsible for the differential expression of genes in DN glomerular cell types, we performed upstream analysis using Ingenuity Pathway Analysis (IPA), which identifies not only transcription factors but also transcription coregulators such as coactivators. Transcriptional regulators commonly changed in single-cell and bulk RNA-seq data are shown in Fig. 3a. Changes in transcriptional regulations were predominant in mesangial cells. Intriguingly, we observed the activation of mechanosensitive transcriptional regulators, including coactivators myocardin-related transcription factor A and B (MRTFA/B) and yes-associated protein 1 (YAP1 or YAP), as well as the transcription factor serum response factor (SRF). Additionally, the gene regulatory network was inferred using SCENIC [24, 25] to determine the cell-type specific activity of transcription factors (Additional file 1: Figure S6b). Consistently, the activity of SRF was high in mesangial cells and moderate or low in podocytes and gECs, respectively (Fig. 3b).

Fig. 3
figure 3

Mechanosensitive transcriptional regulators are activated in DN glomeruli. a Dot plot displaying the commonly changed transcriptional regulations in both single-cell and bulk RNA-seq data estimated by IPA. b Heatmap showing the activity z scores of transcription factors estimated by SCENIC. c Heatmap showing the log2FC (DN vs. Ctrl) values of MRTF transcriptional target genes in mice glomerular cells. Nonsignificant genes are shown in gray (a, c). d Immunofluorescence staining of the mouse kidney paraffin sections showing the expression and localization of MRTFA and MRTFB in control and DN glomeruli. Glomerular endothelial cells were labeled by EMCN (endomucin), podocytes were labeled by NPHS1 (Nephrin), nuclei were counterstained with HOECHST. Arrows indicate mesangial cells, which were negative for EMCN and NPHS1

The MRTF-SRF transcriptional target genes with significant regulations in DN glomerular cell types are shown in Fig. 3c. Notably, the MRTF-SRF-specific targets Acta2 and Tagln [38] were found in both single-cell and bulk RNA-seq data. Enrichment analysis of significantly changed MRTF transcriptional target genes suggested regulations of smooth muscle proliferation, actin cytoskeleton, integrin cell surface interactions, and extracellular matrix organization (Additional file 1: Figure S6c).

MRTF-SRF is an important mechanosensitive transcriptional pathway in cells [39]. MRTF family consists of myocardin, MRTFA (MKL1, MAL), and MRTFB (MKL2). Unlike myocardin, which is limited to myocardial cells and smooth muscle cells, MRTFA and MRTFB are found in various cells and tissues including the kidney [40]. Mechanical stresses lead to the translocation of cytoplasmic MRTF into the nucleus to initiate the SRF-mediated transcriptional responses [40]. We investigated the expression and location of MRTFA and MRTFB in mouse DN by immunofluorescent staining. MRTFA and MRTFB were found in the cytoplasm and nuclei of podocytes and gECs, exhibiting comparable expression levels in both control and DN mice (Fig. 3d). Mesangial cells exhibited a low expression level of MRTFA in control mice, while some of the mesangial cells exhibited nuclear accumulation of MRTFA in DN mice (Fig. 3d). MRTFB expression was absent in most mesangial cells in control mice, but its expression level was clearly increased in DN mice. Importantly, nuclear accumulation of MRTFB in mesangial cells was obvious to detect especially in the region where mesangial expansion was taking place (Fig. 3d).

MRTF transcriptional target genes are upregulated in DN patients with type 2 diabetes

We investigated the expression of MRTF transcriptional target genes in an early DN cohort of Pima Indians from the Gila River Indian Community in Arizona, individuals with type 2 diabetes [31]. mRNA expression analysis of microdissected glomeruli showed that various genes significantly changed in mouse DN were also significantly regulated in early human DN (Fig. 4a). Many of these target genes were positively correlated with mesangial volume (Fig. 4b). Independent single-nucleus transcriptomics data from early human DN [3] showed that these MRTF transcription genes were mostly regulated in mesangial cells (Additional file 1: Figure  S7a). The MRTF transcriptional target genes were mapped in the ERCB dataset. These genes were significantly regulated in glomeruli in the context of DN and other diseases potentially causing mechanical stresses, such as arterial hypertension, as well as various other glomerular diseases with glomerulosclerosis, such as lupus nephritis (caused by systemic lupus erythematosus, SLE) and IgA nephropathy, but not in minimal change disease (MCD) (Fig. 4c).

Fig. 4
figure 4

MRTF transcriptional target genes are upregulated in DN patients with type 2 diabetes. a Heatmap showing the log2FC (DN vs. Ctrl) values of MRTF transcriptional target genes in microdissected glomeruli from Pima early DN patients. b Pearson correlations between MRTF target gene expression levels and mesangial volumes in Pima early DN patients. c Heatmap showing the log2FC (DN vs. Ctrl) values of MRTF transcriptional target genes in microdissected glomeruli from ERCB patient datasets. Nonsignificant genes are shown in gray. DN, diabetic nephropathy; HT, hypertensive nephropathy; IgA, IgA nephropathy; FSGS, focal segmental glomerulosclerosis; MGN, membranous nephropathy; SLE, lupus nephritis; RPGN, ANCA-associated glomerulonephritis; MCD, minimal change disease

Next, we investigated the expression and location of MRTFA and MRTFB in early human DN by immunofluorescent staining. A basal and comparable expression level of MRTFA was detected in all three glomerular cell types in both healthy and early DN kidney biopsy samples (Fig. 5). Consistent with the observations in mouse DN, MRTFB was absent in most mesangial cells from healthy individuals but was accumulated in the nuclei of mesangial cells where glomerular lesions were observed (Fig. 5). Taken together, our immunofluorescent staining analyses on human and mouse DN confirm the activation of MRTFB in mesangial cells during early DN.

Fig. 5
figure 5

Immunofluorescence staining confirms the activation of MRTFB in mesangial cells during early DN. Immunofluorescence staining of the human kidney samples (ctrl: n = 2, DN: n = 5). Representative images showing the expression and localization of MRTFA and MRTFB in control and DN glomeruli. Glomerular endothelial cells were labeled by CD34, podocytes were labeled by NPHS1 (Nephrin), and nuclei were counterstained with HOECHST. Arrows indicate mesangial cells, which were negative for CD34 and NPHS1

Mesangial cells exhibit dominant signaling networks in DN glomeruli

Our results suggested that MRTF-SRF transcriptional regulation was activated in both mouse and human DN, predominantly in mesangial cells. Additionally, we identified various mesangial cell marker genes that encode transmembrane proteins responsible for the perception and transduction of mechanical signals, including mechanosensitive ion channels (MSCs), GPCRs, integrins, and cadherins [39] (Fig. 6a). Furthermore, we performed CellChat [23] to study the changes of intercellular communications (secreted signaling, cell–cell contact, and ECM-receptor) taking place in DN. We compared interaction strength between all kidney cell types in DN versus control mice, which showed dominating interactions of mesangial cells with all other cell types in DN kidneys (Fig. 6b). Additionally, mesangial cells dominated the intercellular communications (Fig. 6c) and interacted as a sender extensively with gECs and podocytes (Fig. 6d) in DN glomeruli. These interactions involved various well-known growth factor signaling pathways associated with the pathogenesis of chronic kidney disease (CKD), such as bone morphogenetic protein (BMP), fibroblast growth factor (FGF), and vascular endothelial growth factor (VEGF), as well as interactions between multiple extracellular matrix proteins and integrins. Notably, the BMP2/4 signaling pathways identified in DN mesangial cell interactions contribute to glomerulosclerosis and tubulointerstitial fibrosis, particularly in DN [41].

Fig. 6
figure 6

Mesangial cells exhibit dominant signaling networks in DN glomeruli. a Dot plot showing mesangial marker genes encoding transmembrane proteins responsible for the perception and transduction of mechanical signals. MSCs, mechanosensitive ion channels; GPCRs, G protein-coupled receptors. b Change of cell–cell interactions in all pairs of kidney cell types in DN vs. Ctrl mice. c Change of cell–cell interactions in all pairs of glomerular cell types in DN vs. Ctrl mice. d Increased cell–cell interactions for mesangial cell signaling to podocytes and glomerular endothelial cell in DN mice

Kidney ex vivo perfusion activates mechanosensitive signaling pathways

Glomerular hyperfiltration as a consequence of increased single-nephron glomerular filtration rate may be a common upstream mechanism contributing to CKD [2]. In diabetes, glomerular hyperfiltration is driven by the sodium–glucose cotransporter 2 (SGLT2) [2]. Hyper-reabsorption of glucose and sodium in the PT ultimately leads to restricted tubuloglomerular feedback, resulting in increased intraglomerular pressure and glomerular hyperfiltration [2]. Consequently, circumferential and axial capillary wall stress, as well as fluid shear stress at the glomerular filtration barrier, is increased [42]. Since the mechanosensitive signaling pathway MRTF-SRF is a central coordinator of fibrosis-relevant mechanoresponse [39], we hypothesize that glomerular hyperfiltration may directly activate MRTF-SRF, leading to cellular changes associated with diabetic glomerulopathy.

We used normothermic machine perfusion to prove this hypothesis. The ex vivo kidney perfusion takes place for 1 to 2 h at 100 mmHg with 100% Gey’s Balanced Salt Solution (GBSS) without the addition of albumin or red blood cells (see the “Methods” section). Under this condition, tubular glomerular feedback is absent and 100 mmHg is a supraphysiologic pressure leading to glomerular hyperfiltration [43].

We generated bulk and single-nucleus transcriptomics data of isolated mouse glomeruli and pig renal cortex tissue, respectively (Fig. 7a, d), and gene expression was compared in perfused versus unperfused conditions. In ex vivo perfused mouse glomeruli, upstream analysis in DEGs suggested activation of MRTFA/B, SRF, and YAP1 (Fig. 7b). Enrichment analysis of significantly changed MRTF transcriptional target genes revealed actin cytoskeleton organization, VEGFA-VEGFR2 signaling pathway, smooth muscle contraction, PDGFRB pathway, and cellular response to external stimulus (Fig. 7c). In ex vivo perfused pig kidney, a total of 11,692 single nuclei were profiled from perfused and unperfused renal cortex tissue. Unsupervised clustering identified 18 clusters (Additional file 1: Figure S7b), which were classified into 12 cell types (Fig. 7e; Additional file 1: Figure S7c) on the basis of cell/nucleus-specific marker genes reported for previous human kidney single-nucleus data (Fig. 7f) [3, 36]. We identified the major glomerular cell types including a stromal cell type (indicated by “STROMA” in figures) defined by the nucleus marker gene MEIS1 and FHL2 (Fig. 7f) [44]. Immunofluorescent staining on the pig renal cortex showed that MEIS1-positive cells were mainly composed of mesangial cells, but also some interstitial cells (Fig. 7g). Upstream analysis on DEGs revealed the activation of MRTFA/B, SRF, and YAP1 in stromal cells, as well as the activation of YAP1 and MRTFB in podocytes and endothelial cells, respectively (Fig. 7h). Significantly regulated MRTF transcriptional target genes were found predominantly in stromal cells (Fig. 7i). These findings support our hypothesis that the activation of mechanosensitive signaling pathways is a common mechanism in response to glomerular hyperfiltration.

Fig. 7
figure 7

Kidney ex vivo perfusion activates mechanosensitive signaling pathways. a Experimental scheme of mouse kidney ex vivo perfusion. b Heatmap showing the log2FC values of MRTF transcriptional target genes in perfused versus control glomeruli. c Top enriched pathways of MRTF transcriptional target genes. d Experimental scheme of pig kidney ex vivo perfusion. e UMAP plot of annotated cell types from snRNA-seq of perfused and control pig kidney tissue. g Immunofluorescence staining against MEIS1 in the pig kidney cortex. Nuclei were counterstained with DAPI. G, glomerulus. Dashed lines indicate the area of glomerulus. Scale bar: 100 μm. f Dot plot displaying defining marker genes for each cell type. STROMA, stromal cells. h Heatmap showing z scores of mechanosensitive transcriptional regulators estimated by IPA. i Dot plot displaying MRTF transcriptional target genes significantly changed in perfused kidney tissue. Nonsignificant regulation or genes are shown in gray g, h

Discussion

BTBR ob/ob is a useful mouse model to investigate cellular changes and molecular mechanisms in early DN. Integrative analysis of single-cell and bulk RNA-seq data led to the identification of the mechanosensitive pathway MRTF-SRF in DN glomeruli. The association of MRTF-SRF activation with human DN was validated in independent patient cohorts at both mRNA and protein levels, and MRTF-SRF is suggested as a common mechanism in response to glomerular hyperfiltration.

Single-cell analysis on BTBR ob/ob kidneys suggested that PT cells and mesangial cells are the most sensitive cells to hyperglycemia and initiate glucose-dependent and glucose-independent downstream pathways in response to stress in the kidneys, respectively. PT cells showed adaptations to glucose metabolic fluxes and upregulation of antioxidant pathways. Increased glucose influx in cells fuels oxidative phosphorylation, leading to the overproduction of superoxide by mitochondria and the generation of reactive oxygen species, which play a central role in initiating diverse pathways responsible for diabetic abnormalities such as cellular dysfunction, inflammation, and fibrosis [1]. Mesangial cells showed the activation of other signaling pathways obviously independent of glucose metabolic flux. This finding is consistent with those of previous studies showing that overexpression of GLUT1 in mesangial cells exposed to high glucose fails to induce TGF-β1 synthesis and oxidative stress [45] and that transgenic mice with GLUT1 overexpression in mesangial cells do not develop pathological phenotypes until 26 weeks of age [46].

Mouse DN glomeruli exhibited shared features with human DN glomeruli. Approximately 70% glomerular DEGs identified in DN mice were also significantly regulated in DN patients. Importantly, MRTF-SRF transcriptional target genes were significantly regulated in both mouse and human DN glomeruli. Notably, mesangial cells responded at the earliest time point and showed a persistent activation of MRTFB compared to podocytes and gECs. Immunofluorescent staining further confirmed the activation of MRTFB in the region where mesangial expansion was taking place. Mechanic stimulation triggers cytoskeleton remodeling, resulting in the activation of mechanosensitive transcriptional pathway MRTF-SRF [39], which plays a key role in organ fibrosis [39, 47]. Upon stimulation, cytoplasmic MRTF translocate to the nucleus and interact with SRF to target the expression of muscle-specific and contractile genes, as well as extracellular matrix genes [48, 49]. Interestingly, single-cell data showed mesangial cell expression of numerous genes that encode transmembrane proteins responsible for the perception and transduction of mechanical signals. DN mesangial cells exhibited a dominant signaling network in the whole kidney, interacted with gECs and podocytes by ECM, and secreted signaling. These results suggest that MRTF-SRF transcriptional regulation in mesangial cells could be an important pathological pathway contributing to diabetic glomerulopathy.

Since the mechanical stress caused by glomerular hyperfiltration is thought to play a critical role in the development of glomerulopathy [42], we speculate that glomerular hyperfiltration, a key event in early DN [1, 2], directly activates mechanosensitive signaling pathways in renal glomeruli, and consequently contribute to diabetic glomerulopathy. This speculation is supported by evidence in our study. First, single-cell and bulk RNA-seq of DN mice showed activation of MRTFA/B in three major glomerular cell types, which was also observed in an early DN cohort. Second, ex vivo kidney perfusion-induced glomerular hyperfiltration showed a prominent activation of MRTFA/B, SRF, and YAP1 in mouse glomeruli. Third, snRNA-seq of perfused pig renal cortex tissue showed activation of these mechanosensitive transcriptional regulators, which were most prominently activated in stromal cells composing mesangial cells.

Many MRTF target genes were significantly regulated not only in DN but also in the context of other glomerular diseases associated with glomerulosclerosis. These findings suggest a common role of MRTF in diabetic glomerulopathy, as well as in glomerulosclerosis generally. Previous studies demonstrated that MRTF-SRF and YAP are essential for podocyte structure and function [50, 51]. Tubulointerstitial fibrosis was diminished in MRTFA-deficient mice and MRTFA was activated in high-fat diet- and streptozotocin-induced DN [52]. In addition, genetic deletion or pharmacological inhibition of MRTF attenuates fibrosis in various pathological conditions, including kidney diseases like diabetic nephropathy, obstructive nephropathy, acute kidney injury, and polycystic kidney disease [47]. Therefore, MRTF is important for glomerular cells and is a potential therapeutic target for DN, as well as glomerulosclerosis.

Finally, while the principle of SGLT2 inhibition has been proven to be the major advance in slowing down CKD, the mechanisms remain unclear [2]. It was proposed that a lowering of the intraglomerular pressure by SGLT2 inhibition might drive the renal protective effects. Our results unravel now mechanosensitive signaling pathways as possible mediators of SGLT2-dependent effects, opening the search for novel druggable downstream targets of SGLT2.

Conclusions

Our study presents a comprehensive single-cell transcriptomic landscape of early DN and reveals cell-specific alterations in gene expression that occur with DN onset and progression. Our intensive analysis of glomerular cell types reveals that mechanosensitive signaling is associated with diabetic glomerulopathy and could play a driving role in response to glomerular hyperfiltration. The MRTF transcription pathway may be part of a common mechanism in the glomerulus in the context of glomerulosclerosis.

Availability of data and materials

Animal data: The single-cell data of mouse kidneys are deposited in Gene Expression Omnibus database (GEO) under accession number GSE218563. The bulk RNA sequencing data of mice ex vivo kidney perfusion and glomerular cell types are available in GEO at GSE218086. The single-nucleus RNA sequencing data of pig ex vivo kidney perfusion are available in GEO at GSE218413. Additional interactive visualization of the mouse kidney single-cell dataset is available through the cellxgene portal located at http://3.139.221.155.

ERCB-KFB data: Gene expression data have been reported in previous studies [29, 53,54,55]. The CEL files of the different diseases and controls can be found at the GEO database (https://www.ncbi.nlm.nih.gov/geo/) under accession numbers GSE32591 [53], GSE35489 [54], GSE37463 [53], GSE47185 [55], GSE99340 [29].

Pima Indians cohort data: Due to privacy protection concerns, individual-level genotype and gene expression data from this study cannot be made publicly available.

Abbreviations

ATAC:

Assay for transposase-accessible chromatin

CAP:

Cold-active protease

CKD:

Chronic kidney disease

Ctrl:

Control

DEG:

Differentially expressed gene

DN:

Diabetic nephropathy

ECM:

Extracellular matrix

ERCB-KFB:

European Renal cDNA Bank-Kröner-Fresenius Biopsy Bank

FACS:

Fluorescence-activated cell sorting

FSGS:

Focal segmental glomerulosclerosis

GBM:

Glomerular basement membrane

GPCR:

G protein-coupled receptor

HT:

Hypertensive nephropathy

IPA:

Ingenuity Pathway Analysis

MCD:

Minimal change disease

MGN:

Membranous nephropathy

MSC:

Mechanosensitive ion channel

PCA:

Principal component analysis

QC:

Quality control

SCENIC:

Single-cell regulatory network inference and clustering

SGLT2:

Sodium–glucose cotransporter 2

SLE:

Systemic lupus erythematosus

UMAP:

Uniform Manifold Approximation and Projection

UMI:

Unique molecular identifier

WT:

Wildtype

References

  1. Thomas MC, Brownlee M, Susztak K, Sharma K, Jandeleit-Dahm KA, Zoungas S, Rossing P, Groop PH, Cooper ME. Diabetic kidney disease. Nat Rev Dis Primers. 2015;1:15018.

    Article  Google Scholar 

  2. Anders HJ, Huber TB, Isermann B, Schiffer M. CKD in diabetes: diabetic kidney disease versus nondiabetic kidney disease. Nat Rev Nephrol. 2018;14:361–77.

    Article  CAS  Google Scholar 

  3. Wilson PC, Wu H, Kirita Y, Uchimura K, Ledru N, Rennke HG, Welling PA, Waikar SS, Humphreys BD. The single-cell transcriptomic landscape of early human diabetic nephropathy. Proc Natl Acad Sci U S A. 2019;116:19619–25.

    Article  CAS  Google Scholar 

  4. Wilson PC, Muto Y, Wu H, Karihaloo A, Waikar SS, Humphreys BD. Multimodal single cell sequencing implicates chromatin accessibility and genetic background in diabetic kidney disease progression. Nat Commun. 2022;13:5253.

    Article  CAS  Google Scholar 

  5. Stefansson VTN, Nair V, Melsom T, Looker HC, Mariani LH, Fermin D, et al. Molecular programs associated with glomerular hyperfiltration in early diabetic kidney disease. Kidney Int. 2022;102(6):1345–58.

    Article  CAS  Google Scholar 

  6. Fu J, Akat KM, Sun Z, Zhang W, Schlondorff D, Liu Z, Tuschl T, Lee K, He JC. Single-cell RNA profiling of glomerular cells shows dynamic changes in experimental diabetic kidney disease. J Am Soc Nephrol. 2019;30:533–45.

    Article  CAS  Google Scholar 

  7. Wu J, Sun Z, Yang S, Fu J, Fan Y, Wang N, Hu J, Ma L, Peng C, Wang Z, et al. Kidney single-cell transcriptome profile reveals distinct response of proximal tubule cells to SGLT2i and ARB treatment in diabetic mice. Mol Ther. 2022;30:1741–53.

    Article  CAS  Google Scholar 

  8. Wu H, Villalobos RG, Yao X, Reilly D, Chen T, Rankin M, Myshkin E, Breyer MD, Humphreys BD. Mapping the single-cell transcriptomic response of murine diabetic kidney disease to therapies. Cell Metab. 2022;34(1064–1078):e1066.

    Google Scholar 

  9. Hudkins KL, Pichaiwong W, Wietecha T, Kowalewska J, Banas MC, Spencer MW, Muhlfeld A, Koelling M, Pippin JW, Shankland SJ, et al. BTBR Ob/Ob mutant mice model progressive diabetic nephropathy. J Am Soc Nephrol. 2010;21:1533–42.

    Article  CAS  Google Scholar 

  10. Ericsson A, Tonelius P, Lal M, Sabirsh A, Bottcher G, William-Olsson L, et al. The effects of dual PPARalpha/gamma agonism compared with ACE inhibition in the BTBRob/ob mouse model of diabetes and diabetic nephropathy. Physiol Rep. 2017;5(5):e13186.

  11. Clee SM, Nadler ST, Attie AD. Genetic and genomic studies of the BTBR ob/ob mouse model of type 2 diabetes. Am J Ther. 2005;12:491–8.

    Article  Google Scholar 

  12. Czogalla J, Schweda F, Loffing J. The mouse isolated perfused kidney technique. J Vis Exp 2016

  13. Breda PC, Wiech T, Meyer-Schwesinger C, Grahammer F, Huber T, Panzer U, Tiegs G, Neumann K. Renal proximal tubular epithelial cells exert immunomodulatory function by driving inflammatory CD4(+) T cell responses. Am J Physiol Renal Physiol. 2019;317:F77–89.

    Article  CAS  Google Scholar 

  14. Menon R, Otto EA, Kokoruda A, Zhou J, Zhang Z, Yoon E, et al. Single-cell analysis of progenitor cell dynamics and lineage specification in the human fetal kidney. Development. 2018;145(16):dev164038.

  15. Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM 3rd, Hao Y, Stoeckius M, Smibert P, Satija R. Comprehensive integration of single-cell data. Cell. 2019;177(1888–1902):e1821.

    Google Scholar 

  16. Park J, Shrestha R, Qiu C, Kondo A, Huang S, Werth M, Li M, Barasch J, Susztak K. Single-cell transcriptomics of the mouse kidney reveals potential cellular targets of kidney disease. Science. 2018;360:758–63.

    Article  CAS  Google Scholar 

  17. Finak G, McDavid A, Yajima M, Deng J, Gersuk V, Shalek AK, Slichter CK, Miller HW, McElrath MJ, Prlic M, et al. MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data. Genome Biol. 2015;16:278.

    Article  Google Scholar 

  18. Young MD, Behjati S. SoupX removes ambient RNA contamination from droplet-based single-cell RNA sequencing data. Gigascience. 2020;9(12):giaa151.

  19. Wolock SL, Lopez R, Klein AM. Scrublet: computational identification of cell doublets in single-cell transcriptomic data. Cell Syst. 2019;8(281–291):e289.

    Google Scholar 

  20. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.

    Article  CAS  Google Scholar 

  21. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.

    Article  CAS  Google Scholar 

  22. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  Google Scholar 

  23. Jin S, Guerrero-Juarez CF, Zhang L, Chang I, Ramos R, Kuan CH, Myung P, Plikus MV, Nie Q. Inference and analysis of cell-cell communication using Cell Chat. Nat Commun. 2021;12:1088.

    Article  CAS  Google Scholar 

  24. Aibar S, Gonzalez-Blas CB, Moerman T, Huynh-Thu VA, Imrichova H, Hulselmans G, Rambow F, Marine JC, Geurts P, Aerts J, et al. SCENIC: single-cell regulatory network inference and clustering. Nat Methods. 2017;14:1083–6.

    Article  CAS  Google Scholar 

  25. Van de Sande B, Flerin C, Davie K, De Waegeneer M, Hulselmans G, Aibar S, Seurinck R, Saelens W, Cannoodt R, Rouchon Q, et al. A scalable SCENIC workflow for single-cell gene regulatory network analysis. Nat Protoc. 2020;15:2247–76.

    Article  Google Scholar 

  26. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, Benner C, Chanda SK. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.

    Article  Google Scholar 

  27. Cohen CD, Frach K, Schlondorff D, Kretzler M. Quantitative gene expression analysis in renal biopsies: a novel protocol for a high-throughput multicenter application. Kidney Int. 2002;61:133–40.

    Article  CAS  Google Scholar 

  28. Neusser MA, Lindenmeyer MT, Kretzler M, Cohen CD. Genomic analysis in nephrology–towards systems biology and systematic medicine? Nephrol Ther. 2008;4:306–11.

    Article  CAS  Google Scholar 

  29. Shved N, Warsow G, Eichinger F, Hoogewijs D, Brandt S, Wild P, Kretzler M, Cohen CD, Lindenmeyer MT. Transcriptome-based network analysis reveals renal cell type-specific dysregulation of hypoxia-associated transcripts. Sci Rep. 2017;7:8576.

    Article  Google Scholar 

  30. Weil EJ, Fufaa G, Jones LI, Lovato T, Lemley KV, Hanson RL, Knowler WC, Bennett PH, Yee B, Myers BD, Nelson RG. Effect of losartan on prevention and progression of early diabetic nephropathy in American Indians with type 2 diabetes. Diabetes. 2013;62:3224–31.

    Article  CAS  Google Scholar 

  31. Nair V, Komorowsky CV, Weil EJ, Yee B, Hodgin J, Harder JL, Godfrey B, Ju W, Boustany-Kari CM, Schwarz M, et al. A molecular morphometric approach to diabetic kidney disease can link structure to function and outcome. Kidney Int. 2018;93:439–49.

    Article  CAS  Google Scholar 

  32. Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98:5116–21.

    Article  CAS  Google Scholar 

  33. Adam M, Potter AS, Potter SS. Psychrophilic proteases dramatically reduce single-cell RNA-seq artifacts: a molecular atlas of kidney development. Development. 2017;144:3625–32.

    CAS  Google Scholar 

  34. Barry DM, McMillan EA, Kunar B, Lis R, Zhang T, Lu T, Daniel E, Yokoyama M, Gomez-Salinero JM, Sureshbabu A, et al. Molecular determinants of nephron vascular specialization in the kidney. Nat Commun. 2019;10:5705.

    Article  CAS  Google Scholar 

  35. He B, Chen P, Zambrano S, Dabaghie D, Hu Y, Moller-Hackbarth K, Unnersjo-Jess D, Korkut GG, Charrin E, Jeansson M, et al. Single-cell RNA sequencing reveals the mesangial identity and species diversity of glomerular cell transcriptomes. Nat Commun. 2021;12:2141.

    Article  CAS  Google Scholar 

  36. Lake BB, Chen S, Hoshi M, Plongthongkum N, Salamon D, Knoten A, Vijayan A, Venkatesh R, Kim EH, Gao D, et al. A single-nucleus RNA-sequencing pipeline to decipher the molecular anatomy and pathophysiology of human kidneys. Nat Commun. 2019;10:2832.

    Article  Google Scholar 

  37. Betz B, Conway BR. An Update on the Use of Animal Models in Diabetic Nephropathy Research. Curr Diab Rep. 2016;16:18.

    Article  Google Scholar 

  38. Foster CT, Gualdrini F, Treisman R. Mutual dependence of the MRTF-SRF and YAP-TEAD pathways in cancer-associated fibroblasts is indirect and mediated by cytoskeletal dynamics. Genes Dev. 2017;31:2361–75.

    Article  CAS  Google Scholar 

  39. Tschumperlin DJ, Ligresti G, Hilscher MB, Shah VH. Mechanosensing and fibrosis. J Clin Invest. 2018;128:74–84.

    Article  Google Scholar 

  40. Esnault C, Stewart A, Gualdrini F, East P, Horswell S, Matthews N, Treisman R. Rho-actin signaling to the MRTF coactivators dominates the immediate transcriptional response to serum in fibroblasts. Genes Dev. 2014;28:943–58.

    Article  CAS  Google Scholar 

  41. Perera N, Ritchie RH, Tate M. The Role of Bone Morphogenetic Proteins in Diabetic Complications. ACS Pharmacol Transl Sci. 2020;3:11–20.

    Article  CAS  Google Scholar 

  42. Kriz W, Lemley KV. A potential role for mechanical forces in the detachment of podocytes and the progression of CKD. J Am Soc Nephrol. 2015;26:258–69.

    Article  Google Scholar 

  43. Schurek HJ, Neumann KH, Schweda F, Czogalla J In A Laboratory Manual of Kidney Perfusion Techniques. Munster Germany. 2017

  44. Avraham S, Korin B, Chung JJ, Oxburgh L, Shaw AS. The mesangial cell - the glomerular stromal cell. Nat Rev Nephrol. 2021;17(12):855–64.

    Article  Google Scholar 

  45. Weigert C, Brodbeck K, Brosius FC 3rd, Huber M, Lehmann R, Friess U, Facchin S, Aulwurm S, Haring HU, Schleicher ED, Heilig CW. Evidence for a novel TGF-beta1-independent mechanism of fibronectin production in mesangial cells overexpressing glucose transporters. Diabetes. 2003;52:527–35.

    Article  CAS  Google Scholar 

  46. Wang Y, Heilig K, Saunders T, Minto A, Deb DK, Chang A, Brosius F, Monteiro C, Heilig CW. Transgenic overexpression of GLUT1 in mouse glomeruli produces renal disease resembling diabetic glomerulosclerosis. Am J Physiol Renal Physiol. 2010;299:F99–111.

    Article  CAS  Google Scholar 

  47. Miranda MZ, Lichner Z, Szaszi K, Kapus A. MRTF: basic biology and role in kidney disease. Int J Mol Sci. 2021;22(11):6040.

    Article  CAS  Google Scholar 

  48. Olson EN, Nordheim A. Linking actin dynamics and gene transcription to drive cellular motile functions. Nat Rev Mol Cell Biol. 2010;11:353–65.

    Article  CAS  Google Scholar 

  49. Gau D, Roy P. SRF'ing and SAP'ing - the role of MRTF proteins in cell migration. J Cell Sci. 2018;131(19):jcs218222.

  50. Schwartzman M, Reginensi A, Wong JS, Basgen JM, Meliambro K, Nicholas SB, D’Agati V, McNeill H, Campbell KN. Podocyte-specific deletion of yes-associated protein causes FSGS and progressive renal failure. J Am Soc Nephrol. 2016;27:216–26.

    Article  CAS  Google Scholar 

  51. Guo B, Lyu Q, Slivano OJ, Dirkx R, Christie CK, Czyzyk J, Hezel AF, Gharavi AG, Small EM, Miano JM. Serum response factor is essential for maintenance of podocyte structure and function. J Am Soc Nephrol. 2018;29:416–22.

    Article  CAS  Google Scholar 

  52. Xu H, Wu X, Qin H, Tian W, Chen J, Sun L, Fang M, Xu Y. Myocardin-related transcription factor A epigenetically regulates renal fibrosis in diabetic nephropathy. J Am Soc Nephrol. 2015;26:1648–60.

    Article  CAS  Google Scholar 

  53. Berthier CC, Bethunaickan R, Gonzalez-Rivera T, Nair V, Ramanujam M, Zhang W, Bottinger EP, Segerer S, Lindenmeyer M, Cohen CD, et al. Cross-species transcriptional network analysis defines shared inflammatory responses in murine and human lupus nephritis. J Immunol. 2012;189:988–1001.

    Article  CAS  Google Scholar 

  54. Reich HN, Tritchler D, Cattran DC, Herzenberg AM, Eichinger F, Boucherot A, Henger A, Berthier CC, Nair V, Cohen CD, et al. A molecular signature of proteinuria in glomerulonephritis. PLoS ONE. 2010;5:e13451.

    Article  Google Scholar 

  55. Ju W, Greene CS, Eichinger F, Nair V, Hodgin JB, Bitzer M, Lee YS, Zhu Q, Kehata M, Li M, et al. Defining cell-type specificity at the transcriptional level in human disease. Genome Res. 2013;23:1862–73.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors wish to thank Alina Borchers, Anja Obser, Silvia Chilla, and Valerie Oberueber, as well as FACS Sorting Core Unit and Microscopy Imaging Facility at University Medical Center Hamburg-Eppendorf, Germany, for the excellent technical support.

Funding

TBH was supported by the DFG (CRC1192, HU 1016/8–2, HU 1016/11–1, HU 1016/ 12–1), by the BMBF (STOP-FSGS- 01GM1901C, NephrESA-031L0191E and DEFEAT PANDEMIcs), by the Else-Kröner Fresenius Foundation (Else Kröner-Promotionskolleg –iPRIME), and by the European Research Council-ERC (grant 616891). This work was further supported by Unicyte AG. This project has received funding from Abbvie and the Innovative Medicines Initiative 2 Joint undertaking under grant agreement No. 115974 (Beat-DKD). This Joint undertaking receives support from the European Union’s Horizon 2020 research and innovation program and EFPIA with JDRF. The work performed at the University of Michigan was supported by George M. O’Brien Michigan Kidney Translational Core Center funded by NIH/NIDDK grant 2P30-DK-081943. The work performed in Phoenix was supported by the Intramural Research Program of the National Institute of Diabetes and Digestive and Kidney Disease and the American Diabetes Association (Clinical Science Award 1–08-CR-42). The ERCB-KFB was supported by the Else Kröner-Fresenius Foundation. We also thank all participating centers of the European Renal cDNA Bank—Kröner-Fresenius biopsy bank (ERCB-KFB) and their patients for their cooperation. Active members at the time of the study are found in N. Shved et al., Scientific reports 7, 8576 (2017). This work was partially funded by the Else Kröner-Fresenius-Stiftung and the Eva Luise und Horst Köhler Stiftung – Project No: 2019_KollegSE.04. Open Access funding enabled and organized by Projekt DEAL.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: SL and TBH. Formal analysis: SL, YZ, and SLU. Investigation: SL, YZ, SLU, TZ, MTL, VN, SEG, GW, JC, HA, SZ, ZL, MS, DD, TW, and CMS. Method development: SL, SLU, CMS, and CK. Resources: MTL, VN, RGN, DF, TW, FG, and MK. Data and sample curation: SL, YZ, SLU, MTL, VN, and SB. Writing (original draft): SL and TBH. Visualization: SL, YZ, and TZ. Supervision and funding acquisition: TBH. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Shuya Liu or Tobias B. Huber.

Ethics declarations

Ethics approval and consent to participate

The study of ERCB-KFB was carried out in accordance with the Declaration of Helsinki. The ERCB-KFB represents a multi-center study capturing renal biopsy tissue for gene expression analysis. Biospecimens were obtained after written consent and approval of the local ethics committees and in the frame of the European Renal cDNA Bank approved by the Specialized Ethics Subcommittee of Internal Medicine of the University Hospital Zurich, Cantonal Ethics Committee Zurich and Ethics Committee of the Medical Faculty of the University of Munich (LMU). The authors declare that the study of Pima Indians cohort was conducted according to the regulations established by the Clinical Research and Ethics Committee and conforms to the principles of the Helsinki Declaration of the World Medical Association. All studies of the diabetic kidney disease population were approved by the Institutional Review Board of the National Institute of Diabetes and Digestive and Kidney Diseases. All participants provided informed consent. The renal biopsy studies were carried out in accordance with the Declaration of Helsinki and were performed in accordance with the law of Hamburg Hospitals (Hamburger Krankenhausgesetz, §12).

All animal experiments were conducted according to the National Institutes of Health Guide for the Care and Use of Laboratory Animals as well as German- and Hamburg law and approved by the veterinary administration of the City of Hamburg (Single-cell analysis license: Ü003-2018; Kidney normothermic machine perfusion license: N002/2020).

Consent for publication

Not applicable.

Competing interests

DD is an employee of Boehringer Ingelheim Pharma GmbH & Co. KG. T.B. Huber reports having consultancy agreements with AstraZeneca, Bayer, Boehringer-Ingelheim, DaVita, Fresenius Medical Care, Novartis, and Retrophin; receiving research funding from Amicus Therapeutics, Fresenius Medical Care; and being on the editorial board of Kidney International and the advisory board of Nature Review Nephrology. The other authors declare that they have no competing interests.

Additional information

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Figure S1. Characterization of BTBR ob/ob podocyte-reporter mice. Figure S2. Overview of the single cells from DN and control kidney samples. Figure S3. Distributions of cells across conditions and individual kidney samples. Figure S4. FACS analysis of podocytes (podo), glomerular endothelial cells (gECs) and mesangial cells (mesan) from the isolated mouse glomeruli. Figure 5S. Bulk RNA-seq results for FACS-sorted glomerular cells. Figure S6. Transcriptional regulations in DN mice. Figure S7. snRNA-seq of early human diabetic nephropathy and ex-vivo-perfused pig kidney tissue. Supplementary methods: FACS-sorted glomerular cells. Urine and serum analysis. Histological and ultrastructural analysis.

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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liu, S., Zhao, Y., Lu, S. et al. Single-cell transcriptomics reveals a mechanosensitive injury signaling pathway in early diabetic nephropathy. Genome Med 15, 2 (2023). https://doi.org/10.1186/s13073-022-01145-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13073-022-01145-4