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.

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 (Grhl2CD-/-), 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 Grhl2CD-/- 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 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.

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.

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.

A bigenic mouse model of FSGS reveals perturbed pathways in podocytes, mesangial cells and endothelial cells.

Focal segmental glomerulosclerosis is a major cause of end stage renal disease. Many patients prove unresponsive to available therapies. An improved understanding of the molecular basis of the disease process could provide insights leading to novel therapeutic approaches. In this study we carried out an RNA-seq analysis of the altered gene expression patterns of podocytes, mesangial cells and glomerular endothelial cells of the bigenic Cd2ap+/-, Fyn-/- mutant mouse model of FSGS. In the podocytes we observed upregulation of many genes related to the Tgfß family/pathway, including Gdnf, Tgfß1, Tgfß2, Snai2, Vegfb, Bmp4, and Tnc. The mutant podocytes also showed upregulation of Acta2, a marker of smooth muscle and associated with myofibroblasts, which are implicated in driving fibrosis. GO analysis of the podocyte upregulated genes identified elevated protein kinase activity, increased expression of growth factors, and negative regulation of cell adhesion, perhaps related to the observed podocyte loss. Both podocytes and mesangial cells showed strong upregulation of aldehyde dehydrogenase genes involved in the synthesis of retinoic acid. Similarly, the Cd2ap+/-, Fyn-/- mesangial cells, as well as podocytes in other genetic models, and the glomeruli of human FSGS patients, all show upregulation of the serine protease Prss23, with the common thread suggesting important functionality. Another gene with strong upregulation in the Cd2ap+/-, Fyn-/- mutant mesangial cells as well as multiple other mutant mouse models of FSGS was thrombospondin, which activates the secreted inactive form of Tgfß. The Cd2ap+/-, Fyn-/- mutant endothelial cells showed elevated expression of genes involved in cell proliferation, angioblast migration, angiogenesis, and neovasculature, all consistent with the formation of new blood vessels in the diseased glomerulus. The resulting global definition of the perturbed molecular pathways in the three major cell types of the mutant glomerulus provide deeper understanding of the molecular pathogenic pathways.

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 Transcriptomics Uncovers Glial Progenitor Diversity and Cell Fate Determinants during Development and Gliomagenesis.

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