Single-cell RNA-seq is a powerful technique. Nevertheless, there are important limitations, including the technical challenges of breaking down an organ or tissue into a single-cell suspension. ...Invariably, this has required enzymatic incubation at 37°C, which can be expected to result in artifactual changes in gene expression patterns. Here, we describe a dissociation method that uses a protease with high activity in the cold, purified from a psychrophilic microorganism. The entire procedure is carried out at 6°C or colder, at which temperature mammalian transcriptional machinery is largely inactive, thereby effectively 'freezing in' the
gene expression patterns. To test this method, we carried out RNA-seq on 20,424 single cells from postnatal day 1 mouse kidneys, comparing the results of the psychrophilic protease method with procedures using 37°C incubation. We show that the cold protease method provides a great reduction in gene expression artifacts. In addition, the results produce a single-cell resolution gene expression atlas of the newborn mouse kidney, an interesting time in development when mature nephrons are present yet nephrogenesis remains extremely active.
Focal segmental glomerulosclerosis is a major cause of end stage renal disease. Many patients prove unresponsive to available therapies. An improved understanding of the molecular basis of the ...disease process could provide insights leading to novel therapeutic approaches. In this study we carried out an RNA-seq analysis of the altered gene expression patterns of podocytes, mesangial cells and glomerular endothelial cells of the bigenic Cd2ap+/-, Fyn-/- mutant mouse model of FSGS. In the podocytes we observed upregulation of many genes related to the Tgfβ family/pathway, including Gdnf, Tgfβ1, Tgfβ2, Snai2, Vegfb, Bmp4, and Tnc. The mutant podocytes also showed upregulation of Acta2, a marker of smooth muscle and associated with myofibroblasts, which are implicated in driving fibrosis. GO analysis of the podocyte upregulated genes identified elevated protein kinase activity, increased expression of growth factors, and negative regulation of cell adhesion, perhaps related to the observed podocyte loss. Both podocytes and mesangial cells showed strong upregulation of aldehyde dehydrogenase genes involved in the synthesis of retinoic acid. Similarly, the Cd2ap+/-, Fyn-/- mesangial cells, as well as podocytes in other genetic models, and the glomeruli of human FSGS patients, all show upregulation of the serine protease Prss23, with the common thread suggesting important functionality. Another gene with strong upregulation in the Cd2ap+/-, Fyn-/- mutant mesangial cells as well as multiple other mutant mouse models of FSGS was thrombospondin, which activates the secreted inactive form of Tgfβ. The Cd2ap+/-, Fyn-/- mutant endothelial cells showed elevated expression of genes involved in cell proliferation, angioblast migration, angiogenesis, and neovasculature, all consistent with the formation of new blood vessels in the diseased glomerulus. The resulting global definition of the perturbed molecular pathways in the three major cell types of the mutant glomerulus provide deeper understanding of the molecular pathogenic pathways.
The developing kidney provides a useful model for study of the principles of organogenesis. In this report we use three independent platforms, Drop-Seq, Chromium 10x Genomics and Fluidigm C1, to ...carry out single cell RNA-Seq (scRNA-Seq) analysis of the E14.5 mouse kidney. Using the software AltAnalyze, in conjunction with the unsupervised approach ICGS, we were unable to identify and confirm the presence of 16 distinct cell populations during this stage of active nephrogenesis. Using a novel integrative supervised computational strategy, we were able to successfully harmonize and compare the cell profiles across all three technological platforms. Analysis of possible cross compartment receptor/ligand interactions identified the nephrogenic zone stroma as a source of GDNF. This was unexpected because the cap mesenchyme nephron progenitors had been thought to be the sole source of GDNF, which is a key driver of branching morphogenesis of the collecting duct system. The expression of Gdnf by stromal cells was validated in several ways, including Gdnf in situ hybridization combined with immunohistochemistry for SIX2, and marker of nephron progenitors, and MEIS1, a marker of stromal cells. Finally, the single cell gene expression profiles generated in this study confirmed and extended previous work showing the presence of multilineage priming during kidney development. Nephron progenitors showed stochastic expression of genes associated with multiple potential differentiation lineages.
•Single cell RNA-seq of over 8000 E14.5 developing kidney cells.•Stromal cells, as well as cap mesenchyme, express GDNF.•Comparison of 10X Genomics, Drop-Seq and Fluidigm 800 cell IFC single cell RNA-Seq.•Multilineage priming at many stages of nephrogenesis.
Current management of AKI, a potentially fatal disorder that can also initiate or exacerbate CKD, is merely supportive. Therefore, deeper understanding of the molecular pathways perturbed in AKI is ...needed to identify targets with potential to lead to improved treatment.
We performed single-cell RNA sequencing (scRNA-seq) with the clinically relevant unilateral ischemia-reperfusion murine model of AKI at days 1, 2, 4, 7, 11, and 14 after AKI onset. Using real-time quantitative PCR, immunofluorescence, Western blotting, and both chromogenic and single-molecule
hybridizations, we validated AKI signatures in multiple experiments.
Our findings show the time course of changing gene expression patterns for multiple AKI stages and all renal cell types. We observed elevated expression of crucial injury response factors-including kidney injury molecule-1 (Kim1), lipocalin 2 (Lcn2), and keratin 8 (Krt8)-and of several novel genes (
,
, and
) not previously examined in kidney pathologies. AKI induced proximal tubule dedifferentiation, with a pronounced nephrogenic signature represented by
and
. Moreover, AKI caused the formation of "mixed-identity cells" (expressing markers of different renal cell types) that are normally seen only during early kidney development. The injured tubules acquired a proinflammatory and profibrotic phenotype; moreover, AKI dramatically modified ligand-receptor crosstalk, with potential pathologic epithelial-to-stromal interactions. Advancing age in AKI onset was associated with maladaptive response and kidney fibrosis.
The scRNA-seq, comprehensive, cell-specific profiles provide a valuable resource for examining molecular pathways that are perturbed in AKI. The results fully define AKI-associated dedifferentiation programs, potential pathologic ligand-receptor crosstalk, novel genes, and the improved injury response in younger mice, and highlight potential targets of kidney injury.
The uterus is a remarkable organ that must guard against infections while maintaining the ability to support growth of a fetus without rejection. The Hoxa10 and Hoxa11 genes have previously been ...shown to play essential roles in uterus development and function. In this report we show that the Hoxa9,10,11, Hoxc9,10,11, Hoxd9,10,11 genes play a redundant role in the formation of uterine glands. In addition, we use single cell RNA-seq to create a high resolution gene expression atlas of the developing wild type mouse uterus. Cell types and subtypes are defined, for example dividing endothelial cells into arterial, venous, capillary, and lymphatic, while epithelial cells separate into luminal and glandular subtypes. Further, a surprising heterogeneity of stromal and myocyte cell types are identified. Transcription factor codes and ligand/receptor interactions are characterized. We also used single cell RNA-seq to globally define the altered gene expression patterns in all developing uterus cell types for two Hox mutants, with 8 or 9 mutant Hox genes. The mutants show a striking disruption of Wnt signaling as well as the Cxcl12/Cxcr4 ligand/receptor axis.
The NIH FACEBASE consortium was established in part to create a central resource for craniofacial researchers. One purpose is to provide a molecular anatomy of craniofacial development. To this end ...we have used a combination of laser capture microdissection and RNA-Seq to define the gene expression programs driving development of the murine palate. We focused on the E14.5 palate, soon after medial fusion of the two palatal shelves. The palate was divided into multiple compartments, including both medial and lateral, as well as oral and nasal, for both the anterior and posterior domains. A total of 25 RNA-Seq datasets were generated. The results provide a comprehensive view of the region specific expression of all transcription factors, growth factors and receptors. Paracrine interactions can be inferred from flanking compartment growth factor/receptor expression patterns. The results are validated primarily through very high concordance with extensive previously published gene expression data for the developing palate. In addition selected immunostain validations were carried out. In conclusion, this report provides an RNA-Seq based atlas of gene expression patterns driving palate development at microanatomic resolution. This FACEBASE resource is designed to promote discovery by the craniofacial research community.
The advent of next-generation sequencing has introduced new opportunities in analyzing gene expression data. Research in systems biology has taken advantage of these opportunities by gleaning ...insights into gene regulatory networks through the analysis of gene association networks. Contrasting networks from different populations can reveal the many different roles genes fill, which can lead to new discoveries in gene function. Pathologies can also arise from aberrations in these gene-gene interactions. Exposing these network irregularities provides a new avenue for understanding and treating diseases. A general framework for integrating known gene regulatory pathways into a differential network analysis between two populations is proposed. The framework importantly allows for any gene-gene association measure to be used, and inference is carried out through permutation testing. A simulation study investigates the performance in identifying differentially connected genes when incorporating known pathways, even if the pathway knowledge is partially inaccurate. Another simulation study compares the general framework with four state-of-the-art methods. Two RNA-seq datasets are analyzed to illustrate the use of this framework in practice. In both examples, the analysis reveals genes and pathways that are known to be biologically significant along with potentially novel findings that may be used to motivate future research.
Progenitor heterogeneity and identities underlying tumor initiation and relapse in medulloblastomas remain elusive. Utilizing single-cell transcriptomic analysis, we demonstrated a developmental ...hierarchy of progenitor pools in Sonic Hedgehog (SHH) medulloblastomas, and identified OLIG2-expressing glial progenitors as transit-amplifying cells at the tumorigenic onset. Although OLIG2+ progenitors become quiescent stem-like cells in full-blown tumors, they are highly enriched in therapy-resistant and recurrent medulloblastomas. Depletion of mitotic Olig2+ progenitors or Olig2 ablation impeded tumor initiation. Genomic profiling revealed that OLIG2 modulates chromatin landscapes and activates oncogenic networks including HIPPO-YAP/TAZ and AURORA-A/MYCN pathways. Co-targeting these oncogenic pathways induced tumor growth arrest. Together, our results indicate that glial lineage-associated OLIG2+ progenitors are tumor-initiating cells during medulloblastoma tumorigenesis and relapse, suggesting OLIG2-driven oncogenic networks as potential therapeutic targets.
Display omitted
•Single-cell-omics reveal OLIG2+ glial progenitors as tumor-initiating cells in MB•OLIG2+ cells are quiescent stem-like in full-blown MB but re-emerge during relapse•Ablation of mitotic OLIG2+ cells or deletion of Olig2 impedes MB tumor initiation•OLIG2 activates HIPPO-YAP and AURORA-A/MYCN oncogenic networks to promote MB growth
Zhang et al. show that OLIG2+ progenitor cells drive Sonic Hedgehog subtype medulloblastoma (SHH-MB) tumorigenesis in mice and are enriched in recurrent or resistant patient SHH-MB. OLIG2 activates HIPPO and AURORA-A/MYCN signaling, and co-inhibition of these pathways arrests SHH-MB growth in mice.
Single-cell transcriptomes from dissociated tissues provide insights into cell types and their gene expression and may harbor additional information on spatial position and the local ...microenvironment. The kidney's cells are embedded into a gradient of increasing tissue osmolality from the cortex to the medulla, which may alter their transcriptomes and provide cues for spatial reconstruction.
Single-cell or single-nuclei mRNA sequencing of dissociated mouse kidneys and of dissected cortex, outer, and inner medulla, to represent the corticomedullary axis, was performed. Computational approaches predicted the spatial ordering of cells along the corticomedullary axis and quantitated expression levels of osmo-responsive genes.
hybridization validated computational predictions of spatial gene-expression patterns. The strategy was used to compare single-cell transcriptomes from wild-type mice to those of mice with a collecting duct-specific knockout of the transcription factor grainyhead-like 2 (Grhl2
), which display reduced renal medullary osmolality.
Single-cell transcriptomics from dissociated kidneys provided sufficient information to approximately reconstruct the spatial position of kidney tubule cells and to predict corticomedullary gene expression. Spatial gene expression in the kidney changes gradually and osmo-responsive genes follow the physiologic corticomedullary gradient of tissue osmolality. Single-nuclei transcriptomes from Grhl2
mice indicated a flattened expression gradient of osmo-responsive genes compared with control mice, consistent with their physiologic phenotype.
Single-cell transcriptomics from dissociated kidneys facilitated the prediction of spatial gene expression along the corticomedullary axis and quantitation of osmotically regulated genes, allowing the prediction of a physiologic phenotype.