Nr0b1 is a negative regulator of Zscan4c in mouse embryonic stem cells

The impact of retrotransposons on zygotic genome activation and the chromatin landscape of early embryos

In mammals, fertilization is followed by extensive reprogramming and reorganization of the chromatin accompanying the transcriptional activation of the embryo. This reprogramming results in blastomeres with the ability to give rise to all cell types and a complete organism, including extra-embryonic tissues, and is known as totipotency. Transcriptional activation occurs in a process known as zygotic genome activation (ZGA) and is tightly linked to the expression of transposable elements, including endogenous retroviruses (ERVs) such as endogenous retrovirus with leucine tRNA primer (ERVL). Recent studies discovered the importance of ERVs in this process, yet the race to decipher the network surrounding these elements is still ongoing, and the molecular mechanism behind their involvement remains a mystery. Amid a recent surge of studies reporting the discovery of various factors and pathways involved in the regulation of ERVs, this review provides an overview of the knowns and unknowns in the field, with a particular emphasis on the chromatin landscape and how ERVs shape preimplantation development in mammals. In so doing, we highlight recent discoveries that have advanced our understanding of how these elements are involved in transforming the quiescent zygote into the most powerful cell type in mammals.

Nr5a2 is dispensable for zygotic genome activation but essential for morula development

Early embryogenesis is driven by transcription factors (TFs) that first activate the zygotic genome and then specify the lineages constituting the blastocyst. Although the TFs specifying the blastocyst's lineages are well characterized, those playing earlier roles remain poorly defined. Using mouse models of the TF Nr5a2, we show that Nr5a2-/- embryos arrest at the early morula stage and exhibit altered lineage specification, frequent mitotic failure, and substantial chromosome segregation defects. Although NR5A2 plays a minor but measurable role during zygotic genome activation, it predominantly acts as a master regulator at the eight-cell stage, controlling expression of lineage-specifying TFs and genes involved in mitosis, telomere maintenance, and DNA repair. We conclude that NR5A2 coordinates proliferation, genome stability, and lineage specification to ensure correct morula development.

The role of Zfp352 in the regulation of transient expression of 2-cell specific genes in mouse embryonic stem cells

Mouse ES cell populations contain a minor sub-population that expresses genes specifically expressed in 2-cell stage embryos. This sub-population consists of 2-cell-gene labeled cells (2CLCs) generated by the transient activation of the 2-cell specific genes initiated by the master regulator, Dux. However, the mechanism regulating the transient expression remains largely unclear. Here we reported a novel function of Zfp352, one of the 2-cell specific genes, in regulating the 2CLC sub-population. Zfp352 encodes zinc-finger transcription factor belonging to the Klf family. Dux transiently activates Zfp352 after the activation of Zscan4c in a subset of the 2CLC subpopulation. Interestingly, in the reporter assay, the transcriptional activation of Zscan4c by Dux is strongly repressed by the co-expression of Zfp352. However, the knockout of Zfp352 resulted in the repression of a subset of the 2-cell-specific genes. These data suggest the dual roles of Zfp352 in regulating the transient activation of the 2-cell-specific genes.

Dex-Benchmark: datasets and code to evaluate algorithms for transcriptomics data analysis

Many tools and algorithms are available for analyzing transcriptomics data. These include algorithms for performing sequence alignment, data normalization and imputation, clustering, identifying differentially expressed genes, and performing gene set enrichment analysis. To make the best choice about which tools to use, objective benchmarks can be developed to compare the quality of different algorithms to extract biological knowledge maximally and accurately from these data. The Dexamethasone Benchmark (Dex-Benchmark) resource aims to fill this need by providing the community with datasets and code templates for benchmarking different gene expression analysis tools and algorithms. The resource provides access to a collection of curated RNA-seq, L1000, and ChIP-seq data from dexamethasone treatment as well as genetic perturbations of its known targets. In addition, the website provides Jupyter Notebooks that use these pre-processed curated datasets to demonstrate how to benchmark the different steps in gene expression analysis. By comparing two independent data sources and data types with some expected concordance, we can assess which tools and algorithms best recover such associations. To demonstrate the usefulness of the resource for discovering novel drug targets, we applied it to optimize data processing strategies for the chemical perturbations and CRISPR single gene knockouts from the L1000 transcriptomics data from the Library of Integrated Network Cellular Signatures (LINCS) program, with a focus on understudied proteins from the Illuminating the Druggable Genome (IDG) program. Overall, the Dex-Benchmark resource can be utilized to assess the quality of transcriptomics and other related bioinformatics data analysis workflows. The resource is available from: https://maayanlab.github.io/dex-benchmark.

Loss of full-length DNA replication regulator Rif1 in two-cell embryos is associated with zygotic transcriptional activation

Rif1 regulates DNA replication timing and double-strand break repair, and its depletion induces transcriptional bursting of two-cell (2C) zygote-specific genes in mouse ES cells. However, how Rif1 regulates zygotic transcription is unclear. We show here that Rif1 depletion promotes the formation of a unique Zscan4 enhancer structure harboring both histone H3 lysine 27 acetylation (H3K27ac) and moderate levels of silencing chromatin mark H3K9me3. Curiously, another enhancer mark H3K4me1 is missing, whereas DNA methylation is still maintained in the structure, which spreads across gene bodies and neighboring regions within the Zscan4 gene cluster. We also found by function analyses of Rif1 domains in ES cells that ectopic expression of Rif1 lacking N-terminal domain results in upregulation of 2C transcripts. This appears to be caused by dominant negative inhibition of endogenous Rif1 protein localization at the nuclear periphery through formation of hetero-oligomers between the N-terminally truncated and endogenous forms. Strikingly, in murine 2C embryos, most of Rif1-derived polypeptides are expressed as truncated forms in soluble nuclear or cytosolic fraction and are likely nonfunctional. Toward the morula stage, the full-length form of Rif1 gradually increased. Our results suggest that the absence of the functional full-length Rif1 due to its instability or alternative splicing and potential inactivation of Rif1 through dominant inhibition by N-terminally truncated Rif1 polypeptides may be involved in 2C-specific transcription program.

nr0b1 (DAX1) loss of function in zebrafish causes hypothalamic defects via abnormal progenitor proliferation and differentiation

The nuclear receptor DAX-1 (encoded by the NR0B1 gene) is presented in the hypothalamic tissues in humans and other vertebrates. Human patients with NR0B1 mutations often have hypothalamic-pituitary defects, but the involvement of NR0B1 in hypothalamic development and function is not well understood. Here, we report the disruption of the nr0b1 gene in zebrafish causes abnormal expression of gonadotropins, a reduction in fertilization rate, and an increase in post-fasting food intake, which is indicative of abnormal hypothalamic functions. We find that loss of nr0b1 increases the number of prodynorphin (pdyn)-expressing neurons but decreases the number of pro-opiomelanocortin (pomcb)-expressing neurons in the zebrafish hypothalamic arcuate region (ARC). Further examination reveals that the proliferation of progenitor cells is reduced in the hypothalamus of nr0b1 mutant embryos accompanying with the decreased expression of genes in the Notch signaling pathway. Additionally, the inhibition of Notch signaling in wild-type embryos increases the number of pdyn neurons, mimicking the nr0b1 mutant phenotype. In contrast, ectopic activation of Notch signaling in nr0b1 mutant embryos decreases the number of pdyn neurons. Taken together, our results suggest that nr0b1 regulates neural progenitor proliferation and maintenance to ensure normal hypothalamic neuron development.

The combined action of Esrrb and Nr5a2 is essential for murine naïve pluripotency

The maintenance of pluripotency in mouse embryonic stem cells (ESCs) is governed by the action of an interconnected network of transcription factors. Among them, only Oct4 and Sox2 have been shown to be strictly required for the self-renewal of ESCs and pluripotency, particularly in culture conditions in which differentiation cues are chemically inhibited. Here, we report that the conjunct activity of two orphan nuclear receptors, Esrrb and Nr5a2, parallels the importance of that of Oct4 and Sox2 in naïve mouse ESCs. By occupying a large common set of regulatory elements, these two factors control the binding of Oct4, Sox2 and Nanog to DNA. Consequently, in their absence the pluripotency network collapses and the transcriptome is substantially deregulated, leading to the differentiation of ESCs. Altogether, this work identifies orphan nuclear receptors, previously thought to be performing supportive functions, as a set of core regulators of naïve pluripotency.

Summary: Esrrb and Nr5a2, two orphan nuclear receptors, are identified as essential regulators of pluripotency in mouse embryonic stem cells.

Pluripotency-State-Dependent Role of Dax1 in Embryonic Stem Cells Self-Renewal

Dax1(also known as Nr0b1) is regarded as an important component of the transcription factor network in mouse embryonic stem cells (ESCs). However, the role and the molecular mechanism of Dax1 in the maintenance of different pluripotency states are poorly understood. Here, we constructed a stable Dax1 knockout (KO) cell line using the CRISPR/Cas9 system to analyze the precise function of Dax1. We reported that 2i/LIF-ESCs had significantly lower Dax1 expression than LIF/serum-ESCs. Dax1KO ES cell lines could be established in 2i/LIF and their pluripotency was confirmed. In contrast, Dax1-null ESCs could not be continuously passaged in LIF/serum due to severe differentiation and apoptosis. In LIF/serum, the activities of the Core module and Myc module were significantly reduced, while the PRC2 module was activated after Dax1KO. The expression of most proapoptotic genes and lineage-commitment genes were drastically increased, while the downregulated expression of antiapoptotic genes and many pluripotency genes was observed. Our research on the pluripotent state-dependent role of Dax1 provides clues to understand the molecular regulation mechanism at different stages of early embryonic development.

Cooperation between HDAC3 and DAX1 mediates lineage restriction of embryonic stem cells

Mouse embryonic stem cells (mESCs) are biased toward producing embryonic rather than extraembryonic endoderm fates. Here, we identify the mechanism of this barrier and report that the histone deacetylase Hdac3 and the transcriptional corepressor Dax1 cooperatively limit the lineage repertoire of mESCs by silencing an enhancer of the extraembryonic endoderm-specifying transcription factor Gata6. This restriction is opposed by the pluripotency transcription factors Nr5a2 and Esrrb, which promote cell type conversion. Perturbation of the barrier extends mESC potency and allows formation of 3D spheroids that mimic the spatial segregation of embryonic epiblast and extraembryonic endoderm in early embryos. Overall, this study shows that transcriptional repressors stabilize pluripotency by biasing the equilibrium between embryonic and extraembryonic lineages that is hardwired into the mESC transcriptional network.

Molecular signature and colony morphology affect in vitro pluripotency of porcine induced pluripotent stem cells

Overall efficiency of cell reprogramming for porcine fibroblasts into induced pluripotent stem cells (iPSCs) is currently poor, and few cell lines have been established. This study examined gene expression during early phase of cellular reprogramming in the relationship to the iPSC colony morphology and in vitro pluripotent characteristics. Fibroblasts were reprogrammed with OCT4, SOX2, KLF4 and c-MYC. Two different colony morphologies referred to either compact (n = 10) or loose (n = 10) colonies were further examined for proliferative activity, gene expression and in vitro pluripotency. A total of 1,697 iPSC-like colonies (2.34%) were observed after gene transduction. The compact colonies contained with tightly packed cells with a distinct-clear border between the colony and feeder cells, while loose colonies demonstrated irregular colony boundary. For quantitative expression of genes responsible for early phase cell reprogramming, the Dppa2 and EpCAM were significantly upregulated while NR0B1 was downregulated in compact colonies compared with loose phenotype (p < .05). Higher proportion of compact iPSC phenotype (5 of 10, 50%) could be maintained in undifferentiated state for more than 50 passages compared unfavourably with loose morphology (3 of 10, 30%). All iPS cell lines obtained from these two types of colony morphologies expressed pluripotent genes and proteins (OCT4, NANOG and E-cadherin). In addition, they could aggregate and form three-dimensional structure of embryoid bodies. However, only compact iPSC colonies differentiated into three germ layers. Molecular signature of early phase of cell reprogramming coupled with primary colony morphology reflected the in vitro pluripotency of porcine iPSCs. These findings can be simply applied for pre-screening selection of the porcine iPSC cell line.

MEAF6 is essential for cell proliferation and plays a role in the assembly of KAT7 complexes

Myst family genes encode lysine acetyltransferases that mainly mediate histone acetylation to control transcription, DNA replication and DNA damage response. They form tetrameric complexes with PHD-finger proteins (Brpfs or Jades) and small non-catalytic subunits Ing4/5 and Meaf6. Although all the components of the complex are well-conserved from yeast to mammals, the function of Meaf6 and its homologs has not been elucidated in any species. Here we revealed the role of Meaf6 utilizing inducible Meaf6 KO ES cells. By elimination of Meaf6, proliferation ceased although histone acetylations were largely unaffected. In the absence of Meaf6, one of the Myst family members Myst2/Kat7 increased the ability to interact with PHD-finger proteins. This study is the first indication of the function of Meaf6, which shows it is not essential for HAT activity but modulates the assembly of the Kat7 complex.

MLL1 Inhibition and Vitamin D Signaling Cooperate to Facilitate the Expanded Pluripotency State

Dynamic establishment of histone modifications in early development coincides with programed cell fate restriction and loss of totipotency beyond the early blastocyst stage. Causal function of histone-modifying enzymes in this process remains to be defined. Here we show that inhibiting histone methyltransferase MLL1 reprograms naive embryonic stem cells (ESCs) to expanded pluripotent stem cells (EPSCs), with differentiation potential toward both embryonic and extraembryonic lineages in vitro and in vivo. MLL1 inhibition or deletion upregulates gene signatures of early blastomere development. The function of MLL1 in restricting induction of EPSCs is mediated partly by Gc, which regulates cellular response to vitamin D signaling. Combined treatment of MLL1 inhibitor and 1α,25-dihydroxyvitamin D₃ (1,25-(OH)₂D₃) cooperatively enhanced functionality of EPSCs, triggering an extended 2C-like state in vitro and robust totipotent-like property in vivo. Our study sheds light on interplay between epigenetics and vitamin D pathway in cell fate determination.

Toll-like receptor mediated inflammation requires FASN-dependent MYD88 palmitoylation

Toll-like receptor (TLR)/myeloid differentiation primary response protein (MYD88) signaling aggravates sepsis by impairing neutrophil migration to infection sites. However, the role of intracellular fatty acids in TLR/MYD88 signaling is unclear. Here, inhibition of fatty acid synthase by C75 improved neutrophil chemotaxis and increased the survival of mice with sepsis in cecal ligation puncture and lipopolysaccharide-induced septic shock models. C75 specifically blocked TLR/MYD88 signaling in neutrophils. Treatment with GSK2194069 that targets a different domain of fatty acid synthase, did not block TLR signaling or MYD88 palmitoylation. De novo fatty acid synthesis and CD36-mediated exogenous fatty acid incorporation contributed to MYD88 palmitoylation. The binding of IRAK4 to the MYD88 intermediate domain and downstream signal activation required MYD88 palmitoylation at cysteine 113. MYD88 was palmitoylated by ZDHHC6, and ZDHHC6 knockdown decreased MYD88 palmitoylation and TLR/MYD88 activation upon lipopolysaccharide stimulus. Thus, intracellular saturated fatty acid-dependent palmitoylation of MYD88 by ZDHHC6 is a therapeutic target of sepsis.

Dppa2 and Dppa4 directly regulate the Dux-driven zygotic transcriptional program

In this study, Eckersley-Maslin et al. investigated the upstream maternal factors that initiate zygotic genome activation (ZGA) in either a Dux-dependent (a transcription factor expressed in the minor wave of ZGA) or Dux-independent manner. They performed a candidate-based overexpression screen, identifying developmental pluripotency-associated 2 (Dppa2) and Dppa4 as positive regulators of 2C-like cells and transcription of ZGA genes, and their results suggest that Dppa2/4 binding to the Dux promoter leads to Dux up-regulation and activation of the 2C-like transcriptional program, which is subsequently reinforced by Zscan4c.

The molecular regulation of zygotic genome activation (ZGA) in mammals remains an exciting area of research. Primed mouse embryonic stem cells contain a rare subset of “2C-like” cells that are epigenetically and transcriptionally similar to the two-cell embryo and thus represent an in vitro approximation for studying ZGA transcription regulation. Recently, the transcription factor Dux, expressed in the minor wave of ZGA, was described to activate many downstream ZGA transcripts. However, it remains unknown what upstream maternal factors initiate ZGA in either a Dux-dependent or Dux-independent manner. Here we performed a candidate-based overexpression screen, identifying, among others, developmental pluripotency-associated 2 (Dppa2) and Dppa4 as positive regulators of 2C-like cells and transcription of ZGA genes. In the germline, promoter DNA demethylation coincides with expression of Dppa2 and Dppa4, which remain expressed until embryonic day 7.5 (E7.5), when their promoters are remethylated. Furthermore, Dppa2 and Dppa4 are also expressed during induced pluripotent stem cell (iPSC) reprogramming at the time that 2C-like transcription transiently peaks. Through a combination of overexpression, knockdown, knockout, and rescue experiments together with transcriptional analyses, we show that Dppa2 and Dppa4 directly regulate the 2C-like cell population and associated transcripts, including Dux and the Zscan4 cluster. Importantly, we teased apart the molecular hierarchy in which the 2C-like transcriptional program is initiated and stabilized. Dppa2 and Dppa4 require Dux to initiate 2C-like transcription, suggesting that they act upstream by directly regulating Dux. Supporting this, ChIP-seq (chromatin immunoprecipitation [ChIP] combined with high-throughput sequencing) analysis revealed that Dppa2 and Dppa4 bind to the Dux promoter and gene body and drive its expression. Zscan4c is also able to induce 2C-like cells in wild-type cells but, in contrast to Dux, can no longer do so in Dppa2/4 double-knockout cells, suggesting that it may act to stabilize rather than drive the transcriptional network. Our findings suggest a model in which Dppa2/4 binding to the Dux promoter leads to Dux up-regulation and activation of the 2C-like transcriptional program, which is subsequently reinforced by Zscan4c.

A FRET biosensor for necroptosis uncovers two different modes of the release of DAMPs

Necroptosis is a regulated form of necrosis that depends on receptor-interacting protein kinase (RIPK)3 and mixed lineage kinase domain-like (MLKL). While danger-associated molecular pattern (DAMP)s are involved in various pathological conditions and released from dead cells, the underlying mechanisms are not fully understood. Here we develop a fluorescence resonance energy transfer (FRET) biosensor, termed SMART (a sensor for MLKL activation by RIPK3 based on FRET). SMART is composed of a fragment of MLKL and monitors necroptosis, but not apoptosis or necrosis. Mechanistically, SMART monitors plasma membrane translocation of oligomerized MLKL, which is induced by RIPK3 or mutational activation. SMART in combination with imaging of the release of nuclear DAMPs and Live-Cell Imaging for Secretion activity (LCI-S) reveals two different modes of the release of High Mobility Group Box 1 from necroptotic cells. Thus, SMART and LCI-S uncover novel regulation of the release of DAMPs during necroptosis.

Necroptotic cells activate MLKL and release inflammatory DAMPs, although the underlying regulatory mechanisms of this process are poorly understood. Here, Murai et al. develop a necroptosis-specific FRET sensor (SMART) that monitors MLKL membrane translocation to identify two modes of DAMP release.

TGMI: an efficient algorithm for identifying pathway regulators through evaluation of triple-gene mutual interaction

Despite their important roles, the regulators for most metabolic pathways and biological processes remain elusive. Presently, the methods for identifying metabolic pathway and biological process regulators are intensively sought after. We developed a novel algorithm called triple-gene mutual interaction (TGMI) for identifying these regulators using high-throughput gene expression data. It first calculated the regulatory interactions among triple gene blocks (two pathway genes and one transcription factor (TF)), using conditional mutual information, and then identifies significantly interacted triple genes using a newly identified novel mutual interaction measure (MIM), which was substantiated to reflect strengths of regulatory interactions within each triple gene block. The TGMI calculated the MIM for each triple gene block and then examined its statistical significance using bootstrap. Finally, the frequencies of all TFs present in all significantly interacted triple gene blocks were calculated and ranked. We showed that the TFs with higher frequencies were usually genuine pathway regulators upon evaluating multiple pathways in plants, animals and yeast. Comparison of TGMI with several other algorithms demonstrated its higher accuracy. Therefore, TGMI will be a valuable tool that can help biologists to identify regulators of metabolic pathways and biological processes from the exploded high-throughput gene expression data in public repositories.

Overlapping functions of Krüppel-like factor family members: targeting multiple transcription factors to maintain the naïve pluripotency of mouse embryonic stem cells

Krüppel-like factors (Klfs) have a pivotal role in maintaining self-renewal of mouse embryonic stem cells (mESCs). The functions of three Klf family members (Klf2, Klf4 and Klf5) have been identified, and are suggested to largely overlap. For further dissection of their functions, we applied an inducible knockout system for these Klf family members and assessed the effects of combinatorial loss of function. As a result, we confirmed that any one of Klf2, Klf4 and Klf5 was sufficient to support self-renewal, whereas the removal of all three compromised it. The activity of any single transcription factor, except for a Klf family member, was not sufficient to restore self-renewal of triple-knockout mESCs. However, some particular combinations of transcription factors were capable of the restoration. The triple-knockout mESCs were successfully captured at primed state. These data indicate that the pivotal function of a Klf family member is transduced into the activation of multiple transcription factors in a naïve-state-specific manner.

The principles that govern transcription factor network functions in stem cells

Tissue-specific transcription factors primarily act to define the phenotype of the cell. The power of a single transcription factor to alter cell fate is often minimal, as seen in gain-of-function analyses, but when multiple transcription factors cooperate synergistically it potentiates their ability to induce changes in cell fate. By contrast, transcription factor function is often dispensable in the maintenance of cell phenotype, as is evident in loss-of-function assays. Why does this phenomenon, commonly known as redundancy, occur? Here, I discuss the role that transcription factor networks play in collaboratively regulating stem cell fate and differentiation by providing multiple explanations for their functional redundancy.

Epigenetic modifications promote the expression of the orphan nuclear receptor NR0B1 in human lung adenocarcinoma cells

The ectopic activation of NR0B1 is involved in the development of some cancers. However, the regulatory mechanisms controlling NR0B1 expression are not well understood. Therefore, the epigenetic modifications promoting NR0B1 activation were examined in this study. NR0B1 protein was detected in cancerous tissues of more than 50% of human lung adenocarcinoma (ADCA) cases and tended to be expressed in low-differentiated cancerous tissues obtained from males. Nevertheless, NR0B1 activation in ADCA has not previously been correlated with DNA demethylation. NR0B1 expression was not detected in 293T cells, although it contains a hypomethylated NR0B1 promoter. Treating 293T cells with a histone deacetylase inhibitor increased acetylated histone H4 binding to the NR0B1 promoter and activated NR0B1 expression. In contrast, treatment with histone methylase inhibitors decreased the methylation of histones H3K9 and H3K27 and slightly induced NR0B1 transcription. Furthermore, the level of acetyl-histone H4 binding to the NR0B1 promoter increased, whereas the occupancy of H3K27me3 was lower in cancerous tissues than in non-cancerous tissues. Similar histone occupancies were confirmed in a comparison of cancerous tissues with strong, moderate and negative NR0B1 expression. In conclusion, this study shows that CpG methylation within the NR0B1 promoter is not involved in the in vivo regulation of NR0B1 expression, whereas the hyperacetylation of histone H4 and the unmethylation of histones H3K9 and H3K27, and their binding to the NR0B1 promoter results in decondensed euchromatin for NR0B1 activation.

MERVL/Zscan4 Network Activation Results in Transient Genome-wide DNA Demethylation of mESCs

Mouse embryonic stem cells are dynamic and heterogeneous. For example, rare cells cycle through a state characterized by decondensed chromatin and expression of transcripts, including the Zscan4 cluster and MERVL endogenous retrovirus, which are usually restricted to preimplantation embryos. Here, we further characterize the dynamics and consequences of this transient cell state. Single-cell transcriptomics identified the earliest upregulated transcripts as cells enter the MERVL/Zscan4 state. The MERVL/Zscan4 transcriptional network was also upregulated during induced pluripotent stem cell reprogramming. Genome-wide DNA methylation and chromatin analyses revealed global DNA hypomethylation accompanying increased chromatin accessibility. This transient DNA demethylation was driven by a loss of DNA methyltransferase proteins in the cells and occurred genome-wide. While methylation levels were restored once cells exit this state, genomic imprints remained hypomethylated, demonstrating a potential global and enduring influence of endogenous retroviral activation on the epigenome.

•

Single-cell transcriptomics reveals dynamics of MERVL/Zscan4 network activation

•

MERVL-LTR transcriptional network is expressed in iPSC reprogramming events

•

Translation block depletes Dnmt proteins, inducing transient global demethylation

•

Passage through the MERVL/Zscan4 state may cause irreversible imprint erasure

Totipotency in the mouse

In mammals, the unicellular zygote starts the process of embryogenesis and differentiates into all types of somatic cells, including both fetal and extraembryonic lineages-in a highly organized manner to eventually give rise to an entire multicellular organism comprising more than 200 different tissue types. This feature is referred to as totipotency. Upon fertilization, oocyte maternal factors epigenetically reprogram the genomes of the terminally differentiated oocyte and spermatozoon and turn the zygote into a totipotent cell. Today, we still do not fully understand the molecular properties of totipotency. In this review, we discuss recent findings on the molecular signature and mechanism of transcriptional regulation networks in the totipotent mouse embryo.

Ground rules of the pluripotency gene regulatory network

Pluripotency is a state that exists transiently in the early embryo and, remarkably, can be recapitulated in vitro by deriving embryonic stem cells or by reprogramming somatic cells to become induced pluripotent stem cells. The state of pluripotency, which is stabilized by an interconnected network of pluripotency-associated genes, integrates external signals and exerts control over the decision between self-renewal and differentiation at the transcriptional, post-transcriptional and epigenetic levels. Recent evidence of alternative pluripotency states indicates the regulatory flexibility of this network. Insights into the underlying principles of the pluripotency network may provide unprecedented opportunities for studying development and for regenerative medicine.

Linking Telomere Regulation to Stem Cell Pluripotency

Embryonic stem cells (ESCs), somatic cell nuclear transfer ESCs, and induced pluripotent stem cells (iPSCs) represent the most studied group of PSCs. Unlimited self-renewal without incurring chromosomal instability and pluripotency are essential for the potential use of PSCs in regenerative therapy. Telomere length maintenance is critical for the unlimited self-renewal, pluripotency, and chromosomal stability of PSCs. While telomerase has a primary role in telomere maintenance, alternative lengthening of telomere pathways involving recombination and epigenetic modifications are also required for telomere length regulation, notably in mouse PSCs. Telomere rejuvenation is part of epigenetic reprogramming to pluripotency. Insights into telomere reprogramming and maintenance in PSCs may have implications for understanding of aging and tumorigenesis. Here, I discuss the link between telomere elongation and homeostasis to the acquisition and maintenance of stem cell pluripotency, and their regulatory mechanisms by epigenetic modifications.

Expression analysis of the endogenous Zscan4 locus and its coding proteins in mouse ES cells and preimplantation embryos

Mouse Zinc finger and SCAN domain containing 4 (Zscan4) is encoded in multiple copies of Zscan4 genes, which are expressed in late two-cell stage preimplantation embryos and in 1–5% of the embryonic stem (ES) cell population at a given time. Due to the highly identical nucleotide sequences of multiple copies of Zscan4 paralogs and pseudogenes in the mouse Zscan4 genomic cluster, previous analyses have been done using exogenous transgenes under the regulation of Zscan4c promoter. In this manuscript, we generated knock-in mouse ES cell lines and mouse lines, in which the expression of endogenous Zscan4c, one of the Zscan4 genes, can be specifically monitored with a green fluorescent protein variant, Emerald. Interestingly, we found that only ∼30% of Zscan4-immunopositive ES cells were Emerald positive, suggesting that even when the Zscan4 locus is active, not all Zscan4 genes are expressed synchronously. We also carried out mass spectrometry of protein complexes associated with endogenous Zscan4 proteins. Taken together, our genetic engineering at an endogenous Zscan4c gene provides the first clue for the expression and function of each gene copy of Zscan4 locus in a physiological context.

MERVL/Zscan4 Network Activation Results in Transient Genome-wide DNA Demethylation of mESCs

Mouse embryonic stem cells are dynamic and heterogeneous. For example, rare cells cycle through a state characterized by decondensed chromatin and expression of transcripts, including the Zscan4 cluster and MERVL endogenous retrovirus, which are usually restricted to preimplantation embryos. Here, we further characterize the dynamics and consequences of this transient cell state. Single-cell transcriptomics identified the earliest upregulated transcripts as cells enter the MERVL/Zscan4 state. The MERVL/Zscan4 transcriptional network was also upregulated during induced pluripotent stem cell reprogramming. Genome-wide DNA methylation and chromatin analyses revealed global DNA hypomethylation accompanying increased chromatin accessibility. This transient DNA demethylation was driven by a loss of DNA methyltransferase proteins in the cells and occurred genome-wide. While methylation levels were restored once cells exit this state, genomic imprints remained hypomethylated, demonstrating a potential global and enduring influence of endogenous retroviral activation on the epigenome.

Esrrb directly binds to Gata6 promoter and regulates its expression with Dax1 and Ncoa3

Estrogen-related receptor beta (Esrrb) is expressed in embryonic stem (ES) cells and is involved in self-renewal ability and pluripotency. Previously, we found that Dax1 is associated with Esrrb and represses its transcriptional activity. Further, the disruption of the Dax1-Esrrb interaction increases the expression of the extra-embryonic endoderm marker Gata6 in ES cells. Here, we investigated the influences of Esrrb and Dax1 on Gata6 expression. Esrrb overexpression in ES cells induced endogenous Gata6 mRNA and Gata6 promoter activity. In addition, the Gata6 promoter was found to contain the Esrrb recognition motifs ERRE1 and ERRE2, and the latter was the responsive element of Esrrb. Associations between ERRE2 and Esrrb were then confirmed by biotin DNA pulldown and chromatin immunoprecipitation assays. Subsequently, we showed that Esrrb activity at the Gata6 promoter was repressed by Dax1, and although Dax1 did not bind to ERRE2, it was associated with Esrrb, which directly binds to ERRE2. In addition, the transcriptional activity of Esrrb was enhanced by nuclear receptor co-activator 3 (Ncoa3), which has recently been shown to be a binding partner of Esrrb. Finally, we showed that Dax1 was associated with Ncoa3 and repressed its transcriptional activity. Taken together, the present study indicates that the Gata6 promoter is activated by Esrrb in association with Ncoa3, and Dax1 inhibited activities of Esrrb and Ncoa3, resulting maintenance of the undifferentiated status of ES cells.

Emergence of undifferentiated colonies from mouse embryonic stem cells undergoing differentiation by retinoic acid treatment

Retinoic acid (RA) is one of the most potent inducers of differentiation of mouse embryonic stem cells (ESCs). However, previous studies show that RA treatment of cells cultured in the presence of a leukemia inhibitory factor (LIF) also result in the upregulation of a gene called Zscan4, whose transient expression is a marker for undifferentiated ESCs. We explored the balance between these two seemingly antagonistic effects of RA. ESCs indeed differentiated in the presence of LIF after RA treatment, but colonies of undifferentiated ESCs eventually emerged from these differentiated cells - even in the presence of RA. These colonies, named secondary colonies, consist of three cell types: typical undifferentiated ESCs expressing pluripotency genes such as Pou5f1, Sox2, and Nanog; cells expressing Zscan4; and endodermal-like cells located at the periphery of the colony. The capacity to form secondary colonies was confirmed for all eight tested ESC lines. Cells from the secondary colonies - after transfer to the standard ESC medium - retained pluripotency, judged by their strong alkaline phosphatase (ALP) staining, typical colony morphology, gene expression profile, stable karyotype, capacity to differentiate into all three germ layers in embryoid body formation assays, and successful contribution to chimeras after injection into blastocysts. Based on flow cytometry analysis (FACS), the proportion of Zscan4-positive cells in secondary colonies was higher than in standard ESC colonies, which may explain the capacity of ESCs to resist the differentiating effects of RA and instead form secondary colonies of undifferentiated ESCs. This hypothesis is supported by cell-lineage tracing analysis, which showed that most cells in the secondary colonies were descendents of cells transiently expressing Zscan4.

Zscan4 Is Activated after Telomere Shortening in Mouse Embryonic Stem Cells

ZSCAN4 is a DNA-binding protein that functions for telomere elongation and genomic stability. In vivo, it is specifically expressed at the two-cell stage during mouse development. In vitro, it is transiently expressed in mouse embryonic stem cells (ESCs), only in 5% of the population at one time. Here we attempted to elucidate when, under what circumstances, Zscan4 is activated in ESCs. Using live cell imaging, we monitored the activity of Zscan4 together with the pluripotency marker Rex1. The lengths of the cell cycles in ESCs were diverse. Longer cell cycles were accompanied by shorter telomeres and higher activation of Zscan4. Since activation of Zscan4 is involved in telomere elongation, we speculate that the extended cell cycles accompanied by Zscan4 activation reflect the time for telomere recovery. Rex1 and Zscan4 did not show any correlation. Taken together, we propose that Zscan4 is activated to recover shortened telomeres during extended cell cycles, irrespective of the pluripotent status.

•

At longer cell cycles, telomeres are shorter

•

Zscan4 is activated when the cell cycles become long

•

After the activation of Zscan4, the next cell cycle becomes short

•

We propose Zscan4 is activated for telomere maintenance irrespective of pluripotency

Tcstv1 and Tcstv3 elongate telomeres of mouse ES cells

Mouse embryonic stem cell (ESC) cultures exhibit a heterogeneous mixture of metastable cells sporadically entering the 2-cell (2C)-embryo-like state, critical for ESC potency. One of 2-cell genes, Zscan4, has been shown to be responsible for telomere maintenance, genomic stability and pluripotency of mouse ESCs. Functions of other 2C-genes in ESCs remain elusive. Here we show that 2C-genes Tcstv1 and Tcstv3 play a role in regulation of telomere lengths. Overexpression or knockdown Tcstv1 and Tcstv3 does not immediately affect proliferation, pluripotency and differentiation in vitro of ESCs. However, ectopic expression of Tcstv1 or Tcstv3 results in telomere elongation, whereas Tcstv1/3 knockdown shortens telomeres of ESCs. Overexpression of Tcstv1 or Tcstv3 does not alter telomere stability. Furthermore, Tcstv1 can increase Zscan4 protein levels and telomere recombination by telomere sister chromatid exchange (T-SCE). Depletion of Tcstv1/3 reduces Zscan4 protein levels. Together, Tcstv1 and Tcstv3 are involved in telomere maintenance that is required for long-term self-renewal of mouse ESCs. Our data also suggests that Tcstv1/3 may co-operate and stabilize Zscan4 protein but the molecular bases remain to be determined.

Lineage Segregation in the Totipotent Embryo

After a spermatozoon enters an oocyte, maternal factors accumulated in the oocyte reprogram the genomes of the terminally differentiated oocyte and spermatozoon epigenetically and turn the zygote into a totipotent cell, with the capacity to differentiate into all types of somatic cells in a highly organized manner and generate the entire organism, a feature referred to as totipotency. Differentiation of the first lineage begins after three cleavages, when the early embryo compacts and becomes polarized, followed by segregation of the first lineages--the inner cell mass (ICM) and the trophectoderm (TE). To date, a full understanding of the molecular mechanisms that underlie the establishment of totipotency and the ICM/TE lineage segregation remains unclear. In this review, we discuss recent findings in the mechanism of transcriptional regulation networks and signaling pathways in the first lineage separation in the totipotent mouse embryo.

LIF signal in mouse embryonic stem cells

Since the establishment of mouse embryonic stem cells (mESCs) in the 1980s, a number of important notions on the self-renewal of pluripotent stem cells in vitro have been found. In serum containing conventional culture, an exogenous cytokine, leukemia inhibitory factor (LIF), is absolutely essential for the maintenance of pluripotency. In contrast, in serum-free culture with simultaneous inhibition of Map-kinase and Gsk3 (so called 2i-culture), LIF is no longer required. However, recent findings also suggest that LIF may have a role not covered by the 2i for the maintenance of naïve pluripotency. These suggest that LIF functions for the maintenance of naïve pluripotency in a context dependent manner. We summarize how LIF-signal pathway is converged to maintain the naïve state of pluripotency.

ETS-related transcription factors ETV4 and ETV5 are involved in proliferation and induction of differentiation-associated genes in embryonic stem (ES) cells

The pluripotency and self-renewal capacity of embryonic stem (ES) cells is regulated by several transcription factors. Here, we show that the ETS-related transcription factors Etv4 and Etv5 (Etv4/5) are specifically expressed in undifferentiated ES cells, and suppression of Oct3/4 results in down-regulation of Etv4/5. Simultaneous deletion of Etv4 and Etv5 (Etv4/5 double knock-out (dKO)) in ES cells resulted in a flat, epithelial cell-like appearance, whereas the morphology changed into compact colonies in a 2i medium (containing two inhibitors for GSK3 and MEK/ERK). Expression levels of self-renewal marker genes, including Oct3/4 and Nanog, were similar between wild-type and dKO ES cells, whereas proliferation of Etv4/5 dKO ES cells was decreased with overexpression of cyclin-dependent kinase inhibitors (p16/p19, p15, and p57). A differentiation assay revealed that the embryoid bodies derived from Etv4/5 dKO ES cells were smaller than the control, and expression of ectoderm marker genes, including Fgf5, Sox1, and Pax3, was not induced in dKO-derived embryoid bodies. Microarray analysis demonstrated that stem cell-related genes, including Tcf15, Gbx2, Lrh1, Zic3, and Baf60c, were significantly repressed in Etv4/5 dKO ES cells. The artificial expression of Etv4 and/or Etv5 in Etv4/5 dKO ES cells induced re-expression of Tcf15 and Gbx2. These results indicate that Etv4 and Etv5, potentially through regulation of Gbx2 and Tcf15, are involved in the ES cell proliferation and induction of differentiation-associated genes in ES cells.