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.

Significant heterogeneities in gene expression among individual cells are typically interrogated using single whole cell approaches. However, tissues that have highly interconnected processes, such as in the brain, present unique challenges. Single-nucleus RNA sequencing (SNS) has emerged as an alternative method of assessing a cell's transcriptome through the use of isolated nuclei. However, studies directly comparing expression data between nuclei and whole cells are lacking. Here, we have characterized nuclear and whole cell transcriptomes in mouse single neurons and provided a normalization strategy to reduce method-specific differences related to the length of genic regions. We confirmed a high concordance between nuclear and whole cell transcriptomes in the expression of cell type and metabolic modeling markers, but less so for a subset of genes associated with mitochondrial respiration. Therefore, our results indicate that single-nucleus transcriptome sequencing provides an effective means to profile cell type expression dynamics in previously inaccessible tissues.

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.

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.

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.

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.

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 developmental origins of most adult stem cells are poorly understood. Here, we report the identification of a transcription factor-RHOX10-critical for the initial establishment of spermatogonial stem cells (SSCs). Conditional loss of the entire 33-gene X-linked homeobox gene cluster that includes Rhox10 causes progressive spermatogenic decline, a phenotype indistinguishable from that caused by loss of only Rhox10. We demonstrate that this phenotype results from dramatically reduced SSC generation. By using a battery of approaches, including single-cell-RNA sequencing (scRNA-seq) analysis, we show that Rhox10 drives SSC generation by promoting pro-spermatogonia differentiation. Rhox10 also regulates batteries of migration genes and promotes the migration of pro-spermatogonia into the SSC niche. The identification of an X-linked homeobox gene that drives the initial generation of SSCs has implications for the evolution of X-linked gene clusters and sheds light on regulatory mechanisms influencing adult stem cell generation in general.

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 human brain has enormously complex cellular diversity and connectivities fundamental to our neural functions, yet difficulties in interrogating individual neurons has impeded understanding of the underlying transcriptional landscape. We developed a scalable approach to sequence and quantify RNA molecules in isolated neuronal nuclei from a postmortem brain, generating 3227 sets of single-neuron data from six distinct regions of the cerebral cortex. Using an iterative clustering and classification approach, we identified 16 neuronal subtypes that were further annotated on the basis of known markers and cortical cytoarchitecture. These data demonstrate a robust and scalable method for identifying and categorizing single nuclear transcriptomes, revealing shared genes sufficient to distinguish previously unknown and orthologous neuronal subtypes as well as regional identity and transcriptomic heterogeneity within the human brain.