Deep sequencing of Phox2a nuclei reveals five classes of anterolateral system neurons.

The anterolateral system (ALS) is a major ascending pathway from the spinal cord that projects to multiple brain areas and underlies the perception of pain, itch, and skin temperature. Despite its importance, our understanding of this system has been hampered by the considerable functional and molecular diversity of its constituent cells. Here, we use fluorescence-activated cell sorting to isolate ALS neurons belonging to the Phox2a-lineage for single-nucleus RNA sequencing. We reveal five distinct clusters of ALS neurons (ALS1-5) and document their laminar distribution in the spinal cord using in situ hybridization. We identify three clusters of neurons located predominantly in laminae I-III of the dorsal horn (ALS1-3) and two clusters with cell bodies located in deeper laminae (ALS4 and ALS5). Our findings reveal the transcriptional logic that underlies ALS neuronal diversity in the adult mouse and uncover the molecular identity of two previously identified classes of projection neurons. We also show that these molecular signatures can be used to target groups of ALS neurons using retrograde viral tracing. Overall, our findings provide a valuable resource for studying somatosensory biology and targeting subclasses of ALS neurons.

Deep sequencing of Phox2a nuclei reveals five classes of anterolateral system neurons.

The anterolateral system (ALS) is a major ascending pathway from the spinal cord that projects to multiple brain areas and underlies the perception of pain, itch and skin temperature. Despite its importance, our understanding of this system has been hampered by the considerable functional and molecular diversity of its constituent cells. Here we use fluorescence-activated cell sorting to isolate ALS neurons belonging to the Phox2a-lineage for single-nucleus RNA sequencing. We reveal five distinct clusters of ALS neurons (ALS1-5) and document their laminar distribution in the spinal cord using in situ hybridization. We identify 3 clusters of neurons located predominantly in laminae I-III of the dorsal horn (ALS1-3) and two clusters with cell bodies located in deeper laminae (ALS4 & ALS5). Our findings reveal the transcriptional logic that underlies ALS neuronal diversity in the adult mouse and uncover the molecular identity of two previously identified classes of projection neurons. We also show that these molecular signatures can be used to target groups of ALS neurons using retrograde viral tracing. Overall, our findings provide a valuable resource for studying somatosensory biology and targeting subclasses of ALS neurons.

Identification of a non-canonical ciliate nuclear genetic code where UAA and UAG code for different amino acids.

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the "universal" genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3'UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids.

Perivascular niche cells sense thrombocytopenia and activate hematopoietic stem cells in an IL-1 dependent manner.

Hematopoietic stem cells (HSCs) residing in specialized niches in the bone marrow are responsible for the balanced output of multiple short-lived blood cell lineages in steady-state and in response to different challenges. However, feedback mechanisms by which HSCs, through their niches, sense acute losses of specific blood cell lineages remain to be established. While all HSCs replenish platelets, previous studies have shown that a large fraction of HSCs are molecularly primed for the megakaryocyte-platelet lineage and are rapidly recruited into proliferation upon platelet depletion. Platelets normally turnover in an activation-dependent manner, herein mimicked by antibodies inducing platelet activation and depletion. Antibody-mediated platelet activation upregulates expression of Interleukin-1 (IL-1) in platelets, and in bone marrow extracellular fluid in vivo. Genetic experiments demonstrate that rather than IL-1 directly activating HSCs, activation of bone marrow Lepr+ perivascular niche cells expressing IL-1 receptor is critical for the optimal activation of quiescent HSCs upon platelet activation and depletion. These findings identify a feedback mechanism by which activation-induced depletion of a mature blood cell lineage leads to a niche-dependent activation of HSCs to reinstate its homeostasis.

Panhematopoietic RNA barcoding enables kinetic measurements of nucleate and anucleate lineages and the activation of myeloid clones following acute platelet depletion.

BACKGROUND: Platelets and erythrocytes constitute over 95% of all hematopoietic stem cell output. However, the clonal dynamics of HSC contribution to these lineages remains largely unexplored. RESULTS: We use lentiviral genetic labeling of mouse hematopoietic stem cells to quantify output from all lineages, nucleate, and anucleate, simultaneously linking these with stem and progenitor cell transcriptomic phenotypes using single-cell RNA-sequencing. We observe dynamic shifts of clonal behaviors through time in same-animal peripheral blood and demonstrate that acute platelet depletion shifts the output of multipotent hematopoietic stem cells to the exclusive production of platelets. Additionally, we observe the emergence of new myeloid-biased clones, which support short- and long-term production of blood cells. CONCLUSIONS: Our approach enables kinetic studies of multi-lineage output in the peripheral blood and transcriptional heterogeneity of individual hematopoietic stem cells. Our results give a unique insight into hematopoietic stem cell reactivation upon platelet depletion and of clonal dynamics in both steady state and under stress.

Single-cell gene and isoform expression analysis reveals signatures of ageing in haematopoietic stem and progenitor cells.

Single-cell approaches have revealed that the haematopoietic hierarchy is a continuum of differentiation, from stem cell to committed progenitor, marked by changes in gene expression. However, many of these approaches neglect isoform-level information and thus do not capture the extent of alternative splicing within the system. Here, we present an integrated short- and long-read single-cell RNA-seq analysis of haematopoietic stem and progenitor cells. We demonstrate that over half of genes detected in standard short-read single-cell analyses are expressed as multiple, often functionally distinct, isoforms, including many transcription factors and key cytokine receptors. We observe global and HSC-specific changes in gene expression with ageing but limited impact of ageing on isoform usage. Integrating single-cell and cell-type-specific isoform landscape in haematopoiesis thus provides a new reference for comprehensive molecular profiling of heterogeneous tissues, as well as novel insights into transcriptional complexity, cell-type-specific splicing events and consequences of ageing.

Single-cell genomics reveals population structures from in vitro evolutionary studies of Salmonella.

Single-cell DNA sequencing has the potential to reveal detailed hierarchical structures in evolving populations of cells. Single cell approaches are increasingly used to study clonal evolution in human ageing and cancer but have not yet been deployed to study evolving clonal microbial populations. Here, we present an approach for single bacterial genomic analysis for in vitro evolution experiments using FACS isolation of individual bacteria followed by whole-genome amplification and sequencing. We apply this to the experimental evolution of a hypermutator strain of Salmonella in response to antibiotic stress (ciprofloxacin). By analysing sequence polymorphisms in individual cells from populations we identified the presence and prevalence of sub-populations which have acquired polymorphisms in genes previously demonstrated to be associated with ciprofloxacin susceptibility. We were also able to identify that the population exposed to antibiotic stress was able to develop resistance whilst maintaining diversity. This population structure could not be resolved from bulk sequence data, and our results show how high-throughput single-cell sequencing can enhance experimental studies of bacterial evolution.

Into the multiverse: advances in single-cell multiomic profiling.

Single-cell transcriptomic approaches have revolutionised the study of complex biological systems, with the routine measurement of gene expression in thousands of cells enabling construction of whole-organism cell atlases. However, the transcriptome is just one layer amongst many that coordinate to define cell type and state and, ultimately, function. In parallel with the widespread uptake of single-cell RNA-seq (scRNA-seq), there has been a rapid emergence of methods that enable multiomic profiling of individual cells, enabling parallel measurement of intercellular heterogeneity in the genome, epigenome, transcriptome, and proteomes. Linking measurements from each of these layers has the potential to reveal regulatory and functional mechanisms underlying cell behaviour in healthy development and disease.

Transcriptional characterization of human megakaryocyte polyploidization and lineage commitment.

BACKGROUND: Megakaryocytes (MKs) originate from cells immuno-phenotypically indistinguishable from hematopoietic stem cells (HSCs), bypassing intermediate progenitors. They mature within the adult bone marrow and release platelets into the circulation. Until now, there have been no transcriptional studies of primary human bone marrow MKs. OBJECTIVES: To characterize MKs and HSCs from human bone marrow using single-cell RNA sequencing, to investigate MK lineage commitment, maturation steps, and thrombopoiesis. RESULTS: We show that MKs at different levels of polyploidization exhibit distinct transcriptional states. Although high levels of platelet-specific gene expression occur in the lower ploidy classes, as polyploidization increases, gene expression is redirected toward translation and posttranslational processing transcriptional programs, in preparation for thrombopoiesis. Our findings are in keeping with studies of MK ultrastructure and supersede evidence generated using in vitro cultured MKs. Additionally, by analyzing transcriptional signatures of a single HSC, we identify two MK-biased HSC subpopulations exhibiting unique differentiation kinetics. We show that human bone marrow MKs originate from these HSC subpopulations, supporting the notion that they display priming for MK differentiation. Finally, to investigate transcriptional changes in MKs associated with stress thrombopoiesis, we analyzed bone marrow MKs from individuals with recent myocardial infarction and found a specific gene expression signature. Our data support the modulation of MK differentiation in this thrombotic state. CONCLUSIONS: Here, we use single-cell sequencing for the first time to characterize the human bone marrow MK transcriptome at different levels of polyploidization and investigate their differentiation from the HSC.

Non-neuronal expression of SARS-CoV-2 entry genes in the olfactory system suggests mechanisms underlying COVID-19-associated anosmia.

Altered olfactory function is a common symptom of COVID-19, but its etiology is unknown. A key question is whether SARS-CoV-2 (CoV-2) - the causal agent in COVID-19 - affects olfaction directly, by infecting olfactory sensory neurons or their targets in the olfactory bulb, or indirectly, through perturbation of supporting cells. Here we identify cell types in the olfactory epithelium and olfactory bulb that express SARS-CoV-2 cell entry molecules. Bulk sequencing demonstrated that mouse, non-human primate and human olfactory mucosa expresses two key genes involved in CoV-2 entry, ACE2 and TMPRSS2. However, single cell sequencing revealed that ACE2 is expressed in support cells, stem cells, and perivascular cells, rather than in neurons. Immunostaining confirmed these results and revealed pervasive expression of ACE2 protein in dorsally-located olfactory epithelial sustentacular cells and olfactory bulb pericytes in the mouse. These findings suggest that CoV-2 infection of non-neuronal cell types leads to anosmia and related disturbances in odor perception in COVID-19 patients.

Defining Cell Identity with Single-Cell Omics.

Cells are a fundamental unit of life, and the ability to study the phenotypes and behaviors of individual cells is crucial to understanding the workings of complex biological systems. Cell phenotypes (epigenomic, transcriptomic, proteomic, and metabolomic) exhibit dramatic heterogeneity between and within the different cell types and states underlying cellular functional diversity. Cell genotypes can also display heterogeneity throughout an organism, in the form of somatic genetic variation-most notably in the emergence and evolution of tumors. Recent technical advances in single-cell isolation and the development of omics approaches sensitive enough to reveal these aspects of cell identity have enabled a revolution in the study of multicellular systems. In this review, we discuss the technologies available to resolve the genomes, epigenomes, transcriptomes, proteomes, and metabolomes of single cells from a wide variety of living systems.

Erratum: Pluripotent state transitions coordinate morphogenesis in mouse and human embryos.

This corrects the article DOI: 10.1038/nature24675.

Pluripotent state transitions coordinate morphogenesis in mouse and human embryos.

The foundations of mammalian development lie in a cluster of embryonic epiblast stem cells. In response to extracellular matrix signalling, these cells undergo epithelialization and create an apical surface in contact with a cavity, a fundamental event for all subsequent development. Concomitantly, epiblast cells transit through distinct pluripotent states, before lineage commitment at gastrulation. These pluripotent states have been characterized at the molecular level, but their biological importance remains unclear. Here we show that exit from an unrestricted naive pluripotent state is required for epiblast epithelialization and generation of the pro-amniotic cavity in mouse embryos. Embryonic stem cells locked in the naive state are able to initiate polarization but fail to undergo lumenogenesis. Mechanistically, exit from naive pluripotency activates an Oct4-governed transcriptional program that results in expression of glycosylated sialomucin proteins and the vesicle tethering and fusion events of lumenogenesis. Similarly, exit of epiblasts from naive pluripotency in cultured human post-implantation embryos triggers amniotic cavity formation and developmental progression. Our results add tissue-level architecture as a new criterion for the characterization of different pluripotent states, and show the relevance of transitions between these states during development of the mammalian embryo.

Proliferation Drives Aging-Related Functional Decline in a Subpopulation of the Hematopoietic Stem Cell Compartment.

Aging of the hematopoietic stem cell (HSC) compartment is characterized by lineage bias and reduced stem cell function, the molecular basis of which is largely unknown. Using single-cell transcriptomics, we identified a distinct subpopulation of old HSCs carrying a p53 signature indicative of stem cell decline alongside pro-proliferative JAK/STAT signaling. To investigate the relationship between JAK/STAT and p53 signaling, we challenged HSCs with a constitutively active form of JAK2 (V617F) and observed an expansion of the p53-positive subpopulation in old mice. Our results reveal cellular heterogeneity in the onset of HSC aging and implicate a role for JAK2V617F-driven proliferation in the p53-mediated functional decline of old HSCs.

Power analysis of single-cell RNA-sequencing experiments.

Single-cell RNA sequencing (scRNA-seq) has become an established and powerful method to investigate transcriptomic cell-to-cell variation, thereby revealing new cell types and providing insights into developmental processes and transcriptional stochasticity. A key question is how the variety of available protocols compare in terms of their ability to detect and accurately quantify gene expression. Here, we assessed the protocol sensitivity and accuracy of many published data sets, on the basis of spike-in standards and uniform data processing. For our workflow, we developed a flexible tool for counting the number of unique molecular identifiers (https://github.com/vals/umis/). We compared 15 protocols computationally and 4 protocols experimentally for batch-matched cell populations, in addition to investigating the effects of spike-in molecular degradation. Our analysis provides an integrated framework for comparing scRNA-seq protocols.

Single-Cell Multiomics: Multiple Measurements from Single Cells.

Single-cell sequencing provides information that is not confounded by genotypic or phenotypic heterogeneity of bulk samples. Sequencing of one molecular type (RNA, methylated DNA or open chromatin) in a single cell, furthermore, provides insights into the cell's phenotype and links to its genotype. Nevertheless, only by taking measurements of these phenotypes and genotypes from the same single cells can such inferences be made unambiguously. In this review, we survey the first experimental approaches that assay, in parallel, multiple molecular types from the same single cell, before considering the challenges and opportunities afforded by these and future technologies.

Single-cell transcriptome analysis of fish immune cells provides insight into the evolution of vertebrate immune cell types.

The immune system of vertebrate species consists of many different cell types that have distinct functional roles and are subject to different evolutionary pressures. Here, we first analyzed conservation of genes specific for all major immune cell types in human and mouse. Our results revealed higher gene turnover and faster evolution of trans-membrane proteins in NK cells compared with other immune cell types, and especially T cells, but similar conservation of nuclear and cytoplasmic protein coding genes. To validate these findings in a distant vertebrate species, we used single-cell RNA sequencing of lck:GFP cells in zebrafish and obtained the first transcriptome of specific immune cell types in a nonmammalian species. Unsupervised clustering and single-cell TCR locus reconstruction identified three cell populations, T cells, a novel type of NK-like cells, and a smaller population of myeloid-like cells. Differential expression analysis uncovered new immune-cell-specific genes, including novel immunoglobulin-like receptors, and neofunctionalization of recently duplicated paralogs. Evolutionary analyses confirmed the higher gene turnover of trans-membrane proteins in NK cells compared with T cells in fish species, suggesting that this is a general property of immune cell types across all vertebrates.

Genome-wide in vivo screen identifies novel host regulators of metastatic colonization.

Metastasis is the leading cause of death for cancer patients. This multi-stage process requires tumour cells to survive in the circulation, extravasate at distant sites, then proliferate; it involves contributions from both the tumour cell and tumour microenvironment ('host', which includes stromal cells and the immune system). Studies suggest the early steps of the metastatic process are relatively efficient, with the post-extravasation regulation of tumour growth ('colonization') being critical in determining metastatic outcome. Here we show the results of screening 810 mutant mouse lines using an in vivo assay to identify microenvironmental regulators of metastatic colonization. We identify 23 genes that, when disrupted in mouse, modify the ability of tumour cells to establish metastatic foci, with 19 of these genes not previously demonstrated to play a role in host control of metastasis. The largest reduction in pulmonary metastasis was observed in sphingosine-1-phosphate (S1P) transporter spinster homologue 2 (Spns2)-deficient mice. We demonstrate a novel outcome of S1P-mediated regulation of lymphocyte trafficking, whereby deletion of Spns2, either globally or in a lymphatic endothelial-specific manner, creates a circulating lymphopenia and a higher percentage of effector T cells and natural killer (NK) cells present in the lung. This allows for potent tumour cell killing, and an overall decreased metastatic burden.

Initial seeding of the embryonic thymus by immune-restricted lympho-myeloid progenitors.

The final stages of restriction to the T cell lineage occur in the thymus after the entry of thymus-seeding progenitors (TSPs). The identity and lineage potential of TSPs remains unclear. Because the first embryonic TSPs enter a non-vascularized thymic rudiment, we were able to directly image and establish the functional and molecular properties of embryonic thymopoiesis-initiating progenitors (T-IPs) before their entry into the thymus and activation of Notch signaling. T-IPs did not include multipotent stem cells or molecular evidence of T cell-restricted progenitors. Instead, single-cell molecular and functional analysis demonstrated that most fetal T-IPs expressed genes of and had the potential to develop into lymphoid as well as myeloid components of the immune system. Moreover, studies of embryos deficient in the transcriptional regulator RBPJ demonstrated that canonical Notch signaling was not involved in pre-thymic restriction to the T cell lineage or the migration of T-IPs.

Separation and parallel sequencing of the genomes and transcriptomes of single cells using G&T-seq.

Parallel sequencing of a single cell's genome and transcriptome provides a powerful tool for dissecting genetic variation and its relationship with gene expression. Here we present a detailed protocol for G&T-seq, a method for separation and parallel sequencing of genomic DNA and full-length polyA(+) mRNA from single cells. We provide step-by-step instructions for the isolation and lysis of single cells; the physical separation of polyA(+) mRNA from genomic DNA using a modified oligo-dT bead capture and the respective whole-transcriptome and whole-genome amplifications; and library preparation and sequence analyses of these amplification products. The method allows the detection of thousands of transcripts in parallel with the genetic variants captured by the DNA-seq data from the same single cell. G&T-seq differs from other currently available methods for parallel DNA and RNA sequencing from single cells, as it involves physical separation of the DNA and RNA and does not require bespoke microfluidics platforms. The process can be implemented manually or through automation. When performed manually, paired genome and transcriptome sequencing libraries from eight single cells can be produced in ∼3 d by researchers experienced in molecular laboratory work. For users with experience in the programming and operation of liquid-handling robots, paired DNA and RNA libraries from 96 single cells can be produced in the same time frame. Sequence analysis and integration of single-cell G&T-seq DNA and RNA data requires a high level of bioinformatics expertise and familiarity with a wide range of informatics tools.