A multiplexed, single-cell sequencing screen identifies compounds that increase neurogenic reprogramming of murine Muller glia

Retinal degeneration in mammals causes permanent loss of vision, due to an inability to regenerate naturally. Some non-mammalian vertebrates show robust regeneration, via Muller glia (MG). We have recently made significant progress in stimulating adult mouse MG to regenerate functional neurons by transgenic expression of the proneural transcription factor Ascl1. While these results showed that MG can serve as an endogenous source of neuronal replacement, the efficacy of this process is limited. With the goal of improving this in mammals, we designed a small molecule screen using sci-Plex, a method to multiplex up to thousands of single nucleus RNA-seq conditions into a single experiment. We used this technology to screen a library of 92 compounds, identified, and validated two that promote neurogenesis in vivo. Our results demonstrate that high-throughput single-cell molecular profiling can substantially improve the discovery process for molecules and pathways that can stimulate neural regeneration and further demonstrate the potential for this approach to restore vision in patients with retinal disease.

Dystrophin deficiency impairs cell junction formation during embryonic myogenesis

Mutations in the DMD gene lead to Duchenne muscular dystrophy, a severe X-linked neuromuscular disorder that manifests itself as young boys acquire motor functions. DMD is typically diagnosed at 2 to 4 years of age, but the absence of dystrophin negatively impacts muscle structure and function before overt symptoms appear in patients, which poses a serious challenge in the optimization of standards of care. In this report, we investigated the early consequences of dystrophin deficiency during skeletal muscle development. We used single-cell transcriptome profiling to characterize the myogenic trajectory of human pluripotent stem cells and showed that DMD cells bifurcate to an alternative branch when they reach the somite stage. Here, dystrophin deficiency was linked to marked dysregulations of cell junction protein families involved in the cell state transitions characteristic of embryonic somitogenesis. Altogether, this work demonstrates that in vitro, dystrophin deficiency has deleterious effects on cell-cell communication during myogenic development, which should be considered in future therapeutic strategies for DMD.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or nuclear factor κB (NF-κB) inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

A single-cell time-lapse of mouse prenatal development from gastrula to birth

The house mouse (Mus musculus) is an exceptional model system, combining genetic tractability with close evolutionary affinity to humans1,2. Mouse gestation lasts only 3 weeks, during which the genome orchestrates the astonishing transformation of a single-cell zygote into a free-living pup composed of more than 500 million cells. Here, to establish a global framework for exploring mammalian development, we applied optimized single-cell combinatorial indexing3 to profile the transcriptional states of 12.4 million nuclei from 83 embryos, precisely staged at 2- to 6-hour intervals spanning late gastrulation (embryonic day 8) to birth (postnatal day 0). From these data, we annotate hundreds of cell types and explore the ontogenesis of the posterior embryo during somitogenesis and of kidney, mesenchyme, retina and early neurons. We leverage the temporal resolution and sampling depth of these whole-embryo snapshots, together with published data4-8 from earlier timepoints, to construct a rooted tree of cell-type relationships that spans the entirety of prenatal development, from zygote to birth. Throughout this tree, we systematically nominate genes encoding transcription factors and other proteins as candidate drivers of the in vivo differentiation of hundreds of cell types. Remarkably, the most marked temporal shifts in cell states are observed within one hour of birth and presumably underlie the massive physiological adaptations that must accompany the successful transition of a mammalian fetus to life outside the womb.

High-capacity sample multiplexing for single cell chromatin accessibility profiling

Single-cell chromatin accessibility has emerged as a powerful means of understanding the epigenetic landscape of diverse tissues and cell types, but profiling cells from many independent specimens is challenging and costly. Here we describe a novel approach, sciPlex-ATAC-seq, which uses unmodified DNA oligos as sample-specific nuclear labels, enabling the concurrent profiling of chromatin accessibility within single nuclei from virtually unlimited specimens or experimental conditions. We first demonstrate our method with a chemical epigenomics screen, in which we identify drug-altered distal regulatory sites predictive of compound- and dose-dependent effects on transcription. We then analyze cell type-specific chromatin changes in PBMCs from multiple donors responding to synthetic and allogeneic immune stimulation. We quantify stimulation-altered immune cell compositions and isolate the unique effects of allogeneic stimulation on chromatin accessibility specific to T-lymphocytes. Finally, we observe that impaired global chromatin decondensation often coincides with chemical inhibition of allogeneic T-cell activation.

Proteostasis governs differential temperature sensitivity across embryonic cell types

Embryonic development is remarkably robust, but temperature stress can degrade its ability to generate animals with invariant anatomy. Phenotypes associated with environmental stress suggest that some cell types are more sensitive to stress than others, but the basis of this sensitivity is unknown. Here, we characterize hundreds of individual zebrafish embryos under temperature stress using whole-animal single-cell RNA sequencing (RNA-seq) to identify cell types and molecular programs driving phenotypic variability. We find that temperature perturbs the normal proportions and gene expression programs of numerous cell types and also introduces asynchrony in developmental timing. The notochord is particularly sensitive to temperature, which we map to a specialized cell type: sheath cells. These cells accumulate misfolded protein at elevated temperature, leading to a cascading structural failure of the notochord and anatomic defects. Our study demonstrates that whole-animal single-cell RNA-seq can identify mechanisms for developmental robustness and pinpoint cell types that constitute key failure points.

Embryo-scale reverse genetics at single-cell resolution

The maturation of single-cell transcriptomic technologies has facilitated the generation of comprehensive cellular atlases from whole embryos1-4. A majority of these data, however, has been collected from wild-type embryos without an appreciation for the latent variation that is present in development. Here we present the 'zebrafish single-cell atlas of perturbed embryos': single-cell transcriptomic data from 1,812 individually resolved developing zebrafish embryos, encompassing 19 timepoints, 23 genetic perturbations and a total of 3.2 million cells. The high degree of replication in our study (eight or more embryos per condition) enables us to estimate the variance in cell type abundance organism-wide and to detect perturbation-dependent deviance in cell type composition relative to wild-type embryos. Our approach is sensitive to rare cell types, resolving developmental trajectories and genetic dependencies in the cranial ganglia neurons, a cell population that comprises less than 1% of the embryo. Additionally, time-series profiling of individual mutants identified a group of brachyury-independent cells with strikingly similar transcriptomes to notochord sheath cells, leading to new hypotheses about early origins of the skull. We anticipate that standardized collection of high-resolution, organism-scale single-cell data from large numbers of individual embryos will enable mapping of the genetic dependencies of zebrafish cell types, while also addressing longstanding challenges in developmental genetics, including the cellular and transcriptional plasticity underlying phenotypic diversity across individuals.

Single-cell, whole-embryo phenotyping of mammalian developmental disorders

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

Transcriptomic profiling of tissue environments critical for post-embryonic patterning and morphogenesis of zebrafish skin

Pigment patterns and skin appendages are prominent features of vertebrate skin. In zebrafish, regularly patterned pigment stripes and an array of calcified scales form simultaneously in the skin during post-embryonic development. Understanding the mechanisms that regulate stripe patterning and scale morphogenesis may lead to the discovery of fundamental mechanisms that govern the development of animal form. To learn about cell types and signaling interactions that govern skin patterning and morphogenesis, we generated and analyzed single-cell transcriptomes of skin from wild-type fish as well as fish having genetic or transgenically induced defects in squamation or pigmentation. These data reveal a previously undescribed population of epidermal cells that express transcripts encoding enamel matrix proteins, suggest hormonal control of epithelial-mesenchymal signaling, clarify the signaling network that governs scale papillae development, and identify a critical role for the hypodermis in supporting pigment cell development. Additionally, these comprehensive single-cell transcriptomic data representing skin phenotypes of biomedical relevance should provide a useful resource for accelerating the discovery of mechanisms that govern skin development and homeostasis.

Single-cell sequencing of individual retinal organoids reveals determinants of cell-fate heterogeneity

With a critical need for more complete in vitro models of human development and disease, organoids hold immense potential. Their complex cellular composition makes single-cell sequencing of great utility; however, the limitation of current technologies to a handful of treatment conditions restricts their use in screens or studies of organoid heterogeneity. Here, we apply sci-Plex, a single-cell combinatorial indexing (sci)-based RNA sequencing (RNA-seq) multiplexing method to retinal organoids. We demonstrate that sci-Plex and 10× methods produce highly concordant cell-class compositions and then expand sci-Plex to analyze the cell-class composition of 410 organoids upon modulation of critical developmental pathways. Leveraging individual organoid data, we develop a method to measure organoid heterogeneity, and we identify that activation of Wnt signaling early in retinal organoid cultures increases retinal cell classes up to 6 weeks later. Our data show sci-Plex's potential to dramatically scale up the analysis of treatment conditions on relevant human models.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or NF-kB inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Single-cell sequencing of individual retinal organoids reveals determinants of cell fate heterogeneity

With a critical need for more complete in vitro models of human development and disease, organoids hold immense potential. Their complex cellular composition makes single-cell sequencing of great utility; however, the limitation of current technologies to a handful of treatment conditions restricts their use in screens or studies of organoid heterogeneity. Here, we apply sci-Plex, a single-cell combinatorial indexing (sci)-based RNA-seq multiplexing method to retinal organoids. We demonstrate that sci-Plex and 10x methods produce highly concordant cell class compositions and then expand sci-Plex to analyze the cell class composition of 410 organoids upon modulation of critical developmental pathways. Leveraging individual organoid data, we develop a method to measure organoid heterogeneity, and we identify that activation of Wnt signaling early in retinal organoid cultures increases retinal cell classes up to six weeks later. Our data show sci-Plex's potential to dramatically scale-up the analysis of treatment conditions on relevant human models.

Predicting cellular responses to complex perturbations in high-throughput screens

Recent advances in multiplexed single-cell transcriptomics experiments facilitate the high-throughput study of drug and genetic perturbations. However, an exhaustive exploration of the combinatorial perturbation space is experimentally unfeasible. Therefore, computational methods are needed to predict, interpret, and prioritize perturbations. Here, we present the compositional perturbation autoencoder (CPA), which combines the interpretability of linear models with the flexibility of deep-learning approaches for single-cell response modeling. CPA learns to in silico predict transcriptional perturbation response at the single-cell level for unseen dosages, cell types, time points, and species. Using newly generated single-cell drug combination data, we validate that CPA can predict unseen drug combinations while outperforming baseline models. Additionally, the architecture's modularity enables incorporating the chemical representation of the drugs, allowing the prediction of cellular response to completely unseen drugs. Furthermore, CPA is also applicable to genetic combinatorial screens. We demonstrate this by imputing in silico 5,329 missing combinations (97.6% of all possibilities) in a single-cell Perturb-seq experiment with diverse genetic interactions. We envision CPA will facilitate efficient experimental design and hypothesis generation by enabling in silico response prediction at the single-cell level and thus accelerate therapeutic applications using single-cell technologies.

A single-cell transcriptional timelapse of mouse embryonic development, from gastrula to pup

The house mouse, Mus musculus, is an exceptional model system, combining genetic tractability with close homology to human biology. Gestation in mouse development lasts just under three weeks, a period during which its genome orchestrates the astonishing transformation of a single cell zygote into a free-living pup composed of >500 million cells. Towards a global framework for exploring mammalian development, we applied single cell combinatorial indexing (sci-*) to profile the transcriptional states of 12.4 million nuclei from 83 precisely staged embryos spanning late gastrulation (embryonic day 8 or E8) to birth (postnatal day 0 or P0), with 2-hr temporal resolution during somitogenesis, 6-hr resolution through to birth, and 20-min resolution during the immediate postpartum period. From these data (E8 to P0), we annotate dozens of trajectories and hundreds of cell types and perform deeper analyses of the unfolding of the posterior embryo during somitogenesis as well as the ontogenesis of the kidney, mesenchyme, retina, and early neurons. Finally, we leverage the depth and temporal resolution of these whole embryo snapshots, together with other published data, to construct and curate a rooted tree of cell type relationships that spans mouse development from zygote to pup. Throughout this tree, we systematically nominate sets of transcription factors (TFs) and other genes as candidate drivers of the in vivo differentiation of hundreds of mammalian cell types. Remarkably, the most dramatic shifts in transcriptional state are observed in a restricted set of cell types in the hours immediately following birth, and presumably underlie the massive changes in physiology that must accompany the successful transition of a placental mammal to extrauterine life.

Transcriptional Activation of Regenerative Hematopoiesis via Vascular Niche Sensing

Transition between activation and quiescence programs in hematopoietic stem and progenitor cells (HSC/HSPCs) is perceived to be governed intrinsically and by microenvironmental co-adaptation. However, HSC programs dictating both transition and adaptability, remain poorly defined. Single cell multiome analysis divulging differential transcriptional activity between distinct HSPC states, indicated for the exclusive absence of Fli-1 motif from quiescent HSCs. We reveal that Fli-1 activity is essential for HSCs during regenerative hematopoiesis. Fli-1 directs activation programs while manipulating cellular sensory and output machineries, enabling HSPCs co-adoptability with a stimulated vascular niche. During regenerative conditions, Fli-1 presets and enables propagation of niche-derived Notch1 signaling. Constitutively induced Notch1 signaling is sufficient to recuperate functional HSC impairments in the absence of Fli-1. Applying FLI-1 modified-mRNA transduction into lethargic adult human mobilized HSPCs, enables their vigorous niche-mediated expansion along with superior engraftment capacities. Thus, decryption of stem cell activation programs offers valuable insights for immune regenerative medicine.

Modulation of FGF pathway signaling and vascular differentiation using designed oligomeric assemblies

Growth factors and cytokines signal by binding to the extracellular domains of their receptors and drive association and transphosphorylation of the receptor intracellular tyrosine kinase domains, initiating downstream signaling cascades. To enable systematic exploration of how receptor valency and geometry affects signaling outcomes, we designed cyclic homo-oligomers with up to 8 subunits using repeat protein building blocks that can be modularly extended. By incorporating a de novo designed fibroblast growth-factor receptor (FGFR) binding module into these scaffolds, we generated a series of synthetic signaling ligands that exhibit potent valency- and geometry-dependent Ca2+ release and MAPK pathway activation. The high specificity of the designed agonists reveal distinct roles for two FGFR splice variants in driving endothelial and mesenchymal cell fates during early vascular development. The ability to incorporate receptor binding domains and repeat extensions in a modular fashion makes our designed scaffolds broadly useful for probing and manipulating cellular signaling pathways.

High-Capacity Sample Multiplexing for Single Cell Chromatin Accessibility Profiling

Single-cell chromatin accessibility has emerged as a powerful means of understanding the epigenetic landscape of diverse tissues and cell types, but profiling cells from many independent specimens is challenging and costly. Here we describe a novel approach, sciPlex-ATAC-seq, which uses unmodified DNA oligos as sample-specific nuclear labels, enabling the concurrent profiling of chromatin accessibility within single nuclei from virtually unlimited specimens or experimental conditions. We first demonstrate our method with a chemical epigenomics screen, in which we identify drug-altered distal regulatory sites predictive of compound- and dose-dependent effects on transcription. We then analyze cell type-specific chromatin changes in PBMCs from multiple donors responding to synthetic and allogeneic immune stimulation. We quantify stimulation-altered immune cell compositions and isolate the unique effects of allogeneic stimulation on chromatin accessibility specific to T-lymphocytes. Finally, we observe that impaired global chromatin decondensation often coincides with chemical inhibition of allogeneic T-cell activation.

Pulmonary osteoclast-like cells in silica induced pulmonary fibrosis

The pathophysiology of silicosis is poorly understood, limiting development of therapies for those who have been exposed to the respirable particle. We explored the mechanisms of silica-induced pulmonary fibrosis in a mouse model using multiple modalities including whole-lung single-nucleus RNA sequencing. These analyses revealed that in addition to pulmonary inflammation and fibrosis, intratracheal silica challenge induced osteoclast-like differentiation of alveolar macrophages and recruited monocytes, driven by induction of the osteoclastogenic cytokine, receptor activator of nuclear factor-κB ligand (RANKL) in pulmonary lymphocytes and alveolar type II cells. Furthermore, anti-RANKL monoclonal antibody treatment suppressed silica-induced osteoclast-like differentiation in the lung and attenuated silica-induced pulmonary fibrosis. We conclude that silica induces osteoclast-like differentiation of distinct recruited and tissue resident monocyte populations, leading to progressive lung injury, likely due to sustained elaboration of bone resorbing proteases and hydrochloric acid. Interrupting osteoclast-like differentiation may therefore constitute a promising avenue for moderating lung damage in silicosis.

Single-cell transcriptomic profiling of the zebrafish inner ear reveals molecularly distinct hair cell and supporting cell subtypes

A major cause of human deafness and vestibular dysfunction is permanent loss of the mechanosensory hair cells of the inner ear. In non-mammalian vertebrates such as zebrafish, regeneration of missing hair cells can occur throughout life. While a comparative approach has the potential to reveal the basis of such differential regenerative ability, the degree to which the inner ears of fish and mammals share common hair cells and supporting cell types remains unresolved. Here, we perform single-cell RNA sequencing of the zebrafish inner ear at embryonic through adult stages to catalog the diversity of hair cells and non-sensory supporting cells. We identify a putative progenitor population for hair cells and supporting cells, as well as distinct hair and supporting cell types in the maculae versus cristae. The hair cell and supporting cell types differ from those described for the lateral line system, a distributed mechanosensory organ in zebrafish in which most studies of hair cell regeneration have been conducted. In the maculae, we identify two subtypes of hair cells that share gene expression with mammalian striolar or extrastriolar hair cells. In situ hybridization reveals that these hair cell subtypes occupy distinct spatial domains within the three macular organs, the utricle, saccule, and lagena, consistent with the reported distinct electrophysiological properties of hair cells within these domains. These findings suggest that primitive specialization of spatially distinct striolar and extrastriolar hair cells likely arose in the last common ancestor of fish and mammals. The similarities of inner ear cell type composition between fish and mammals validate zebrafish as a relevant model for understanding inner ear-specific hair cell function and regeneration.

Optimized single-nucleus transcriptional profiling by combinatorial indexing

Single-cell combinatorial indexing RNA sequencing (sci-RNA-seq) is a powerful method for recovering gene expression data from an exponentially scalable number of individual cells or nuclei. However, sci-RNA-seq is a complex protocol that has historically exhibited variable performance on different tissues, as well as lower sensitivity than alternative methods. Here, we report a simplified, optimized version of the sci-RNA-seq protocol with three rounds of split-pool indexing that is faster, more robust and more sensitive and has a higher yield than the original protocol, with reagent costs on the order of 1 cent per cell or less. The total hands-on time from nuclei isolation to final library preparation takes 2-3 d, depending on the number of samples sharing the experiment. The improvements also allow RNA profiling from tissues rich in RNases like older mouse embryos or adult tissues that were problematic for the original method. We showcase the optimized protocol via whole-organism analysis of an E16.5 mouse embryo, profiling ~380,000 nuclei in a single experiment. Finally, we introduce a 'Tiny-Sci' protocol for experiments in which input material is very limited.

Single-cell analyses reveal early thymic progenitors and pre-B cells in zebrafish

The zebrafish has proven to be a valuable model organism for studying hematopoiesis, but relatively little is known about zebrafish immune cell development and functional diversity. Elucidating key aspects of zebrafish lymphocyte development and exploring the breadth of effector functions would provide valuable insight into the evolution of adaptive immunity. We performed single-cell RNA sequencing on ∼70,000 cells from the zebrafish marrow and thymus to establish a gene expression map of zebrafish immune cell development. We uncovered rich cellular diversity in the juvenile and adult zebrafish thymus, elucidated B- and T-cell developmental trajectories, and transcriptionally characterized subsets of hematopoietic stem and progenitor cells and early thymic progenitors. Our analysis permitted the identification of two dendritic-like cell populations and provided evidence in support of the existence of a pre-B cell state. Our results provide critical insights into the landscape of zebrafish immunology and offer a foundation for cellular and genetic studies.

The continuum of Drosophila embryonic development at single-cell resolution

Drosophila melanogaster is a powerful, long-standing model for metazoan development and gene regulation. We profiled chromatin accessibility in almost 1 million and gene expression in half a million nuclei from overlapping windows spanning the entirety of embryogenesis. Leveraging developmental asynchronicity within embryo collections, we applied deep neural networks to infer the age of each nucleus, resulting in continuous, multimodal views of molecular and cellular transitions in absolute time. We identify cell lineages; infer their developmental relationships; and link dynamic changes in enhancer usage, transcription factor (TF) expression, and the accessibility of TFs' cognate motifs. With these data, the dynamics of enhancer usage and gene expression can be explored within and across lineages at the scale of minutes, including for precise transitions like zygotic genome activation.

Nuclear oligo hashing improves differential analysis of single-cell RNA-seq

Single-cell RNA sequencing (scRNA-seq) offers a high-resolution molecular view into complex tissues, but suffers from high levels of technical noise which frustrates efforts to compare the gene expression programs of different cell types. "Spike-in" RNA standards help control for technical variation in scRNA-seq, but using them with recently developed, ultra-scalable scRNA-seq methods based on combinatorial indexing is not feasible. Here, we describe a simple and cost-effective method for normalizing transcript counts and subtracting technical variability that improves differential expression analysis in scRNA-seq. The method affixes a ladder of synthetic single-stranded DNA oligos to each cell that appears in its RNA-seq library. With improved normalization we explore chemical perturbations with broad or highly specific effects on gene regulation, including RNA pol II elongation, histone deacetylation, and activation of the glucocorticoid receptor. Our methods reveal that inhibiting histone deacetylation prevents cells from executing their canonical program of changes following glucocorticoid stimulation.

Engineering a niche supporting hematopoietic stem cell development using integrated single-cell transcriptomics

Hematopoietic stem cells (HSCs) develop from hemogenic endothelium within embryonic arterial vessels such as the aorta of the aorta-gonad-mesonephros region (AGM). To identify the signals responsible for HSC formation, here we use single cell RNA-sequencing to simultaneously analyze the transcriptional profiles of AGM-derived cells transitioning from hemogenic endothelium to HSCs, and AGM-derived endothelial cells which provide signals sufficient to support HSC maturation and self-renewal. Pseudotemporal ordering reveals dynamics of gene expression during the hemogenic endothelium to HSC transition, identifying surface receptors specifically expressed on developing HSCs. Transcriptional profiling of niche endothelial cells identifies corresponding ligands, including those signaling to Notch receptors, VLA-4 integrin, and CXCR4, which, when integrated in an engineered platform, are sufficient to support the generation of engrafting HSCs. These studies provide a transcriptional map of the signaling interactions necessary for the development of HSCs and advance the goal of engineering HSCs for therapeutic applications.

Systematic reconstruction of cellular trajectories across mouse embryogenesis

Mammalian embryogenesis is characterized by rapid cellular proliferation and diversification. Within a few weeks, a single-cell zygote gives rise to millions of cells expressing a panoply of molecular programs. Although intensively studied, a comprehensive delineation of the major cellular trajectories that comprise mammalian development in vivo remains elusive. Here, we set out to integrate several single-cell RNA-sequencing (scRNA-seq) datasets that collectively span mouse gastrulation and organogenesis, supplemented with new profiling of ~150,000 nuclei from approximately embryonic day 8.5 (E8.5) embryos staged in one-somite increments. Overall, we define cell states at each of 19 successive stages spanning E3.5 to E13.5 and heuristically connect them to their pseudoancestors and pseudodescendants. Although constructed through automated procedures, the resulting directed acyclic graph (TOME (trajectories of mammalian embryogenesis)) is largely consistent with our contemporary understanding of mammalian development. We leverage TOME to systematically nominate transcription factors (TFs) as candidate regulators of each cell type's specification, as well as 'cell-type homologs' across vertebrate evolution.

A murine model of hereditary pulmonary alveolar proteinosis caused by homozygous Csf2ra gene disruption

Hereditary pulmonary alveolar proteinosis (hPAP) is a rare disorder caused by recessive mutations in GM-CSF receptor subunit α/β genes (CSF2RA/CSF2RB, respectively) characterized by impaired GM-CSF-dependent surfactant clearance by alveolar macrophages (AMs) resulting in alveolar surfactant accumulation and hypoxemic respiratory failure. Because hPAP is caused by CSF2RA mutations in most patients, we created an animal model of hPAP caused by Csf2ra gene disruption (Csf2ra^-/- mice) and evaluated the effects on AMs and lungs. Macrophages from Csf2ra^-/- mice were unable to bind and clear GM-CSF, did not exhibit GM-CSF signaling, and had functional defects in phagocytosis, cholesterol clearance, and surfactant clearance. Csf2ra^-/- mice developed a time-dependent, progressive lung disease similar to hPAP in children caused by CSF2RA mutations with respect to the clinical, physiological, histopathological, biochemical abnormalities, biomarkers of PAP lung disease, and clinical course. In contrast, Csf2ra^+/- mice had functionally normal AMs and no lung disease. Pulmonary macrophage transplantation (PMT) without myeloablation resulted in long-term engraftment, restoration of GM-CSF responsiveness to AMs, and a safe and durable treatment effect that lasted for the duration of the experiment (6 mo). Results demonstrate that homozygous (but not heterozygous) Csf2ra gene ablation caused hPAP identical to hPAP in children with CSF2RA mutations, identified AMs as the cellular site of hPAP pathogenesis in Csf2ra^-/- mice, and have implications for preclinical studies supporting the translation of PMT as therapy of hPAP in humans.

Airway epithelial interferon response to SARS-CoV-2 is inferior to rhinovirus and heterologous rhinovirus infection suppresses SARS-CoV-2 replication

INTRODUCTION

Common alphacoronaviruses and human rhinoviruses (HRV) induce type I and III interferon (IFN) responses important to limiting viral replication in the airway epithelium. In contrast, highly pathogenic betacoronaviruses including SARS-CoV-2 may evade or antagonize RNA-induced IFN I/III responses.

METHODS

In airway epithelial cells (AECs) from children and older adults we compared IFN I/III responses to SARS-CoV-2 and HRV-16, and assessed whether pre-infection with HRV-16, or pretreatment with recombinant IFN-β or IFN-λ, modified SARS-CoV-2 replication. Bronchial AECs from children (ages 6-18 yrs.) and older adults (ages 60-75 yrs.) were differentiated ex vivo to generate organotypic cultures. In a biosafety level 3 (BSL-3) facility, cultures were infected with SARS-CoV-2 or HRV-16, and RNA and protein was harvested from cell lysates 96 hrs. following infection and supernatant was collected 48 and 96 hrs. following infection. In additional experiments cultures were pre-infected with HRV-16, or pre-treated with recombinant IFN-β1 or IFN-λ2 before SARS-CoV-2 infection.

RESULTS

Despite significant between-donor heterogeneity SARS-CoV-2 replicated 100 times more efficiently than HRV-16. IFNB1, INFL2, and CXCL10 gene expression and protein production following HRV-16 infection was significantly greater than following SARS-CoV-2. IFN gene expression and protein production were inversely correlated with SARS-CoV-2 replication. Treatment of cultures with recombinant IFNβ1 or IFNλ2, or pre-infection of cultures with HRV-16, markedly reduced SARS-CoV-2 replication.

DISCUSSION

In addition to marked between-donor heterogeneity in IFN responses and viral replication, SARS-CoV-2 elicits a less robust IFN response in primary AEC cultures than does rhinovirus, and heterologous rhinovirus infection, or treatment with recombinant IFN-β1 or IFN-λ2, markedly reduces SARS-CoV-2 replication.

Single cell biology-a Keystone Symposia report

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

Single-cell landscape of nuclear configuration and gene expression during stem cell differentiation and X inactivation

BACKGROUND

Mammalian development is associated with extensive changes in gene expression, chromatin accessibility, and nuclear structure. Here, we follow such changes associated with mouse embryonic stem cell differentiation and X inactivation by integrating, for the first time, allele-specific data from these three modalities obtained by high-throughput single-cell RNA-seq, ATAC-seq, and Hi-C.

RESULTS

Allele-specific contact decay profiles obtained by single-cell Hi-C clearly show that the inactive X chromosome has a unique profile in differentiated cells that have undergone X inactivation. Loss of this inactive X-specific structure at mitosis is followed by its reappearance during the cell cycle, suggesting a "bookmark" mechanism. Differentiation of embryonic stem cells to follow the onset of X inactivation is associated with changes in contact decay profiles that occur in parallel on both the X chromosomes and autosomes. Single-cell RNA-seq and ATAC-seq show evidence of a delay in female versus male cells, due to the presence of two active X chromosomes at early stages of differentiation. The onset of the inactive X-specific structure in single cells occurs later than gene silencing, consistent with the idea that chromatin compaction is a late event of X inactivation. Single-cell Hi-C highlights evidence of discrete changes in nuclear structure characterized by the acquisition of very long-range contacts throughout the nucleus. Novel computational approaches allow for the effective alignment of single-cell gene expression, chromatin accessibility, and 3D chromosome structure.

CONCLUSIONS

Based on trajectory analyses, three distinct nuclear structure states are detected reflecting discrete and profound simultaneous changes not only to the structure of the X chromosomes, but also to that of autosomes during differentiation. Our study reveals that long-range structural changes to chromosomes appear as discrete events, unlike progressive changes in gene expression and chromatin accessibility.

Multipotent progenitors and hematopoietic stem cells arise independently from hemogenic endothelium in the mouse embryo

During embryogenesis, waves of hematopoietic progenitors develop from hemogenic endothelium (HE) prior to the emergence of self-renewing hematopoietic stem cells (HSCs). Although previous studies have shown that yolk-sac-derived erythromyeloid progenitors and HSCs emerge from distinct populations of HE, it remains unknown whether the earliest lymphoid-competent progenitors, multipotent progenitors, and HSCs originate from common HE. In this study, we demonstrate by clonal assays and single-cell transcriptomics that rare HE with functional HSC potential in the early murine embryo are distinct from more abundant HE with multilineage hematopoietic potential that fail to generate HSCs. Specifically, HSC-competent HE are characterized by expression of CXCR4 surface marker and by higher expression of genes tied to arterial programs regulating HSC dormancy and self-renewal. Taken together, these findings suggest a revised model of developmental hematopoiesis in which the initial populations of multipotent progenitors and HSCs arise independently from HE with distinct phenotypic and transcriptional properties.

Embryo-scale, single-cell spatial transcriptomics

Spatial patterns of gene expression manifest at scales ranging from local (e.g., cell-cell interactions) to global (e.g., body axis patterning). However, current spatial transcriptomics methods either average local contexts or are restricted to limited fields of view. Here, we introduce sci-Space, which retains single-cell resolution while resolving spatial heterogeneity at larger scales. Applying sci-Space to developing mouse embryos, we captured approximate spatial coordinates and whole transcriptomes of about 120,000 nuclei. We identify thousands of genes exhibiting anatomically patterned expression, leverage spatial information to annotate cellular subtypes, show that cell types vary substantially in their extent of spatial patterning, and reveal correlations between pseudotime and the migratory patterns of differentiating neurons. Looking forward, we anticipate that sci-Space will facilitate the construction of spatially resolved single-cell atlases of mammalian development.

Neural G0: a quiescent-like state found in neuroepithelial-derived cells and glioma

Single‐cell RNA sequencing has emerged as a powerful tool for resolving cellular states associated with normal and maligned developmental processes. Here, we used scRNA‐seq to examine the cell cycle states of expanding human neural stem cells (hNSCs). From these data, we constructed a cell cycle classifier that identifies traditional cell cycle phases and a putative quiescent‐like state in neuroepithelial‐derived cell types during mammalian neurogenesis and in gliomas. The Neural G0 markers are enriched with quiescent NSC genes and other neurodevelopmental markers found in non‐dividing neural progenitors. Putative glioblastoma stem‐like cells were significantly enriched in the Neural G0 cell population. Neural G0 cell populations and gene expression are significantly associated with less aggressive tumors and extended patient survival for gliomas. Genetic screens to identify modulators of Neural G0 revealed that knockout of genes associated with the Hippo/Yap and p53 pathways diminished Neural G0 in vitro, resulting in faster G1 transit, down‐regulation of quiescence‐associated markers, and loss of Neural G0 gene expression. Thus, Neural G0 represents a dynamic quiescent‐like state found in neuroepithelial‐derived cells and gliomas.

A human cell atlas of fetal gene expression

The gene expression program underlying the specification of human cell types is of fundamental interest. We generated human cell atlases of gene expression and chromatin accessibility in fetal tissues. For gene expression, we applied three-level combinatorial indexing to >110 samples representing 15 organs, ultimately profiling ~4 million single cells. We leveraged the literature and other atlases to identify and annotate hundreds of cell types and subtypes, both within and across tissues. Our analyses focused on organ-specific specializations of broadly distributed cell types (such as blood, endothelial, and epithelial), sites of fetal erythropoiesis (which notably included the adrenal gland), and integration with mouse developmental atlases (such as conserved specification of blood cells). These data represent a rich resource for the exploration of in vivo human gene expression in diverse tissues and cell types.

A human cell atlas of fetal chromatin accessibility

The chromatin landscape underlying the specification of human cell types is of fundamental interest. We generated human cell atlases of chromatin accessibility and gene expression in fetal tissues. For chromatin accessibility, we devised a three-level combinatorial indexing assay and applied it to 53 samples representing 15 organs, profiling ~800,000 single cells. We leveraged cell types defined by gene expression to annotate these data and cataloged hundreds of thousands of candidate regulatory elements that exhibit cell type-specific chromatin accessibility. We investigated the properties of lineage-specific transcription factors (such as POU2F1 in neurons), organ-specific specializations of broadly distributed cell types (such as blood and endothelial), and cell type-specific enrichments of complex trait heritability. These data represent a rich resource for the exploration of in vivo human gene regulation in diverse tissues and cell types.

Sci-fate characterizes the dynamics of gene expression in single cells

Gene expression programs change over time, differentiation and development, and in response to stimuli. However, nearly all techniques for profiling gene expression in single cells do not directly capture transcriptional dynamics. In the present study, we present a method for combined single-cell combinatorial indexing and messenger RNA labeling (sci-fate), which uses combinatorial cell indexing and 4-thiouridine labeling of newly synthesized mRNA to concurrently profile the whole and newly synthesized transcriptome in each of many single cells. We used sci-fate to study the cortisol response in >6,000 single cultured cells. From these data, we quantified the dynamics of the cell cycle and glucocorticoid receptor activation, and explored their intersection. Finally, we developed software to infer and analyze cell-state transitions. We anticipate that sci-fate will be broadly applicable to quantitatively characterize transcriptional dynamics in diverse systems.

A-kinase-anchoring protein 1 (dAKAP1)-based signaling complexes coordinate local protein synthesis at the mitochondrial surface

Compartmentalization of macromolecules is a ubiquitous molecular mechanism that drives numerous cellular functions. The appropriate organization of enzymes in space and time enables the precise transmission and integration of intracellular signals. Molecular scaffolds constrain signaling enzymes to influence the regional modulation of these physiological processes. Mitochondrial targeting of protein kinases and protein phosphatases provides a means to locally control the phosphorylation status and action of proteins on the surface of this organelle. Dual-specificity protein kinase A anchoring protein 1 (dAKAP1) is a multivalent binding protein that targets protein kinase A (PKA), RNAs, and other signaling enzymes to the outer mitochondrial membrane. Many AKAPs recruit a diverse set of binding partners that coordinate a broad range of cellular processes. Here, results of MS and biochemical analyses reveal that dAKAP1 anchors additional components, including the ribonucleoprotein granule components La-related protein 4 (LARP4) and polyadenylate-binding protein 1 (PABPC1). Local translation of mRNAs at organelles is a means to spatially control the synthesis of proteins. RNA-Seq data demonstrate that dAKAP1 binds mRNAs encoding proteins required for mitochondrial metabolism, including succinate dehydrogenase. Functional studies suggest that the loss of dAKAP1-RNA interactions reduces mitochondrial electron transport chain activity. Hence, dAKAP1 plays a previously unappreciated role as a molecular interface between second messenger signaling and local protein synthesis machinery.

Trajectory analysis quantifies transcriptional plasticity during macrophage polarization

In recent years, macrophages have been shown to be tremendously plastic in both in vitro and in vivo settings; however, it remains unclear whether macrophages retain any persistent memory of past polarization states which may then impact their future repolarization to new states. Here, we perform deep transcriptomic profiling at high temporal resolution as macrophages are polarized with cytokines that drive them into "M1" and "M2" molecular states. We find through trajectory analysis of their global transcriptomic profiles that macrophages which are first polarized to M1 or M2 and then subsequently repolarized demonstrate little to no memory of their polarization history. We observe complete repolarization both from M1 to M2 and vice versa, and we find that macrophage transcriptional phenotypes are defined by the current cell microenvironment, rather than an amalgamation of past and present states.

Wnt Signaling Separates the Progenitor and Endocrine Compartments during Pancreas Development

In vitro differentiation of pluripotent cells into β cells is a promising alternative to cadaveric-islet transplantation as a cure for type 1 diabetes (T1D). During the directed differentiation of human embryonic stem cells (hESCS) by exogenous factors, numerous genes that affect the differentiation process are turned on and off autonomously. Manipulating these reactions could increase the efficiency of differentiation and provide a more complete control over the final composition of cell populations. To uncover in vitro autonomous responses, we performed single-cell RNA sequencing on hESCs as they differentiate in spherical clusters. We observed that endocrine cells and their progenitors exist beside one another in separate compartments that activate distinct genetic pathways. WNT pathway inhibition in the endocrine domain of the differentiating clusters reveals a necessary role for the WNT inhibitor APC during islet formation in vivo. Accordingly, WNT inhibition in vitro causes an increase in the proportion of differentiated endocrine cells.

High-throughput, microscope-based sorting to dissect cellular heterogeneity

Microscopy is a powerful tool for characterizing complex cellular phenotypes, but linking these phenotypes to genotype or RNA expression at scale remains challenging. Here, we present Visual Cell Sorting, a method that physically separates hundreds of thousands of live cells based on their visual phenotype. Automated imaging and phenotypic analysis directs selective illumination of Dendra2, a photoconvertible fluorescent protein expressed in live cells; these photoactivated cells are then isolated using fluorescence-activated cell sorting. First, we use Visual Cell Sorting to assess hundreds of nuclear localization sequence variants in a pooled format, identifying variants that improve nuclear localization and enabling annotation of nuclear localization sequences in thousands of human proteins. Second, we recover cells that retain normal nuclear morphologies after paclitaxel treatment, and then derive their single-cell transcriptomes to identify pathways associated with paclitaxel resistance in cancers. Unlike alternative methods, Visual Cell Sorting depends on inexpensive reagents and commercially available hardware. As such, it can be readily deployed to uncover the relationships between visual cellular phenotypes and internal states, including genotypes and gene expression programs.

SUN-LB55 Connect-seq to Superimpose Molecular on Anatomical Neural Circuit Maps

Animals exhibit instinctive behavioral and physiological responses to a variety of stressors to overcome danger and restore homeostasis. The physiological response to stress is governed by hypothalamic corticotropin-releasing hormone (CRH) neurons which regulate the hypothalamic-pituitary-adrenal axis to control blood levels of stress hormones. At present, the neural circuits and signaling mechanisms through which different stress signals are transmitted to CRH neurons are poorly understood. Here, we devised a new method, termed “Connect-Seq,” which couples single-cell transcriptomics and retrograde viral tracing to define the molecular identities of individual neurons in neural circuits. As a proof of concept, using Connect-Seq, we profiled single-cell transcriptomes of 124 brain neurons upstream of CRH neurons and identified subpopulations that are likely to communicate stress-related signals to CRH neurons. Analyses of single-cell transcriptomes for ‘fast-acting’ neurotransmitters revealed subsets of upstream neurons that expressed markers of inhibitory GABAergic neurons or excitatory glutamatergic neurons. Further analyses showed a number of other neuromodulators/neurotransmitters in upstream neurons, including acetylcholine, dopamine, histamine, and 43 different neuropeptides, each expressed in individual neurons or subsets of neurons. These findings reveal extreme molecular heterogeneity among upstream neurons and suggest the upstream neurons use diverse neurochemical messengers to transmit signals to CRH neurons. Many neurons coexpressed different neurotransmitters/neuromodulators, suggesting the co-release of neurochemical messengers. Dual labeling of brain sections verified expression of specific neuromodulators in virus-infected neurons upstream of CRH neurons in selected brain areas. Our results indicate that Connect-Seq can be applied to genetically dissect neural circuits and uncover molecular identities of neurons upstream of specific neuronal types of known function. Molecular markers identified in those neurons lay a foundation for the application of cell-specific genetic tools to investigate the functions and physiological significance of diverse neuronal subsets within complex neural circuits.

IL-2 enhances ex vivo-expanded regulatory T-cell persistence after adoptive transfer

As regulatory T cell (Treg) adoptive therapy continues to develop clinically, there is a need to determine which immunomodulatory agents pair most compatibly with Tregs to enable persistence and stabilize suppressor function. Prior work has shown that mechanistic target of rapamycin inhibition can increase the stability of thymic Tregs. In this study, we investigated the transcriptomic signatures of ex vivo-expanded Tregs after adoptive transfer in the setting of clinically relevant immunosuppression using a nonhuman primate (NHP) model as a prelude to future transplant studies. Here, we found that adding interleukin-2 (IL-2) to rapamycin in vivo supported a logarithmic increase in the half-life of adoptively transferred carboxyfluorescein diacetate succinimidyl ester-labeled, autologous NHP Tregs, effectively doubling the number of cells in the peripheral blood Treg compartment compared with Treg infusion when rapamycin was given alone. Using single-cell transcriptomics, we found that transferred ex vivo-expanded Tregs initially exhibit a gene expression signature consistent with an activated state. Moreover, those cells with the highest levels of activation also expressed genes associated with p53-mediated apoptosis. In contrast, transferred Tregs interrogated at day +20 posttransfer demonstrated a gene signature more similar to published profiles of resting Tregs. Together, these preclinical data further support combining IL-2 and rapamycin in vivo as adjunctive therapy for ex vivo-expanded adoptively transferred Tregs and suggest that the activation status of ex vivo-expanded Tregs is critical to their persistence.

The promise of single-cell genomics in plants

Single-cell genomic approaches have the potential to revolutionize the study of plant systems. Here, we highlight newly developed techniques to analyze transcriptomes at single-cell resolution. We focus on the rigorous standards necessary to generate and compare these data sets introducing analysis methods that can be applied to interpret their results. Lastly, we discuss the inherent limitations of single-cell studies and address future directions for plant single-cell genomics.

Inferring Causal Gene Regulatory Networks from Coupled Single-Cell Expression Dynamics Using Scribe

Here, we present Scribe (https://github.com/aristoteleo/Scribe-py), a toolkit for detecting and visualizing causal regulatory interactions between genes and explore the potential for single-cell experiments to power network reconstruction. Scribe employs restricted directed information to determine causality by estimating the strength of information transferred from a potential regulator to its downstream target. We apply Scribe and other leading approaches for causal network reconstruction to several types of single-cell measurements and show that there is a dramatic drop in performance for "pseudotime"-ordered single-cell data compared with true time-series data. We demonstrate that performing causal inference requires temporal coupling between measurements. We show that methods such as "RNA velocity" restore some degree of coupling through an analysis of chromaffin cell fate commitment. These analyses highlight a shortcoming in experimental and computational methods for analyzing gene regulation at single-cell resolution and suggest ways of overcoming it.

Dimensionality reduction by UMAP to visualize physical and genetic interactions

Dimensionality reduction is often used to visualize complex expression profiling data. Here, we use the Uniform Manifold Approximation and Projection (UMAP) method on published transcript profiles of 1484 single gene deletions of Saccharomyces cerevisiae. Proximity in low-dimensional UMAP space identifies groups of genes that correspond to protein complexes and pathways, and finds novel protein interactions, even within well-characterized complexes. This approach is more sensitive than previous methods and should be broadly useful as additional transcriptome datasets become available for other organisms.

Massively multiplex chemical transcriptomics at single-cell resolution

High-throughput chemical screens typically use coarse assays such as cell survival, limiting what can be learned about mechanisms of action, off-target effects, and heterogeneous responses. Here, we introduce "sci-Plex," which uses "nuclear hashing" to quantify global transcriptional responses to thousands of independent perturbations at single-cell resolution. As a proof of concept, we applied sci-Plex to screen three cancer cell lines exposed to 188 compounds. In total, we profiled ~650,000 single-cell transcriptomes across ~5000 independent samples in one experiment. Our results reveal substantial intercellular heterogeneity in response to specific compounds, commonalities in response to families of compounds, and insight into differential properties within families. In particular, our results with histone deacetylase inhibitors support the view that chromatin acts as an important reservoir of acetate in cancer cells.

The human body at cellular resolution: the NIH Human Biomolecular Atlas Program

Transformative technologies are enabling the construction of three-dimensional maps of tissues with unprecedented spatial and molecular resolution. Over the next seven years, the NIH Common Fund Human Biomolecular Atlas Program (HuBMAP) intends to develop a widely accessible framework for comprehensively mapping the human body at single-cell resolution by supporting technology development, data acquisition, and detailed spatial mapping. HuBMAP will integrate its efforts with other funding agencies, programs, consortia, and the biomedical research community at large towards the shared vision of a comprehensive, accessible three-dimensional molecular and cellular atlas of the human body, in health and under various disease conditions.

A lineage-resolved molecular atlas of C. elegans embryogenesis at single-cell resolution

Caenorhabditis elegans is an animal with few cells but a wide diversity of cell types. In this study, we characterize the molecular basis for their specification by profiling the transcriptomes of 86,024 single embryonic cells. We identify 502 terminal and preterminal cell types, mapping most single-cell transcriptomes to their exact position in C. elegans' invariant lineage. Using these annotations, we find that (i) the correlation between a cell's lineage and its transcriptome increases from middle to late gastrulation, then falls substantially as cells in the nervous system and pharynx adopt their terminal fates; (ii) multilineage priming contributes to the differentiation of sister cells at dozens of lineage branches; and (iii) most distinct lineages that produce the same anatomical cell type converge to a homogenous transcriptomic state.

A pooled single-cell genetic screen identifies regulatory checkpoints in the continuum of the epithelial-to-mesenchymal transition

Integrating single-cell trajectory analysis with pooled genetic screening could reveal the genetic architecture that guides cellular decisions in development and disease. We applied this paradigm to probe the genetic circuitry that controls epithelial-to-mesenchymal transition (EMT). We used single-cell RNA sequencing to profile epithelial cells undergoing a spontaneous spatially determined EMT in the presence or absence of transforming growth factor-β. Pseudospatial trajectory analysis identified continuous waves of gene regulation as opposed to discrete 'partial' stages of EMT. KRAS was connected to the exit from the epithelial state and the acquisition of a fully mesenchymal phenotype. A pooled single-cell CRISPR-Cas9 screen identified EMT-associated receptors and transcription factors, including regulators of KRAS, whose loss impeded progress along the EMT. Inhibiting the KRAS effector MEK and its upstream activators EGFR and MET demonstrates that interruption of key signaling events reveals regulatory 'checkpoints' in the EMT continuum that mimic discrete stages, and reconciles opposing views of the program that controls EMT.

A CLK3-HMGA2 Alternative Splicing Axis Impacts Human Hematopoietic Stem Cell Molecular Identity throughout Development

While gene expression dynamics have been extensively cataloged during hematopoietic differentiation in the adult, less is known about transcriptome diversity of human hematopoietic stem cells (HSCs) during development. To characterize transcriptional and post-transcriptional changes in HSCs during development, we leveraged high-throughput genomic approaches to profile miRNAs, lincRNAs, and mRNAs. Our findings indicate that HSCs manifest distinct alternative splicing patterns in key hematopoietic regulators. Detailed analysis of the splicing dynamics and function of one such regulator, HMGA2, identified an alternative isoform that escapes miRNA-mediated targeting. We further identified the splicing kinase CLK3 that, by regulating HMGA2 splicing, preserves HMGA2 function in the setting of an increase in let-7 miRNA levels, delineating how CLK3 and HMGA2 form a functional axis that influences HSC properties during development. Collectively, our study highlights molecular mechanisms by which alternative splicing and miRNA-mediated post-transcriptional regulation impact the molecular identity and stage-specific developmental features of human HSCs.