Heterogeneity in Oct4 and Sox2 Targets Biases Cell Fate in 4-Cell Mouse Embryos.

The major and essential objective of pre-implantation development is to establish embryonic and extra-embryonic cell fates. To address when and how this fundamental process is initiated in mammals, we characterize transcriptomes of all individual cells throughout mouse pre-implantation development. This identifies targets of master pluripotency regulators Oct4 and Sox2 as being highly heterogeneously expressed between blastomeres of the 4-cell embryo, with Sox21 showing one of the most heterogeneous expression profiles. Live-cell tracking demonstrates that cells with decreased Sox21 yield more extra-embryonic than pluripotent progeny. Consistently, decreasing Sox21 results in premature upregulation of the differentiation regulator Cdx2, suggesting that Sox21 helps safeguard pluripotency. Furthermore, Sox21 is elevated following increased expression of the histone H3R26-methylase CARM1 and is lowered following CARM1 inhibition, indicating the importance of epigenetic regulation. Therefore, our results indicate that heterogeneous gene expression, as early as the 4-cell stage, initiates cell-fate decisions by modulating the balance of pluripotency and differentiation.

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.

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.

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.

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.

The earliest thymic T cell progenitors sustain B cell and myeloid lineage potential.

The stepwise commitment from hematopoietic stem cells in the bone marrow to T lymphocyte-restricted progenitors in the thymus represents a paradigm for understanding the requirement for distinct extrinsic cues during different stages of lineage restriction from multipotent to lineage-restricted progenitors. However, the commitment stage at which progenitors migrate from the bone marrow to the thymus remains unclear. Here we provide functional and molecular evidence at the single-cell level that the earliest progenitors in the neonatal thymus had combined granulocyte-monocyte, T lymphocyte and B lymphocyte lineage potential but not megakaryocyte-erythroid lineage potential. These potentials were identical to those of candidate thymus-seeding progenitors in the bone marrow, which were closely related at the molecular level. Our findings establish the distinct lineage-restriction stage at which the T cell lineage-commitment process transits from the bone marrow to the remote thymus.

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.

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.

Resolving early mesoderm diversification through single-cell expression profiling.

In mammals, specification of the three major germ layers occurs during gastrulation, when cells ingressing through the primitive streak differentiate into the precursor cells of major organ systems. However, the molecular mechanisms underlying this process remain unclear, as numbers of gastrulating cells are very limited. In the mouse embryo at embryonic day 6.5, cells located at the junction between the extra-embryonic region and the epiblast on the posterior side of the embryo undergo an epithelial-to-mesenchymal transition and ingress through the primitive streak. Subsequently, cells migrate, either surrounding the prospective ectoderm contributing to the embryo proper, or into the extra-embryonic region to form the yolk sac, umbilical cord and placenta. Fate mapping has shown that mature tissues such as blood and heart originate from specific regions of the pre-gastrula epiblast, but the plasticity of cells within the embryo and the function of key cell-type-specific transcription factors remain unclear. Here we analyse 1,205 cells from the epiblast and nascent Flk1(+) mesoderm of gastrulating mouse embryos using single-cell RNA sequencing, representing the first transcriptome-wide in vivo view of early mesoderm formation during mammalian gastrulation. Additionally, using knockout mice, we study the function of Tal1, a key haematopoietic transcription factor, and demonstrate, contrary to previous studies performed using retrospective assays, that Tal1 knockout does not immediately bias precursor cells towards a cardiac fate.

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.

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.

Parallel single-cell sequencing links transcriptional and epigenetic heterogeneity.

We report scM&T-seq, a method for parallel single-cell genome-wide methylome and transcriptome sequencing that allows for the discovery of associations between transcriptional and epigenetic variation. Profiling of 61 mouse embryonic stem cells confirmed known links between DNA methylation and transcription. Notably, the method revealed previously unrecognized associations between heterogeneously methylated distal regulatory elements and transcription of key pluripotency genes.

Platelet-biased stem cells reside at the apex of the haematopoietic stem-cell hierarchy.

The blood system is maintained by a small pool of haematopoietic stem cells (HSCs), which are required and sufficient for replenishing all human blood cell lineages at millions of cells per second throughout life. Megakaryocytes in the bone marrow are responsible for the continuous production of platelets in the blood, crucial for preventing bleeding--a common and life-threatening side effect of many cancer therapies--and major efforts are focused at identifying the most suitable cellular and molecular targets to enhance platelet production after bone marrow transplantation or chemotherapy. Although it has become clear that distinct HSC subsets exist that are stably biased towards the generation of lymphoid or myeloid blood cells, we are yet to learn whether other types of lineage-biased HSC exist or understand their inter-relationships and how differently lineage-biased HSCs are generated and maintained. The functional relevance of notable phenotypic and molecular similarities between megakaryocytes and bone marrow cells with an HSC cell-surface phenotype remains unclear. Here we identify and prospectively isolate a molecularly and functionally distinct mouse HSC subset primed for platelet-specific gene expression, with enhanced propensity for short- and long-term reconstitution of platelets. Maintenance of platelet-biased HSCs crucially depends on thrombopoietin, the primary extrinsic regulator of platelet development. Platelet-primed HSCs also frequently have a long-term myeloid lineage bias, can self-renew and give rise to lymphoid-biased HSCs. These findings show that HSC subtypes can be organized into a cellular hierarchy, with platelet-primed HSCs at the apex. They also demonstrate that molecular and functional priming for platelet development initiates already in a distinct HSC population. The identification of a platelet-primed HSC population should enable the rational design of therapies enhancing platelet output.

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.

G&T-seq: parallel sequencing of single-cell genomes and transcriptomes.

The simultaneous sequencing of a single cell's genome and transcriptome offers a powerful means to dissect genetic variation and its effect on gene expression. Here we describe G&T-seq, a method for separating and sequencing genomic DNA and full-length mRNA from single cells. By applying G&T-seq to over 220 single cells from mice and humans, we discovered cellular properties that could not be inferred from DNA or RNA sequencing alone.

Dynamic spatio-temporal contribution of single ß5t+ cortical epithelial precursors to the thymus medulla.

Intrathymic T-cell development is critically dependent on cortical and medullary thymic epithelial cells (TECs). Both epithelial subsets originate during early thymus organogenesis from progenitor cells that express the thymoproteasome subunit ß5t, a typical feature of cortical TECs. Using in vivo lineage fate mapping, we demonstrate in mice that ß5t(+) TEC progenitors give rise to the medullary TEC compartment early in life but significantly limit their contribution once the medulla has completely formed. Lineage-tracing studies at single cell resolution demonstrate for young mice that the postnatal medulla is expanded from individual ß5t(+) cortical progenitors located at the cortico-medullary junction. These results therefore not only define a developmental window during which the expansion of medulla is efficiently enabled by progenitors resident in the thymic cortex, but also reveal the spatio-temporal dynamics that control the growth of the thymic medulla.

Population and single-cell genomics reveal the Aire dependency, relief from Polycomb silencing, and distribution of self-antigen expression in thymic epithelia.

Promiscuous gene expression (PGE) by thymic epithelial cells (TEC) is essential for generating a diverse T cell antigen receptor repertoire tolerant to self-antigens, and thus for avoiding autoimmunity. Nevertheless, the extent and nature of this unusual expression program within TEC populations and single cells are unknown. Using deep transcriptome sequencing of carefully identified mouse TEC subpopulations, we discovered a program of PGE that is common between medullary (m) and cortical TEC, further elaborated in mTEC, and completed in mature mTEC expressing the autoimmune regulator gene (Aire). TEC populations are capable of expressing up to 19,293 protein-coding genes, the highest number of genes known to be expressed in any cell type. Remarkably, in mouse mTEC, Aire expression alone positively regulates 3980 tissue-restricted genes. Notably, the tissue specificities of these genes include known targets of autoimmunity in human AIRE deficiency. Led by the observation that genes induced by Aire expression are generally characterized by a repressive chromatin state in somatic tissues, we found these genes to be strongly associated with H3K27me3 marks in mTEC. Our findings are consistent with AIRE targeting and inducing the promiscuous expression of genes previously epigenetically silenced by Polycomb group proteins. Comparison of the transcriptomes of 174 single mTEC indicates that genes induced by Aire expression are transcribed stochastically at low cell frequency. Furthermore, when present, Aire expression-dependent transcript levels were 16-fold higher, on average, in individual TEC than in the mTEC population.

Single cell genomics: advances and future perspectives.

Advances in whole-genome and whole-transcriptome amplification have permitted the sequencing of the minute amounts of DNA and RNA present in a single cell, offering a window into the extent and nature of genomic and transcriptomic heterogeneity which occurs in both normal development and disease. Single-cell approaches stand poised to revolutionise our capacity to understand the scale of genomic, epigenomic, and transcriptomic diversity that occurs during the lifetime of an individual organism. Here, we review the major technological and biological breakthroughs achieved, describe the remaining challenges to overcome, and provide a glimpse into the promise of recent and future developments.