Cell cycle dynamics of mouse embryonic stem cells in the ground state and during transition to formative pluripotency

Single-cell nascent RNA sequencing unveils coordinated global transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1,2. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations3. However, fundamental questions about the temporal regulation of transcription and enhancer-gene coordination remain unanswered, primarily because of the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq-a new single-cell nascent RNA sequencing assay that uses click chemistry-and unveil coordinated transcription throughout the genome. We demonstrate the episodic nature of transcription and the co-transcription of functionally related genes. scGRO-seq can estimate burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells and can leverage replication-dependent non-polyadenylated histone gene transcription to elucidate cell cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq enables the identification of networks of enhancers and genes. Our results suggest that the bursting of transcription at super-enhancers precedes bursting from associated genes. By imparting insights into the dynamic nature of global transcription and the origin and propagation of transcription signals, we demonstrate the ability of scGRO-seq to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

A 3D micro-printed single cell micro-niche with asymmetric niche signals - An in vitro model for asymmetric cell division study

Intricate microenvironment signals orchestrate to affect cell behavior and fate during tissue morphogenesis. However, the underlying mechanisms on how specific local niche signals influence cell behavior and fate are not fully understood, owing to the lack of in vitro platform able to precisely, quantitatively, spatially, and independently manipulate individual niche signals. Here, microarrays of protein-based 3D single cell micro-niche (3D-SCμN), with precisely engineered biophysical and biochemical niche signals, are micro-printed by a multiphoton microfabrication and micropatterning technology. Mouse embryonic stem cell (mESC) is used as the model cell to study how local niche signals affect stem cell behavior and fate. By precisely engineering the internal microstructures of the 3D SCμNs, we demonstrate that the cell division direction can be controlled by the biophysical niche signals, in a cell shape-independent manner. After confining the cell division direction to a dominating axis, single mESCs are exposed to asymmetric biochemical niche signals, specifically, cell-cell adhesion molecule on one side and extracellular matrix on the other side. We demonstrate that, symmetry-breaking (asymmetric) niche signals successfully trigger cell polarity formation and bias the orientation of asymmetric cell division, the mitosis process resulting in two daughter cells with differential fates, in mESCs.

Transitions in development - an interview with Ariel Waisman

Ariel Waisman is a CONICET Junior Researcher in the Laboratory of Applied Research in Neurosciences at FLENI in Buenos Aires, Argentina. Ariel's group studies gene regulatory networks in human pluripotent stem cells to address mechanisms of development and cardiac differentiation, among other topics. We spoke to Ariel over Teams to learn more about his career path, the research interests in his group, and the challenges faced by researchers in Argentina and the Global South.

The transcription factor OCT6 promotes the dissolution of the naïve pluripotent state by repressing Nanog and activating a formative state gene regulatory network

In the mouse embryo, the transition from the preimplantation to the postimplantation epiblast is governed by changes in the gene regulatory network (GRN) that lead to transcriptional, epigenetic, and functional changes. This transition can be faithfully recapitulated in vitro by the differentiation of mouse embryonic stem cells (mESCs) to epiblast-like cells (EpiLCs), that reside in naïve and formative states of pluripotency, respectively. However, the GRN that drives this conversion is not fully elucidated. Here we demonstrate that the transcription factor OCT6 is a key driver of this process. Firstly, we show that Oct6 is not expressed in mESCs but is rapidly induced as cells exit the naïve pluripotent state. By deleting Oct6 in mESCs, we find that knockout cells fail to acquire the typical morphological changes associated with the formative state when induced to differentiate. Additionally, the key naïve pluripotency TFs Nanog, Klf2, Nr5a2, Prdm14, and Esrrb were expressed at higher levels than in wild-type cells, indicating an incomplete dismantling of the naïve pluripotency GRN. Conversely, premature expression of Oct6 in naïve cells triggered a rapid morphological transformation mirroring differentiation, that was accompanied by the upregulation of the endogenous Oct6 as well as the formative genes Sox3, Zic2/3, Foxp1, Dnmt3A and FGF5. Strikingly, we found that OCT6 represses Nanog in a bistable manner and that this regulation is at the transcriptional level. Moreover, our findings also reveal that Oct6 is repressed by NANOG. Collectively, our results establish OCT6 as a key TF in the dissolution of the naïve pluripotent state and support a model where Oct6 and Nanog form a double negative feedback loop which could act as an important toggle mediating the transition to the formative state.

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.

Dynamic Roles of Signaling Pathways in Maintaining Pluripotency of Mouse and Human Embryonic Stem Cells

Culturing of mouse and human embryonic stem cells (ESCs) in vitro was a major breakthrough in the field of stem cell biology. These models gained popularity very soon mainly due to their pluripotency. Evidently, the ESCs of mouse and human origin share typical phenotypic responses due to their pluripotent nature, such as self-renewal capacity and potency. The conserved network of core transcription factors regulates these responses. However, significantly different signaling pathways and upstream transcriptional networks regulate expression and activity of these core pluripotency factors in ESCs of both the species. In fact, ample evidence shows that a pathway, which maintains pluripotency in mouse ESCs, promotes differentiation in human ESCs. In this review, we discuss the role of canonical signaling pathways implicated in regulation of pluripotency and differentiation particularly in mouse and human ESCs. We believe that understanding these distinct and at times-opposite mechanisms-is critical for the progress in the field of stem cell biology and regenerative medicine.

Single-cell nascent RNA sequencing using click-chemistry unveils coordinated transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1-5. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations6-9. However, fundamental questions in the temporal regulation of transcription and enhancer-gene synchrony remain unanswered primarily due to the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq - a novel single-cell nascent RNA sequencing assay using click-chemistry - and unveil the coordinated transcription throughout the genome. scGRO-seq demonstrates the episodic nature of transcription, and estimates burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells. It reveals the co-transcription of functionally related genes and leverages the replication-dependent non-polyadenylated histone genes transcription to elucidate cell-cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq identifies networks of enhancers and genes and indicates that the bursting of transcription at super-enhancers precedes the burst from associated genes. By imparting insights into the dynamic nature of transcription and the origin and propagation of transcription signals, scGRO-seq demonstrates its unique ability to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

The dynamical organization of the core pluripotency transcription factors responds to differentiation cues in early S-phase

DNA replication in stem cells is a major challenge for pluripotency preservation and cell fate decisions. This process involves massive changes in the chromatin architecture and the reorganization of many transcription-related molecules in different spatial and temporal scales. Pluripotency is controlled by the master transcription factors (TFs) OCT4, SOX2 and NANOG that partition into condensates in the nucleus of embryonic stem cells. These condensates are proposed to play relevant roles in the regulation of gene expression and the maintenance of pluripotency. Here, we asked whether the dynamical distribution of the pluripotency TFs changes during the cell cycle, particularly during DNA replication. Since the S phase is considered to be a window of opportunity for cell fate decisions, we explored if differentiation cues in G1 phase trigger changes in the distribution of these TFs during the subsequent S phase. Our results show a spatial redistribution of TFs condensates during DNA replication which was not directly related to chromatin compaction. Additionally, fluorescence fluctuation spectroscopy revealed TF-specific, subtle changes in the landscape of TF-chromatin interactions, consistent with their particularities as key players of the pluripotency network. Moreover, we found that differentiation stimuli in the preceding G1 phase triggered a relatively fast and massive reorganization of pluripotency TFs in early-S phase. Particularly, OCT4 and SOX2 condensates dissolved whereas the lifetimes of TF-chromatin interactions increased suggesting that the reorganization of condensates is accompanied with a change in the landscape of TF-chromatin interactions. Notably, NANOG showed impaired interactions with chromatin in stimulated early-S cells in line with its role as naïve pluripotency TF. Together, these findings provide new insights into the regulation of the core pluripotency TFs during DNA replication of embryonic stem cells and highlight their different roles at early differentiation stages.

Transcriptional Fluctuations Govern the Serum-Dependent Cell Cycle Duration Heterogeneities in Mammalian Cells

Mammalian cells exhibit a high degree of intercellular variability in cell cycle period and phase durations. However, the factors orchestrating the cell cycle duration heterogeneities remain unclear. Herein, by combining cell cycle network-based mathematical models with live single-cell imaging studies under varied serum conditions, we demonstrate that fluctuating transcription rates of cell cycle regulatory genes across cell lineages and during cell cycle progression in mammalian cells majorly govern the robust correlation patterns of cell cycle period and phase durations among sister, cousin, and mother-daughter lineage pairs. However, for the overall cellular population, alteration in the serum level modulates the fluctuation and correlation patterns of cell cycle period and phase durations in a correlated manner. These heterogeneities at the population level can be fine-tuned under limited serum conditions by perturbing the cell cycle network using a p38-signaling inhibitor without affecting the robust lineage-level correlations. Overall, our approach identifies transcriptional fluctuations as the key controlling factor for the cell cycle duration heterogeneities and predicts ways to reduce cell-to-cell variabilities by perturbing the cell cycle network regulations.

Dynamical modeling of the H3K27 epigenetic landscape in mouse embryonic stem cells

The Polycomb system via the methylation of the lysine 27 of histone H3 (H3K27) plays central roles in the silencing of many lineage-specific genes during development. Recent experimental evidence suggested that the recruitment of histone modifying enzymes like the Polycomb repressive complex 2 (PRC2) at specific sites and their spreading capacities from these sites are key to the establishment and maintenance of a proper epigenomic landscape around Polycomb-target genes. Here, to test whether such mechanisms, as a minimal set of qualitative rules, are quantitatively compatible with data, we developed a mathematical model that can predict the locus-specific distributions of H3K27 modifications based on previous biochemical knowledge. Within the biological context of mouse embryonic stem cells, our model showed quantitative agreement with experimental profiles of H3K27 acetylation and methylation around Polycomb-target genes in wild-type and mutants. In particular, we demonstrated the key role of the reader-writer module of PRC2 and of the competition between the binding of activating and repressing enzymes in shaping the H3K27 landscape around transcriptional start sites. The predicted dynamics of establishment and maintenance of the repressive trimethylated H3K27 state suggest a slow accumulation, in perfect agreement with experiments. Our approach represents a first step towards a quantitative description of PcG regulation in various cellular contexts and provides a generic framework to better characterize epigenetic regulation in normal or disease situations.

Transcriptional heterogeneity and cell cycle regulation as central determinants of primitive endoderm priming

During embryonic development cells acquire identity as they proliferate, implying that an intrinsic facet of cell fate choice requires coupling lineage decisions to cell division. How is the cell cycle regulated to promote or suppress heterogeneity and differentiation? We explore this question combining time lapse imaging with single-cell RNA-seq in the contexts of self-renewal, priming, and differentiation of mouse embryonic stem cells (ESCs) towards the Primitive Endoderm (PrE) lineage. Since ESCs are derived from the inner cell mass (ICM) of the mammalian blastocyst, ESCs in standard culture conditions are transcriptionally heterogeneous containing dynamically interconverting subfractions primed for either of the two ICM lineages, Epiblast and PrE. Here, we find that differential regulation of cell cycle can tip the balance between these primed populations, such that naïve ESC culture promotes Epiblast-like expansion and PrE differentiation stimulates the selective survival and proliferation of PrE-primed cells. In endoderm differentiation, this change is accompanied by a counter-intuitive increase in G1 length, also observed in vivo. While fibroblast growth factor/extracellular signal-regulated kinase (FGF/ERK) signalling is a key regulator of ESC differentiation and PrE specification, we find it is not just responsible for ESCs heterogeneity, but also the inheritance of similar cell cycles between sisters and cousins. Taken together, our results indicate a tight relationship between transcriptional heterogeneity and cell cycle regulation in lineage specification, with primed cell populations providing a pool of flexible cell types that can be expanded in a lineage-specific fashion while allowing plasticity during early determination.

Esrrb is a cell-cycle-dependent associated factor balancing pluripotency and XEN differentiation

Cell cycle and differentiation decisions are linked; however, the underlying principles that drive these decisions are unclear. Here, we combined cell-cycle reporter system and single-cell RNA sequencing (scRNA-seq) profiling to study the transcriptomes of embryonic stem cells (ESCs) in the context of cell-cycle states and differentiation. By applying retinoic acid, to G1 and G2/M ESCs, we show that, while both populations can differentiate toward epiblast stem cells (EpiSCs), only G2/M ESCs could differentiate into extraembryonic endoderm cells. We identified Esrrb, a pluripotency factor that is upregulated during G2/M, as a driver of extraembryonic endoderm stem cell (XEN) differentiation. Furthermore, enhancer chromatin states based on wild-type (WT) and ESRRB knockout (KO) ESCs show association of ESRRB with XEN poised enhancers. G1 cells overexpressing Esrrb allow ESCs to produce XENs, while ESRRB-KO ESCs lost their potential to differentiate into XEN. Overall, this study reveals a vital link between Esrrb and cell-cycle states during the exit from pluripotency.

A data-driven method to learn a jump diffusion process from aggregate biological gene expression data

Dynamic models of gene expression are urgently required. In this paper, we describe the time evolution of gene expression by learning a jump diffusion process to model the biological process directly. Our algorithm needs aggregate gene expression data as input and outputs the parameters of the jump diffusion process. The learned jump diffusion process can predict population distributions of gene expression at any developmental stage, obtain long-time trajectories for individual cells, and offer a novel approach to computing RNA velocity. Moreover, it studies biological systems from a stochastic dynamic perspective. Gene expression data at a time point, which is a snapshot of a cellular process, is treated as an empirical marginal distribution of a stochastic process. The Wasserstein distance between the empirical distribution and predicted distribution by the jump diffusion process is minimized to learn the dynamics. For the learned jump diffusion process, its trajectories correspond to the development process of cells, the stochasticity determines the heterogeneity of cells, its instantaneous rate of state change can be taken as "RNA velocity", and the changes in scales and orientations of clusters can be noticed too. We demonstrate that our method can recover the underlying nonlinear dynamics better compared to previous parametric models and the diffusion processes driven by Brownian motion for both synthetic and real world datasets. Our method is also robust to perturbations of data because the computation involves only population expectations.

Imaging developmental cell cycles

The last decade has seen a major expansion in development of live biosensors, the tools needed to genetically encode them into model organisms, and the microscopic techniques used to visualize them. When combined, these offer us powerful tools with which to make fundamental discoveries about complex biological processes. In this review, we summarize the availability of biosensors to visualize an essential cellular process, the cell cycle, and the techniques for single-cell tracking and quantification of these reporters. We also highlight studies investigating the connection of cellular behavior to the cell cycle, particularly through live imaging, and anticipate exciting discoveries with the combination of these technologies in developmental contexts.

Embryonic diapause in mammals and dormancy in embryonic stem cells with the European roe deer as experimental model

In species displaying embryonic diapause, the developmental pace of the embryo is either temporarily and reversibly halted or largely reduced. Only limited knowledge on its regulation and the inhibition of cell proliferation extending pluripotency is available. In contrast with embryos from other diapausing species that reversibly halt during diapause, embryos of the roe deer Capreolus capreolus slowly proliferate over a period of 4-5 months to reach a diameter of approximately 4mm before elongation. The diapausing roe deer embryos present an interesting model species for research on preimplantation developmental progression. Based on our and other research, we summarise the available knowledge and indicate that the use of embryonic stem cells (ESCs) would help to increase our understanding of embryonic diapause. We report on known molecular mechanisms regulating embryonic diapause, as well as cellular dormancy of pluripotent cells. Further, we address the promising application of ESCs to study embryonic diapause, and highlight the current knowledge on the cellular microenvironment regulating embryonic diapause and cellular dormancy.

Opposing Effects of Cohesin and Transcription on CTCF Organization Revealed by Super-resolution Imaging

CCCTC-binding factor (CTCF) and cohesin play critical roles in organizing mammalian genomes into topologically associating domains (TADs). Here, by combining genetic engineering with quantitative super-resolution stimulated emission depletion (STED) microscopy, we demonstrate that in living cells, CTCF forms clusters typically containing 2-8 molecules. A fraction of CTCF clusters, enriched for those with ≥3 molecules, are coupled with cohesin complexes with a characteristic physical distance suggestive of a defined molecular interaction. Acute degradation of the cohesin unloader WAPL or transcriptional inhibition (TI) result in increased CTCF clustering. Furthermore, the effect of TI on CTCF clusters is alleviated by the acute loss of the cohesin subunit SMC3. Our study provides quantitative characterization of CTCF clusters in living cells, uncovers the opposing effects of cohesin and transcription on CTCF clustering, and highlights the power of quantitative super-resolution microscopy as a tool to bridge the gap between biochemical and genomic methodologies in chromatin research.

SMC5/6 is required for replication fork stability and faithful chromosome segregation during neurogenesis

Mutations of SMC5/6 components cause developmental defects, including primary microcephaly. To model neurodevelopmental defects, we engineered a mouse wherein Smc5 is conditionally knocked out (cKO) in the developing neocortex. Smc5 cKO mice exhibited neurodevelopmental defects due to neural progenitor cell (NPC) apoptosis, which led to reduction in cortical layer neurons. Smc5 cKO NPCs formed DNA bridges during mitosis and underwent chromosome missegregation. SMC5/6 depletion triggers a CHEK2-p53 DNA damage response, as concomitant deletion of the Trp53 tumor suppressor or Chek2 DNA damage checkpoint kinase rescued Smc5 cKO neurodevelopmental defects. Further assessment using Smc5 cKO and auxin-inducible degron systems demonstrated that absence of SMC5/6 leads to DNA replication stress at late-replicating regions such as pericentromeric heterochromatin. In summary, SMC5/6 is important for completion of DNA replication prior to entering mitosis, which ensures accurate chromosome segregation. Thus, SMC5/6 functions are critical in highly proliferative stem cells during organism development.

Morphine leads to global genome changes in H3K27me3 levels via a Polycomb Repressive Complex 2 (PRC2) self-regulatory mechanism in mESCs

Background

Environmentally induced epigenetic changes can lead to health problems or disease, but the mechanisms involved remain unclear. Morphine can pass through the placental barrier leading to abnormal embryo development. However, the mechanism by which morphine causes these effects and how they sometimes persist into adulthood is not well known. To unravel the morphine-induced chromatin alterations involved in aberrant embryo development, we explored the role of the H3K27me3/PRC2 repressive complex in gene expression and its transmission across cellular generations in response to morphine.

Results

Using mouse embryonic stem cells as a model system, we found that chronic morphine treatment induces a global downregulation of the histone modification H3K27me3. Conversely, ChIP-Seq showed a remarkable increase in H3K27me3 levels at specific genomic sites, particularly promoters, disrupting selective target genes related to embryo development, cell cycle and metabolism. Through a self-regulatory mechanism, morphine downregulated the transcription of PRC2 components responsible for H3K27me3 by enriching high H3K27me3 levels at the promoter region. Downregulation of PRC2 components persisted for at least 48 h (4 cell cycles) following morphine removal, though promoter H3K27me3 levels returned to control levels.

Conclusions

Morphine induces targeting of the PRC2 complex to selected promoters, including those of PRC2 components, leading to characteristic changes in gene expression and a global reduction in H3K27me3. Following morphine removal, enhanced promoter H3K27me3 levels revert to normal sooner than global H3K27me3 or PRC2 component transcript levels. We suggest that H3K27me3 is involved in initiating morphine-induced changes in gene expression, but not in their maintenance.

Graphic abstract

Model of Polycomb repressive complex 2 (PRC2) and H3K27me3 alterations induced by chronic morphine exposure. Morphine induces H3K27me3 enrichment at promoters of genes encoding core members of the PRC2 complex and is associated with their transcriptional downregulation.

AICAR Stimulates the Pluripotency Transcriptional Complex in Embryonic Stem Cells Mediated by PI3K, GSK3β, and β-Catenin

Pluripotent stem cells maintain the property of self-renewal and differentiate into all cell types under clear environments. Though the gene regulatory mechanism for pluripotency has been investigated in recent years, it is still not completely understood. Here, we show several signaling pathways involved in the maintenance of pluripotency. To investigate whether AMPK is involved in maintaining the pluripotency in mouse embryonic stem cells (mESCs) and elucidating the possible molecular mechanisms, implicated D3 and R1/E mESC lines were used in this study. Cells were cultured in the absence or presence of LIF and treated with 1 mM and 0.5 mM 5-aminoimidazole-4-carboxamide-1-β-d-ribofuranoside (AICAR), 2 mM metformin, compound C, and the PI3K inhibitor LY294002 for 24, 72, and 120 h. The levels of Nanog, Oct3/4, and REX1 and Brachyury, Notch2, and Gata4 mRNAs and Nanog or OCT3/4 protein levels were analyzed. Alkaline phosphatase and the cellular cycle were determined. The pGSK3β, GSK3β, p-β-catenin, and β-catenin protein levels were also investigated. We found that AMPK activators such as AICAR and metformin increase mRNA expression of pluripotency markers and decrease mRNA expression of differentiation markers in R1/E and D3 ES cells. AICAR increases phosphatase activity and arrests the cellular cycle in the G1 phase in these cells. We describe that AICAR effects were mediated by AMPK activation using a chemical inhibitor or by silencing this gene. AICAR effects were also mediated by PI3K, GSK3β, and β-catenin in R1/E ES cells. According to our findings, we provide a mechanism by which AICAR increases and maintains a pluripotency state through enhanced Nanog expression, involving AMPK/PI3K and p-GSK3β Ser21/9 pathways backing up the AICAR function as a potential target for this drug controlling pluripotency. The highlights of this study are that AICAR (5-aminoimidazole-4-carboxamied-1-b-riboside), an AMP protein kinase (AMPK) activator, blocks the ESC differentiation and AMPK is a key enzyme for pluripotency and shows valuable data to clarify the molecular pluripotency mechanism.