Defining an essential transcription factor program for naïve pluripotency

Overexpression of Nuclear Receptor 5A1 Induces and Maintains an Intermediate State of Conversion between Primed and Naive Pluripotency

Naive and primed human pluripotent stem cells (hPSCs) have provided useful insights into the regulation of pluripotency. However, the molecular mechanisms regulating naive conversion remain elusive. Here, we report intermediate naive conversion induced by overexpressing nuclear receptor 5A1 (NR5A1) in hPSCs. The cells displayed some naive features, such as clonogenicity, glycogen synthase kinase 3β, and mitogen-activated protein kinase (MAPK) independence, expression of naive-associated genes, and two activated X chromosomes, but lacked others, such as KLF17 expression, transforming growth factor β independence, and imprinted gene demethylation. Notably, NR5A1 negated MAPK activation by fibroblast growth factor 2, leading to cell-autonomous self-renewal independent of MAPK inhibition. These phenotypes may be associated with naive conversion, and were regulated by a DPPA2/4-dependent pathway that activates the selective expression of naive-associated genes. This study increases our understanding of the mechanisms regulating the conversion from primed to naive pluripotency.

Landscape inferred from gene expression data governs pluripotency in embryonic stem cells

Embryonic stem cells (ESCs) can differentiate into diverse cell types and have the ability of self-renewal. Therefore, the study of cell fate decisions on embryonic stem cells has far-reaching significance for regenerative medicine and other biomedical fields. Mathematical models have been used to study emryonic stem cell differentiation. However, the underlying mechanisms of cell differentiation and lineage reprogramming remain to be elucidated. Especially, how to integrate the computational models with quantitative experimental data is still challenging. In this work, we developed a data-constrained modelling approach, and established a model of mouse embryonic stem cells. We used the truncated moment equations (TME) method to quantify the potential landscape of the ESC network. We identified two attractors on the landscape, which represent the embryonic stem cell (ESC) state and differentiated cell (DC) state, respectively, and quantified high dimensional biological paths for differentiation and reprogramming process. Through identifying the optimal combinations of gene targets based on a landscape control strategy, we offered some predictions about the key regulatory factors that govern the differentiation and reprogramming in ESCs.

Synthetic and genomic regulatory elements reveal aspects of cis-regulatory grammar in mouse embryonic stem cells

In embryonic stem cells (ESCs), a core transcription factor (TF) network establishes the gene expression program necessary for pluripotency. To address how interactions between four key TFs contribute to cis-regulation in mouse ESCs, we assayed two massively parallel reporter assay (MPRA) libraries composed of binding sites for SOX2, POU5F1 (OCT4), KLF4, and ESRRB. Comparisons between synthetic cis-regulatory elements and genomic sequences with comparable binding site configurations revealed some aspects of a regulatory grammar. The expression of synthetic elements is influenced by both the number and arrangement of binding sites. This grammar plays only a small role for genomic sequences, as the relative activities of genomic sequences are best explained by the predicted occupancy of binding sites, regardless of binding site identity and positioning. Our results suggest that the effects of transcription factor binding sites (TFBS) are influenced by the order and orientation of sites, but that in the genome the overall occupancy of TFs is the primary determinant of activity.

Transcription factors are proteins that flip genetic switches; their role is to control when and where genes are active. They do this by binding to short stretches of DNA called cis-regulatory sequences. Each sequence can have several binding sites for different transcription factors, but it is largely unclear whether the transcription factors binding to the same regulatory sequence actually work together.

It is possible that each transcription factor may work independently and there only needs to be critical mass of transcription factors bound to throw the genetic switch. If this is the case, the most important features of a cis-regulatory sequence should be the number of binding sites it contains, and how tightly the transcription factors bind to those sites. The more transcription factors and the more strongly they bind, the more active the gene should be. An alternative option is that certain transcription factors may work better together, enhancing each other's effects such that the total effect is more than the sum of its parts. If this is true, the order, orientation and spacing of the binding sites within a sequence should matter more than the number.

One way to investigate to distinguish between these possibilities is to study mouse embryonic stem cells, which have a core set of four transcription factors. Looking directly at a real genome, however, can be confusing and it is difficult to measure the effects of different cis-regulatory sequences because genes differ in so many other ways. To tackle this problem, King et al. created a synthetic set of cis-regulatory sequences based on the four core transcription factors found in mouse stem cells.

The synthetic set had every combination of two, three or four of the binding sites, with each site either facing forwards or backwards along the DNA strand. King et al. attached each of the synthetic cis-regulatory sequences to a reporter gene to find out how well each sequence performed. This revealed that the cis-regulatory sequences with the most binding sites and the tightest binding affinities work best, suggesting that transcription factors mainly work independently.

There was evidence of some interaction between some transcription factors, because, of the synthetic sequences with four binding sites, some worked better than others, and there were patterns in the most effective binding site combinations. However, these effects were small and when King et al. went on to test sequences from the real mouse genome, the most important factor by far was the number of binding sites.

Synthetic libraries of DNA sequences allow researchers to examine gene regulation more clearly than is possible in real genomes. Yet this approach does have its limitations and it is impossible to capture every type of cis-regulatory sequence in one library. The next step to extend this work is to combine the two approaches, taking sequences from the real genome and manipulating them one by one. This could help to unravel the rules that govern how cis-regulatory sequences work in real cells.

Automated inference of gene regulatory networks using explicit regulatory modules

Gene regulatory networks are a popular tool for modelling important biological phenomena, such as cell differentiation or oncogenesis. Efficient identification of the causal connections between genes, their products and regulating transcription factors, is key to understanding how defects in their function may trigger diseases. Modelling approaches should keep up with the ever more detailed descriptions of the biological phenomena at play, as provided by new experimental findings and technical improvements. In recent years, we have seen great improvements in mapping of specific binding sites of many transcription factors to distinct regulatory regions. Recent gene regulatory network models use binding measurements; but usually only to define gene-to-gene interactions, ignoring regulatory module structure. Moreover, current huge amount of transcriptomic data, and exploration of all possible cis-regulatory arrangements which can lead to the same transcriptomic response, makes manual model building both tedious and time-consuming. In our paper, we propose a method to specify possible regulatory connections in a given Boolean network, based on transcription factor binding evidence. This is implemented by an algorithm which expands a regular Boolean network model into a "cis-regulatory" Boolean network model. This expanded model explicitly defines regulatory regions as additional nodes in the network, and adds new, valuable biological insights to the system dynamics. The expanded model can automatically be compared with expression data. And, for each node, a regulatory function, consistent with the experimental data, can be found. The resulting models are usually more constrained (by biologically-motivated metadata), and can then be inspected in in silico simulations. The fully automated method for model identification has been implemented in Python, and the expansion algorithm in R. The method resorts to the Z3 Satisfiability Modulo Theories (SMT) solver, and is similar to the RE:IN application (Yordanov et al., 2016). It is available on https://github.com/regulomics/expansion-network.

Nicotinamide adenine dinucleotide induces a bivalent metabolism and maintains pluripotency in human embryonic stem cells

Nicotinamide adenine dinucleotide (NAD⁺ ) and its precursor metabolites are emerging as important regulators of both cell metabolism and cell state. Interestingly, the role of NAD⁺ in human embryonic stem cell (hESC) metabolism and the regulation of pluripotent cell state is unresolved. Here we show that NAD⁺ simultaneously increases hESC mitochondrial oxidative metabolism and partially suppresses glycolysis and stimulates amino acid turnover, doubling the consumption of glutamine. Concurrent with this metabolic remodeling, NAD⁺ increases hESC pluripotent marker expression and proliferation, inhibits BMP4-induced differentiation and reduces global histone 3 lysine 27 trimethylation, plausibly inducing an intermediate naïve-to-primed bivalent metabolism and pluripotent state. Furthermore, maintenance of NAD⁺ recycling via malate aspartate shuttle activity is identified as an absolute requirement for hESC self-renewal, responsible for 80% of the oxidative capacity of hESC mitochondria. Our findings implicate NAD⁺ in the regulation of cell state, suggesting that the hESC pluripotent state is dependent upon cellular NAD⁺ .

scMAGeCK links genotypes with multiple phenotypes in single-cell CRISPR screens

We present scMAGeCK, a computational framework to identify genomic elements associated with multiple expression-based phenotypes in CRISPR/Cas9 functional screening that uses single-cell RNA-seq as readout. scMAGeCK outperforms existing methods, identifies genes and enhancers with known and novel functions in cell proliferation, and enables an unbiased construction of genotype-phenotype network. Single-cell CRISPR screening on mouse embryonic stem cells identifies key genes associated with different pluripotency states. Applying scMAGeCK on multiple datasets, we identify key factors that improve the power of single-cell CRISPR screening. Collectively, scMAGeCK is a novel tool to study genotype-phenotype relationships at a single-cell level.

Phosphorylation of TFCP2L1 by CDK1 is required for stem cell pluripotency and bladder carcinogenesis

Molecular programs involved in embryogenesis are frequently upregulated in oncogenic dedifferentiation and metastasis. However, their precise roles and regulatory mechanisms remain elusive. Here, we showed that CDK1 phosphorylation of TFCP2L1, a pluripotency-associated transcription factor, orchestrated pluripotency and cell-cycling in embryonic stem cells (ESCs) and was aberrantly activated in aggressive bladder cancers (BCs). In murine ESCs, the protein interactome and transcription targets of Tfcp2l1 indicated its involvement in cell cycle regulation. Tfcp2l1 was phosphorylated at Thr177 by Cdk1, which affected ESC cell cycle progression, pluripotency, and differentiation. The CDK1-TFCP2L1 pathway was activated in human BC cells, stimulating their proliferation, self-renewal, and invasion. Lack of TFCP2L1 phosphorylation impaired the tumorigenic potency of BC cells in a xenograft model. In patients with BC, high co-expression of TFCP2L1 and CDK1 was associated with unfavorable clinical characteristics including tumor grade, lymphovascular and muscularis propria invasion, and distant metastasis and was an independent prognostic factor for cancer-specific survival. These findings demonstrate the molecular and clinical significance of CDK1-mediated TFCP2L1 phosphorylation in stem cell pluripotency and in the tumorigenic stemness features associated with BC progression.

Machine learning applications in drug development

Due to the huge amount of biological and medical data available today, along with well-established machine learning algorithms, the design of largely automated drug development pipelines can now be envisioned. These pipelines may guide, or speed up, drug discovery; provide a better understanding of diseases and associated biological phenomena; help planning preclinical wet-lab experiments, and even future clinical trials. This automation of the drug development process might be key to the current issue of low productivity rate that pharmaceutical companies currently face. In this survey, we will particularly focus on two classes of methods: sequential learning and recommender systems, which are active biomedical fields of research.

PRDM1 controls the sequential activation of neural, neural crest and sensory progenitor determinants

During early embryogenesis, the ectoderm is rapidly subdivided into neural, neural crest and sensory progenitors. How the onset of lineage determinants and the loss of pluripotency markers are temporally and spatially coordinated in vivo is still debated. Here, we identify a crucial role for the transcription factor PRDM1 in the orderly transition from epiblast to defined neural lineages in chick. PRDM1 is initially expressed broadly in the entire epiblast, but becomes gradually restricted as cell fates are specified. We find that PRDM1 is required for the loss of some pluripotency markers and the onset of neural, neural crest and sensory progenitor specifier genes. PRDM1 directly activates their expression by binding to their promoter regions and recruiting the histone demethylase Kdm4a to remove repressive histone marks. However, once neural lineage determinants become expressed, they in turn repress PRDM1, whereas prolonged PRDM1 expression inhibits neural, neural crest and sensory progenitor genes, suggesting that its downregulation is necessary for cells to maintain their identity. Therefore, PRDM1 plays multiple roles during ectodermal cell fate allocation.

Dynamic regulation of chromatin accessibility by pluripotency transcription factors across the cell cycle

The pioneer activity of transcription factors allows for opening of inaccessible regulatory elements and has been extensively studied in the context of cellular differentiation and reprogramming. In contrast, the function of pioneer activity in self-renewing cell divisions and across the cell cycle is poorly understood. Here we assessed the interplay between OCT4 and SOX2 in controlling chromatin accessibility of mouse embryonic stem cells. We found that OCT4 and SOX2 operate in a largely independent manner even at co-occupied sites, and that their cooperative binding is mostly mediated indirectly through regulation of chromatin accessibility. Controlled protein degradation strategies revealed that the uninterrupted presence of OCT4 is required for post-mitotic re-establishment and interphase maintenance of chromatin accessibility, and that highly OCT4-bound enhancers are particularly vulnerable to transient loss of OCT4 expression. Our study sheds light on the constant pioneer activity required to maintain the dynamic pluripotency regulatory landscape in an accessible state.

Recent insights into the naïve state of human pluripotency and its applications

The past decade has seen significant interest in the isolation of pluripotent stem cells corresponding to various stages of mammalian embryonic development. Two distinct and well-defined pluripotent states can be derived from mouse embryos: "naïve" pluripotent cells with properties of pre-implantation epiblast, and "primed" pluripotent cells, resembling post-implantation epiblast. Prompted by the successful interconversion between these two stem cell states in the mouse system, several groups have devised strategies for inducing a naïve state of pluripotency in human pluripotent stem cells. Here, we review recent insights into the naïve state of human pluripotency, focusing on two methods that confer defining transcriptomic and epigenomic signatures of the pre-implantation embryo. The isolation of naïve human pluripotent stem cells offers a window into early developmental mechanisms that cannot be adequately modeled in primed cells, such as X chromosome reactivation, metabolic reprogramming, and the regulation of hominid-specific transposable elements. We outline key unresolved questions regarding naïve human pluripotency, including its extrinsic and intrinsic control mechanisms, potential for embryonic and extraembryonic differentiation, and general utility as a model system for human development and disease.

A gene regulatory network controls the balance between mesendoderm and ectoderm at pluripotency exit

During embryogenesis, differentiation of pluripotent cells into somatic cell types depends both on signaling cues and intrinsic gene expression programs. While the molecular underpinnings of pluripotency are well mapped, much less is known on how mouse embryonic stem cells (mESCs) differentiate. Using RNA-Seq profiling during specification to the three germ layers, we showed that mESCs switched on condition-specific gene expression programs from the onset of the differentiation procedure and that primed pluripotency did not constitute an obligatory intermediate state. After inferring the gene network controlling mESC differentiation, we tested the role of the highly connected nodes by deleting them in a triple knock-in Sox1-Brachyury-Eomes mESC line reporting on ectoderm, mesoderm, and endoderm fates. This led to the identification of regulators of mESC differentiation that acted at several levels: Sp1 as a global break on differentiation, Nr5a2 controlling ectoderm specification, and notably Fos:Jun and Zfp354c as opposite switches between ectoderm and mesendoderm fate.

Regulation of Cell Fate Decisions in Early Mammalian Embryos

Early embryogenesis is characterized by the segregation of cell lineages that fulfill critical roles in the establishment of pregnancy and development of the fetus. The formation of the blastocyst marks the emergence of extraembryonic precursors, needed for implantation, and of pluripotent cells, which differentiate toward the major lineages of the adult organism. The coordinated emergence of these cell types shows that these processes are broadly conserved in mammals. However, developmental heterochrony and changes in gene regulatory networks highlight unique evolutionary adaptations that may explain the diversity in placentation and in the mechanisms controlling pluripotency in mammals. The incorporation of new technologies, including single-cell omics, imaging, and gene editing, is instrumental for comparative embryology. Broadening the knowledge of mammalian embryology will provide new insights into the mechanisms driving evolution and development. This knowledge can be readily translated into biomedical and biotechnological applications in humans and livestock, respectively.

The developmental origins of sex-biased expression in cardiac development

BACKGROUND

Expression patterns between males and females vary in every adult tissue, even in organs with no conspicuous dimorphisms such as the heart. While studies of male and female differences have traditionally focused on the influence of sex hormones, these do not account for all the differences at the molecular and epigenetic levels. We previously reported that a substantial number of genes were differentially expressed in male and female mouse embryonic stem (ES) cells and revealed dose-dependent enhancer activity in response to Prdm14, a key pluripotency factor expressed more highly in female ES cells. In this work, we investigated the role of Prdm14 in establishing sex-specific gene expression networks. We surveyed the sex-specific landscape in early embryogenesis with special reference to cardiac development. We generated sex-specific co-expression networks from mouse ES cells, examined the presence of sex-specific chromatin domains, and analyzed previously published datasets from different developmental time points to characterize how sex-biased gene expression waxes and wanes to evaluate whether sex-biased networks are detectable throughout heart development.

RESULTS

We performed ChIP-seq on male and female mouse ES cells to determine differences in chromatin status. Our study reveals sex-biased histone modifications, underscoring the potential for the sex chromosome complement to prime the genome differently in early development with consequences for later expression biases. Upon differentiation of ES cells to cardiac precursors, we found sex-biased expression of key transcription and epigenetic factors, some of which persisted from the undifferentiated state. Using network analyses, we also found that Prdm14 plays a prominent role in regulating a subset of dimorphic expression patterns. To determine whether sex-biased expression is present throughout cardiogenesis, we re-analyzed data from two published studies that sampled the transcriptomes of mouse hearts from 8.5 days post-coitum embryos to neonates and adults. We found sex-biased expression at every stage in heart development, and interestingly, identified a subset of genes that exhibit the same bias across multiple cardiogenic stages.

CONCLUSIONS

Overall, our results support the existence of sexually dimorphic gene expression profiles and regulatory networks at every stage of cardiac development, some of which may be established in early embryogenesis and epigenetically perpetuated.

Deriving rabbit embryonic stem cells by small molecule inhibitors

We previously developed pluripotent rabbit embryonic stem cells (rbES) using a culture system supplemented with basic fibroblast growth factor (bFGF) and leukemia inhibitory factor (LIF), noggin and Y-27632 (referred to as iFLY). In present work, we explored multiple approaches to enhance the chance of deriving domed pluripotent rbES cells by inhibition of MEK, GSK, and PKC signaling pathways. Domed stated rbES were derived in defined medium supplemented with 15% KOSR, 103 IU/mL mouse LIF, 10 ng/mL bFGF and three inhibitors to the MEK (PD0325901, 1 µM), GSK3 (CHIR99021, 3 µM) and PKC (Gö6983, 5 µM) (3i). Domed rbES were passaged every 3-4 days till passage 3-4 for the designated experiments. We showed that bFGF and LIF are indispensable for the derivation and maintenance of rbES; whereas the 3i medium containing inhibitors to the MEK (PD0325901), GSK3 (CHIR99021) and PKC (Gö6983) were necessary for deriving domed rbES. Domed rbES possessed naïve ES markers as Rex1 and ERAS in addition to Oct4, Klf4, Sox 2 and c-myc by RT-PCR. Domed rbES showed positive staining for Rex1, Fgf4, Klf4, Nanog and Oct4 by immunofluorescence chemistry. Further deleting either one factor in 3i medium as CHIR99021, PD0325901, Gö6983 or bFGF resulted in disappearing of domed rbES colonies. The optimal concentrations of 3i contained 0.75 µM PD0325901, 2.25 µM CHIR99021, and 4.5 µM Gö6983. Our work, in combination of different inhibitors for deriving rabbit ES, supports that the network of signal pathways plays an important role in ES self-renew, propagation and maintenance, and sheds light on deriving authentic properties of rbES in an important yet understudied model animal species.

Stem cell bioengineering: building from stem cell biology

New fundamental discoveries in stem cell biology have yielded potentially transformative regenerative therapeutics. However, widespread implementation of stem-cell-derived therapeutics remains sporadic. Barriers that impede the development of these therapeutics can be linked to our incomplete understanding of how the regulatory networks that encode stem cell fate govern the development of the complex tissues and organs that are ultimately required for restorative function. Bioengineering tools, strategies and design principles represent core components of the stem cell bioengineering toolbox. Applied to the different layers of complexity present in stem-cell-derived systems - from gene regulatory networks in single stem cells to the systemic interactions of stem-cell-derived organs and tissues - stem cell bioengineering can address existing challenges and advance regenerative medicine and cellular therapies.

Stem Cell Differentiation as a Non-Markov Stochastic Process

Pluripotent stem cells can self-renew in culture and differentiate along all somatic lineages in vivo. While much is known about the molecular basis of pluripotency, the mechanisms of differentiation remain unclear. Here, we profile individual mouse embryonic stem cells as they progress along the neuronal lineage. We observe that cells pass from the pluripotent state to the neuronal state via an intermediate epiblast-like state. However, analysis of the rate at which cells enter and exit these observed cell states using a hidden Markov model indicates the presence of a chain of unobserved molecular states that each cell transits through stochastically in sequence. This chain of hidden states allows individual cells to record their position on the differentiation trajectory, thereby encoding a simple form of cellular memory. We suggest a statistical mechanics interpretation of these results that distinguishes between functionally distinct cellular “macrostates” and functionally similar molecular “microstates” and propose a model of stem cell differentiation as a non-Markov stochastic process.

•

We profile individual stem cells as they differentiate along the neural lineage

•

Regulatory network changes and increased cell variability accompany differentiation

•

Analysis of dynamics with a hidden Markov model reveals unobserved molecular states

•

We propose a model of stem cell differentiation as a non-Markov stochastic process

Glutamine independence is a selectable feature of pluripotent stem cells

Most rapidly proliferating mammalian cells rely on the oxidation of exogenous glutamine to support cell proliferation. We previously found that culture of mouse embryonic stem cells (ESCs) in the presence of inhibitors against MEK and GSK3β to maintain pluripotency reduces cellular reliance on glutamine for tricarboxylic acid (TCA) cycle anaplerosis, enabling ESCs to proliferate in the absence of exogenous glutamine. Here we show that reduced dependence on exogenous glutamine is a generalizable feature of pluripotent stem cells. Enhancing self-renewal, through either overexpression of pluripotency-associated transcription factors or altered signal transduction, decreases the utilization of glutamine-derived carbons in the TCA cycle. As a result, cells with the highest potential for self-renewal can be enriched by transient culture in glutamine-deficient media. During pluripotent cell culture or reprogramming to pluripotency, transient glutamine withdrawal selectively leads to the elimination of non-pluripotent cells. These data reveal that reduced dependence on glutamine anaplerosis is an inherent feature of self-renewing pluripotent stem cells and reveal a simple, non-invasive mechanism to select for mouse and human pluripotent stem cells within a heterogeneous population during both ESC passage and induced pluripotent cell reprogramming.

Epiblast Formation by TEAD-YAP-Dependent Expression of Pluripotency Factors and Competitive Elimination of Unspecified Cells

The epiblast is a pluripotent cell population first formed in preimplantation embryos, and its quality is important for proper development. Here, we examined the mechanisms of epiblast formation and found that the Hippo pathway transcription factor TEAD and its coactivator YAP regulate expression of pluripotency factors. After specification of the inner cell mass, YAP accumulates in the nuclei and activates TEAD. TEAD activity is required for strong expression of pluripotency factors and is variable in the forming epiblast. Cells showing low TEAD activity are eliminated from the epiblast through cell competition. Pluripotency factor expression and MYC control cell competition downstream of TEAD activity. Cell competition eliminates unspecified cells and is required for proper organization of the epiblast. These results suggest that induction of pluripotency factors by TEAD activity and elimination of unspecified cells via cell competition ensure the production of an epiblast with naive pluripotency.

Codon usage optimization in pluripotent embryonic stem cells

BACKGROUND

The uneven use of synonymous codons in the transcriptome regulates the efficiency and fidelity of protein translation rates. Yet, the importance of this codon bias in regulating cell state-specific expression programmes is currently debated. Here, we ask whether different codon usage controls gene expression programmes in self-renewing and differentiating embryonic stem cells.

RESULTS

Using ribosome and transcriptome profiling, we identify distinct codon signatures during human embryonic stem cell differentiation. We find that cell state-specific codon bias is determined by the guanine-cytosine (GC) content of differentially expressed genes. By measuring the codon frequencies at the ribosome active sites interacting with transfer RNAs (tRNA), we further discover that self-renewing cells optimize translation of codons that depend on the inosine tRNA modification in the anticodon wobble position. Accordingly, inosine levels are highest in human pluripotent embryonic stem cells. This effect is conserved in mice and is independent of the differentiation stimulus.

CONCLUSIONS

We show that GC content influences cell state-specific mRNA levels, and we reveal how translational mechanisms based on tRNA modifications change codon usage in embryonic stem cells.

Complementary Activity of ETV5, RBPJ, and TCF3 Drives Formative Transition from Naive Pluripotency

The gene regulatory network (GRN) of naive mouse embryonic stem cells (ESCs) must be reconfigured to enable lineage commitment. TCF3 sanctions rewiring by suppressing components of the ESC transcription factor circuitry. However, TCF3 depletion only delays and does not prevent transition to formative pluripotency. Here, we delineate additional contributions of the ETS-family transcription factor ETV5 and the repressor RBPJ. In response to ERK signaling, ETV5 switches activity from supporting self-renewal and undergoes genome relocation linked to commissioning of enhancers activated in formative epiblast. Independent upregulation of RBPJ prevents re-expression of potent naive factors, TBX3 and NANOG, to secure exit from the naive state. Triple deletion of Etv5, Rbpj, and Tcf3 disables ESCs, such that they remain largely undifferentiated and locked in self-renewal, even in the presence of differentiation stimuli. Thus, genetic elimination of three complementary drivers of network transition stalls developmental progression, emulating environmental insulation by small-molecule inhibitors.

The gene regulatory network of mESC differentiation: a benchmark for reverse engineering methods

A large body of data have accumulated that characterize the gene regulatory network of stem cells. Yet, a comprehensive and integrative understanding of this complex network is lacking. Network reverse engineering methods that use transcriptome data to derive these networks may help to uncover the topology in an unbiased way. Many methods exist that use co-expression to reconstruct networks. However, it remains unclear how these methods perform in the context of stem cell differentiation, as most systematic assessments have been made for regulatory networks of unicellular organisms. Here, we report a systematic benchmark of different reverse engineering methods against functional data. We show that network pruning is critical for reconstruction performance. We also find that performance is similar for algorithms that use different co-expression measures, i.e. mutual information or correlation. In addition, different methods yield very different network topologies, highlighting the challenge of interpreting these resulting networks as a whole.This article is part of the theme issue 'Designer human tissue: coming to a lab near you'.

Epigenetic modulation of a hardwired 3D chromatin landscape in two naive states of pluripotency

The mechanisms underlying enhancer activation and the extent to which enhancer-promoter rewiring contributes to spatiotemporal gene expression are not well understood. Using integrative and time-resolved analyses we show that the extensive transcriptome and epigenome resetting during the conversion between 'serum' and '2i' states of mouse embryonic stem cells (ESCs) takes place with minimal enhancer-promoter rewiring that becomes more evident in primed-state pluripotency. Instead, differential gene expression is strongly linked to enhancer activation via H3K27ac. Conditional depletion of transcription factors and allele-specific enhancer analysis reveal an essential role for Esrrb in H3K27 acetylation and activation of 2i-specific enhancers. Restoration of a polymorphic ESRRB motif using CRISPR-Cas9 in a hybrid ESC line restores ESRRB binding and enhancer H3K27ac in an allele-specific manner but has no effect on chromatin interactions. Our study shows that enhancer activation in serum- and 2i-ESCs is largely driven by transcription factor binding and epigenetic marking in a hardwired network of chromatin interactions.

Defined conditions for propagation and manipulation of mouse embryonic stem cells

The power of mouse embryonic stem (ES) cells to colonise the developing embryo has revolutionised mammalian developmental genetics and stem cell research. This power is vulnerable, however, to the cell culture environment, deficiencies in which can lead to cellular heterogeneity, adaptive phenotypes, epigenetic aberrations and genetic abnormalities. Here, we provide detailed methodologies for derivation, propagation, genetic modification and primary differentiation of ES cells in 2i or 2i+LIF media without serum or undefined serum substitutes. Implemented diligently, these procedures minimise variability and deviation, thereby improving the efficiency, reproducibility and biological validity of ES cell experimentation.

A TRIM71 binding long noncoding RNA Trincr1 represses FGF/ERK signaling in embryonic stem cells

Long noncoding RNAs (lncRNAs) have emerged as important components of gene regulatory network in embryonic stem cells (ESCs). However, the function and molecular mechanism of lncRNAs are still largely unknown. Here we identifies Trincr1 (TRIM71 interacting long noncoding RNA 1) lncRNA that regulates the FGF/ERK signaling and self-renewal of ESCs. Trincr1 is exported by THOC complex to cytoplasm where it binds and represses TRIM71, leading to the downregulation of SHCBP1 protein. Knocking out Trincr1 leads to the upregulation of phosphorylated ERK and ERK pathway target genes and the decrease of ESC self-renewal, while knocking down Trim71 completely rescues the defects of Trincr1 knockout. Furthermore, ectopic expression of Trincr1 represses FGF/ERK signaling and the self-renewal of neural progenitor cells (NPCs). Together, this study highlights lncRNA as an important player in cell signaling network to coordinate cell fate specification.

FGF signaling through ERK is known to promote the differentiation of embryonic stem cells (ES cells). Here, the authors demonstrate that the lncRNA Trincr1 binds and represses TRIM71 in ES cells, leading to downregulation of SHCBP1 protein, the reduction of FGF/ERK signaling and the promotion of self-renewal.

Dynamic m6A mRNA methylation reveals the role of METTL3-m6A-CDCP1 signaling axis in chemical carcinogenesis

N6-methyladenosine (m⁶A) is the most abundant internal modification in mammalian mRNAs. Despite its functional importance in various physiological events, the role of m⁶A in chemical carcinogenesis remains largely unknown. Here we profiled the dynamic m⁶A mRNA modification during cellular transformation induced by chemical carcinogens and identified a subset of cell transformation-related, concordantly modulated m⁶A sites. Notably, the increased m⁶A in 3'-UTR mRNA of oncogene CDCP1 was found in malignant transformed cells. Mechanistically, the m⁶A methyltransferase METTL3 and demethylases ALKBH5 mediate the m⁶A modification in 3'-UTR of CDCP1 mRNA. METTL3 and m⁶A reader YTHDF1 preferentially recognize m⁶A residues on CPCP1 3'-UTR and promote CDCP1 translation. We further showed that METTL3 and CDCP1 are upregulated in the bladder cancer patient samples and the expression of METTL3 and CDCP1 is correlated with the progression status of the bladder cancers. Inhibition of the METTL3-m⁶A-CDCP1 axis resulted in decreased growth and progression of chemical-transformed cells and bladder cancer cells. Most importantly, METTL3-m⁶A-CDCP1 axis has synergistic effect with chemical carcinogens in promoting malignant transformation of uroepithelial cells and bladder cancer tumorigenesis in vitro and in vivo. Taken together, our results identify dynamic m⁶A modification in chemical-induced malignant transformation and provide insight into critical roles of the METTL3-m⁶A-CDCP1 axis in chemical carcinogenesis.

The transcription factor TFCP2L1 induces expression of distinct target genes and promotes self-renewal of mouse and human embryonic stem cells

TFCP2L1 (transcription factor CP2-like 1) is a transcriptional regulator critical for maintaining mouse and human embryonic stem cell (ESC) pluripotency. However, the direct TFCP2L1 target genes are uncharacterized. Here, using gene overexpression, immunoblotting, quantitative real-time PCR, ChIP, and reporter gene assays, we show that TFCP2L1 primarily induces estrogen-related receptor β (Esrrb) expression that supports mouse ESC identity and also selectively enhances Kruppel-like factor 4 (Klf4) expression and thereby promotes human ESC self-renewal. Specifically, we found that in mouse ESCs, TFCP2L1 binds directly to the Esrrb gene promoter and regulates its transcription. Esrrb knockdown impaired Tfcp2l1's ability to induce interleukin 6 family cytokine (leukemia inhibitory factor)-independent ESC self-renewal and to reprogram epiblast stem cells to naïve pluripotency. Conversely, Esrrb overexpression blocked differentiation induced by Tfcp2l1 down-regulation. Moreover, we identified Klf4 as a direct TFCP2L1 target in human ESCs, bypassing the requirement for activin A and basic fibroblast growth factor in short-term human ESC self-renewal. Enforced Klf4 expression recapitulated the self-renewal-promoting effect of Tfcp2l1, whereas Klf4 knockdown eliminated these effects and caused loss of colony-forming capability. These findings indicate that TFCP2L1 functions differently in naïve and primed pluripotency, insights that may help elucidate the different states of pluripotency.

Computational Strategies for Niche-Dependent Cell Conversion to Assist Stem Cell Therapy

The field of regenerative medicine has blossomed in recent decades. However, the ultimate goal of tissue regeneration - replacing damaged or aged cells with healthy functioning cells - still faces a number of challenges. In particular, better understanding of the role of the cellular niche in shaping stem cell phenotype and conversion would aid in improving current protocols for stem cell therapies. In this regard, the implementation of novel computational approaches that consider the niche effect on stem cells would be valuable. Here we discuss current problems in stem cell transplantation and rejuvenation, and we propose computational strategies to control niche-dependent cell conversion to overcome them.

Folliculin regulates mTORC1/2 and WNT pathways in early human pluripotency

To reveal how cells exit human pluripotency, we designed a CRISPR-Cas9 screen exploiting the metabolic and epigenetic differences between naïve and primed pluripotent cells. We identify the tumor suppressor, Folliculin(FLCN) as a critical gene required for the exit from human pluripotency. Here we show that FLCN Knock-out (KO) hESCs maintain the naïve pluripotent state but cannot exit the state since the critical transcription factor TFE3 remains active in the nucleus. TFE3 targets up-regulated in FLCN KO exit assay are members of Wnt pathway and ESRRB. Treatment of FLCN KO hESC with a Wnt inhibitor, but not ESRRB/FLCN double mutant, rescues the cells, allowing the exit from the naïve state. Using co-immunoprecipitation and mass spectrometry analysis we identify unique FLCN binding partners. The interactions of FLCN with components of the mTOR pathway (mTORC1 and mTORC2) reveal a mechanism of FLCN function during exit from naïve pluripotency.

The pathways involved in exit from pluripotency in human embryonic stem cells are poorly understood. Here, the authors performed a CRISPR-based screen to identify genes that promote exit from naïve pluripotency and find a role for folliculin (FLCN) by regulating the mTOR and Wnt pathways.

Embryonic stem cells become mechanoresponsive upon exit from ground state of pluripotency

Stem cell fate decisions are driven by a broad array of signals, both chemical and mechanical. Although much progress has been made in our understanding of the impact of chemical signals on cell fate choice, much less is known about the role and influence of mechanical signalling, particularly in embryonic stem (ES) cells. Many studies use substrates with different stiffness to study mechanical signalling, but changing substrate stiffness can induce secondary effects which are difficult to disentangle from the direct effects of forces/mechanical signals. To probe the direct impact of mechanical stress on cells, we developed an adaptable cell substrate stretcher to exert specific, reproducible forces on cells. Using this device to test the response of ES cells to tensile strain, we found that cells experienced a transient influx of calcium followed by an upregulation of the so-called immediate and early genes. On longer time scales, however, ES cells in ground state conditions were largely insensitive to mechanical stress. Nonetheless, as ES cells exited the ground state, their susceptibility to mechanical signals increased, resulting in broad transcriptional changes. Our findings suggest that exit from ground state of pluripotency is unaffected by mechanical signals, but that these signals could become important during the next stage of lineage specification. A better understanding of this process could improve our understanding of cell fate choice in early development and improve protocols for differentiation guided by mechanical cues.

Machine Learning of Stem Cell Identities From Single-Cell Expression Data via Regulatory Network Archetypes

The molecular regulatory network underlying stem cell pluripotency has been intensively studied, and we now have a reliable ensemble model for the "average" pluripotent cell. However, evidence of significant cell-to-cell variability suggests that the activity of this network varies within individual stem cells, leading to differential processing of environmental signals and variability in cell fates. Here, we adapt a method originally designed for face recognition to infer regulatory network patterns within individual cells from single-cell expression data. Using this method we identify three distinct network configurations in cultured mouse embryonic stem cells-corresponding to naïve and formative pluripotent states and an early primitive endoderm state-and associate these configurations with particular combinations of regulatory network activity archetypes that govern different aspects of the cell's response to environmental stimuli, cell cycle status and core information processing circuitry. These results show how variability in cell identities arise naturally from alterations in underlying regulatory network dynamics and demonstrate how methods from machine learning may be used to better understand single cell biology, and the collective dynamics of cell communities.

Frontiers of Pluripotency

Humans develop from a unique group of pluripotent cells in early embryos that can produce all cells of the human body. While pluripotency is only transiently manifest in the embryo, scientists have identified conditions that sustain pluripotency indefinitely in the laboratory. Pluripotency is not a monolithic entity, however, but rather comprises a spectrum of different cellular states. Questions regarding the scientific value of examining the continuum of pluripotent stem (PS) cell states have gained increased significance in light of attempts to generate interspecies chimeras between humans and animals. In this chapter, I review our ever-evolving understanding of the continuum of pluripotency. Historically, the discovery of two different PS cell states in mice fostered a general conception of pluripotency comprised of two distinct attractor states: naïve and primed. Naïve pluripotency has been defined by competence to form germline chimeras and governance by unique KLF-based transcription factor (TF) circuitry, whereas primed state is distinguished by an inability to generate chimeras and alternative TF regulation. However, the discovery of many alternative PS cell states challenges the concept of pluripotency as a binary property. Moreover, it remains unclear whether the current molecular criteria used to classify human naïve-like pluripotency also identify human chimera-competent PS cells. Therefore, I examine the pluripotency continuum more closely in light of recent advances in PS cell research and human interspecies chimera research.

Automated Formal Reasoning to Uncover Molecular Programs of Self-Renewal

The Reasoning Engine for Interaction Networks (RE:IN) is a tool that was developed initially for the study of pluripotency in mouse embryonic stem cells. A set of critical factors that regulate the pluripotent state had been identified experimentally, but it was not known how these genes interacted to stabilize self-renewal or commit the cell to differentiation. The methodology encapsulated in RE:IN enabled the exploration of a space of possible network interaction models, allowing for uncertainty in whether individual interactions exist between the pluripotency factors. This concept of an "abstract" network was combined with automated reasoning that allows the user to eliminate models that are inconsistent with experimental observations. The tool generalizes beyond the study of stem cell decision-making, allowing for the study of interaction networks more broadly across biology.

[Fifty shades of pluripotency]

Since the derivation of the first pluripotent embryonic stem cell lines in mice in the early 1980s, a plethora of lines has been obtained from various mammalian species including rodents, lagomorphs and primates. These lines are distinguished by their molecular and functional characteristics and correspond to the different pluripotency states observed in the developing embryo between the "blastocyst" and "gastrula" stages. These cell lines are positioned along a gradient, or continuum of pluripotency, the ends of which are epitomized by the naïve and primed states, respectively. Conventional human pluripotent stem cells self-renew in the primed state of pluripotency (ie, at the bottom of the gradient), a position that is undoubtedly the cause of their natural instability. Recent studies aim to generate naive human pluripotent stem cells (at the top of the gradient). The importance of this research in the perspective of medical applications will be discussed.

Evidence of Extracellular Vesicles Biogenesis and Release in Mouse Embryonic Stem Cells

Extracellular vesicles (EVs) released by mouse embryonic stem cells (mESCs) are considered a source of bioactive molecules that modulate their microenvironment by acting on intercellular communication. Either intracellular endosomal machinery or their derived EVs have been considered a relevant system of signal circuits processing. Herein, we show that these features are found in mESCs. Ultrastructural analysis revealed structures and organelles of the endosomal system such as coated pits and endocytosis-related vesicles, prominent rough endoplasmic reticulum and Golgi apparatus, and multivesicular bodies (MVBs) containing either few or many intraluminal vesicles (ILVs) that could be released as exosomes to extracellular milieu. Besides, budding vesicles shed from the plasma membrane to the extracellular space is suggestive of microvesicle biogenesis in mESCs. mESCs and mouse blastocyst express specific markers of the Endosomal Sorting Complex Required for Transport (ESCRT) system. Ultrastructural analysis and Nanoparticle Tracking Analysis (NTA) of isolated EVs revealed a heterogeneous population of exosomes and microvesicles released by mESCs. These vesicles contain Wnt10b and the Notch ligand Delta-like 4 (DLL4) and also the co-chaperone stress inducible protein 1 (STI1) and its partner Hsp90. Wnt10b and Dll4 colocalize with EVs biogenesis markers in mESCs. Overall, the present study supports the function of the mESCs endocytic network and their EVs as players in stem cell biology.

MTCH2-mediated mitochondrial fusion drives exit from naïve pluripotency in embryonic stem cells

The role of mitochondria dynamics and its molecular regulators remains largely unknown during naïve-to-primed pluripotent cell interconversion. Here we report that mitochondrial MTCH2 is a regulator of mitochondrial fusion, essential for the naïve-to-primed interconversion of murine embryonic stem cells (ESCs). During this interconversion, wild-type ESCs elongate their mitochondria and slightly alter their glutamine utilization. In contrast, MTCH2^-/- ESCs fail to elongate their mitochondria and to alter their metabolism, maintaining high levels of histone acetylation and expression of naïve pluripotency markers. Importantly, enforced mitochondria elongation by the pro-fusion protein Mitofusin (MFN) 2 or by a dominant negative form of the pro-fission protein dynamin-related protein (DRP) 1 is sufficient to drive the exit from naïve pluripotency of both MTCH2^-/- and wild-type ESCs. Taken together, our data indicate that mitochondria elongation, governed by MTCH2, plays a critical role and constitutes an early driving force in the naïve-to-primed pluripotency interconversion of murine ESCs.

A common molecular logic determines embryonic stem cell self-renewal and reprogramming

During differentiation and reprogramming, new cell identities are generated by reconfiguration of gene regulatory networks. Here, we combined automated formal reasoning with experimentation to expose the logic of network activation during induction of naïve pluripotency. We find that a Boolean network architecture defined for maintenance of naïve state embryonic stem cells (ESC) also explains transcription factor behaviour and potency during resetting from primed pluripotency. Computationally identified gene activation trajectories were experimentally substantiated at single‐cell resolution by RT–qPCR. Contingency of factor availability explains the counterintuitive observation that Klf2, which is dispensable for ESC maintenance, is required during resetting. We tested 124 predictions formulated by the dynamic network, finding a predictive accuracy of 77.4%. Finally, we show that this network explains and predicts experimental observations of somatic cell reprogramming. We conclude that a common deterministic program of gene regulation is sufficient to govern maintenance and induction of naïve pluripotency. The tools exemplified here could be broadly applied to delineate dynamic networks underlying cell fate transitions.

Causal Transcription Regulatory Network Inference Using Enhancer Activity as a Causal Anchor

Transcription control plays a crucial role in establishing a unique gene expression signature for each of the hundreds of mammalian cell types. Though gene expression data have been widely used to infer cellular regulatory networks, existing methods mainly infer correlations rather than causality. We developed statistical models and likelihood-ratio tests to infer causal gene regulatory networks using enhancer RNA (eRNA) expression information as a causal anchor and applied the framework to eRNA and transcript expression data from the FANTOM Consortium. Predicted causal targets of transcription factors (TFs) in mouse embryonic stem cells, macrophages and erythroblastic leukaemia overlapped significantly with experimentally-validated targets from ChIP-seq and perturbation data. We further improved the model by taking into account that some TFs might act in a quantitative, dosage-dependent manner, whereas others might act predominantly in a binary on/off fashion. We predicted TF targets from concerted variation of eRNA and TF and target promoter expression levels within a single cell type, as well as across multiple cell types. Importantly, TFs with high-confidence predictions were largely different between these two analyses, demonstrating that variability within a cell type is highly relevant for target prediction of cell type-specific factors. Finally, we generated a compendium of high-confidence TF targets across diverse human cell and tissue types.

Canonical and single-cell Hi-C reveal distinct chromatin interaction sub-networks of mammalian transcription factors

Background

Transcription factor (TF) binding to regulatory DNA sites is a key determinant of cell identity within multi-cellular organisms and has been studied extensively in relation to site affinity and chromatin modifications. There has been a strong focus on the inference of TF-gene regulatory networks and TF-TF physical interaction networks. Here, we present a third type of TF network, the spatial network of co-localized TF binding sites within the three-dimensional genome.

Results

Using published canonical Hi-C data and single-cell genome structures, we assess the spatial proximity of a genome-wide array of potential TF-TF co-localizations in human and mouse cell lines. For individual TFs, the abundance of occupied binding sites shows a positive correspondence with their clustering in three dimensions, and this is especially apparent for weak TF binding sites and at enhancer regions. An analysis between different TF proteins identifies significantly proximal pairs, which are enriched in reported physical interactions. Furthermore, clustering of different TFs based on proximity enrichment identifies two partially segregated co-localization sub-networks, involving different TFs in different cell types. Using data from both human lymphoblastoid cells and mouse embryonic stem cells, we find that these sub-networks are enriched within, but not exclusive to, different chromosome sub-compartments that have been identified previously in Hi-C data.

Conclusions

This suggests that the association of TFs within spatial networks is closely coupled to gene regulatory networks. This applies to both differentiated and undifferentiated cells and is a potential causal link between lineage-specific TF binding and chromosome sub-compartment segregation.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1558-2) contains supplementary material, which is available to authorized users.

Tracing the transitions from pluripotency to germ cell fate with CRISPR screening

Early mammalian development entails transit through naive pluripotency towards post-implantation epiblast, which subsequently gives rise to primordial germ cells (PGC), the founding germline population. To investigate these cell fate transitions, we developed a compound-reporter to track cellular identity in a model of PGC specification (PGC-like cells; PGCLC), and coupled it with genome-wide CRISPR screening. We identify key genes both for exit from pluripotency and for acquisition of PGC fate, and characterise a central role for the transcription regulators Nr5a2 and Zfp296 in germline ontogeny. Abrogation of these genes results in widespread activation (Nr5a2^−/−) or inhibition (Zfp296^−/−) of WNT pathway factors in PGCLC. This leads to aberrant upregulation of the somatic programme or failure to activate germline genes, respectively, and consequently loss of germ cell identity. Our study places Zfp296 and Nr5a2 as key components of an expanded PGC gene regulatory network, and outlines a transferable strategy for identifying critical regulators of complex cell fate decisions.

Primordial Germ Cell-Like Cells (PGCLCs) are an in vitro model for primordial germ cell development. Here, the authors couple a novel compound reporter with CRISPR screening to identify key genes for exit from pluripotency and acquisition of PGCLC fate; specifically identifying Nr5a2 and Zfp296.

Synthesizing Signaling Pathways from Temporal Phosphoproteomic Data

We present a method for automatically discovering signaling pathways from time-resolved phosphoproteomic data. The Temporal Pathway Synthesizer (TPS) algorithm uses constraint-solving techniques first developed in the context of formal verification to explore paths in an interaction network. It systematically eliminates all candidate structures for a signaling pathway where a protein is activated or inactivated before its upstream regulators. The algorithm can model more than one hundred thousand dynamic phosphosites and can discover pathway members that are not differentially phosphorylated. By analyzing temporal data, TPS defines signaling cascades without needing to experimentally perturb individual proteins. It recovers known pathways and proposes pathway connections when applied to the human epidermal growth factor and yeast osmotic stress responses. Independent kinase mutant studies validate predicted substrates in the TPS osmotic stress pathway.

Quantitative subcellular proteomics using SILAC reveals enhanced metabolic buffering in the pluripotent ground state

The ground state of pluripotency is defined as a minimal unrestricted epigenetic state as present in the Inner Cell Mass. Mouse embryonic stem cells (ESCs) grown in a defined serum-free medium with two kinase inhibitors ("2i ESCs") have been postulated to reflect ground-state pluripotency, whereas ESCs grown in the presence of serum ("serum ESCs") share more similarities with post-implantation epiblast cells. Pluripotency results from an intricate interplay between cytoplasmic, nuclear and chromatin-associated proteins. Here, we perform quantitative subcellular proteomics to gain insight in the molecular mechanisms sustaining the pluripotent states reflected by 2i and serum ESCs. We describe a full SILAC workflow and quality controls for proteomic comparison of 2i and serum ESCs, allowing subcellular proteomics of the cytoplasm, nucleoplasm and chromatin. The obtained quantitative information revealed increased levels of naïve pluripotency factors on the chromatin of 2i ESCs. Surprisingly, the cytoplasmic proteome suggests that 2i and serum ESCs utilize distinct metabolic programs, which include upregulation of free radical buffering by the glutathione pathway in 2i ESCs. Through induction of intracellular radicals, we show that the altered metabolic environment renders 2i ESCs less sensitive to oxidative stress. Altogether, this work provides novel insights into the proteomic landscape underlying ground state pluripotency.

LIF-dependent survival of embryonic stem cells is regulated by a novel palmitoylated Gab1 signalling protein

The cytokine leukaemia inhibitory factor (LIF) promotes self-renewal of mouse embryonic stem cells (ESCs) through activation of the transcription factor Stat3. However, the contribution of other ancillary pathways stimulated by LIF in ESCs, such as the MAPK and PI3K pathways, is less well understood. We show here that naive-type mouse ESCs express high levels of a novel effector of the MAPK and PI3K pathways. This effector is an isoform of the Gab1 (Grb2-associated binder protein 1) adaptor protein that lacks the N-terminal pleckstrin homology (PH) membrane-binding domain. Although not essential for rapid unrestricted growth of ESCs under optimal conditions, the novel Gab1 variant (Gab1β) is required for LIF-mediated cell survival under conditions of limited nutrient availability. This enhanced survival is absolutely dependent upon a latent palmitoylation site that targets Gab1β directly to ESC membranes. These results show that constitutive association of Gab1 with membranes through a novel mechanism promotes LIF-dependent survival of murine ESCs in nutrient-poor conditions.

Transition state characteristics during cell differentiation

Models describing the process of stem-cell differentiation are plentiful, and may offer insights into the underlying mechanisms and experimentally observed behaviour. Waddington's epigenetic landscape has been providing a conceptual framework for differentiation processes since its inception. It also allows, however, for detailed mathematical and quantitative analyses, as the landscape can, at least in principle, be related to mathematical models of dynamical systems. Here we focus on a set of dynamical systems features that are intimately linked to cell differentiation, by considering exemplar dynamical models that capture important aspects of stem cell differentiation dynamics. These models allow us to map the paths that cells take through gene expression space as they move from one fate to another, e.g. from a stem-cell to a more specialized cell type. Our analysis highlights the role of the transition state (TS) that separates distinct cell fates, and how the nature of the TS changes as the underlying landscape changes-change that can be induced by e.g. cellular signaling. We demonstrate that models for stem cell differentiation may be interpreted in terms of either a static or transitory landscape. For the static case the TS represents a particular transcriptional profile that all cells approach during differentiation. Alternatively, the TS may refer to the commonly observed period of heterogeneity as cells undergo stochastic transitions.

Single-cell RNA-sequencing reveals the existence of naive and primed pluripotency in pre-implantation rhesus monkey embryos

Naive pluripotency exists in epiblast cells of mouse pre-implantation embryos. However, whether the naive pluripotency is transient or nonexistent in primate embryos remains unclear. Using RNA-seq in single blastomeres from 16-cell embryos through to hatched blastocysts of rhesus monkey, we constructed the lineage segregation roadmap in which the specification of trophectoderm, epiblast, and primitive endoderm is initiated simultaneously at the early blastocyst stage. Importantly, we uncovered the existence of distinct pluripotent states in monkey pre-implantation embryos. At the early- and middle-blastocyst stages, the epiblast cells have the transcriptome features of naive pluripotency, whereas they display a continuum of primed pluripotency characteristics at the late and hatched blastocyst stages. Moreover, we identified potential regulators that might play roles in the transition from naive to primed pluripotency. Thus, our study suggests the transient existence of naive pluripotency in primates and proposes an ideal time window for derivation of primate embryonic stem cells with naive pluripotency.

Negative feedback via RSK modulates Erk-dependent progression from naïve pluripotency

Mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase (ERK) signalling is implicated in initiation of embryonic stem (ES) cell differentiation. The pathway is subject to complex feedback regulation. Here, we examined the ERK-responsive phosphoproteome in ES cells and identified the negative regulator RSK1 as a prominent target. We used CRISPR/Cas9 to create combinatorial mutations in RSK family genes. Genotypes that included homozygous null mutations in Rps6ka1, encoding RSK1, resulted in elevated ERK phosphorylation. These RSK-depleted ES cells exhibit altered kinetics of transition into differentiation, with accelerated downregulation of naïve pluripotency factors, precocious expression of transitional epiblast markers and early onset of lineage specification. We further show that chemical inhibition of RSK increases ERK phosphorylation and expedites ES cell transition without compromising multilineage potential. These findings demonstrate that the ERK activation profile influences the dynamics of pluripotency progression and highlight the role of signalling feedback in temporal control of cell state transitions.