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.

Automated Purification of DNA Origami with SPRI Beads.

DNA origami synthesis is a well-established technique with wide-ranging applications. In most cases, the synthesized origami must be purified to remove excess materials such as DNA oligos and other functional molecules. While several purification techniques are routinely used, all have limitations, and cannot be integrated with robotic systems. Here the use of solid-phase reversible immobilization (SPRI) beads as a scalable, high-throughput, and automatable method to purify DNA origami is demonstrated. Not only can this method remove unreacted oligos and biomolecules with yields comparable to existing methods while maintaining the high structural integrity of the origami, but it can also be integrated into an automated workflow to purify simultaneously large numbers and quantities of samples. It is envisioned that the SPRI beads purification method will improve the scalability of DNA nanostructures synthesis both for research and commercial applications.

High expression of the MADS-box gene VRT2 increases the number of rudimentary basal spikelets in wheat.

Spikelets are the fundamental building blocks of Poaceae inflorescences, and their development and branching patterns determine the various inflorescence architectures and grain yield of grasses. In wheat (Triticum aestivum), the central spikelets produce the most and largest grains, while spikelet size gradually decreases acropetally and basipetally, giving rise to the characteristic lanceolate shape of wheat spikes. The acropetal gradient corresponds with the developmental age of spikelets; however, the basal spikelets are developed first, and the cause of their small size and rudimentary development is unclear. Here, we adapted G&T-seq, a low-input transcriptomics approach, to characterize gene expression profiles within spatial sections of individual spikes before and after the establishment of the lanceolate shape. We observed larger differences in gene expression profiles between the apical, central, and basal sections of a single spike than between any section belonging to consecutive developmental time points. We found that SHORT VEGETATIVE PHASE MADS-box transcription factors, including VEGETATIVE TO REPRODUCTIVE TRANSITION 2 (VRT-A2), are expressed highest in the basal section of the wheat spike and display the opposite expression gradient to flowering E-class SEPALLATA1 genes. Based on multi-year field trials and transgenic lines, we show that higher expression of VRT-A2 in the basal sections of the spike is associated with increased numbers of rudimentary basal spikelets. Our results, supported by computational modeling, suggest that the delayed transition of basal spikelets from vegetative to floral developmental programs results in the lanceolate shape of wheat spikes. This study highlights the value of spatially resolved transcriptomics to gain insights into developmental genetics pathways of grass inflorescences.

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.

ROS-mediated PI3K activation drives mitochondrial transfer from stromal cells to hematopoietic stem cells in response to infection.

Hematopoietic stem cells (HSCs) undergo rapid expansion in response to stress stimuli. Here we investigate the bioenergetic processes which facilitate the HSC expansion in response to infection. We find that infection by Gram-negative bacteria drives an increase in mitochondrial mass in mammalian HSCs, which results in a metabolic transition from glycolysis toward oxidative phosphorylation. The initial increase in mitochondrial mass occurs as a result of mitochondrial transfer from the bone marrow stromal cells (BMSCs) to HSCs through a reactive oxygen species (ROS)-dependent mechanism. Mechanistically, ROS-induced oxidative stress regulates the opening of connexin channels in a system mediated by phosphoinositide 3-kinase (PI3K) activation, which allows the mitochondria to transfer from BMSCs into HSCs. Moreover, mitochondria transfer from BMSCs into HSCs, in the response to bacterial infection, occurs before the HSCs activate their own transcriptional program for mitochondrial biogenesis. Our discovery demonstrates that mitochondrial transfer from the bone marrow microenvironment to HSCs is an early physiologic event in the mammalian response to acute bacterial infection and results in bioenergetic changes which underpin emergency granulopoiesis.

Acute myeloid leukemia induces protumoral p16INK4a-driven senescence in the bone marrow microenvironment.

Acute myeloid leukemia (AML) is an age-related disease that is highly dependent on the bone marrow (BM) microenvironment. With increasing age, tissues accumulate senescent cells, characterized by an irreversible arrest of cell proliferation and the secretion of a set of proinflammatory cytokines, chemokines, and growth factors, collectively known as the senescence-associated secretory phenotype (SASP). Here, we report that AML blasts induce a senescent phenotype in the stromal cells within the BM microenvironment and that the BM stromal cell senescence is driven by p16INK4a expression. The p16INK4a-expressing senescent stromal cells then feed back to promote AML blast survival and proliferation via the SASP. Importantly, selective elimination of p16INK4a+ senescent BM stromal cells in vivo improved the survival of mice with leukemia. Next, we find that the leukemia-driven senescent tumor microenvironment is caused by AML-induced NOX2-derived superoxide. Finally, using the p16-3MR mouse model, we show that by targeting NOX2 we reduced BM stromal cell senescence and consequently reduced AML proliferation. Together, these data identify leukemia-generated NOX2-derived superoxide as a driver of protumoral p16INK4a-dependent senescence in BM stromal cells. Our findings reveal the importance of a senescent microenvironment for the pathophysiology of leukemia. These data now open the door to investigate drugs that specifically target the "benign" senescent cells that surround and support AML.

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.