A Single-Cell Transcriptome Atlas of the Mouse Glomerulus

Glomerular filtration rate and sexual dimorphism: lessons from animal and human studies

PURPOSE OF REVIEW

Integrating sex-based analyses is becoming a key point in the new recommendations, particularly in nephrology.

RECENT FINDINGS

Whereas single nephron glomerular filtration rate (GFR) is not different between men and women, male sex is associated, after multiple adjustments, with a higher number of nephrons. However, after indexation to body surface area, measured GFR (mGFR) in healthy potential kidney donors is not different between men and women between 20 and 50 years of age. After 50 years, mGFR decline seems faster in women than in men, which is concordant with the protective role of estrogens on renal function, as demonstrated in animal and some human studies. Conversely, testosterone has a detrimental effect on renal function. Of note, although testosterone has been shown to increase the kidney volume of the remnant kidney after a unilateral nephrectomy in animal models, this may generate deleterious hyperfiltration in the longer term.

SUMMARY

Taken together, these data highlight the impact of sex on GFR, notably through sexual hormones whose receptors are expressed in glomerular cells.

Upregulation of Piezo2 and increased extracellular matrix protein in diabetic kidney disease mice

Mechanical forces such as glomerular hyperfiltration are crucial in the pathogenesis and progression of diabetic kidney disease. Piezo2 is a mechanosensitive cation channel and plays a major role in various biological and pathophysiological phenomena. We previously reported Piezo2 expression in mouse and rat kidneys and its alteration by dehydration and hypertension. To elucidate the alteration of Piezo2 expression and its consequences in a mouse model of diabetic kidney disease, we used high salt-fed male KK-Ay mice, an accelerated genetic model of diabetic kidney disease. KK-Ay mice exhibited marked obesity, hyperglycemia, elevated blood pressure, higher creatinine clearance, and overt albuminuria. Histopathological analysis revealed glomerular hypertrophy, mesangial expansion, macrophage infiltration, tubular vacuolization, and interstitial fibrosis. The mRNA and protein expression analyses revealed (1) increased fibronectin protein expression in the glomerular areas, (2) upregulated Piezo2 expression in the glomerular mesangial cells and interstitial region, (3) increased Piezo2 and the fibronectin-coding gene Fn1 mRNA, and (4) a strong correlation of Piezo2 expression with that of Fn1 in the kidneys of diabetic kidney disease mice. Piezo2 upregulation and fibronectin accumulation were alleviated by an angiotensin II receptor blocker. In accordance with these in vivo results, in vitro study demonstrated that Piezo2 overexpression increased fibronectin production in HEK293T cells. In conclusion, we demonstrated that Piezo2 expression was upregulated in glomerular mesangial cells in a mouse model of diabetic kidney disease. Our results suggest that Piezo2 contributes to the progression of diabetic kidney disease by mediating glomerular fibronectin production, leading to glomerulosclerosis. Hyperfiltration is crucial in the pathogenesis of diabetic kidney disease. We showed that Piezo2 expression is upregulated in mesangial cells of diabetic kidney disease mice with glomerular fibronectin accumulation. Piezo2 overexpression increased fibronectin production in HEK293T cells. Piezo2 may contribute to diabetic kidney disease progression by mediating glomerular fibronectin production.

Single-cell atlas of human pancreatic islet and acinar endothelial cells in health and diabetes

Characterization of the vascular heterogeneity within the pancreas has previously been lacking. Here, we develop strategies to enrich islet-specific endothelial cells (ISECs) and acinar-specific endothelial cells (ASECs) from three human pancreases and corroborate these findings with three published pancreatic datasets. Single-cell RNA sequencing reveals the unique molecular signatures of ISECs, including structural genes COL13A1, ESM1, PLVAP, UNC5B, and LAMA4, angiocrine genes KDR, THBS1, BMPs and CXCR4, and metabolic genes ACE, PASK and F2RL3. ASECs display distinct signatures including GPIHBP1, CCL14, CD74, AQP1, KLF4, and KLF2, which may manage the inflammatory and metabolic needs of the exocrine pancreas. Ligand-receptor analysis suggests ISECs and ASECs interact with LUM+ fibroblasts and RGS5+ pericytes and smooth muscle cells via VEGF-A:VEGFR2, CXCL12:CXCR4, and LIF:LIFR pathways. Comparative expression and immunohistochemistry indicate disruption of endothelial-expressed CD74, ESM1, PLVAP, THBD, VWA1, and VEGF-A cross-talk among vascular and other cell types in diabetes. Thus, our data provide a single-cell vascular atlas of human pancreas, enabling deeper understanding of pancreatic pathophysiology in health and disease.

TRAIL induces podocyte PANoptosis via death receptor 5 in diabetic kidney disease

Podocytes can undergo PANoptosis (apoptosis, pyroptosis, and necroptosis). Diabetic kidney disease (DKD) is the leading cause of kidney failure, and podocyte loss is a major event leading to the progression of DKD. Here, we compared single cell RNA sequencing (scRNA-seq) data between three normal and three DKD human kidney samples and found a significant increase of TNFSF10 and TNFRSF10B expression in podocytes of patients with DKD. Tumor necrosis factor (TNF)-related apoptosis-inducing ligand (TRAIL), coded by TNFSF10, belongs to the TNF superfamily members and TNFRSF10B codes for death receptor 5 (DR5). We confirmed that expression of TRAIL and DR5 increased in podocytes of patients with DKD and correlated with the severity of DKD. In vitro, TNF-α stimulated TRAIL and DR5 expression in cultured human podocytes. Silence of TRAIL or DR5 by small interfering RNA alleviated TNF-α-stimulated podocytes PANoptosis, while overexpression of TRAIL, treatment with recombinant human TRAIL (rh-TRAIL) or the DR5 activator (Bioymifi) enhanced podocytes PANoptosis. In vivo, podocyte-specific deletion of TNFSF10 or TNFRSF10B alleviated podocyte and glomerular injury in high fat diet and streptozotocin-induced obese diabetic mice and was associated with decreased podocyte PANoptosis. Conversely, the induction of TNFSF10 overexpression specifically in podocytes exacerbated albuminuria and kidney injury in diabetic mice with increased podocyte PANoptosis. Additionally, administration of soluble DR5-Fc, an inhibitor of DR5, resulted in a marked reduction in albuminuria and glomerular injury in BTBR ob/ob mice. Our findings suggest a critical autocrine role of TRAIL/DR5 in inducing podocyte injury in DKD via activation of PANoptosis.

Single-cell RNA sequencing in diabetic kidney disease: a literature review

Diabetic kidney disease (DKD) is the leading cause of end-stage renal disease (ESRD), and its pathogenesis has not been clarified. Current research suggests that DKD involves multiple cell types and extra-renal factors, and it is particularly important to clarify the pathogenesis and identify new therapeutic targets. Single-cell RNA sequencing (scRNA-seq) technology is high-throughput sequencing of the transcriptomes of individual cells at the single-cell level, which is an effective technology for exploring the development of diseases by comparing genetic information, reflecting the differences in genetic information between cells, and identifying different cell subpopulations. Accumulating evidence supports the role of scRNA-seq in revealing the pathogenesis of diabetes and strengthening our understanding of the molecular mechanisms of DKD. We reviewed the scRNA-seq data this time. Then, we analyzed and discussed the applications of scRNA-seq technology in DKD research, including annotation of cell types, identification of novel cell types (or subtypes), identification of intercellular communication, analysis of cell differentiation trajectories, gene expression detection, and analysis of gene regulatory networks, and lastly, we explored the future perspectives of scRNA-seq technology in DKD research.

Loss of 3-O-sulfotransferase enzymes, Hs3st3a1 and Hs3st3b1, reduces kidney and glomerular size and disrupts glomerular architecture

Heparan sulfate (HS) is an important component of the kidney anionic filtration barrier, the glomerular basement membrane (GBM). HS chains attached to proteoglycan protein cores are modified by sulfotransferases in a highly ordered series of biosynthetic steps resulting in immense structural diversity due to negatively charged sulfate modifications. 3-O-sulfation is the least abundant modification generated by a family of seven isoforms but creates the most highly sulfated HS domains. We analyzed the kidney phenotypes in the Hs3st3a1, Hs3st3b1 and Hs3st6 -knockout (KO) mice, the isoforms enriched in kidney podocytes. Individual KO mice show no overt kidney phenotype, although Hs3st3b1 kidneys were smaller than wildtype (WT). Furthermore, Hs3st3a1-/-; Hs3st3b1-/- double knockout (DKO) kidneys were smaller but also had a reduction in glomerular size relative to wildtype (WT). Mass spectrometry analysis of kidney HS showed reduced 3-O-sulfation in Hs3st3a1-/- and Hs3st3b1-/-, but not in Hs3st6-/- kidneys. Glomerular HS showed reduced HS staining and reduced ligand-and-carbohydrate engagement (LACE) assay, a tool that detects changes in binding of growth factor receptor-ligand complexes to HS. Interestingly, DKO mice have increased levels of blood urea nitrogen, although no differences were detected in urinary levels of albumin, creatinine and nephrin. Finally, transmission electron microscopy showed irregular and thickened GBM and podocyte foot process effacement in the DKO compared to WT. Together, our data suggest that loss of 3-O-HS domains disrupts the kidney glomerular architecture without affecting the glomerular filtration barrier and overall kidney function.

Single-cell transcriptomics identifies aberrant glomerular angiogenic signalling in the early stages of WT1 kidney disease

WT1 encodes a podocyte transcription factor whose variants can cause an untreatable glomerular disease in early childhood. Although WT1 regulates many podocyte genes, it is poorly understood which of them are initiators in disease and how they subsequently influence other cell-types in the glomerulus. We hypothesised that this could be resolved using single-cell RNA sequencing (scRNA-seq) and ligand-receptor analysis to profile glomerular cell-cell communication during the early stages of disease in mice harbouring an orthologous human mutation in WT1 (Wt1R394W/+). Podocytes were the most dysregulated cell-type in the early stages of Wt1R394W/+ disease, with disrupted angiogenic signalling between podocytes and the endothelium, including the significant downregulation of transcripts for the vascular factors Vegfa and Nrp1. These signalling changes preceded glomerular endothelial cell loss in advancing disease, a feature also observed in biopsy samples from human WT1 glomerulopathies. Addition of conditioned medium from murine Wt1R394W/+ primary podocytes to wild-type glomerular endothelial cells resulted in impaired endothelial looping and reduced vascular complexity. Despite the loss of key angiogenic molecules in Wt1R394W/+ podocytes, the pro-vascular molecule adrenomedullin was upregulated in Wt1R394W/+ podocytes and plasma and its further administration was able to rescue the impaired looping observed when glomerular endothelium was exposed to Wt1R394W/+ podocyte medium. In comparative analyses, adrenomedullin upregulation was part of a common injury signature across multiple murine and human glomerular disease datasets, whilst other gene changes were unique to WT1 disease. Collectively, our study describes a novel role for altered angiogenic signalling in the initiation of WT1 glomerulopathy. We also identify adrenomedullin as a proangiogenic factor, which despite being upregulated in early injury, offers an insufficient protective response due to the wider milieu of dampened vascular signalling that results in endothelial cell loss in later disease. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Piezo ion channels: long-sought-after mechanosensors mediating hypertension and hypertensive nephropathy

Recent advances in mechanobiology and the discovery of mechanosensitive ion channels have opened a new era of research on hypertension and related diseases. Piezo1 and Piezo2, first reported in 2010, are regarded as bona fide mechanochannels that mediate various biological and pathophysiological phenomena in multiple tissues and organs. For example, Piezo channels have pivotal roles in blood pressure control, triggering shear stress-induced nitric oxide synthesis and vasodilation, regulating baroreflex in the carotid sinus and aorta, and releasing renin from renal juxtaglomerular cells. Herein, we provide an overview of recent literature on the roles of Piezo channels in the pathogenesis of hypertension and related kidney damage, including our experimental data on the involvement of Piezo1 in podocyte injury and that of Piezo2 in renin expression and renal fibrosis in animal models of hypertensive nephropathy. The mechanosensitive ion channels Piezo1 and Piezo2 play various roles in the pathogenesis of systemic hypertension by acting on vascular endothelial cells, baroreceptors in the carotid artery and aorta, and the juxtaglomerular apparatus. Piezo channels also contribute to hypertensive nephropathy by acting on mesangial cells, podocytes, and perivascular mesenchymal cells.

Lessons from Klotho mouse models to understand mineral homeostasis

AIM

Klotho, a key component of the endocrine fibroblast growth factor receptor-fibroblast growth factor axis, is a multi-functional protein that impacts renal electrolyte handling. The physiological significance of Klotho will be highlighted in the regulation of calcium, phosphate, and potassium metabolism.

METHODS

In this review, we compare several murine models with different renal targeted deletions of Klotho and the insights into the molecular and physiological function that these models offer.

RESULTS

In vivo, Klotho deficiency is associated with severely impaired mineral metabolism, with consequences on growth, longevity and disease development. Additionally, we explore the perspectives of Klotho in renal pathology and vascular events, as well as potential Klotho treatment options.

CONCLUSION

This comprehensive review emphasizes the use of Klotho to shed light on deciphering the renal molecular in vivo mechanisms in electrolyte handling, as well as novel therapeutic interventions.

Generating kidney organoids based on developmental nephrology

Over the past decade, the induction protocols for the two types of kidney organoids (nephron organoids and ureteric bud organoids) from pluripotent stem cells (PSCs) have been established based on the knowledge gained in developmental nephrology. Kidney organoids are now used for disease modeling and drug screening, but they also have potential as tools for clinical transplantation therapy. One of the options to achieve this goal would be to assemble multiple renal progenitor cells (nephron progenitor, ureteric bud, stromal progenitor) to reproduce the organotypic kidney structure from PSCs. At least from mouse PSCs, all the three progenitors have been induced and assembled into such "higher order" kidney organoids. We will provide an overview of the developmental nephrology required for the induction of renal progenitors and discuss recent advances and remaining challenges of kidney organoids for clinical transplantation therapy.

Identification of novel therapeutic targets for chronic kidney disease and kidney function by integrating multi-omics proteome with transcriptome

BACKGROUND

Chronic kidney disease (CKD) is a progressive disease for which there is no effective cure. We aimed to identify potential drug targets for CKD and kidney function by integrating plasma proteome and transcriptome.

METHODS

We designed a comprehensive analysis pipeline involving two-sample Mendelian randomization (MR) (for proteins), summary-based MR (SMR) (for mRNA), and colocalization (for coding genes) to identify potential multi-omics biomarkers for CKD and combined the protein-protein interaction, Gene Ontology (GO), and single-cell annotation to explore the potential biological roles. The outcomes included CKD, extensive kidney function phenotypes, and different CKD clinical types (IgA nephropathy, chronic glomerulonephritis, chronic tubulointerstitial nephritis, membranous nephropathy, nephrotic syndrome, and diabetic nephropathy).

RESULTS

Leveraging pQTLs of 3032 proteins from 3 large-scale GWASs and corresponding blood- and tissue-specific eQTLs, we identified 32 proteins associated with CKD, which were validated across diverse CKD datasets, kidney function indicators, and clinical types. Notably, 12 proteins with prior MR support, including fibroblast growth factor 5 (FGF5), isopentenyl-diphosphate delta-isomerase 2 (IDI2), inhibin beta C chain (INHBC), butyrophilin subfamily 3 member A2 (BTN3A2), BTN3A3, uromodulin (UMOD), complement component 4A (C4a), C4b, centrosomal protein of 170 kDa (CEP170), serologically defined colon cancer antigen 8 (SDCCAG8), MHC class I polypeptide-related sequence B (MICB), and liver-expressed antimicrobial peptide 2 (LEAP2), were confirmed. To our knowledge, 20 novel causal proteins have not been previously reported. Five novel proteins, namely, GCKR (OR 1.17, 95% CI 1.10-1.24), IGFBP-5 (OR 0.43, 95% CI 0.29-0.62), sRAGE (OR 1.14, 95% CI 1.07-1.22), GNPTG (OR 0.90, 95% CI 0.86-0.95), and YOD1 (OR 1.39, 95% CI 1.18-1.64,) passed the MR, SMR, and colocalization analysis. The other 15 proteins were also candidate targets (GATM, AIF1L, DQA2, PFKFB2, NFATC1, activin AC, Apo A-IV, MFAP4, DJC10, C2CD2L, TCEA2, HLA-E, PLD3, AIF1, and GMPR1). These proteins interact with each other, and their coding genes were mainly enrichment in immunity-related pathways or presented specificity across tissues, kidney-related tissue cells, and kidney single cells.

CONCLUSIONS

Our integrated analysis of plasma proteome and transcriptome data identifies 32 potential therapeutic targets for CKD, kidney function, and specific CKD clinical types, offering potential targets for the development of novel immunotherapies, combination therapies, or targeted interventions.

Mechanosensing by Piezo1 and its implications in the kidney

Piezo1 is an essential mechanosensitive transduction ion channel in mammals. Its unique structure makes it capable of converting mechanical cues into electrical and biological signals, modulating biological and (patho)physiological processes in a wide variety of cells. There is increasing evidence demonstrating that the piezo1 channel plays a vital role in renal physiology and disease conditions. This review summarizes the current evidence on the structure and properties of Piezo1, gating modulation, and pharmacological characteristics, with special focus on the distribution and (patho)physiological significance of Piezo1 in the kidney, which may provide insights into potential treatment targets for renal diseases involving this ion channel.

Macrophages communicate with mesangial cells through the CXCL12/DPP4 axis in lupus nephritis pathogenesis

Lupus nephritis (LN) occurs in 50% of cases of systemic lupus erythematosus (SLE) and is one of the most serious complications that can occur during lupus progression. Mesangial cells (MCs) are intrinsic cells in the kidney that can regulate capillary blood flow, phagocytose apoptotic cells, and secrete vasoactive substances and growth factors. Previous studies have shown that various types of inflammatory cells can activate MCs for hyperproliferation, leading to disruption of the filtration barrier and impairment of renal function in LN. Here, we characterized the heterogeneity of kidney cells of LN mice by single-nucleus RNA sequencing (snRNA-seq) and revealed the interaction between macrophages and MCs through the CXC motif chemokine ligand 12 (CXCL12)/dipeptidyl peptidase 4 (DPP4) axis. In culture, macrophages modulated the proliferation and migration of MCs through this ligand-receptor interaction. In LN mice, treatment with linagliptin, a DPP4 inhibitor, effectively inhibited MC proliferation and reduced urinary protein levels. Together, our findings indicated that targeting the CXCL12/DPP4 axis with linagliptin treatment may serve as a novel strategy for the treatment of LN via the CXCL12/DPP4 axis.

Single-Cell RNA Sequencing Identifies Response of Renal Lymphatic Endothelial Cells to Acute Kidney Injury

SIGNIFICANCE STATEMENT

The renal lymphatic vasculature and the lymphatic endothelial cells that make up this network play important immunomodulatory roles during inflammation. How lymphatics respond to AKI may affect AKI outcomes. The authors used single-cell RNA sequencing to characterize mouse renal lymphatic endothelial cells in quiescent and cisplatin-injured kidneys. Lymphatic endothelial cell gene expression changes were confirmed in ischemia-reperfusion injury and in cultured lymphatic endothelial cells, validating renal lymphatic endothelial cells single-cell RNA sequencing data. This study is the first to describe renal lymphatic endothelial cell heterogeneity and uncovers molecular pathways demonstrating lymphatic endothelial cells regulate the local immune response to AKI. These findings provide insights into previously unidentified molecular pathways for lymphatic endothelial cells and roles that may serve as potential therapeutic targets in limiting the progression of AKI.

BACKGROUND

The inflammatory response to AKI likely dictates future kidney health. Lymphatic vessels are responsible for maintaining tissue homeostasis through transport and immunomodulatory roles. Owing to the relative sparsity of lymphatic endothelial cells in the kidney, past sequencing efforts have not characterized these cells and their response to AKI.

METHODS

Here, we characterized murine renal lymphatic endothelial cell subpopulations by single-cell RNA sequencing and investigated their changes in cisplatin AKI 72 hours postinjury. Data were processed using the Seurat package. We validated our findings by quantitative PCR in lymphatic endothelial cells isolated from both cisplatin-injured and ischemia-reperfusion injury, by immunofluorescence, and confirmation in in vitro human lymphatic endothelial cells.

RESULTS

We have identified renal lymphatic endothelial cells and their lymphatic vascular roles that have yet to be characterized in previous studies. We report unique gene changes mapped across control and cisplatin-injured conditions. After AKI, renal lymphatic endothelial cells alter genes involved in endothelial cell apoptosis and vasculogenic processes as well as immunoregulatory signaling and metabolism. Differences between injury models were also identified with renal lymphatic endothelial cells further demonstrating changed gene expression between cisplatin and ischemia-reperfusion injury models, indicating the renal lymphatic endothelial cell response is both specific to where they lie in the lymphatic vasculature and the kidney injury type.

CONCLUSIONS

In this study, we uncover lymphatic vessel structural features of captured populations and injury-induced genetic changes. We further determine that lymphatic endothelial cell gene expression is altered between injury models. How lymphatic endothelial cells respond to AKI may therefore be key in regulating future kidney disease progression.

Sepsis induces heterogeneous transcription of coagulation- and inflammation-associated genes in renal microvasculature

BACKGROUND

Acute kidney injury (AKI) in sepsis patients increases patient mortality. Endothelial cells are important players in the pathophysiology of sepsis-associated AKI (SA-AKI), yet knowledge regarding their spatiotemporal involvement in coagulation disbalance and leukocyte recruitment is lacking. This study investigated the identity and kinetics of responses of different microvascular compartments in kidney cortex in response to SA-AKI.

METHODS

Laser microdissected arterioles, glomeruli, peritubular capillaries, and postcapillary venules from kidneys of mice subjected to cecal ligation and puncture (CLP) were analyzed using RNA sequencing. Differential expression and pathway enrichment analyses identified genes involved in coagulation and inflammation. A selection of these genes was evaluated by RT-qPCR in microvascular compartments of renal biopsies from patients with SA-AKI. The role of two identified genes in lipopolysaccharide-induced endothelial coagulation and inflammatory activation were determined in vitro in HUVEC using siRNA-based gene silencing.

RESULTS

CLP-sepsis in mice induced altered expression of approximately 400 genes in the renal microvasculature, with microvascular compartments exhibiting unique spatiotemporal responses. In mice, changes in gene expression related to coagulation and inflammation were most extensive in glomeruli at early and intermediate time points, with high induction of Plat, Serpine1, Thbd, Icam1, Stat3, and Ifitm3. In human SA-AKI, PROCR and STAT3 were induced in postcapillary venules, while SERPINE1 expression was diminished. IFITM3 was increased in arterioles and glomeruli. In vitro studies revealed that STAT3 and IFITM3 partly control endothelial coagulation and inflammatory activation.

CONCLUSION

Renal microvascular compartments in mice and humans exhibited heterogeneous changes in coagulation- and inflammation-related gene expression in response to SA-AKI. Additional research should aim at understanding the functional consequences of the here described heterogeneous microvascular responses to establish the usefulness of identified genes as therapeutic targets in SA-AKI.

Neuronally differentiated macula densa cells regulate tissue remodeling and regeneration in the kidney

Tissue regeneration is limited in several organs, including the kidney, contributing to the high prevalence of kidney disease globally. However, evolutionary and physiological adaptive responses and the presence of renal progenitor cells suggest an existing remodeling capacity. This study uncovered endogenous tissue remodeling mechanisms in the kidney that were activated by the loss of body fluid and salt and regulated by a unique niche of a minority renal cell type called the macula densa (MD). Here, we identified neuronal differentiation features of MD cells that sense the local and systemic environment and secrete angiogenic, growth, and extracellular matrix remodeling factors, cytokines and chemokines, and control resident progenitor cells. Serial intravital imaging, MD nerve growth factor receptor and Wnt mouse models, and transcriptome analysis revealed cellular and molecular mechanisms of these MD functions. Human and therapeutic translation studies illustrated the clinical potential of MD factors, including CCN1, as a urinary biomarker and therapeutic target in chronic kidney disease. The concept that a neuronally differentiated key sensory and regulatory cell type responding to organ-specific physiological inputs controls local progenitors to remodel or repair tissues may be applicable to other organs and diverse tissue-regenerative therapeutic strategies.

CDSKNNXMBD: a novel clustering framework for large-scale single-cell data based on a stable graph structure

BACKGROUND

Accurate and efficient cell grouping is essential for analyzing single-cell transcriptome sequencing (scRNA-seq) data. However, the existing clustering techniques often struggle to provide timely and accurate cell type groupings when dealing with datasets with large-scale or imbalanced cell types. Therefore, there is a need for improved methods that can handle the increasing size of scRNA-seq datasets while maintaining high accuracy and efficiency.

METHODS

We propose CDSKNNXMBD (Community Detection based on a Stable K-Nearest Neighbor Graph Structure), a novel single-cell clustering framework integrating partition clustering algorithm and community detection algorithm, which achieves accurate and fast cell type grouping by finding a stable graph structure.

RESULTS

We evaluated the effectiveness of our approach by analyzing 15 tissues from the human fetal atlas. Compared to existing methods, CDSKNN effectively counteracts the high imbalance in single-cell data, enabling effective clustering. Furthermore, we conducted comparisons across multiple single-cell datasets from different studies and sequencing techniques. CDSKNN is of high applicability and robustness, and capable of balancing the complexities of across diverse types of data. Most importantly, CDSKNN exhibits higher operational efficiency on datasets at the million-cell scale, requiring an average of only 6.33 min for clustering 1.46 million single cells, saving 33.3% to 99% of running time compared to those of existing methods.

CONCLUSIONS

The CDSKNN is a flexible, resilient, and promising clustering tool that is particularly suitable for clustering imbalanced data and demonstrates high efficiency on large-scale scRNA-seq datasets.

Crosstalk among podocytes, glomerular endothelial cells and mesangial cells in diabetic kidney disease: an updated review

Diabetic kidney disease (DKD) is a long-term and serious complication of diabetes that affects millions of people worldwide. It is characterized by proteinuria, glomerular damage, and renal fibrosis, leading to end-stage renal disease, and the pathogenesis is complex and involves multiple cellular and molecular mechanisms. Among three kinds of intraglomerular cells including podocytes, glomerular endothelial cells (GECs) and mesangial cells (MCs), the alterations in one cell type can produce changes in the others. The cell-to-cell crosstalk plays a crucial role in maintaining the glomerular filtration barrier (GFB) and homeostasis. In this review, we summarized the recent advances in understanding the pathological changes and interactions of these three types of cells in DKD and then focused on the signaling pathways and factors that mediate the crosstalk, such as angiopoietins, vascular endothelial growth factors, transforming growth factor-β, Krüppel-like factors, retinoic acid receptor response protein 1 and exosomes, etc. Furthermore, we also simply introduce the application of the latest technologies in studying cell interactions within glomerular cells and new promising mediators for cell crosstalk in DKD. In conclusion, this review provides a comprehensive and updated overview of the glomerular crosstalk in DKD and highlights its importance for the development of novel intervention approaches.

Glucocorticoid- and pioglitazone-induced proteinuria reduction in experimental NS both correlate with glomerular ECM modulation

Idiopathic nephrotic syndrome (NS) is a common glomerular disease. Although glucocorticoids (GC) are the primary treatment, the PPARγ agonist pioglitazone (Pio) also reduces proteinuria in patients with NS and directly protects podocytes from injury. Because both drugs reduce proteinuria, we hypothesized these effects result from overlapping transcriptional patterns. Systems biology approaches compared glomerular transcriptomes from rats with PAN-induced NS treated with GC vs. Pio and identified 29 commonly regulated genes-of-interest, primarily involved in extracellular matrix (ECM) remodeling. Correlation with clinical idiopathic NS patient datasets confirmed glomerular ECM dysregulation as a potential mechanism of injury. Cellular deconvolution in silico revealed GC- and Pio-induced amelioration of altered genes primarily within podocytes and mesangial cells. While validation studies are indicated, these analyses identified molecular pathways involved in the early stages of NS (prior to scarring), suggesting that targeting glomerular ECM dysregulation may enable a future non-immunosuppressive approach for proteinuria reduction in idiopathic NS.

Organotypic heterogeneity in microvascular endothelial cell responses in sepsis-a molecular treasure trove and pharmacological Gordian knot

In the last decades, it has become evident that endothelial cells (ECs) in the microvasculature play an important role in the pathophysiology of sepsis-associated multiple organ dysfunction syndrome (MODS). Studies on how ECs orchestrate leukocyte recruitment, control microvascular integrity and permeability, and regulate the haemostatic balance have provided a wealth of knowledge and potential molecular targets that could be considered for pharmacological intervention in sepsis. Yet, this information has not been translated into effective treatments. As MODS affects specific vascular beds, (organotypic) endothelial heterogeneity may be an important contributing factor to this lack of success. On the other hand, given the involvement of ECs in sepsis, this heterogeneity could also be leveraged for therapeutic gain to target specific sites of the vasculature given its full accessibility to drugs. In this review, we describe current knowledge that defines heterogeneity of organ-specific microvascular ECs at the molecular level and elaborate on studies that have reported EC responses across organ systems in sepsis patients and animal models of sepsis. We discuss hypothesis-driven, single-molecule studies that have formed the basis of our understanding of endothelial cell engagement in sepsis pathophysiology, and include recent studies employing high-throughput technologies. The latter deliver comprehensive data sets to describe molecular signatures for organotypic ECs that could lead to new hypotheses and form the foundation for rational pharmacological intervention and biomarker panel development. Particularly results from single cell RNA sequencing and spatial transcriptomics studies are eagerly awaited as they are expected to unveil the full spatiotemporal signature of EC responses to sepsis. With increasing awareness of the existence of distinct sepsis subphenotypes, and the need to develop new drug regimen and companion diagnostics, a better understanding of the molecular pathways exploited by ECs in sepsis pathophysiology will be a cornerstone to halt the detrimental processes that lead to MODS.

GPR56 Promotes Diabetic Kidney Disease Through eNOS Regulation in Glomerular Endothelial Cells

Although glomerular endothelial dysfunction is well recognized as contributing to the pathogenesis of diabetic kidney disease (DKD), the molecular pathways contributing to DKD pathogenesis in glomerular endothelial cells (GECs) are only partially understood. To uncover pathways that are differentially regulated in early DKD that may contribute to disease pathogenesis, we recently conducted a transcriptomic analysis of isolated GECs from diabetic NOS3-null mice. The analysis identified several potential mediators of early DKD pathogenesis, one of which encoded an adhesion G protein-coupled receptor-56 (GPR56), also known as ADGRG1. Enhanced glomerular expression of GPR56 was observed in human diabetic kidneys, which was negatively associated with kidney function. Using cultured mouse GECs, we observed that GPR56 expression was induced with exposure to advanced glycation end products, as well as in high-glucose conditions, and its overexpression resulted in decreased phosphorylation and expression of endothelial nitric oxide synthase (eNOS). This effect on eNOS by GPR56 was mediated by coupling of Gα12/13-RhoA pathway activation and Gαi-mediated cAMP/PKA pathway inhibition. The loss of GPR56 in mice led to a significant reduction in diabetes-induced albuminuria and glomerular injury, which was associated with reduced oxidative stress and restoration of eNOS expression in GECs. These findings suggest that GPR56 promotes DKD progression mediated, in part, through enhancing glomerular endothelial injury and dysfunction.

ARTICLE HIGHLIGHTS

The unique structural and functional characteristics of glomerular endothelial cell fenestrations and their potential as a therapeutic target in kidney disease

Glomerular endothelial cell (GEnC) fenestrations are a critical component of the glomerular filtration barrier. Their unique nondiaphragmed structure is key to their function in glomerular hydraulic permeability, and their aberration in disease can contribute to loss of glomerular filtration function. This review provides a comprehensive update of current understanding of the regulation and biogenesis of fenestrae. We consider diseases in which GEnC fenestration loss is recognized or may play a role and discuss methods with potential to facilitate the study of these critical structures. Literature is drawn from GEnCs as well as other fenestrated cell types such as liver sinusoidal endothelial cells that most closely parallel GEnCs.

Unique miRNome and transcriptome profiles underlie microvascular heterogeneity in mouse kidney

Endothelial cells in blood vessels in the kidney exert different functions depending on the (micro)vascular bed they are located in. The present study aimed to investigate microRNA and mRNA transcription patterns that underlie these differences. We zoomed in on microvascular compartments in the mouse renal cortex by laser microdissecting the microvessels prior to small RNA- and RNA-sequencing analyses. By these means, we characterized microRNA and mRNA transcription profiles of arterioles, glomeruli, peritubular capillaries, and postcapillary venules. Quantitative RT-PCR, in situ hybridization, and immunohistochemistry were used to validate sequencing results. Unique microRNA and mRNA transcription profiles were found in all microvascular compartments, with dedicated marker microRNAs and mRNAs showing enriched transcription in a single microvascular compartment. In situ hybridization validated the localization of microRNAs mmu-miR-140-3p in arterioles, mmu-miR-322-3p in glomeruli, and mmu-miR-451a in postcapillary venules. Immunohistochemical staining showed that von Willebrand factor protein was mainly expressed in arterioles and postcapillary venules, whereas GABRB1 expression was enriched in glomeruli, and IGF1 was enriched in postcapillary venules. More than 550 compartment-specific microRNA-mRNA interaction pairs were identified that carry functional implications for microvascular behavior. In conclusion, our study identified unique microRNA and mRNA transcription patterns in microvascular compartments of the mouse kidney cortex that underlie microvascular heterogeneity. These patterns provide important molecular information for future studies into differential microvascular engagement in health and disease.NEW & NOTEWORTHY Renal endothelial cells display a high level of heterogeneity depending on the (micro)vascular bed they reside in. The molecular basis contributing to these differences is poorly understood yet of high importance to increase understanding of microvascular engagement in the kidney in health and disease. This report describes m(icro)RNA expression profiles of microvascular beds in the mouse renal cortex and uncovers microvascular compartment-specific m(icro)RNAs and miRNA-mRNA pairs, thereby revealing important molecular mechanisms underlying renal microvascular heterogeneity.

Advances and challenges toward developing kidney organoids for clinical applications

Kidney organoids have enabled modeling of human development and disease. While methods of generating the nephron lineage are well established, new protocols to induce another lineage, the ureteric bud/collecting duct, have been reported in the past 5 years. Many reports have described modeling of various hereditary kidney diseases, with polycystic kidney disease serving as the archetypal disease, by using patient-derived or genome-edited kidney organoids. The generation of more organotypic kidneys is also becoming feasible. In this review, I also discuss the significant challenges for more sophisticated disease modeling and for realizing the ambitious goal of generating transplantable synthetic kidneys.

Identification of renal cyst cells of type I Nephronophthisis by single-nucleus RNA sequencing

Background: Nephronophthisis (NPH) is the most common genetic cause of end-stage renal disease (ESRD) in childhood, and NPHP1 is the major pathogenic gene. Cyst formation at the corticomedullary junction is a pathological feature of NPH, but the mechanism underlying cystogenesis is not well understood. The isolation and identification of cystic cell subpopulation could help to identify their origins and provide vital clues to the mechanisms underlying cystogenesis in NPH. Methods: Single-nucleus RNA sequencing (snRNA-seq) was performed to produce an atlas of NPHP1 renal cells. Kidney samples were collected from WT (Nphp1 +/+) mice and NPHP1 (Nphp1 del2-20/del2-20) model mice. Results: A comprehensive atlas of the renal cellular landscape in NPHP1 was generated, consisting of 14 basic renal cell types as well as a subpopulation of DCT cells that was overrepresented in NPHP1 kidneys compared to WT kidneys. GO analysis revealed significant downregulation of genes associated with tubular development and kidney morphogenesis in this subpopulation. Furthermore, the reconstruction of differentiation trajectories of individual cells within this subpopulation confirmed that a specific group of cells in NPHP1 mice become arrested at an early stage of differentiation and proliferate to form cysts. We demonstrate that Niban1 is a specific molecular marker of cystic cells in both mice and human NPHP1. Conclusion: In summary, we report a novel subpopulation of DCT cells, marked by Niban1, that are classified as cystic cells in the NPHP1 mice kidney. These results offer fresh insights into the cellular and molecular basis of cystogenesis in NPH.

Single-cell transcriptomic profiles in the pathophysiology within the microenvironment of early diabetic kidney disease

Diabetic kidney disease (DKD) is the leading cause of end-stage kidney disease, resulting in a huge socio-economic impact. Kidney is a highly complex organ and the pathogenesis underlying kidney organization involves complex cell-to-cell interaction within the heterogeneous kidney milieu. Advanced single-cell RNA sequencing (scRNA-seq) could reveal the complex architecture and interaction with the microenvironment in early DKD. We used scRNA-seq to investigate early changes in the kidney of db/m mice and db/db mice at the 14th week. Uniform Manifold Approximation and Projection were applied to classify cells into different clusters at a proper resolution. Weighted gene co-expression network analysis was used to identify the key molecules specifically expressed in kidney tubules. Information of cell-cell communication within the kidney was obtained using receptor-ligand pairing resources. In vitro model, human subjects, and co-detection by indexing staining were used to identify the pathophysiologic role of the hub genes in DKD. Among four distinct subsets of the proximal tubule (PT), lower percentages of proliferative PT and PT containing AQP4 expression (PT^AQP4+) in db/db mice induced impaired cell repair activity and dysfunction of renin-angiotensin system modulation in early DKD. We found that ferroptosis was involved in DKD progression, and ceruloplasmin acted as a central regulator of the induction of ferroptosis in PT^AQP4+. In addition, lower percentages of thick ascending limbs and collecting ducts with impaired metabolism function were also critical pathogenic features in the kidney of db/db mice. Secreted phosphoprotein 1 (SPP1) mediated pathogenic cross-talk in the tubular microenvironment, as validated by a correlation between urinary SPP1/Cr level and tubular injury. Finally, mesangial cell-derived semaphorin 3C (SEMA3C) further promoted endothelium-mesenchymal transition in glomerular endothelial cells through NRP1 and NRP2, and urinary SEMA3C/Cr level was positively correlated with glomerular injury. These data identified the hub genes involved in pathophysiologic changes within the microenvironment of early DKD.

Single cell landscape of parietal epithelial cells in healthy and diseased states

The biology and diversity of glomerular parietal epithelial cells (PECs) are important for understanding podocyte regeneration and crescent formation. Although protein markers have revealed the morphological heterogeneity of PECs, the molecular characteristics of PEC subpopulations remain largely unknown. Here, we performed a comprehensive analysis of PECs using single-cell RNA sequencing (scRNA-seq) data. Our analysis identified five distinct PEC subpopulations: PEC-A1, PEC-A2, PEC-A3, PEC-A4 and PEC-B. Among these subpopulations, PEC- A1 and PEC-A2 were characterized as podocyte progenitors while PEC-A4 represented tubular progenitors. Further dynamic signaling network analysis indicated that activation of PEC-A4 and the proliferation of PEC-A3 played pivotal roles in crescent formation. Analyses suggested that upstream signals released by podocytes, immune cells, endothelial cells and mesangial cells serve as pathogenic signals and may be promising intervention targets in crescentic glomerulonephritis. Pharmacological blockade of two such pathogenic signaling targets, proteins Mif and Csf1r, reduced hyperplasia of the PECs and crescent formation in anti-glomerular basement membrane glomerulonephritis murine models. Thus, our study demonstrates that scRNA-seq-based analysis provided valuable insights into the pathology and therapeutic strategies for crescentic glomerulonephritis.

Single-cell RNA sequencing identifies response of renal lymphatic endothelial cells to acute kidney injury

The inflammatory response to acute kidney injury (AKI) likely dictates future renal health. Lymphatic vessels are responsible for maintaining tissue homeostasis through transport and immunomodulatory roles. Due to the relative sparsity of lymphatic endothelial cells (LECs) in the kidney, past sequencing efforts have not characterized these cells and their response to AKI. Here we characterized murine renal LEC subpopulations by single-cell RNA sequencing and investigated their changes in cisplatin AKI. We validated our findings by qPCR in LECs isolated from both cisplatin-injured and ischemia reperfusion injury, by immunofluorescence, and confirmation in in vitro human LECs. We have identified renal LECs and their lymphatic vascular roles that have yet to be characterized in previous studies. We report unique gene changes mapped across control and cisplatin injured conditions. Following AKI, renal LECs alter genes involved endothelial cell apoptosis and vasculogenic processes as well as immunoregulatory signaling and metabolism. Differences between injury models are also identified with renal LECs further demonstrating changed gene expression between cisplatin and ischemia reperfusion injury models, indicating the renal LEC response is both specific to where they lie in the lymphatic vasculature and the renal injury type. How LECs respond to AKI may therefore be key in regulating future kidney disease progression.

Integrated single-nucleus sequencing and spatial architecture analysis identified distinct injured-proximal tubular types in calculi rats

BACKGROUND

Urolithiasis with high prevalence and recurrence rate, has impacts on kidney injury in patients, becomes a socioeconomic and healthcare problem in worldwide. However, the biology of kidney with crystal formation and proximal tubular injury remains essentially unclear. The present study aims to evaluate the cell biology and immune-communications in urolithiasis mediated kidney injury, to provide new insights in the kidney stone treatment and prevention.

RESULTS

We identified 3 distinct injured-proximal tubular cell types based on the differentially expression injury markers (Havcr1 and lcn2) and functional solute carriers (slc34a3, slc22a8, slc38a3 and slc7a13), and characterized 4 main immune cell types in kidney and one undefined cell population, where F13a1^+/high/CD163^+/high monocyte & macrophage and Sirpa/Fcgr1a/Fcgr2a^+/high granulocyte were the most enriched. We performed intercellular crosstalk analysis based on the snRNA-seq data and explored the potential immunomodulation of calculi stone formation, and founded that the interaction between ligand Gas6 and its receptors (Gas6-Axl, Gas6-Mertk) was specifically observed in the injured-PT1 cells, but not injured-PT2 and -PT3 cells. The interaction of Ptn-Plxnb2 was only observed between the injured-PT3 cells and its receptor enriched cells.

CONCLUSIONS

Present study comprehensively characterized the gene expression profile in the calculi rat kidney at single nucleus level, identified novel marker genes for all cell types of rat kidney, and determined 3 distinct sub-population of injured-PT clusters, as well as intercellular communication between injured-PTs and immune cells. Our collection of data provides a reliable resource and reference for studies on renal cell biology and kidney disease.

Assessing the effects of aging on the renal endothelial cell landscape using single-cell RNA sequencing

Endothelial cells (ECs) with senescence-associated secretory phenotypes (SASP) have been identified as a key mechanism of aging that contributes to various age-related kidney diseases. In this study, we used single-cell RNA sequencing (scRNA-seq) to create a transcriptome atlas of murine renal ECs and identify transcriptomic changes that occur during aging. We identified seven different subtypes of renal ECs, with glomerular ECs and angiogenic ECs being the most affected by senescence. We confirmed our scRNA-seq findings by using double immunostaining for an EC marker (CD31) and markers of specialized EC phenotypes. Our analysis of the dynamics of capillary lineage development revealed a chronic state of inflammation and compromised glomerular function as prominent aging features. Additionally, we observed an elevated pro-inflammatory and pro-coagulant microenvironment in aged glomerular ECs, which may contribute to age-related glomerulosclerosis and renal fibrosis. Through intercellular communication analysis, we also identified changes in signaling involved in immune regulation that may contribute to a hostile microenvironment for renal homeostasis and function. Overall, our findings provide new insights into the mechanisms of aging in the renal endothelium and may pave the way for the discovery of diagnostic biomarkers and therapeutic interventions against age-related kidney diseases.

The tissue-specific transcriptional landscape underlines the involvement of endothelial cells in health and disease

Endothelial cells (ECs) that line vascular and lymphatic vessels are being increasingly recognized as important to organ function in health and disease. ECs participate not only in the trafficking of gases, metabolites, and cells between the bloodstream and tissues but also in the angiocrine-based induction of heterogeneous parenchymal cells, which are unique to their specific tissue functions. The molecular mechanisms regulating EC heterogeneity between and within different tissues are modeled during embryogenesis and become fully established in adults. Any changes in adult tissue homeostasis induced by aging, stress conditions, and various noxae may reshape EC heterogeneity and induce specific transcriptional features that condition a functional phenotype. Heterogeneity is sustained via specific genetic programs organized through the combinatory effects of a discrete number of transcription factors (TFs) that, at the single tissue-level, constitute dynamic networks that are post-transcriptionally and epigenetically regulated. This review is focused on outlining the TF-based networks involved in EC specialization and physiological and pathological stressors thought to modify their architecture.

Single-Cell Transcriptome Sequencing Reveals Molecular Mechanisms of Renal Injury in Essential Hypertension

INTRODUCTION

Hypertensive nephropathy is characterized by glomerular and tubulointerstitial damage, but we know little about changes in cell-specific gene expression in the early stages of hypertensive kidney injury, which usually has no obvious pathological changes.

METHODS

We performed unbiased single-cell RNA sequencing of rat kidney samples from hypertensive kidney injury to generate 10,602 single-cell transcriptomes from 2 control and 2 early stage hypertensive kidney injury samples.

RESULTS

All major cell types of the kidney were represented in the final dataset. Side-by-side comparisons showed that cell type-specific changes in gene expression are critical for functional impairment of glomeruli and tubules and activation of immune cells. In particular, we found a significantly reduced gene expression profile of maintaining vascular integrity in glomerular cells of hypertensive kidney injury. Meanwhile, the expression of genes associated with oxidative stress injury and fibrosis in the renal tubules and collecting ducts was elevated, but the degree of tubular cells response to injury differed between parts. We also found a signature of immune cell infiltration in hypertensive kidney injury.

CONCLUSION

Exploring the changes of gene expression in hypertension-injured kidneys may be helpful to identify the early biomarkers and signal pathways of this disease. Our data provide rich resources for understanding the pathogenesis of hypertensive renal injury and formulating effective treatment strategies.

Identification of IgA autoantibodies targeting mesangial cells redefines the pathogenesis of IgA nephropathy

Immunoglobulin A (IgA) nephropathy (IgAN) is the most common type of primary glomerulonephritis, often progressing to renal failure. IgAN is triggered by IgA deposition in the glomerular mesangium by an undefined mechanism. Here, we show that grouped ddY (gddY) mice, a spontaneous IgAN model, produce serum IgA against mesangial antigens, including βII-spectrin. Most patients with IgAN also have serum anti-βII-spectrin IgA. As in patients with IgAN, IgA⁺ plasmablasts accumulate in the kidneys of gddY mice. IgA antibodies cloned from the plasmablasts carry substantial V-region mutations and bind to βII-spectrin and the surface of mesangial cells. These IgAs recognize transfected and endogenous βII-spectrin exposed on the surface of embryonic kidney-derived cells. Last, we demonstrate that the cloned IgA can bind selectively to glomerular mesangial regions in situ. The identification of IgA autoantibody and its antigen in IgAN provides key insights into disease onset and redefines IgAN as a tissue-specific autoimmune disease.

Spatial transcriptomics using multiplexed deterministic barcoding in tissue

Spatially resolved transcriptomics of tissue sections enables advances in fundamental and applied biomedical research. Here, we present Multiplexed Deterministic Barcoding in Tissue (xDBiT) to acquire spatially resolved transcriptomes of nine tissue sections in parallel. New microfluidic chips were developed to spatially encode mRNAs over a total tissue area of 1.17 cm² with a 50 µm resolution. Optimization of the biochemical protocol increased read and gene counts per spot by one order of magnitude compared to previous reports. Furthermore, the introduction of alignment markers allowed seamless registration of images and spatial transcriptomic spots. Together with technological advances, we provide an open-source computational pipeline to prepare raw sequencing data for downstream analysis. The functionality of xDBiT was demonstrated by acquiring 16 spatially resolved transcriptomic datasets from five different murine organs, including the cerebellum, liver, kidney, spleen, and heart. Factor analysis and deconvolution of spatial transcriptomes allowed for in-depth characterization of the murine kidney.

Deciphering endothelial heterogeneity in health and disease at single-cell resolution: progress and perspectives

Endothelial cells (ECs) constitute the inner lining of vascular beds in mammals and are crucial for homeostatic regulation of blood vessel physiology, but also play a key role in pathogenesis of many diseases, thereby representing realistic therapeutic targets. However, it has become evident that ECs are heterogeneous, encompassing several subtypes with distinct functions, which makes EC targeting and modulation in diseases challenging. The rise of the new single-cell era has led to an emergence of studies aimed at interrogating transcriptome diversity along the vascular tree, and has revolutionized our understanding of EC heterogeneity from both a physiological and pathophysiological context. Here, we discuss recent landmark studies aimed at teasing apart the heterogeneous nature of ECs. We cover driving (epi)genetic, transcriptomic, and metabolic forces underlying EC heterogeneity in health and disease, as well as current strategies used to combat disease-enriched EC phenotypes, and propose strategies to transcend largely descriptive heterogeneity towards prioritization and functional validation of therapeutically targetable drivers of EC diversity. Lastly, we provide an overview of the most recent advances and hurdles in single EC OMICs.

Esm-1 mediates transcriptional polarization associated with diabetic kidney disease

Background

Esm-1, endothelial cell-specific molecule-1, is a susceptibility gene for diabetic kidney disease (DKD) and is a cytokine- and glucose-regulated, secreted proteoglycan, that is notably expressed in kidney and attenuates inflammation and albuminuria. Esm1 has restricted expression at the vascular tip during development but little is known about its expression pattern in mature tissues, and its precise effects in diabetes.

Methods

We utilized publicly available single-cell RNA sequencing data to explore the characteristics of Esm1 expression in 27,786 renal endothelial cells obtained from four adult human and three mouse databases. We validated our findings using bulk transcriptome data from an additional 20 healthy subjects and 41 patients with DKD and using RNAscope. Using correlation matrices, we relate Esm1 expression to the glomerular transcriptome and evaluated these matrices with systemic over-expression of Esm-1.

Results

In both mice and humans, Esm1 is expressed in a subset of all renal endothelial cell types and represents a minority of glomerular endothelial cells. In patients, Esm1 (+) cells exhibit a highly conserved enrichment for blood vessel development genes. With diabetes, these cells are fewer in number and profoundly shift expression to reflect chemotaxis pathways. Analysis of these gene sets highlight candidate genes such as Igfbp5 for cross talk between cell types. We also find that diabetes induces correlations in the expression of large clusters of genes, within cell type-enriched transcripts. Esm1 significantly correlates with a majority genes within these clusters, delineating a glomerular transcriptional polarization reflected by the magnitude of Esm1 deficiency. In diabetic mice, these gene clusters link Esm1 expression to albuminuria, and over-expression of Esm-1 reverses the expression pattern in many of these genes.

Conclusions

A comprehensive analysis of single cell and bulk transcriptomes demonstrates that diabetes correlates with lower Esm1 expression and with changes in the functional characterization of Esm1 (+) cells. Esm1 is both a marker for glomerular transcriptional polarization, and a mediator that re-orients the transcriptional program in DKD.

Vascular endothelial cell development and diversity

Vascular endothelial cells form the inner layer of blood vessels where they have a key role in the development and maintenance of the functional circulatory system and provide paracrine support to surrounding non-vascular cells. Technical advances in the past 5 years in single-cell genomics and in in vivo genetic labelling have facilitated greater insights into endothelial cell development, plasticity and heterogeneity. These advances have also contributed to a new understanding of the timing of endothelial cell subtype differentiation and its relationship to the cell cycle. Identification of novel tissue-specific gene expression patterns in endothelial cells has led to the discovery of crucial signalling pathways and new interactions with other cell types that have key roles in both tissue maintenance and disease pathology. In this Review, we describe the latest findings in vascular endothelial cell development and diversity, which are often supported by large-scale, single-cell studies, and discuss the implications of these findings for vascular medicine. In addition, we highlight how techniques such as single-cell multimodal omics, which have become increasingly sophisticated over the past 2 years, are being utilized to study normal vascular physiology as well as functional perturbations in disease.

Soluble Klotho protects against glomerular injury through regulation of ER stress response

αKlotho (Klotho) has well established renoprotective effects; however, the molecular pathways mediating its glomerular protection remain incompletely understood. Recent studies have reported that Klotho is expressed in podocytes and protects glomeruli through auto- and paracrine effects. Here, we examined renal expression of Klotho in detail and explored its protective effects in podocyte-specific Klotho knockout mice, and by overexpressing human Klotho in podocytes and hepatocytes. We demonstrate that Klotho is not significantly expressed in podocytes, and transgenic mice with either a targeted deletion or overexpression of Klotho in podocytes lack a glomerular phenotype and have no altered susceptibility to glomerular injury. In contrast, mice with hepatocyte-specific overexpression of Klotho have high circulating levels of soluble Klotho, and when challenged with nephrotoxic serum have less albuminuria and less severe kidney injury compared to wildtype mice. RNA-seq analysis suggests an adaptive response to increased endoplasmic reticulum stress as a putative mechanism of action. To evaluate the clinical relevance of our findings, the results were validated in patients with diabetic nephropathy, and in precision cut kidney slices from human nephrectomies. Together, our data reveal that the glomeruloprotective effects of Klotho is mediated via endocrine actions, which increases its therapeutic potential for patients with glomerular diseases.

RefTM: reference-guided topic modeling of single-cell chromatin accessibility data

Single-cell analysis is a valuable approach for dissecting the cellular heterogeneity, and single-cell chromatin accessibility sequencing (scCAS) can profile the epigenetic landscapes for thousands of individual cells. It is challenging to analyze scCAS data, because of its high dimensionality and a higher degree of sparsity compared with scRNA-seq data. Topic modeling in single-cell data analysis can lead to robust identification of the cell types and it can provide insight into the regulatory mechanisms. Reference-guided approach may facilitate the analysis of scCAS data by utilizing the information in existing datasets. We present RefTM (Reference-guided Topic Modeling of single-cell chromatin accessibility data), which not only utilizes the information in existing bulk chromatin accessibility and annotated scCAS data, but also takes advantage of topic models for single-cell data analysis. RefTM simultaneously models: (1) the shared biological variation among reference data and the target scCAS data; (2) the unique biological variation in scCAS data; (3) other variations from known covariates in scCAS data.

Single-cell transcriptomics: A new tool for studying diabetic kidney disease

The kidney is a complex organ comprising various functional partitions and special cell types that play important roles in maintaining homeostasis in the body. Diabetic kidney disease (DKD) is the leading cause of end-stage renal disease and is an independent risk factor for cardiovascular diseases. Owing to the complexity and heterogeneity of kidney structure and function, the mechanism of DKD development has not been fully elucidated. Single-cell sequencing, including transcriptomics, epigenetics, metabolomics, and proteomics etc., is a powerful technology that enables the analysis of specific cell types and states, specifically expressed genes or pathways, cell differentiation trajectories, intercellular communication, and regulation or co-expression of genes in various diseases. Compared with other omics, RNA sequencing is a more developed technique with higher utilization of tissues or samples. This article reviewed the application of single-cell transcriptomics in the field of DKD and highlighted the key signaling pathways in specific tissues or cell types involved in the occurrence and development of DKD. The comprehensive understanding of single-cell transcriptomics through single-cell RNA-seq and single-nucleus RNA-seq will provide us new insights into the pathogenesis and treatment strategy of various diseases including DKD.

Protocol for Classification Single-Cell PBMC Types from Pathological Samples Using Supervised Machine Learning

Peripheral blood mononuclear cells (PBMC) are mixed subpopulations of blood cells composed of five cell types. PBMC are widely used in the study of the immune system, infectious diseases, cancer, and vaccine development. Single-cell transcriptomics (SCT) allows the labeling of cell types by gene expression patterns from biological samples. Classifying cells into cell types and states is essential for single-cell analyses, especially in the classification of diseases and the assessment of therapeutic interventions, and for many secondary analyses. Most of the classification of cell types from SCT data use unsupervised clustering or a combination of unsupervised and supervised methods including manual correction. In this chapter, we describe a protocol that uses supervised machine learning (ML) methods with SCT data for the classification of PBMC cell types in samples representing pathological states. This protocol has three parts: (1) data preprocessing, (2) labeling of reference PBMC SCT datasets and training supervised ML models, and (3) labeling new PBMC datasets from disease samples. This protocol enables building classification models that are of high accuracy and efficiency. Our example focuses on 10× Genomics technology but applies to datasets from other SCT platforms.

Paraoxonase 2 (PON2) Deficiency Reproduces Lipid Alterations of Diabetic and Inflammatory Glomerular Disease and Affects TRPC6 Signaling

Diabetes and inflammatory diseases are associated with an altered cellular lipid composition due to lipid peroxidation. The pathogenic potential of these lipid alterations in glomerular kidney diseases remains largely obscure as suitable cell culture and animal models are lacking. In glomerular disease, a loss of terminally differentiated glomerular epithelial cells called podocytes refers to irreversible damage. Podocytes are characterized by a complex ramified cellular architecture and highly active transmembrane signaling. Alterations in lipid composition in states of disease have been described in podocytes but the pathophysiologic mechanisms mediating podocyte damage are unclear. In this study, we employ a genetic deletion of the anti-oxidative, lipid-modifying paraoxonase 2 enzyme (PON2) as a model to study altered cellular lipid composition and its effects on cellular signaling in glomerular disease. PON2 deficiency reproduces features of an altered lipid composition of glomerular disease, characterized by an increase in ceramides and cholesterol. PON2 knockout mice are more susceptible to glomerular damage in models of aggravated oxidative stress such as adriamycin-induced nephropathy. Voltage clamp experiments in cultured podocytes reveal a largely increased TRPC6 conductance after a membrane stretch in PON2 deficiency. Correspondingly, a concomitant knockout of TRPC6 and PON2 partially rescues the aggravated glomerular phenotype of a PON2 knockout in the adriamycin model. This study establishes PON2 deficiency as a model to investigate the pathophysiologic mechanisms of podocyte dysfunction related to alterations in the lipid composition, as seen in diabetic and inflammatory glomerular disease. Expanding the knowledge on these routes and options of intervention could lead to novel treatment strategies for glomerular disease.

Blood pressure and the kidney cortex transcriptome response to high-sodium diet challenge in female nonhuman primates

Blood pressure (BP) is influenced by genetic variation and sodium intake with sex-specific differences; however, studies to identify renal molecular mechanisms underlying the influence of sodium intake on BP in nonhuman primates (NHP) have focused on males. To address the gap in our understanding of molecular mechanisms regulating BP in female primates, we studied sodium-naïve female baboons (n = 7) fed a high-sodium (HS) diet for 6 wk. We hypothesized that in female baboons variation in renal transcriptional networks correlates with variation in BP response to a high-sodium diet. BP was continuously measured for 64-h periods throughout the study by implantable telemetry devices. Sodium intake, blood samples for clinical chemistries, and ultrasound-guided kidney biopsies were collected before and after the HS diet for RNA-Seq and bioinformatic analyses. We found that on the LS diet but not the HS diet, sodium intake and serum 17 β-estradiol concentration correlated with BP. Furthermore, kidney transcriptomes differed by diet-unbiased weighted gene coexpression network analysis revealed modules of genes correlated with BP on the HS diet but not the LS diet. Our results showed variation in BP on the HS diet correlated with variation in novel kidney gene networks regulated by ESR1 and MYC; i.e., these regulators have not been associated with BP regulation in male humans or rodents. Validation of the mechanisms underlying regulation of BP-associated gene networks in female NHP will inform better therapies toward greater precision medicine for women.

Early Molecular Events Mediating Loss of Aquaporin-2 during Ureteral Obstruction in Rats

BACKGROUND

Ureteral obstruction is marked by disappearance of the vasopressin-dependent water channel aquaporin-2 (AQP2) in the renal collecting duct and polyuria upon reversal. Most studies of unilateral ureteral obstruction (UUO) models have examined late time points, obscuring the early signals that trigger loss of AQP2.

METHODS

We performed RNA-Seq on microdissected rat cortical collecting ducts (CCDs) to identify early signaling pathways after establishment of UUO.

RESULTS

Vasopressin V2 receptor (AVPR2) mRNA was decreased 3 hours after UUO, identifying one cause of AQP2 loss. Collecting duct principal cell differentiation markers were lost, including many not regulated by vasopressin. Immediate early genes in CCDs were widely induced 3 hours after UUO, including Myc, Atf3, and Fos (confirmed at the protein level). Simultaneously, expression of NF-κB signaling response genes known to repress Aqp2 increased. RNA-Seq for CCDs at an even earlier time point (30 minutes) showed widespread mRNA loss, indicating a "stunned" profile. Immunocytochemical labeling of markers of mRNA-degrading P-bodies DDX6 and 4E-T indicated an increase in P-body formation within 30 minutes.

CONCLUSIONS

Immediately after establishment of UUO, collecting ducts manifest a stunned state with broad disappearance of mRNAs. Within 3 hours, there is upregulation of immediate early and inflammatory genes and disappearance of the V2 vasopressin receptor, resulting in loss of AQP2 (confirmed by lipopolysaccharide administration). The inflammatory response seen rapidly after UUO establishment may be relevant to both UUO-induced polyuria and long-term development of fibrosis in UUO kidneys.

Integrative transcriptome analysis reveals TEKT2 and PIAS2 involvement in diabetic nephropathy

Cell heterogeneity has impeded the accurate interpretation of the bulk transcriptome data from patients with diabetic nephropathy (DN). We performed an analysis by integrating bulk and single-cell transcriptome datasets to uncover novel mechanisms leading to DN, especially in the podocytes. Microdissected glomeruli and tubules transcriptome datasets were selected from Gene Expression Omnibus (GEO). Then the consistency between datasets was evaluated. The analysis of the bulk dataset and single-nucleus RNA dataset was integrated to reveal the cell type-specific responses to DN. The candidate genes were validated in kidney tissues from DN patients and diabetic mice. We compared 4 glomerular and 4 tubular datasets and found considerable discrepancies among datasets regarding the deferentially expressed genes (DEGs), involved signaling pathways, and the hallmark enrichment profiles. Deconvolution of the bulk data revealed that the variations in cell-type proportion contributed greatly to this discrepancy. The integrative analysis uncovered that the dysregulation of spermatogenesis-related genes, including TEKT2 and PIAS2, was involved in the development of DN. Importantly, the mRNA level of TEKT2 was negatively correlated with the mRNA levels of NPHS1 (r = -.66, p < .0001) and NPHS2 (r = -.85, p < .0001) in human diabetic glomeruli. Immunostaining confirmed that the expression of TEKT2 and PIAS2 were up-regulated in podocytes of DN patients and diabetic mice. Knocking down TEKT2 resisted high glucose-induced cytoskeletal remodeling and down-regulation of NPHS1 protein in the cultured podocyte. In conclusion, the integrative strategy can help us efficiently use the publicly available transcriptomics resources. Using this approach and combining it with classical research methods, we identified TEKT2 and PIAS2, two spermatogenesis-related genes involved in the pathogenesis of DN. Furthermore, TEKT2 is involved in this pathogenesis by regulating the podocyte cytoskeleton.

Bioinformatics analysis of potential pathogenesis and risk genes of immunoinflammation-promoted renal injury in severe COVID-19

Renal injury secondary to COVID-19 is an important factor for the poor prognosis of COVID-19 patients. The pathogenesis of renal injury caused by aberrant immune inflammatory of COVID-19 remains unclear. In this study, a total of 166 samples from 4 peripheral blood transcriptomic datasets of COVID-19 patients were integrated. By using the weighted gene co-expression network (WGCNA) algorithm, we identified key genes for mild, moderate, and severe COVID-19. Subsequently, taking these genes as input genes, we performed Short Time-series Expression Miner (STEM) analysis in a time consecutive ischemia-reperfusion injury (IRI) -kidney dataset to identify genes associated with renal injury in COVID-19. The results showed that only in severe COVID-19 there exist a small group of genes associated with the progression of renal injury. Gene enrichment analysis revealed that these genes are involved in extensive immune inflammation and cell death-related pathways. A further protein-protein interaction (PPI) network analysis screened 15 PPI-hub genes: ALOX5, CD38, GSF3R, LGR, RPR1, HCK, ITGAX, LYN, MAPK3, NCF4, SELP, SPI1, WAS, TLR2 and TLR4. Single-cell sequencing analysis indicated that PPI-hub genes were mainly distributed in neutrophils, macrophages, and dendritic cells. Intercellular ligand-receptor analysis characterized the activated ligand-receptors between these immune cells and parenchyma cells in depth. And KEGG enrichment analysis revealed that viral protein interaction with cytokine and cytokine receptor, necroptosis, and Toll-like receptor signaling pathway may be potentially essential for immune cell infiltration leading to COVID-19 renal injury. Finally, we validated the expression pattern of PPI-hub genes in an independent data set by random forest. In addition, we found that the high expression of these genes was correlated with a low glomerular filtration rate. Including them as risk genes in lasso regression, we constructed a Nomogram model for predicting severe COVID-19. In conclusion, our study explores the pathogenesis of renal injury promoted by immunoinflammatory in severe COVID-19 and extends the clinical utility of its key genes.

Systemic gene therapy with thymosin β4 alleviates glomerular injury in mice

Plasma ultrafiltration in the kidney occurs across glomerular capillaries, which are surrounded by epithelial cells called podocytes. Podocytes have a unique shape maintained by a complex cytoskeleton, which becomes disrupted in glomerular disease resulting in defective filtration and albuminuria. Lack of endogenous thymosin β4 (TB4), an actin sequestering peptide, exacerbates glomerular injury and disrupts the organisation of the podocyte actin cytoskeleton, however, the potential of exogenous TB4 therapy to improve podocyte injury is unknown. Here, we have used Adriamycin (ADR), a toxin which injures podocytes and damages the glomerular filtration barrier leading to albuminuria in mice. Through interrogating single-cell RNA-sequencing data of isolated glomeruli we demonstrate that ADR injury results in reduced levels of podocyte TB4. Administration of an adeno-associated viral vector encoding TB4 increased the circulating level of TB4 and prevented ADR-induced podocyte loss and albuminuria. ADR injury was associated with disorganisation of the podocyte actin cytoskeleton in vitro, which was ameliorated by treatment with exogenous TB4. Collectively, we propose that systemic gene therapy with TB4 prevents podocyte injury and maintains glomerular filtration via protection of the podocyte cytoskeleton thus presenting a novel treatment strategy for glomerular disease.

Consensus draft of the native mouse podocyte-ome

The podocyte is a key cell in maintaining renal filtration barrier integrity. Several recent studies have analyzed the entity of genome-coded molecules in the podocyte at deep resolution. This avenue of "podocyte-ome" research was enabled by a variety of techniques, including single-cell transcriptomics, FACS-sorting with and without genetically encoded markers, and deep acquisition of proteomics. However, data across various omics studies are not well-integrated with each other. Here, we aim to establish a common, simplified knowledgebase for the mouse "podocyte-ome" by integrating bulk RNA sequencing and bulk proteomics of sorted podocytes and single cell transcriptomics. Three datasets of each omics type from different laboratories, respectively, were integrated, visualized and bioinformatically analyzed. The procedure sheds light on conserved processes of podocytes, but also on limitations and specific features of the used technologies. High expression of glycan GPI anchor synthesis and turnover, and retinol metabolism was identified as a relatively understudied feature of podocytes, while there are both podocyte-enriched and podocyte-depleted actin binding molecules. We compiled aggregated data in an application that illustrates the features of the dataset and allows for exploratory analyses through individual gene query of podocyte identity in absolute and relative quantification towards other glomerular cell types, keywords, GO-terms and gene set enrichments. This consensus draft is a first step towards common molecular omics knowledge of kidney cells.

Identification, discrimination and heterogeneity of fibroblasts

Fibroblasts, the principal cell type of connective tissue, secrete extracellular matrix components during tissue development, homeostasis, repair and disease. Despite this crucial role, the identification and distinction of fibroblasts from other cell types are challenging and laden with caveats. Rapid progress in single-cell transcriptomics now yields detailed molecular portraits of fibroblasts and other cell types in our bodies, which complement and enrich classical histological and immunological descriptions, improve cell class definitions and guide further studies on the functional heterogeneity of cell subtypes and states, origins and fates in physiological and pathological processes. In this review, we summarize and discuss recent advances in the understanding of fibroblast identification and heterogeneity and how they discriminate from other cell types.

In this review, the authors look at how recent progress in single-cell transcriptomics complement and enrich the classical, largely morphological, portraits of fibroblasts. The detailed molecular information now available provides new insights into fibroblast identity, heterogeneity and function.

Reduced Glomerular Filtration in Diabetes Is Attributable to Loss of Density and Increased Resistance of Glomerular Endothelial Cell Fenestrations

BACKGROUND

Glomerular endothelial cell (GEnC) fenestrations are recognized as an essential component of the glomerular filtration barrier, yet little is known about how they are regulated and their role in disease.

METHODS

We comprehensively characterized GEnC fenestral and functional renal filtration changes including measurement of glomerular Kf and GFR in diabetic mice (BTBR ob-/ob- ). We also examined and compared human samples. We evaluated Eps homology domain protein-3 (EHD3) and its association with GEnC fenestrations in diabetes in disease samples and further explored its role as a potential regulator of fenestrations in an in vitro model of fenestration formation using b.End5 cells.

RESULTS

Loss of GEnC fenestration density was associated with decreased filtration function in diabetic nephropathy. We identified increased diaphragmed fenestrations in diabetes, which are posited to increase resistance to filtration and further contribute to decreased GFR. We identified decreased glomerular EHD3 expression in diabetes, which was significantly correlated with decreased fenestration density. Reduced fenestrations in EHD3 knockdown b.End5 cells in vitro further suggested a mechanistic role for EHD3 in fenestration formation.

CONCLUSIONS

This study demonstrates the critical role of GEnC fenestrations in renal filtration function and suggests EHD3 may be a key regulator, loss of which may contribute to declining glomerular filtration function through aberrant GEnC fenestration regulation. This points to EHD3 as a novel therapeutic target to restore filtration function in disease.