Systematic repression of transcription factors reveals limited patterns of gene expression changes in ES cells

Atlas of regulated target genes of transcription factors (ART-TF) in human ES cells

Background

Transcription factors (TFs) play central roles in maintaining “stemness” of embryonic stem (ES) cells and their differentiation into several hundreds of adult cell types. The regulatory competence of TFs is routinely assessed by detecting target genes to which they bind. However, these data do not indicate which target genes are activated, repressed, or not affected by the change of TF abundance. There is a lack of large-scale studies that compare the genome binding of TFs with the expression change of target genes after manipulation of each TF.

Results

In this paper we associated human TFs with their target genes by two criteria: binding to genes, evaluated from published ChIP-seq data (n = 1868); and change of target gene expression shortly after induction of each TF in human ES cells. Lists of direction- and strength-specific regulated target genes are generated for 311 TFs (out of 351 TFs tested) with expected proportion of false positives less than or equal to 0.30, including 63 new TFs not present in four existing databases of target genes. Our lists of direction-specific targets for 152 TFs (80.0%) are larger that in the TRRUST database. In average, 30.9% of genes that respond greater than or equal to twofold to the induction of TFs are regulated targets. Regulated target genes indicate that the majority of TFs are either strong activators or strong repressors, whereas sets of genes that responded greater than or equal to twofold to the induction of TFs did not show strong asymmetry in the direction of expression change. The majority of human TFs (82.1%) regulated their target genes primarily via binding to enhancers. Repression of target genes is more often mediated by promoter-binding than activation of target genes. Enhancer-promoter loops are more abundant among strong activator and repressor TFs.

Conclusions

We developed an atlas of regulated targets of TFs (ART-TF) in human ES cells by combining data on TF binding with data on gene expression change after manipulation of individual TFs. Sets of regulated gene targets were identified with a controlled rate of false positives. This approach contributes to the understanding of biological functions of TFs and organization of gene regulatory networks. This atlas should be a valuable resource for ES cell-based regenerative medicine studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04924-3.

Dynamic regulatory module networks for inference of cell type-specific transcriptional networks

Changes in transcriptional regulatory networks can significantly alter cell fate. To gain insight into transcriptional dynamics, several studies have profiled bulk multi-omic data sets with parallel transcriptomic and epigenomic measurements at different stages of a developmental process. However, integrating these data to infer cell type-specific regulatory networks is a major challenge. We present dynamic regulatory module networks (DRMNs), a novel approach to infer cell type-specific cis-regulatory networks and their dynamics. DRMN integrates expression, chromatin state, and accessibility to predict cis-regulators of context-specific expression, where context can be cell type, developmental stage, or time point, and uses multitask learning to capture network dynamics across linearly and hierarchically related contexts. We applied DRMNs to study regulatory network dynamics in three developmental processes, each showing different temporal relationships and measuring a different combination of regulatory genomic data sets: cellular reprogramming, liver dedifferentiation, and forward differentiation. DRMN identified known and novel regulators driving cell type-specific expression patterns, showing its broad applicability to examine dynamics of gene regulatory networks from linearly and hierarchically related multi-omic data sets.

Arp2/3 complex activity is necessary for mouse ESC differentiation, times formative pluripotency, and enables lineage specification

Mouse embryonic stem cells (mESCs), a model for differentiation into primed epiblast-like cells (EpiLCs), have revealed transcriptional and epigenetic control of early embryonic development. The control and significance of morphological changes, however, remain less defined. We show marked changes in morphology and actin architectures during differentiation that depend on Arp2/3 complex but not formin activity. Inhibiting Arp2/3 complex activity pharmacologically or genetically does not block exit from naive pluripotency, but attenuates increases in EpiLC markers. We find that inhibiting Arp2/3 complex activity delays formative pluripotency and causes globally defective lineage specification as indicated by RNA sequencing, with significant effects on TBX3-depedendent transcriptional programs. We also identify two previously unreported indicators of mESC differentiation, namely, MRTF and FHL2, which have inverse Arp2/3 complex-dependent nuclear translocation. Our findings on Arp2/3 complex activity in differentiation and the established role of formins in EMT indicate that these two actin nucleators regulate distinct modes of epithelial plasticity.

Transcriptome analysis reveals high tumor heterogeneity with respect to re-activation of stemness and proliferation programs

Significant alterations in signaling pathways and transcriptional regulatory programs together represent major hallmarks of many cancers. These, among all, include the reactivation of stemness, which is registered by the expression of pathways that are active in the embryonic stem cells (ESCs). Here, we assembled gene sets that reflect the stemness and proliferation signatures and used them to analyze a large panel of RNA-seq data from The Cancer Genome Atlas (TCGA) Consortium in order to specifically assess the expression of stemness-related and proliferation-related genes across a collection of different tumor types. We introduced a metric that captures the collective similarity of the expression profile of a tumor to that of ESCs, which showed that stemness and proliferation signatures vary greatly between different tumor types. We also observed a high degree of intertumoral heterogeneity in the expression of stemness- and proliferation-related genes, which was associated with increased hazard ratios in a fraction of tumors and mirrored by high intratumoral heterogeneity and a remarkable stemness capacity in metastatic lesions across cancer cells in single cell RNA-seq datasets. Taken together, these results indicate that the expression of stemness signatures is highly heterogeneous and cannot be used as a universal determinant of cancer. This calls into question the universal validity of diagnostic tests that are based on stem cell markers.

Cell fate conversion prediction by group sparse optimization method utilizing single-cell and bulk OMICs data

Cell fate conversion by overexpressing defined factors is a powerful tool in regenerative medicine. However, identifying key factors for cell fate conversion requires laborious experimental efforts; thus, many of such conversions have not been achieved yet. Nevertheless, cell fate conversions found in many published studies were incomplete as the expression of important gene sets could not be manipulated thoroughly. Therefore, the identification of master transcription factors for complete and efficient conversion is crucial to render this technology more applicable clinically. In the past decade, systematic analyses on various single-cell and bulk OMICs data have uncovered numerous gene regulatory mechanisms, and made it possible to predict master gene regulators during cell fate conversion. By virtue of the sparse structure of master transcription factors and the group structure of their simultaneous regulatory effects on the cell fate conversion process, this study introduces a novel computational method predicting master transcription factors based on group sparse optimization technique integrating data from multi-OMICs levels, which can be applicable to both single-cell and bulk OMICs data with a high tolerance of data sparsity. When it is compared with current prediction methods by cross-referencing published and validated master transcription factors, it possesses superior performance. In short, this method facilitates fast identification of key regulators, give raise to the possibility of higher successful conversion rate and in the hope of reducing experimental cost.

The SAM domain-containing protein 1 (SAMD1) acts as a repressive chromatin regulator at unmethylated CpG islands

CpG islands (CGIs) are key regulatory DNA elements at most promoters, but how they influence the chromatin status and transcription remains elusive. Here, we identify and characterize SAMD1 (SAM domain-containing protein 1) as an unmethylated CGI-binding protein. SAMD1 has an atypical winged-helix domain that directly recognizes unmethylated CpG-containing DNA via simultaneous interactions with both the major and the minor groove. The SAM domain interacts with L3MBTL3, but it can also homopolymerize into a closed pentameric ring. At a genome-wide level, SAMD1 localizes to H3K4me3-decorated CGIs, where it acts as a repressor. SAMD1 tethers L3MBTL3 to chromatin and interacts with the KDM1A histone demethylase complex to modulate H3K4me2 and H3K4me3 levels at CGIs, thereby providing a mechanism for SAMD1-mediated transcriptional repression. The absence of SAMD1 impairs ES cell differentiation processes, leading to misregulation of key biological pathways. Together, our work establishes SAMD1 as a newly identified chromatin regulator acting at unmethylated CGIs.

An Esrrb and Nanog Cell Fate Regulatory Module Controlled by Feed Forward Loop Interactions

Cell fate decisions during development are governed by multi-factorial regulatory mechanisms including chromatin remodeling, DNA methylation, binding of transcription factors to specific loci, RNA transcription and protein synthesis. However, the mechanisms by which such regulatory “dimensions” coordinate cell fate decisions are currently poorly understood. Here we quantified the multi-dimensional molecular changes that occur in mouse embryonic stem cells (mESCs) upon depletion of Estrogen related receptor beta (Esrrb), a key pluripotency regulator. Comparative analyses of expression changes subsequent to depletion of Esrrb or Nanog, indicated that a system of interlocked feed-forward loops involving both factors, plays a central part in regulating the timing of mESC fate decisions. Taken together, our meta-analyses support a hierarchical model in which pluripotency is maintained by an Oct4-Sox2 regulatory module, while the timing of differentiation is regulated by a Nanog-Esrrb module.

Metabolic control of DNA methylation in naive pluripotent cells

Naive epiblast and embryonic stem cells (ESCs) give rise to all cells of adults. Such developmental plasticity is associated with genome hypomethylation. Here, we show that LIF-Stat3 signaling induces genomic hypomethylation via metabolic reconfiguration. Stat3-/- ESCs show decreased α-ketoglutarate production from glutamine, leading to increased Dnmt3a and Dnmt3b expression and DNA methylation. Notably, genome methylation is dynamically controlled through modulation of α-ketoglutarate availability or Stat3 activation in mitochondria. Alpha-ketoglutarate links metabolism to the epigenome by reducing the expression of Otx2 and its targets Dnmt3a and Dnmt3b. Genetic inactivation of Otx2 or Dnmt3a and Dnmt3b results in genomic hypomethylation even in the absence of active LIF-Stat3. Stat3-/- ESCs show increased methylation at imprinting control regions and altered expression of cognate transcripts. Single-cell analyses of Stat3-/- embryos confirmed the dysregulated expression of Otx2, Dnmt3a and Dnmt3b as well as imprinted genes. Several cancers display Stat3 overactivation and abnormal DNA methylation; therefore, the molecular module that we describe might be exploited under pathological conditions.

Identification of a dynamic gene regulatory network required for pluripotency factor-induced reprogramming of mouse fibroblasts and hepatocytes

The generation of induced pluripotent stem cells (iPSCs) from somatic cells provides an excellent model to study mechanisms of transcription factor-induced global alterations of the epigenome and genome function. Here, we have investigated the early transcriptional events of cellular reprogramming triggered by the co-expression of Oct4, Sox2, Klf4, and c-Myc (OSKM) in mouse embryonic fibroblasts (MEFs) and mouse hepatocytes (mHeps). In this analysis, we identified a gene regulatory network composed of nine transcriptional regulators (9TR; Cbfa2t3, Gli2, Irf6, Nanog, Ovol1, Rcan1, Taf1c, Tead4, and Tfap4), which are directly targeted by OSKM, in vivo. Functional studies using single and double shRNA knockdowns of any of these factors caused disruption of the network and dramatic reductions in reprogramming efficiency, indicating that this network is essential for the induction and establishment of pluripotency. We demonstrate that the stochastic co-expression of 9TR network components occurs in a remarkably small number of cells, approximating the percentage of terminally reprogrammed cells as a result of dynamic molecular events. Thus, the early DNA-binding patterns of OSKM and the subsequent probabilistic co-expression of essential 9TR components in subpopulations of cells undergoing reprogramming steer the reconstruction of a gene regulatory network marking the transition to pluripotency.

Perturbation-based gene regulatory network inference to unravel oncogenic mechanisms

The gene regulatory network (GRN) of human cells encodes mechanisms to ensure proper functioning. However, if this GRN is dysregulated, the cell may enter into a disease state such as cancer. Understanding the GRN as a system can therefore help identify novel mechanisms underlying disease, which can lead to new therapies. To deduce regulatory interactions relevant to cancer, we applied a recent computational inference framework to data from perturbation experiments in squamous carcinoma cell line A431. GRNs were inferred using several methods, and the false discovery rate was controlled by the NestBoot framework. We developed a novel approach to assess the predictiveness of inferred GRNs against validation data, despite the lack of a gold standard. The best GRN was significantly more predictive than the null model, both in cross-validated benchmarks and for an independent dataset of the same genes under a different perturbation design. The inferred GRN captures many known regulatory interactions central to cancer-relevant processes in addition to predicting many novel interactions, some of which were experimentally validated, thus providing mechanistic insights that are useful for future cancer research.

Transcriptional network dynamics during the progression of pluripotency revealed by integrative statistical learning

The developmental potential of cells, termed pluripotency, is highly dynamic and progresses through a continuum of naive, formative and primed states. Pluripotency progression of mouse embryonic stem cells (ESCs) from naive to formative and primed state is governed by transcription factors (TFs) and their target genes. Genomic techniques have uncovered a multitude of TF binding sites in ESCs, yet a major challenge lies in identifying target genes from functional binding sites and reconstructing dynamic transcriptional networks underlying pluripotency progression. Here, we integrated time-resolved ‘trans-omic’ datasets together with TF binding profiles and chromatin conformation data to identify target genes of a panel of TFs. Our analyses revealed that naive TF target genes are more likely to be TFs themselves than those of formative TFs, suggesting denser hierarchies among naive TFs. We also discovered that formative TF target genes are marked by permissive epigenomic signatures in the naive state, indicating that they are poised for expression prior to the initiation of pluripotency transition to the formative state. Finally, our reconstructed transcriptional networks pinpointed the precise timing from naive to formative pluripotency progression and enabled the spatiotemporal mapping of differentiating ESCs to their in vivo counterparts in developing embryos.

Linc1557 is critical for the initiation of embryonic stem cell differentiation by directly targeting the LIF/STAT3 signaling pathway

Embryonic stem cells (ESCs) have self-renewal and multi-lineage differentiation potential and perform critical functions in development and biomedicine. Several long noncoding RNAs (lncRNAs) have been reported as key regulators of stem cell pluripotency and differentiation. However, the function and regulatory mechanism of lncRNAs during the initiation of ESC differentiation remains unclear. Here, we found that linc1557 was highly expressed in mouse ESCs and required for the initiation of ESC differentiation. Knockdown of linc1557 increased the expression and phosphorylation levels of signal transducer and activator of transcription 3 (STAT3), a key factor in the leukemia inhibitory factor (LIF)/STAT3 signaling pathway. Furthermore, we found that linc1557 directly bound to Stat3 mRNA and affected its stability. The differentially expressed transcriptome after linc1557 knockdown in ESCs was involved primarily in multicellular organism development and cell differentiation as similar to that after Stat3 knockdown. Moreover, either knockdown of Stat3 or addition of a LIF/STAT3 signaling inhibitor rescued the suppressive effects of linc1557 knockdown on the initiation of mouse ESC differentiation. These findings not only elucidated the critical function of linc1557 in the initiation of mouse ESC differentiation but also clarified that its specific mechanism as directly affecting Stat3 mRNA stability, which enhanced the understanding of the lncRNA-mediated regulatory mechanism for mRNA stability and key signaling pathways in ESC pluripotency and differentiation.

Integrative and perturbation-based analysis of the transcriptional dynamics of TGFβ/BMP system components in transition from embryonic stem cells to neural progenitors

Cooperative actions of extrinsic signals and cell‐intrinsic transcription factors alter gene regulatory networks enabling cells to respond appropriately to environmental cues. Signaling by transforming growth factor type β (TGFβ) family ligands (eg, bone morphogenetic proteins [BMPs] and Activin/Nodal) exerts cell‐type specific and context‐dependent transcriptional changes, thereby steering cellular transitions throughout embryogenesis. Little is known about coordinated regulation and transcriptional interplay of the TGFβ system. To understand intrafamily transcriptional regulation as part of this system's actions during development, we selected 95 of its components and investigated their mRNA‐expression dynamics, gene‐gene interactions, and single‐cell expression heterogeneity in mouse embryonic stem cells transiting to neural progenitors. Interrogation at 24 hour intervals identified four types of temporal gene transcription profiles that capture all stages, that is, pluripotency, epiblast formation, and neural commitment. Then, between each stage we performed esiRNA‐based perturbation of each individual component and documented the effect on steady‐state mRNA levels of the remaining 94 components. This exposed an intricate system of multilevel regulation whereby the majority of gene‐gene interactions display a marked cell‐stage specific behavior. Furthermore, single‐cell RNA‐profiling at individual stages demonstrated the presence of detailed co‐expression modules and subpopulations showing stable co‐expression modules such as that of the core pluripotency genes at all stages. Our combinatorial experimental approach demonstrates how intrinsically complex transcriptional regulation within a given pathway is during cell fate/state transitions.

The gene regulatory network of mESC differentiation: a benchmark for reverse engineering methods

A large body of data have accumulated that characterize the gene regulatory network of stem cells. Yet, a comprehensive and integrative understanding of this complex network is lacking. Network reverse engineering methods that use transcriptome data to derive these networks may help to uncover the topology in an unbiased way. Many methods exist that use co-expression to reconstruct networks. However, it remains unclear how these methods perform in the context of stem cell differentiation, as most systematic assessments have been made for regulatory networks of unicellular organisms. Here, we report a systematic benchmark of different reverse engineering methods against functional data. We show that network pruning is critical for reconstruction performance. We also find that performance is similar for algorithms that use different co-expression measures, i.e. mutual information or correlation. In addition, different methods yield very different network topologies, highlighting the challenge of interpreting these resulting networks as a whole.This article is part of the theme issue 'Designer human tissue: coming to a lab near you'.

Somatic Pluripotent Genes in Tissue Repair, Developmental Disease, and Cancer

Embryonic stem cells possess the ability to differentiate into all cell types of the body. This pliable developmental state is achieved by the function of a series of pluripotency factors, classically identified as OCT4, SOX2, and NANOG. These pluripotency factors are responsible for activating the larger pluripotency networks and the self-renewal programs which give ES cells their unique characteristics. However, during differentiation pluripotency networks become downregulated as cells achieve greater lineage specification and exit the cell cycle. Typically the repression of pluripotency is viewed as a positive factor to ensure the fidelity of cellular identity by restricting cellular pliancy. Consistent with this view, the expression of pluripotency factors is greatly restricted in somatic cells. However, there are examples whereby cells either maintain or reactivate pluripotency factors to preserve the increased potential for the healing of wounds or tissue homeostasis. Additionally there are many examples where these pluripotency factors become reactivated in a variety of human pathologies, particularly cancer. In this review, we will summarize the somatic repression of pluripotency factors, their role in tissue homeostasis and wound repair, and the human diseases that are associated with pluripotency factor misregulation with an emphasis on their role in the etiology of multiple cancers.

EPOP Interacts with Elongin BC and USP7 to Modulate the Chromatin Landscape

Gene regulatory networks are pivotal for many biological processes. In mouse embryonic stem cells (mESCs), the transcriptional network can be divided into three functionally distinct modules: Polycomb, Core, and Myc. The Polycomb module represses developmental genes, while the Myc module is associated with proliferative functions, and its mis-regulation is linked to cancer development. Here, we show that, in mESCs, the Polycomb repressive complex 2 (PRC2)-associated protein EPOP (Elongin BC and Polycomb Repressive Complex 2-associated protein; a.k.a. C17orf96, esPRC2p48, and E130012A19Rik) co-localizes at chromatin with members of the Myc and Polycomb module. EPOP interacts with the transcription elongation factor Elongin BC and the H2B deubiquitinase USP7 to modulate transcriptional processes in mESCs similar to MYC. EPOP is commonly upregulated in human cancer, and its loss impairs the proliferation of several human cancer cell lines. Our findings establish EPOP as a transcriptional modulator, which impacts both Polycomb and active gene transcription in mammalian cells.

The antiproliferative ELF2 isoform, ELF2B, induces apoptosis in vitro and perturbs early lymphocytic development in vivo

Background

ELF2 (E74-like factor 2) also known as NERF (new Ets-related factor), a member of the Ets family of transcription factors, regulates genes important in B and T cell development, cell cycle progression, and angiogenesis. Conserved ELF2 isoforms, ELF2A, and ELF2B, arising from alternative promoter usage can exert opposing effects on target gene expression. ELF2A activates, whilst ELF2B represses, gene expression, and the balance of expression between these isoforms may be important in maintaining normal cellular function.

Methods

We compared the function of ELF2 isoforms ELF2A and ELF2B with other ELF subfamily proteins ELF1 and ELF4 in primary and cancer cell lines using proliferation, colony-forming, cell cycle, and apoptosis assays. We further examined the role of ELF2 isoforms in haemopoietic development using a Rag1^-/-murine bone marrow reconstitution model.

Results

ELF2B overexpression significantly reduced cell proliferation and clonogenic capacity, minimally disrupted cell cycle kinetics, and induced apoptosis. In contrast, ELF2A overexpression only marginally reduced clonogenic capacity with little effect on proliferation, cell cycle progression, or apoptosis. Deletion of the N-terminal 19 amino acids unique to ELF2B abrogated the antiproliferative and proapoptotic functions of ELF2B thereby confirming its crucial role. Mice expressing Elf2a or Elf2b in haemopoietic cells variously displayed perturbations in the pre-B cell stage and multiple stages of T cell development. Mature B cells, T cells, and myeloid cells in steady state were unaffected, suggesting that the main role of ELF2 is restricted to the early development of B and T cells and that compensatory mechanisms exist. No differences in B and T cell development were observed between ELF2 isoforms.

Conclusions

We conclude that ELF2 isoforms are important regulators of cellular proliferation, cell cycle progression, and apoptosis. In respect to this, ELF2B acts in a dominant negative fashion compared to ELF2A and as a putative tumour suppressor gene. Given that these cellular processes are critical during haemopoiesis, we propose that the regulatory interplay between ELF2 isoforms contributes substantially to early B and T cell development.

Electronic supplementary material

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

Interactomic analysis of REST/NRSF and implications of its functional links with the transcription suppressor TRIM28 during neuronal differentiation

RE-1 silencing transcription factor (REST) is a transcriptional repressor that regulates gene expression by binding to repressor element 1. However, despite its critical function in physiology, little is known about its interaction proteins. Here we identified 204 REST-interacting proteins using affinity purification and mass spectrometry. The interactome included proteins associated with mRNA processing/splicing, chromatin organization, and transcription. The interactions of these REST-interacting proteins, which included TRIM28, were confirmed by co-immunoprecipitation and immunocytochemistry, respectively. Gene Ontology (GO) analysis revealed that neuronal differentiation-related GO terms were enriched among target genes that were co-regulated by REST and TRIM28, while the level of CTNND2 was increased by the knockdown of REST and TRIM28. Consistently, the level of CTNND2 increased while those of REST and TRIM28 decreased during neuronal differentiation in the primary neurons, suggesting that CTNND2 expression may be co-regulated by both. Furthermore, neurite outgrowth was increased by depletion of REST or TRIM28, implying that reduction of both REST and TRIM28 could promote neuronal differentiation via induction of CTNND2 expression. In conclusion, our study of REST reveals novel interacting proteins which could be a valuable resource for investigating unidentified functions of REST and also suggested functional links between REST and TRIM28 during neuronal development.

A Dynamic Role of TBX3 in the Pluripotency Circuitry

Pluripotency represents a cell state comprising a fine-tuned pattern of transcription factor activity required for embryonic stem cell (ESC) self-renewal. TBX3 is the earliest expressed member of the T-box transcription factor family and is involved in maintenance and induction of pluripotency. Hence, TBX3 is believed to be a key member of the pluripotency circuitry, with loss of TBX3 coinciding with loss of pluripotency. We report a dynamic expression of TBX3 in vitro and in vivo using genetic reporter tools tracking TBX3 expression in mouse ESCs (mESCs). Low TBX3 levels are associated with reduced pluripotency, resembling the more mature epiblast. Notably, TBX3-low cells maintain the intrinsic capability to switch to a TBX3-high state and vice versa. Additionally, we show TBX3 to be dispensable for induction and maintenance of naive pluripotency as well as for germ cell development. These data highlight novel facets of TBX3 action in mESCs.

Past Roadblocks and New Opportunities in Transcription Factor Network Mapping

One of the principal mechanisms by which cells differentiate and respond to changes in external signals or conditions is by changing the activity levels of transcription factors (TFs). This changes the transcription rates of target genes via the cell's TF network, which ultimately contributes to reconfiguring cellular state. Since microarrays provided our first window into global cellular state, computational biologists have eagerly attacked the problem of mapping TF networks, a key part of the cell's control circuitry. In retrospect, however, steady-state mRNA abundance levels were a poor substitute for TF activity levels and gene transcription rates. Likewise, mapping TF binding through chromatin immunoprecipitation proved less predictive of functional regulation and less amenable to systematic elucidation of complete networks than originally hoped. This review explains these roadblocks and the current, unprecedented blossoming of new experimental techniques built on second-generation sequencing, which hold out the promise of rapid progress in TF network mapping.

Generation and gene expression profiling of 48 transcription-factor-inducible mouse embryonic stem cell lines

Mouse embryonic stem cells (ESCs) can differentiate into a wide range – and possibly all cell types in vitro, and thus provide an ideal platform to study systematically the action of transcription factors (TFs) in cell differentiation. Previously, we have generated and analyzed 137 TF-inducible mouse ESC lines. As an extension of this “NIA Mouse ESC Bank,” we generated and characterized 48 additional mouse ESC lines, in which single TFs in each line could be induced in a doxycycline-controllable manner. Together, with the previous ESC lines, the bank now comprises 185 TF-manipulable ESC lines (>10% of all mouse TFs). Global gene expression (transcriptome) profiling revealed that the induction of individual TFs in mouse ESCs for 48 hours shifts their transcriptomes toward specific differentiation fates (e.g., neural lineages by Myt1 Isl1, and St18; mesodermal lineages by Pitx1, Pitx2, Barhl2, and Lmx1a; white blood cells by Myb, Etv2, and Tbx6, and ovary by Pitx1, Pitx2, and Dmrtc2). These data also provide and lists of inferred target genes of each TF and possible functions of these TFs. The results demonstrate the utility of mouse ESC lines and their transcriptome data for understanding the mechanism of cell differentiation and the function of TFs.

Single-Cell Analyses of ESCs Reveal Alternative Pluripotent Cell States and Molecular Mechanisms that Control Self-Renewal

Analyses of gene expression in single mouse embryonic stem cells (mESCs) cultured in serum and LIF revealed the presence of two distinct cell subpopulations with individual gene expression signatures. Comparisons with published data revealed that cells in the first subpopulation are phenotypically similar to cells isolated from the inner cell mass (ICM). In contrast, cells in the second subpopulation appear to be more mature. Pluripotency Gene Regulatory Network (PGRN) reconstruction based on single-cell data and published data suggested antagonistic roles for Oct4 and Nanog in the maintenance of pluripotency states. Integrated analyses of published genomic binding (ChIP) data strongly supported this observation. Certain target genes alternatively regulated by OCT4 and NANOG, such as Sall4 and Zscan10, feed back into the top hierarchical regulator Oct4. Analyses of such incoherent feedforward loops with feedback (iFFL-FB) suggest a dynamic model for the maintenance of mESC pluripotency and self-renewal.

•

Mouse embryonic stem cells grown on serum and LIF contain two subpopulations of cells

•

Oct4 and Nanog alternatively regulate a class of pluripotency genes

•

We demonstrate stabilization of Oct4 concentration and pluripotency via feedback control

•

The “state exchange” model explains self-renewal

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.

Systematic identification of gene family regulators in mouse and human embryonic stem cells

Pluripotent self-renewing embryonic stem cells (ESCs) have been the focus of a growing number of high-throughput experiments, revealing the genome-wide locations of hundreds of transcription factors and histone modifications. While most of these datasets were used in a specific context, all datasets combined offer a comprehensive view of chromatin characteristics and regulatory elements that govern cell states. Here, using hundreds of datasets in ESCs, we generated colocalization maps of chromatin proteins and modifications, and built a discovery pipeline for regulatory proteins of gene families. By comparing genome-wide binding data with over-expression and knockdown analysis of hundreds of genes, we discovered that the pluripotency-related factor NR5A2 separates mitochondrial from cytosolic ribosomal genes, regulating their expression. We further show that genes with a common chromatin profile are enriched for distinct Gene Ontology (GO) categories. Our approach can be generalized to reveal common regulators of any gene group; discover novel gene families, and identify common genomic elements based on shared chromatin features.

Emerging Modeling Concepts and Solutions in Stem Cell Research

Modern stem cell research, as well as other fields of contemporary biology involves quantitative sciences in many ways. Identifying candidates for key differentiation or reprogramming factors, tracing global transcriptome changes, or finding drugs is now broadly involves bioinformatics and biostatistics. However, the next key step, understanding the underlying reasons and establishing causal links leading to differentiation or reprogramming requires qualitative and quantitative biological models describing complex biological systems. Currently, quantitative modeling is a challenging science, capable to deliver rather modest results or predictions. What model types are the most popular and what features of stem cell behavior they are capturing? What new insights do we expect from the computational modeling of stem cells in the foreseeable future? Current review attempts to approach these essential questions by considering published quantitative models and solutions emerging in the area of stem cell research.

Messenger RNA profile analysis deciphers new Esrrb responsive genes in prostate cancer cells

Background

Orphan nuclear receptor estrogen related receptor β (Esrrb or ERRβ) is well known in stem cells and early embryonic development. However, little is known about its function in cancer.

Method

We investigated the mRNA profile alterations induced by Esrrb expression and its synthetic ligand DY131 in human prostate cancer DU145 cells via RNA-Seq analysis.

Results

We distinguished 67 mRNAs differentially expressed by Esrrb alone. Although DY131 alone did not change any gene, treatment of DY131 in the presence of Esrrb altered 1161 mRNAs. These observations indicated Esrrb had both ligand-independent and ligand-dependent activity. When Esrrb was expressed, DY131 treatment further regulated 15 Esrrb-altered mRNAs. DY131 acted as an antagonist for 11 of 15 mRNAs (wdr52, f13a1, pxdn, spns2, loc100506599, tagln, loc441454, tkel1, sema3f, zcwpw2, sdc2) and as an agonist for 4 of the 15 mRNAs (rarres3, oasl, padi2, ddx60). Gene ontology analyses showed altered genes are related to transcription and translation regulation, cell proliferation and apoptosis regulation, and cellular metabolism.

Conclusion

Our results characterized mRNA profiles in DU145 prostate cancer cells driven by Esrrb expression and Esrrb ligand DY131, and provided multiple markers to characterize Esrrb’s function in Esrrb research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12867-015-0049-1) contains supplementary material, which is available to authorized users.

Genes targeted by the Hedgehog-signaling pathway can be regulated by Estrogen related receptor β

BACKGROUND

Nuclear receptor family member, Estrogen related receptor β, and the Hedgehog signal transduction pathway are both reported to relate to tumorigenesis and induced pluripotent stem cell reprogramming. We hypothesize that Estrogen related receptor β can modulate the Hedgehog signaling pathway and affect Hedgehog driven downstream gene expression.

RESULTS

We established an estrogen related receptor β-expressing Hedgehog-responsive NIH3T3 cell line by Esrrb transfection, and performed mRNA profiling using RNA-Seq after Hedgehog ligand conditioned medium treatment. Esrrb expression altered 171 genes, while Hedgehog signaling activation alone altered 339 genes. Additionally, estrogen related receptor β expression in combination with Hedgehog signaling activation affects a group of 109 Hedgehog responsive mRNAs, including Hsd11b1, Ogn, Smoc2, Igf1, Pdcd4, Igfbp4, Stmn1, Hp, Hoxd8, Top2a, Tubb4b, Sfrp2, Saa3, Prl2c3 and Dpt.

CONCLUSIONS

We conclude that Estrogen related receptor β is capable of interacting with Hh-signaling downstream targets. Our results suggest a new level of regulation of Hedgehog signaling by Estrogen related receptor β, and indicate modulation of Estrogen related receptor β can be a new strategy to regulate various functions driven by the Hedgehog signaling pathway.

Systems Analyses Reveal Shared and Diverse Attributes of Oct4 Regulation in Pluripotent Cells

We combine a genome-scale RNAi screen in mouse epiblast stem cells (EpiSCs) with genetic interaction, protein localization, and "protein-level dependency" studies-a systematic technique that uncovers post-transcriptional regulation-to delineate the network of factors that control the expression of Oct4, a key regulator of pluripotency. Our data signify that there are similarities, but also fundamental differences in Oct4 regulation in EpiSCs versus embryonic stem cells (ESCs). Through multiparametric data analyses, we predict that Tox4 is associating with the Paf1C complex, which maintains cell identity in both cell types, and validate that this protein-protein interaction exists in ESCs and EpiSCs. We also identify numerous knockdowns that increase Oct4 expression in EpiSCs, indicating that, in stark contrast to ESCs, Oct4 is under active repressive control in EpiSCs. These studies provide a framework for better understanding pluripotency and for dissecting the molecular events that govern the transition from the pre-implantation to the post-implantation state.

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.

ChIP-Array 2: integrating multiple omics data to construct gene regulatory networks

Transcription factors (TFs) play an important role in gene regulation. The interconnections among TFs, chromatin interactions, epigenetic marks and cis-regulatory elements form a complex gene transcription apparatus. Our previous work, ChIP-Array, combined TF binding and transcriptome data to construct gene regulatory networks (GRNs). Here we present an enhanced version, ChIP-Array 2, to integrate additional types of omics data including long-range chromatin interaction, open chromatin region and histone modification data to dissect more comprehensive GRNs involving diverse regulatory components. Moreover, we substantially extended our motif database for human, mouse, rat, fruit fly, worm, yeast and Arabidopsis, and curated large amount of omics data for users to select as input or backend support. With ChIP-Array 2, we compiled a library containing regulatory networks of 18 TFs/chromatin modifiers in mouse embryonic stem cell (mESC). The web server and the mESC library are publicly free and accessible athttp://jjwanglab.org/chip-array.

Inferring dynamic gene regulatory networks in cardiac differentiation through the integration of multi-dimensional data

Background

Decoding the temporal control of gene expression patterns is key to the understanding of the complex mechanisms that govern developmental decisions during heart development. High-throughput methods have been employed to systematically study the dynamic and coordinated nature of cardiac differentiation at the global level with multiple dimensions. Therefore, there is a pressing need to develop a systems approach to integrate these data from individual studies and infer the dynamic regulatory networks in an unbiased fashion.

Results

We developed a two-step strategy to integrate data from (1) temporal RNA-seq, (2) temporal histone modification ChIP-seq, (3) transcription factor (TF) ChIP-seq and (4) gene perturbation experiments to reconstruct the dynamic network during heart development. First, we trained a logistic regression model to predict the probability (LR score) of any base being bound by 543 TFs with known positional weight matrices. Second, four dimensions of data were combined using a time-varying dynamic Bayesian network model to infer the dynamic networks at four developmental stages in the mouse [mouse embryonic stem cells (ESCs), mesoderm (MES), cardiac progenitors (CP) and cardiomyocytes (CM)]. Our method not only infers the time-varying networks between different stages of heart development, but it also identifies the TF binding sites associated with promoter or enhancers of downstream genes.

The LR scores of experimentally verified ESCs and heart enhancers were significantly higher than random regions (p <10⁻¹⁰⁰), suggesting that a high LR score is a reliable indicator for functional TF binding sites. Our network inference model identified a region with an elevated LR score approximately −9400 bp upstream of the transcriptional start site of Nkx2-5, which overlapped with a previously reported enhancer region (−9435 to −8922 bp). TFs such as Tead1, Gata4, Msx2, and Tgif1 were predicted to bind to this region and participate in the regulation of Nkx2-5 gene expression. Our model also predicted the key regulatory networks for the ESC-MES, MES-CP and CP-CM transitions.

Conclusion

We report a novel method to systematically integrate multi-dimensional -omics data and reconstruct the gene regulatory networks. This method will allow one to rapidly determine the cis-modules that regulate key genes during cardiac differentiation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0460-0) contains supplementary material, which is available to authorized users.

GATA2 as a potential metastasis-driving gene in prostate cancer

Effective treatment for metastatic prostate cancer is critically needed. The present study was aimed at identifying metastasis-driving genes as potential targets for therapy (oncotargets). A differential gene expression profile of metastatic LTL-313H and non-metastatic LTL-313B prostate cancer tissue xenografts, derived from one patient's specimen, was subjected to integrative analysis using the Ingenuity Upstream Regulator Analysis tool. Six candidate master regulatory genes were identified, including GATA2, a gene encoding a pioneer factor, a special transcription factor facilitating the recruitment of additional transcription factors. Elevated GATA2 expression in metastatic prostate cancer tissues correlated with poor patient prognosis. Furthermore, GATA2 gene silencing in human prostate cancer LNCaP cells led to a marked reduction in cell migration, tissue invasion, focal adhesion disassembly and to a dramatic change in cell transcriptomes, indicating that GATA2 plays a critical role in prostate cancer metastasis. As such, GATA2 could represent a prostate cancer metastasis-driving gene and a potential target for therapy of metastatic prostate cancer.

SOX9 accelerates ESC differentiation to three germ layer lineages by repressing SOX2 expression through P21 (WAF1/CIP1)

Upon removal of culture conditions that maintain an undifferentiated state, mouse embryonic stem cells (ESCs) differentiate into various cell types. Differentiation can be facilitated by forced expression of certain transcription factors (TFs), each of which can generally specify a particular developmental lineage. We previously established 137 mouse ESC lines, each of which carried a doxycycline-controllable TF. Among them, Sox9 has unique capacity: its forced expression accelerates differentiation of mouse ESCs into cells of all three germ layers. With the additional use of specific culture conditions, overexpression of Sox9 facilitated the generation of endothelial cells, hepatocytes and neurons from ESCs. Furthermore, Sox9 action increases formation of p21 (WAF1/CIP1), which then binds to the SRR2 enhancer of pluripotency marker Sox2 and inhibits its expression. Knockdown of p21 abolishes inhibition of Sox2 and Sox9-accelerated differentiation, and reduction of Sox2 2 days after the beginning of ESC differentiation can comparably accelerate mouse ESC formation of cells of three germ layers. These data implicate the involvement of the p21-Sox2 pathway in the mechanism of accelerated ESC differentiation by Sox9 overexpression. The molecular cascade could be among the first steps to program ESC differentiation.

Primate-specific endogenous retrovirus-driven transcription defines naive-like stem cells

Naive embryonic stem cells hold great promise for research and therapeutics as they have broad and robust developmental potential. While such cells are readily derived from mouse blastocysts it has not been possible to isolate human equivalents easily, although human naive-like cells have been artificially generated (rather than extracted) by coercion of human primed embryonic stem cells by modifying culture conditions or through transgenic modification. Here we show that a sub-population within cultures of human embryonic stem cells (hESCs) and induced pluripotent stem cells (hiPSCs) manifests key properties of naive state cells. These naive-like cells can be genetically tagged, and are associated with elevated transcription of HERVH, a primate-specific endogenous retrovirus. HERVH elements provide functional binding sites for a combination of naive pluripotency transcription factors, including LBP9, recently recognized as relevant to naivety in mice. LBP9-HERVH drives hESC-specific alternative and chimaeric transcripts, including pluripotency-modulating long non-coding RNAs. Disruption of LBP9, HERVH and HERVH-derived transcripts compromises self-renewal. These observations define HERVH expression as a hallmark of naive-like hESCs, and establish novel primate-specific transcriptional circuitry regulating pluripotency.

Dax1 and Nanog act in parallel to stabilize mouse embryonic stem cells and induced pluripotency

Nanog expression is heterogeneous and dynamic in embryonic stem cells (ESCs). However, the mechanism for stabilizing pluripotency during the transitions between Nanog(high) and Nanog(low) states is not well understood. Here we report that Dax1 acts in parallel with Nanog to regulate mouse ESC (mESCs) identity. Dax1 stable knockdown mESCs are predisposed towards differentiation but do not lose pluripotency, whereas Dax1 overexpression supports LIF-independent self-renewal. Although partially complementary, Dax1 and Nanog function independently and cannot replace one another. They are both required for full reprogramming to induce pluripotency. Importantly, Dax1 is indispensable for self-renewal of Nanog(low) mESCs. Moreover, we report that Dax1 prevents extra-embryonic endoderm (ExEn) commitment by directly repressing Gata6 transcription. Dax1 may also mediate inhibition of trophectoderm differentiation independent or as a downstream effector of Oct4. These findings establish a basal role of Dax1 in maintaining pluripotency during the state transition of mESCs and somatic cell reprogramming.

BRD4 regulates Nanog expression in mouse embryonic stem cells and preimplantation embryos

Bromodomain-containing protein 4 (BRD4) is an important epigenetic reader implicated in the pathogenesis of a number of different cancers and other diseases. Brd4-null mouse embryos die shortly after implantation and are compromised in their ability to maintain the inner cell mass, which gives rise to embryonic stem cells (ESCs). Here we report that BRD4 regulates expression of the pluripotency factor Nanog in mouse ESCs and preimplantation embryos, as well as in human ESCs and embryonic cancer stem cells. Inhibition of BRD4 function using a chemical inhibitor, small interfering RNAs, or a dominant-negative approach suppresses Nanog expression, and abolishes the self-renewal ability of ESCs. We also find that BRD4 associates with BRG1 (brahma-related gene 1, aka Smarca4 (SWI/SNF-related, matrix-associated, actin-dependent regulator of chromatin, subfamily a, member 4)), a key regulator of ESC self-renewal and pluripotency, in the Nanog regulatory regions to regulate Nanog expression. Our study identifies Nanog as a novel BRD4 target gene, providing new insights for the biological function of BRD4 in stem cells and mouse embryos. Knowledge gained from these non-cancerous systems will facilitate future investigations of how Brd4 dysfunction leads to cancers.

Chromatin properties of regulatory DNA probed by manipulation of transcription factors

Transcription factors (TFs) bind to DNA and regulate the transcription of nearby genes. However, only a small fraction of TF binding sites have such regulatory effects. Here we search for the predictors of functional binding sites by carrying out a systematic computational screening of a variety of contextual factors (histone modifications, nuclear lamin-bindings, and cofactor bindings). We used regression analysis to test if contextual factors are associated with upregulation or downregulation of neighboring genes following the induction or knockdown of the 9 TFs in mouse embryonic stem (ES) cells. Functional TF binding sites appeared to be either active (i.e., bound by P300, CHD7, mediator, cohesin, and SWI/SNF) or repressed (i.e., with H3K27me3 histone marks and bound by Polycomb factors). Active binding sites mediated the downregulation of nearby genes upon knocking down the activating TFs or inducing repressors. Repressed TF binding sites mediated the upregulation of nearby genes (e.g., poised developmental regulators) upon inducing TFs. In addition, repressed binding sites mediated repressive effects of TFs, identified by the downregulation of target genes after the induction of TFs or by the upregulation of target genes after the knockdown of TFs. The contextual factors associated with functions of DNA-bound TFs were used to improve the identification of candidate target genes regulated by TFs.

Defining an essential transcription factor program for naïve pluripotency

The gene regulatory circuitry through which pluripotent embryonic stem (ES) cells choose between self-renewal and differentiation appears vast and has yet to be distilled into an executive molecular program. We developed a data-constrained, computational approach to reduce complexity and to derive a set of functionally validated components and interaction combinations sufficient to explain observed ES cell behavior. This minimal set, the simplest version of which comprises only 16 interactions, 12 components, and three inputs, satisfies all prior specifications for self-renewal and furthermore predicts unknown and nonintuitive responses to compound genetic perturbations with an overall accuracy of 70%. We propose that propagation of ES cell identity is not determined by a vast interactome but rather can be explained by a relatively simple process of molecular computation.

A genome scale overexpression screen to reveal drug activity in human cells

Target identification is a critical step in the lengthy and expensive process of drug development. Here, we describe a genome-wide screening platform that uses systematic overexpression of pooled human ORFs to understand drug mode-of-action and resistance mechanisms. We first calibrated our screen with the well-characterized drug methotrexate. We then identified new genes involved in the bioactivity of diverse drugs including antineoplastic agents and biologically active molecules. Finally, we focused on the transcription factor RHOXF2 whose overexpression conferred resistance to DNA damaging agents. This approach represents an orthogonal method for functional screening and, to our knowledge, has never been reported before.

Decoding the Pluripotency Network: The Emergence of New Transcription Factors

Since the successful isolation of mouse and human embryonic stem cells (ESCs) in the past decades, massive investigations have been conducted to dissect the pluripotency network that governs the ability of these cells to differentiate into all cell types. Beside the core Oct4-Sox2-Nanog circuitry, accumulating regulators, including transcription factors, epigenetic modifiers, microRNA and signaling molecules have also been found to play important roles in preserving pluripotency. Among the various regulations that orchestrate the cellular pluripotency program, transcriptional regulation is situated in the central position and appears to be dominant over other regulatory controls. In this review, we would like to summarize the recent advancements in the accumulating findings of new transcription factors that play a critical role in controlling both pluripotency network and ESC identity.

Identification of transcription factors for lineage-specific ESC differentiation

A network of transcription factors (TFs) determines cell identity, but identity can be altered by overexpressing a combination of TFs. However, choosing and verifying combinations of TFs for specific cell differentiation have been daunting due to the large number of possible combinations of ∼2,000 TFs. Here, we report the identification of individual TFs for lineage-specific cell differentiation based on the correlation matrix of global gene expression profiles. The overexpression of identified TFs—Myod1, Mef2c, Esx1, Foxa1, Hnf4a, Gata2, Gata3, Myc, Elf5, Irf2, Elf1, Sfpi1, Ets1, Smad7, Nr2f1, Sox11, Dmrt1, Sox9, Foxg1, Sox2, or Ascl1—can direct efficient, specific, and rapid differentiation into myocytes, hepatocytes, blood cells, and neurons. Furthermore, transfection of synthetic mRNAs of TFs generates their appropriate target cells. These results demonstrate both the utility of this approach to identify potent TFs for cell differentiation, and the unanticipated capacity of single TFs directly guides differentiation to specific lineage fates.

•

Lineage-determining single TFs are identified based on the correlation matrix

•

A proof of concept is demonstrated for ESC differentiation by 21 TFs

•

TFs orchestrate global gene expression changes via direct binding to target genes

•

Transfections of synthetic TF mRNAs generate desired differentiated cells

A competitive protein interaction network buffers Oct4-mediated differentiation to promote pluripotency in embryonic stem cells

The dynamic competition for complex formation between the pluripotency network components Oct4, Nanog, Tcf3, and β-catenin prevents embryonic stem cell differentiation by controlling the levels of free Oct4.

Pluripotency is defined by the ratios between the levels of pluripotency factors rather than by their absolute levels.

Competition between different protein complexes involving Nanog, Oct4, Tcf3, and β-catenin can account for the ratios associated with pluripotency.

The unstable pluripotency of Nanog mutant cells was shown to depend on the interactions between Oct4 and β-catenin.

The function of the protein competition network is to control the levels of free Oct4, which are balanced by Nanog and β-catenin in embryonic stem cells.

Pluripotency in embryonic stem cells is maintained through the activity of a small set of transcription factors centred around Oct4 and Nanog, which control the expression of ‘self-renewal' and ‘differentiation' genes. Here, we combine single-cell quantitative immunofluorescence microscopy and gene expression analysis, together with theoretical modelling, to investigate how the activity of those factors is regulated. We uncover a key role for post-translational regulation in the maintenance of pluripotency, which complements the well-established transcriptional regulatory layer. Specifically, we find that the activity of a network of protein complexes involving Nanog, Oct4, Tcf3, and β-catenin suffices to account for the behavior of ES cells under different conditions. Our results suggest that the function of the network is to buffer the transcriptional activity of Oct4, which appears to be the main determinant to exit pluripotency. The protein network explains the mechanisms underlying the gain and loss of function in different mutants, and brings us closer to a full understanding of the molecular basis of pluripotency.

Identification of the missing pluripotency mediator downstream of leukaemia inhibitory factor

Self-renewal of pluripotent mouse embryonic stem (ES) cells is sustained by the cytokine leukaemia inhibitory factor (LIF) acting through the transcription factor Stat3. Several targets of Stat3 have previously been identified, most notably the reprogramming factor Klf4. However, such factors are neither required nor sufficient for the potent effect of LIF. We took advantage of Stat3 null ES cells to confirm that Stat3 mediates the self-renewal response to LIF. Through comparative transcriptome analysis intersected with genome location data, we arrived at a set of candidate transcription factor effectors. Among these, Tfcp2l1 (also known as Crtr-1) was most abundant. Constitutive expression of Tfcp2l1 at levels similar to those induced by LIF effectively substituted for LIF or Stat3 in sustaining clonal self-renewal and pluripotency. Conversely, knockdown of Tfcp2l1 profoundly compromised responsiveness to LIF. We further found that Tfcp2l1 is both necessary and sufficient to direct molecular reprogramming of post-implantation epiblast stem cells to naïve pluripotency. These results establish Tfcp2l1 as the principal bridge between LIF/Stat3 input and the transcription factor core of naïve pluripotency.