An atlas of healthy and injured cell states and niches in the human kidney.

Understanding kidney disease relies on defining the complexity of cell types and states, their associated molecular profiles and interactions within tissue neighbourhoods1. Here we applied multiple single-cell and single-nucleus assays (>400,000 nuclei or cells) and spatial imaging technologies to a broad spectrum of healthy reference kidneys (45 donors) and diseased kidneys (48 patients). This has provided a high-resolution cellular atlas of 51 main cell types, which include rare and previously undescribed cell populations. The multi-omic approach provides detailed transcriptomic profiles, regulatory factors and spatial localizations spanning the entire kidney. We also define 28 cellular states across nephron segments and interstitium that were altered in kidney injury, encompassing cycling, adaptive (successful or maladaptive repair), transitioning and degenerative states. Molecular signatures permitted the localization of these states within injury neighbourhoods using spatial transcriptomics, while large-scale 3D imaging analysis (around 1.2 million neighbourhoods) provided corresponding linkages to active immune responses. These analyses defined biological pathways that are relevant to injury time-course and niches, including signatures underlying epithelial repair that predicted maladaptive states associated with a decline in kidney function. This integrated multimodal spatial cell atlas of healthy and diseased human kidneys represents a comprehensive benchmark of cellular states, neighbourhoods, outcome-associated signatures and publicly available interactive visualizations.

Comparative cellular analysis of motor cortex in human, marmoset and mouse.

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch-seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

A multimodal and integrated approach to interrogate human kidney biopsies with rigor and reproducibility: guidelines from the Kidney Precision Medicine Project.

Comprehensive and spatially mapped molecular atlases of organs at a cellular level are a critical resource to gain insights into pathogenic mechanisms and personalized therapies for diseases. The Kidney Precision Medicine Project (KPMP) is an endeavor to generate three-dimensional (3-D) molecular atlases of healthy and diseased kidney biopsies by using multiple state-of-the-art omics and imaging technologies across several institutions. Obtaining rigorous and reproducible results from disparate methods and at different sites to interrogate biomolecules at a single-cell level or in 3-D space is a significant challenge that can be a futile exercise if not well controlled. We describe a "follow the tissue" pipeline for generating a reliable and authentic single-cell/region 3-D molecular atlas of human adult kidney. Our approach emphasizes quality assurance, quality control, validation, and harmonization across different omics and imaging technologies from sample procurement, processing, storage, shipping to data generation, analysis, and sharing. We established benchmarks for quality control, rigor, reproducibility, and feasibility across multiple technologies through a pilot experiment using common source tissue that was processed and analyzed at different institutions and different technologies. A peer review system was established to critically review quality control measures and the reproducibility of data generated by each technology before their being approved to interrogate clinical biopsy specimens. The process established economizes the use of valuable biopsy tissue for multiomics and imaging analysis with stringent quality control to ensure rigor and reproducibility of results and serves as a model for precision medicine projects across laboratories, institutions and consortia.

Mapping human tissues with highly multiplexed RNA in situ hybridization.

In situ transcriptomic techniques promise a holistic view of tissue organization and cell-cell interactions. There has been a surge of multiplexed RNA in situ mapping techniques but their application to human tissues has been limited due to their large size, general lower tissue quality and high autofluorescence. Here we report DART-FISH, a padlock probe-based technology capable of profiling hundreds to thousands of genes in centimeter-sized human tissue sections. We introduce an omni-cell type cytoplasmic stain that substantially improves the segmentation of cell bodies. Our enzyme-free isothermal decoding procedure allows us to image 121 genes in large sections from the human neocortex in <10 h. We successfully recapitulated the cytoarchitecture of 20 neuronal and non-neuronal subclasses. We further performed in situ mapping of 300 genes on a diseased human kidney, profiled >20 healthy and pathological cell states, and identified diseased niches enriched in transcriptionally altered epithelial cells and myofibroblasts.

The chromatin landscape of healthy and injured cell types in the human kidney.

There is a need to define regions of gene activation or repression that control human kidney cells in states of health, injury, and repair to understand the molecular pathogenesis of kidney disease and design therapeutic strategies. Comprehensive integration of gene expression with epigenetic features that define regulatory elements remains a significant challenge. We measure dual single nucleus RNA expression and chromatin accessibility, DNA methylation, and H3K27ac, H3K4me1, H3K4me3, and H3K27me3 histone modifications to decipher the chromatin landscape and gene regulation of the kidney in reference and adaptive injury states. We establish a spatially-anchored epigenomic atlas to define the kidney's active, silent, and regulatory accessible chromatin regions across the genome. Using this atlas, we note distinct control of adaptive injury in different epithelial cell types. A proximal tubule cell transcription factor network of ELF3, KLF6, and KLF10 regulates the transition between health and injury, while in thick ascending limb cells this transition is regulated by NR2F1. Further, combined perturbation of ELF3, KLF6, and KLF10 distinguishes two adaptive proximal tubular cell subtypes, one of which manifested a repair trajectory after knockout. This atlas will serve as a foundation to facilitate targeted cell-specific therapeutics by reprogramming gene regulatory networks.

Context-dependent enhancement of induced pluripotent stem cell reprogramming by silencing Puma.

Reprogramming of the somatic state to pluripotency can be induced by a defined set of transcription factors including Oct3/4, Sox2, Klf4, and c-Myc [Cell 2006;126:663-676]. These induced pluripotent stem cells (iPSCs) hold great promise in human therapy and disease modeling. However, tumor suppressive activities of p53, which are necessary to prevent persistence of DNA damage in mammalian cells, have proven a serious impediment to formation of iPSCs [Nat Methods 2011;8:409-412]. We examined the requirement for downstream p53 activities in suppressing efficiency of reprogramming as well as preventing persistence of DNA damage into the early iPSCs. We discovered that the majority of the p53 activation occurred through early reprogramming-induced DNA damage with the activated expression of the apoptotic inducer Puma and the cell cycle inhibitor p21. While Puma deficiency increases reprogramming efficiency only in the absence of c-Myc, double deficiency of Puma and p21 has achieved a level of efficiency that exceeded that of p53 deficiency alone. We further demonstrated that, in both the presence and absence of p21, Puma deficiency was able to prevent any increase in persistent DNA damage in early iPSCs. This may be due to a compensatory cellular senescent response to reprogramming-induced DNA damage in pre-iPSCs. Therefore, our findings provide a potentially safe approach to enhance iPSC derivation by transiently silencing Puma and p21 without compromising genomic integrity.

Rab11 regulates planar polarity and migratory behavior of multiciliated cells in Xenopus embryonic epidermis.

BACKGROUND: Xenopus embryonic skin is composed of the superficial layer with defined apicobasal polarity and the inner layer lacking the apical domain. Multiciliated cells (MCCs) originate in the inner layer of the epidermal ectoderm and subsequently migrate to the surface. How MCCs acquire the apicobasal polarity and intercalate into the superficial layer during neurulation is largely unknown. As Rab11-dependent vesicle trafficking has been implicated in ciliary membrane assembly and in apical domain formation in epithelial cells, we assessed the involvement of Rab11 in MCC development. RESULTS: Here we report that Rab11 is specifically enriched and becomes apically polarized in skin MCCs. Interference with Rab11 function by overexpression of a dominant negative mutant or injection of a specific morpholino oligonucleotide inhibited MCC intercalation into the superficial layer. Dominant negative Rab11-expressing MCC precursors revealed intrinsic apicobasal polarity, characterized by the apical domain, which is not normally observed in inner layer cells. Despite the presence of the apical domain, the cells with inhibited Rab11 function were randomly oriented relative to the plane of the tissue, thereby demonstrating a defect in planar polarity. CONCLUSIONS: These results establish a requirement for Rab11 in MCC development and support a two-step model, in which the initial polarization of MCC precursors is critical for their integration into the superficial cell layer.

Mapping Human Tissues with Highly Multiplexed RNA in situ Hybridization.

In situ transcriptomic techniques promise a holistic view of tissue organization and cell-cell interactions. Recently there has been a surge of multiplexed RNA in situ techniques but their application to human tissues and clinical biopsies has been limited due to their large size, general lower tissue quality and high background autofluorescence. Here we report DART-FISH, a versatile padlock probe-based technology capable of profiling hundreds to thousands of genes in centimeter-sized human tissue sections at cellular resolution. We introduced an omni-cell type cytoplasmic stain, dubbed RiboSoma that substantially improves the segmentation of cell bodies. We developed a computational decoding-by-deconvolution workflow to extract gene spots even in the presence of optical crowding. Our enzyme-free isothermal decoding procedure allowed us to image 121 genes in a large section from the human neocortex in less than 10 hours, where we successfully recapitulated the cytoarchitecture of 20 neuronal and non-neuronal subclasses. Additionally, we demonstrated the detection of transcripts as short as 461 nucleotides, including neuropeptides and discovered new cortical layer markers. We further performed in situ mapping of 300 genes on a diseased human kidney, profiled >20 healthy and pathological cell states, and identified diseased niches enriched in transcriptionally altered epithelial cells and myofibroblasts.

A spatially anchored transcriptomic atlas of the human kidney papilla identifies significant immune injury in patients with stone disease.

Kidney stone disease causes significant morbidity and increases health care utilization. In this work, we decipher the cellular and molecular niche of the human renal papilla in patients with calcium oxalate (CaOx) stone disease and healthy subjects. In addition to identifying cell types important in papillary physiology, we characterize collecting duct cell subtypes and an undifferentiated epithelial cell type that was more prevalent in stone patients. Despite the focal nature of mineral deposition in nephrolithiasis, we uncover a global injury signature characterized by immune activation, oxidative stress and extracellular matrix remodeling. We also identify the association of MMP7 and MMP9 expression with stone disease and mineral deposition, respectively. MMP7 and MMP9 are significantly increased in the urine of patients with CaOx stone disease, and their levels correlate with disease activity. Our results define the spatial molecular landscape and specific pathways contributing to stone-mediated injury in the human papilla and identify associated urinary biomarkers.

The chromatin landscape of healthy and injured cell types in the human kidney.

There is a need to define regions of gene activation or repression that control human kidney cells in states of health, injury, and repair to understand the molecular pathogenesis of kidney disease and design therapeutic strategies. However, comprehensive integration of gene expression with epigenetic features that define regulatory elements remains a significant challenge. We measured dual single nucleus RNA expression and chromatin accessibility, DNA methylation, and H3K27ac, H3K4me1, H3K4me3, and H3K27me3 histone modifications to decipher the chromatin landscape and gene regulation of the kidney in reference and adaptive injury states. We established a comprehensive and spatially-anchored epigenomic atlas to define the kidney's active, silent, and regulatory accessible chromatin regions across the genome. Using this atlas, we noted distinct control of adaptive injury in different epithelial cell types. A proximal tubule cell transcription factor network of ELF3 , KLF6 , and KLF10 regulated the transition between health and injury, while in thick ascending limb cells this transition was regulated by NR2F1 . Further, combined perturbation of ELF3 , KLF6 , and KLF10 distinguished two adaptive proximal tubular cell subtypes, one of which manifested a repair trajectory after knockout. This atlas will serve as a foundation to facilitate targeted cell-specific therapeutics by reprogramming gene regulatory networks.

A reference tissue atlas for the human kidney.

Kidney Precision Medicine Project (KPMP) is building a spatially specified human kidney tissue atlas in health and disease with single-cell resolution. Here, we describe the construction of an integrated reference map of cells, pathways, and genes using unaffected regions of nephrectomy tissues and undiseased human biopsies from 56 adult subjects. We use single-cell/nucleus transcriptomics, subsegmental laser microdissection transcriptomics and proteomics, near-single-cell proteomics, 3D and CODEX imaging, and spatial metabolomics to hierarchically identify genes, pathways, and cells. Integrated data from these different technologies coherently identify cell types/subtypes within different nephron segments and the interstitium. These profiles describe cell-level functional organization of the kidney following its physiological functions and link cell subtypes to genes, proteins, metabolites, and pathways. They further show that messenger RNA levels along the nephron are congruent with the subsegmental physiological activity. This reference atlas provides a framework for the classification of kidney disease when multiple molecular mechanisms underlie convergent clinical phenotypes.

Scalable dual-omics profiling with single-nucleus chromatin accessibility and mRNA expression sequencing 2 (SNARE-seq2).

Comprehensive characterization of cellular heterogeneity and the underlying regulatory landscapes of tissues and organs requires a highly robust and scalable method to acquire matched RNA and chromatin accessibility profiles on the same cells. Here, we describe a single-nucleus chromatin accessibility and mRNA expression sequencing 2 (SNARE-seq2) assay, implemented with cellular combinatorial indexing. This method involves tagmentation within permeabilized and fixed single-nucleus isolates to capture accessible chromatin (AC) regions, followed by the capture and reverse transcription of RNA transcripts. Through combinatorial split pool ligations, cDNA and AC within each single nucleus become appended with a common cell barcode combination. The captured cDNA and AC are then co-amplified before splitting and enrichment into single-nucleus RNA and single-nucleus AC sequencing libraries. This protocol is compatible with both nuclei and whole cells and can be completed in 3.5 d. SNARE-seq2 permits robust generation of high-quality, joint single-cell RNA and AC sequencing libraries from hundreds of thousands of single cells per experiment.

The role of the NMD factor UPF3B in olfactory sensory neurons.

The UPF3B-dependent branch of the nonsense-mediated RNA decay (NMD) pathway is critical for human cognition. Here, we examined the role of UPF3B in the olfactory system. Single-cell RNA-sequencing (scRNA-seq) analysis demonstrated considerable heterogeneity of olfactory sensory neuron (OSN) cell populations in wild-type (WT) mice, and revealed that UPF3B loss influences specific subsets of these cell populations. UPF3B also regulates the expression of a large cadre of antimicrobial genes in OSNs, and promotes the selection of specific olfactory receptor (Olfr) genes for expression in mature OSNs (mOSNs). RNA-seq and Ribotag analyses identified classes of mRNAs expressed and translated at different levels in WT and Upf3b-null mOSNs. Integrating multiple computational approaches, UPF3B-dependent NMD target transcripts that are candidates to mediate the functions of NMD in mOSNs were identified in vivo. Together, our data provides a valuable resource for the olfactory field and insights into the roles of NMD in vivo.

High-throughput sequencing of the transcriptome and chromatin accessibility in the same cell.

Single-cell RNA sequencing can reveal the transcriptional state of cells, yet provides little insight into the upstream regulatory landscape associated with open or accessible chromatin regions. Joint profiling of accessible chromatin and RNA within the same cells would permit direct matching of transcriptional regulation to its outputs. Here, we describe droplet-based single-nucleus chromatin accessibility and mRNA expression sequencing (SNARE-seq), a method that can link a cell's transcriptome with its accessible chromatin for sequencing at scale. Specifically, accessible sites are captured by Tn5 transposase in permeabilized nuclei to permit, within many droplets in parallel, DNA barcode tagging together with the mRNA molecules from the same cells. To demonstrate the utility of SNARE-seq, we generated joint profiles of 5,081 and 10,309 cells from neonatal and adult mouse cerebral cortices, respectively. We reconstructed the transcriptome and epigenetic landscapes of major and rare cell types, uncovered lineage-specific accessible sites, especially for low-abundance cells, and connected the dynamics of promoter accessibility with transcription level during neurogenesis.

Cellular Recruitment by Podocyte-Derived Pro-migratory Factors in Assembly of the Human Renal Filter.

Analysis of kidney disease-causing genes and pathology resulting from systemic diseases highlight the importance of the kidney's filtering system, the renal corpuscles. To elucidate the developmental processes that establish the renal corpuscle, we performed single-nucleus droplet-based sequencing of the human fetal kidney. This enabled the identification of nephron, interstitial, and vascular cell types that together generate the renal corpuscles. Trajectory analysis identified transient developmental gene expression, predicting precursors or mature podocytes express FBLN2, BMP4, or NTN4, in conjunction with recruitment, differentiation, and modeling of vascular and mesangial cell types into a functional filter. In vitro studies provide evidence that these factors exhibit angiogenic or mesangial recruiting and inductive properties consistent with a key organizing role for podocyte precursors in kidney development. Together these studies define a spatiotemporal developmental program for the primary filtration unit of the human kidney and provide novel insights into cell interactions regulating co-assembly of constituent cell types.

A single-nucleus RNA-sequencing pipeline to decipher the molecular anatomy and pathophysiology of human kidneys.

Defining cellular and molecular identities within the kidney is necessary to understand its organization and function in health and disease. Here we demonstrate a reproducible method with minimal artifacts for single-nucleus Droplet-based RNA sequencing (snDrop-Seq) that we use to resolve thirty distinct cell populations in human adult kidney. We define molecular transition states along more than ten nephron segments spanning two major kidney regions. We further delineate cell type-specific expression of genes associated with chronic kidney disease, diabetes and hypertension, providing insight into possible targeted therapies. This includes expression of a hypertension-associated mechano-sensory ion channel in mesangial cells, and identification of proximal tubule cell populations defined by pathogenic expression signatures. Our fully optimized, quality-controlled transcriptomic profiling pipeline constitutes a tool for the generation of healthy and diseased molecular atlases applicable to clinical samples.

Integrative single-cell analysis of transcriptional and epigenetic states in the human adult brain.

Detailed characterization of the cell types in the human brain requires scalable experimental approaches to examine multiple aspects of the molecular state of individual cells, as well as computational integration of the data to produce unified cell-state annotations. Here we report improved high-throughput methods for single-nucleus droplet-based sequencing (snDrop-seq) and single-cell transposome hypersensitive site sequencing (scTHS-seq). We used each method to acquire nuclear transcriptomic and DNA accessibility maps for >60,000 single cells from human adult visual cortex, frontal cortex, and cerebellum. Integration of these data revealed regulatory elements and transcription factors that underlie cell-type distinctions, providing a basis for the study of complex processes in the brain, such as genetic programs that coordinate adult remyelination. We also mapped disease-associated risk variants to specific cellular populations, which provided insights into normal and pathogenic cellular processes in the human brain. This integrative multi-omics approach permits more detailed single-cell interrogation of complex organs and tissues.

Requirement of Pygopus 2 in breast cancer.

The development of novel therapeutic strategies for breast cancer requires the identification of molecular targets involved in malignancy. Human Pygopus (Pygo)-1 and -2 are recently discovered components of the Wnt signaling pathway required for beta-Catenin/Tcf dependent transcription in embryos and colorectal cancer cells, but the role of these proteins in malignant cell growth and survival has not yet been determined. We report the expression and requirement for proliferation of hPygo2 in breast cancer cells. hPygo2 protein was overexpressed in malignant breast tumors and in the nuclei of five breast cancer cell lines, but was not expressed in the nuclei of non-malignant breast cells. Phosphorothioated antisense oligonucleotides were used to specifically knockdown expression hPygo2 in Mcf-7 and MDA-MB-468 cell lines. hPygo2 was required for the growth, in tissue culture and anchorage-independent assays, of both cell lines and for the expression of the Wnt target gene Cyclin D1. We conclude that hPygo2 is highly expressed in, and required for the growth of breast carcinoma cells.