Search

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.

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.

Chick Chorioallantoic Membrane (CAM) assay for UM-SCC1 cells.

Cell attachment and migration assays for UM-SCC1 cells.

DNA-PK and Rac1 activity promote UM-SCC47 cell invasion.

RNA sequencing data. Gene ontologies and the associated genes based on biological processes with the most significance based on p-values and corrected q-values.

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. Clin Cancer Res; 22(8); 2062–73. ©2015 AACR.

Generation and characterization of FA-deficient (HPV positive and HPV negative) HNSCC cell models.

Legends

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.

Background Epidemiological data associate high levels of combustion-derived particulate matter (PM) with deleterious respiratory outcomes, but the mechanism underlying those outcomes remains elusive. It has been acknowledged by the World Health Organization that PM exposure contributes to more than 4.2 million all-cause mortalities worldwide each year. Current literature demonstrates that PM exacerbates respiratory diseases, impairs lung function, results in chronic respiratory illnesses, and is associated with increased mortality. The proposed mechanisms revolve around oxidative stress and inflammation promoting pulmonary physiological remodeling. However, our previous data found that PM is capable of inducing T helper cell 17 (Th17) immune responses via aryl hydrocarbon receptor (Ahr) activation, which was associated with neutrophilic invasion characteristic of steroid insensitive asthma. Methods In the present study, we utilized a combination of microarray and single cell RNA sequencing data to analyze the immunological landscape in mouse lungs following acute exposure to combustion derived particulate matter. Results We present data that suggest epithelial cells produce specific cytokines in the aryl hydrocarbon receptor (Ahr) pathway that inform dendritic cells to initiate the production of pathogenic T helper (eTh17) cells. Using single-cell RNA sequencing analysis, we observed that upon exposure epithelial cells acquire a transcriptomic profile indicative of increased Il-17 signaling, Ahr activation, Egfr signaling, and T cell receptor and co-stimulatory signaling pathways. Epithelial cells further showed, Ahr activation is brought on by Ahr/ARNT nuclear translocation and activation of tyrosine kinase c-src, Egfr, and subsequently Erk1/2 pathways. Conclusions Collectively, our data corroborates that PM initiates an eTh17 specific inflammatory response causing neutrophilic asthma through pathways in epithelial, dendritic, and T cells that promote eTh17 differentiation during initial PM exposure.

BackgroundPersistent kidney fibroblast activation and tubular epithelial cell (TEC) injury are key contributors to CKD. However, transcriptional and cellular identities of advanced kidney disease, along with renal fibroblast specific markers and molecular targets contributing to persistent tubular injury, remain elusive.MethodsWe performed single-cell RNA sequencing with two clinically relevant murine kidney fibrosis models. Day 28 post-injury was chosen to ensure advanced fibrotic disease. Identified gene expression signatures were validated using multiple quantitative molecular analyses.ResultsWe revealed comprehensive single cell transcriptomic profiles of two independent kidney fibrosis models compared to normal control. Both models exhibited key CKD characteristics including renal blood flow decline, inflammatory expansion and proximal tubular loss. We identified novel populations including “secretory”, “migratory” and “contractile” activated fibroblasts, specifically labelled by newly identified fibroblast-specific Gucy1a3 expression. Fibrotic kidneys elicited elevated embryonic and pro-fibrotic signaling, including separate “Embryonic” and “Pro-fibrotic” TEC clusters. Also, fibrosis caused enhanced cell-to-cell crosstalk, particularly between activated fibroblasts and pro-fibrotic TECs. Analysis of factors mediating mesenchymal phenotype in the injured epithelium identified persistent elevation of Ahnak, previously reported in AKI, in both CKD models. AHNAK knockdown in primary human renal proximal tubular epithelial cells induced a pro-fibrotic phenotype and exacerbated TGFβ response via p38, p42/44, pAKT, BMP and MMP signaling.ConclusionsOur study comprehensively examined kidney fibrosis in two independent models at the singe-cell resolution, providing a valuable resource for the field. Moreover, we newly identified Gucy1a3 as a kidney activated fibroblast specific marker and validated AHNAK as a putative disease target.Significance StatementMechanistic understanding of kidney fibrosis is principal for mechanistic understanding and developing targeted strategies against CKD. However, specific markers and molecular targets of key effector cells - activated kidney fibroblasts and injured tubular epithelial cells - remain elusive. Here, we created comprehensive single cell transcriptomic profiles of two clinically relevant kidney fibrosis models. We revealed “secretory”, “contractile” and “migratory” fibroblasts and identified Gucy1a3 as a novel marker selectively labelling all three populations. We revealed that kidney fibrosis elicited remarkable epithelial-to-stromal crosstalk and pro-fibrotic signaling in the tubular cells. Moreover, we mechanistically validated AHNAK as a putative novel kidney injury target in a primary human in vitro model of epithelial-to-mesenchymal transition. Our findings advance understanding of and targeted intervention in fibrotic kidney disease.

Notch signaling in the renal interstitium is known to be required for the formation of mesangial cells and Ren1 (Renin)-expressing cells. However, little is known about how interstitial Notch signaling affects nephron development. We found that blocking Notch signaling in the renal interstitium in mice caused developmental arrest of proximal tubules accompanied by defective formation of mesangial cells. We examined the interstitial Pdgfrb mutant kidney which exhibits a similar mesangial cell defect and found that the Pdgfrb mutant kidney showed normal proximal tubule development, suggesting that the absence of mesangial cells was not the cause of defective proximal tubule development. Our single cell RNA-seq analysis of the interstitial Rbpj mutant kidney showed that a subset of proximal tubule genes were downregulated in the mutant kidney and that Gata3 was downregulated in the mutant interstitium during the development of Ren1-expressing cells. We found that deleting Gata3 in the interstitium caused the loss of Renin and the developmental arrest of proximal tubules, phenocopying the interstitial Notch/Rbpj mutants. Our results suggest that interstitial Notch signaling regulates the development of proximal tubules via the Gata3-Renin axis in the mouse kidney.

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 naive 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.

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

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.

Background Rhabdoid tumors (RTs) arise within (atypical teratoid/rhabdoid tumor [AT/RT]) or outside the brain (extra [e]CNS-RT) and are driven mainly by inactivation of the SWItch/sucrose nonfermentable (SWI/SNF) complex subunit SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily B member 1 (SMARCB1). A pathognomonic hallmark of RTs is heterogeneous multilineage differentiation, including anomalous neuronal differentiation in some eCNS-RTs. Because remodeling of the SWI/SNF complex regulates differentiation, we hypothesized that SWI/SNF Brahma-associated factors (BAF) and polybromo-associated BAF (PBAF) complex heterogeneity are related to both multilineage differentiation and clinical outcome. Methods We performed an integrated analysis of SWI/SNF complex alterations in the developing kidney and cerebellum (most common regions of RT origin) in comparison to eCNS-RT (n = 14) and AT/RT (n = 25) tumors. RT samples were interrogated using immunohistochemistry, DNA methylation, and gene expression analyses. Results The SWI/SNF BAF paralogs actin-like protein (ACTL)6A and ACTL6B were expressed in a mutually exclusive manner in the developing cerebellum and kidney. In contrast, a subset of eCNS-RTs lost mutual exclusivity and coexpressed both subunits. These tumors showed aberrant DNA methylation of genes that regulate neuronal and renal development and demonstrated immunohistochemical evidence of neuronal differentiation. In addition, low expression of the PBAF subunit polybromo-1 (PBRM1) identified a group of AT/RTs in younger children with better overall prognosis. PBRM1-low AT/RT and eCNS-RTs showed altered DNA methylation and gene expression in immune-related genes. PBRM1 knockdown resulted in lowering immunosuppressive cytokines, and PBRM1 levels in tumor samples showed an inverse relationship with cluster of differentiation (CD)8 cytotoxic T-cell infiltration. Conclusions Heterogeneity in SWI/SNF BAF (ACTL6A/ACTL6B) and PBAF (PBRM1) subunits is related to histogenesis, contributes to the immune microenvironment and prognosis in RTs, and may inform opportunities to develop immunotherapies.

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.

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

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.

BACKGROUND & AIMSThe enteric nervous system (ENS) coordinates essential intestinal functions through the concerted action of diverse enteric neurons (EN). 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.METHODSHuman myenteric ganglia and adjacent smooth muscle were isolated by laser-capture microdissection for RNA-Seq. 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-Seq (snRNA-Seq), totaling over four billion reads and 25,208 neurons. Neuronal subtypes were defined using mouse snRNA-Seq data. Comparative informatics between human and mouse datasets identified shared EN subtype markers, which were visualized in situ using hybridization chain reaction (HCR).RESULTSSeveral 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 one or more marker genes.CONCLUSIONSOur 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 GI studies.Graphical Abstract

RationaleWhile animal model studies have extensively defined mechanisms controlling cell diversity in the developing mammalian lung, the limited data available from late stage human lung development represents a significant knowledge gap. 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.MethodsSingle cell RNA sequencing generated transcriptional profiles of 5500 cells obtained from two one-day old human lungs (born at gestational ages of 39 and 31 weeks) from the LungMAP Human Tissue Core Biorepository at the University of Rochester. Frozen single cell isolates 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. Single cell sequence data from 32000 postnatal day 1, 3, 7 and 10 mouse lung (n = 2 at each time point) cells generated by the LungMAP Research Center at Cincinnati Children’s Hospital and Medical Center, using Dropseq platform, was integrated with the human data. In situ hybridization was used to confirm the spatial location of cellular phenotypes.ResultsTranscriptional interrogation of donor newborn human lung cells identified distinct clusters representing multiple populations of epithelial, endothelial, fibroblasts, pericytes, smooth muscle, and immune cells and signature genes for each of these populations were identified. Computational integration of newborn human and postnatal mouse lung development cellular transcriptomes facilitated the identification of distinct cellular lineages among all the major cell types. Integration of the human and mouse cellular transcriptomes also demonstrated cell type-specific differences in developmental states of the newborn human lung cells. In particular, matrix fibroblasts could be separated into those representative of younger cells (n=393), or older cells (n=158).This is the first comprehensive molecular map of the cellular landscape of neonatal human lung, including biomarkers for cells at distinct states of development. Our results indicate that integrated single cell RNA profiling of human and mouse lung will help identify common and species-specific mechanisms of lung development and respiratory disease.

Microglia are heterogeneous and ubiquitous 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 (scRNA-seq) analysis to compare the transcriptomes of microglia in adult mouse (C57/Bl) brains and spinal cords to identify microglial subtypes in these CNS compartments. Cortical microglia from 2-month mice consisted of a predominant population of the homeostatic subtype (HOM-M) and a small population (4%) of the inflammatory subtype (IFLAM-M), while spinal microglia consisted of 55% HOM-M and 45% IFLAM-M subtype. Comparison of cortical and spinal microglia at 2, 4 and 8 months revealed consistently a higher composition of the IFLAM-M subtype in the spinal cord. At 8-month, cortical microglia differentiated a small new subtype with interferon response phenotypes (INF-M), while spinal microglia polarized toward a proinflammatory phenotype, as indicated by the increase of microglia expressing IL-1β. To further characterize the differential plasticity of cortical and spinal microglial heterogeneity, we determined the microglial transcriptomes from HIV-1 gp 120 transgenic (Tg) mice, a model of HIV-associated neurological disorders. Compared with wild-type (Wt) cortical microglia, the gp120Tg cortical microglia had three new subtypes, with signatures of interferon I response (INF-M), cell proliferation (PLF-M), and myelination or demyelination (MYE-M) respectively; while INF-M and PLF-M subtypes presented at all ages, the MYE-M only at 4-month. In contrast, only the INF-M subtype was observed in the spinal microglia from 2- and 4-month gp120tg mice. Bioinformatic analysis of regulated molecular pathways of individual microglial subtypes indicated that gp120 more severely impaired the biological function of microglia in cortices than in the spinal cord. The results collectively reveal differential heterogeneity and plasticity of cortical and spinal microglia, and suggest functional differentiation of microglia in different CNS compartments.

SummaryAcute kidney injury (AKI) is a rapid decline of renal function, with an incidence of up to 67% of intensive care unit patients. Current treatments are merely supportive, emphasizing the need for deeper understanding that could lead to improved therapies. We used single cell RNA sequencing,in situhybridization and protein expression analyses to create comprehensive renal cell specific transcriptional profiles of multiple AKI stages. We revealed that AKI induces marked dedifferentiation, renal developmental gene activation and mixed identities in injured renal tubules. Moreover, we identified potential pathologic crosstalk between epithelial and stromal cells, and several novel genes involved in AKI. We also demonstrated the definitive effects of age on AKI outcome, and showed that renal developmental genes hold a potential as novel AKI markers. Moreover, our study provides the resource power which will aid in unraveling the molecular genetics of AKI.

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.

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.

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.

‘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.

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, while 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. Moreover, Gsx2 physically interacts with the basic helix-loop-helix (bHLH) domain of Ascl1, and DNA binding assays demonstrated that this interaction interferes with Ascl1's ability 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.

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.

The identity and degree of heterogeneity of glial progenitors and their contributions to brain tumor malignancy remain elusive. By applying lineage-targeted single-cell transcriptomics, we uncover an unanticipated diversity of glial progenitor pools with unique molecular identities in developing brain. Our analysis identifies distinct transitional intermediate states and their divergent developmental trajectories in astroglial and oligodendroglial lineages. Moreover, intersectional analysis uncovers analogous intermediate progenitors during brain tumorigenesis, wherein oligodendrocyte-progenitor intermediates are abundant, hyper-proliferative, and progressively reprogrammed toward a stem-like state susceptible to further malignant transformation. Similar actively cycling intermediate progenitors are prominent components in human gliomas with distinct driver mutations. We further unveil lineage-driving networks underlying glial fate specification and identify Zfp36l1 as necessary for oligodendrocyte-astrocyte lineage transition and glioma growth. Together, our results resolve the dynamic repertoire of common and divergent glial progenitors during development and tumorigenesis and highlight Zfp36l1 as a molecular nexus for balancing glial cell-fate decision and controlling gliomagenesis.

The urothelium is an epithelial barrier lining the bladder that protects against infection, fluid exchange and damage from toxins. The nuclear receptor Pparg promotes urothelial differentiation in vitro, and Pparg mutations are associated with bladder cancer. However, the function of Pparg in the healthy urothelium is unknown. Here we show that Pparg is critical in urothelial cells for mitochondrial biogenesis, cellular differentiation and regulation of inflammation in response to urinary tract infection (UTI). Superficial cells, which are critical for maintaining the urothelial barrier, fail to mature in Pparg mutants and basal cells undergo squamous-like differentiation. Pparg mutants display persistent inflammation after UTI, and Nf-KB, which is transiently activated in response to infection in the wild type urothelium, persists for months. Our observations suggest that in addition to its known roles in adipogegnesis and macrophage differentiation, that Pparg-dependent transcription plays a role in the urothelium controlling mitochondrial function development and regeneration., The nuclear receptor Pparg regulates urothelial differentiation in vitro but its role in healthy urothelium is unclear. Here, the authors show that PPAR gamma mediates urothelial development during both homeostasis (via mitochondrial function) and following infection, via an inflammatory response.

SUMMARY 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 mesenchymeosteoblast 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., Graphical Abstract, In Brief Calvarial suture dysgenesis is a significant human pathology. Holmeset al. transcriptionally profile the frontal suture in normal mice and mouse models of frontal suture dysgenesis at multiple embryonic ages using bulk and single-cell RNA-seq. This comprehensive dataset allows high-resolution exploration of suture development at the transcriptional and subpopulation levels.

It is now a routine to carry out single-cell RNA-Seq to define the gene expression patterns of thousands of cells, thereby revolutionizing many areas of research. Projects are underway to use these techniques to create an atlas of the expressed genes in all cell types of the human body. Here we describe cold-active protease methods for single-cell dissociation of organs and tissues that better preserve the in vivo gene expression patterns.

The respiratory system undergoes a diversity of structural, biochemical, and functional changes necessary for adaptation to air breathing at birth. To identify the heterogeneity of pulmonary cell types and dynamic changes in gene expression mediating adaptation to respiration, here we perform single cell RNA analyses of mouse lung on postnatal day 1. Using an iterative cell type identification strategy we unbiasedly identify the heterogeneity of murine pulmonary cell types. We identify distinct populations of epithelial, endothelial, mesenchymal, and immune cells, each containing distinct subpopulations. Furthermore we compare temporal changes in RNA expression patterns before and after birth to identify signaling pathways selectively activated in specific pulmonary cell types, including activation of cell stress and the unfolded protein response during perinatal adaptation of the lung. The present data provide a single cell view of the adaptation to air breathing after birth., The respiratory system is transformed in terms of functional change at birth to adapt to breathing air. Here, the authors examine the molecular changes behind the first breath in the mouse by Drop-seq based RNA sequencing, identifying activation of the unfolded protein response as a perinatal adaptation of the lung.