Defining an essential transcription factor program for naïve pluripotency

The regulatory architecture of the primed pluripotent cell state

Despite extensive research, the gene regulatory architecture governing mammalian cell states remains poorly understood. Here we present an integrative systems biology approach to elucidate the network architecture of primed state pluripotency. Using an unbiased methodology, we identified and experimentally confirmed 132 transcription factors as master regulators (MRs) of mouse epiblast stem cell (EpiSC) pluripotency, many of which were further validated by CRISPR-mediated functional assays. To assemble a comprehensive regulatory network, we silenced each of the 132 MRs to assess their effects on the other MRs and their transcriptional targets, yielding a network of 1273 MR → MR interactions. Network architecture analyses revealed four functionally distinct MR modules (communities), and identified key Speaker and Mediator MRs based on their hierarchical rank and centrality. Our findings elucidate the de-centralized logic of a "communal interaction" model in which the balanced activities of four MR communities maintain primed state pluripotency.

Optimized reporters for multiplexed detection of transcription factor activity

In any given cell type, dozens of transcription factors (TFs) act in concert to control the activity of the genome by binding to specific DNA sequences in regulatory elements. Despite their considerable importance, we currently lack simple tools to directly measure the activity of many TFs in parallel. Massively parallel reporter assays (MPRAs) allow the detection of TF activities in a multiplexed fashion; however, we lack basic understanding to rationally design sensitive reporters for many TFs. Here, we use an MPRA to systematically optimize transcriptional reporters for 86 TFs and evaluate the specificity of all reporters across a wide array of TF perturbation conditions. We thus identified critical TF reporter design features and obtained highly sensitive and specific reporters for 62 TFs, many of which outperform available reporters. The resulting collection of "prime" TF reporters can be used to uncover TF regulatory networks and to illuminate signaling pathways. A record of this paper's transparent peer review process is included in the supplemental information.

Propagating pluripotency - The conundrum of self-renewal

The discovery of mouse embryonic stem cells in 1981 transformed research in mammalian developmental biology and functional genomics. The subsequent generation of human pluripotent stem cells (PSCs) and the development of molecular reprogramming have opened unheralded avenues for drug discovery and cell replacement therapy. Here, I review the history of PSCs from the perspective that long-term self-renewal is a product of the in vitro signaling environment, rather than an intrinsic feature of embryos. I discuss the relationship between pluripotent states captured in vitro to stages of epiblast in the embryo and suggest key considerations for evaluation of PSCs. A remaining fundamental challenge is to determine whether naïve pluripotency can be propagated from the broad range of mammals by exploiting common principles in gene regulatory architecture.

OneSC: a computational platform for recapitulating cell state transitions

MOTIVATION

Computational modeling of cell state transitions has been a great interest of many in the field of developmental biology, cancer biology, and cell fate engineering because it enables performing perturbation experiments in silico more rapidly and cheaply than could be achieved in a lab. Recent advancements in single-cell RNA-sequencing (scRNA-seq) allow the capture of high-resolution snapshots of cell states as they transition along temporal trajectories. Using these high-throughput datasets, we can train computational models to generate in silico "synthetic" cells that faithfully mimic the temporal trajectories.

RESULTS

Here we present OneSC, a platform that can simulate cell state transitions using systems of stochastic differential equations govern by a regulatory network of core transcription factors (TFs). Different from many current network inference methods, OneSC prioritizes on generating Boolean network that produces faithful cell state transitions and terminal cell states that mimic real biological systems. Applying OneSC to real data, we inferred a core TF network using a mouse myeloid progenitor scRNA-seq dataset and showed that the dynamical simulations of that network generate synthetic single-cell expression profiles that faithfully recapitulate the four myeloid differentiation trajectories going into differentiated cell states (erythrocytes, megakaryocytes, granulocytes, and monocytes). Finally, through the in silico perturbations of the mouse myeloid progenitor core network, we showed that OneSC can accurately predict cell fate decision biases of TF perturbations that closely match with previous experimental observations.

AVAILABILITY AND IMPLEMENTATION

OneSC is implemented as a Python package on GitHub (https://github.com/CahanLab/oneSC) and on Zenodo (https://zenodo.org/records/14052421).

Approaches for Benchmarking Single-Cell Gene Regulatory Network Methods

Gene regulatory networks are powerful tools for modeling genetic interactions that control the expression of genes driving cell differentiation, and single-cell sequencing offers a unique opportunity to build these networks with high-resolution genomic data. There are many proposed computational methods to build these networks using single-cell data, and different approaches are used to benchmark these methods. However, a comprehensive discussion specifically focusing on benchmarking approaches is missing. In this article, we lay the GRN terminology, present an overview of common gold-standard studies and data sets, and define the performance metrics for benchmarking network construction methodologies. We also point out the advantages and limitations of different benchmarking approaches, suggest alternative ground truth data sets that can be used for benchmarking, and specify additional considerations in this context.

Large-scale analysis of the integration of enhancer-enhancer signals by promoters

Genes are often regulated by multiple enhancers. It is poorly understood how the individual enhancer activities are combined to control promoter activity. Anecdotal evidence has shown that enhancers can combine sub-additively, additively, synergistically, or redundantly. However, it is not clear which of these modes are more frequent in mammalian genomes. Here, we systematically tested how pairs of enhancers activate promoters using a three-way combinatorial reporter assay in mouse embryonic stem cells. By assaying about 69,000 enhancer-enhancer-promoter combinations we found that enhancer pairs generally combine near-additively. This behaviour was conserved across seven developmental promoters tested. Surprisingly, these promoters scale the enhancer signals in a non-linear manner that depends on promoter strength. A housekeeping promoter showed an overall different response to enhancer pairs, and a smaller dynamic range. Thus, our data indicate that enhancers mostly act additively, but promoters transform their collective effect non-linearly.

GPR160 regulates the self-renewal and pluripotency of mouse embryonic stem cells via JAK1/STAT3 signaling pathway

G-protein-coupled receptors (GPCRs) are the largest family of transmembrane receptors and regulate various physiological and pathological processes. Despite extensive studies, the roles of GPCRs in mouse embryonic stem cells (mESCs) remain poorly understood. Here, we show that GPR160, a class A member of GPCRs, is dramatically downregulated concurrent with mESC differentiation into embryoid bodies in vitro. Knockdown of Gpr160 leads to downregulation of the expression of pluripotency-associated transcription factors and upregulation of the expression of lineage markers, accompanying with the arrest of the mESC cell-cycle in the G0/G1 phase. RNA-seq analysis shows that GPR160 participates in the JAK/STAT signaling pathway crucial for maintaining ESC stemness, and the knockdown of Gpr160 results in the downregulation of STAT3 phosphorylation level, which in turn is partially rescued by colivelin, a STAT3 activator. Consistent with these observations, GPR160 physically interacts with JAK1, and cooperates with leukemia inhibitory factor receptor (LIFR) and gp130 to activate the STAT3 pathway. In summary, our results suggest that GPR160 regulates mESC self-renewal and pluripotency by interacting with the JAK1-LIFR-gp130 complex to mediate the JAK1/STAT3 signaling pathway.

FoxO transcription factors actuate the formative pluripotency specific gene expression programme

Naïve pluripotency is sustained by a self-reinforcing gene regulatory network (GRN) comprising core and naïve pluripotency-specific transcription factors (TFs). Upon exiting naïve pluripotency, embryonic stem cells (ESCs) transition through a formative post-implantation-like pluripotent state, where they acquire competence for lineage choice. However, the mechanisms underlying disengagement from the naïve GRN and initiation of the formative GRN are unclear. Here, we demonstrate that phosphorylated AKT acts as a gatekeeper that prevents nuclear localisation of FoxO TFs in naïve ESCs. PTEN-mediated reduction of AKT activity upon exit from naïve pluripotency allows nuclear entry of FoxO TFs, enforcing a cell fate transition by binding and activating formative pluripotency-specific enhancers. Indeed, FoxO TFs are necessary and sufficient for the activation of the formative pluripotency-specific GRN. Our work uncovers a pivotal role for FoxO TFs in establishing formative post-implantation pluripotency, a critical early embryonic cell fate transition.

Dissecting SOX2 expression and function reveals an association with multiple signaling pathways during embryonic development and in cancer progression

SRY (Sex Determining Region) box 2 (SOX2) is an essential transcription factor that plays crucial roles in activating genes involved in pre- and post-embryonic development, adult tissue homeostasis, and lineage specifications. SOX2 maintains the self-renewal property of stem cells and is involved in the generation of induced pluripotency stem cells. SOX2 protein contains a particular high-mobility group domain that enables SOX2 to achieve the capacity to participate in a broad variety of functions. The information about the involvement of SOX2 with gene regulatory elements, signaling networks, and microRNA is gradually emerging, and the higher expression of SOX2 is functionally relevant to various cancer types. SOX2 facilitates the oncogenic phenotype via cellular proliferation and enhancement of invasive tumor properties. Evidence are accumulating in favor of three dimensional (higher order) folding of chromatin and epigenetic control of the SOX2 gene by chromatin modifications, which implies that the expression level of SOX2 can be modulated by epigenetic regulatory mechanisms, specifically, via DNA methylation and histone H3 modification. In view of this, and to focus further insights into the roles SOX2 plays in physiological functions, involvement of SOX2 during development, precisely, the advances of our knowledge in pre- and post-embryonic development, and interactions of SOX2 in this scenario with various signaling pathways in tumor development and cancer progression, its potential as a therapeutic target against many cancers are summarized and discussed in this article.

Establishment of porcine embryonic stem cells in simplified serum free media and feeder free expansion

BACKGROUND

The establishment of stable porcine embryonic stem cells (pESCs) can contribute to basic and biomedical research, including comparative developmental biology, as well as assessing the safety of stem cell-based therapies. Despite these advantages, most pESCs obtained from in vitro blastocysts require complex media and feeder layers, making routine use, genetic modification, and differentiation into specific cell types difficult. We aimed to establish pESCs with a single cell-passage ability, high proliferative potency, and stable in long-term culture from in vitro-derived blastocysts using a simplified serum-free medium.

METHODS

We evaluated the establishment efficiency of pESCs from in vitro blastocysts using various basal media (DMEM/F10 (1:1), DMEM/F12, and a-MEM) and factors (FGF2, IWR-1, CHIR99021, and WH-4-023). The pluripotency and self-renewal capacity of the established pESCs were analyzed under feeder or feeder-free conditions. Ultimately, we developed a simplified culture medium (FIW) composed of FGF2, IWR-1, and WH-4-023 under serum-free conditions.

RESULTS

The pESC-FIW lines were capable of single-cell passaging with short cell doubling times and expressed the pluripotency markers POU5F1, SOX2, and NANOG, as well as cell surface markers SSEA1, SSEA4, and TRA-1-60. pESC-FIW showed a stable proliferation rate and normal karyotype, even after 50 passages. Transcriptome analysis revealed that pESC-FIW were similar to reported pESC maintained in complex media and showed gastrulating epiblast cell characteristics. pESC-FIW were maintained for multiple passages under feeder-free conditions on fibronectin-coated plates using mTeSR™, a commercial medium used for feeder-free culture, exhibiting characteristics similar to those observed under feeder conditions.

CONCLUSIONS

These results indicated that inhibition of WNT and SRC was sufficient to establish pESCs capable of single-cell passaging and feeder-free expansion under serum-free conditions. The easy maintenance of pESCs facilitates their application in gene editing technology for agriculture and biomedicine, as well as lineage commitment studies.

Optimized reporters for multiplexed detection of transcription factor activity

In any given cell type, dozens of transcription factors (TFs) act in concert to control the activity of the genome by binding to specific DNA sequences in regulatory elements. Despite their considerable importance in determining cell identity and their pivotal role in numerous disorders, we currently lack simple tools to directly measure the activity of many TFs in parallel. Massively parallel reporter assays (MPRAs) allow the detection of TF activities in a multiplexed fashion; however, we lack basic understanding to rationally design sensitive reporters for many TFs. Here, we use an MPRA to systematically optimize transcriptional reporters for 86 TFs and evaluate the specificity of all reporters across a wide array of TF perturbation conditions. We thus identified critical TF reporter design features and obtained highly sensitive and specific reporters for 60 TFs, many of which outperform available reporters. The resulting collection of "prime" TF reporters can be used to uncover TF regulatory networks and to illuminate signaling pathways.

ERK signalling eliminates Nanog and maintains Oct4 to drive the formative pluripotency transition

Naïve epiblast cells in the embryo and pluripotent stem cells in vitro undergo developmental progression to a formative state competent for lineage specification. During this transition, transcription factors and chromatin are rewired to encode new functional features. Here, we examine the role of mitogen-activated protein kinase (ERK1/2) signalling in pluripotent state transition. We show that a primary consequence of ERK activation in mouse embryonic stem cells is elimination of Nanog, which precipitates breakdown of the naïve state gene regulatory network. Variability in pERK dynamics results in heterogeneous loss of Nanog and metachronous state transition. Knockdown of Nanog allows exit without ERK activation. However, transition to formative pluripotency does not proceed and cells collapse to an indeterminate identity. This outcome is due to failure to maintain expression of the central pluripotency factor Oct4. Thus, during formative transition ERK signalling both dismantles the naïve state and preserves pluripotency. These results illustrate how a single signalling pathway can both initiate and secure transition between cell states.

Mechanisms of Embryonic Stem Cell Pluripotency Maintenance and Their Application in Livestock and Poultry Breeding

Embryonic stem cells (ESCs) are remarkably undifferentiated cells that originate from the inner cell mass of the blastocyst. They possess the ability to self-renew and differentiate into multiple cell types, making them invaluable in diverse applications such as disease modeling and the creation of transgenic animals. In recent years, as agricultural practices have evolved from traditional to biological breeding, it has become clear that pluripotent stem cells (PSCs), either ESCs or induced pluripotent stem cells (iPSCs), are optimal for continually screening suitable cellular materials. However, the technologies for long-term in vitro culture or establishment of cell lines for PSCs in livestock are still immature, and research progress is uneven, which poses challenges for the application of PSCs in various fields. The establishment of a robust in vitro system for these cells is critically dependent on understanding their pluripotency maintenance mechanisms. It is believed that the combined effects of pluripotent transcription factors, pivotal signaling pathways, and epigenetic regulation contribute to maintaining their pluripotent state, forming a comprehensive regulatory network. This article will delve into the primary mechanisms underlying the maintenance of pluripotency in PSCs and elaborate on the applications of PSCs in the field of livestock.

OneSC: A computational platform for recapitulating cell state transitions

Computational modelling of cell state transitions has been a great interest of many in the field of developmental biology, cancer biology and cell fate engineering because it enables performing perturbation experiments in silico more rapidly and cheaply than could be achieved in a wet lab. Recent advancements in single-cell RNA sequencing (scRNA-seq) allow the capture of high-resolution snapshots of cell states as they transition along temporal trajectories. Using these high-throughput datasets, we can train computational models to generate in silico 'synthetic' cells that faithfully mimic the temporal trajectories. Here we present OneSC, a platform that can simulate synthetic cells across developmental trajectories using systems of stochastic differential equations govern by a core transcription factors (TFs) regulatory network. Different from the current network inference methods, OneSC prioritizes on generating Boolean network that produces faithful cell state transitions and steady cell states that mimic real biological systems. Applying OneSC to real data, we inferred a core TF network using a mouse myeloid progenitor scRNA-seq dataset and showed that the dynamical simulations of that network generate synthetic single-cell expression profiles that faithfully recapitulate the four myeloid differentiation trajectories going into differentiated cell states (erythrocytes, megakaryocytes, granulocytes and monocytes). Finally, through the in-silico perturbations of the mouse myeloid progenitor core network, we showed that OneSC can accurately predict cell fate decision biases of TF perturbations that closely match with previous experimental observations.

A gray box framework that optimizes a white box logical model using a black box optimizer for simulating cellular responses to perturbations

Predicting cellular responses to perturbations requires interpretable insights into molecular regulatory dynamics to perform reliable cell fate control, despite the confounding non-linearity of the underlying interactions. There is a growing interest in developing machine learning-based perturbation response prediction models to handle the non-linearity of perturbation data, but their interpretation in terms of molecular regulatory dynamics remains a challenge. Alternatively, for meaningful biological interpretation, logical network models such as Boolean networks are widely used in systems biology to represent intracellular molecular regulation. However, determining the appropriate regulatory logic of large-scale networks remains an obstacle due to the high-dimensional and discontinuous search space. To tackle these challenges, we present a scalable derivative-free optimizer trained by meta-reinforcement learning for Boolean network models. The logical network model optimized by the trained optimizer successfully predicts anti-cancer drug responses of cancer cell lines, while simultaneously providing insight into their underlying molecular regulatory mechanisms.

Enhancer-promoter interactions are reconfigured through the formation of long-range multiway hubs as mouse ES cells exit pluripotency

Enhancers bind transcription factors, chromatin regulators, and non-coding transcripts to modulate the expression of target genes. Here, we report 3D genome structures of single mouse ES cells as they are induced to exit pluripotency and transition through a formative stage prior to undergoing neuroectodermal differentiation. We find that there is a remarkable reorganization of 3D genome structure where inter-chromosomal intermingling increases dramatically in the formative state. This intermingling is associated with the formation of a large number of multiway hubs that bring together enhancers and promoters with similar chromatin states from typically 5-8 distant chromosomal sites that are often separated by many Mb from each other. In the formative state, genes important for pluripotency exit establish contacts with emerging enhancers within these multiway hubs, suggesting that the structural changes we have observed may play an important role in modulating transcription and establishing new cell identities.

Beyond pathogens: the intriguing genetic legacy of endogenous retroviruses in host physiology

The notion that viruses played a crucial role in the evolution of life is not a new concept. However, more recent insights suggest that this perception might be even more expansive, highlighting the ongoing impact of viruses on host evolution. Endogenous retroviruses (ERVs) are considered genomic remnants of ancient viral infections acquired throughout vertebrate evolution. Their exogenous counterparts once infected the host's germline cells, eventually leading to the permanent endogenization of their respective proviruses. The success of ERV colonization is evident so that it constitutes 8% of the human genome. Emerging genomic studies indicate that endogenous retroviruses are not merely remnants of past infections but rather play a corollary role, despite not fully understood, in host genetic regulation. This review presents some evidence supporting the crucial role of endogenous retroviruses in regulating host genetics. We explore the involvement of human ERVs (HERVs) in key physiological processes, from their precise and orchestrated activities during cellular differentiation and pluripotency to their contributions to aging and cellular senescence. Additionally, we discuss the costs associated with hosting a substantial amount of preserved viral genetic material.

Mitotic bookmarking redundancy by nuclear receptors in pluripotent cells

Mitotic bookmarking transcription factors (TFs) are thought to mediate rapid and accurate reactivation after mitotic gene silencing. However, the loss of individual bookmarking TFs often leads to the deregulation of only a small proportion of their mitotic targets, raising doubts on the biological significance and importance of their bookmarking function. Here we used targeted proteomics of the mitotic bookmarking TF ESRRB, an orphan nuclear receptor, to discover a large redundancy in mitotic binding among members of the protein super-family of nuclear receptors. Focusing on the nuclear receptor NR5A2, which together with ESRRB is essential in maintaining pluripotency in mouse embryonic stem cells, we demonstrate conjoint bookmarking activity of both factors on promoters and enhancers of a large fraction of active genes, particularly those most efficiently reactivated in G1. Upon fast and simultaneous degradation of both factors during mitotic exit, hundreds of mitotic targets of ESRRB/NR5A2, including key players of the pluripotency network, display attenuated transcriptional reactivation. We propose that redundancy in mitotic bookmarking TFs, especially nuclear receptors, confers robustness to the reestablishment of gene regulatory networks after mitosis.

The Nucleosome Remodelling and Deacetylation complex coordinates the transcriptional response to lineage commitment in pluripotent cells

As cells exit the pluripotent state and begin to commit to a specific lineage they must activate genes appropriate for that lineage while silencing genes associated with pluripotency and preventing activation of lineage-inappropriate genes. The Nucleosome Remodelling and Deacetylation (NuRD) complex is essential for pluripotent cells to successfully undergo lineage commitment. NuRD controls nucleosome density at regulatory sequences to facilitate transcriptional responses, and also has been shown to prevent unscheduled transcription (transcriptional noise) in undifferentiated pluripotent cells. How these activities combine to ensure cells engage a gene expression program suitable for successful lineage commitment has not been determined. Here, we show that NuRD is not required to silence all genes. Rather, it restricts expression of genes primed for activation upon exit from the pluripotent state, but maintains them in a transcriptionally permissive state in self-renewing conditions, which facilitates their subsequent activation upon exit from naïve pluripotency. We further show that NuRD coordinates gene expression changes, which acts to maintain a barrier between different stable states. Thus NuRD-mediated chromatin remodelling serves multiple functions, including reducing transcriptional noise, priming genes for activation and coordinating the transcriptional response to facilitate lineage commitment.

Comparative transcriptomic analysis of Illumina and MGI next-generation sequencing platforms using RUNX3- and ZBTB46-instructed embryonic stem cells

Introduction: We have previously observed phenotypic and developmental changes upon the ectopic expression of the RUNX3 or the ZBTB46 transcription factors in mouse embryonic stem cell (ESC) derived progenitors. In this study, we evaluated the gene expression profiles of the RUNX3- and the ZBTB46-instructed murine ESCs with RNA-seq testing two next-generation sequencing technologies. Methods: We compared the DNA nanoball-based DNBSEQ G400 sequencer (MGI) with the bridge-PCR-based NextSeq 500 instrument (Illumina) for RNA sequencing. Moreover, we also compared two types of MGI sequencing reagents (Standard versus Hot-massive parallel sequencing (MPS)) with the DNBSEQ G400. Results: We observed that both sequencing platforms showed comparable levels of quality, sequencing uniformity, and gene expression profiles. For example, highly overlapping RUNX3- and ZBTB46-regulated gene lists were obtained from both sequencing datasets. Moreover, we observed that the Standard and the Hot-MPS-derived RUNX3- and ZBTB46-regulated gene lists were also considerably overlapped. This transcriptome analysis also helped us to identify differently expressed genes in the presence of the transgenic RUNX3 or ZBTB46. For example, we found that Gzmb, Gzmd, Gzme, Gdf6, and Ccr7 genes were robustly upregulated upon the forced expression of Runx3; on the other hand, Gpx2, Tdpoz4, and Arg2 were induced alongside the ectopic expression of Zbtb46. Discussion: Similar gene expression profile and greatly overlapping RUNX3- and ZBTB46-regulated gene sets were detected with both DNA sequencing platforms. Our analyses demonstrate that both sequencing technologies are suitable for transcriptome profiling and target gene selection. These findings suggest that DNBSEQ G400 represents a cost-effective alternative sequencing platform for gene expression monitoring. Moreover, this analysis provides a resource for exploration of the RUNX3- and ZBTB46-dependent gene regulatory networks.

Human Stem Cells for Ophthalmology: Recent Advances in Diagnostic Image Analysis and Computational Modelling

Purpose of Review

To explore the advances and future research directions in image analysis and computational modelling of human stem cells (hSCs) for ophthalmological applications.

Recent Findings

hSCs hold great potential in ocular regenerative medicine due to their application in cell-based therapies and in disease modelling and drug discovery using state-of-the-art 2D and 3D organoid models. However, a deeper characterisation of their complex, multi-scale properties is required to optimise their translation to clinical practice. Image analysis combined with computational modelling is a powerful tool to explore mechanisms of hSC behaviour and aid clinical diagnosis and therapy.

Summary

Many computational models draw on a variety of techniques, often blending continuum and discrete approaches, and have been used to describe cell differentiation and self-organisation. Machine learning tools are having a significant impact in model development and improving image classification processes for clinical diagnosis and treatment and will be the focus of much future research.

Inferring single-cell gene regulatory network by non-redundant mutual information

Gene regulatory network plays a crucial role in controlling the biological processes of living creatures. Deciphering the complex gene regulatory networks from experimental data remains a major challenge in system biology. Recent advances in single-cell RNA sequencing technology bring massive high-resolution data, enabling computational inference of cell-specific gene regulatory networks (GRNs). Many relevant algorithms have been developed to achieve this goal in the past years. However, GRN inference is still less ideal due to the extra noises involved in pseudo-time information and large amounts of dropouts in datasets. Here, we present a novel GRN inference method named Normi, which is based on non-redundant mutual information. Normi manipulates these problems by employing a sliding size-fixed window approach on the entire trajectory and conducts average smoothing strategy on the gene expression of the cells in each window to obtain representative cells. To further alleviate the impact of dropouts, we utilize the mixed KSG estimator to quantify the high-order time-delayed mutual information among genes, then filter out the redundant edges by adopting Max-Relevance and Min Redundancy algorithm. Moreover, we determined the optimal time delay for each gene pair by distance correlation. Normi outperforms other state-of-the-art GRN inference methods on both simulated data and single-cell RNA sequencing (scRNA-seq) datasets, demonstrating its superiority in robustness. The performance of Normi in real scRNA-seq data further reveals its ability to identify the key regulators and crucial biological processes.

NAT10-mediated N4-acetylcytidine mRNA modification regulates self-renewal in human embryonic stem cells

NAT10-catalyzed N4-acetylcytidine (ac4C) has emerged as a vital post-transcriptional modulator on the coding transcriptome by promoting mRNA stability. However, its role in mammalian development remains unclear. Here, we found that NAT10 expression positively correlates with pluripotency in vivo and in vitro. High throughput ac4C-targeted RNA immunoprecipitation sequencing (ac4C-RIP-seq), NaCNBH3-based chemical ac4C sequencing (ac4C-seq) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) assays revealed noticeable ac4C modifications in transcriptome of hESCs, among which transcripts encoding core pluripotency transcription factors are favorable targets of ac4C modification. Further validation assays demonstrate that genetic inactivation of NAT10, the ac4C writer enzyme, led to ac4C level decrease on target genes, promoted the core pluripotency regulator OCT4 (POU5F1) transcript decay, and finally impaired self-renewal and promoted early differentiation in hESCs. Together, our work presented here elucidates a previously unrecognized interconnectivity between the core pluripotent transcriptional network for the maintenance of human ESC self-renewal and NAT10-catalyzed ac4C RNA epigenetic modification.

Inhibition of Wnt activity improves peri-implantation development of somatic cell nuclear transfer embryos

Somatic cell nuclear transfer (SCNT) can reprogram differentiated somatic cells into totipotency. Although pre-implantation development of SCNT embryos has greatly improved, most SCNT blastocysts are still arrested at the peri-implantation stage, and the underlying mechanism remains elusive. Here, we develop a 3D in vitro culture system for SCNT peri-implantation embryos and discover that persistent Wnt signals block the naïve-to-primed pluripotency transition of epiblasts with aberrant H3K27me3 occupancy, which in turn leads to defects in epiblast transformation events and subsequent implantation failure. Strikingly, manipulating Wnt signals can attenuate the pluripotency transition and H3K27me3 deposition defects in epiblasts and achieve up to a 9-fold increase in cloning efficiency. Finally, single-cell RNA-seq analysis reveals that Wnt inhibition markedly enhances the lineage developmental trajectories of SCNT blastocysts during peri-implantation development. Overall, these findings reveal diminished potentials of SCNT blastocysts for lineage specification and validate a critical peri-implantation barrier for SCNT embryos.

Going circular: history, present, and future of circRNAs in cancer

To date, thousands of highly abundant and conserved single-stranded RNA molecules shaped into ring structures (circRNAs) have been identified. CircRNAs are multifunctional molecules that have been shown to regulate gene expression transcriptionally and post-transcriptionally and exhibit distinct tissue- and development-specific expression patterns associated with a variety of normal and disease conditions, including cancer pathogenesis. Over the past years, due to their intrinsic stability and resistance to ribonucleases, particular attention has been drawn to their use as reliable diagnostic and prognostic biomarkers in cancer diagnosis, treatment, and prevention. However, there are some critical caveats to their utility in the clinic. Their circular shape limits their annotation and a complete functional elucidation is lacking. This makes their detection and biomedical application still challenging. Herein, we review the current knowledge of circRNA biogenesis and function, and of their involvement in tumorigenesis and potential utility in cancer-targeted therapy.

The Reasoning Engine: A Satisfiability Modulo Theories-Based Framework for Reasoning About Discrete Biological Models

We present a framework called the Reasoning Engine, which implements Satisfiability Modulo Theories (SMT)-based methods within a unified computational environment to address diverse biological analysis problems. The Reasoning Engine was used to reproduce results from key scientific studies, as well as supporting new research in stem cell biology. The framework utilizes an intermediate language for encoding partially specified discrete dynamical systems, which bridges the gap between high-level domain-specific languages and low-level SMT solvers. We provide this framework as open source together with various biological case studies, illustrating the synthesis, enumeration, optimization, and reasoning over models consistent with experimental observations to reveal novel biological insights.

Collective behavior and self-organization in neural rosette morphogenesis

Neural rosettes develop from the self-organization of differentiating human pluripotent stem cells. This process mimics the emergence of the embryonic central nervous system primordium, i.e., the neural tube, whose formation is under close investigation as errors during such process result in severe diseases like spina bifida and anencephaly. While neural tube formation is recognized as an example of self-organization, we still do not understand the fundamental mechanisms guiding the process. Here, we discuss the different theoretical frameworks that have been proposed to explain self-organization in morphogenesis. We show that an explanation based exclusively on stem cell differentiation cannot describe the emergence of spatial organization, and an explanation based on patterning models cannot explain how different groups of cells can collectively migrate and produce the mechanical transformations required to generate the neural tube. We conclude that neural rosette development is a relevant experimental 2D in-vitro model of morphogenesis because it is a multi-scale self-organization process that involves both cell differentiation and tissue development. Ultimately, to understand rosette formation, we first need to fully understand the complex interplay between growth, migration, cytoarchitecture organization, and cell type evolution.

Parallel Genome-Wide CRISPR Screens to Identify State-Dependent Self-Renewal Regulators of Mouse Embryonic Stem Cells

The pluripotency of embryonic stem cells (ESCs) is more accurately viewed as a continuous developmental process rather than a fixed state. However, the factors that play general or state-specific roles in regulating self-renewal in different pluripotency states remain poorly defined. In this study, parallel genome-wide CRISPR/Cas9 knockout (KO) screens were applied in ESCs cultured in the serum plus LIF (SL) and in the 2i plus LIF (2iL) conditions. The candidate genes were classified into seven groups based on their positive or negative effects on self-renewal, and whether this effect was general or state-specific for ESCs under SL and 2iL culture conditions. We characterized the expression and function of genes in these seven groups. The loss of function of novel pluripotent candidate genes Usp28, Zfp598, and Zfp296 was further evaluated in mouse ESCs. Consistent with our screen, the knockout of Usp28 promotes the proliferation of SL-ESCs and 2iL-ESCs, whereas Zfp598 is indispensable for the self-renewal of ESCs under both culture conditions. The cell phenotypes of Zfp296 KO ESCs under SL and 2iL culture conditions were different. Our work provided a valuable resource for dissecting the molecular regulation of ESC self-renewal in different pluripotency states.

Plakoglobin is a mechanoresponsive regulator of naive pluripotency

Biomechanical cues are instrumental in guiding embryonic development and cell differentiation. Understanding how these physical stimuli translate into transcriptional programs will provide insight into mechanisms underlying mammalian pre-implantation development. Here, we explore this type of regulation by exerting microenvironmental control over mouse embryonic stem cells. Microfluidic encapsulation of mouse embryonic stem cells in agarose microgels stabilizes the naive pluripotency network and specifically induces expression of Plakoglobin (Jup), a vertebrate homolog of β-catenin. Overexpression of Plakoglobin is sufficient to fully re-establish the naive pluripotency gene regulatory network under metastable pluripotency conditions, as confirmed by single-cell transcriptome profiling. Finally, we find that, in the epiblast, Plakoglobin was exclusively expressed at the blastocyst stage in human and mouse embryos - further strengthening the link between Plakoglobin and naive pluripotency in vivo. Our work reveals Plakoglobin as a mechanosensitive regulator of naive pluripotency and provides a paradigm to interrogate the effects of volumetric confinement on cell-fate transitions.

STAT3 and HIF1α cooperatively mediate the transcriptional and physiological responses to hypoxia

STAT3 and HIF1α are two fundamental transcription factors involved in many merging processes, like angiogenesis, metabolism, and cell differentiation. Notably, under pathological conditions, the two factors have been shown to interact genetically, but both the molecular mechanisms underlying such interactions and their relevance under physiological conditions remain unclear. In mouse embryonic stem cells (ESCs) we manage to determine the specific subset of hypoxia-induced genes that need STAT3 to be properly transcribed and, among them, fundamental genes like Vegfa, Hk1, Hk2, Pfkp and Hilpda are worth mentioning. Unexpectedly, we also demonstrated that the absence of STAT3 does not affect the expression of Hif1α mRNA nor the stabilization of HIF1α protein, but the STAT3-driven regulation of the hypoxia-dependent subset of gene could rely on the physical interaction between STAT3 and HIF1α. To further elucidate the physiological roles of this STAT3 non-canonical nuclear activity, we used a CRISPR/Cas9 zebrafish stat3 knock-out line. Notably, hypoxia-related fluorescence of the hypoxia zebrafish reporter line (HRE:mCherry) cannot be induced when Stat3 is not active and, while Stat3 Y705 phosphorylation seems to have a pivotal role in this process, S727 does not affect the Stat3-dependent hypoxia response. Hypoxia is fundamental for vascularization, angiogenesis and immune cells mobilization; all processes that, surprisingly, cannot be induced by low oxygen levels when Stat3 is genetically ablated. All in all, here we report the specific STAT3/HIF1α-dependent subset of genes in vitro and, for the first time with an in vivo model, we determined some of the physiological roles of STAT3-hypoxia crosstalk.

The Role of ATP-Binding Cassette Proteins in Stem Cell Pluripotency

Pluripotent stem cells (PSCs) are highly proliferative cells that can self-renew indefinitely in vitro. Upon receiving appropriate signals, PSCs undergo differentiation and can generate every cell type in the body. These unique properties of PSCs require specific gene expression patterns that define stem cell identity and dynamic regulation of intracellular metabolism to support cell growth and cell fate transitions. PSCs are prone to DNA damage due to elevated replicative and transcriptional stress. Therefore, mechanisms to prevent deleterious mutations in PSCs that compromise stem cell function or increase the risk of tumor formation from becoming amplified and propagated to progenitor cells are essential for embryonic development and for using PSCs including induced PSCs (iPSCs) as a cell source for regenerative medicine. In this review, we discuss the role of the ATP-binding cassette (ABC) superfamily in maintaining PSC homeostasis, and propose how their activities can influence cellular signaling and stem cell fate decisions. Finally, we highlight recent discoveries that not all ABC family members perform only canonical metabolite and peptide transport functions in PSCs; rather, they can participate in diverse cellular processes from genome surveillance to gene transcription and mRNA translation, which are likely to maintain the pristine state of PSCs.

Esrrb guides naive pluripotent cells through the formative transcriptional programme

During embryonic development, naive pluripotent epiblast cells transit to a formative state. The formative epiblast cells form a polarized epithelium, exhibit distinct transcriptional and epigenetic profiles and acquire competence to differentiate into all somatic and germline lineages. However, we have limited understanding of how the transition to a formative state is molecularly controlled. Here we used murine embryonic stem cell models to show that ESRRB is both required and sufficient to activate formative genes. Genetic inactivation of Esrrb leads to illegitimate expression of mesendoderm and extra-embryonic markers, impaired formative expression and failure to self-organize in 3D. Functionally, this results in impaired ability to generate formative stem cells and primordial germ cells in the absence of Esrrb. Computational modelling and genomic analyses revealed that ESRRB occupies key formative genes in naive cells and throughout the formative state. In so doing, ESRRB kickstarts the formative transition, leading to timely and unbiased capacity for multi-lineage differentiation.

Identifying strengths and weaknesses of methods for computational network inference from single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) offers unparalleled insight into the transcriptional programs of different cellular states by measuring the transcriptome of thousands of individual cells. An emerging problem in the analysis of scRNA-seq is the inference of transcriptional gene regulatory networks and a number of methods with different learning frameworks have been developed to address this problem. Here, we present an expanded benchmarking study of eleven recent network inference methods on seven published scRNA-seq datasets in human, mouse, and yeast considering different types of gold standard networks and evaluation metrics. We evaluate methods based on their computing requirements as well as on their ability to recover the network structure. We find that, while most methods have a modest recovery of experimentally derived interactions based on global metrics such as Area Under the Precision Recall curve, methods are able to capture targets of regulators that are relevant to the system under study. Among the top performing methods that use only expression were SCENIC, PIDC, MERLIN or Correlation. Addition of prior biological knowledge and the estimation of transcription factor activities resulted in the best overall performance with the Inferelator and MERLIN methods that use prior knowledge outperforming methods that use expression alone. We found that imputation for network inference did not improve network inference accuracy and could be detrimental. Comparisons of inferred networks for comparable bulk conditions showed that the networks inferred from scRNA-seq datasets are often better or at par with the networks inferred from bulk datasets. Our analysis should be beneficial in selecting methods for network inference. At the same time, this highlights the need for improved methods and better gold standards for regulatory network inference from scRNAseq datasets.

Zeb2 DNA-Binding Sites in Neuroprogenitor Cells Reveal Autoregulation and Affirm Neurodevelopmental Defects, Including in Mowat-Wilson Syndrome

Functional perturbation and action mechanism studies have shown that the transcription factor Zeb2 controls cell fate decisions, differentiation, and/or maturation in multiple cell lineages in embryos and after birth. In cultured embryonic stem cells (ESCs), Zeb2’s mRNA/protein upregulation is necessary for the exit from primed pluripotency and for entering general and neural differentiation. We edited mouse ESCs to produce Flag-V5 epitope-tagged Zeb2 protein from one endogenous allele. Using chromatin immunoprecipitation coupled with sequencing (ChIP-seq), we mapped 2432 DNA-binding sites for this tagged Zeb2 in ESC-derived neuroprogenitor cells (NPCs). A new, major binding site maps promoter-proximal to Zeb2 itself. The homozygous deletion of this site demonstrates that autoregulation of Zeb2 is necessary to elicit the appropriate Zeb2-dependent effects in ESC-to-NPC differentiation. We have also cross-referenced all the mapped Zeb2 binding sites with previously obtained transcriptome data from Zeb2 perturbations in ESC-derived NPCs, GABAergic interneurons from the ventral forebrain of mouse embryos, and stem/progenitor cells from the post-natal ventricular-subventricular zone (V-SVZ) in mouse forebrain, respectively. Despite the different characteristics of each of these neurogenic systems, we found interesting target gene overlaps. In addition, our study also contributes to explaining developmental disorders, including Mowat-Wilson syndrome caused by ZEB2 deficiency, and also other monogenic syndromes.

Chemical conversion of human conventional PSCs to TSCs following transient naive gene activation

In human embryos, naive pluripotent cells of the inner cell mass (ICM) generate epiblast, primitive endoderm and trophectoderm (TE) lineages, whence trophoblast cells derive. In vitro, naive pluripotent stem cells (PSCs) retain this potential and efficiently generate trophoblast stem cells (TSCs), while conventional PSCs form TSCs at low efficiency. Transient histone deacetylase and MEK inhibition combined with LIF stimulation is used to chemically reset conventional to naive PSCs. Here, we report that chemical resetting induces the expression of both naive and TSC markers and of placental imprinted genes. A modified chemical resetting protocol allows for the fast and efficient conversion of conventional PSCs into TSCs, entailing shutdown of pluripotency genes and full activation of the trophoblast master regulators, without induction of amnion markers. Chemical resetting generates a plastic intermediate state, characterised by co-expression of naive and TSC markers, after which cells steer towards one of the two fates in response to the signalling environment. The efficiency and rapidity of our system will be useful to study cell fate transitions and to generate models of placental disorders.

Transcription factor SOX15 regulates stem cell pluripotency and promotes neural fate during differentiation by activating the neurogenic gene Hes5

SOX2 and SOX15 are Sox family transcription factors enriched in embryonic stem cells (ESCs). The role of SOX2 in activating gene expression programs essential for stem cell self-renewal and acquisition of pluripotency during somatic cell reprogramming is well-documented. However, the contribution of SOX15 to these processes is unclear and often presumed redundant with SOX2 largely because overexpression of SOX15 can partially restore self-renewal in SOX2-deficient ESCs. Here, we show that SOX15 contributes to stem cell maintenance by cooperating with ESC-enriched transcriptional coactivators to ensure optimal expression of pluripotency-associated genes. We demonstrate that SOX15 depletion compromises reprogramming of fibroblasts to pluripotency which cannot be compensated by SOX2. Ectopic expression of SOX15 promotes the reversion of a postimplantation, epiblast stem cell state back to a preimplantation, ESC-like identity even though SOX2 is expressed in both cell states. We also uncover a role of SOX15 in lineage specification, by showing that loss of SOX15 leads to defects in commitment of ESCs to neural fates. SOX15 promotes neural differentiation by binding to and activating a previously uncharacterized distal enhancer of a key neurogenic regulator, Hes5. Together, these findings identify a multifaceted role of SOX15 in induction and maintenance of pluripotency and neural differentiation.

Label-free and non-destructive identification of naïve and primed embryonic stem cells based on differences in cellular metabolism

Pluripotent stem cells (PSCs) exist in naïve or primed states based on their origin. For in vitro culture, these PSCs require different supplements and growth factors. However, owing to their similar phenotypic features, identifying both cell types without harming cellular functions is challenging. This study reports an electrochemical method that enables simple, label-free, and non-destructive detection of naïve embryonic stem cells (ESCs) derived from mouse ESCs, based on the differences in cellular metabolism. Two major metabolic pathways to generate adenosine triphosphate (ATP)-glycolysis and oxidative phosphorylation (OXPHOS)-were blocked, and it was found that mitochondrial energy generation is the origin of the strong electrochemical signals of naïve ESCs. The number of ESCs is quantified when mixed with primed ESCs or converted from naïve-primed switchable metastable ESCs. The mouse PSCs derived from doxycycline-inducible mouse embryonic fibroblasts (MEFs) are also sensitively identified among other cell types such as unconverted MEFs and primed PSCs. The developed sensing platform operates in a non-invasive and label-free manner. Thus, it can be useful in the development of stem cell-derived therapeutics.

Molecular versatility during pluripotency progression

A challenge during development is to ensure lineage segregation while preserving plasticity. Using pluripotency progression as a paradigm, we review how developmental transitions are coordinated by redeployments, rather than global resettings, of cellular components. We highlight how changes in response to extrinsic cues (FGF, WNT, Activin/Nodal, Netrin-1), context- and stoichiometry-dependent action of transcription factors (Oct4, Nanog) and reconfigurations of epigenetic regulators (enhancers, promoters, TrxG, PRC) may confer robustness to naïve to primed pluripotency transition. We propose the notion of Molecular Versatility to regroup mechanisms by which molecules are repurposed to exert different, sometimes opposite, functions in close stem cell configurations.

Requirement for STAT3 and its target, TFCP2L1, in self-renewal of naïve pluripotent stem cells in vivo and in vitro

We previously demonstrated gradual loss of epiblast during diapause in embryos lacking components of the LIF/IL6 receptor. Here, we explore the requirement for the downstream signalling transducer andactivator of transcription STAT3 and its target, TFCP2L1, in maintenance of naïve pluripotency. Unlike conventional markers, such as NANOG, which remains high in epiblast until implantation, both STAT3 and TFCP2L1 proteins decline during blastocyst expansion, but intensify in the embryonic region after induction of diapause, as observed visually and confirmed using our image-analysis pipeline, consistent with our previous transcriptional expression data. Embryos lacking STAT3 or TFCP2L1 underwent catastrophic loss of most of the inner cell mass during the first few days of diapause, indicating involvement of signals in addition to LIF/IL6 for sustaining naïve pluripotency in vivo. By blocking MEK/ERK signalling from the morula stage, we could derive embryonic stem cells with high efficiency from STAT3 null embryos, but not those lacking TFCP2L1, suggesting a hitherto unknown additional role for this essential STAT3 target in transition from embryo to embryonic stem cells in vitro. This article has an associated First Person interview with the first author of the paper.

circTMEM181 upregulates ARHGAP29 to inhibit hepatocellular carcinoma migration and invasion by sponging miR-519a-5p

AIM

Circular RNAs (circRNAs) are a novel class of noncoding RNAs and are conserved in various species. Although numerous circRNAs have been identified, their role in cancer remains unclear.

METHODS

The expression of circTMEM181 in 90 paired human hepatocellular carcinoma (HCC) and adjacent nontumor tissues were detected using quantitative reverse transcription-polymerase chain reaction. Transwell assay was performed for functional analysis of HCC cell migration and invasion. Luciferase reporter assay was used to verify the combination of circTMEM181 and miR-519a-5p.

RESULTS

In this study, we identified a novel circRNA, named circTMEM181, was downregulated in HCC tissues. Decreased expression of circTMEM181 was associated with shorter overall survival of patients with HCC. CircTMEM181 overexpression reduced HCC cell migration and invasion abilities, while circTMEM181 knockdown increased cell motility. Mechanically, circTMEM181 could directly bind to miR-519a-5p and subsequently upregulate ARHGAP29 protein expression.

CONCLUSION

These data provide the first evidence of clinical significance and function of circTMEM181, and suggest the circTMEM181/miR-519a-5p/ARHGAP29 axis in HCC development.

Virtual cells in a virtual microenvironment recapitulate early development-like patterns in human pluripotent stem cell colonies

The mechanism by which morphogenetic signals engage the regulatory networks responsible for early embryonic tissue patterning is incompletely understood. Here, we developed a minimal gene regulatory network (GRN) model of human pluripotent stem cell (hPSC) lineage commitment and embedded it into "cellular" agents that respond to a dynamic morphogenetic signaling microenvironment. Simulations demonstrated that GRN wiring had significant non-intuitive effects on tissue pattern order, composition, and dynamics. Experimental perturbation of GRN connectivities supported model predictions and demonstrated the role of OCT4 as a master regulator of peri-gastrulation fates. Our so-called GARMEN strategy provides a multiscale computational platform to understand how single-cell-based regulatory interactions scale to tissue domains. This foundation provides new opportunities to simulate the impact of network motifs on normal and aberrant tissue development.

Transient upregulation of IRF1 during exit from naive pluripotency confers viral protection

Stem cells intrinsically express a subset of genes which are normally associated with interferon stimulation and the innate immune response. However, the expression of these interferon-stimulated genes (ISG) in stem cells is independent from external stimuli such as viral infection. Here, we show that the interferon regulatory factor 1, Irf1, is directly controlled by the murine formative pluripotency gene regulatory network and transiently upregulated during the transition from naive to formative pluripotency. IRF1 binds to regulatory regions of a conserved set of ISGs and is required for their faithful expression upon exit from naive pluripotency. We show that in the absence of IRF1, cells exiting the naive pluripotent stem cell state are more susceptible to viral infection. Irf1 therefore acts as a link between the formative pluripotency network, regulation of innate immunity genes, and defense against viral infections during formative pluripotency.

Major transcriptomic, epigenetic and metabolic changes underlie the pluripotency continuum in rabbit preimplantation embryos

Despite the growing interest in the rabbit model for developmental and stem cell biology, the characterization of embryos at the molecular level is still poorly documented. We conducted a transcriptome analysis of rabbit preimplantation embryos from E2.7 (morula stage) to E6.6 (early primitive streak stage) using bulk and single-cell RNA-sequencing. In parallel, we studied oxidative phosphorylation and glycolysis, and analysed active and repressive epigenetic modifications during blastocyst formation and expansion. We generated a transcriptomic, epigenetic and metabolic map of the pluripotency continuum in rabbit preimplantation embryos, and identified novel markers of naive pluripotency that might be instrumental for deriving naive pluripotent stem cell lines. Although the rabbit is evolutionarily closer to mice than to primates, we found that the transcriptome of rabbit epiblast cells shares common features with those of humans and non-human primates.

Chimpanzee and pig-tailed macaque iPSCs: Improved culture and generation of primate cross-species embryos

As our closest living relatives, non-human primates uniquely enable explorations of human health, disease, development, and evolution. Considerable effort has thus been devoted to generating induced pluripotent stem cells (iPSCs) from multiple non-human primate species. Here, we establish improved culture methods for chimpanzee (Pan troglodytes) and pig-tailed macaque (Macaca nemestrina) iPSCs. Such iPSCs spontaneously differentiate in conventional culture conditions, but can be readily propagated by inhibiting endogenous WNT signaling. As a unique functional test of these iPSCs, we injected them into the pre-implantation embryos of another non-human species, rhesus macaques (Macaca mulatta). Ectopic expression of gene BCL2 enhances the survival and proliferation of chimpanzee and pig-tailed macaque iPSCs within the pre-implantation embryo, although the identity and long-term contribution of the transplanted cells warrants further investigation. In summary, we disclose transcriptomic and proteomic data, cell lines, and cell culture resources that may be broadly enabling for non-human primate iPSCs research.

The Divergent Pluripotent States in Mouse and Human Cells

Pluripotent stem cells (PSCs), which can self-renew and give rise to all cell types in all three germ layers, have great potential in regenerative medicine. Recent studies have shown that PSCs can have three distinct but interrelated pluripotent states: naive, formative, and primed. The PSCs of each state are derived from different stages of the early developing embryo and can be maintained in culture by different molecular mechanisms. In this review, we summarize the current understanding on features of the three pluripotent states and review the underlying molecular mechanisms of maintaining their identities. Lastly, we discuss the interrelation and transition among these pluripotency states. We believe that comprehending the divergence of pluripotent states is essential to fully harness the great potential of stem cells in regenerative medicine.

CELLoGeNe - An energy landscape framework for logical networks controlling cell decisions

Experimental and computational efforts are constantly made to elucidate mechanisms controlling cell fate decisions during development and reprogramming. One powerful computational method is to consider cell commitment and reprogramming as movements in an energy landscape. Here, we develop Computation of Energy Landscapes of Logical Gene Networks (CELLoGeNe), which maps Boolean implementation of gene regulatory networks (GRNs) into energy landscapes. CELLoGeNe removes inadvertent symmetries in the energy landscapes normally arising from standard Boolean operators. Furthermore, CELLoGeNe provides tools to visualize and stochastically analyze the shapes of multi-dimensional energy landscapes corresponding to epigenetic landscapes for development and reprogramming. We demonstrate CELLoGeNe on two GRNs governing different aspects of induced pluripotent stem cells, identifying experimentally validated attractors and revealing potential reprogramming roadblocks. CELLoGeNe is a general framework that can be applied to various biological systems offering a broad picture of intracellular dynamics otherwise inaccessible with existing methods.

•

CELLoGeNe – Computation of Energy Landscapes of Logical Gene Networks

•

Cell states as landscape attractors

•

Maintenance and acquisition of cell pluripotency applications

•

Single cell stochastic landscape navigation and visualization tool

Cell cycle gene regulation dynamics revealed by RNA velocity and deep-learning

Despite the fact that the cell cycle is a fundamental process of life, a detailed quantitative understanding of gene regulation dynamics throughout the cell cycle is far from complete. Single-cell RNA-sequencing (scRNA-seq) technology gives access to these dynamics without externally perturbing the cell. Here, by generating scRNA-seq libraries in different cell systems, we observe cycling patterns in the unspliced-spliced RNA space of cell cycle-related genes. Since existing methods to analyze scRNA-seq are not efficient to measure cycling gene dynamics, we propose a deep learning approach (DeepCycle) to fit these patterns and build a high-resolution map of the entire cell cycle transcriptome. Characterizing the cell cycle in embryonic and somatic cells, we identify major waves of transcription during the G1 phase and systematically study the stages of the cell cycle. Our work will facilitate the study of the cell cycle in multiple cellular models and different biological contexts.

Single-cell RNA-sequencing technology gives access to cell cycle dynamics without externally perturbing the cell. Here the authors present DeepCycle,a robust deep learning method to infer the cell cycle state in single cells from scRNA-seq data.

Cnot8 eliminates naïve regulation networks and is essential for naïve-to-formative pluripotency transition

Mammalian early epiblasts at different phases are characterized by naïve, formative, and primed pluripotency states, involving extensive transcriptome changes. Here, we report that deadenylase Cnot8 of Ccr4-Not complex plays essential roles during the transition from naïve to formative state. Knock out (KO) Cnot8 resulted in early embryonic lethality in mice, but Cnot8 KO embryonic stem cells (ESCs) could be established. Compared with the cells differentiated from normal ESCs, Cnot8 KO cells highly expressed a great many genes during their differentiation into the formative state, including several hundred naïve-like genes enriched in lipid metabolic process and gene expression regulation that may form the naïve regulation networks. Knockdown expression of the selected genes of naïve regulation networks partially rescued the differentiation defects of Cnot8 KO ESCs. Cnot8 depletion led to the deadenylation defects of its targets, increasing their poly(A) tail lengths and half-life, eventually elevating their expression levels. We further found that Cnot8 was involved in the clearance of targets through its deadenylase activity and the binding of Ccr4-Not complex, as well as the interacting with Tob1 and Pabpc1. Our results suggest that Cnot8 eliminates naïve regulation networks through mRNA clearance, and is essential for naïve-to-formative pluripotency transition.

Smarcb1 Loss Results in a Deregulation of esBAF Binding and Impacts the Expression of Neurodevelopmental Genes

The murine esBAF complex plays a major role in the regulation of gene expression during stem cell development and differentiation. As one of its core subunits, Smarcb1 is indispensable for its function and its loss is connected to neurodevelopmental disorders and participates in the carcinogenesis of entities such as rhabdoid tumours. We explored how Smarcb1 regulates gene programs in murine embryonic stem cells (mESC) and in this way orchestrates differentiation. Our data underline the importance of Smarcb1 expression and function for the development of the nervous system along with basic cellular functions, such as cell adhesion and cell organisation. Using ChIP-seq, we were able to portray the consequences of Smarcb1 knockdown (kd) for the binding of esBAF and PRC2 as well as its influence on histone marks H3K27me3, H3K4me3 and H3K27ac. Their signals are changed in gene and enhancer regions of genes connected to nervous system development and offers a plausible explanation for changes in gene expression. Further, we describe a group of genes that are, despite increased BAF binding, suppressed after Smarcb1 kd by mechanisms independent of PRC2 function.

Modeling the C. elegans germline stem cell genetic network using automated reasoning

Computational methods and tools are a powerful complementary approach to experimental work for studying regulatory interactions in living cells and systems. We demonstrate the use of formal reasoning methods as applied to the Caenorhabditis elegans germ line, which is an accessible system for stem cell research. The dynamics of the underlying genetic networks and their potential regulatory interactions are key for understanding mechanisms that control cellular decision-making between stem cells and differentiation. We model the “stem cell fate” versus entry into the “meiotic development” pathway decision circuit in the young adult germ line based on an extensive study of published experimental data and known/hypothesized genetic interactions. We apply a formal reasoning framework to derive predictive networks for control of differentiation. Using this approach we simultaneously specify many possible scenarios and experiments together with potential genetic interactions, and synthesize genetic networks consistent with all encoded experimental observations. In silico analysis of knock-down and overexpression experiments within our model recapitulate published phenotypes of mutant animals and can be applied to make predictions on cellular decision-making. A methodological contribution of this work is demonstrating how to effectively model within a formal reasoning framework a complex genetic network with a wealth of known experimental data and constraints. We provide a summary of the steps we have found useful for the development and analysis of this model and can potentially be applicable to other genetic networks. This work also lays a foundation for developing realistic whole tissue models of the C. elegans germ line where each cell in the model will execute a synthesized genetic network.