Urinary Single-Cell Profiling Captures the Cellular Diversity of the Kidney

Amin Abedini, Yuan O Zhu, Shatakshee Chatterjee, Gabor Halasz, Kishor Devalaraja-Narashimha, Rojesh Shrestha, Michael S Balzer, Jihwan Park, Tong Zhou, Ziyuan Ma, Katie Marie Sullivan, Hailong Hu, Xin Sheng, Hongbo Liu, Yi Wei, Carine M Boustany-Kari, Uptal Patel, Salem Almaani, Matthew Palmer, Raymond Townsend, Shira Blady, Jonathan Hogan, Lori Morton, Katalin Susztak, TRIDENT Study Investigators, Katalin Susztak, Raymond Townsend, Shira Blady, Matthew Palmer, Carine Boustany, Richard Urquhart, Paolo Guarnieri, Lea Sarov-Blat, Erding Hu, Lori Morton, Kishor Devalaraja, Uptal Patel, Shawn Badal, John Liles, Jonathan Rosen, Anil Karihaloo, Randy Luciano, Jonathan Hogan, Amy Mottl, Shweta Bansal, Salem Almaani, Christos Argyropoulos, Kirk Campbell, Tamara Isakova, Oliver Lenz, Harold Szerlip, Matthias Kretzler, Pietro Canetta, Jeffery Schelling, Rupali Avasare, Frank Brosius, Michael Ross, Nelson Kopyt, James Tumlin, Julia Scialla, Richard Lafayette, Manisha Singh, Yan Zhong, Amin Abedini, Yuan O Zhu, Shatakshee Chatterjee, Gabor Halasz, Kishor Devalaraja-Narashimha, Rojesh Shrestha, Michael S Balzer, Jihwan Park, Tong Zhou, Ziyuan Ma, Katie Marie Sullivan, Hailong Hu, Xin Sheng, Hongbo Liu, Yi Wei, Carine M Boustany-Kari, Uptal Patel, Salem Almaani, Matthew Palmer, Raymond Townsend, Shira Blady, Jonathan Hogan, Lori Morton, Katalin Susztak, TRIDENT Study Investigators, Katalin Susztak, Raymond Townsend, Shira Blady, Matthew Palmer, Carine Boustany, Richard Urquhart, Paolo Guarnieri, Lea Sarov-Blat, Erding Hu, Lori Morton, Kishor Devalaraja, Uptal Patel, Shawn Badal, John Liles, Jonathan Rosen, Anil Karihaloo, Randy Luciano, Jonathan Hogan, Amy Mottl, Shweta Bansal, Salem Almaani, Christos Argyropoulos, Kirk Campbell, Tamara Isakova, Oliver Lenz, Harold Szerlip, Matthias Kretzler, Pietro Canetta, Jeffery Schelling, Rupali Avasare, Frank Brosius, Michael Ross, Nelson Kopyt, James Tumlin, Julia Scialla, Richard Lafayette, Manisha Singh, Yan Zhong

Abstract

Background: Microscopic analysis of urine sediment is probably the most commonly used diagnostic procedure in nephrology. The urinary cells, however, have not yet undergone careful unbiased characterization.

Methods: Single-cell transcriptomic analysis was performed on 17 urine samples obtained from five subjects at two different occasions, using both spot and 24-hour urine collection. A pooled urine sample from multiple healthy individuals served as a reference control. In total 23,082 cells were analyzed. Urinary cells were compared with human kidney and human bladder datasets to understand similarities and differences among the observed cell types.

Results: Almost all kidney cell types can be identified in urine, such as podocyte, proximal tubule, loop of Henle, and collecting duct, in addition to macrophages, lymphocytes, and bladder cells. The urinary cell-type composition was subject specific and reasonably stable using different collection methods and over time. Urinary cells clustered with kidney and bladder cells, such as urinary podocytes with kidney podocytes, and principal cells of the kidney and urine, indicating their similarities in gene expression.

Conclusions: A reference dataset for cells in human urine was generated. Single-cell transcriptomics enables detection and quantification of almost all types of cells in the kidney and urinary tract.

Keywords: RNA sequencing; diabetic kidney disease; single-cell transcriptomics; urine.

Copyright © 2021 by the American Society of Nephrology.

Figures

Figure 1.
Figure 1.
Study design and urine characteristics. (A) Subjects with diabetes underwent diagnostic kidney biopsy. Tissue samples analyzed and scored by light and electron microscopy and immunofluorescence. Subjects were followed for 18 months. Spot and 24-hour urine samples were collected 1 month apart and analyzed. Ten spot urine samples from ten healthy individuals as control were obtained and pooled for the single-cell RNA-seq study. (B) Sample level, urinary cell number per ml urine (y-axis), x-axis is each patient, collection method, and time. (C) Sample level, urinary cell viability (y-axis), x-axis is each patient, collection method, and time. (D) Total captured cell number per sample, each column represents one sample, lighter and darker shades of the same color represent spot urine sample and the darker color shows the 24-hour urine collection. V1 and V2 represent visit 1 and visit 2, respectively, for each subject. (E) Images of periodic acid–Schiff-stained kidney section of each kidney biopsy. Scale bars are shown at the left bottom of each image represent 200 µm. DM, diabetes mellitus; F/U, follow up.
Figure 2.
Figure 2.
Single-cell survey of the human urine. (A) Uniform Manifold Approximation and Projection (UMAP) dimension reduction of 23,089 urinary cells identified in 17 DKD and one pooled control urine samples. Clusters 0, 2, 8, 9, 10, 11, 15, 16, 19, and 20 were called Epi. (B) Bubble dot plots of the top cluster specific genes. The size of the dot indicates expression percentage and the darkness of the color indicates average expression. Epi, variety epithelial cells; Epi-PLAT+, PLAT-positive cells; Epi-KRT1+, KRT1-positive cells; umbrella, umbrella cells; endo, endothelial cells; podo, podocytes; LOH, loop of Henle; fibro, fibroblast; CD, collecting duct principal cell; macro, macrophages; lympho, lymphocytes; mes, mesenchymal cells.
Figure 3.
Figure 3.
Distribution of cells in control and DKD samples. (A) Subject-level feature plot of cells over the Uniform Manifold Approximation and Projection (UMAP) of all cells (gray). All cells analyzed for a single subject (visit 1 [V1] and V2, spot and 24-hour) were combined and shown in red. (B) The fractions of cells in each cell cluster in control and DKD urine samples. Darker color indicates higher percentage of cells. (C) The fractions of cells in each cell cluster in each sample. Darker color indicates a higher percentage of cells. (D) The stability of urinary cell number on the basis of the collection method. The collection method (spot versus 24-hour) level feature models (colored red) over all cells in the UMAP (gray) (upper panel). Cumulative urine cell fractions of spot and 24-hour urine samples (the two visits were combined) (lower panel). (E) The stability of urinary cell number on the basis of time. V1 or V2 level feature models (colored red) over all cells in the UMAP (gray) (upper panel). Cumulative urine cell fractions of visits 1 and 2 samples (combined spot and 24-hour collection) (lower panel). Epi, variety epithelial cells; Epi-PLAT+, PLAT-positive cells; Epi-KRT1+, KRT1-positive cells; LOH, loop of Henle; CD, collecting duct principal cell.
Figure 4.
Figure 4.
Integration of urine single cell with human kidney single nuclei datasets. (A) Uniform Manifold Approximation and Projection (UMAP) of Harmony-based integration of urinary and kidney cells. (B) UMAP of Harmony-based integration of urinary and kidney cells colored by the sample of origin. Blue indicates the urine origin and pink shows the cells originated from kidney. Each cluster is labelled with the cells of origin. (C) The percent of cells in each integrated (urine and kidney) cluster (y-axis) that came from each original cluster (x-axis). (D) Bubble dot plots of the top cell-type–specific differentially expressed genes in the integrated clusters of urine and kidney samples. The size of the dot indicates the expression percentage and the darkness of the color indicates average expression. (E) The fraction of cells in each integrated cluster and their sample of origin. The data are colored by the percentage of cell. Epi, variety epithelial cells; Epi-PLAT+, PLAT-positive cells; Epi-KRT1+, KRT1-positive cells; umbrella, umbrella cells; endo, endothelial cells; EC, endothelial cells; LOH, loop of Henle; CD PC, collecting duct principal cells; CD ICA, collecting duct intercalated cells A; CD ICB, collecting duct intercalated cells B; DCT, distal convoluted tubule; mesangial, mesangial cells; podo, podocytes; PEC, parietal epithelial cells; macro, macrophages; lympho, lymphocytes; U, urine; K, kidney; leuko, leukocyte; mes, mesenchymal cells.
Figure 5.
Figure 5.
Integration of urine single-cell, kidney single-nucleus, and bladder single-cell datasets. (A) Uniform Manifold Approximation and Projection (UMAP) of Harmony-based integration of urinary, bladder, and kidney cells. (B) UMAP of Harmony-based integration of urinary bladder and kidney cells colored by the sample of origin. Blue, green, and pink indicates the urine, kidney, and bladder origins, respectively. The origins of the cells are written in each plot. (C) The percent of cells in each integrated (kidney, urine, and bladder) cluster (y-axis) that originated from the original clusters (x-axis). (D) Bubble dot plots of the top cell-type–specific differentially expressed genes in the integrated clusters of urine, kidney, and bladder samples. The size of the dot indicates expression percentage and the darkness of the color indicates average expression. (E) The fraction of cells in each integrated cluster and their sample of origin. Epi, variety epithelial cells; Epi-PLAT+, PLAT-positive cells; Epi-KRT1+, KRT1-positive cells; Epi-TNNT1+, TNNT1-positive cells; PEC, parietal epithelial cells; LOH, loop of Henle; DCT, distal convoluted tubule; CD PC, collecting duct principal cells; CD ICA, collecting duct intercalated cells A; CD ICB, collecting duct intercalated cells B; SMC, smooth muscle cell; basal, basal cells; inter, intermediate cells; umbrella, umbrella cells; endo, endothelial cells; podo, podocyte; fibro, fibroblast; SMC, smooth muscle cell; macro, macrophages; mono, monocyte; lympho, lymphocytes; U, urine; K, kidney; B, bladder; leuko, leukocyte.
Figure 6.
Figure 6.
Urinary single-cell profiling provides a read-out for kidney disease genes and drug targets. (A) Urinary single cell-type specific expression enrichment of monogenic nephrotic syndrome genes. (B) Urinary single cell-type–specific expression enrichment of CKD-GWAS–nominated genes. (C) Urinary single cell-type–specific expression enrichment genes related to nephrolithiasis. (D) Urinary single cell-type–specific expression enrichment of Food and Drug Administration–approved drug target genes. Mean gene expression values were calculated in each cell type cluster. The color scheme of the heatmap is on the basis of z score distribution. In each heatmap, x-axis represents the urinary cell clusters and y-axis shows the genes. Epi, variety epithelial cells; Epi-PLAT+, PLAT-positive cells; Epi-KRT1+, KRT1-positive cells; PT, proximal tubule; LOH, loop of Henle; CD, collecting duct principal cells.

References

    1. Cavanaugh C, Perazella MA: Urine sediment examination in the diagnosis and management of kidney disease: Core curriculum 2019. Am J Kidney Dis 73: 258–272, 2019
    1. Oliveira Arcolino F, Tort Piella A, Papadimitriou E, Bussolati B, Antonie DJ, Murray P, et al. .: Human urine as a noninvasive source of kidney cells. Stem Cells Int 2015: 362562, 2015
    1. Racusen LC, Fivush BA, Andersson H, Gahl WA: Culture of renal tubular cells from the urine of patients with nephropathic cystinosis. J Am Soc Nephrol 1: 1028–1033, 1991
    1. Dörrenhaus A, Müller JI, Golka K, Jedrusik P, Schulze H, Föllmann W: Cultures of exfoliated epithelial cells from different locations of the human urinary tract and the renal tubular system. Arch Toxicol 74: 618–626, 2000
    1. Inoue CN, Sunagawa N, Morimoto T, Ohnuma S, Katsushima F, Nishio T, et al. .: Reconstruction of tubular structures in three-dimensional collagen gel culture using proximal tubular epithelial cells voided in human urine. In Vitro Cell Dev Biol Anim 39: 364–367, 2003
    1. Al-Malki AL: Assessment of urinary osteopontin in association with podocyte for early predication of nephropathy in diabetic patients. Dis Markers 2014: 493736, 2014
    1. Nakamura T, Ushiyama C, Suzuki S, Hara M, Shimada N, Ebihara I, et al. .: Urinary excretion of podocytes in patients with diabetic nephropathy. Nephrol Dial Transplant 15: 1379–1383, 2000
    1. Detrisac CJ, Mayfield RK, Colwell JA, Garvin AJ, Sens DA: In vitro culture of cells exfoliated in the urine by patients with diabetes mellitus. J Clin Invest 71: 170–173, 1983
    1. Bharadwaj S, Liu G, Shi Y, Wu R, Yang B, He T, et al. .: Multipotential differentiation of human urine-derived stem cells: Potential for therapeutic applications in urology. Stem Cells 31: 1840–1856, 2013
    1. Lang R, Liu G, Shi Y, Bharadwaj S, Leng X, Zhou X, et al. .: Self-renewal and differentiation capacity of urine-derived stem cells after urine preservation for 24 hours. PLoS One 8: e53980, 2013
    1. Da Sacco S, Sedrakyan S, Boldrin F, Giuliani S, Parnigotto P, Habibian R, et al. .: Human amniotic fluid as a potential new source of organ specific precursor cells for future regenerative medicine applications. J Urol 183: 1193–1200, 2010
    1. Rahman MS, Wruck W, Spitzhorn L-S, Nguyen L, Bohndorf M, Martins S, et al. .: The fGf, tGfβ and Wnt axis modulate self-renewal of human SIX2+ urine derived renal progenitor cells. Sci Rep 10: 1–16, 2020
    1. Zhang Y, McNeill E, Tian H, Soker S, Andersson K-E, Yoo JJ, et al. .: Urine derived cells are a potential source for urological tissue reconstruction. J Urol 180: 2226–2233, 2008
    1. Arazi A, Rao DA, Berthier CC, Davidson A, Liu Y, Hoover PJ, et al. .; Accelerating Medicines Partnership in SLE network: The immune cell landscape in kidneys of patients with lupus nephritis [published correction appears in Nat Immunol 20: 1404, 2019 10.1038/s41590-019-0473-3]. Nat Immunol 20: 902–914, 2019
    1. Wang YJ, Kaestner KH: Single-cell RNA-seq of the pancreatic islets––a promise not yet fulfilled? Cell Metab 29: 539–544, 2019
    1. Park J, Liu CL, Kim J, Susztak K: Understanding the kidney one cell at a time. Kidney Int 96: 862–870, 2019
    1. Park J, Shrestha R, Qiu C, Kondo A, Huang S, Werth M, et al. .: Single-cell transcriptomics of the mouse kidney reveals potential cellular targets of kidney disease. Science 360: 758–763, 2018
    1. Menon R, Otto EA, Sealfon R, Nair V, Wong AK, Theesfeld CL, et al. .: SARS-CoV-2 receptor networks in diabetic and COVID-19-associated kidney disease. Kidney Int 98: 1502–1518, 2020
    1. Townsend RR, Guarnieri P, Argyropoulos C, Blady S, Boustany-Kari CM, Devalaraja-Narashimha K, et al. .; TRIDENT Study Investigators: Rationale and design of the transformative research in diabetic nephropathy (TRIDENT) study [published correction appears in Kidney Int 97: 809, 2020 10.1016/j.kint.2020.02.005]. Kidney Int 97: 10–13, 2020
    1. Butler A, Hoffman P, Smibert P, Papalexi E, Satija R: Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol 36: 411–420, 2018
    1. Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, et al. .: Fast, sensitive and accurate integration of single-cell data with Harmony. Nat Methods 16: 1289–1296, 2019
    1. Young MD, Behjati S: SoupX removes ambient RNA contamination from droplet based single cell RNA sequencing data. BioRxiv 303727, 2020. 10.1101/303727
    1. McGinnis CS, Murrow LM, Gartner ZJ: DoubletFinder: Doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell Systems 8: 329–337.e4, 2019
    1. Wilson PC, Wu H, Kirita Y, Uchimura K, Ledru N, Rennke HG, et al. .: The single-cell transcriptomic landscape of early human diabetic nephropathy. Proc Natl Acad Sci U S A 116: 19619–19625, 2019
    1. Yu Z, Liao J, Chen Y, Zou C, Zhang H, Cheng J, et al. .: Single-cell transcriptomic map of the human and mouse bladders. J Am Soc Nephrol 30: 2159–2176, 2019
    1. Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM 3rd, et al. .: Comprehensive integration of single-cell data. Cell 177: 1888–1902.e21, 2019
    1. Wu H, Uchimura K, Donnelly EL, Kirita Y, Morris SA, Humphreys BD: Comparative analysis and refinement of human PSC-derived kidney organoid differentiation with single-cell transcriptomics. Cell Stem Cell 23: 869–881.e8, 2018
    1. Nacken W, Roth J, Sorg C, Kerkhoff C: S100A9/S100A8: Myeloid representatives of the S100 protein family as prominent players in innate immunity. Microsc Res Tech 60: 569–580, 2003
    1. Menon R, Otto EA, Hoover P, Eddy S, Mariani L, Godfrey B, et al. .; Nephrotic Syndrome Study Network (NEPTUNE): Single cell transcriptomics identifies focal segmental glomerulosclerosis remission endothelial biomarker. JCI Insight 5: e133267, 2020
    1. Dumas SJ, Meta E, Borri M, Goveia J, Rohlenova K, Conchinha NV, et al. .: Single-cell RNA sequencing reveals renal endothelium heterogeneity and metabolic adaptation to water deprivation. J Am Soc Nephrol 31: 118–138, 2020
    1. Pattaro C, Teumer A, Gorski M, Chu AY, Li M, Mijatovic V, et al. .; ICBP Consortium; AGEN Consortium; CARDIOGRAM; CHARGe-Heart Failure Group; ECHOGen Consortium: Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function. Nat Commun 7: 10023, 2016
    1. Qiu C, Huang S, Park J, Park Y, Ko Y-A, Seasock MJ, et al. .: Renal compartment–specific genetic variation analyses identify new pathways in chronic kidney disease. Nature Medicine 24: 1721–1731, 2018
    1. Wu H, Malone AF, Donnelly EL, Kirita Y, Uchimura K, Ramakrishnan SM, et al. .: Single-cell transcriptomics of a human kidney allograft biopsy specimen defines a diverse inflammatory response. J Am Soc Nephrol 29: 2069–2080, 2018
    1. Edeling M, Ragi S, Huang S, Pavenstädt H, Susztak K: Developmental signalling pathways in renal fibrosis: the roles of Notch, Wnt and Hedgehog. Nature Reviews Nephrology 12: 426, 2016

Source: PubMed

Patrocinadores y Colaboradores

Condiciones médicas

Intervenciones de drogas

3
Suscribir