Inductive and Selective Effects of GSK3 and MEK Inhibition on Nanog Heterogeneity in Embryonic Stem Cells

Gene-specific RNA homeostasis revealed by perturbation of coactivator complexes

Transcript buffering entails the reciprocal modulation of mRNA synthesis and degradation rates to maintain stable RNA levels under varying cellular conditions. Current research supports a global, non-sequence-specific connection between mRNA synthesis and degradation, but the underlying mechanisms are still unclear. In this study, we investigated changes in RNA metabolism following acute depletion of TIP60/KAT5, the acetyltransferase subunit of the NuA4 transcriptional coactivator complex, in mouse embryonic stem cells. By combining RNA sequencing of nuclear, cytoplasmic, and newly synthesised transcript fractions with biophysical modelling, we demonstrate that TIP60 predominantly enhances transcription of numerous genes, while a smaller set of genes undergoes TIP60-dependent transcriptional repression. Surprisingly, transcription changes caused by TIP60 depletion were offset by corresponding changes in RNA nuclear export and cytoplasmic stability, indicating gene-specific buffering mechanisms. Similarly, disruption of the unrelated ATAC coactivator complex also resulted in gene-specific transcript buffering. These findings reveal that transcript buffering functions at a gene-specific level and suggest that cells dynamically adjust RNA splicing, export, and degradation in response to individual RNA synthesis alterations, thereby sustaining cellular homeostasis.

Combinatorial quantification of 5mC and 5hmC at individual CpG dyads and the transcriptome in single cells reveals modulators of DNA methylation maintenance fidelity

Inheritance of 5-methylcytosine from one cell generation to the next by DNA methyltransferase 1 (DNMT1) plays a key role in regulating cellular identity. While recent work has shown that the activity of DNMT1 is imprecise, it remains unclear how the fidelity of DNMT1 is tuned in different genomic and cell state contexts. Here we describe Dyad-seq, a method to quantify the genome-wide methylation status of cytosines at the resolution of individual CpG dinucleotides to find that the fidelity of DNMT1-mediated maintenance methylation is related to the local density of DNA methylation and the landscape of histone modifications. To gain deeper insights into methylation/demethylation turnover dynamics, we first extended Dyad-seq to quantify all combinations of 5-methylcytosine and 5-hydroxymethylcytosine at individual CpG dyads. Next, to understand how cell state transitions impact maintenance methylation, we scaled the method down to jointly profile genome-wide methylation levels, maintenance methylation fidelity and the transcriptome from single cells (scDyad&T-seq). Using scDyad&T-seq, we demonstrate that, while distinct cell states can substantially impact the activity of the maintenance methylation machinery, locally there exists an intrinsic relationship between DNA methylation density, histone modifications and DNMT1-mediated maintenance methylation fidelity that is independent of cell state.

The BAF chromatin remodeler synergizes with RNA polymerase II and transcription factors to evict nucleosomes

Chromatin accessibility is a hallmark of active transcription and entails ATP-dependent nucleosome remodeling, which is carried out by complexes such as Brahma-associated factor (BAF). However, the mechanistic links between transcription, nucleosome remodeling and chromatin accessibility are unclear. Here, we used a chemical-genetic approach coupled with time-resolved chromatin profiling to dissect the interplay between RNA Polymerase II (RNAPII), BAF and DNA-sequence-specific transcription factors in mouse embryonic stem cells. We show that BAF dynamically unwraps and evicts nucleosomes at accessible chromatin regions, while RNAPII promoter-proximal pausing stabilizes BAF chromatin occupancy and enhances ATP-dependent nucleosome eviction by BAF. We find that although RNAPII and BAF dynamically probe both transcriptionally active and Polycomb-repressed genomic regions, pluripotency transcription factor chromatin binding confers locus specificity for productive chromatin remodeling and nucleosome eviction by BAF. Our study suggests a paradigm for how functional synergy between dynamically acting chromatin factors regulates locus-specific nucleosome organization and chromatin accessibility.

An improved rhythmicity analysis method using Gaussian Processes detects cell-density dependent circadian oscillations in stem cells

MOTIVATION

Detecting oscillations in time series remains a challenging problem even after decades of research. In chronobiology, rhythms (for instance in gene expression, eclosion, egg-laying, and feeding) tend to be low amplitude, display large variations amongst replicates, and often exhibit varying peak-to-peak distances (non-stationarity). Most currently available rhythm detection methods are not specifically designed to handle such datasets, and are also limited by their use of P-values in detecting oscillations.

RESULTS

We introduce a new method, ODeGP (Oscillation Detection using Gaussian Processes), which combines Gaussian Process regression and Bayesian inference to incorporate measurement errors, non-uniformly sampled data, and a recently developed non-stationary kernel to improve detection of oscillations. By using Bayes factors, ODeGP models both the null (non-rhythmic) and the alternative (rhythmic) hypotheses, thus providing an advantage over P-values. Using synthetic datasets, we first demonstrate that ODeGP almost always outperforms eight commonly used methods in detecting stationary as well as non-stationary symmetric oscillations. Next, by analyzing existing qPCR datasets, we demonstrate that our method is more sensitive compared to the existing methods at detecting weak and noisy oscillations. Finally, we generate new qPCR data on mouse embryonic stem cells. Surprisingly, we discover using ODeGP that increasing cell-density results in rapid generation of oscillations in the Bmal1 gene, thus highlighting our method's ability to discover unexpected and new patterns. In its current implementation, ODeGP is meant only for analyzing single or a few time-trajectories, not genome-wide datasets.

AVAILABILITY AND IMPLEMENTATION

ODeGP is available at https://github.com/Shaonlab/ODeGP.

Embryonic stem cell ERK, AKT, plus STAT3 response dynamics combinatorics are heterogeneous but NANOG state independent

Signaling is central in cell fate regulation, and relevant information is encoded in its activity over time (i.e., dynamics). However, simultaneous dynamics quantification of several pathways in single mammalian stem cells has not yet been accomplished. Here we generate mouse embryonic stem cell (ESC) lines simultaneously expressing fluorescent reporters for ERK, AKT, and STAT3 signaling activity, which all control pluripotency. We quantify their single-cell dynamics combinations in response to different self-renewal stimuli and find striking heterogeneity for all pathways, some dependent on cell cycle but not pluripotency states, even in ESC populations currently assumed to be highly homogeneous. Pathways are mostly independently regulated, but some context-dependent correlations exist. These quantifications reveal surprising single-cell heterogeneity in the important cell fate control layer of signaling dynamics combinations and raise fundamental questions about the role of signaling in (stem) cell fate control.

Combinatorial quantification of 5mC and 5hmC at individual CpG dyads and the transcriptome in single cells reveals modulators of DNA methylation maintenance fidelity

Transmission of 5-methylcytosine (5mC) from one cell generation to the next plays a key role in regulating cellular identity in mammalian development and diseases. While recent work has shown that the activity of DNMT1, the protein responsible for the stable inheritance of 5mC from mother to daughter cells, is imprecise; it remains unclear how the fidelity of DNMT1 is tuned in different genomic and cell state contexts. Here we describe Dyad-seq, a method that combines enzymatic detection of modified cytosines with nucleobase conversion techniques to quantify the genome-wide methylation status of cytosines at the resolution of individual CpG dinucleotides. We find that the fidelity of DNMT1-mediated maintenance methylation is directly related to the local density of DNA methylation, and for genomic regions that are lowly methylated, histone modifications can dramatically alter the maintenance methylation activity. Further, to gain deeper insights into the methylation and demethylation turnover dynamics, we extended Dyad-seq to quantify all combinations of 5mC and 5-hydroxymethylcytosine (5hmC) at individual CpG dyads to show that TET proteins preferentially hydroxymethylate only one of the two 5mC sites in a symmetrically methylated CpG dyad rather than sequentially convert both 5mC to 5hmC. To understand how cell state transitions impact DNMT1-mediated maintenance methylation, we scaled the method down and combined it with the measurement of mRNA to simultaneously quantify genome-wide methylation levels, maintenance methylation fidelity and the transcriptome from the same cell (scDyad&T-seq). Applying scDyad&T-seq to mouse embryonic stem cells transitioning from serum to 2i conditions, we observe dramatic and heterogenous demethylation and the emergence of transcriptionally distinct subpopulations that are closely linked to the cell-to-cell variability in loss of DNMT1-mediated maintenance methylation activity, with regions of the genome that escape 5mC reprogramming retaining high levels of maintenance methylation fidelity. Overall, our results demonstrate that while distinct cell states can substantially impact the genome-wide activity of the DNA methylation maintenance machinery, locally there exists an intrinsic relationship between DNA methylation density, histone modifications and DNMT1-mediated maintenance methylation fidelity that is independent of cell state.

An improved rhythmicity analysis method using Gaussian Processes detects cell-density dependent circadian oscillations in stem cells

Detecting oscillations in time series remains a challenging problem even after decades of research. In chronobiology, rhythms in time series (for instance gene expression, eclosion, egg-laying and feeding) datasets tend to be low amplitude, display large variations amongst replicates, and often exhibit varying peak-to-peak distances (non-stationarity). Most currently available rhythm detection methods are not specifically designed to handle such datasets. Here we introduce a new method, ODeGP ( O scillation De tection using G aussian P rocesses), which combines Gaussian Process (GP) regression with Bayesian inference to provide a flexible approach to the problem. Besides naturally incorporating measurement errors and non-uniformly sampled data, ODeGP uses a recently developed kernel to improve detection of non-stationary waveforms. An additional advantage is that by using Bayes factors instead of p-values, ODeGP models both the null (non-rhythmic) and the alternative (rhythmic) hypotheses. Using a variety of synthetic datasets we first demonstrate that ODeGP almost always outperforms eight commonly used methods in detecting stationary as well as non-stationary oscillations. Next, on analyzing existing qPCR datasets that exhibit low amplitude and noisy oscillations, we demonstrate that our method is more sensitive compared to the existing methods at detecting weak oscillations. Finally, we generate new qPCR time-series datasets on pluripotent mouse embryonic stem cells, which are expected to exhibit no oscillations of the core circadian clock genes. Surprisingly, we discover using ODeGP that increasing cell density can result in the rapid generation of oscillations in the Bmal1 gene, thus highlighting our method’s ability to discover unexpected patterns. In its current implementation, ODeGP (available as an R package) is meant only for analyzing single or a few time-trajectories, not genome-wide datasets.

RNA Polymerase II, the BAF remodeler and transcription factors synergize to evict nucleosomes

Chromatin accessibility is a hallmark of active transcription and requires ATP-dependent nucleosome remodeling by Brahma-Associated Factor (BAF). However, the mechanistic link between transcription, nucleosome remodeling, and chromatin accessibility is unclear. Here, we used a chemical-genetic approach to dissect the interplay between RNA Polymerase II (RNAPII), BAF, and DNA-sequence-specific transcription factors (TFs) in mouse embryonic stem cells. By time-resolved chromatin profiling with acute transcription block at distinct stages, we show that RNAPII promoter-proximal pausing stabilizes BAF chromatin occupancy and enhances nucleosome eviction by BAF. We find that RNAPII and BAF probe both transcriptionally active and Polycomb-repressed genomic regions and provide evidence that TFs capture transient site exposure due to nucleosome unwrapping by BAF to confer locus specificity for persistent chromatin remodeling. Our study reveals the mechanistic basis of cell-type-specific chromatin accessibility. We propose a new paradigm for how functional synergy between dynamically acting chromatin factors regulates nucleosome organization.

H3K9 tri-methylation at Nanog times differentiation commitment and enables the acquisition of primitive endoderm fate

Mouse embryonic stem cells have an inherent propensity to explore gene regulatory states associated with either self-renewal or differentiation. This property depends on ERK, which downregulates pluripotency genes such as Nanog. Here, we aimed at identifying repressive histone modifications that would mark Nanog for inactivation in response to ERK activity. We found that the transcription factor ZFP57, which binds methylated DNA to nucleate heterochromatin, is recruited upstream of Nanog, within a region enriched for histone H3 lysine 9 tri-methylation (H3K9me3). Whereas before differentiation H3K9me3 at Nanog depends on ERK, in somatic cells it becomes independent of ERK. Moreover, the loss of H3K9me3 at Nanog, induced by deleting the region or by knocking out DNA methyltransferases or Zfp57, is associated with reduced heterogeneity of NANOG, delayed commitment into differentiation and impaired ability to acquire a primitive endoderm fate. Hence, a network axis centred on DNA methylation, ZFP57 and H3K9me3 links Nanog regulation to ERK activity for the timely establishment of new cell identities. We suggest that establishment of irreversible H3K9me3 at specific master regulators allows the acquisition of particular cell fates during differentiation.

Summary: A regulatory axis integrating ERK, ZFP57, DNA and H3K9 methylation underlies the transition of Nanog expression from heterogeneous and dynamic to irreversibly silenced, enabling differentiation commitment and primitive endoderm specification.

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.

Epiblast inducers capture mouse trophectoderm stem cells in vitro and pattern blastoids for implantation in utero

The embryo instructs the allocation of cell states to spatially regulate functions. In the blastocyst, patterning of trophoblast (TR) cells ensures successful implantation and placental development. Here, we defined an optimal set of molecules secreted by the epiblast (inducers) that captures in vitro stable, highly self-renewing mouse trophectoderm stem cells (TESCs) resembling the blastocyst stage. When exposed to suboptimal inducers, these stem cells fluctuate to form interconvertible subpopulations with reduced self-renewal and facilitated differentiation, resembling peri-implantation cells, known as TR stem cells (TSCs). TESCs have enhanced capacity to form blastoids that implant more efficiently in utero due to inducers maintaining not only local TR proliferation and self-renewal, but also WNT6/7B secretion that stimulates uterine decidualization. Overall, the epiblast maintains sustained growth and decidualization potential of abutting TR cells, while, as known, distancing imposed by the blastocyst cavity differentiates TR cells for uterus adhesion, thus patterning the essential functions of implantation.

Bioimaging approaches for quantification of individual cell behavior during cell fate decisions

Tracking individual cells has allowed a new understanding of cellular behavior in human health and disease by adding a dynamic component to the already complex heterogeneity of single cells. Technically, despite countless advances, numerous experimental variables can affect data collection and interpretation and need to be considered. In this review, we discuss the main technical aspects and biological findings in the analysis of the behavior of individual cells. We discuss the most relevant contributions provided by these approaches in clinically relevant human conditions like embryo development, stem cells biology, inflammation, cancer and microbiology, along with the cellular mechanisms and molecular pathways underlying these conditions. We also discuss the key technical aspects to be considered when planning and performing experiments involving the analysis of individual cells over long periods. Despite the challenges in automatic detection, features extraction and long-term tracking that need to be tackled, the potential impact of single-cell bioimaging is enormous in understanding the pathogenesis and development of new therapies in human pathophysiology.

Hypoblast Formation in Bovine Embryos Does Not Depend on NANOG

The role of the pluripotency factor NANOG during the second embryonic lineage differentiation has been studied extensively in mouse, although species-specific differences exist. To elucidate the role of NANOG in an alternative model organism, we knocked out NANOG in fibroblast cells and produced bovine NANOG-knockout (KO) embryos via somatic cell nuclear transfer (SCNT). At day 8, NANOG-KO blastocysts showed a decreased total cell number when compared to controls from SCNT (NT Ctrl). The pluripotency factors OCT4 and SOX2 as well as the hypoblast (HB) marker GATA6 were co-expressed in all cells of the inner cell mass (ICM) and, in contrast to mouse Nanog-KO, expression of the late HB marker SOX17 was still present. We blocked the MEK-pathway with a MEK 1/2 inhibitor, and control embryos showed an increase in NANOG positive cells, but SOX17 expressing HB precursor cells were still present. NANOG-KO together with MEK-inhibition was lethal before blastocyst stage, similarly to findings in mouse. Supplementation of exogenous FGF4 to NANOG-KO embryos did not change SOX17 expression in the ICM, unlike mouse Nanog-KO embryos, where missing SOX17 expression was completely rescued by FGF4. We conclude that NANOG mediated FGF/MEK signaling is not required for HB formation in the bovine embryo and that another—so far unknown—pathway regulates HB differentiation.

Fine-tuning Nanog expression heterogeneity in embryonic stem cells by regulating a Nanog transcript-specific microRNA

In embryonic stem cells (ESCs), the transcription factor Nanog maintains the stemness of ESCs despite exhibiting heterogeneous expression patterns under varied culture conditions. Efficient fine-tuning of Nanog expression heterogeneity could enable ESC proliferation and differentiation along specific lineages to be regulated. Herein, by employing a stochastic modeling approach, we show that Nanog expression heterogeneity can be controlled by modulating the regulatory features of a Nanog transcript-specific microRNA, mir-296. We demonstrate how and why the extent of origin-dependent fluctuations in Nanog expression level can be altered by varying either the binding efficiency of the microRNA-mRNA complex or the expression level of mir-296. Moreover, our model makes experimentally feasible and insightful predictions to maneuver Nanog expression heterogeneity explicitly to achieve cell-type-specific differentiation of ESCs.

2i, or Not 2i: The Soliloquy of Nanog-Negative Mouse Embryonic Stem Cells

In this issue of Stem Cell Reports, Hastreiter et al. (2018) use continuous time-lapse imaging of mouse embryonic stem cells to investigate how the inhibition of GSK3b and MEK/ERK (2i) leads to homogeneous expression of the transcription factor Nanog. They show that both induction of Nanog expression and selection against cells expressing low levels of Nanog contribute to the homogeneous appearance of 2i cultures.

Restoration of primary cilia in obese adipose-derived mesenchymal stem cells by inhibiting Aurora A or extracellular signal-regulated kinase

BACKGROUND

Obesity impairs a variety of cell types including adipose-derived mesenchymal stem cells (ASCs). ASCs are indispensable for tissue homeostasis/repair, immunomodulation, and cell renewal. It has been demonstrated that obese ASCs are defective in differentiation, motility, immunomodulation, and replication. We have recently reported that some of these defects are linked to impaired primary cilia, which are unable to properly convey and coordinate a variety of signaling pathways. We hypothesized that the rescue of the primary cilium in obese ASCs would restore their functional properties.

METHODS

Obese ASCs derived from subcutaneous and visceral adipose tissues were treated with a specific inhibitor against Aurora A or with an inhibitor against extracellular signal-regulated kinase 1/2 (Erk1/2). Multiple molecular and cellular assays were performed to analyze the altered functionalities and their involved pathways.

RESULTS

The treatment with low doses of these inhibitors extended the length of the primary cilium, restored the invasion and migration potential, and improved the differentiation capacity of obese ASCs. Associated with enhanced differentiation ability, the cells displayed an increased expression of self-renewal/stemness-related genes like SOX2, OCT4, and NANOG, mediated by reduced active glycogen synthase kinase 3 β (GSK3β).

CONCLUSION

This work describes a novel phenomenon whereby the primary cilium of obese ASCs is rescuable by the low-dose inhibition of Aurora A or Erk1/2, restoring functional ASCs with increased stemness. These cells might be able to improve tissue homeostasis in obese patients and thereby ameliorate obesity-associated diseases. Additionally, these functionally restored obese ASCs could be useful for novel autologous mesenchymal stem cell-based therapies.

Extensive epigenetic and transcriptomic variability between genetically identical human B-lymphoblastoid cells with implications in pharmacogenomics research

Genotyped human B-lymphoblastoid cell lines (LCLs) are widely used models in mapping quantitative trait loci for chromatin features, gene expression, and drug response. The extent of genotype-independent functional genomic variability of the LCL model, although largely overlooked, may inform association study design. In this study, we use flow cytometry, chromatin immunoprecipitation sequencing and mRNA sequencing to study surface marker patterns, quantify genome-wide chromatin changes (H3K27ac) and transcriptome variability, respectively, among five isogenic LCLs derived from the same individual. Most of the studied LCLs were non-monoclonal and had mature B cell phenotypes. Strikingly, nearly one-fourth of active gene regulatory regions showed significantly variable H3K27ac levels, especially enhancers, among which several were classified as clustered enhancers. Large, contiguous genomic regions showed signs of coordinated activity change. Regulatory differences were mirrored by mRNA expression changes, preferentially affecting hundreds of genes involved in specialized cellular processes including immune and drug response pathways. Differential expression of DPYD, an enzyme involved in 5-fluorouracil (5-FU) catabolism, was associated with variable LCL growth inhibition mediated by 5-FU. The extent of genotype-independent functional genomic variability might highlight the need to revisit study design strategies for LCLs in pharmacogenomics.

Esrrb extinction triggers dismantling of naïve pluripotency and marks commitment to differentiation

Self-renewal of embryonic stem cells (ESCs) cultured in LIF/fetal calf serum (FCS) is incomplete with some cells initiating differentiation. While this is reflected in heterogeneous expression of naive pluripotency transcription factors (TFs), the link between TF heterogeneity and differentiation is not fully understood. Here, we purify ESCs with distinct TF expression levels from LIF/FCS cultures to uncover early events during commitment from naïve pluripotency. ESCs carrying fluorescent Nanog and Esrrb reporters show Esrrb downregulation only in Nanoglow cells. Independent Esrrb reporter lines demonstrate that Esrrbnegative ESCs cannot effectively self-renew. Upon Esrrb loss, pre-implantation pluripotency gene expression collapses. ChIP-Seq identifies different regulatory element classes that bind both OCT4 and NANOG in Esrrbpositive cells. Class I elements lose NANOG and OCT4 binding in Esrrbnegative ESCs and associate with genes expressed preferentially in naïve ESCs. In contrast, Class II elements retain OCT4 but not NANOG binding in ESRRB-negative cells and associate with more broadly expressed genes. Therefore, mechanistic differences in TF function act cumulatively to restrict potency during exit from naïve pluripotency.