The integration of single cell transcriptome and chromatin accessibility datasets enables a deeper understanding of cell heterogeneity. We performed single nucleus ATAC (snATAC-seq) and RNA ...(snRNA-seq) sequencing to generate paired, cell-type-specific chromatin accessibility and transcriptional profiles of the adult human kidney. We demonstrate that snATAC-seq is comparable to snRNA-seq in the assignment of cell identity and can further refine our understanding of functional heterogeneity in the nephron. The majority of differentially accessible chromatin regions are localized to promoters and a significant proportion are closely associated with differentially expressed genes. Cell-type-specific enrichment of transcription factor binding motifs implicates the activation of NF-κB that promotes VCAM1 expression and drives transition between a subpopulation of proximal tubule epithelial cells. Our multi-omics approach improves the ability to detect unique cell states within the kidney and redefines cellular heterogeneity in the proximal tubule and thick ascending limb.
A challenge for single-cell genomic studies in kidney and other solid tissues is generating a high-quality single-cell suspension that contains rare or difficult-to-dissociate cell types and is free ...of both RNA degradation and artifactual transcriptional stress responses.
We compared single-cell RNA sequencing (scRNA-seq) using the DropSeq platform with single-nucleus RNA sequencing (snRNA-seq) using sNuc-DropSeq, DroNc-seq, and 10X Chromium platforms on adult mouse kidney. We validated snRNA-seq on fibrotic kidney from mice 14 days after unilateral ureteral obstruction (UUO) surgery.
A total of 11,391 transcriptomes were generated in the comparison phase. We identified ten clusters in the scRNA-seq dataset, but glomerular cell types were absent, and one cluster consisted primarily of artifactual dissociation
induced stress response genes. By contrast, snRNA-seq from all three platforms captured a diversity of kidney cell types that were not represented in the scRNA-seq dataset, including glomerular podocytes, mesangial cells, and endothelial cells. No stress response genes were detected. Our snRNA-seq protocol yielded 20-fold more podocytes compared with published scRNA-seq datasets (2.4% versus 0.12%, respectively). Unexpectedly, single-cell and single-nucleus platforms had equivalent gene detection sensitivity. For validation, analysis of frozen day 14 UUO kidney revealed rare juxtaglomerular cells, novel activated proximal tubule and fibroblast cell states, and previously unidentified tubulointerstitial signaling pathways.
snRNA-seq achieves comparable gene detection to scRNA-seq in adult kidney, and it also has substantial advantages, including reduced dissociation bias, compatibility with frozen samples, elimination of dissociation-induced transcriptional stress responses, and successful performance on inflamed fibrotic kidney.
The underlying cellular events driving kidney fibrogenesis and metabolic dysfunction are incompletely understood. Here, we employed single-cell combinatorial indexing RNA sequencing to analyze 24 ...mouse kidneys from two fibrosis models. We profiled 309,666 cells in one experiment, representing 50 cell types/states encompassing epithelial, endothelial, immune, and stromal populations. Single-cell analysis identified diverse injury states of the proximal tubule, including two distinct early-phase populations with dysregulated lipid and amino acid metabolism, respectively. Lipid metabolism was defective in the chronic phase but was transiently activated in the very early stages of ischemia-induced injury, where we discovered increased lipid deposition and increased fatty acid β-oxidation. Perilipin 2 was identified as a surface marker of intracellular lipid droplets, and its knockdown in vitro disrupted cell energy state maintenance during lipid accumulation. Surveying epithelial cells across nephron segments identified shared and unique injury responses. Stromal cells exhibited high heterogeneity and contributed to fibrogenesis by epithelial-stromal crosstalk.
Display omitted
•sci-RNA-seq3 transcriptionally profiles 309,666 cells from 24 kidneys without batch effects•Two injured proximal tubule cell states with distinct metabolic profiles revealed•Transiently activated lipid metabolism and PLIN2+ lipid droplets appear in early IRI•Nephron epithelia possess both shared and segment-specific injury and repair responses
Li et al. profile the full time courses of mouse kidney fibrogenesis using single-cell combinatorial indexing RNA sequencing. They describe diverse injury states of proximal tubular cells, including one cell state with enhanced lipid metabolism at an early phase of ischemia-induced injury. This single-cell atlas defines kidney epithelial injury responses in fibrosis.
Recent techniques for single-cell RNA sequencing (scRNA-seq) at high throughput are leading to profound new discoveries in biology. The ability to generate vast amounts of transcriptomic data at ...cellular resolution represents a transformative advance, allowing the identification of novel cell types, states, and dynamics. In this review, we summarize the development of scRNA-seq methodologies and highlight their advantages and drawbacks. We discuss available software tools for analyzing scRNA-Seq data and summarize current computational challenges. Finally, we outline ways in which this powerful technology might be applied to discovery research in kidney development and disease.
Diabetic nephropathy is characterized by damage to both the glomerulus and tubulointerstitium, but relatively little is known about accompanying cell-specific changes in gene expression. We performed ...unbiased single-nucleus RNA sequencing (snRNA-seq) on cryopreserved human diabetic kidney samples to generate 23,980 single-nucleus transcriptomes from 3 control and 3 early diabetic nephropathy samples. All major cell types of the kidney were represented in the final dataset. Side-by-side comparison demonstrated cell-type–specific changes in gene expression that are important for ion transport, angiogenesis, and immune cell activation. In particular, we show that the diabetic thick ascending limb, late distal convoluted tubule, and principal cells all adopt a gene expression signature consistent with increased potassium secretion, including alterations in Na⁺/K⁺-ATPase, WNK1, mineralocorticoid receptor, and NEDD4L expression, as well as decreased paracellular calcium and magnesium reabsorption. We also identify strong angiogenic signatures in glomerular cell types, proximal convoluted tubule, distal convoluted tubule, and principal cells. Taken together, these results suggest that increased potassium secretion and angiogenic signaling represent early kidney responses in human diabetic nephropathy.
Maximizing the potential of human kidney organoids for drug testing and regenerative medicine and to model development and disease requires addressing cell immaturity, the lack of a mature collecting ...system, and off-target cell types. By independently generating two kidney progenitor cell populations—metanephric mesenchyme and ureteric bud (UB)-like cells—we could generate kidney organoids with a collecting system. We also identify the hormones aldosterone and arginine vasopressin (AVP) as critical to promote differentiation of collecting duct cell types including both principal cells (PCs) and intercalated cells (ICs). The resulting PCs express aquaporin-2 (AQP2) protein, which undergoes translocation to the apical membrane after vasopressin or forskolin stimulation. By single-cell RNA sequencing (scRNA-seq), we demonstrate improved proximal tubule maturation and reduced off-target cell populations. We also show appropriate downregulation of progenitor cell types, improved modeling of tubular injury, the presence of urothelium (Uro), and the ability of Notch pathway modulation to regulate PC:IC ratios during organoid development.
Display omitted
•Combining differentiated progenitors leads to kidney organoids with collecting duct•Aldosterone and vasopressin drive principal and intercalated cell differentiation•Organoids show improved maturation and model tubular injury•Notch regulates principal:intercalated cell ratios
Uchimura et al. combine independently differentiated metanephric mesenchyme-like and ureteric bud-like progenitors to generate human kidney organoids with a collecting system. Hormones aldosterone and arginine vasopressin drive principal and intercalated cell maturation, and Notch signaling could regulate cell ratios. Organoids also showed improved maturation and injury modeling.
Single-cell genomics techniques are revolutionizing our ability to characterize complex tissues. By contrast, the techniques used to analyze renal biopsy specimens have changed little over several ...decades. We tested the hypothesis that single-cell RNA-sequencing can comprehensively describe cell types and states in a human kidney biopsy specimen.
We generated 8746 single-cell transcriptomes from a healthy adult kidney and a single kidney transplant biopsy core by single-cell RNA-sequencing. Unsupervised clustering analysis of the biopsy specimen was performed to identify 16 distinct cell types, including all of the major immune cell types and most native kidney cell types, in this biopsy specimen, for which the histologic read was mixed rejection.
Monocytes formed two subclusters representing a nonclassical CD16+ group and a classic CD16- group expressing dendritic cell maturation markers. The presence of both monocyte cell subtypes was validated by staining of independent transplant biopsy specimens. Comparison of healthy kidney epithelial transcriptomes with biopsy specimen counterparts identified novel segment-specific proinflammatory responses in rejection. Endothelial cells formed three distinct subclusters: resting cells and two activated endothelial cell groups. One activated endothelial cell group expressed Fc receptor pathway activation and Ig internalization genes, consistent with the pathologic diagnosis of antibody-mediated rejection. We mapped previously defined genes that associate with rejection outcomes to single cell types and generated a searchable online gene expression database.
We present the first step toward incorporation of single-cell transcriptomics into kidney biopsy specimen interpretation, describe a heterogeneous immune response in mixed rejection, and provide a searchable resource for the scientific community.
Single-cell sequencing technologies have advanced our understanding of kidney biology and disease, but the loss of spatial information in these datasets hinders our interpretation of intercellular ...communication networks and regional gene expression patterns. New spatial transcriptomic sequencing platforms make it possible to measure the topography of gene expression at genome depth.
We optimized and validated a female bilateral ischemia-reperfusion injury model. Using the 10× Genomics Visium Spatial Gene Expression solution, we generated spatial maps of gene expression across the injury and repair time course, and applied two open-source computational tools, Giotto and SPOTlight, to increase resolution and measure cell-cell interaction dynamics.
An ischemia time of 34 minutes in a female murine model resulted in comparable injury to 22 minutes for males. We report a total of 16,856 unique genes mapped across our injury and repair time course. Giotto, a computational toolbox for spatial data analysis, enabled increased resolution mapping of genes and cell types. Using a seeded nonnegative matrix regression (SPOTlight) to deconvolute the dynamic landscape of cell-cell interactions, we found that injured proximal tubule cells were characterized by increasing macrophage and lymphocyte interactions even 6 weeks after injury, potentially reflecting the AKI to CKD transition.
In this transcriptomic atlas, we defined region-specific and injury-induced loss of differentiation markers and their re-expression during repair, as well as region-specific injury and repair transcriptional responses. Lastly, we created an interactive data visualization application for the scientific community to explore these results (http://humphreyslab.com/SingleCell/).
Diabetic kidney disease (DKD) occurs in ∼40% of patients with diabetes and causes kidney failure, cardiovascular disease, and premature death. We analyzed the response of a murine DKD model to five ...treatment regimens using single-cell RNA sequencing (scRNA-seq). Our atlas of ∼1 million cells revealed a heterogeneous response of all kidney cell types both to DKD and its treatment. Both monotherapy and combination therapies targeted differing cell types and induced distinct and non-overlapping transcriptional changes. The early effects of sodium-glucose cotransporter-2 inhibitors (SGLT2i) on the S1 segment of the proximal tubule suggest that this drug class induces fasting mimicry and hypoxia responses. Diabetes downregulated the spliceosome regulator serine/arginine-rich splicing factor 7 (Srsf7) in proximal tubule that was specifically rescued by SGLT2i. In vitro proximal tubule knockdown of Srsf7 induced a pro-inflammatory phenotype, implicating alternative splicing as a driver of DKD and suggesting SGLT2i regulation of proximal tubule alternative splicing as a potential mechanism of action for this drug class.
Abstract
The proximal tubule is a key regulator of kidney function and glucose metabolism. Diabetic kidney disease leads to proximal tubule injury and changes in chromatin accessibility that modify ...the activity of transcription factors involved in glucose metabolism and inflammation. Here we use single nucleus RNA and ATAC sequencing to show that diabetic kidney disease leads to reduced accessibility of glucocorticoid receptor binding sites and an injury-associated expression signature in the proximal tubule. We hypothesize that chromatin accessibility is regulated by genetic background and closely-intertwined with metabolic memory, which pre-programs the proximal tubule to respond differently to external stimuli. Glucocorticoid excess has long been known to increase risk for type 2 diabetes, which raises the possibility that glucocorticoid receptor inhibition may mitigate the adverse metabolic effects of diabetic kidney disease.
PDF
