We used a single cell RNA-seq strategy to create an atlas of gene expression patterns in the developing kidney. At several stages of kidney development, histologically uniform populations of cells give rise to multiple distinct lineages. We performed single cell RNA-seq analysis of total mouse kidneys at E11.5 and E12.5, as well as the renal vesicles at P4. We define an early stage of progenitor cell induction driven primarily by gene repression. Surprising stochastic expression of marker genes associated with differentiated cell types was observed in E11.5 progenitors. We provide a global view of the polarized gene expression already present in the renal vesicle, the first epithelial precursor of the nephron. We show that Hox gene readthrough transcripts can be spliced to produce intergenic homeobox swaps. We also identify a surprising number of genes with partially degraded noncoding RNA. Perhaps most interesting, at early developmental times single cells often expressed genes related to several developmental pathways. This provides powerful evidence that initial organogenesis involves a process of multilineage priming. This is followed by a combination of gene repression, which turns off the genes associated with most possible lineages, and the activation of increasing numbers of genes driving the chosen developmental direction.

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.

A blastocyst will implant only when the uterus becomes receptive. Following attachment, luminal epithelial cells undergo degeneration at the site of the blastocyst. Although many genes critical for uterine receptivity are primarily regulated by ovarian hormones, Kruppel-like factor 5 (KLF5), a zinc finger-containing transcription factor, is persistently expressed in epithelial cells independently of ovarian hormones. Loss of uterine Klf5 causes female infertility due to defective implantation. Cox2 is normally expressed in the luminal epithelium and stroma at the site of blastocyst attachment, but luminal epithelial COX2 expression is absent with loss of Klf5. This is associated with the retention of the epithelium around the implantation chamber with arrested embryonic growth. These results suggest that Klf5 is indispensable for normal implantation.decidualization | fertility | pregnancy