Single cell lineage analysis of mouse embryonic stem cells at the exit from pluripotency

TOBF1 modulates mouse embryonic stem cell fate through regulating alternative splicing of pluripotency genes

Embryonic stem cells (ESCs) can undergo lineage-specific differentiation, giving rise to different cell types that constitute an organism. Although roles of transcription factors and chromatin modifiers in these cells have been described, how the alternative splicing (AS) machinery regulates their expression has not been sufficiently explored. Here, we show that the long non-coding RNA (lncRNA)-associated protein TOBF1 modulates the AS of transcripts necessary for maintaining stem cell identity in mouse ESCs. Among the genes affected is serine/arginine splicing factor 1 (SRSF1), whose AS leads to global changes in splicing and expression of a large number of downstream genes involved in the maintenance of ESC pluripotency. By overlaying information derived from TOBF1 chromatin occupancy, the distribution of its pluripotency-associated OCT-SOX binding motifs, and transcripts undergoing differential expression and AS upon its knockout, we describe local nuclear territories where these distinct events converge. Collectively, these contribute to the maintenance of mouse ESC identity.

Development of a novel embryonic germline gene-related prognostic model of lung adenocarcinoma

Background

Emerging evidence implicates the correlation of embryonic germline genes with the tumor progress and patient's outcome. However, the prognostic value of these genes in lung adenocarcinoma (LUAD) has not been fully studied. Here we systematically evaluated this issue, and constructed a novel signature and a nomogram associated with embryonic germline genes for predicting the outcomes of lung adenocarcinoma.

Methods

The LUAD cohorts retrieved from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database were used as training set and testing set, respectively. The embryonic germline genes were downloaded from the website https://venn.lodder.dev. Then, the differentially expressed embryonic germline genes (DEGGs) between the tumor and normal samples were identified by limma package. The functional enrichment and pathway analyses were also performed by clusterProfiler package. The prognostic model was constructed by the least absolute shrinkage and selection operator (LASSO)-Cox regression method. Survival and Receiver Operating Characteristic (ROC) analyses were performed to validate the model using training set and four testing GEO datasets. Finally, a prognostic nomogram based on the signature genes was constructed using multivariate regression method.

Results

Among the identified 269 DEGGs, 249 were up-regulated and 20 were down-regulated. GO and KEGG analyses revealed that these DEGGs were mainly enriched in the process of cell proliferation and DNA damage repair. Then, 103 DEGGs with prognostic value were identified by univariate Cox regression and further filtered by LASSO method. The resulting sixteen DEGGs were included in step multivariate Cox regression and an eleven embryonic germline gene related signature (EGRS) was constructed. The model could robustly stratify the LUAD patients into high-risk and low-risk groups in both training and testing sets, and low-risk patients had much better outcomes. The multi-ROC analysis also showed that the EGRS model had the best predictive efficacy compared with other common clinicopathological factors. The EGRS model also showed robust predictive ability in four independent external datasets, and the area under curve (AUC) was 0.726 (GSE30219), 0.764 (GSE50081), 0.657 (GSE37745) and 0.668 (GSE72094). More importantly, the expression level of some genes in EGRS has a significant correlation with the progression of LUAD clinicopathology, suggesting these genes might play an important role in the progression of LUAD. Finally, based on EGRS genes, we built and calibrated a nomogram for conveniently evaluating patients' outcomes.

Nanog regulates Pou3f1 expression at the exit from pluripotency during gastrulation

Pluripotency is regulated by a network of transcription factors that maintain early embryonic cells in an undifferentiated state while allowing them to proliferate. NANOG is a critical factor for maintaining pluripotency and its role in primordial germ cell differentiation has been well described. However, Nanog is expressed during gastrulation across all the posterior epiblast, and only later in development is its expression restricted to primordial germ cells. In this work, we unveiled a previously unknown mechanism by which Nanog specifically represses genes involved in anterior epiblast lineage. Analysis of transcriptional data from both embryonic stem cells and gastrulating mouse embryos revealed Pou3f1 expression to be negatively correlated with that of Nanog during the early stages of differentiation. We have functionally demonstrated Pou3f1 to be a direct target of NANOG by using a dual transgene system for the controlled expression of Nanog Use of Nanog null ES cells further demonstrated a role for Nanog in repressing a subset of anterior neural genes. Deletion of a NANOG binding site (BS) located nine kilobases downstream of the transcription start site of Pou3f1 revealed this BS to have a specific role in the regionalization of the expression of this gene in the embryo. Our results indicate an active role of Nanog inhibiting neural regulatory networks by repressing Pou3f1 at the onset of gastrulation.This article has an associated First Person interview with the joint first authors of the paper.

Developmental Xist induction is mediated by enhanced splicing

X-inactive-specific transcript (Xist) is a long noncoding RNA (lncRNA) essential for inactivating one of the two X chromosomes in mammalian females. Random X chromosome inactivation is mediated by Xist RNA expressed from the inactive X chromosome. We found that Xist RNA is unspliced in naïve embryonic stem (ES) cells. Upon differentiation, Xist splicing becomes efficient across all exons independent of transcription, suggesting interdependent or coordinated removal of Xist introns. In female cells with mutated polypyrimidine tract binding protein 1 (Ptbp1), differentiation fails to substantially upregulate mature Xist RNA because of a defect in Xist splicing. We further found both Xist129 and XistCAS RNA are unspliced in Mus musculus 129SvJ/Mus castaneous (CAS) hybrid female ES cells. Upon differentiation, Xist129 exhibits a higher splicing efficiency than XistCAS, likely contributing to preferential inhibition of the X129 chromosome. Single cell analysis shows that the allelic choice of Xist splicing is linked to the inactive X chromosome. We conclude post-transcriptional control of Xist RNA splicing is an essential regulatory step of Xist induction. Our studies shed light on the developmental roles of splicing for nuclear-retained Xist lncRNA and suggest inefficient Xist splicing is an additional fail-safe mechanism to prevent Xist activity in ES cells.

Histone Acetyltransferase KAT2A Stabilizes Pluripotency with Control of Transcriptional Heterogeneity

Cell fate transitions in mammalian stem cell systems have often been associated with transcriptional heterogeneity; however, existing data have failed to establish a functional or mechanistic link between the two phenomena. Experiments in unicellular organisms support the notion that transcriptional heterogeneity can be used to facilitate adaptability to environmental changes and have identified conserved chromatin‐associated factors that modulate levels of transcriptional noise. Herein, we show destabilization of pluripotency‐associated gene regulatory networks through increased transcriptional heterogeneity of mouse embryonic stem cells in which paradigmatic histone acetyl‐transferase, and candidate noise modulator, Kat2a (yeast orthologue Gcn5), have been inhibited. Functionally, network destabilization associates with reduced pluripotency and accelerated mesendodermal differentiation, with increased probability of transitions into lineage commitment. Thus, we show evidence of a relationship between transcriptional heterogeneity and cell fate transitions through manipulation of the histone acetylation landscape of mouse embryonic stem cells, suggesting a general principle that could be exploited in other normal and malignant stem cell fate transitions. stem cells2018;36:1828–11

Network Features and Dynamical Landscape of Naive and Primed Pluripotency

Although the broad and unique differentiation potential of pluripotent stem cells relies on a complex transcriptional network centered around Oct4, Sox2, and Nanog, two well-distinct pluripotent states, called "naive" and "primed", have been described in vitro and markedly differ in their developmental potential, their expression profiles, their signaling requirements, and their reciprocal conversion. Aiming to determine the key features that segregate and coordinate these two states, data-driven optimization of network models is performed to identify relevant parameter regimes and reduce network complexity to its core structure. Decision dynamics of optimized networks is characterized by signal-dependent multistability and strongly asymmetric transitions among naive, primed, and nonpluripotent states. Further model perturbation and reduction approaches reveal that such a dynamical landscape of pluripotency involves a functional partitioning of the regulatory network. Specifically, two overlapping positive feedback modules, Klf4/Esrrb/Nanog and Oct4/Nanog, stabilize the naive or the primed state, respectively. In turn, their incoherent feedforward and negative feedback coupling mediated by the Erk/Gsk3 module is critical for robust segregation and sequential progression between naive and primed states before irreversible exit from pluripotency.

Dynamics of lineage commitment revealed by single-cell transcriptomics of differentiating embryonic stem cells

Gene expression heterogeneity in the pluripotent state of mouse embryonic stem cells (mESCs) has been increasingly well-characterized. In contrast, exit from pluripotency and lineage commitment have not been studied systematically at the single-cell level. Here we measure the gene expression dynamics of retinoic acid driven mESC differentiation from pluripotency to lineage commitment, using an unbiased single-cell transcriptomics approach. We find that the exit from pluripotency marks the start of a lineage transition as well as a transient phase of increased susceptibility to lineage specifying signals. Our study reveals several transcriptional signatures of this phase, including a sharp increase of gene expression variability and sequential expression of two classes of transcriptional regulators. In summary, we provide a comprehensive analysis of the exit from pluripotency and lineage commitment at the single cell level, a potential stepping stone to improved lineage manipulation through timing of differentiation cues.

Commitment to different fates by differentiating pluripotent cells depends upon integration of external and internal signals. Here the authors analyse the entry of mouse embryonic stem cells into retinoic acid-mediated differentiation using single cell transcriptomics with high temporal resolution.

Zeb2 Regulates Cell Fate at the Exit from Epiblast State in Mouse Embryonic Stem Cells

In human embryonic stem cells (ESCs) the transcription factor Zeb2 regulates neuroectoderm versus mesendoderm formation, but it is unclear how Zeb2 affects the global transcriptional regulatory network in these cell‐fate decisions. We generated Zeb2 knockout (KO) mouse ESCs, subjected them as embryoid bodies (EBs) to neural and general differentiation and carried out temporal RNA‐sequencing (RNA‐seq) and reduced representation bisulfite sequencing (RRBS) analysis in neural differentiation. This shows that Zeb2 acts preferentially as a transcriptional repressor associated with developmental progression and that Zeb2 KO ESCs can exit from their naïve state. However, most cells in these EBs stall in an early epiblast‐like state and are impaired in both neural and mesendodermal differentiation. Genes involved in pluripotency, epithelial‐to‐mesenchymal transition (EMT), and DNA‐(de)methylation, including Tet1, are deregulated in the absence of Zeb2. The observed elevated Tet1 levels in the mutant cells and the knowledge of previously mapped Tet1‐binding sites correlate with loss‐of‐methylation in neural‐stimulating conditions, however, after the cells initially acquired the correct DNA‐methyl marks. Interestingly, cells from such Zeb2 KO EBs maintain the ability to re‐adapt to 2i + LIF conditions even after prolonged differentiation, while knockdown of Tet1 partially rescues their impaired differentiation. Hence, in addition to its role in EMT, Zeb2 is critical in ESCs for exit from the epiblast state, and links the pluripotency network and DNA‐methylation with irreversible commitment to differentiation. Stem Cells2017;35:611–625

Transition states and cell fate decisions in epigenetic landscapes

Waddington's epigenetic landscape is an abstract metaphor frequently used to represent the relationship between gene activity and cell fates during development. Over the past few years, it has become a useful framework for interpreting results from single-cell transcriptomics experiments. It has led to the proposal that, during fate transitions, cells experience smooth, continuous progressions of global transcriptional activity, which can be captured by (pseudo)temporal dynamics. Here, focusing strictly on the fate decision events, we suggest an alternative view: that fate transitions occur in a discontinuous, stochastic manner whereby signals modulate the probability of the transition events.

Visualization of the Epiblast and Visceral Endodermal Cells Using Fgf5-P2A-Venus BAC Transgenic Mice and Epiblast Stem Cells

Fibroblast growth factor 5 (Fgf5) has been widely used as a marker for the epiblast in the postimplantation embryo and epiblast stem cells (mEpiSCs) in the mouse, making it valuable for study of differentiation of various tissues and epiblast cells in vivo and in vitro. Here, we report for the first time the generation of Fgf5-P2A-Venus BAC transgenic (Tg) mice and show that the BAC Tg can recapitulate endogenous Fgf5 expression in epiblast and visceral endodermal cells of E6.5 and 7.5 embryos. We also show that Fgf5-P2A-Venus BAC Tg mEpiSCs in the undifferentiated state expressed abundant Venus, and upon reprogramming into naïve state, Venus was suppressed. Furthermore, while most Tg mEpiSCs expressed Venus abundantly, surprisingly the Tg mEpiSCs contained a minor subpopulation of Venus-negative cells that were capable of conversion to Venus-positive cells, indicating that even Fgf5 expression shows dynamic heterogeneity in mEpiSCs. Taken together, Fgf5-P2A-Venus BAC Tg mice and mEpiSCs generated in this study will be useful for developmental biology as well as stem cell biology research.

Inhibiting Smad2/3 signaling in pluripotent mouse embryonic stem cells enhances endoderm formation by increasing transcriptional priming of lineage-specifying target genes

BACKGROUND

Pluripotent embryonic stem cells (ESCs) offer great potential for regenerative medicine. However, efficient in vitro generation of specific desired cell types is still a challenge. We previously established that Smad2/3 signaling, essential for endoderm formation, regulates target gene expression by counteracting epigenetic repression mediated by Polycomb Repressive Complex 2 (PRC2). Although this mechanism has been demonstrated during differentiation and reprogramming, little is known of its role in pluripotent cells.

RESULTS

Chromatin immunoprecipitation-deep sequencing of undifferentiated mouse ESCs inhibited for Smad2/3 signaling identified Prdm14, important for protecting pluripotency, as a target gene. Although Prdm14 accumulates the normally repressive PRC2 deposited histone modification H3K27me3 under these conditions, surprisingly, expression increases. Analysis indicates that increased H3K27me3 leads to increased binding of PRC2 accessory component Jarid2 and recruitment of RNA polymerase II. Similar increases were found at the Nodal endoderm target gene Eomes but it remained unexpressed in pluripotent cells as normal. Upon differentiation, however, Eomes expression was significantly higher than in cells that had not been inhibited for signaling before differentiation. In addition, endoderm formation was markedly increased.

CONCLUSIONS

Blocking Smad2/3 signaling in pluripotent stem cells results in epigenetic changes that enhance the capacity for endoderm differentiation. Developmental Dynamics 245:807-815, 2016. © 2016 Wiley Periodicals, Inc.

BMP-SMAD Signaling Regulates Lineage Priming, but Is Dispensable for Self-Renewal in Mouse Embryonic Stem Cells

Naive mouse embryonic stem cells (mESCs) are in a metastable state and fluctuate between inner cell mass- and epiblast-like phenotypes. Here, we show transient activation of the BMP-SMAD signaling pathway in mESCs containing a BMP-SMAD responsive reporter transgene. Activation of the BMP-SMAD reporter transgene in naive mESCs correlated with lower levels of genomic DNA methylation, high expression of 5-methylcytosine hydroxylases Tet1/2 and low levels of DNA methyltransferases Dnmt3a/b. Moreover, naive mESCs, in which the BMP-SMAD reporter transgene was activated, showed higher resistance to differentiation. Using double Smad1;Smad5 knockout mESCs, we showed that BMP-SMAD signaling is dispensable for self-renewal in both naive and ground state. These mutant mESCs were still pluripotent, but they exhibited higher levels of DNA methylation than their wild-type counterparts and had a higher propensity to differentiate. We showed that BMP-SMAD signaling modulates lineage priming in mESCs, by transiently regulating the enzymatic machinery responsible for DNA methylation.

•

BMP-SMAD signaling in mESCs is more prominent in naive than ground state

•

BMP-SMAD signaling is dispensable for pluripotency in mESCs

•

BMP-SMAD signaling facilitates lineage priming in mESCs

•

BMP-SMAD signaling regulates Dnmt3b and hence levels of DNA methylation

Heterogeneous lineage marker expression in naive embryonic stem cells is mostly due to spontaneous differentiation

Populations of cultured mouse embryonic stem cells (ESCs) exhibit a subfraction of cells expressing uncharacteristically low levels of pluripotency markers such as Nanog. Yet, the extent to which individual Nanog-negative cells are differentiated, both from ESCs and from each other, remains unclear. Here, we show the transcriptome of Nanog-negative cells exhibits expression of classes of genes associated with differentiation that are not yet active in cells exposed to differentiation conditions for one day. Long non-coding RNAs, however, exhibit more changes in expression in the one-day-differentiated cells than in Nanog-negative cells. These results are consistent with the concept that Nanog-negative cells may contain subpopulations of both lineage-primed and differentiated cells. Single cell analysis showed that Nanog-negative cells display substantial and coherent heterogeneity in lineage marker expression in progressively nested subsets of cells exhibiting low levels of Nanog, then low levels of Oct4, and then a set of lineage markers, which express intensely in a small subset of these more differentiated cells. Our results suggest that the observed enrichment of lineage-specific marker gene expression in Nanog-negative cells is associated with spontaneous differentiation of a subset of these cells rather than the more random expression that may be associated with reversible lineage priming.

From blastocyst to gastrula: gene regulatory networks of embryonic stem cells and early mouse embryogenesis

To date, many regulatory genes and signalling events coordinating mammalian development from blastocyst to gastrulation stages have been identified by mutational analyses and reverse-genetic approaches, typically on a gene-by-gene basis. More recent studies have applied bioinformatic approaches to generate regulatory network models of gene interactions on a genome-wide scale. Such models have provided insights into the gene networks regulating pluripotency in embryonic and epiblast stem cells, as well as cell-lineage determination in vivo. Here, we review how regulatory networks constructed for different stem cell types relate to corresponding networks in vivo and provide insights into understanding the molecular regulation of the blastocyst-gastrula transition.

The nature of embryonic stem cells

Mouse embryonic stem (ES) cells perpetuate in vitro the broad developmental potential of naïve founder cells in the preimplantation embryo. ES cells self-renew relentlessly in culture but can reenter embryonic development seamlessly, differentiating on schedule to form all elements of the fetus. Here we review the properties of these remarkable cells. Arising from the stability, homogeneity, and equipotency of ES cells, we consider the concept of a pluripotent ground state. We evaluate the authenticity of ES cells in relation to cells in the embryo and examine their utility for dissecting mechanisms that confer pluripotency and that execute fate choice. We summarize current knowledge of the transcription factor circuitry that governs the ES cell state and discuss the opportunity to expose molecular logic further through iterative computational modeling and experimentation. Finally, we present a perspective on unresolved questions, including the challenge of deriving ground state pluripotent stem cells from non-rodent species.

β-Catenin Regulates Primitive Streak Induction through Collaborative Interactions with SMAD2/SMAD3 and OCT4

Canonical Wnt and Nodal signaling are both required for induction of the primitive streak (PS), which guides organization of the early embryo. The Wnt effector β-catenin is thought to function in these early lineage specification decisions via transcriptional activation of Nodal signaling. Here, we demonstrate a broader role for β-catenin in PS formation by analyzing its genome-wide binding in a human embryonic stem cell model of PS induction. β-catenin occupies regulatory regions in numerous PS and neural crest genes, and direct interactions between β-catenin and the Nodal effectors SMAD2/SMAD3 are required at these regions for PS gene activation. Furthermore, OCT4 binding in proximity to these sites is likewise required for PS induction, suggesting a collaborative interaction between β-catenin and OCT4. Induction of neural crest genes by β-catenin is repressed by SMAD2/SMAD3, ensuring proper lineage specification. This study provides mechanistic insight into how Wnt signaling controls early cell lineage decisions.

Brachyury cooperates with Wnt/β-catenin signalling to elicit primitive-streak-like behaviour in differentiating mouse embryonic stem cells

Background

The formation of the primitive streak is the first visible sign of gastrulation, the process by which the three germ layers are formed from a single epithelium during early development. Embryonic stem cells (ESCs) provide a good system for understanding the molecular and cellular events associated with these processes. Previous work, both in embryos and in culture, has shown how converging signals from both nodal/TGFβR and Wnt/β-catenin signalling pathways specify cells to adopt a primitive-streak-like fate and direct them to undertake an epithelial-to-mesenchymal transition (EMT). However, many of these approaches have relied on genetic analyses without taking into account the temporal progression of events within single cells. In addition, it is still unclear to what extent events in the embryo are able to be reproduced in culture.

Results

Here, we combine flow cytometry and a quantitative live single-cell imaging approach to demonstrate how the controlled differentiation of mouse ESCs towards a primitive streak fate in culture results in cells displaying many of the characteristics observed during early mouse development including transient brachyury expression, EMT and increased motility. We also find that the EMT initiates the process, and this is both fuelled and terminated by the action of brachyury, whose expression is dependent on the EMT and β-catenin activity.

Conclusions

As a consequence of our analysis, we propose that a major output of brachyury expression is in controlling the velocity of the cells that are transiting out of the primitive streak.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-014-0063-7) contains supplementary material, which is available to authorized users.