Resolution of cell fate decisions revealed by single-cell gene expression analysis from zygote to blastocyst

3D computational reconstruction of tissues with hollow spherical morphologies using single-cell gene expression data

Single-cell gene expression analysis has contributed to a better understanding of the transcriptional heterogeneity in a variety of model systems, including those used in research in developmental, cancer and stem cell biology. Nowadays, technological advances facilitate the generation of large gene expression data sets in high-throughput format. Strategies are needed to pertinently visualize this information in a tissue structure-related context, so as to improve data analysis and aid the drawing of meaningful conclusions. Here we describe an approach that uses spatial properties of the tissue source to enable the reconstruction of hollow sphere-shaped tissues and organs from single-cell gene expression data in 3D space. To demonstrate our method, we used cells of the mouse otocyst and the renal vesicle as examples. This protocol presents a straightforward computational expression analysis workflow, and it is implemented on the MATLAB and R statistical computing and graphics software platforms. Hands-on time for typical experiments can be <1 h using a standard desktop PC or Mac.

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.

The Roles of the Stem Cell-Controlling Sox2 Transcription Factor: from Neuroectoderm Development to Alzheimer's Disease?

Sox2 is a component of the core transcriptional regulatory network which maintains the totipotency of the cells during embryonic preimplantation period, the pluripotency of embryonic stem cells, and the multipotency of neural stem cells. This maintenance is controlled by internal loops between Sox2 and other transcription factors of the core such as Oct4, Nanog, Dax1, and Klf4, downstream proteins of extracellular ligands, epigenetic modifiers, and miRNAs. As Sox2 plays an important role in the balance between stem cells maintenance and commitment to differentiated lineages throughout the lifetime, it is supposed that Sox2 could regulate stem cells aging processes. In this review, we provide an update concerning the involvement of Sox2 in neurogenesis during normal aging and discuss its possible role in Alzheimer's disease.

Identification of cell types from single-cell transcriptomes using a novel clustering method

MOTIVATION

The recent advance of single-cell technologies has brought new insights into complex biological phenomena. In particular, genome-wide single-cell measurements such as transcriptome sequencing enable the characterization of cellular composition as well as functional variation in homogenic cell populations. An important step in the single-cell transcriptome analysis is to group cells that belong to the same cell types based on gene expression patterns. The corresponding computational problem is to cluster a noisy high dimensional dataset with substantially fewer objects (cells) than the number of variables (genes).

RESULTS

In this article, we describe a novel algorithm named shared nearest neighbor (SNN)-Cliq that clusters single-cell transcriptomes. SNN-Cliq utilizes the concept of shared nearest neighbor that shows advantages in handling high-dimensional data. When evaluated on a variety of synthetic and real experimental datasets, SNN-Cliq outperformed the state-of-the-art methods tested. More importantly, the clustering results of SNN-Cliq reflect the cell types or origins with high accuracy.

AVAILABILITY AND IMPLEMENTATION

The algorithm is implemented in MATLAB and Python. The source code can be downloaded at http://bioinfo.uncc.edu/SNNCliq.

Bioinformatics approaches to single-blastomere transcriptomics

The totipotent zygote gives rise to cells with differing identities during mouse preimplantation development. Many studies have focused on analyzing the spatio-temporal dependencies during these lineage decision processes and much has been learnt by tracing transgenic marker gene expression up to the blastocyst stage and by analyzing the effects of genetic manipulations (knockout/ overexpression) on embryo development. However, until recently, it has not been possible to get broader overviews on the gene expression networks that distinguish one cell from the other within the same embryo. With the advent of whole genome amplification methodology and microfluidics-based quantitative RT-PCR it became possible to generate transcriptomes of single cells. Here we review the current state of the art of single-cell transcriptomics applied to mouse preimplantation embryo blastomeres and summarize findings made by pioneering studies in recent years. Furthermore we use the PluriNetWork and ExprEssence to investigate cell transitions based on published data.

Reprogramming the genome to totipotency in mouse embryos

Despite investigative interest, the artificial derivation of pluripotent stem cells remains inefficient and incomplete reprogramming hinders its potential as a reliable tool in regenerative medicine. By contrast, fusion of terminally differentiated gametes at fertilization activates efficient epigenetic reprogramming to ensure totipotency of early embryos. Understanding the epigenetic mechanisms required for the transition from the fertilized egg to the embryo can improve efforts to reprogram differentiated cells to pluripotent/totipotent cells for therapeutic use. We review recent discoveries that are providing insight into the molecular mechanisms required for epigenetic reprogramming to totipotency in vivo.

Identification of cardiovascular lineage descendants at single-cell resolution

The transcriptional profiles of cardiac cells derived from murine embryos and from mouse embryonic stem cells (mESCs) have primarily been studied within a cell population. However, the characterization of gene expression in these cells at a single-cell level might demonstrate unique variations that cannot be appreciated within a cell pool. In this study, we aimed to establish a single-cell quantitative PCR platform and perform side-by-side comparison between cardiac progenitor cells (CPCs) and cardiomyocytes (CMs) derived from mESCs and mouse embryos. We first generated a reference map for cardiovascular single cells through quantifying lineage-defining genes for CPCs, CMs, smooth muscle cells (SMCs), endothelial cells (EDCs), fibroblasts and mESCs. This panel was then applied against single embryonic day 10.5 heart cells to demonstrate its ability to identify each endocardial cell and chamber-specific CM. In addition, we compared the gene expression profile of embryo- and mESC-derived CPCs and CMs at different developmental stages and showed that mESC-derived CMs are phenotypically similar to embryo-derived CMs up to the neonatal stage. Furthermore, we showed that single-cell expression assays coupled with time-lapse microscopy can resolve the identity and the lineage relationships between progenies of single cultured CPCs. With this approach, we found that mESC-derived Nkx2-5(+) CPCs preferentially become SMCs or CMs, whereas single embryo-derived Nkx2-5(+) CPCs represent two phenotypically distinct subpopulations that can become either EDCs or CMs. These results demonstrate that multiplex gene expression analysis in single cells is a powerful tool for examining the unique behaviors of individual embryo- or mESC-derived cardiac cells.

Quantification of cell identity from single-cell gene expression profiles

The definition of cell identity is a central problem in biology. While single-cell RNA-seq provides a wealth of information regarding cell states, better methods are needed to map their identity, especially during developmental transitions. Here, we use repositories of cell type-specific transcriptomes to quantify identities from single-cell RNA-seq profiles, accurately classifying cells from Arabidopsis root tips and human glioblastoma tumors. We apply our approach to single cells captured from regenerating roots following tip excision. Our technique exposes a previously uncharacterized transient collapse of identity distant from the injury site, demonstrating the biological relevance of a quantitative cell identity index.

The human protein PRR14 tethers heterochromatin to the nuclear lamina during interphase and mitotic exit

The nuclear lamina is a protein meshwork that lies under the inner nuclear membrane of metazoan cells. One function of the nuclear lamina is to organize heterochromatin at the inner nuclear periphery. However, very little is known about how heterochromatin attaches to the nuclear lamina and how such attachments are restored at mitotic exit. Here, we show that a previously unstudied human protein, PRR14, functions to tether heterochromatin to the nuclear periphery during interphase, through associations with heterochromatin protein 1 (HP1) and the nuclear lamina. During early mitosis, PRR14 is released from the nuclear lamina and chromatin and remains soluble. Strikingly, at the onset of anaphase, PRR14 is incorporated rapidly into chromatin through HP1 binding. Finally, in telophase, PRR14 relocalizes to the reforming nuclear lamina. This stepwise reassembly of PRR14 suggests a function in the selection of HP1-bound heterochromatin for reattachment to the nuclear lamina as cells exit mitosis.

Assessment of cell lineages and cell death in blastocysts by immunostaining

During the last decade it has been shown that most mammalian blastocysts consisted of three cell lineages. Immunofluorescence with multiple antibodies enables to identify each cell type allowing an easy detection of eventual defects. It is complementary to RT-PCR experiments as this technique allows to look at cell position and to analyze and count the proportions between the different cell types. Thus after any kind of embryo manipulation such as nuclear transfer (NT), the analysis of the three cell lineages by immunofluorescence will provide criteria for good or poor development.

Ex Uno Plures: Molecular Designs for Embryonic Pluripotency

Pluripotent cells in embryos are situated near the apex of the hierarchy of developmental potential. They are capable of generating all cell types of the mammalian body proper. Therefore, they are the exemplar of stem cells. In vivo, pluripotent cells exist transiently and become expended within a few days of their establishment. Yet, when explanted into artificial culture conditions, they can be indefinitely propagated in vitro as pluripotent stem cell lines. A host of transcription factors and regulatory genes are now known to underpin the pluripotent state. Nonetheless, how pluripotent cells are equipped with their vast multilineage differentiation potential remains elusive. Consensus holds that pluripotency transcription factors prevent differentiation by inhibiting the expression of differentiation genes. However, this does not explain the developmental potential of pluripotent cells. We have presented another emergent perspective, namely, that pluripotency factors function as lineage specifiers that enable pluripotent cells to differentiate into specific lineages, therefore endowing pluripotent cells with their multilineage potential. Here we provide a comprehensive overview of the developmental biology, transcription factors, and extrinsic signaling associated with pluripotent cells, and their accompanying subtypes, in vitro heterogeneity and chromatin states. Although much has been learned since the appreciation of mammalian pluripotency in the 1950s and the derivation of embryonic stem cell lines in 1981, we will specifically emphasize what currently remains unclear. However, the view that pluripotency factors capacitate differentiation, recently corroborated by experimental evidence, might perhaps address the long-standing question of how pluripotent cells are endowed with their multilineage differentiation potential.

Establishment of a rapid drug screening system based on embryonic stem cells

Embryonic stem (ES) cells have the capacity for self-renewal and differentiation into three germ layers following formation of embryonic bodies (EB). To investigate toxicity of pharmaceutical compounds, five toxic chemicals, indomethacin, dexamethasone, hydroxyurea, 5-fluorouracil, and cytosine arabinoside were applied in mouse ES cells during formation of EBs. Using microscopic evaluation, the size of EBs was reduced in a dose-dependent manner by treatment with pharmaceutical chemicals. While apoptosis-related proteins, cleaved caspase-3 and PARP, were decreased in compound-exposed EBs, necrosis-related protein (Hmgb1) was present in culture media of EBs, indicating that detection of Hmgb1 can result in activation of necrosis by pharmaceutical compounds. While pharmaceutical compounds impaired the differentiation of mES cells linked with spontaneous apoptotic cell death, it was determined that cytotoxic cell damage is necrosis-dependent in mES cells. In addition, an apoptotic transcript (Noxa mRNA) in toxicant-exposed EBs was decreased in parallel with apoptosis-related proteins. Following impairment of apoptosis, differentiation-related markers including un-differentiation (Sox2), endoderm (Hnf4), mesoderm (Bmp4), and ectoderm (Pax6) also fluctuated by treatment with pharmaceutical compounds. Taken together, the data imply that exposure to pharmaceutical compounds results in increased cell death hindering the spontaneous apoptosis of cells to undergo differentiation. Using both characteristics of ES cells like self-renewal or cellular pluripotency and potentials of ES cells for evaluation in toxicity of various compounds, the current study was conducted for establishment of a novel drug screening system beyond hidden virtues of the well-known chemicals.

BMP signalling regulates the pre-implantation development of extra-embryonic cell lineages in the mouse embryo

Pre-implantation development requires the specification and organization of embryonic and extra-embryonic lineages. The separation of these lineages takes place when asymmetric divisions generate inside and outside cells that differ in polarity, position and fate. Here we assess the global transcriptional identities of these precursor cells to gain insight into the molecular mechanisms regulating lineage segregation. Unexpectedly, this reveals that complementary components of the BMP signalling pathway are already differentially expressed after the first wave of asymmetric divisions. We investigate the role of BMP signalling by expressing dominant negative forms of Smad4 and Bmpr2, by down-regulating the pathway using RNAi against BMP ligands and by applying three different BMP inhibitors at distinct stages. This reveals that BMP signalling regulates the correct development of both extra-embryonic lineages, primitive endoderm and trophectoderm, but not the embryonic lineage, prior to implantation. Together these findings indicate multiple roles of BMP signalling in the early mouse embryo.

Single-cell technologies sharpen up mammalian stem cell research

Analysis of the mechanisms underlying cell fates requires the molecular quantification of cellular features. Classical techniques use population average readouts at single time points. However, these approaches mask cellular heterogeneity and dynamics and are limited for studying rare and heterogeneous cell populations like stem cells. Techniques for single-cell analyses, ideally allowing non-invasive quantification of molecular dynamics and cellular behaviour over time, are required for studying stem cells. Here, we review the development and application of these techniques to stem cell research.

Bifurcation analysis of single-cell gene expression data reveals epigenetic landscape

We present single-cell clustering using bifurcation analysis (SCUBA), a novel computational method for extracting lineage relationships from single-cell gene expression data and modeling the dynamic changes associated with cell differentiation. SCUBA draws techniques from nonlinear dynamics and stochastic differential equation theories, providing a systematic framework for modeling complex processes involving multilineage specifications. By applying SCUBA to analyze two complementary, publicly available datasets we successfully reconstructed the cellular hierarchy during early development of mouse embryos, modeled the dynamic changes in gene expression patterns, and predicted the effects of perturbing key transcriptional regulators on inducing lineage biases. The results were robust with respect to experimental platform differences between RT-PCR and RNA sequencing. We selectively tested our predictions in Nanog mutants and found good agreement between SCUBA predictions and the experimental data. We further extended the utility of SCUBA by developing a method to reconstruct missing temporal-order information from a typical single-cell dataset. Analysis of a hematopoietic dataset suggests that our method is effective for reconstructing gene expression dynamics during human B-cell development. In summary, SCUBA provides a useful single-cell data analysis tool that is well-suited for the investigation of developmental processes.

Characterizing heterogeneity in leukemic cells using single-cell gene expression analysis

BACKGROUND

A fundamental challenge for cancer therapy is that each tumor contains a highly heterogeneous cell population whose structure and mechanistic underpinnings remain incompletely understood. Recent advances in single-cell gene expression profiling have created new possibilities to characterize this heterogeneity and to dissect the potential intra-cancer cellular hierarchy.

RESULTS

Here, we apply single-cell analysis to systematically characterize the heterogeneity within leukemic cells using the MLL-AF9 driven mouse model of acute myeloid leukemia. We start with fluorescence-activated cell sorting analysis with seven surface markers, and extend by using a multiplexing quantitative polymerase chain reaction approach to assay the transcriptional profile of a panel of 175 carefully selected genes in leukemic cells at the single-cell level. By employing a set of computational tools we find striking heterogeneity within leukemic cells. Mapping to the normal hematopoietic cellular hierarchy identifies two distinct subtypes of leukemic cells; one similar to granulocyte/monocyte progenitors and the other to macrophage and dendritic cells. Further functional experiments suggest that these subtypes differ in proliferation rates and clonal phenotypes. Finally, co-expression network analysis reveals similarities as well as organizational differences between leukemia and normal granulocyte/monocyte progenitor networks.

CONCLUSIONS

Overall, our single-cell analysis pinpoints previously uncharacterized heterogeneity within leukemic cells and provides new insights into the molecular signatures of acute myeloid leukemia.

Molecular regulation of effector and memory T cell differentiation

Immunological memory is a cardinal feature of adaptive immunity and an important goal of vaccination strategies. Here we highlight advances in the understanding of the diverse T lymphocyte subsets that provide acute and long-term protection from infection. These include new insights into the transcription factors, and the upstream 'pioneering' factors that regulate their accessibility to key sites of gene regulation, as well as metabolic regulators that contribute to the differentiation of effector and memory subsets; ontogeny and defining characteristics of tissue-resident memory lymphocytes; and origins of the remarkable heterogeneity exhibited by activated T cells. Collectively, these findings underscore progress in delineating the underlying pathways that control diversification in T cell responses but also reveal gaps in the knowledge, as well as the challenges that arise in the application of this knowledge to rationally elicit desired T cell responses through vaccination and immunotherapy.

Identifying functional gene regulatory network phenotypes underlying single cell transcriptional variability

Recent analysis of single-cell transcriptomic data has revealed a surprising organization of the transcriptional variability pervasive across individual neurons. In response to distinct combinations of synaptic input-type, a new organization of neuronal subtypes emerged based on transcriptional states that were aligned along a gradient of correlated gene expression. Individual neurons traverse across these transcriptional states in response to cellular inputs. However, the regulatory network interactions driving these changes remain unclear. Here we present a novel fuzzy logic-based approach to infer quantitative gene regulatory network models from highly variable single-cell gene expression data. Our approach involves developing an a priori regulatory network that is then trained against in vivo single-cell gene expression data in order to identify causal gene interactions and corresponding quantitative model parameters. Simulations of the inferred gene regulatory network response to experimentally observed stimuli levels mirrored the pattern and quantitative range of gene expression across individual neurons remarkably well. In addition, the network identification results revealed that distinct regulatory interactions, coupled with differences in the regulatory network stimuli, drive the variable gene expression patterns observed across the neuronal subtypes. We also identified a key difference between the neuronal subtype-specific networks with respect to negative feedback regulation, with the catecholaminergic subtype network lacking such interactions. Furthermore, by varying regulatory network stimuli over a wide range, we identified several cases in which divergent neuronal subtypes could be driven towards similar transcriptional states by distinct stimuli operating on subtype-specific regulatory networks. Based on these results, we conclude that heterogenous single-cell gene expression profiles should be interpreted through a regulatory network modeling perspective in order to separate the contributions of network interactions from those of cellular inputs.

Single-cell transcriptional analysis to uncover regulatory circuits driving cell fate decisions in early mouse development

MOTIVATION

Transcriptional regulatory networks controlling cell fate decisions in mammalian embryonic development remain elusive despite a long time of research. The recent emergence of single-cell RNA profiling technology raises hope for new discovery. Although experimental works have obtained intriguing insights into the mouse early development, a holistic and systematic view is still missing. Mathematical models of cell fates tend to be concept-based, not designed to learn from real data. To elucidate the regulatory mechanisms behind cell fate decisions, it is highly desirable to synthesize the data-driven and knowledge-driven modeling approaches.

RESULTS

We propose a novel method that integrates the structure of a cell lineage tree with transcriptional patterns from single-cell data. This method adopts probabilistic Boolean network (PBN) for network modeling, and genetic algorithm as search strategy. Guided by the 'directionality' of cell development along branches of the cell lineage tree, our method is able to accurately infer the regulatory circuits from single-cell gene expression data, in a holistic way. Applied on the single-cell transcriptional data of mouse preimplantation development, our algorithm outperforms conventional methods of network inference. Given the network topology, our method can also identify the operational interactions in the gene regulatory network (GRN), corresponding to specific cell fate determination. This is one of the first attempts to infer GRNs from single-cell transcriptional data, incorporating dynamics of cell development along a cell lineage tree.

AVAILABILITY AND IMPLEMENTATION

Implementation of our algorithm is available from the authors upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Quantitative single-cell approaches to stem cell research

Understanding the molecular control of cell fates is central to stem cell research. Such insight requires quantification of molecular and cellular behavior at the single-cell level. Recent advances now permit high-throughput molecular readouts from single cells as well as continuous, noninvasive observation of cell behavior over time. Here, we review current state-of-the-art approaches used to query stem cell fate at the single-cell level, including advances in lineage tracing, time-lapse imaging, and molecular profiling. We also offer our perspective on the advantages and drawbacks of available approaches, key technical limitations, considerations for data interpretation, and future innovation.

Defining midbrain dopaminergic neuron diversity by single-cell gene expression profiling

Effective approaches to neuropsychiatric disorders require detailed understanding of the cellular composition and circuitry of the complex mammalian brain. Here, we present a paradigm for deconstructing the diversity of neurons defined by a specific neurotransmitter using a microfluidic dynamic array to simultaneously evaluate the expression of 96 genes in single neurons. With this approach, we successfully identified multiple molecularly distinct dopamine neuron subtypes and localized them in the adult mouse brain. To validate the anatomical and functional correlates of molecular diversity, we provide evidence that one Vip+ subtype, located in the periaqueductal region, has a discrete projection field within the extended amygdala. Another Aldh1a1+ subtype, located in the substantia nigra, is especially vulnerable in the 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) model of Parkinson's disease. Overall, this rapid, cost-effective approach enables the identification and classification of multiple dopamine neuron subtypes, with distinct molecular, anatomical, and functional properties.

The transcription repressors Bach2 and Bach1 promote B cell development by repressing the myeloid program

Mature lymphoid cells express the transcription repressor Bach2, which imposes regulation on humoral and cellular immunity. Here we found critical roles for Bach2 in the development of cells of the B lineage, commencing from the common lymphoid progenitor (CLP) stage, with Bach1 as an auxiliary. Overexpression of Bach2 in pre-pro-B cells deficient in the transcription factor EBF1 and single-cell analysis of CLPs revealed that Bach2 and Bach1 repressed the expression of genes important for myeloid cells ('myeloid genes'). Bach2 and Bach1 bound to presumptive regulatory regions of the myeloid genes. Bach2(hi) CLPs showed resistance to myeloid differentiation even when cultured under myeloid conditions. Our results suggest that Bach2 functions with Bach1 and EBF1 to promote B cell development by repressing myeloid genes in CLPs.

HIPPO pathway members restrict SOX2 to the inner cell mass where it promotes ICM fates in the mouse blastocyst

Pluripotent epiblast (EPI) cells, present in the inner cell mass (ICM) of the mouse blastocyst, are progenitors of both embryonic stem (ES) cells and the fetus. Discovering how pluripotency genes regulate cell fate decisions in the blastocyst provides a valuable way to understand how pluripotency is normally established. EPI cells are specified by two consecutive cell fate decisions. The first decision segregates ICM from trophectoderm (TE), an extraembryonic cell type. The second decision subdivides ICM into EPI and primitive endoderm (PE), another extraembryonic cell type. Here, we investigate the roles and regulation of the pluripotency gene Sox2 during blastocyst formation. First, we investigate the regulation of Sox2 patterning and show that SOX2 is restricted to ICM progenitors prior to blastocyst formation by members of the HIPPO pathway, independent of CDX2, the TE transcription factor that restricts Oct4 and Nanog to the ICM. Second, we investigate the requirement for Sox2 in cell fate specification during blastocyst formation. We show that neither maternal (M) nor zygotic (Z) Sox2 is required for blastocyst formation, nor for initial expression of the pluripotency genes Oct4 or Nanog in the ICM. Rather, Z Sox2 initially promotes development of the primitive endoderm (PE) non cell-autonomously via FGF4, and then later maintains expression of pluripotency genes in the ICM. The significance of these observations is that 1) ICM and TE genes are spatially patterned in parallel prior to blastocyst formation and 2) both the roles and regulation of Sox2 in the blastocyst are unique compared to other pluripotency factors such as Oct4 or Nanog.

Pluripotent stem cells can give rise to any cell type in the body, making them an attractive tool for regenerative medicine. Pluripotent stem cells can be derived from the mammalian embryo at the blastocyst stage or they can be created from mature adult cells by reprogramming. During reprogramming, SOX2 helps establish pluripotency, but it is not clear how SOX2 establishes pluripotency in the blastocyst. We evaluated where SOX2 is present, how SOX2 is regulated, and where SOX2 is active during blastocyst formation. Our data show that the roles and the regulation of SOX2 are unique compared to other pluripotency/reprogramming factors, such as OCT4 and NANOG. SOX2 marks pluripotent cells earlier than do other factors, but does not regulate pluripotency until several days later. Rather, the earlier role of SOX2 is to help establish the yolk sac lineage, which is essential for gestation.

Time-variant clustering model for understanding cell fate decisions

Both spatial characteristics and temporal features are often the subjects of concern in physical, social, and biological studies. This work tackles the clustering problems for time course data in which the cluster number and clustering structure change with respect to time, dubbed time-variant clustering. We developed a hierarchical model that simultaneously clusters the objects at every time point and describes the relationships of the clusters between time points. The hidden layer of this model is a generalized form of branching processes. A reversible-jump Markov Chain Monte Carlo method was implemented for model inference, and a feature selection procedure was developed. We applied this method to explore an open question in preimplantation embryonic development. Our analyses using single-cell gene expression data suggested that the earliest cell fate decision could start at the 4-cell stage in mice, earlier than the commonly thought 8- to 16-cell stage. These results together with independent experimental data from single-cell RNA-seq provided support against a prevailing hypothesis in mammalian development.

BMP signaling is required for cell cleavage in preimplantation-mouse embryos

The mechanisms regulating cell division during development of the mouse pre-implantation embryo are poorly understood. We have investigated whether bone morphogenetic protein (BMP) signaling is involved in controlling cell cycle during mouse pre-implantation development. We mapped and quantitated the dynamic activities of BMP signaling through high-resolution immunofluorescence imaging combined with a 3D segmentation method. Immunostaining for phosphorylated Smad1/5/8 shows that BMP signaling is activated in mouse embryos as early as the 4-cell stage, and becomes spatially restricted by late blastocyst stage. Perturbation of BMP signaling in preimplantation mouse embryos, whether by treatment with a small molecule inhibitor, with Noggin protein, or by overexpression of a dominant-negative BMP receptor, indicates that BMPs regulate cell cleavage up to the morula stage. These results indicate that BMP signaling is active during mouse pre-implantation development and is required for cell cleavage in preimplantation mouse embryos.

Dynamic heterogeneity and DNA methylation in embryonic stem cells

Cell populations can be strikingly heterogeneous, composed of multiple cellular states, each exhibiting stochastic noise in its gene expression. A major challenge is to disentangle these two types of variability and to understand the dynamic processes and mechanisms that control them. Embryonic stem cells (ESCs) provide an ideal model system to address this issue because they exhibit heterogeneous and dynamic expression of functionally important regulatory factors. We analyzed gene expression in individual ESCs using single-molecule RNA-FISH and quantitative time-lapse movies. These data discriminated stochastic switching between two coherent (correlated) gene expression states and burst-like transcriptional noise. We further showed that the “2i” signaling pathway inhibitors modulate both types of variation. Finally, we found that DNA methylation plays a key role in maintaining these metastable states. Together, these results show how ESC gene expression states and dynamics arise from a combination of intrinsic noise, coherent cellular states, and epigenetic regulation.

•

smFISH in ESCs reveals two transcriptional states and highly stochastic expression

•

Live-cell expression dynamics reveal the in situ transition rates between states

•

DNA methylation regulates state-switching dynamics

•

“2i” signaling inhibitors impact both gene expression noise and state transitions

Chromatin dynamics in the regulation of cell fate allocation during early embryogenesis

Following fertilization, gametes undergo epigenetic reprogramming in order to revert to a totipotent state. How embryonic cells subsequently acquire their fate and the role of chromatin dynamics in this process are unknown. Genetic and experimental embryology approaches have identified some of the players and morphological changes that are involved in early mammalian development, but the exact events underlying cell fate allocation in single embryonic cells have remained elusive. Experimental and technological advances have recently provided novel insights into chromatin dynamics and nuclear architecture in single cells; these insights have reshaped our understanding of the mechanisms underlying cell fate allocation and plasticity in early mammalian development.

Gata6, Nanog and Erk signaling control cell fate in the inner cell mass through a tristable regulatory network

During blastocyst formation, inner cell mass (ICM) cells differentiate into either epiblast (Epi) or primitive endoderm (PrE) cells, labeled by Nanog and Gata6, respectively, and organized in a salt-and-pepper pattern. Previous work in the mouse has shown that, in absence of Nanog, all ICM cells adopt a PrE identity. Moreover, the activation or the blockade of the Fgf/RTK pathway biases cell fate specification towards either PrE or Epi, respectively. We show that, in absence of Gata6, all ICM cells adopt an Epi identity. Furthermore, the analysis of Gata6+/− embryos reveals a dose-sensitive phenotype, with fewer PrE-specified cells. These results and previous findings have enabled the development of a mathematical model for the dynamics of the regulatory network that controls ICM differentiation into Epi or PrE cells. The model describes the temporal dynamics of Erk signaling and of the concentrations of Nanog, Gata6, secreted Fgf4 and Fgf receptor 2. The model is able to recapitulate most of the cell behaviors observed in different experimental conditions and provides a unifying mechanism for the dynamics of these developmental transitions. The mechanism relies on the co-existence between three stable steady states (tristability), which correspond to ICM, Epi and PrE cells, respectively. Altogether, modeling and experimental results uncover novel features of ICM cell fate specification such as the role of the initial induction of a subset of cells into Epi in the initiation of the salt-and-pepper pattern, or the precocious Epi specification in Gata6+/− embryos.

A new toolbox for assessing single cells

Unprecedented access to the biology of single cells is now feasible, enabled by recent technological advancements that allow us to manipulate and measure sparse samples and achieve a new level of resolution in space and time. This review focuses on advances in tools to study single cells for specific areas of biology. We examine both mature and nascent techniques to study single cells at the genomics, transcriptomics, and proteomics level. In addition, we provide an overview of tools that are well suited for following biological responses to defined perturbations with single-cell resolution. Techniques to analyze and manipulate single cells through soluble and chemical ligands, the microenvironment, and cell-cell interactions are provided. For each of these topics, we highlight the biological motivation, applications, methods, recent advances, and opportunities for improvement. The toolbox presented in this review can function as a starting point for the design of single-cell experiments.

Single-cell analyses of regulatory network perturbations using enhancer-targeting TALEs suggest novel roles for PU.1 during haematopoietic specification

Transcription factors (TFs) act within wider regulatory networks to control cell identity and fate. Numerous TFs, including Scl (Tal1) and PU.1 (Spi1), are known regulators of developmental and adult haematopoiesis, but how they act within wider TF networks is still poorly understood. Transcription activator-like effectors (TALEs) are a novel class of genetic tool based on the modular DNA-binding domains of Xanthomonas TAL proteins, which enable DNA sequence-specific targeting and the manipulation of endogenous gene expression. Here, we report TALEs engineered to target the PU.1-14kb and Scl+40kb transcriptional enhancers as efficient new tools to perturb the expression of these key haematopoietic TFs. We confirmed the efficiency of these TALEs at the single-cell level using high-throughput RT-qPCR, which also allowed us to assess the consequences of both PU.1 activation and repression on wider TF networks during developmental haematopoiesis. Combined with comprehensive cellular assays, these experiments uncovered novel roles for PU.1 during early haematopoietic specification. Finally, transgenic mouse studies confirmed that the PU.1-14kb element is active at sites of definitive haematopoiesis in vivo and PU.1 is detectable in haemogenic endothelium and early committing blood cells. We therefore establish TALEs as powerful new tools to study the functionality of transcriptional networks that control developmental processes such as early haematopoiesis.

Spontaneous reactivation of clusters of X-linked genes is associated with the plasticity of X-inactivation in mouse trophoblast stem cells

Random epigenetic silencing of the X-chromosome in somatic tissues of female mammals equalizes the dosage of X-linked genes between the sexes. Unlike this form of X-inactivation that is essentially irreversible, the imprinted inactivation of the paternal X, which characterizes mouse extra-embryonic tissues, appears highly unstable in the trophoblast giant cells of the placenta. Here, we wished to determine whether such instability is already present in placental progenitor cells prior to differentiation toward lineage-specific cell types. To this end, we analyzed the behavior of a GFP transgene on the paternal X both in vivo and in trophoblast stem (TS) cells derived from the trophectoderm of XX(GFP) blastocysts. Using single-cell studies, we show that not only the GFP transgene but also a large number of endogenous genes on the paternal X are subject to orchestrated cycles of reactivation/de novo inactivation in placental progenitor cells. This reversal of silencing is associated with local losses of histone H3 lysine 27 trimethylation extending over several adjacent genes and with the topological relocation of the hypomethylated loci outside of the nuclear compartment of the inactive X. The "reactivated" state is maintained through several cell divisions. Our study suggests that this type of "metastable epigenetic" states may underlie the plasticity of TS cells and predispose specific genes to relaxed regulation in specific subtypes of placental cells.

Temporal dynamics and transcriptional control using single-cell gene expression analysis

Background

Changes in environmental conditions lead to expression variation that manifest at the level of gene regulatory networks. Despite a strong understanding of the role noise plays in synthetic biological systems, it remains unclear how propagation of expression heterogeneity in an endogenous regulatory network is distributed and utilized by cells transitioning through a key developmental event.

Results

Here we investigate the temporal dynamics of a single-cell transcriptional network of 45 transcription factors in THP-1 human myeloid monocytic leukemia cells undergoing differentiation to macrophages. We systematically measure temporal regulation of expression and variation by profiling 120 single cells at eight distinct time points, and infer highly controlled regulatory modules through which signaling operates with stochastic effects. This reveals dynamic and specific rewiring as a cellular strategy for differentiation. The integration of both positive and negative co-expression networks further identifies the proto-oncogene MYB as a network hinge to modulate both the pro- and anti-differentiation pathways.

Conclusions

Compared to averaged cell populations, temporal single-cell expression profiling provides a much more powerful technique to probe for mechanistic insights underlying cellular differentiation. We believe that our approach will form the basis of novel strategies to study the regulation of transcription at a single-cell level.

A close look at the mammalian blastocyst: epiblast and primitive endoderm formation

During early development, the mammalian embryo undergoes a series of profound changes that lead to the formation of two extraembryonic tissues--the trophectoderm and the primitive endoderm. These tissues encapsulate the pluripotent epiblast at the time of implantation. The current model proposes that the formation of these lineages results from two consecutive binary cell fate decisions. The first controls the formation of the trophectoderm and the inner cell mass, and the second controls the formation of the primitive endoderm and the epiblast within the inner cell mass. While early mammalian embryos develop with extensive plasticity, the embryonic pattern prior to implantation is remarkably reproducible. Here, we review the molecular mechanisms driving the cell fate decision between primitive endoderm and epiblast in the mouse embryo and integrate data from recent studies into the current model of the molecular network regulating the segregation between these lineages and their subsequent differentiation.

Cell fate regulation during preimplantation development: a view of adhesion-linked molecular interactions

In the developmental process of the early mammalian embryo, it is crucial to understand how the identical cells in the early embryo later develop different fates. Along with existing models, many recently discovered molecular, cellular and developmental factors play roles in cell position, cell polarity and transcriptional networks in cell fate regulation during preimplantation. A structuring process known as compaction provides the "start signal" for cells to differentiate and orchestrates the developmental cascade. The proper intercellular junctional complexes assembled between blastomeres act as a conducting mechanism governing cellular diversification. Here, we provide an overview of the diversification process during preimplantation development as it relates to intercellular junctional complexes. We also evaluate transcriptional differences between embryonic lineages according to cell- cell adhesion and the contributions of adhesion to lineage commitment. These series of processes indicate that proper cell fate specification in the early mammalian embryo depends on junctional interactions and communication, which play essential roles during early morphogenesis.

Induction of a human pluripotent state with distinct regulatory circuitry that resembles preimplantation epiblast

Human embryonic stem cells (hESCs) are derived from the inner cell mass of the blastocyst. Despite sharing the common property of pluripotency, hESCs are notably distinct from epiblast cells of the preimplantation blastocyst. Here we use a combination of three small-molecule inhibitors to sustain hESCs in a LIF signaling-dependent hESC state (3iL hESCs) with elevated expression of NANOG and epiblast-enriched genes such as KLF4, DPPA3, and TBX3. Genome-wide transcriptome analysis confirms that the expression signature of 3iL hESCs shares similarities with native preimplantation epiblast cells. We also show that 3iL hESCs have a distinct epigenetic landscape, characterized by derepression of preimplantation epiblast genes. Using genome-wide binding profiles of NANOG and OCT4, we identify enhancers that contribute to rewiring of the regulatory circuitry. In summary, our study identifies a distinct hESC state with defined regulatory circuitry that will facilitate future analysis of human preimplantation embryogenesis and pluripotency.

Cell fate inclination within 2-cell and 4-cell mouse embryos revealed by single-cell RNA sequencing

It remains an open question when and how the first cell fate decision is made in mammals. Using deep single-cell RNA-seq of matched sister blastomeres, we report highly reproducible inter-blastomere differences among 10 2-cell and five 4-cell mouse embryos. Inter-blastomere gene expression differences dominated between-embryo differences and noise, and were sufficient to cluster sister blastomeres into distinct groups. Dozens of protein-coding genes exhibited reproducible bimodal expression in sister blastomeres, which cannot be explained by random fluctuations. The protein expression of one gene out of four of these bimodal genes tested, Gadd45a, exhibited clear inter-blastomeric contrasts. We traced some of the bimodal mRNA expressions to embryonic genome activation, and others to blastomere-specific RNA depletion. Inter-blastomere differences created coexpression gene networks that were much stronger and larger than those that can possibly be created by random noise. The highly correlated gene pairs at the 4-cell stage overlapped with those showing the same directions of differential expression between inner cell mass (ICM) and trophectoderm (TE). These data substantiate the hypothesis of inter-blastomere differences in 2- and 4-cell mouse embryos, and associate these differences with ICM/TE differences.

Low-coverage single-cell mRNA sequencing reveals cellular heterogeneity and activated signaling pathways in developing cerebral cortex

Large-scale surveys of single-cell gene expression have the potential to reveal rare cell populations and lineage relationships, but require efficient methods for cell capture and mRNA sequencing^1–4. Although cellular barcoding strategies allow parallel sequencing of single cells at ultra-low depths⁵, the limitations of shallow sequencing have not been directly investigated. By capturing 301 single cells from 11 populations using microfluidics and analyzing single-cell transcriptomes across downsampled sequencing depths, we demonstrate that shallow single-cell mRNA sequencing (~50,000 reads per cell) is sufficient for unbiased cell-type classification and biomarker identification. In developing cortex we identify diverse cell types including multiple progenitor and neuronal subtypes, and we identify EGR1 and FOS as previously unreported candidate targets of Notch signaling in human but not mouse radial glia. Our strategy establishes an efficient method for unbiased analysis and comparison of cell populations from heterogeneous tissue by microfluidic single-cell capture and low-coverage sequencing of many cells.

Lytic gene expression is frequent in HSV-1 latent infection and correlates with the engagement of a cell-intrinsic transcriptional response

Herpes simplex viruses (HSV) are significant human pathogens that provide one of the best-described examples of viral latency and reactivation. HSV latency occurs in sensory neurons, being characterized by the absence of virus replication and only fragmentary evidence of protein production. In mouse models, HSV latency is especially stable but the detection of some lytic gene transcription and the ongoing presence of activated immune cells in latent ganglia have been used to suggest that this state is not entirely quiescent. Alternatively, these findings can be interpreted as signs of a low, but constant level of abortive reactivation punctuating otherwise silent latency. Using single cell analysis of transcription in mouse dorsal root ganglia, we reveal that HSV-1 latency is highly dynamic in the majority of neurons. Specifically, transcription from areas of the HSV genome associated with at least one viral lytic gene occurs in nearly two thirds of latently-infected neurons and more than half of these have RNA from more than one lytic gene locus. Further, bioinformatics analyses of host transcription showed that progressive appearance of these lytic transcripts correlated with alterations in expression of cellular genes. These data show for the first time that transcription consistent with lytic gene expression is a frequent event, taking place in the majority of HSV latently-infected neurons. Furthermore, this transcription is of biological significance in that it influences host gene expression. We suggest that the maintenance of HSV latency involves an active host response to frequent viral activity.

Primary herpes simplex virus (HSV) infections are characterized by acute disease that resolves rapidly, but the virus persists in a latent form in sensory neurons that can be a source of renewed disease. Analyzing gene expression in single mouse neurons harboring latent HSV, we show directly that HSV latency is dynamic and heterogeneous. HSV lytic gene transcripts were frequently detected in latently infected neurons and often in combinations. Expression of selected cellular anti-viral and survival genes showed that transcriptional profiles differed between latently infected and uninfected neurons from the same ganglia. The pattern of host gene expression also differed between latently infected neurons that were and were not experiencing HSV lytic gene expression. Our study suggests that HSV latency is characterized by very frequent switching on of lytic genes and a rapid response by the host, presumably to halt progression to reactivation.

Single-cell analysis reveals functionally distinct classes within the planarian stem cell compartment

Planarians are flatworms capable of regenerating any missing body region. This capacity is mediated by neoblasts, a proliferative cell population that contains pluripotent stem cells. Although population-based studies have revealed many neoblast characteristics, whether functionally distinct classes exist within this population is unclear. Here, we used high-dimensional single-cell transcriptional profiling from over a thousand individual neoblasts to directly compare gene expression fingerprints during homeostasis and regeneration. We identified two prominent neoblast classes that we named ζ (zeta) and σ (sigma). Zeta-neoblasts encompass specified cells that give rise to an abundant postmitotic lineage, including epidermal cells, and are not required for regeneration. By contrast, sigma-neoblasts proliferate in response to injury, possess broad lineage capacity, and can give rise to zeta-neoblasts. These findings indicate that planarian neoblasts comprise two major and functionally distinct cellular compartments.

Increased MAP kinase inhibition enhances epiblast-specific gene expression in bovine blastocysts

Mammalian blastocysts comprise three distinct lineages, namely, trophoblast, hypoblast, and epiblast, which develop into fetal placenta, extraembryonic yolk sac, and embryo proper, respectively. Pluripotent embryonic stem cells, capable of forming all adult cell types, can only be derived from the epiblast. In mouse and rat, this process is promoted by the double inhibition (2i) of mitogen-activated protein kinase kinase (MAP2K), which antagonizes FGF signaling, and glycogen synthase kinase 3 (GSK3), which stimulates the WNT pathway. We investigated variations of the 2i treatment on lineage segregation and pluripotency-related gene expression in bovine blastocysts. In vitro-fertilized embryos were cultured either in the presence of inhibitors of GSK3 (3 μM CHIR) and MAP2K (0.4 vs. 10 μM PD0325901, designated 2i and 2i+, respectively) or in 2i/2i+ with FGFR inhibitor (0.1 μM PD173074, designated 3i [2i and PD173074] and 3i+ [2i+ and PD173074]). Compared with 2i, both 2i+ and 3i+ potentiated the improvement in blastocyst morphology. Using an automated platform for multiplexed digital mRNA profiling, we simultaneously counted transcripts of 76 candidate genes in bovine blastocysts treated with multiple kinase inhibitors. We show that 2i+ medium specifically increased FGF4 and NANOG while reducing PDGFRalpha and SOX17 levels. The shift from a hypoblast to an epiblast gene expression signature was confirmed by quantitative PCR. A wide range of functionally related genes, including candidates involved in DNA methylation, were not significantly changed. This well-defined 2i+ effect was not observed after pharmacologically inhibiting FGF receptor or related MAP kinases (p38, JNK, and ERK5). In summary, our data suggest that increased MAP2K inhibition exerts its pluripotency-promoting effects through as yet unidentified signals.

Single-cell analysis of the transcriptome and its application in the characterization of stem cells and early embryos

Cellular heterogeneity within a cell population is a common phenomenon in multicellular organisms, tissues, cultured cells, and even FACS-sorted subpopulations. Important information may be masked if the cells are studied as a mass. Transcriptome profiling is a parameter that has been intensively studied, and relatively easier to address than protein composition. To understand the basis and importance of heterogeneity and stochastic aspects of the cell function and its mechanisms, it is essential to examine transcriptomes of a panel of single cells. High-throughput technologies, starting from microarrays and now RNA-seq, provide a full view of the expression of transcriptomes but are limited by the amount of RNA for analysis. Recently, several new approaches for amplification and sequencing the transcriptome of single cells or a limited low number of cells have been developed and applied. In this review, we summarize these major strategies, such as PCR-based methods, IVT-based methods, phi29-DNA polymerase-based methods, and several other methods, including their principles, characteristics, advantages, and limitations, with representative applications in cancer stem cells, early development, and embryonic stem cells. The prospects for development of future technology and application of transcriptome analysis in a single cell are also discussed.

Lineage-specific regulation of imprinted X inactivation in extraembryonic endoderm stem cells

Background

Silencing of the paternal X chromosome (Xp), a phenomenon known as imprinted X-chromosome inactivation (I-XCI), characterises, amongst mouse extraembryonic lineages, the primitive endoderm and the extraembryonic endoderm (XEN) stem cells derived from it.

Results

Using a combination of chromatin immunoprecipitation characterisation of histone modifications and single-cell expression studies, we show that whilst the Xp in XEN cells, like the inactive X chromosome in other cell types, globally accumulates the repressive histone mark H3K27me3, a large number of Xp genes locally lack H3K27me3 and escape from I-XCI. In most cases this escape is specific to the XEN cell lineage. Importantly, the degree of escape and the genes concerned remain unchanged upon XEN conversion into visceral endoderm, suggesting stringent control of I-XCI in XEN derivatives. Surprisingly, chemical inhibition of EZH2, a member of the Polycomb repressive complex 2 (PRC2), and subsequent loss of H3K27me3 on the Xp, do not drastically perturb the pattern of silencing of Xp genes in XEN cells.

Conclusions

The observations that we report here suggest that the maintenance of gene expression profiles of the inactive Xp in XEN cells involves a tissue-specific mechanism that acts partly independently of PRC2 catalytic activity.

Transcription factor heterogeneity in pluripotent stem cells: a stochastic advantage

When pluripotent cells are exposed to a uniform culture environment they routinely display heterogeneous gene expression. Aspects of this heterogeneity, such as Nanog expression, are linked to differences in the propensity of individual cells to either self-renew or commit towards differentiation. Recent findings have provided new insight into the underlying causes of this heterogeneity, which we summarise here using Nanog, a key regulator of pluripotency, as a model gene. We discuss the role of transcription factor heterogeneity in facilitating the intrinsically dynamic and stochastic nature of the pluripotency network, which in turn provides a potential benefit to a population of cells that needs to balance cell fate decisions.

The ability of inner-cell-mass cells to self-renew as embryonic stem cells is acquired following epiblast specification

The precise relationship of embryonic stem cells (ESCs) to cells in the mouse embryo remains controversial. We present transcriptional and functional data to identify the embryonic counterpart of ESCs. Marker profiling shows that ESCs are distinct from early inner cell mass (ICM) and closely resemble pre-implantation epiblast. A characteristic feature of mouse ESCs is propagation without ERK signalling. Single-cell culture reveals that cell-autonomous capacity to thrive when the ERK pathway is inhibited arises late during blastocyst development and is lost after implantation. The frequency of deriving clonal ESC lines suggests that all E4.5 epiblast cells can become ESCs. We further show that ICM cells from early blastocysts can progress to ERK independence if provided with a specific laminin substrate. These findings suggest that formation of the epiblast coincides with competence for ERK-independent self-renewal in vitro and consequent propagation as ESC lines.