Single-Cell Transcriptomic Profiling Identifies Molecular Phenotypes of Newborn Human Lung Cells

While animal model studies have extensively defined the mechanisms controlling cell diversity in the developing mammalian lung, there exists a significant knowledge gap with regards to late-stage human lung development. The NHLBI Molecular Atlas of Lung Development Program (LungMAP) seeks to fill this gap by creating a structural, cellular and molecular atlas of the human and mouse lung. Transcriptomic profiling at the single-cell level created a cellular atlas of newborn human lungs. Frozen single-cell isolates obtained from two newborn human lungs from the LungMAP Human Tissue Core Biorepository, were captured, and library preparation was completed on the Chromium 10X system. Data was analyzed in Seurat, and cellular annotation was performed using the ToppGene functional analysis tool. Transcriptional interrogation of 5500 newborn human lung cells identified distinct clusters representing multiple populations of epithelial, endothelial, fibroblasts, pericytes, smooth muscle, immune cells and their gene signatures. Computational integration of data from newborn human cells and with 32,000 cells from postnatal days 1 through 10 mouse lungs generated by the LungMAP Cincinnati Research Center facilitated the identification of distinct cellular lineages among all the major cell types. Integration of the newborn human and mouse cellular transcriptomes also demonstrated cell type-specific differences in maturation states of newborn human lung cells. Specifically, newborn human lung matrix fibroblasts could be separated into those representative of younger cells (n = 393), or older cells (n = 158). Cells with each molecular profile were spatially resolved within newborn human lung tissue. This is the first comprehensive molecular map of the cellular landscape of neonatal human lung, including biomarkers for cells at distinct states of maturity.

Single-cell sequencing dissects the transcriptional identity of activated fibroblasts and identifies novel persistent distal tubular injury patterns in kidney fibrosis

Examining kidney fibrosis is crucial for mechanistic understanding and developing targeted strategies against chronic kidney disease (CKD). Persistent fibroblast activation and tubular epithelial cell (TEC) injury are key CKD contributors. However, cellular and transcriptional landscapes of CKD and specific activated kidney fibroblast clusters remain elusive. Here, we analyzed single cell transcriptomic profiles of two clinically relevant kidney fibrosis models which induced robust kidney parenchymal remodeling. We dissected the molecular and cellular landscapes of kidney stroma and newly identified three distinctive fibroblast clusters with "secretory", "contractile" and "vascular" transcriptional enrichments. Also, both injuries generated failed repair TECs (frTECs) characterized by decline of mature epithelial markers and elevation of stromal and injury markers. Notably, frTECs shared transcriptional identity with distal nephron segments of the embryonic kidney. Moreover, we identified that both models exhibited robust and previously unrecognized distal spatial pattern of TEC injury, outlined by persistent elevation of renal TEC injury markers including Krt8 and Vcam1, while the surviving proximal tubules (PTs) showed restored transcriptional signature. We also found that long-term kidney injuries activated a prominent nephrogenic signature, including Sox4 and Hox gene elevation, which prevailed in the distal tubular segments. Our findings might advance understanding of and targeted intervention in fibrotic kidney disease.

Single-cell sequencing dissects the transcriptional identity of activated fibroblasts and identifies novel persistent distal tubular injury patterns in kidney fibrosis

Examining kidney fibrosis is crucial for mechanistic understanding and developing targeted strategies against chronic kidney disease (CKD). Persistent fibroblast activation and tubular epithelial cell (TEC) injury are key CKD contributors. However, cellular and transcriptional landscapes of CKD and specific activated kidney fibroblast clusters remain elusive. Here, we analyzed single cell transcriptomic profiles of two clinically relevant kidney fibrosis models which induced robust kidney parenchymal remodeling. We dissected the molecular and cellular landscapes of kidney stroma and newly identified three distinctive fibroblast clusters with "secretory", "contractile" and "vascular" transcriptional enrichments. Also, both injuries generated failed repair TECs (frTECs) characterized by decline of mature epithelial markers and elevation of stromal and injury markers. Notably, frTECs shared transcriptional identity with distal nephron segments of the embryonic kidney. Moreover, we identified that both models exhibited robust and previously unrecognized distal spatial pattern of TEC injury, outlined by persistent elevation of renal TEC injury markers including Krt8, while the surviving proximal tubules (PTs) showed restored transcriptional signature. Furthermore, we found that long-term kidney injuries activated a prominent nephrogenic signature, including Sox4 and Hox gene elevation, which prevailed in the distal tubular segments. Our findings might advance understanding of and targeted intervention in fibrotic kidney disease.

Single cell transcriptomic analysis of HPV16-infected epithelium identifies a keratinocyte subpopulation implicated in cancer

Persistent HPV16 infection is a major cause of the global cancer burden. The viral life cycle is dependent on the differentiation program of stratified squamous epithelium, but the landscape of keratinocyte subpopulations which support distinct phases of the viral life cycle has yet to be elucidated. Here, single cell RNA sequencing of HPV16 infected compared to uninfected organoids identifies twelve distinct keratinocyte populations, with a subset mapped to reconstruct their respective 3D geography in stratified squamous epithelium. Instead of conventional terminally differentiated cells, an HPV-reprogrammed keratinocyte subpopulation (HIDDEN cells) forms the surface compartment and requires overexpression of the ELF3/ESE-1 transcription factor. HIDDEN cells are detected throughout stages of human carcinogenesis including primary human cervical intraepithelial neoplasias and HPV positive head and neck cancers, and a possible role in promoting viral carcinogenesis is supported by TCGA analyses. Single cell transcriptome information on HPV-infected versus uninfected epithelium will enable broader studies of the role of individual keratinocyte subpopulations in tumor virus infection and cancer evolution.

IL-1/MyD88-Dependent G-CSF and IL-6 Secretion Mediates Postburn Anemia

The anemia of critical illness (ACI) is a nearly universal pathophysiological consequence of burn injury and a primary reason burn patients require massive quantities of transfused blood. Inflammatory processes are expected to drive postburn ACI and prevent meaningful erythropoietic stimulation through iron or erythropoietin supplementation, but to this day no specific inflammatory pathways have been identified as a critical mechanism. In this study, we examined whether secretion of G-CSF and IL-6 mediates distinct features of postburn ACI and interrogated inflammatory mechanisms that could be responsible for their secretion. Our analysis of mouse and human skin samples identified the burn wound as a primary source of G-CSF and IL-6 secretion. We show that G-CSF and IL-6 are secreted independently through an IL-1/MyD88-dependent mechanism, and we ruled out TLR2 and TLR4 as critical receptors. Our results indicate that IL-1/MyD88-dependent G-CSF secretion plays a key role in impairing medullary erythropoiesis and IL-6 secretion plays a key role in limiting the access of erythroid cells to iron. Importantly, we found that IL-1α/β neutralizing Abs broadly attenuated features of postburn ACI that could be attributed to G-CSF or IL-6 secretion and rescued deficits of circulating RBC counts, hemoglobin, and hematocrit caused by burn injury. We conclude that wound-based IL-1/MyD88 signaling mediates postburn ACI through induction of G-CSF and IL-6 secretion.

Fibroblast inflammatory priming determines regenerative versus fibrotic skin repair in reindeer

Adult mammalian skin wounds heal by forming fibrotic scars. We report that full-thickness injuries of reindeer antler skin (velvet) regenerate, whereas back skin forms fibrotic scar. Single-cell multi-omics reveal that uninjured velvet fibroblasts resemble human fetal fibroblasts, whereas back skin fibroblasts express inflammatory mediators mimicking pro-fibrotic adult human and rodent fibroblasts. Consequently, injury elicits site-specific immune responses: back skin fibroblasts amplify myeloid infiltration and maturation during repair, whereas velvet fibroblasts adopt an immunosuppressive phenotype that restricts leukocyte recruitment and hastens immune resolution. Ectopic transplantation of velvet to scar-forming back skin is initially regenerative, but progressively transitions to a fibrotic phenotype akin to the scarless fetal-to-scar-forming transition reported in humans. Skin regeneration is diminished by intensifying, or enhanced by neutralizing, these pathologic fibroblast-immune interactions. Reindeer represent a powerful comparative model for interrogating divergent wound healing outcomes, and our results nominate decoupling of fibroblast-immune interactions as a promising approach to mitigate scar.

Cross talk between LAM cells and fibroblasts may influence alveolar epithelial cell behavior in lymphangioleiomyomatosis

Lymphangioleiomyomatosis (LAM) is a female-specific cystic lung disease in which tuberous sclerosis complex 2 (TSC2)-deficient LAM cells, LAM-associated fibroblasts (LAFs), and other cell types infiltrate the lungs. LAM lesions can be associated with type II alveolar epithelial (AT2) cells. We hypothesized that the behavior of AT2 cells in LAM is influenced locally by LAFs. We tested this hypothesis in the patient samples and in vitro. In human LAM lung, nodular AT2 cells show enhanced proliferation when compared with parenchymal AT2 cells, demonstrated by increased Ki67 expression. Furthermore, nodular AT2 cells express proteins associated with epithelial activation in other disease states including matrix metalloproteinase 7, and fibroblast growth factor 7 (FGF7). In vitro, LAF-conditioned medium is mitogenic and positively chemotactic for epithelial cells, increases the rate of epithelial repair, and protects against apoptosis. In vitro, LAM patient-derived TSC2 null cells cocultured with LAFs upregulate LAF expression of the epithelial chemokine and mitogen FGF7, a potential mediator of fibroblast-epithelial cross talk, in a mechanistic target of rapamycin (mTOR)-dependent manner. In a novel in vitro model of LAM, ex vivo cultured LAM lung-derived microtissues promote both epithelial migration and adhesion. Our findings suggest that AT2 cells in LAM display a proliferative, activated phenotype and fibroblast accumulation following LAM cell infiltration into the parenchyma contributes to this change in AT2 cell behavior. Fibroblast-derived FGF7 may contribute to the cross talk between LAFs and hyperplastic epithelium in vivo, but does not appear to be the main driver of the effects of LAFs on epithelial cells in vitro.

HoxD transcription factors define monosynaptic sensory-motor specificity in the developing spinal cord

The specificity of monosynaptic connections between proprioceptive sensory neurons and their recipient spinal motor neurons depends on multiple factors, including motor neuron positioning and dendrite morphology, axon projection patterns of proprioceptive sensory neurons in the spinal cord, and the ligand-receptor molecules involved in cell-to-cell recognition. However, with few exceptions, the transcription factors engaged in this process are poorly characterized. Here, we show that members of the HoxD family of transcription factors play a crucial role in the specificity of monosynaptic sensory-motor connections. Mice lacking Hoxd9, Hoxd10 and Hoxd11 exhibit defects in locomotion but have no obvious defects in motor neuron positioning or dendrite morphology through the medio-lateral and rostro-caudal axes. However, we found that quadriceps motor neurons in these mice show aberrant axon development and receive inappropriate inputs from proprioceptive sensory axons innervating the obturator muscle. These genetic studies demonstrate that the HoxD transcription factors play an integral role in the synaptic specificity of monosynaptic sensory-motor connections in the developing spinal cord.

The Rhesus Macaque Serves As a Model for Human Lateral Branch Nephrogenesis

BACKGROUND

Most nephrons are added in late gestation. Truncated extrauterine nephrogenesis in premature infants results in fewer nephrons and significantly increased risk for CKD in adulthood. To overcome the ethical and technical difficulties associated with studies of late-gestation human fetal kidney development, third-trimester rhesus macaques served as a model to understand lateral branch nephrogenesis (LBN) at the molecular level.

METHODS

Immunostaining and 3D rendering assessed morphology. Single-cell (sc) and single-nucleus (sn) RNA-Seq were performed on four cortically enriched fetal rhesus kidneys of 129-131 days gestational age (GA). An integrative bioinformatics strategy was applied across single-cell modalities, species, and time. RNAScope validation studies were performed on human archival tissue.

RESULTS

Third-trimester rhesus kidney undergoes human-like LBN. scRNA-Seq of 23,608 cells revealed 37 transcriptionally distinct cell populations, including naïve nephron progenitor cells (NPCs), with the prior noted marker genes CITED1, MEOX1, and EYA1 (c25). These same populations and markers were reflected in snRNA-Seq of 5972 nuclei. Late-gestation rhesus NPC markers resembled late-gestation murine NPC, whereas early second-trimester human NPC markers aligned to midgestation murine NPCs. New, age-specific rhesus NPCs (SHISA8) and ureteric buds (POU3F4 and TWIST) predicted markers were verified in late-gestation human archival samples.

CONCLUSIONS

Rhesus macaque is the first model of bona fide LBN, enabling molecular studies of late gestation, human-like nephrogenesis. These molecular findings support the hypothesis that aging nephron progenitors have a distinct molecular signature and align to their earlier human counterparts, with unique markers highlighting LBN-specific progenitor maturation.

Single-cell RNA-seq analysis reveals compartment-specific heterogeneity and plasticity of microglia

Microglia are ubiquitous central nervous system (CNS)-resident macrophages that maintain homeostasis of neural tissues and protect them from pathogen attacks. Yet, their differentiation in different compartments remains elusive. We performed single-cell RNA-seq to compare microglial subtypes in the cortex and the spinal cord. A multi-way comparative analysis was carried out on samples from C57/BL and HIV gp120 transgenic mice at two, four, and eight months of age. The results revealed overlapping but distinct microglial populations in the cortex and the spinal cord. The differential heterogeneity of microglia in these CNS regions was further suggested by their disparity of plasticity in response to life span progression and HIV-1 pathogenic protein gp120. Our findings indicate that microglia in different CNS compartments are adapted to their local environments to fulfill region-specific biological functions.

Kidney Single-cell Transcriptomes Predict Spatial Corticomedullary Gene Expression and Tissue Osmolality Gradients

BACKGROUND

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.

METHODS

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. In situ 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^CD-/-), which display reduced renal medullary osmolality.

RESULTS

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^CD-/- mice indicated a flattened expression gradient of osmo-responsive genes compared with control mice, consistent with their physiologic phenotype.

CONCLUSIONS

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.

Combinatorial Transcriptional Profiling of Mouse and Human Enteric Neurons Identifies Shared and Disparate Subtypes In Situ

BACKGROUND & AIMS

The enteric nervous system (ENS) coordinates essential intestinal functions through the concerted action of diverse enteric neurons (ENs). However, integrated molecular knowledge of EN subtypes is lacking. To compare human and mouse ENs, we transcriptionally profiled healthy ENS from adult humans and mice. We aimed to identify transcripts marking discrete neuron subtypes and visualize conserved EN subtypes for humans and mice in multiple bowel regions.

METHODS

Human myenteric ganglia and adjacent smooth muscle were isolated by laser-capture microdissection for RNA sequencing. Ganglia-specific transcriptional profiles were identified by computationally subtracting muscle gene signatures. Nuclei from mouse myenteric neurons were isolated and subjected to single-nucleus RNA sequencing, totaling more than 4 billion reads and 25,208 neurons. Neuronal subtypes were defined using mouse single-nucleus RNA sequencing data. Comparative informatics between human and mouse data sets identified shared EN subtype markers, which were visualized in situ using hybridization chain reaction.

RESULTS

Several EN subtypes in the duodenum, ileum, and colon are conserved between humans and mice based on orthologous gene expression. However, some EN subtype-specific genes from mice are expressed in completely distinct morphologically defined subtypes in humans. In mice, we identified several neuronal subtypes that stably express gene modules across all intestinal segments, with graded, regional expression of 1 or more marker genes.

CONCLUSIONS

Our combined transcriptional profiling of human myenteric ganglia and mouse EN provides a rich foundation for developing novel intestinal therapeutics. There is congruency among some EN subtypes, but we note multiple species differences that should be carefully considered when relating findings from mouse ENS research to human gastrointestinal studies.

Single-Cell Transcriptomic Analysis Identifies a Unique Pulmonary Lymphangioleiomyomatosis Cell

Rationale: Lymphangioleiomyomatosis (LAM) is a metastatic neoplasm of reproductive-age women associated with mutations in tuberous sclerosis complex genes. LAM causes cystic remodeling of the lung and progressive respiratory failure. The sources and cellular characteristics of LAM cells underlying disease pathogenesis remain elusive.Objectives: Identification and characterization of LAM cells in human lung and uterus using a single-cell approach.Methods: Single-cell and single-nuclei RNA sequencing on LAM (n = 4) and control (n = 7) lungs, immunofluorescence confocal microscopy, ELISA, and aptamer proteomics were used to identify and validate LAMCORE cells and secreted biomarkers, predict cellular origins, and define molecular and cellular networks in LAM.Measurements and Main Results: A unique cell type termed LAMCORE was identified, which was distinct from, but closely related to, lung mesenchymal cells. LAMCORE cells expressing signature genes included known LAM markers such as PMEL, FIGF, CTSK, and MLANA and novel biomarkers validated by aptamer screening, ELISA, and immunofluorescence microscopy. LAM cells in lung and uterus are morphologically indistinguishable and share similar gene expression profiles and biallelic TSC2 mutations, supporting a potential uterine origin for the LAMCORE cell. Effects of LAM on resident pulmonary cell types indicated recruitment and activation of lymphatic endothelial cells.Conclusions: A unique population of LAMCORE cells was identified in lung and uterus of patients with LAM, sharing close transcriptomic identity. LAM cell selective markers, secreted biomarkers, and the predicted cellular molecular features provide new insights into the signaling and transcriptional programs that may serve as diagnostic markers and therapeutic targets to influence the pathogenesis of LAM.

Single-Cell Profiling of AKI in a Murine Model Reveals Novel Transcriptional Signatures, Profibrotic Phenotype, and Epithelial-to-Stromal Crosstalk

BACKGROUND

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.

METHODS

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 in situ hybridizations, we validated AKI signatures in multiple experiments.

RESULTS

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 (Ahnak, Sh3bgrl3, and Col18a1) not previously examined in kidney pathologies. AKI induced proximal tubule dedifferentiation, with a pronounced nephrogenic signature represented by Sox4 and Cd24a. 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.

CONCLUSIONS

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.

Integrated Transcriptome and Network Analysis Reveals Spatiotemporal Dynamics of Calvarial Suturogenesis

Craniofacial abnormalities often involve sutures, the growth centers of the skull. To characterize the organization and processes governing their development, we profile the murine frontal suture, a model for sutural growth and fusion, at the tissue- and single-cell level on embryonic days (E)16.5 and E18.5. For the wild-type suture, bulk RNA sequencing (RNA-seq) analysis identifies mesenchyme-, osteogenic front-, and stage-enriched genes and biological processes, as well as alternative splicing events modifying the extracellular matrix. Single-cell RNA-seq analysis distinguishes multiple subpopulations, of which five define a mesenchyme-osteoblast differentiation trajectory and show variation along the anteroposterior axis. Similar analyses of in vivo mouse models of impaired frontal suturogenesis in Saethre-Chotzen and Apert syndromes, Twist1+/- and Fgfr2+/S252W, demonstrate distinct transcriptional changes involving angiogenesis and ribogenesis, respectively. Co-expression network analysis reveals gene expression modules from which we validate key driver genes regulating osteoblast differentiation. Our study provides a global approach to gain insights into suturogenesis.

Physical interactions between Gsx2 and Ascl1 balance progenitor expansion versus neurogenesis in the mouse lateral ganglionic eminence

The Gsx2 homeodomain transcription factor promotes neural progenitor identity in the lateral ganglionic eminence (LGE), despite upregulating the neurogenic factor Ascl1. How this balance in maturation is maintained is unclear. Here, we show that Gsx2 and Ascl1 are co-expressed in subapical progenitors that have unique transcriptional signatures in LGE ventricular zone (VZ) cells. Moreover, whereas Ascl1 misexpression promotes neurogenesis in dorsal telencephalic progenitors, the co-expression of Gsx2 with Ascl1 inhibits neurogenesis. Using luciferase assays, we found that Gsx2 reduces the ability of Ascl1 to activate gene expression in a dose-dependent and DNA binding-independent manner. Furthermore, Gsx2 physically interacts with the basic helix-loop-helix (bHLH) domain of Ascl1, and DNA-binding assays demonstrated that this interaction interferes with the ability of Ascl1 to bind DNA. Finally, we modified a proximity ligation assay for tissue sections and found that Ascl1-Gsx2 interactions are enriched within LGE VZ progenitors, whereas Ascl1-Tcf3 (E-protein) interactions predominate in the subventricular zone. Thus, Gsx2 contributes to the balance between progenitor maintenance and neurogenesis by physically interacting with Ascl1, interfering with its DNA binding and limiting neurogenesis within LGE progenitors.

Single-Cell Transcriptomics in Medulloblastoma Reveals Tumor-Initiating Progenitors and Oncogenic Cascades during Tumorigenesis and Relapse

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.

Single-cell RNA sequencing for the study of development, physiology and disease

An ongoing technological revolution is continually improving our ability to carry out very high-resolution studies of gene expression patterns. Current technology enables the global gene expression profiles of single cells to be defined, facilitating dissection of heterogeneity in cell populations that was previously hidden. In contrast to gene expression studies that use bulk RNA samples and provide only a virtual average of the diverse constituent cells, single-cell studies enable the molecular distinction of all cell types within a complex population mix, such as a tumour or developing organ. For instance, single-cell gene expression profiling has contributed to improved understanding of how histologically identical, adjacent cells make different differentiation decisions during development. Beyond development, single-cell gene expression studies have enabled the characteristics of previously known cell types to be more fully defined and facilitated the identification of novel categories of cells, contributing to improvements in our understanding of both normal and disease-related physiological processes and leading to the identification of new treatment approaches. Although limitations remain to be overcome, technology for the analysis of single-cell gene expression patterns is improving rapidly and beginning to provide a detailed atlas of the gene expression patterns of all cell types in the human body.

The Pediatric Cell Atlas: Defining the Growth Phase of Human Development at Single-Cell Resolution

Single-cell gene expression analyses of mammalian tissues have uncovered profound stage-specific molecular regulatory phenomena that have changed the understanding of unique cell types and signaling pathways critical for lineage determination, morphogenesis, and growth. We discuss here the case for a Pediatric Cell Atlas as part of the Human Cell Atlas consortium to provide single-cell profiles and spatial characterization of gene expression across human tissues and organs. Such data will complement adult and developmentally focused HCA projects to provide a rich cytogenomic framework for understanding not only pediatric health and disease but also environmental and genetic impacts across the human lifespan.

Integrating gene regulatory pathways into differential network analysis of gene expression data

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.

Single cell RNA-seq study of wild type and Hox9,10,11 mutant developing uterus

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.

Maturation of heart valve cell populations during postnatal remodeling

Heart valve cells mediate extracellular matrix (ECM) remodeling during postnatal valve leaflet stratification, but phenotypic and transcriptional diversity of valve cells in development is largely unknown. Single cell analysis of mouse heart valve cells was used to evaluate cell heterogeneity during postnatal ECM remodeling and leaflet morphogenesis. The transcriptomic analysis of single cells from postnatal day (P)7 and P30 murine aortic (AoV) and mitral (MV) heart valves uncovered distinct subsets of melanocytes, immune and endothelial cells present at P7 and P30. By contrast, interstitial cell populations are different from P7 to P30. P7 valve leaflets exhibit two distinct collagen- and glycosaminoglycan-expressing interstitial cell clusters, and prevalent ECM gene expression. At P30, four interstitial cell clusters are apparent with leaflet specificity and differential expression of complement factors, ECM proteins and osteogenic genes. This initial transcriptomic analysis of postnatal heart valves at single cell resolution demonstrates that subpopulations of endothelial and immune cells are relatively constant throughout postnatal development, but interstitial cell subpopulations undergo changes in gene expression and cellular functions in primordial and mature valves.

Hnf4a deletion in the mouse kidney phenocopies Fanconi renotubular syndrome

Different nephron tubule segments perform distinct physiological functions, collectively acting as a blood filtration unit. Dysfunction of the proximal tubule segment can lead to Fanconi renotubular syndrome (FRTS), with major symptoms such as excess excretion of water, glucose, and phosphate in the urine. It has been shown that a mutation in HNF4A is associated with FRTS in humans and that Hnf4a is expressed specifically in proximal tubules in adult rat nephrons. However, little is known about the role of Hnf4a in nephrogenesis. Here, we found that Hnf4a is expressed in both presumptive and differentiated proximal tubules in the developing mouse kidney. We show that Hnf4a is required for the formation of differentiated proximal tubules but is dispensable for the formation of presumptive proximal tubules. Furthermore, we show that loss of Hnf4a decreased the expression of proximal tubule-specific genes. Adult Hnf4a mutant mice presented with FRTS-like symptoms, including polyuria, polydipsia, glycosuria, and phosphaturia. Analysis of the adult Hnf4a mutant kidney also showed proximal tubule dysgenesis and nephrocalcinosis. Our results demonstrate the critical role of Hnf4a in proximal tubule development and provide mechanistic insight into the etiology of FRTS.

Disruption of Hox9,10,11 function results in cellular level lineage infidelity in the kidney

Hox genes are important regulators of development. The 39 mammalian Hox genes have considerable functional overlap, greatly confounding their study. In this report, we generated mice with multiple combinations of paralogous and flanking Abd-B Hox gene mutations to investigate functional redundancies in kidney development. The resulting mice developed a number of kidney abnormalities, including hypoplasia, agenesis, and severe cysts, with distinct Hox functions observed in early metanephric kidney formation and nephron progenitor maintenance. Most surprising, however, was that extensive removal of Hox shared function in these kidneys resulted in cellular level lineage infidelity. Strikingly, mutant nephron tubules consisted of intermixed cells with proximal tubule, loop of Henle, and collecting duct identities, with some single cells expressing markers associated with more than one nephron segment. These results indicate that Hox genes are required for proper lineage selection/maintenance and full repression of genes involved in cell fate restriction in the developing kidney.

Reduced Abd-B Hox function during kidney development results in lineage infidelity

Hox genes can function as key drivers of segment identity, with Hox mutations in Drosophila often resulting in dramatic homeotic transformations. In addition, however, they can serve other essential functions. In mammals, the study of Hox gene roles in development is complicated by the presence of four Hox clusters with a total of 39 genes showing extensive functional overlap. In this study, in order to better understand shared core Hox functions, we examined kidney development in mice with frameshift mutations of multiple Abd-B type Hox genes. The resulting phenotypes included dramatically reduced branching morphogenesis of the ureteric bud, premature depletion of nephron progenitors and abnormal development of the stromal compartment. Most unexpected, however, we also observed a cellular level lineage infidelity in nephron segments. Scattered cells within the proximal tubules, for example, expressed genes normally expressed only in collecting ducts. Multiple combinations of inappropriate nephron segment specific marker expression were found. In some cases, cells within a tubule showed incorrect identity, while in other cases cells showed ambiguous character, with simultaneous expression of genes associated with more than one nephron segment. These results give evidence that Hox genes have an overlapping core function at the cellular level in driving and/or maintaining correct differentiation decisions.

Cross-platform single cell analysis of kidney development shows stromal cells express Gdnf

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.

Activated Hedgehog-GLI Signaling Causes Congenital Ureteropelvic Junction Obstruction

Intrinsic ureteropelvic junction obstruction is the most common cause of congenital hydronephrosis, yet the underlying pathogenesis is undefined. Hedgehog proteins control morphogenesis by promoting GLI-dependent transcriptional activation and inhibiting the formation of the GLI3 transcriptional repressor. Hedgehog regulates differentiation and proliferation of ureteric smooth muscle progenitor cells during murine kidney-ureter development. Histopathologic findings of smooth muscle cell hypertrophy and stroma-like cells, consistently observed in obstructing tissue at the time of surgical correction, suggest that Hedgehog signaling is abnormally regulated during the genesis of congenital intrinsic ureteropelvic junction obstruction. Here, we demonstrate that constitutively active Hedgehog signaling in murine intermediate mesoderm-derived renal progenitors results in hydronephrosis and failure to develop a patent pelvic-ureteric junction. Tissue obstructing the ureteropelvic junction was marked as early as E13.5 by an ectopic population of cells expressing Ptch2, a Hedgehog signaling target. Constitutive expression of GLI3 repressor in Ptch1-deficient mice rescued ectopic Ptch2 expression and obstructive hydronephrosis. Whole transcriptome analysis of isolated Ptch2⁺ cells revealed coexpression of genes characteristic of stromal progenitor cells. Genetic lineage tracing indicated that stromal cells blocking the ureteropelvic junction were derived from intermediate mesoderm-derived renal progenitors and were distinct from the smooth muscle or epithelial lineages. Analysis of obstructive ureteric tissue resected from children with congenital intrinsic ureteropelvic junction obstruction revealed a molecular signature similar to that observed in Ptch1-deficient mice. Together, these results demonstrate a Hedgehog-dependent mechanism underlying mammalian intrinsic ureteropelvic junction obstruction.

Repression of Interstitial Identity in Nephron Progenitor Cells by Pax2 Establishes the Nephron-Interstitium Boundary during Kidney Development

The kidney contains the functional units, the nephrons, surrounded by the renal interstitium. Previously we discovered that, once Six2-expressing nephron progenitor cells and Foxd1-expressing renal interstitial progenitor cells form at the onset of kidney development, descendant cells from these populations contribute exclusively to the main body of nephrons and renal interstitial tissues, respectively, indicating a lineage boundary between the nephron and renal interstitial compartments. Currently it is unclear how lineages are regulated during kidney organogenesis. We demonstrate that nephron progenitor cells lacking Pax2 fail to differentiate into nephron cells but can switch fates into renal interstitium-like cell types. These data suggest that Pax2 function maintains nephron progenitor cells by repressing a renal interstitial cell program. Thus, the lineage boundary between the nephron and renal interstitial compartments is maintained by the Pax2 activity in nephron progenitor cells during kidney organogenesis.

Psychrophilic proteases dramatically reduce single-cell RNA-seq artifacts: a molecular atlas of kidney development

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 in vivo 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.

Lung Gene Expression Analysis (LGEA): an integrative web portal for comprehensive gene expression data analysis in lung development

'LungGENS', our previously developed web tool for mapping single-cell gene expression in the developing lung, has been well received by the pulmonary research community. With continued support from the 'LungMAP' consortium, we extended the scope of the LungGENS database to accommodate transcriptomics data from pulmonary tissues and cells from human and mouse at different stages of lung development. Lung Gene Expression Analysis (LGEA) web portal is an extended version of LungGENS useful for the analysis, display and interpretation of gene expression patterns obtained from single cells, sorted cell populations and whole lung tissues. The LGEA web portal is freely available at http://research.cchmc.org/pbge/lunggens/mainportal.html.

Intrinsic Age-Dependent Changes and Cell-Cell Contacts Regulate Nephron Progenitor Lifespan

During fetal development, nephrons of the metanephric kidney form from a mesenchymal progenitor population that differentiates en masse before or shortly after birth. We explored intrinsic and extrinsic mechanisms controlling progenitor lifespan in a transplantation assay that allowed us to compare engraftment of old and young progenitors into the same young niche. The progenitors displayed an age-dependent decrease in proliferation and concomitant increase in niche exit rates. Single-cell transcriptome profiling revealed progressive age-dependent changes, with heterogeneity increasing in older populations. Age-dependent elevation in mTor and reduction in Fgf20 could contribute to increased exit rates. Importantly, 30% of old progenitors remained in the niche for up to 1 week post engraftment, a net gain of 50% to their lifespan, but only if surrounded by young neighbors. We provide evidence in support of a model in which intrinsic age-dependent changes affect inter-progenitor interactions that drive cessation of nephrogenesis.

Uterine inactivation of muscle segment homeobox (Msx) genes alters epithelial cell junction proteins during embryo implantation

Embryo implantation requires that the uterus differentiate into the receptive state. Failure to attain uterine receptivity will impede blastocyst attachment and result in a compromised pregnancy. The molecular mechanism by which the uterus transitions from the prereceptive to the receptive stage is complex, involving an intricate interplay of various molecules. We recently found that mice with uterine deletion ofMsxgenes (Msx1(d/d)/Msx2(d/d)) are infertile because of implantation failure associated with heightened apicobasal polarity of luminal epithelial cells during the receptive period. However, information on Msx's roles in regulating epithelial polarity remains limited. To gain further insight, we analyzed cell-type-specific gene expression by RNA sequencing of separated luminal epithelial and stromal cells by laser capture microdissection fromMsx1(d/d)/Msx2(d/d)and floxed mouse uteri on d 4 of pseudopregnancy. We found that claudin-1, a tight junction protein, and small proline-rich (Sprr2) protein, a major component of cornified envelopes in keratinized epidermis, were substantially up-regulated inMsx1(d/d)/Msx2(d/d)uterine epithelia. These factors also exhibited unique epithelial expression patterns at the implantation chamber (crypt) inMsx1(f/f)/Msx2(f/f)females; the patterns were lost inMsx1(d/d)/Msx2(d/d)epithelia on d 5, suggesting important roles during implantation. The results suggest thatMsxgenes play important roles during uterine receptivity including modulation of epithelial junctional activity.-Sun, X., Park, C. B., Deng, W., Potter, S. S., Dey, S. K. Uterine inactivation of muscle segment homeobox (Msx) genes alters epithelial cell junction proteins during embryo implantation.

Defects in the Fanconi Anemia Pathway in Head and Neck Cancer Cells Stimulate Tumor Cell Invasion through DNA-PK and Rac1 Signaling

PURPOSE

Head and neck squamous cell carcinoma (HNSCC) remains a devastating disease, and Fanconi anemia (FA) gene mutations and transcriptional repression are common. Invasive tumor behavior is associated with poor outcome, but relevant pathways triggering invasion are poorly understood. There is a significant need to improve our understanding of genetic pathways and molecular mechanisms driving advanced tumor phenotypes, to develop tailored therapies. Here we sought to investigate the phenotypic and molecular consequences of FA pathway loss in HNSCC cells.

EXPERIMENTAL DESIGN

Using sporadic HNSCC cell lines with and without FA gene knockdown, we sought to characterize the phenotypic and molecular consequences of FA deficiency. FA pathway inactivation was confirmed by the detection of classic hallmarks of FA following exposure to DNA cross-linkers. Cells were subjected to RNA sequencing with qRT-PCR validation, followed by cellular adhesion and invasion assays in the presence and absence of DNA-dependent protein kinase (DNA-PK) and Rac1 inhibitors.

RESULTS

We demonstrate that FA loss in HNSCC cells leads to cytoskeletal reorganization and invasive tumor cell behavior in the absence of proliferative gains. We further demonstrate that cellular invasion following FA loss is mediated, at least in part, through NHEJ-associated DNA-PK and downstream Rac1 GTPase activity.

CONCLUSIONS

These findings demonstrate that FA loss stimulates HNSCC cell motility and invasion, and implicate a targetable DNA-PK/Rac1 signaling axis in advanced tumor phenotypes.

SINCERA: A Pipeline for Single-Cell RNA-Seq Profiling Analysis

A major challenge in developmental biology is to understand the genetic and cellular processes/programs driving organ formation and differentiation of the diverse cell types that comprise the embryo. While recent studies using single cell transcriptome analysis illustrate the power to measure and understand cellular heterogeneity in complex biological systems, processing large amounts of RNA-seq data from heterogeneous cell populations creates the need for readily accessible tools for the analysis of single-cell RNA-seq (scRNA-seq) profiles. The present study presents a generally applicable analytic pipeline (SINCERA: a computational pipeline for SINgle CEll RNA-seq profiling Analysis) for processing scRNA-seq data from a whole organ or sorted cells. The pipeline supports the analysis for: 1) the distinction and identification of major cell types; 2) the identification of cell type specific gene signatures; and 3) the determination of driving forces of given cell types. We applied this pipeline to the RNA-seq analysis of single cells isolated from embryonic mouse lung at E16.5. Through the pipeline analysis, we distinguished major cell types of fetal mouse lung, including epithelial, endothelial, smooth muscle, pericyte, and fibroblast-like cell types, and identified cell type specific gene signatures, bioprocesses, and key regulators. SINCERA is implemented in R, licensed under the GNU General Public License v3, and freely available from CCHMC PBGE website, https://research.cchmc.org/pbge/sincera.html.

Key pathways regulated by HoxA9,10,11/HoxD9,10,11 during limb development

Background

The 39 mammalian Hox genes show problematic patterns of functional overlap. In order to more fully define the developmental roles of Hox genes it is necessary to remove multiple combinations of paralogous and flanking genes. In addition, the downstream molecular pathways regulated by Hox genes during limb development remain incompletely delineated.

Results

In this report we examine limb development in mice with frameshift mutations in six Hox genes, Hoxa9,10,11 and Hoxd9,10,11. The mice were made with a novel recombineering method that allows the simultaneous targeting of frameshift mutations into multiple flanking genes. The Hoxa9,10,11^−/−/Hoxd9,10,11^−/− mutant mice show a reduced ulna and radius that is more severe than seen in Hoxa11^−/−/Hoxd11^−/− mice, indicating a minor role for the flanking Hox9,10 genes in zeugopod development, as well as their primary function in stylopod development. The mutant mice also show severe reduction of Shh expression in the zone of polarizing activity, and decreased Fgf8 expression in the apical ectodermal ridge, thereby better defining the roles of these specific Hox genes in the regulation of critical signaling centers during limb development. Importantly, we also used laser capture microdissection coupled with RNA-Seq to characterize the gene expression programs in wild type and mutant limbs. Resting, proliferative and hypertrophic compartments of E15.5 forelimb zeugopods were examined. The results provide an RNA-Seq characterization of the progression of gene expression patterns during normal endochondral bone formation. In addition of the Hox mutants showed strongly altered expression of Pknox2, Zfp467, Gdf5, Bmpr1b, Dkk3, Igf1, Hand2, Shox2, Runx3, Bmp7 and Lef1, all of which have been previously shown to play important roles in bone formation.

Conclusions

The recombineering based frameshift mutation of the six flanking and paralogous Hoxa9,10,11 and Hoxd9,10,11 genes provides a resource for the analysis of their overlapping functions. Analysis of the Hoxa9,10,11^−/−/Hoxd9,10,11^−/− mutant limbs confirms and extends the results of previous studies using mice with Hox mutations in single paralogous groups or with entire Hox cluster deletions. The RNA-Seq analysis of specific compartments of the normal and mutant limbs defines the multiple key perturbed pathways downstream of these Hox genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-015-0078-5) contains supplementary material, which is available to authorized users.

Pathogenic pathways are activated in each major cell type of the glomerulus in the Cd2ap mutant mouse model of focal segmental glomerulosclerosis

Background

Mutations in several genes expressed in podocytes, including Cd2ap, have been associated with focal segmental glomerulosclerosis in humans. Mutant mouse models provide an opportunity to better understand the molecular pathology that drives these diseases.

Methods

In this report we use a battery of transgenic-GFP mice to facilitate the purification of all three major cell types of the glomerulus from Cd2ap mutant mice. Both microarrays and RNA-seq were used to characterize the gene expression profiles of the podocytes, mesangial cells and endothelial cells, providing a global dual platform cross-validating dataset.

Results

The mesangial cells showed increased expression of profibrotic factors, including thrombospondin, Tgfb2 and Tgfb3, as well as the angiogenesis factor Vegf. They also showed upregulation of protective genes, including Aldh1a2, involved in retinoic acid synthesis and Decorin, a Tgfb antagonist. Of interest, the mesangial cells also showed significant expression of Wt1, which has generally been considered podocyte specific. The Cd2ap mutant podocytes showed upregulation of proteases as well as genes involved in muscle and vasculature development and showed a very strong gene expression signature indicating programmed cell death. Endothelial cells showed increased expression of the leukocyte adhesion associated factors Vcam1 and Sele, as well as Midkine (promoting angiogenesis), endothelin and many genes responsive to cytokines and interferons.

Conclusions

This study provides a comprehensive analysis of the changing properties of the three cell types of the glomerulus in Cd2ap mutants, identifying activated and repressed pathways and responsible genes, thereby delivering a deeper molecular understanding of this genetic disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12882-015-0063-z) contains supplementary material, which is available to authorized users.

Single cell dissection of early kidney development: multilineage priming

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 read-through 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.

Building an atlas of gene expression driving kidney development: pushing the limits of resolution

Changing gene expression patterns is the essential driver of developmental processes. Growth factors, micro-RNAs, long intergenic noncoding RNAs, and epigenetic marks, such as DNA methylation and histone modifications, all work by impacting gene expression. The key features of developing cells, including their ability to communicate with others, are defined primarily by their gene-expression profiles. It is therefore clear that a gene-expression atlas of the developing kidney can provide a useful tool for the developmental nephrology research community. Toward this end, the GenitoUrinary Development Molecular Anatomy Project (GUDMAP) consortium has worked to create an atlas of the changing gene-expression patterns that drive kidney development. In this article, the global gene-expression profiling strategies of GUDMAP are reviewed. The initial work used laser-capture microdissection to purify multiple compartments of the developing kidney, including cap mesenchyme, renal vesicle, S-shaped bodies, proximal tubules, and more, which were then gene-expression profiled using microarrays. Resolution of the atlas was then improved by using transgenic mice with specific cell types labeled with green fluorescent protein (GFP), allowing their purification and profiling. In addition, RNA-Seq replaced microarrays. Currently, the atlas is being pushed to the single-cell resolution using microfluidic approaches that allow high-throughput RNA-Seq analysis of hundreds of individual cells. Results can identify novel types of cells and define interesting heterogeneities present within cell populations.

SP8 regulates signaling centers during craniofacial development

Much of the bone, cartilage and smooth muscle of the vertebrate face is derived from neural crest (NC) cells. During craniofacial development, the anterior neural ridge (ANR) and olfactory pit (OP) signaling centers are responsible for driving the outgrowth, survival, and differentiation of NC populated facial prominences, primarily via FGF. While much is known about the functional importance of signaling centers, relatively little is understood of how these signaling centers are made and maintained. In this report we describe a dramatic craniofacial malformation in mice mutant for the zinc finger transcription factor gene Sp8. At E14.5 they show facial prominences that are reduced in size and underdeveloped, giving an almost faceless phenotype. At later times they show severe midline defects, excencephaly, hyperterlorism, cleft palate, and a striking loss of many NC and paraxial mesoderm derived cranial bones. Sp8 expression was primarily restricted to the ANR and OP regions during craniofacial development. Analysis of an extensive series of conditional Sp8 mutants confirmed the critical role of Sp8 in signaling centers, and not directly in the NC and paraxial mesoderm cells. The NC cells of the Sp8 mutants showed increased levels of apoptosis and decreased cell proliferation, thereby explaining the reduced sizes of the facial prominences. Perturbed gene expression in the Sp8 mutants was examined by laser capture microdissection coupled with microarrays, as well as in situ hybridization and immunostaining. The most dramatic differences included striking reductions in Fgf8 and Fgf17 expression in the ANR and OP signaling centers. We were also able to achieve genetic and pharmaceutical partial rescue of the Sp8 mutant phenotype by reducing Sonic Hedgehog (SHH) signaling. These results show that Sp8 primarily functions to promote Fgf expression in the ANR and OP signaling centers that drive the survival, proliferation, and differentiation of the NC and paraxial mesoderm that make the face.

Recombineering-based dissection of flanking and paralogous Hox gene functions in mouse reproductive tracts

Hox genes are key regulators of development. In mammals, the study of these genes is greatly confounded by their large number, overlapping functions and interspersed shared enhancers. Here, we describe the use of a novel recombineering strategy to introduce simultaneous frameshift mutations into the flanking Hoxa9, Hoxa10 and Hoxa11 genes, as well as their paralogs on the HoxD cluster. The resulting Hoxa9,10,11 mutant mice displayed dramatic synergistic homeotic transformations of the reproductive tracts, with the uterus anteriorized towards oviduct and the vas deferens anteriorized towards epididymis. The Hoxa9,10,11 mutant mice also provided a genetic setting that allowed the discovery of Hoxd9,10,11 redundant reproductive tract patterning function. Both shared and distinct Hox functions were defined. Hoxd9,10,11 play a crucial role in the regulation of uterine immune function. Non-coding non-polyadenylated RNAs were among the key Hox targets, with dramatic downregulation in mutants. We observed Hox cross-regulation of transcription and splicing. In addition, we observed a surprising anti-dogmatic apparent posteriorization of the uterine epithelium. In caudal regions of the uterus, the normal simple columnar epithelium flanking the lumen was replaced by a pseudostratified transitional epithelium, normally found near the more posterior cervix. These results identify novel molecular functions of Hox genes in the development of the male and female reproductive tracts.

Changes in the gene expression programs of renal mesangial cells during diabetic nephropathy

Background

Diabetic nephropathy is the leading cause of end stage renal disease. All three cell types of the glomerulus, podocytes, endothelial cells and mesangial cells, play important roles in diabetic nephropathy. In this report we used Meis1-GFP transgenic mice to purify mesangial cells from normal mice and from db/db mice, which suffer diabetic nephropathy. The purpose of the study is to better define the unique character of normal mesangial cells, and to characterize their pathogenic and protective responses during diabetic nephropathy.

Methods

Comprehensive gene expression states of the normal and diseased mesangial cells were defined with microarrays. By comparing the gene expression profiles of mesangial cells with those of multiple other renal cell types, including podocytes, endothelial cells and renal vesicles, it was possible to better define their exceptional nature, which includes smooth muscle, phagocytic and neuronal traits.

Results

The complete set of mesangial cell expressed transcription factors, growth factors and receptors were identified. In addition, the analysis of the mesangial cells from diabetic nephropathy mice characterized their changes in gene expression. Molecular functions and biological processes specific to diseased mesangial cells were characterized, identifying genes involved in extracellular matrix, cell division, vasculogenesis, and growth factor modulation. Selected gene changes considered of particular importance to the disease process were validated and localized within the glomuerulus by immunostaining. For example, thrombospondin, a key mediator of TGFβ signaling, was upregulated in the diabetic nephropathy mesangial cells, likely contributing to fibrosis. On the other hand the decorin gene was also upregulated, and expression of this gene has been strongly implicated in the reduction of TGFβ induced fibrosis.

Conclusions

The results provide an important complement to previous studies examining mesangial cells grown in culture. The remarkable qualities of the mesangial cell are more fully defined in both the normal and diabetic nephropathy diseased state. New gene expression changes and biological pathways are discovered, yielding a deeper understanding of the diabetic nephropathy pathogenic process, and identifying candidate targets for the development of novel therapies.

RNA-Seq defines novel genes, RNA processing patterns and enhancer maps for the early stages of nephrogenesis: Hox supergenes

During kidney development the cap mesenchyme progenitor cells both self renew and differentiate into nephrons. The balance between renewal and differentiation determines the final nephron count, which is of considerable medical importance. An important goal is to create a precise genetic definition of the early differentiation of cap mesenchyme progenitors. We used RNA-Seq to transcriptional profile the cap mesenchyme progenitors and their first epithelial derivative, the renal vesicles. The results provide a global view of the changing gene expression program during this key period, defining expression levels for all transcription factors, growth factors, and receptors. The RNA-Seq was performed using two different biochemistries, with one examining only polyadenylated RNA and the other total RNA. This allowed the analysis of noncanonical transcripts, which for many genes were more abundant than standard exonic RNAs. Since a large fraction of enhancers are now known to be transcribed the results also provide global maps of potential enhancers. Further, the RNA-Seq data defined hundreds of novel splice patterns and large numbers of new genes. Particularly striking was the extensive sense/antisense transcription and changing RNA processing complexities of the Hox clusters.

Laser capture

This chapter describes detailed methods used for laser capture microdissection (LCM) of discrete subpopulations of cells. Topics covered include preparing tissue blocks, cryostat sectioning, processing slides, performing the LCM, and purification of RNA from LCM samples. Notes describe the fine points of each operation, which can often mean the difference between success and failure.

Genes that confer the identity of the renin cell

Renin-expressing cells modulate BP, fluid-electrolyte homeostasis, and kidney development, but remarkably little is known regarding the genetic regulatory network that governs the identity of these cells. Here we compared the gene expression profiles of renin cells with most cells in the kidney at various stages of development as well as after a physiologic challenge known to induce the transformation of arteriolar smooth muscle cells into renin-expressing cells. At all stages, renin cells expressed a distinct set of genes characteristic of the renin phenotype, which was vastly different from other cell types in the kidney. For example, cells programmed to exhibit the renin phenotype expressed Akr1b7, and maturing cells expressed angiogenic factors necessary for the development of the kidney vasculature and RGS (regulator of G-protein signaling) genes, suggesting a potential relationship between renin cells and pericytes. Contrary to the plasticity of arteriolar smooth muscle cells upstream from the glomerulus, which can transiently acquire the embryonic phenotype in the adult under physiologic stress, the adult juxtaglomerular cell always possessed characteristics of both smooth muscle and renin cells. Taken together, these results identify the gene expression profile of renin-expressing cells at various stages of maturity, and suggest that juxtaglomerular cells maintain properties of both smooth muscle and renin-expressing cells, likely to allow the rapid control of body fluids and BP through both contractile and endocrine functions.

Defining the molecular character of the developing and adult kidney podocyte

Background

The podocyte is a remarkable cell type, which encases the capillaries of the kidney glomerulus. Although mesodermal in origin it sends out axonal like projections that wrap around the capillaries. These extend yet finer projections, the foot processes, which interdigitate, leaving between them the slit diaphragms, through which the glomerular filtrate must pass. The podocytes are a subject of keen interest because of their key roles in kidney development and disease.

Methodology/Principal Findings

In this report we identified and characterized a novel transgenic mouse line, MafB-GFP, which specifically marked the kidney podocytes from a very early stage of development. These mice were then used to facilitate the fluorescent activated cell sorting based purification of podocytes from embryos at E13.5 and E15.5, as well as adults. Microarrays were then used to globally define the gene expression states of podocytes at these different developmental stages. A remarkable picture emerged, identifying the multiple sets of genes that establish the neuronal, muscle, and phagocytic properties of podocytes. The complete combinatorial code of transcription factors that create the podocyte was characterized, and the global lists of growth factors and receptors they express were defined.

Conclusions/Significance

The complete molecular character of the in vivo podocyte is established for the first time. The active molecular functions and biological processes further define their unique combination of features. The results provide a resource atlas of gene expression patterns of developing and adult podocytes that will help to guide further research of these incredible cells.

The GUDMAP database--an online resource for genitourinary research

The GenitoUrinary Development Molecular Anatomy Project (GUDMAP) is an international consortium working to generate gene expression data and transgenic mice. GUDMAP includes data from large-scale in situ hybridisation screens (wholemount and section) and microarray gene expression data of microdissected, laser-captured and FACS-sorted components of the developing mouse genitourinary (GU) system. These expression data are annotated using a high-resolution anatomy ontology specific to the developing murine GU system. GUDMAP data are freely accessible at www.gudmap.org via easy-to-use interfaces. This curated, high-resolution dataset serves as a powerful resource for biologists, clinicians and bioinformaticians interested in the developing urogenital system. This paper gives examples of how the data have been used to address problems in developmental biology and provides a primer for those wishing to use the database in their own research.

Microdissection of the gene expression codes driving nephrogenesis

The kidney represents an excellent model system for learning the principles of organogenesis. It is intermediate in complexity, and employs many commonly used developmental processes. As such, kidney development has been the subject of intensive study, using a variety of techniques, including in situ hybridization, organ culture and gene targeting, revealing many critical genes and pathways. Nevertheless, proper organogenesis requires precise patterns of cell type specific differential gene expression, involving very large numbers of genes. This review is focused on the use of global profiling technologies to create an atlas of gene expression codes driving development of different mammalian kidney compartments. Such an atlas allows one to select a gene of interest, and to determine its expression level in each element of the developing kidney, or to select a structure of interest, such as the renal vesicle, and to examine its complete gene expression state. Novel component specific molecular markers are identified, and the changing waves of gene expression that drive nephrogenesis are defined. As the tools continue to improve for the purification of specific cell types and expression profiling of even individual cells it is possible to predict an atlas of gene expression during kidney development that extends to single cell resolution.

Microarrays and RNA-Seq identify molecular mechanisms driving the end of nephron production

Background

The production of nephrons suddenly ends in mice shortly after birth when the remaining cells of the multi-potent progenitor mesenchyme begin to differentiate into nephrons. We exploited this terminal wave of nephron production using both microarrays and RNA-Seq to serially evaluate gene transcript levels in the progenitors. This strategy allowed us to define the changing gene expression states following induction and the onset of differentiation after birth.

Results

Microarray and RNA-Seq studies of the progenitors detected a change in the expression profiles of several classes of genes early after birth. One functional class, a class of genes associated with cellular proliferation, was activated. Analysis of proliferation with a nucleotide analog demonstrated in vivo that entry into the S-phase of the cell cycle preceded increases in transcript levels of genetic markers of differentiation. Microarrays and RNA-Seq also detected the onset of expression of markers of differentiation within the population of progenitors prior to detectable Six2 repression. Validation by in situ hybridization demonstrated that the markers were expressed in a subset of Six2 expressing progenitors. Finally, the studies identified a third set of genes that provide indirect evidence of an altered cellular microenvironment of the multi-potential progenitors after birth.

Conclusions

These results demonstrate that Six2 expression is not sufficient to suppress activation of genes associated with growth and differentiation of nephrons. They also better define the sequence of events after induction and suggest mechanisms contributing to the rapid end of nephron production after birth in mice.

Analysis of early nephron patterning reveals a role for distal RV proliferation in fusion to the ureteric tip via a cap mesenchyme-derived connecting segment

While nephron formation is known to be initiated by a mesenchyme-to-epithelial transition of the cap mesenchyme to form a renal vesicle (RV), the subsequent patterning of the nephron and fusion with the ureteric component of the kidney to form a patent contiguous uriniferous tubule has not been fully characterized. Using dual section in situ hybridization (SISH)/immunohistochemistry (IHC) we have revealed distinct distal/proximal patterning of Notch, BMP and Wnt pathway components within the RV stage nephron. Quantitation of mitoses and Cyclin D1 expression indicated that cell proliferation was higher in the distal RV, reflecting the differential developmental programs of the proximal and distal populations. A small number of RV genes were also expressed in the early connecting segment of the nephron. Dual ISH/IHC combined with serial section immunofluorescence and 3D reconstruction revealed that fusion occurs between the late RV and adjacent ureteric tip via a process that involves loss of the intervening ureteric epithelial basement membrane and insertion of cells expressing RV markers into the ureteric tip. Using Six2-eGFPCre x R26R-lacZ mice, we demonstrate that these cells are derived from the cap mesenchyme and not the ureteric epithelium. Hence, both nephron patterning and patency are evident at the late renal vesicle stage.