Single-cell transcriptomics of the mouse kidney reveals potential cellular targets of kidney disease - PubMed (original) (raw)
Single-cell transcriptomics of the mouse kidney reveals potential cellular targets of kidney disease
Jihwan Park et al. Science. 2018.
Abstract
Our understanding of kidney disease pathogenesis is limited by an incomplete molecular characterization of the cell types responsible for the organ's multiple homeostatic functions. To help fill this knowledge gap, we characterized 57,979 cells from healthy mouse kidneys by using unbiased single-cell RNA sequencing. On the basis of gene expression patterns, we infer that inherited kidney diseases that arise from distinct genetic mutations but share the same phenotypic manifestation originate from the same differentiated cell type. We also found that the collecting duct in kidneys of adult mice generates a spectrum of cell types through a newly identified transitional cell. Computational cell trajectory analysis and in vivo lineage tracing revealed that intercalated cells and principal cells undergo transitions mediated by the Notch signaling pathway. In mouse and human kidney disease, these transitions were shifted toward a principal cell fate and were associated with metabolic acidosis.
Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.
Figures
Fig. 1.. Cell diversity in mouse kidney cells delineated by single cell transcriptomic analysis.
(A) Unsupervised clustering demonstrates 16 distinct cell types shown in a tSNE map. Left panels are subclusters of clusters 1, 3, and 7. Percentage of assigned cell types are summarized in the right panel. Endo: containing endothelial, vascular, and descending loop of Henle, Podo: podocyte, PT: proximal tubule, LOH: ascending loop of Henle, DCT: distal convoluted tubule, CD-PC: collecting duct principal cell, CD-IC: CD intercalated cell, CD-Trans: CD transitional cell, Fib: fibroblast, Macro: macrophage, Neutro: neutrophil, NK: natural killer cell. (B and C) Violin plots showing the expression level of representative marker genes across the 16 main clusters. Y-axis is log scale normalized read count. (C) Cluster1 (from panel A left) separates into endothelial cells, pericytes/vascular smooth muscle cells and descending loop of Henle (DLH). Cluster3 (proximal tubules) separates into S1, S2 and S3 segments or PCT (proximal convoluted tubules) and PST (proximal straight tubules). Cluster 7 intercalated cells separates into A and B type intercalated cells.
Fig. 2.. Discrete human disease phenotypes are due to mutations in single specific cell types
Single cell-type specific average expression of human (A) monogenic disease genes and (B) complex trait genes identified by GWAS. Mean expression values of the genes were calculated in each cluster. The color scheme is based on z-score distribution; the map shows genes with z-score>2. In the heatmap, each row represents one gene and each column is single cell type (defined in Fig. 1). The full list of cell types and genes are shown in fig. S11 and 12.
Fig. 3.. Identification of a transitional cell type and a conversion process in the collecting duct.
(A) The expression level of marker genes across the 16 clusters. Y-axis is log scale normalized read count. (B) Gene expression levels in PC (Aquaporin 2, Aqp2), IC (H+ATPase, Atp6v1g3,) and in Transitional cells (Syt7) demonstrated by a tSNE plot. (C) Representative immunofluorescence images of AQP2 (PC marker), ATP6V1B1 (IC marker) and DAPI in the kidney collecting duct. Arrow demonstrates transitory PC/IC cell type expressing AQP2 and ATP6V1B1. (D) Heatmap showing the expression level of differentially expressed genes in collecting duct cell types. Color scheme is based on z-score distribution. (E) Venn diagram showing the overlaps of differentially expressed genes between PC, IC and the novel cell type. (F) Immunofluorescence staining for PARM1 (IC specific) and AQP2 (upper panel) or ATP6V1B1 (lower panel) in the kidney collecting duct. A “double positive” cell is shown at the arrow (G) Ordering single cells along a cell conversion trajectory using Monocle. Three collecting duct cell clusters were used for ordering and plotted on low dimensional space with different colors. The tSNE plots next to the trajectory map show differentially expressed genes in the corresponding cell lineages. (H) Aqp2CremT/mG mouse model used for lineage tracing of AQP2 positive cells. Immunofluorescence staining for GFP, ATP6V1B1 and AQP2. Right panel is the quantification of GFP positive cells (mean ± SD). n=3. Note that AQP2 driven GFP (white) is found in PC (red + white) and IC (green + white) and in transitional cells (not shown) (I) Atp6CremT/mG mouse model used for lineage tracing of ATP6ase positive cells. Immunofluorescence staining for GFP, ATP6V1B1 and AQP2 in Atp6CremT/mG mice. Note that ATP6V1B1 driven GFP (white) is found in PC (red+ white) and IC (green+ white) and in transitional cells (red+ green+ white).
Fig. 4.. The IC-to- PC cell transition is driven by Notch ligand and receptor expression.
(A) Transcriptional profiles demonstrating the spectrum of expression of Notch genes in the collecting duct. Cells are ordered in pseudotime and color represents expression levels. (B) Double immunofluorescence staining for AQP2 (red) and JAG1 (green) in the kidney collecting duct. (C) Generation of mice with inducible expression of Notch (ICN1) in kidney tubules. Excess AQP2+ cells and conversely reciprocally decreased ATP6V1B1+ positive cells are found in Pax8rtTA/NICD tubules (mean ± SD). n=3. Asterisk represents significant difference; p-value < 0.01. (D) In silico deconvolution of mouse kidney bulk RNA profiling data. Wild type (WT) and Pax8rtTA/NICD samples were used for analysis. (E) Immunofluorescence quantification of cells labeled with AQP2 and ATP6V1B1 in control and FA mice (mean ± SD). n=3. Asterisk represents significant differences; p-value < 0.01. (F) In silico deconvolution of mouse kidney bulk RNA profiling. Control and kidney samples from FA injected mice were used for analysis. (G) Immunofluorescence staining for AQP2 and ATP6V1B1 in control, Pax8rtTA/NICD and FA collecting ducts. Note the abundance of AQP2+ cells in Pax8rtTA/NICD and FA and conversely the disappearance of ATP6V1B1+ cells (H) In silico deconvolution of bulk RNA profiling data from normal and CKD human samples (n=91). The histological fibrosis scores and HES1 expression levels for the corresponding samples are also shown. (I) Total serum bicarbonate level in control and FA induced kidney fibrosis (mean ± SD). n=5.
Comment in
- Whole-kidney single-cell transcriptomics identifies new cell types.
Heintze JM. Heintze JM. Nat Rev Nephrol. 2018 Jun;14(6):353. doi: 10.1038/s41581-018-0013-7. Nat Rev Nephrol. 2018. PMID: 29695753 No abstract available. - Mapping kidney cellular complexity.
Humphreys BD. Humphreys BD. Science. 2018 May 18;360(6390):709-710. doi: 10.1126/science.aat7271. Science. 2018. PMID: 29773732 No abstract available. - Single-Cell RNA Sequencing of the Adult Mouse Kidney: From Molecular Cataloging of Cell Types to Disease-Associated Predictions.
Lindström NO, De Sena Brandine G, Ransick A, McMahon AP. Lindström NO, et al. Am J Kidney Dis. 2019 Jan;73(1):140-142. doi: 10.1053/j.ajkd.2018.07.002. Epub 2018 Sep 18. Am J Kidney Dis. 2019. PMID: 30241960 Free PMC article. No abstract available.
Similar articles
- Elf5 is a principal cell lineage specific transcription factor in the kidney that contributes to Aqp2 and Avpr2 gene expression.
Grassmeyer J, Mukherjee M, deRiso J, Hettinger C, Bailey M, Sinha S, Visvader JE, Zhao H, Fogarty E, Surendran K. Grassmeyer J, et al. Dev Biol. 2017 Apr 1;424(1):77-89. doi: 10.1016/j.ydbio.2017.02.007. Epub 2017 Feb 17. Dev Biol. 2017. PMID: 28215940 Free PMC article. - Transcriptomes of major renal collecting duct cell types in mouse identified by single-cell RNA-seq.
Chen L, Lee JW, Chou CL, Nair AV, Battistone MA, Păunescu TG, Merkulova M, Breton S, Verlander JW, Wall SM, Brown D, Burg MB, Knepper MA. Chen L, et al. Proc Natl Acad Sci U S A. 2017 Nov 14;114(46):E9989-E9998. doi: 10.1073/pnas.1710964114. Epub 2017 Oct 31. Proc Natl Acad Sci U S A. 2017. PMID: 29089413 Free PMC article. - RNA sequencing of kidney distal tubule cells reveals multiple mediators of chronic aldosterone action.
Poulsen SB, Limbutara K, Fenton RA, Pisitkun T, Christensen BM. Poulsen SB, et al. Physiol Genomics. 2018 May 1;50(5):343-354. doi: 10.1152/physiolgenomics.00084.2017. Epub 2018 Mar 9. Physiol Genomics. 2018. PMID: 29521601 Free PMC article. - Single-Cell Transcriptomics Meets Lineage Tracing.
Kester L, van Oudenaarden A. Kester L, et al. Cell Stem Cell. 2018 Aug 2;23(2):166-179. doi: 10.1016/j.stem.2018.04.014. Epub 2018 May 10. Cell Stem Cell. 2018. PMID: 29754780 Review. - Single-cell RNA sequencing for the study of development, physiology and disease.
Potter SS. Potter SS. Nat Rev Nephrol. 2018 Aug;14(8):479-492. doi: 10.1038/s41581-018-0021-7. Nat Rev Nephrol. 2018. PMID: 29789704 Free PMC article. Review.
Cited by
- Super-resolved highly multiplexed immunofluorescence imaging for precise protein localization and podocyte ultrastructure.
Siegerist F, Kitzel S, Telli N, Dikou JS, Drenić V, Chadjichristos CE, Chatziantoniou C, Endlich N. Siegerist F, et al. J Cell Mol Med. 2024 Sep;28(18):e70066. doi: 10.1111/jcmm.70066. J Cell Mol Med. 2024. PMID: 39334561 Free PMC article. - Introducing synthetic thermostable RNase inhibitors to single-cell RNA-seq.
Noble JC, Lentini A, Hagemann-Jensen M, Sandberg R, Reinius B. Noble JC, et al. Nat Commun. 2024 Sep 27;15(1):8373. doi: 10.1038/s41467-024-52717-4. Nat Commun. 2024. PMID: 39333520 Free PMC article. - Nucleoporin-associated steroid-resistant nephrotic syndrome.
Yao L, Li Y, Wang P, Xu C, Yu Z. Yao L, et al. Pediatr Nephrol. 2024 Sep 27. doi: 10.1007/s00467-024-06494-3. Online ahead of print. Pediatr Nephrol. 2024. PMID: 39331077 Review. - nsDCC: dual-level contrastive clustering with nonuniform sampling for scRNA-seq data analysis.
Wang L, Li W, Zhou F, Yu K, Feng C, Zhao D. Wang L, et al. Brief Bioinform. 2024 Sep 23;25(6):bbae477. doi: 10.1093/bib/bbae477. Brief Bioinform. 2024. PMID: 39327063 Free PMC article. - Epigenomic, cistromic, and transcriptomic profiling of primary kidney tubular cells.
Liu Z, Zhang L, Chen Y. Liu Z, et al. J Biol Methods. 2024 Jul 10;11(2):e99010015. doi: 10.14440/jbm.2024.0009. eCollection 2024. J Biol Methods. 2024. PMID: 39323486 Free PMC article.
References
- Kriz W, Bankir L, A standard nomenclature for structures of the kidney. The Renal Commission of the International Union of Physiological Sciences (IUPS). Kidney international 33, 1–7 (1988). - PubMed
- Rozenblatt-Rosen O, Stubbington MJT, Regev A, Teichmann SA, The Human Cell Atlas: from vision to reality. Nature 550, 451–453 (2017). - PubMed
Publication types
MeSH terms
Substances
Grants and funding
- DP3 DK108220/DK/NIDDK NIH HHS/United States
- R01 DK087635/DK/NIDDK NIH HHS/United States
- R01 DK092684/DK/NIDDK NIH HHS/United States
- R01 DK073462/DK/NIDDK NIH HHS/United States
- U54 DK104309/DK/NIDDK NIH HHS/United States
- UG3 DK114926/DK/NIDDK NIH HHS/United States
- R01 DK076077/DK/NIDDK NIH HHS/United States
- R01 DK105821/DK/NIDDK NIH HHS/United States
LinkOut - more resources
Full Text Sources
Other Literature Sources
Medical
Miscellaneous