Resolving Cell Fate: Experimental, computational and mathematical methods in single cell transcriptomic analysis

Deeper understanding of the embryological origins of tissues and organs is likely to provide insights into novel clinically relevant preventative and therapeutic strategies. To do so effectively and at a large scale so as to have clinical significance requires an exhaustive and meticulously accurate knowledge of normal morphological development and the underlying molecular pathways. A variety of lineage tracing and classical molecular biology techniques have led to key insights but emerging technologies now offer the possibility of more fine grained and precise measurements at the level of the single cell rather than ensembles of heterogeneous cells. In this study the rapidly developing technology of single cell RNA sequencing was combined with the development of state of the art computational methods to study murine gastrulation and early organogenesis at a hitherto unprecedented granularity. A comprehensive analysis of FLK1+ cells harvested from gastrulating murine embryos was performed. In-silico reconstruction of pseudo-temporal and pseudo-spacial relations were demonstrated. Using prior knowledge of known driver genes deeper substructure was revealed in clusters that were initially defined by unsupervised algorithms, illustrating that careful implementation of supervised approaches can outperform naïve unsupervised methods. In this way multiple members of the leukotriene branch of the arachidonic acid pathway were found to be enriched within a subset of the endothelial cluster that has a molecular signature consistent with yolk sac derived definitive wave haematopoiesis. This was subsequently validated in an in-vitro embryonic stem cell differentiation colony assay. By developing a novel adaptation to tSNE dimensionality reduction that now allows new data points to be mapped backed to previously calculated embeddings, the FLK1+ data set became a reference to which cells from other experiments were mapped. The Tal1-/- knockout mutant was characterised, reaffirming the known phenotype with complete failure of embryonic haematopoiesis. Analysing the Tal1-/- endothelial cluster shows definitively in the in-vivo organism no activation of an alternative cardiac fate programme as was previously postulated from an in-vitro model system. Additionally tSNE mapped Brachyury cells into the void in the FLK1+ dataset and identified a node like population. Loss of spacial context remains an Achilles’ heel of single cell protocols. Generation of the single cell suspensions leads to disruption of cellular contacts and loss of any spacial information. A novel method is described which uses tSNE of bulk spacial data to pre- condition a tSNE map upon which single cells can then be computationally positioned reconstructing spacial context. To model gene interactions and perturbations from single-cell data, a hybrid feed-forward deep neural network was trained on branching pseudo-temporally arranged single cell qPCR data of in-vivo wild-type murine developmental haematopoiesis. Strikingly despite the model never having ‘seen’ a mutant, in-silico gene perturbations in the deep neural network are able to faithfully reproduce the Tal1-/- phenotype. In summary the use of single cell transcriptomics to probe early murine embryology com- bined with development of new methods has uncovered a novel pathway in embryonic haematopoietic development and allowed in-silico reconstruction of a short period of early embryonic haematopoiesis. Critically these methods have broad application within the fields of developmental and stem cell biology.Wellcome Trust GRANT_NUMBER: 103392/Z/13/A GRANT_NUMBER: 103392/Z/13/

Neuro-mesodermal progenitors (NMPs): a comparative study between pluripotent stem cells and embryo-derived populations.

The mammalian embryo's caudal lateral epiblast (CLE) harbours bipotent progenitors, called neural mesodermal progenitors (NMPs), that contribute to the spinal cord and the paraxial mesoderm throughout axial elongation. Here, we performed a single cell analysis of different in vitro NMP populations produced either from embryonic stem cells (ESCs) or epiblast stem cells (EpiSCs) and compared them with E8.25 CLE mouse embryos. In our analysis of this region, our findings challenge the notion that NMPs can be defined by the exclusive co-expression of Sox2 and T at mRNA level. We analyse the in vitro NMP-like populations using a purpose-built support vector machine (SVM) based on the embryo CLE and use it as a classification model to compare the in vivo and in vitro populations. Our results show that NMP differentiation from ESCs leads to heterogeneous progenitor populations with few NMP-like cells, as defined by the SVM algorithm, whereas starting with EpiSCs yields a high proportion of cells with the embryo NMP signature. We find that the population from which the Epi-NMPs are derived in culture contains a node-like population, which suggests that this population probably maintains the expression of T in vitro and thereby a source of NMPs. In conclusion, differentiation of EpiSCs into NMPs reproduces events in vivo and suggests a sequence of events for the emergence of the NMP population

The pluripotency factor NANOG controls primitive hematopoiesis and directly regulates Tal1

Progenitors of the first hematopoietic cells in the mouse arise in the early embryo from Brachyury-positive multipotent cells in the posterior-proximal region of the epiblast, but the mechanisms that specify primitive blood cells are still largely unknown. Pluripotency factors maintain uncommitted cells of the blastocyst and embryonic stem cells in the pluripotent state. However, little is known about the role played by these factors during later development, despite being expressed in the postimplantation epiblast. Using a dual transgene system for controlled expression at postimplantation stages, we found that Nanog blocks primitive hematopoiesis in the gastrulating embryo, resulting in a loss of red blood cells and downregulation of erythropoietic genes. Accordingly, Nanog-deficient embryonic stem cells are prone to erythropoietic differentiation. Moreover, Nanog expression in adults prevents the maturation of erythroid cells. By analysis of previous data for NANOG binding during stem cell differentiation and CRISPR/Cas9 genome editing, we found that Tal1 is a direct NANOG target. Our results show that Nanog regulates primitive hematopoiesis by directly repressing critical erythroid lineage specifiers.This work was supported by the Spanish government (grant BFU2014-54608-P and BFU2017-84914-P to MM; grants RYC-2011-09209 and BFU-2012-35892 to JI). The Gottgens and Nichols laboratories are supported by core funding from the Wellcome Trust and MRC to the Wellcome and MRC Cambridge Stem Cell Institute. The CNIC is supported by the Spanish Ministry of Science, Innovation and Universities (MINECO) and the Pro CNIC Foundation, and is a Severo Ochoa Center of Excellence (SEV-2015-0505)S

Nanog regulates Pou3f1 expression at the exit from pluripotency during gastrulation.

Pluripotency is regulated by a network of transcription factors that maintain early embryonic cells in an undifferentiated state while allowing them to proliferate. NANOG is a critical factor for maintaining pluripotency and its role in primordial germ cell differentiation has been well described. However, Nanog is expressed during gastrulation across all the posterior epiblast, and only later in development is its expression restricted to primordial germ cells. In this work, we unveiled a previously unknown mechanism by which Nanog specifically represses genes involved in anterior epiblast lineage. Analysis of transcriptional data from both embryonic stem cells and gastrulating mouse embryos revealed Pou3f1 expression to be negatively correlated with that of Nanog during the early stages of differentiation. We have functionally demonstrated Pou3f1 to be a direct target of NANOG by using a dual transgene system for the controlled expression of Nanog Use of Nanog null ES cells further demonstrated a role for Nanog in repressing a subset of anterior neural genes. Deletion of a NANOG binding site (BS) located nine kilobases downstream of the transcription start site of Pou3f1 revealed this BS to have a specific role in the regionalization of the expression of this gene in the embryo. Our results indicate an active role of Nanog inhibiting neural regulatory networks by repressing Pou3f1 at the onset of gastrulation.This article has an associated First Person interview with the joint first authors of the paper.This work was funded by the Spanish government [grant BFU2017-84914-P to M.M.]. The Gottgens laboratory is supported by core funding from the Wellcome Trust and Medical Research Council to the Wellcome and Medical Research Council Cambridge Stem Cell Institute. The CNIC is supported by the Instituto de Salud Carlos III (ISCIII), the Ministerio de Ciencia, Innovación y Universidades (MCNU) and the Pro CNIC Foundation, and is a Severo Ochoa Center of Excellence [SEV-2015-0505]

Decoding the regulatory network of early blood development from single-cell gene expression measurements.

Reconstruction of the molecular pathways controlling organ development has been hampered by a lack of methods to resolve embryonic progenitor cells. Here we describe a strategy to address this problem that combines gene expression profiling of large numbers of single cells with data analysis based on diffusion maps for dimensionality reduction and network synthesis from state transition graphs. Applying the approach to hematopoietic development in the mouse embryo, we map the progression of mesoderm toward blood using single-cell gene expression analysis of 3,934 cells with blood-forming potential captured at four time points between E7.0 and E8.5. Transitions between individual cellular states are then used as input to develop a single-cell network synthesis toolkit to generate a computationally executable transcriptional regulatory network model of blood development. Several model predictions concerning the roles of Sox and Hox factors are validated experimentally. Our results demonstrate that single-cell analysis of a developing organ coupled with computational approaches can reveal the transcriptional programs that underpin organogenesis.We thank J. Downing (St. Jude Children's Research Hospital, Memphis, TN, USA) for the Runx1-ires-GFP mouse. Research in the authors' laboratory is supported by the Medical Research Council, Biotechnology and Biological Sciences Research Council, Leukaemia and Lymphoma Research, the Leukemia and Lymphoma Society, Microsoft Research and core support grants by the Wellcome Trust to the Cambridge Institute for Medical Research and Wellcome Trust - MRC Cambridge Stem Cell Institute. V.M. is supported by a Medical Research Council Studentship and Centenary Award and S.W. by a Microsoft Research PhD Scholarship.This is the accepted manuscript for a paper published in Nature Biotechnology 33, 269–276 (2015) doi:10.1038/nbt.315

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.We thank M. de Bruijn, A. Martinez-Arias, J. Nichols and C. Mulas for discussion, the Cambridge Institute for Medical Research Flow Cytometry facility for their expertise in single-cell index sorting, and S. Lorenz from the Sanger Single Cell Genomics Core for supervising purification of Tal1−/− sequencing libraries. ChIP-seq reads were processed by R. Hannah. Research in the authors’ laboratories is supported by the Medical Research Council, Cancer Research UK, the Biotechnology and Biological Sciences Research Council, Bloodwise, the Leukemia and Lymphoma Society, and the Sanger-EBI Single Cell Centre, and by core support grants from the Wellcome Trust to the Cambridge Institute for Medical Research and Wellcome Trust - MRC Cambridge Stem Cell Institute and by core funding from Cancer Research UK and the European Molecular Biology Laboratory. Y.T. was supported by a fellowship from the Japan Society for the Promotion of Science. W.J. is a Wellcome Trust Clinical Research Fellow. A.S. is supported by the Sanger-EBI Single Cell Centre. This work was funded as part of Wellcome Trust Strategic Award 105031/D/14/Z ‘Tracing early mammalian lineage decisions by single-cell genomics’ awarded to W. Reik, S. Teichmann, J. Nichols, B. Simons, T. Voet, S. Srinivas, L. Vallier, B. Göttgens and J. Marioni.This is the author accepted manuscript. The final version is available from Nature Publishing Group via http://dx.doi.org/10.1038/nature1863

Defining murine organogenesis at single-cell resolution reveals a role for the leukotriene pathway in regulating blood progenitor formation.

During gastrulation, cell types from all three germ layers are specified and the basic body plan is established 1 . However, molecular analysis of this key developmental stage has been hampered by limited cell numbers and a paucity of markers. Single-cell RNA sequencing circumvents these problems, but has so far been limited to specific organ systems 2 . Here, we report single-cell transcriptomic characterization of >20,000 cells immediately following gastrulation at E8.25 of mouse development. We identify 20 major cell types, which frequently contain substructure, including three distinct signatures in early foregut cells. Pseudo-space ordering of somitic progenitor cells identifies dynamic waves of transcription and candidate regulators, which are validated by molecular characterization of spatially resolved regions of the embryo. Within the endothelial population, cells that transition from haemogenic endothelial to erythro-myeloid progenitors specifically express Alox5 and its co-factor Alox5ap, which control leukotriene production. Functional assays using mouse embryonic stem cells demonstrate that leukotrienes promote haematopoietic progenitor cell generation. Thus, this comprehensive single-cell map can be exploited to reveal previously unrecognized pathways that contribute to tissue development