# Regulation of Cyclin E by transcription factors of the naïve pluripotency network in mouse embryonic stem cells

Fabrice Gonnot, Diana Langer, Pierre-Yves Bourillot, Nathalie Doerflinger, Pierre Savatier

► To cite this version:

Fabrice Gonnot, Diana Langer, Pierre-Yves Bourillot, Nathalie Doerflinger, Pierre Savatier. Regulation of Cyclin E by transcription factors of the naïve pluripotency network in mouse embryonic stem cells. Cell Cycle, 2019, 18 (20), pp.2697-2712. 10.1080/15384101.2019.1656475. hal-02322526

HAL Id: hal-02322526

https://hal.science/hal-02322526v1

Submitted on 3 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH PAPER

<figure><img src="image_2.png" /></figure>

OPEN ACCESS

<figure><img src="image_3.png" /></figure>

Check for updates

# Regulation of Cyclin E by transcription factors of the naïve pluripotency network in mouse embryonic stem cells

Fabrice Gonnot, Diana Langer, Pierre-Yves Bourillot, Nathalie Doerflinger, and Pierre Savatier

Stem Cell and Brain Research Institute, Univ Lyon, Université Lyon 1, Inserm, Bron, France

## ABSTRACT

Continuous, non-cell cycle-dependent expression of cyclin E is a characteristic feature of mouse embryonic stem cells (mESCs). We studied the 5' regulatory region of *Cyclin E*, also known as *Ccne1*, and identified binding sites for transcription factors of the naïve pluripotency network, including *Esrrb*, *Klf4*, and *Tfcp2l1* within 1 kilobase upstream of the transcription start site. Luciferase assay and chromatin immunoprecipitation-quantitative polymerase chain reaction (ChiP-qPCR) study highlighted one binding site for *Esrrb* that is essential to transcriptional activity of the promoter region, and three binding sites for *Klf4* and *Tfcp2l1*. Knockdown of *Esrrb*, *Klf4*, and *Tfcp2l1* reduced *Cyclin E* expression whereas overexpression of *Esrrb* and *Klf4* increased it, indicating a strong correlation between the expression level of these factors and that of cyclin E. We observed that cyclin E overexpression delays differentiation induced by *Esrrb* depletion, suggesting that cyclin E is an important target of *Esrrb* for differentiation blockade. We observed that mESCs express a low level of *miR-15a* and that transfection of a *miR-15a* mimic decreases *Cyclin E* mRNA level. These results lead to the conclusion that the high expression level of *Cyclin E* in mESCs can be attributed to transcriptional activation by *Esrrb* as well as to the absence of its negative regulator, *miR-15a*.

## ARTICLE HISTORY

Received 22 February 2019

Revised 26 July 2019

Accepted 2 August 2019

## KEYWORDS

Embryonic stem cell;
pluripotent; cell-cycle; Cyclin E

## Introduction

Cyclin E is a regulatory subunit of cyclin-dependent kinase (Cdk) 2 involved in many cellular processes including cell cycle progression, replication complex assembly, centrosome cycle, and epigenetic regulation. Its expression is regulated at both the transcriptional and protein level to achieve a timely control of cell division in connection with cell environment and fate decision. Deregulated expression of cyclin E has been shown to play a key role in tumorigenesis [1,2]. Transcription of the *Cyclin E* gene, also known as *Ccne1*, is activated during the G₁ phase and depends on mitogenic input, which is integrated through E2F and Myc transcription factors [3,4]. E2F activity is regulated by phosphorylation of the retinoblastoma (Rb) protein in response to cyclin D/Cdk4 and cyclin D/Cdk6 kinase activity [5,6]. The miRNA *miR-15a* was shown to act as a negative regulator of *Cyclin E* in somatic cells [7,8]. Since, both *Cyclin E* and *mir-15a* are direct transcriptional targets of E2F, it raises the possibility that E2F, *miR-15a*, and cyclin E constitute

a feed-forward loop that modulates E2F activity and cell-cycle progression [8].

There is a growing body of evidence showing that the cell cycle of mouse embryonic stem cells (mESCs) lacks some of the regulatory pathways that operate in somatic cells [9–11]. These include extensive phosphorylation of the Rb family proteins despite little cyclin D/Cdk4 kinase activity [12], $p16$<sup>ink4a</sup>-resistant residual cyclin D3/Cdk6 kinase activity [13], and lack of functional Chk/ $p53/p21$<sup>cip1</sup> and Chk/Cdc25A pathways resulting in the absence of the DNA damage checkpoint in the G₁ phase [14–16]. A key feature of the pluripotent stem cell cycle is the constitutive activity of Cdk2 due to seemingly continuous expression of both cyclin E and A throughout the cell cycle [17,18] in addition to low expression levels of the Cdk2 inhibitors $p21$<sup>cip1</sup>, $p27$<sup>kip1</sup>, and $p57$<sup>kip2</sup> [12,17]. In a previous report, we showed that cyclin E partially rescues mESC differentiation induced by leukemia inhibitory factor (LIF) starvation, suggesting that cyclin E participates in the regulation of pluripotency [19]. It was established

## CONTACT

<figure><img src="image_4.png" /></figure>

Pierre Savatier pierre.savatier@inserm.fr

<figure><img src="image_5.png" /></figure>

Supplemental data for this article can be accessed on the publisher's website.

© 2019 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.

This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives License (http://creativecommons.org/licenses/by-nc-nd/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited, and is not altered, transformed, or built upon in any way.
<figure><img src="image_6.png" /></figure>

that cyclin E:Cdk2 complexes phosphorylate and thereby stabilize the core pluripotency factors Nanog, Sox2, and Oct4 [20]. These findings point to a connection between the cell cycle machinery regulating G₁/S phase transition and the core pluripotency network [21].

In this context, it is important to understand how *Cyclin E* is transcriptionally regulated in pluripotent stem cells. We hypothesized that the transcription factors of the naïve pluripotency network would participate in the transcriptional regulation of *Cyclin E*. These factors include the cardinal pluripotency factors Oct4, Sox2, and Nanog, as well as the ancillary transcription factors Klf2, Klf4, Klf5, Esrrb, Tbx3, Gbx2, Nr0b1, and Tfcp2l1, all of which have been shown to sustain the naïve state of pluripotency in mice [22–34]. The present study points to Esrrb as a transcriptional activator and miR-15a as a negative regulator of *Cyclin E* in mESCs.

# Material and methods

## In silico analysis

Published data were obtained from NCBI Gene Expression Omnibus (http://www.ncbi.nlm.nih. gov/geo) and analyzed using UCSC Genome Browser [35; http://genome.ucsc.edu]. DNAse I hypersensitive sites, were identified from GSM1003830 (DNAseDgf on mESC-CJ7), GSM1014154 (DNAseHS on mESC-E14), and GSM1014187 (DNAseHS on mESC-CJ7) datasets. Histone marks were identified from GSM769008 (H3K4me3 on mESC-Bruce4), GSM1000089 (H3K27me3 on mESC-Bruce4) and GSM1000124 (H3K4me3 on mESC-E14) datasets. ChIP-seq data were from GSM288345 (Nanog), GSM288346 (Oct4), GSM288347 (Sox2), GSM288349 (E2f1), GSM288350 (Tfcp2l1), GSM288353 (Stat3), GSM288354 (Klf4), GSM288355 (Esrrb), and GSM288356 (c-Myc) compendiums [36], and GSM470523 (Nr5a2) [37] and GSM1208217 (Klf4) [38]. Several resources were used to predict the transcription factor binding site (TFBS)'s relative scores on the genomic sequence upstream of the *Ccne1* gene, downloaded from the Ensembl database (genome assembly GRCm38/mm10, December 2011). They include JASPAR [39; http://jaspar.genereg.net], TRANSFAC 7.0 public

by BIOBASE [40; http://www.gene-regulation. com], MAPPER2 [41; http://genome.ufl.edu/map perdb], cisRED mouse v4 [42; http://www.cisred. org/mouse4], UniPROBE [43; http://the_brain. bwh.harvard.edu/uniprobe], MotifViz [44; http:// biowulf.bu.edu/MotifViz] and CONSITE [45; http://consite.genereg.net]. A transcription factor and DNA sequence matching degree greater than 80% was considered as a putative TFBS.

## Quantitative real-time PCR (qRT-PCR)

Total RNA was isolated from cell pellets using TRIzol (Ambion) according to the manufacturer's protocol and reverse-transcribed using a High-Capacity RNA-to-cDNA kit (Applied Biosystems). For microRNAs reverse-transcription, a stem-loop primer specific to each miRNA was used. Real-time PCR was performed using the StepOnePlus real-time PCR system (Applied Biosystems) and Fast SYBR Green Master Mix (Applied Biosystems) according to the manufacturer's instructions. The relative quantitation of gene expression was calculated using StepOne Software 2.3 (Applied Biosystems). Expression of the target genes was normalized to those of the mouse *β-actin* gene (*Actb*) or to the mouse *sno234* RNA for miRNA. Primers are listed in Table S1.

## ChIP-PCR

ChIP for Esrrb, Klf4, and Tfcp2l1 was performed on E14Tg2a mESCs using previously described protocols [46]. In brief, 10⁷ cells were cross-linked with 1% formaldehyde for 15 min. Chromatin was sonicated to a length of less than 400 bp, and subsequently immunoprecipitated with 5 µg of anti-Esrrb (Perseus, pp-H6705-00), anti-Klf4 (Stemgent, 09-0021), and anti-Tfcp2l1 (AbCam, ab123354). DNA fragments encompassing binding sites for Esrrb, Klf4, and Tfcp2l1 in the P region of *Cyclin E* and the *Nanog* promoters were subsequently amplified by qPCR. A 3' untranslated region of the *Cyclin E* gene lacking putative binding sites for Esrrb, Klf4, and Tfcp2l1 was used as negative control. Primers are listed in Table S2. ChIP-qPCR data obtained for each specific antibody were normalized using the percent input method that normalizes according to the
<figure><img src="image_7.png" /></figure>

amount of chromatin input. The percentage value for each sample was calculated based on the equation as follows: % Input = 100 x [primer pair efficiency]^(Ct[adjusted input] – Ct[IP]). The “% Input” value represents the enrichment of factor on specific region.

## Plasmid constructs

Regions I (P) (1.5 kb), II (PE) (1.5 kb), and III (DE) (1.7 kb) of the *Cyclin E* gene 5' flanking sequence were synthesized by GeneArt (Invitrogen) with appropriate restriction sites at both ends and subcloned into the pMA plasmid to generate *pMA-P*, *pMA-PE*, and *pMA-DE* plasmids, respectively (Table S3). A 1,512 base pair (bp) BglII-HindIII fragment encompassing region P was prepared from *pMA-P* and subcloned between BglII and Hind III in *pGL4.10[luc2]* (Promega) to generate *pGL4.10-P*. A 1,506 bp MluI fragment encompassing region PE was prepared from *pMA-PE* and subcloned into the MluI site in *pGL4.10-P* to generate *pGL4.10-P+ PE*. A 1,706 bp EcoRV-BglII fragment was subcloned between EcoRV and BglII sites in *pGL4.10-P* and *pGL4.10-P+ PE* to generate *pGL4.10-P+ DE* and *pGL4.10-P+ PE+DE*, respectively.

For site-directed mutagenesis of *Esrrb*, *Klf4*, and *Tfcp2l1* binding sites, the *pGL4.10-P* plasmid was mutated by PCR using mutant primers and Q5 site-directed mutagenesis kit (New England Biolabs, E0554) according to the manufacturer's instructions. Mutant primers were designed to modify 10-14 bp encompassing *Esrrb*, *Klf4*, and *Tfcp2l1* binding sites into 10-14 bp sequences with low binding scores using JASPAR software (Table S4).

The shRNA sequences used to knock down *Esrrb*, *Klf4*, and *Tfcp2l1* are described in Table S5. An expression cassette containing the *eGFP* reporter was synthesized by GeneArt (Invitrogen) with MfeI restriction sites at either ends and subcloned into pMA to generate *pMA-eGFP-miR30* plasmid as previously described [47]. shRNA sequences for *Esrrb*, *Klf4*, *Tfcp2l1*, and a control sequence were subsequently introduced between XhoI and EcoRI restriction sites in *pMA-eGFP-miR30* to generate *pMA-eGFP-miR30-shEsrrb*, *pMA-eGFP-miR30-shKlf4*, *pMA-eGFP-miR30-shTfcp2l1*, and *pMA-eGFP-miR30-shControl*, respectively. MfeI fragments

containing the seven *eGFP-miR30-shRNA* sequences were subsequently subcloned into the EcoRI site in *pBS31* [48] to generate *pBS31-TetON-shEsrrb#1*, *pBS31-TetON-shEsrrb#2*, *pBS31-TetON-shKlf4#1*, *pBS31-TetON-shKlf4#2*, *pBS31-TetON-shTfcp2l1#1*, *pBS31-TetON-shTfcp2l1#2*, and *pBS31-TetON-shControl*, respectively.

For cDNA overexpression, the coding sequences of mouse *Esrrb* (*NM_011934*), *Klf4* (*NM_010637*), and *Tfcp2l1* (*NM_023755*) were amplified by PCR with Q5 DNA polymerase (NEB, M0493) and primers containing an EcoRI restriction site and subsequently subcloned into the unique EcoRI restriction site in pBS31 [48] to generate *pBS31-TetON-Esrrb*, *pBS31-TetONhKlf4*, and *pBS31-TetON-Tfcp2l1*. For *Tfcp2l1*, before pBS31 subcloning, the PCR product was cloned into pJET1.2 with the CloneJET PCR Cloning Kit (Thermo Scientific, K1231), according to the manufacturer's instructions, and a silent mutation of the EcoRI site, present at position 523-528 of the sequence was performed (GAATTC in to GAGTTC) according to the previously described protocol.

## Cell culture, generation of stable transfectants, and colony assay

Parental E14Tg2a, E14Tg2a-Fucci, KH2 [47], and EKOiE [49] mESC lines were routinely cultured on 0.1% gelatin-coated dishes in Glasgow's Modified Eagle Medium supplemented with 10% fetal calf serum, 100 µM nonessential amino acids, 100 U/ mL penicillin, 100 µg/mL streptomycin, 2 mM L-glutamine, 1 mM sodium pyruvate, 100 µM 2-mercaptoethanol, and 1000 U/mL LIF. Routine culture of EKOiE mESCs included 1 µg/mL Doxycycline [49]. Differentiation of E14Tg2a-Fucci cells was induced after withdrawal of LIF for 5 days. Protocol for routine culture of epiblast stem cells (EpiSC) is described elsewhere [50].

For generation of stable transfectants, 1 × 10⁶ KH2 mESCs were electroporated with 5 µg of *pBS31-TetON* plasmid and 5 µg of *pCAGgs-FLPe* plasmid using the Neon system (Invitrogen) with two impulses (20 ms, 1300 volts). After 48 h, stable transfectants were selected with 40 µg/mL Hygromycin B (Roche, 10,843,555,001). cDNA and shRNA expression were induced by Doxycycline (Sigma, D9891) at concentrations ranging from 0.1
<figure><img src="image_8.png" /></figure>

to 1.0 µg/mL. EKOiE-MycCyclin E and EKOiE-control cells were generated by electroporating EKOiE cells with 10 µg of *pCAGgs*-MycCyclin E plasmid [19] followed by Hygromycin selection.

For colony assays, mESCs were plated at a density of $10^3$ cells per gelatin-coated 100-mm tissue culture dish in complete mESC medium. Cells were exposed to the medium without doxycycline for 1 to 7 days. The protocol for *in situ* detection of alkaline phosphatase activity is described elsewhere [19].

## Infection with lentivirus vectors, flow cytometry, cell transfection and luciferase assay

Production of simian immunodeficiency virus (SIV)-derived lentivectors expressing the Fucci reporters *mKO2:Cdt1*, *mAG:Geminin* and infection of mESCs are described elsewhere [19]. Cells were either analyzed using a LSRFortessa X-20 (Becton-Dickinson), or sorted with a FACS Aria cell sorter (Becton-Dickinson) as described in [19]. For transient expression assay, $5 \times 10^4$ E14Tg2a cells were transfected with 100 ng of reporter plasmids (*pGL4.10*[luc2] and their derivatives) and 1 ng of *pGL4.70*[hRluc] control plasmid using Lipofectamine 2000 (Invitrogen) in 96-well plates. Luciferase activity was measured after 48 h using the Dual-Glo Luciferase Assay System kit (Promega) and the GloMax Multi Detection System (Promega). For microRNAs mimics transient expression assay, 25 mM of miR-15a-5p (Ambion, 4,464,066 - MC10235) and miR-1 positive control (Ambion, 4,464,065) mimics were transfected and the gene expression was measured after 48 h.

## Immunoblotting and immunolabeling

For immunoblotting, frozen cell pellets were lysed in RIPA buffer complemented with protease and phosphatase inhibitors. Protein lysates were then cleared by centrifugation (17,000 × g for 20 min). After SDS-PAGE and electroblotting on polyvinylidene fluoride, the membranes were incubated with specific primary antibodies (mouse anti-Esrrb, Perseus PP-H6705-00; rabbit anti-Klf4, Santa Cruz, sc-20,691; rabbit anti-Tfcp2l1, AbCam, ab123354; anti-βactin, Sigma, A3854). Blots were incubated with

horseradish peroxidase (HRP)-coupled sheep anti-mouse IgG (GE Healthcare, NA931VS) and (HRP)-coupled goat anti-rabbit IgG (GE Healthcare, NA934VS), and developed with Clarity Western ECL Substrate (Bio-Rad, 1,705,060).

For immunolabeling, cells were fixed with 2% paraformaldehyde in PBS at 4°C for 20 min, and permeabilized in Tris-buffered saline (TBS; 50 mM Tris [pH 7.6], 0.9% NaCl, and 0.2% Triton X-100). The cells were then incubated overnight at 4°C with primary antibodies [anti-cyclin E1 rabbit polyclonal, Santa Cruz, sc-481 (1/100 dilution)]; anti-mKO2 mouse IgG1 (1/200 dilution), Clinisciences, M168-3M; anti-Esrrb mouse IgG2a (1/500 dilution), Perseus, PP-H6705-00; anti-Klf4 mouse IgG1β (1/100 dilution), Stemgent, 09-0021; anti-Oct4 rabbit polyclonal (1/300 dilution), Santa Cruz, sc-9081; anti-Oct4 mouse IgG2b (1/300 dilution), Santa Cruz, sc-5279; anti-SOX2 mouse IgG2a (1/50), R&D Systems, MAB2018]. After three rinses (10 min each) with TBS, the cells were incubated with fluorochrome-conjugated secondary antibody [Alexa Fluor 488-conjugated donkey anti-rabbit IgG [H + L], (1/500 dilution), Life Technologies, A21206; Alexa Fluor 647-conjugated donkey anti-mouse IgG [H + L], (1/400 dilution), Life Technologies, A31571] at room temperature for 1 h. The cells were examined under confocal imaging (DM 6000 CS SP5; Leica). Acquisitions were performed using an oil immersion objective (40×/1.25 0.75, PL APO HCX; Leica).

## Results

### Cell cycle expression patterns of cyclins and pluripotency factors

Expression patterns of *Cyclin E* (*Ccne1*), *Cyclin E2* (*Ccne2*), and *Cyclin A* (*Ccna1*) were examined during differentiation of mESCs into embryoid bodies, and compared with that of transcription factors implicated in the regulation of pluripotency including *Oct4*, *Sox2*, *Nanog*, *Esrrb*, *Klf2*, *Klf4*, *Klf5*, and *Tfcp2l1*. Only *Cyclin E* mRNA decreased during differentiation concomitantly with mRNA of all the transcription factors analyzed (Figure 1(a)). We next examined the expression of pattern for *Cyclin E*, *Cyclin E2*, and *Cyclin A*, and for 10 pluripotency genes (*Oct4*, *Sox2*, *Nanog*, *Esrrb*, *Klf2*, *Klf4*,
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 1.** Transcript levels for G₁ cyclins and pluripotency factors during cell cycle progression. (a) Gene expression levels measured by qRT-PCR in E14Tg2a mESCs, before and after differentiation to embryoid bodies, after normalization to β-actin (*Actb*) and levels measured on day 0. (b) E14Tg2a–Fucci mESCs, expressing mKO2–hCdt1 and mAG–hGeminin, before and after withdrawal of LIF for 5 days. Top panel: representative fluorescence image. Mid panel: flow cytometry analysis of E14Tg2a–Fucci mESCs showing the distribution of mKO2 (−) mAG(−), mKO2(+) mAG(−), mKO2(−) mAG(+) and mKO2(−) mAG(++) cells [19]. Lower panel: cell population histogram of E14Tg2a-Fucci mESCs showing the DNA content of mKO2(−) mAG(−), mKO2(+) mAG(−), mKO2(−) mAG(+) and mKO2(−) mAG(++) cells after propidium iodide staining (Scale bar = 20 µm). (c) Gene expression levels measured by qRT-PCR after fluorescence-activated cell sorting (FACS) of E14Tg2a–Fucci mESCs in three distinct fractions corresponding to cells in the G₁ [mKO2(−) mAG(−)], S [mKO2(−) mAG(+)], and G₂ [mKO2(−) mAG(++)] phases, respectively. Transcript levels are normalized to β-actin (*Actb*) and the level measured in the total population prior to FACS. (d) Gene expression levels measured by qRT-PCR after FACS of LIF-deprived E14Tg2a–Fucci cells in four distinct fractions corresponding to cells in the early G₁ [mKO2(−) mAG(−)], late G₁ [mKO2(+) mAG(−)], S [mKO2(−) mAG(+)], and G₂ [mKO2(−) mAG(++)] phases, respectively. Expression levels are normalized to β-actin (*Actb*) and the level measured in the total population prior to FACS. (e) Gene expression levels measured by qRT-PCR in E14Tg2a mESCs and EpiSC, normalized to β-actin (*Actb*) and expression measured in E14Tg2a mESCs. (a, c–d) Means and standard deviations (SD) calculated from three independent experiments are shown. (c–d) Statistical analysis was performed using Bonferroni's multiple comparison test (*p* < 0.05, **p** < 0.01, ***p** < 0.001, ****p** < 0.0001).</figcaption></figure>
<figure><img src="image_11.png" /></figure>

Gbx2, Tfcp2l1, Stat3 and Nr5a2) in each phase of the cell cycle in mESCs using the Fucci reporter to sort cells out according to their position in the cell cycle (Figure 1(b)). No significant variation in mRNA levels for Cyclin E and for the 10 pluripotency genes was observed between mAG(-)/mKO2 (-) (G₁ phase), mAG(+)/mKO2(-) (S phase), and mAG(++)/mKO2(-) (G2 phase) (Figure 1(c)). Only Cyclin E2 showed lower expression in the G2 phase. Note that the mAG(-)/mKO2(+) fraction (cells in the late G₁ phase) were excluded from analysis as most of them displayed low or no expression of Oct4 (Suppl. Figure 1A), and therefore have spontaneously committed to differentiation. This could explain why Klf2 expression was higher in the three fractions (G1, S and G2 phases) than in total cells. Co-expression of cyclin E and the pluripotency regulators Oct4, Esrrb and Klf4 was confirmed by immunofluorescence analysis. Rare Oct4-negative, Esrrb-negative, Sox2-negative, and Klf4-negative cells showed low cyclin E content in line with downregulation of Cyclin E transcripts observed during controlled differentiation (Suppl. Figure 1B). After differentiation induced by withdrawal of LIF for 5 days, cyclin E transcripts displayed a somatic-like pattern of cell cycle expression, showing higher levels in G1 and S phase with respect to the G2 phase (Figure 1(d)). Examination of transcripts levels in a whole cell population of EpiSC revealed both a strong reduction for all naïve pluripotency markers and a 70% reduction for Cyclin E (Figure 1(e)). Taken together these results indicate that naïve ES cells express Cyclin E mRNAs at a constant level throughout the cell-cycle. Cyclin E expression is downregulated in pluripotent stem cells in the primed state of pluripotency. It is further reduced in differentiated cells, where it resumes cyclic expression.

## Mapping of cis-regulatory elements in the 5' flanking region of Cyclin E

We examined the transcriptional regulation of Cyclin E by transcription factors of the naïve pluripotency network. The distribution of DNAse I hypersensitivity sites was analyzed over a 15 kb region encompassing the Cyclin E transcription start site in two mESC lines, CJ7 and E14Tg2a, using data available in the mouse ENCODE data

base. We identified three regions hypersensitive to DNAse I in the 5' flanking region: region I between the transcription start site and -1.2 kb, region II between -4.5 and -6 kb, and region III between -9.8 and -11.5 kb (Figure 2(a)). In both Bruce4 and E14Tg2a mESCs, region I and II displayed a strong and moderate enrichment in H3K4me3 histone marks, respectively. The transcription-promoting activity of the three regions was analyzed using a luciferase assay after transfection into E14Tg2a mESCs (Figure 2(b)). Region I showed high transcriptional activity, region II enhanced this activity by a factor of two, and region III reduced it by a factor of two. Regions I, II, and III were thereafter called “promoter” (P), “Proximal Enhancer” (PE), and “Distal regulator Element” (DE), respectively. After transfection into a mouse fibroblast (STO) cell line, luciferase activity was dramatically reduced compared with that observed in mESCs, in line with the reduced transcript levels observed after mESC differentiation. Moreover, in contrast to the situation observed in mESCs, the PE region had no enhancer activity on transcription initiated from the P region.

Using the JASPAR database [39], putative binding sites for the pluripotency regulators Oct4, Sox2, Esrrb, Tfcp2l1, Klf4, Klf5, Tcf3, and STAT3 were identified in the P and PE regions (Suppl. Figure 2). Binding sites for E2F1 and Nr5a2 were also identified in the P region as previously reported [3,4,51]. Binding of transcription factors to their respective sites was analyzed from published ChIP-seq data [36,38]. We identified a strong enrichment in P region-specific sequences after chromatin immunoprecipitation with Esrrb, Tfcp2l1, Klf4, and E2F antibodies, but not with Oct4, Sox2, Nanog, Stat3, Nr5a2 and cMyc antibodies (Figure 2(c)). Based on these results, all subsequent analyses focused on the role of Esrrb, Tfcp2l1, and Klf4 in the transcriptional regulation of Cyclin E via the P region.

## Esrrb, Tfcp2l1, and Klf4 binding sites in the promoter region of Cyclin E

The P region contains two binding sites for Klf4 at positions -9/-19 and -185/-195 with respect to transcription start site (relative scores of 98.4% and 95.7%,
<figure><img src="image_12.png" /><figcaption>a</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_15.png" /></figure>

**Figure 2.** In silico analysis and transcriptional activity of *Cyclin E* promoter region. (a) Mapping of DNAse I hypersensitive sites, H3K4me3 and H3K27me3 binding in the 5' flanking and coding regions of *Cyclin E* identified from ENCODE database [GSM1003830 (DNAseDgf on mESC-CJ7), GSM1014154 (DNAseHS on mESC-E14), and GSM1014187 (DNAseHS on mESC-CJ7), GSM769008 (H3K4me3 on mESC-Bruce4), GSM1000089 (H3K27me3 on mESC-Bruce4), and GSM1000124 (H3K4me3 on mESC-E14) datasets]. (b) Luciferase assay to measure the transcriptional activity of regions P (I), PE (II), and DE (III) in E14Tg2a cells (top panel) and STO cells (bottom panel). Mean and standard deviations (SD) calculated from three independent experiments are shown and two-way Welch test analysis was used to assess significance (ns p ≥ 0.05, \*p < 0.05, \*\*p < 0.01, \*\*\*p < 0.001, \*\*\*\*p < 0.0001). (c) Binding of transcription factors to P (I), PE (II), and DE (III) regions of the *Cyclin E* 5' flanking sequence identified from published previously ChIP-seq data [36,37] [GSM288349 (E2f1), GSM470523 (Nr5a2), GSM288355 (Esrrb), GSM288350 (Tfcp2l1), GSM288354 (Klf4), GSM288345 (Nanog), GSM288346 (Oct4), GSM288347 (Sox2), GSM288353 (Stat3), GSM288356 (c-Myc), and GSM1208217 (Klf4) datasets].

respectively), one binding site for Tfcp2l1 at position -390/-404 (relative score of 90.2%), and two binding sites for Esrrb at positions -538/-548 and -870/-880 (relative scores of 97.4% and 96.3%, respectively) (Figure 3(a)). The role of the five binding sites in *Cyclin E* transcription was explored by site-directed

mutagenesis and analysis of P region transcriptional activity in a luciferase assay (Figure 3(b)). Mutation of Klf4 binding sites had no significant effect on transcriptional activity. Mutation of the Tfcp2l1 and the proximal Esrrb binding sites had a moderate effect on transcription of luciferase (reduction of
<figure><img src="image_16.png" /><figcaption>a</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_18.png" /></figure>

Luciferase assay

<figure><img src="image_19.png" /><figcaption>c</figcaption></figure>

Luciferase assay

<figure><img src="image_20.png" /><figcaption>d</figcaption></figure>

ChIP-seq

<figure><img src="image_21.png" /><figcaption>e</figcaption></figure>

STO

ChIP-qPCR Esrrb

<figure><img src="image_22.png" /></figure>

INPUT IP

<figure><img src="image_23.png" /></figure>

ChIP-qPCR Tfcp2l1

<figure><img src="image_24.png" /></figure>

INPUT IP

<figure><img src="image_25.png" /></figure>

ChIP-qPCR Klf4

<figure><img src="image_26.png" /></figure>

INPUT IP

<figure><img src="image_27.png" /></figure>

**Figure 3.** In silico analysis of the P region of the *Cyclin E* gene. (a) Mapping of Esrrb, Klf4, Tfcp2l1, E2F, and Nr5a2 putative transcription factor binding sites as determined by Transcription Factor Binding Sites annotation (TFBS) public databases as well as published data. (b) Luciferase assay to measure the transcriptional activity of the P region after disruptive mutations of the putative Esrrb, Klf4, and Tfcp2l1 binding sites in E14Tg2a cells (top panel) and STO cells (bottom panel). (c) Luciferase assay to measure the transcriptional activity of the P region of Cyclin E after transfection of pGL4.10-P in E14Tg2a cells, EpiSC and STO cells. (d) Binding of Esrrb, Klf4, and Tfcp2l1 to the P region of *Cyclin E* and to the promoter of *Nanog* determined by ChIP-seq (data mining from [36] and [38] [GSM288350 (Tfcp2l1), GSM288355 (Esrrb), and GSM1208217 (Klf4) datasets]). Graph represents the maximum enrichment value for each region of interest. (e) Binding of Esrrb, Klf4, and Tfcp2l1 to the P region of *Cyclin E* and to the promoter of *Nanog* determined by ChIP-qPCR. (b,c) Means and standard deviations (SD) calculated from at least three independent experiments are shown and two-way Welch test analysis was used to assess significance (ns p ≥ 0.05, *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001).

transcriptional activity to 20% and 58% relative to wild type P region, respectively). In contrast, mutation of the distal Esrrb binding site reduced transcriptional activity to 99% of the wild type sequence (Figure 3(b)). None of these mutations significantly altered the

transcriptional activity of the P region when transfected into STO fibroblast cells. We concluded that the distal Esrrb binding site was essential to the transcriptional activity of the promoter element of *Cyclin E* in mESCs, and the proximal Esrrb site and the
<figure><img src="image_28.png" /></figure>

Tfcp2l1 site play ancillary roles. Luciferase activity associated with pGL4.10-P was also strongly reduced in EpiSC as compared to ES cells (Figure 3(c)), suggesting that the *Cyclin E* promoter region is less active in the primed state than in the naïve state of pluripotency.

Binding of Esrrb, Tfcp2l1, and Klf4 to the P region of *Cyclin E* was measured using previously published ChIP-seq data [36]. We observed a high occupancy of Esrrb binding sites at positions -538/-548 and -870/-880 (Figure 3(d); Suppl. Figure 3). For comparison, the Esrrb binding site found at position -89/-97 in the *Nanog* promoter [52] showed only a low occupancy. In contrast, we observed a low occupancy of Tfcp2l1 and Klf4 binding sites at position -390/-404, -9/-19 and -185/-195 in the *Cyclin E* promoter, compared with a high occupancy at position -4714/-4890, +43/-137 and -4242/4442 in the *Nanog* promoter [34,53]. These data indicate a strong interaction of Esrrb and a much weaker interaction of Tfcp2l1 and Klf4 to their respective predicted binding sites in the *Cyclin E* promoter. Importantly, these results could be corroborated by ChIP-qPCR (Figure 3(e)). Strong enrichment in PCR fragments encompassing the Esrrb binding sites in the *Cyclin E* promoter were observed compared with moderate enrichment in fragments encompassing the Klf4 binding sites and weak enrichment in a fragment encompassing the Tfcp2l1 binding site. As the two Esrrb binding sites are located in close proximity, we could not assess whether Esrrb shows a similar affinity to the proximal as well as the distal predicted binding site. Altogether, these results indicate a significant regulatory role of Esrrb to the *Cyclin E* promoter, in line with the results of the transcriptional activity of mutant promoters.

## Regulation of Cyclin E by Esrrb, Tfcp2l1, and Klf4

To further substantiate the implication of Esrrb, Tfcp2l1 and Klf4 in the transcriptional regulation of *Cyclin E*, their expression was knocked down by means of doxycyclin-induced expression of two independent shRNAs (Sh#1 and Sh#2) in KH2 mESCs [48]. The expression levels of *Esrrb*, *Klf4*, and *Tfcp2l1* could be reduced to less than 10% after 48 h induction with doxycycline for sh#1 (Figure 4(a)). *Cyclin E* transcript levels showed a correlated decrease compared to control cells after knockdown of *Esrrb* (77%), *Tfcp2l1* (69%) and *Klf4* (74%), respectively. In

mESCs expressing Sh#2, doxycycline treatment resulted in a substantially smaller reduction of *Esrrb*, *Tfcp2l1*, and *Klf4* transcript levels than Sh#1 (i.e. to 25%, 63%, and 30% of their original levels, respectively). Nevertheless, the *Cyclin E* transcript level was significantly reduced (65%, 61%, and 24%, respectively). No alteration of *Esrrb*, *Tfcp2l1*, *Klf4*, or *Cyclin E* transcript levels was observed in control cells expressing Sh-Control. In accordance with the reduced transcripts, Cyclin E protein levels decreased to 25% and 50% of the level measured in control cells after knockdown of *Esrrb* and *Tfcp2l1*, respectively (Figure 4(b)). mESCs expressing Sh-Esrrb and Sh-Tfcp2l1 showed only minor alterations in the expression of other pluripotency markers including *Nanog*, *Oct4*, *Sox2*, *Klf2*, *Klf5*, *Tbx3*, *Nr0b1*, and *Zfp42*, indicating that the observed reduction of *Cyclin E* transcript levels after *Esrrb* and *Tfcp2l1* knockdown is not a consequence of differentiation (Suppl. Figure 4). In contrast, *Klf4* knockdown resulted in a significant attenuation of most of the aforementioned pluripotency regulators. As the Klf4 binding sites are apparently not involved in the transcriptional regulation of *Cyclin E*, the observed downregulation of *Cyclin E* in *shKlf4* mESCs might be explained by a substantial rate of spontaneous differentiation. In the next step, Esrrb, Klf4, and Tfcp2l1 were overexpressed using a doxycycline-inducible vector system in KH2 mESCs (Figure 4(c)). This resulted in a 2.7- and 2.2-fold increase in *Cyclin E* transcript levels after Esrrb and Klf4 overexpression, respectively. No significant increase was observed after Tfcp2l1 overexpression. Overall, these results were corroborated by results obtained at the protein level. Taking expected variations between biological replicates into account, cyclin E1 protein level increased after induction of Esrrb, Klf4 and Tfcp2l1 with 0.1 µg/mL doxycycline (Figure 4(d)). Induction with 1 µg/mL doxycycline did not further increase *Ccne1* mRNA and cyclin E1 protein levels. This suggests that *Ccne1* expression has reached its maximum level with only 0.1 µg/mL doxycycline and higher expression of Esrrb and Tfcp2l1 failed to further activate the *Ccne1* promoter. For Klf4 overexpression, we observed maximal activation of Klf4 with 0.1 µg/mL doxycycline at both RNA and protein levels, consistent with our previous observations indicating that strong overexpression of Klf4 is detrimental to self-renewal of mESCs [38]. For Tfcp2l1 overexpression, Cyclin E increase was
<figure><img src="image_29.png" /><figcaption>a</figcaption></figure>

<figure><img src="image_30.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_32.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_33.png" /></figure>

**Figure 4.** Regulation of *Cyclin E* expression by *Esrrb*, *Klf4*, and *Tfcp2l1*. (a) Doxycycline-induced expression of *sh-Esrrb#1*, *sh-Esrrb#2*, *sh-Klf4#1*, *sh-Klf4#2*, *sh-Tfcp2l1#1*, *sh-Tfcp2l1#2*, and *sh-Control* in KH2 mESCs. Expression of the indicated genes is measured by qRT-PCR after 48 h of treatment with doxycycline at 0.1 and 1.0 µg/mL. All gene expression values are normalized to the value measured in the absence of doxycycline (0 µg/mL). (b) Doxycycline-induced expression of *sh-Esrrb#1*, *sh-Klf4#1*, and *sh-Tfcp2l1#1* in KH2 mESCs. Expression of the indicated genes was detected by immunoblotting after 48 h of treatment with doxycycline at 0.1 and 1.0 µg/mL and compared with the absence of doxycycline (0 µg/mL). Histograms show means of 2 independent experiments obtained with 2 different shRNAs. (c,d) Doxycycline-induced expression of *Esrrb*, *Klf4*, and *Tfcp2l1* mouse cDNAs in KH2 mESCs. Expression of the indicated genes was measured by qRT-PCR (c) and immunoblotting (d) after 48 h of treatment with doxycycline at 0.1 and 1.0 µg/mL, and compared with the absence of doxycycline (0 µg/mL). (A,C,D) Means and standard deviations (SD) calculated from three independent experiments are shown and two-way Welch test analysis was used to assess significance (*p* < 0.05, **p** < 0.01, ***p** < 0.001, ****p < 0.0001).
<figure><img src="image_34.png" /></figure>

observed only at the protein level, which may suggest a post-transcriptional regulation as well. Together, these results indicate a correlation between the expression levels of Esrrb and Cyclin E.

## Partial rescue of Esrrb knockdown-induced differentiation by Cyclin E

Several studies have pointed to Esrrb as an inducer of somatic cell reprogramming and mESC self-renewal [29,30,54]. The role of Esrrb in the regulation of Cyclin E therefore prompted us to investigate the capacity of cyclin E to oppose Esrrb knockdown-induced differentiation. We therefore used EKOiE cells expressing doxycycline-regulated Esrrb cDNA in an Esrrb-null background [49]. After doxycycline deprivation for 48 h, Esrrb was undetectable. Consequently, a twofold reduction of both Cyclin E mRNA and cyclin E protein was observed (Figure 5(a)). Expression of Tfcp2l1 and Klf4 were also reduced, which may have contributed to the downregulation of Cyclin E (Figure 5(b)). When doxycycline-deprived EKOiE cells were supplemented with doxycycline for 48 h, they restored expression of Esrrb, Tfcp2l1, Klf4, and Cyclin E to their original levels. Expression of Oct4 and Sox2 remained unchanged in this experimental setting. Next, EKOiE cells were transfected with a plasmid carrying a Myc-tagged rat cyclin E cDNA, or with an empty plasmid [55] (Figure 5(c)). EKOiE-Myccyclin E and EKOiE-control cell populations were analyzed using a colony-forming assay to assess the balance between self-renewal and differentiation (Figure 5(d)). Withdrawal of doxycycline resulted in a gradual increase of the proportion of mixed and differentiated colonies indicating that the loss of Esrrb disrupted self-renewal as previously reported [30]. In the presence of doxycyline, no difference was observed between EKOiE-Myccyclin E and control cells regarding the proportion of undifferentiated, mixed, and differentiated colonies, suggesting that enforced expression of cyclin E has no observable effect on self-renewal in the presence of Esrrb. In contrast, the withdrawal of doxycycline for 3, 5, and 7 days, leads to a significantly increased proportion of mixed and undifferentiated colonies in EKOiE-Myccyclin E cells when compared to control cells. These results strongly suggest that

enforced expression of cyclin E opposes mESC differentiation induced by downregulation of Esrrb.

## Regulation of Cyclin E transcript level by mir-15a

MicroRNA miR-15a was shown to act as a negative regulator of Cyclin E in somatic cells [7,8]. We therefore asked if the decrease in Cyclin E transcript levels observed during ES cell differentiation could also be a result of a rise in miR-15a levels. A qRT-PCR analysis revealed that changes in miR-15a expression mirrored those of Cyclin E between 1 and 9 days of differentiation, suggesting a cross-regulation (Figure 6(a)). We analyzed the miR-15a level in EpiSC and observed a 1.75-fold increase compared to mESCs (Figure 6(b)). This result is consistent with the 70% reduction of Cyclin E transcripts observed in EpiSC versus mESCs (Figure 1(e)). To demonstrate the regulation of Cyclin E expression levels by miR-15a, E14Tg2a mESCs were transfected with a miR-15a-5p mimic, resulting in a 58% reduction of Cyclin E transcripts levels after 48 hours of culture. Nanog transcript levels were also decreased, albeit to a lesser extend (41%) (Figure 6(c)). Together, these results suggest that the high level of cyclin E observed in mouse ESCs could at least in part be attributed to low miR-15a expression.

## Discussion

In somatic cells, expression of cyclin E varies during the cell cycle, reaching a maximum at the G₁/S phase boundary. Repression of the cyclin E gene during G₂/M and the early G₁ phase of the cell cycle are mediated through the assembly of a multiprotein complex containing hypo-phosphorylated Rb, a histone deacetylase as well as the SWI/SNF chromatin remodeling complex, which bind to the cyclin E promoter in order to silence transcription. Transcriptional activation of the cyclin E gene during progression through the G₁ phase depends on the activity of cyclin D/Cdk4 and cyclin D/Cdk6 complexes, which phosphorylate and inactivate retinoblastoma protein (pRb) leading to the release of the repressor proteins of the E2F-transcription factor family [56]. In mESCs, cell cycle regulation of Cyclin E seems to obey a different rule. Despite very low
<figure><img src="image_35.png" /><figcaption>a</figcaption></figure>

<figure><img src="image_36.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_37.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_38.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_39.png" /></figure>

**Figure 5.** Rescue of Esrrb knockdown-induced differentiation by Cyclin E. (a) Expression of Esrrb and Cyclin E in EKOiE cells, before and after withdrawal of doxycycline for 48 h, and analyzed by immunoblotting. Histograms represent the quantification of Esrrb and Cyclin E after normalization to β-actin and to expression in + DOX condition. (b) Expression of *Esrrb*, *Oct4*, *Sox2*, *Nanog*, *Cyclin E* (*Ccne1*), *Klf4*, and *Tfcp2l1* in EKOiE cells (+ LIF + DOX), after withdrawal of doxycycline for 48 h [+ LIF -DOX (48 h)], and after addition of doxycycline for 48 h [+ LIF -DOX (48 h) + DOX (48 h)], normalized to expression measured in control cells [+ LIF + DOX]. Three independent experiments are shown and two-way Welch test analysis was used to assess significance (*p* < 0.05; **p** < 0.01. ***p* < 0.001). (c) qRT-PCR analysis of *Cyclin E* mRNA levels in EKOiE-Myccyclin E and control EKOiE cells using myc-tagged *Cyclin E*-specific primers. (d) Histogram for the proportion of undifferentiated (U), mixed (M), and differentiated (D) colonies in EKOiE cells expressing Myccyclin E and control vector cells before (+ DOX) and after withdrawal of doxycycline (- DOX) for 3, 5 and 7 days. Three independent experiments are shown and two-way Welch test analysis was used to assess significance (**p* < 0.01). Right panel: representative pictures of undifferentiated (U), mixed (M), and differentiated (D) colonies after staining for stem cell-specific alkaline phosphatase activity.

levels of D-type cyclins [12,17,57], the G₁-specific
hypo-phosphorylated form of Rb is almost absent
[17,58], and essentially all E2F transcription factors
are free from Rb proteins and bind the cyclin
E promoter to stimulate transcription, regardless of
the cell's growth cycle position [18]. The theory of
a continuous and uniform expression of *Cyclin E* was
challenged in a recent study using mESCs synchro-
nized with nocodazole, showing that *Cyclin E*
transcripts indeed reached a maximum during G₁
phase [59]. Strikingly, a similar cell cycle-dependent
pattern was observed for Nanog and Esrrb. In sharp
contrast, we did not observe a similar pattern using
non-synchronized mESC-Fucci, which raises the
question of whether chemically-synchronized mESC

can regulate gene expression after release from the
mitotic block.

Esrrb, Tfcp2l1, and Klf4 are three transcription
factors implicated in the control of naive pluripotency.
Esrrb is a direct target gene of both Nanog and the
GSK3/β-catenin/Tcf3 pathway [29,30], and both
Tfcp2l1 and Klf4 are direct target genes of the LIF/
STAT3 signaling pathway [23,25,33]. In the present
study, we showed that *Cyclin E* is a direct target of
these three transcription factors, pointing to regula-
tion of cyclin E expression by the gene regulatory
circuitry that controls naive pluripotency. Among
the three factors studied, Esrrb seems to play a major
role as shown, first, by the dramatic reduction of
transcriptional activity from the promoter region
<figure><img src="image_40.png" /></figure>

<figure><img src="image_41.png" /></figure>

<figure><img src="image_42.png" /></figure>

<figure><img src="image_43.png" /><figcaption>**Figure 6.** Regulation of *Cyclin E* transcript level by miR-15a. (a) RNA expression levels of *Cyclin E* and miR-15a measured by qRT-PCR in E14Tg2a before and after differentiation to embryoid bodies (1–9 days), after normalization to β-actin (*Actb*) and levels measured on day 0. (b) Expression levels of miR-15a measured by qRT-PCR in E14Tg2a and EpiSC after normalization to β-actin (*Actb*). (c) *Cyclin E* (*Ccne1*, *Nanog*, *Oct4* and *Sox2* transcript levels measured by qRT-PCR 48 hours after transfection of 25 mM miR-15a-5p mimic in E14Tg2a mESCs, normalized to β-actin (*Actb*) and levels measured on untransfected cells. Transfection of a miR-1 mimic and qRT-PCR analysis of *Twf1* transcripts is used as control. (a–c) Three independent experiments are shown and two-way Welch test analysis was used to assess significance (ns p ≥ 0.05, \*p < 0.05, \*\*p < 0.01, \*\*\*p < 0.001).</figcaption></figure>

lacking the distal Esrrb binding site, and second, by the downregulation of *Cyclin E* after Esrrb expression has been knocked down or turned off. Moreover, Esrrb overexpression resulted in a 4-fold increase of the steady-state level of *Cyclin E* RNA, further supporting the link between Esrrb and the regulation of *Cyclin E* expression. Tfcp2l1 contributes to the transcriptional regulation of *Cyclin E* in conjunction with Esrrb. The proximal Esrrb binding site overlaps with the binding site for Nr5a2 (LRH-1), a transcription factor involved in the control of naïve pluripotency [60]. In addition, Nr5a2 has been shown to regulate expression of *Cyclin E* in conjunction with β-catenin in intestinal crypt cells [51]. In mESCs, mutation of the Esrrb/Nr5a2 binding site has only a minor effect on the transcriptional activity of the promoter region and ChIP-seq studies revealed no binding of Nr5a2 to the promoter region of *Cyclin E*, strongly suggesting that Nr5a2 plays no role in *Cyclin E* transcriptional regulation in mESCs.

We showed that microRNA miR-15a is a negative regulator of Cyclin E in mESCs, in line with its function in somatic cells [7,8]. Thus,

we propose a regulatory model of *Cyclin E* expression, in which the elevated level of *Cyclin E* transcripts observed in mESCs results from both a transcriptional activation by Esrrb and a lack of negative regulation by miR-15a. Differentiation would trigger both the down-regulation of Esrrb and the elevation of miR-15a, resulting in a rapid drop in *Cyclin E* transcript level. Interestingly, miR-15a is a key direct transcriptional target of E2F in somatic cells [8]. Low expression of miR-15a in mESCs and its up-regulation during differentiation strongly suggests that E2F activity is very low in mESCs, and it is only restored after exit from naïve pluripotency.

We showed that mouse EpiSCs, which epitomize the primed state of pluripotency, express *Cyclin E* at a much lower level as compared to the naïve mESCs. We know little about cell cycle regulation in mouse EpiSCs. However, human ESCs, i.e. the human counterparts of mouse EpiSCs, seem to exhibit a somatic-like cell-cycle regulation. In particular, it was observed that *Cyclin E* mRNA levels increase sharply at the G1/S transition and that the regulation of *Cyclin E* mRNA
<figure><img src="image_44.png" /></figure>

expression levels involves the activation of MEK/ERK pathway and the transcription factors c-Myc and E2F [61]. These observations suggest that the transition from naïve- to primed-state pluripotency is accompanied by the loss of the regulation by transcription factors of the naïve pluripotency network and the gain of the regulation by c-Myc and E2F.

## Acknowledgments

This work was supported by the Infrastructure Nationale en Biologie et Santé INGESTEM (ANR-11-0009), the IHU-B CESAME (ANR-10-IBHU-003), the LabEx REVIVE (ANR-10-LABX-0073), the LabEx "DEVweCAN" (ANR-10-LABX-0061), the LabEx "CORTEX" (ANR-11-LABX-0042), the University of Lyon, program 'Investissements d'Avenir' (ANR-11-IDEX-0007), and the Fondation pour la Recherche Médicale [DEQ20170336757] The authors would like to thank ENAGO (www.enago.com) for the English language review.

## Disclosure statement

No potential conflict of interest was reported by the authors.

## Funding

This work was supported by the Agence Nationale de la Recherche [ANR-10-LABX-0061]; Agence Nationale de la Recherche [ANR-11-LABX-0042]; Agence Nationale de la Recherche [ANR-10-IBHU-003]; Agence Nationale de la Recherche [ANR-11-INBS-0009]; Agence Nationale de la Recherche [ANR-11-IDEX-0007]; Fondation pour la Recherche Médicale [DEQ20170336757].

## ORCID

Pierre-Yves Bourillot http://orcid.org/0000-0003-1657-8555

Pierre Savatier http://orcid.org/0000-0002-8944-9853

## References

[1] Tat Siu K, Rosner MR, Minella AC. An integrated view of cyclin E function and regulation. Cell Cycle. 2012;11:57-64.

[2] Bendris N, Lemmers B, Blanchard JM. Cell cycle, cytoskeleton dynamics and beyond: the many functions of cyclins and CDK inhibitors. Cell Cycle. 2015;14:1786-1798.

[3] Ohtani K, DeGregori J, Nevins JR. Regulation of the cyclin E gene by transcription factor E2F1. Proc Natl Acad Sci U S A. 1995;92:12146-12150.

[4] Botz J, Zerfass-Thome K, Spitkovsky D, et al. Cell cycle regulation of the murine cyclin E gene depends on an E2F binding site in the promoter. Mol Cell Biol. 1996;16:3401-3409.

[5] Harbour JW, Dean DC. The Rb/E2F pathway: expanding roles and emerging paradigms. Genes Dev. 2000;14:2393-2409.

[6] Nevins JR. The Rb/E2F pathway and cancer. Hum Mol Genet. 2001;10:699-703.

[7] Liu Q, Fu H, Sun F, et al. miR-16 family induces cell cycle arrest by regulating multiple cell cycle genes. Nucleic Acids Res. 2008;36:5391-5404.

[8] Ofir M, Hacohen D, Ginsberg D. MiR-15 and miR-16 are direct transcriptional targets of E2F1 that limit E2F-induced proliferation by targeting cyclin E. Mol Cancer Res. 2011;9:440-447.

[9] Burdon T, Smith A, Savatier P. Signalling, cell cycle and pluripotency in embryonic stem cells. Trends Cell Biol. 2002;12:432-438.

[10] Savatier P, Malashicheva A Cell-cycle control in embryonic stem cells. Handbook of stem cells. Embryonic Stem Cells. 2004;2: 53-62.

[11] Abdelalim EM. Molecular mechanisms controlling the cell cycle in embryonic stem cells. Stem Cell Rev. 2013;9:764-773.

[12] Savatier P, Lapillonne H, van Grunsven LA, et al. Withdrawal of differentiation inhibitory activity/leukemia inhibitory factor up-regulates D-type cyclins and cyclin-dependent kinase inhibitors in mouse embryonic stem cells. Oncogene. 1996;12:309-322.

[13] Faast R, White J, Cartwright P, et al. Cdk6-cyclin D3 activity in murine ES cells is resistant to inhibition by p16(INK4a). Oncogene. 2004;23:491-502.

[14] Aladjem MI, Spike BT, Rodewald LW, et al. ES cells do not activate p53-dependent stress responses and undergo p53-independent apoptosis in response to DNA damage. Curr Biol. 1998;8:145-155.

[15] Hong Y, Stambrook PJ. Restoration of an absent G1 arrest and protection from apoptosis in embryonic stem cells after ionizing radiation. Proc Natl Acad Sci U S A. 2004;101:14443-14448.

[16] van der Laan S, Tsanov N, Crozet C, et al. High dub3 expression in mouse ESCs couples the G1/S checkpoint to pluripotency. Mol Cell. 2013;52:366-379.

[17] Stead E, White J, Faast R, et al. Pluripotent cell division cycles are driven by ectopic Cdk2, cyclin A/E and E2F activities. Oncogene. 2002;21:8320-8333.

[18] White J, Stead E, Faast R, et al. Developmental activation of the Rb-E2F pathway and establishment of cell cycle regulated Cdk activity during embryonic stem cell differentiation. Mol Biol Cell. 2005;16:2018-2027.

[19] Coronado D, Godet M, Bourillot PY, et al. A short G1 phase is an intrinsic determinant of naive embryonic stem cell pluripotency. Stem Cell Res. 2013;10:118-131.

[20] Liu L, Michowski W, Inuzuka H, et al. G1 cyclins link proliferation, pluripotency and differentiation of embryonic stem cells. Nat Cell Biol. 2017;19:177-188.
<figure><img src="image_45.png" /></figure>

[21] Soufi A, Dalton S. Cycling through developmental decisions: how cell cycle dynamics control pluripo-tency, differentiation and reprogramming. Development. 2016;143:4301-4311.

[22] Qiu D, Ye S, Ruiz B, et al. Klf2 and Tfcp2l1, two Wnt/beta-catenin targets, act synergistically to induce and maintain naive pluripotency. Stem Cell Reports. 2015;5:314-322.

[23] Hall J, Guo G, Wray J, et al. Oct4 and LIF/Stat3 additively induce Kruppel factors to sustain embryonic stem cell self-renewal. Cell Stem Cell. 2009;5:597-609.

[24] Yeo JC, Jiang J, Tan ZY, et al. Klf2 is an essential factor that sustains ground state pluripotency. Cell Stem Cell. 2014;14:864-872.

[25] Bourillot PY, Aksoy I, Schreiber V, et al. Novel STAT3 target genes exert distinct roles in the inhibition of mesoderm and endoderm differentiation in coopera-tion with Nanog. Stem Cells. 2009;27:1760-1771.

[26] Bourillot PY, Savatier P. Kruppel-like transcription factors and control of pluripotency. BMC Biol. 2010;8:125.

[27] Ema M, Mori D, Niwa H, et al. Kruppel-like factor 5 is essential for blastocyst development and the normal self-renewal of mouse ESCs. Cell Stem Cell. 2008;3:555-567.

[28] Parisi S, Passaro F, Aloia L, et al. Klf5 is involved in self-renewal of mouse embryonic stem cells. J Cell Sci. 2008;121:2629-2634.

[29] Martello G, Sugimoto T, Diamanti E, et al. Esrrb is a pivotal target of the gsk3/tcf3 axis regulating embryonic stem cell self-renewal. Cell Stem Cell. 2012;11:491-504.

[30] Festuccia N, Osorno R, Halbritter F, et al. Esrrb is a direct nanog target gene that can substitute for nanog function in pluripotent cells. Cell Stem Cell. 2012;11:477-490.

[31] Khalfallah O, Rouleau M, Barbry P, et al. Dax-1 knock-down in mouse embryonic stem cells induces loss of pluripotency and multilineage differentiation. Stem Cells. 2009;27:1529-1537.

[32] Tai CI, Ying QL. Gbx2, a LIF/Stat3 target, promotes reprogramming to and retention of the pluripotent ground state. J Cell Sci. 2013;126:1093-1098.

[33] Martello G, Bertone P, Smith A. Identification of the missing pluripotency mediator downstream of leukaemia inhibitory factor. Embo J. 2013;32:2561-2574.

[34] Ye S, Li P, Tong C, et al. Embryonic stem cell self-renewal pathways converge on the transcription factor Tfcp2l1. Embo J. 2013;32:2548-2560.

[35] Kent WJ, Sugnet CW, Furey TS, et al. The human genome browser at UCSC. Genome Res. 2002;12:996-1006.

[36] Chen X, Xu H, Yuan P, et al. Integration of external signaling pathways with the core transcriptional network in embryonic stem cells. Cell. 2008;133:1106-1117.

[37] Heng JC, Feng B, Han J, et al. The nuclear receptor Nr5a2 can replace Oct4 in the reprogramming of murine somatic cells to pluripotent cells. Cell Stem Cell. 2010;6:167-174.

[38] Aksoy I, Giudice V, Delahaye E, et al. Klf4 and Klf5 differentially inhibit mesoderm and endoderm differentiation in embryonic stem cells. Nat Commun. 2014;5:3719-3733.

[39] Sandelin A, Alkema W, Engstrom P, et al. JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004;32:D91-D4.

[40] Wingender E, Chen X, Fricke E, et al. The TRANSFAC system on gene expression regulation. Nucleic Acids Res. 2001;29:281-283.

[41] Riva A. The MAPPER2 Database: a multi-genome catalog of putative transcription factor binding sites. Nucleic Acids Res. 2012;40:D155-D61.

[42] Robertson G, Bilenky M, Lin K, et al. cisRED: a database system for genome-scale computational discovery of regulatory elements. Nucleic Acids Res. 2006;34:D68-D73.

[43] Newburger DE, Bulyk ML. UniPROBE: an online database of protein binding microarray data on protein-DNA interactions. Nucleic Acids Res. 2009;37:D77-D82.

[44] Fu Y, Frith MC, Haverty PM, et al. MotifViz: an analysis and visualization tool for motif discovery. Nucleic Acids Res. 2004;32:W420-W3.

[45] Sandelin A, Wasserman WW, Lenhard B. ConSite: web-based prediction of regulatory elements using cross-species comparison. Nucleic Acids Res. 2004;32: W249-W52.

[46] Aksoy I, Jauch R, Chen J, et al. Oct4 switches partnering from Sox2 to Sox17 to reinterpret the enhancer code and specify endoderm. Embo J. 2013;32:938-953.

[47] Premsrirut PK, Dow LE, Kim SY, et al. A rapid and scalable system for studying gene function in mice using conditional RNA interference. Cell. 2011;145:145-158.

[48] Beard C, Hochedlinger K, Plath K, et al. Efficient method to generate single-copy transgenic mice by site-specific integration in embryonic stem cells. Genesis. 2006;44:23-28.

[49] Festuccia N, Dubois A, Vandormael-Pournin S, et al. Mitotic binding of Esrrb marks key regulatory regions of the pluripotency network. Nat Cell Biol. 2016;18:1139-1148.

[50] Osteil P, Studdert J, Wilkie E, et al. Generation of genome-edited mouse epiblast stem cells via a detour through ES cell-chimeras. Differentiation. 2016;91:119-125.

[51] Botrugno OA, Fayard E, Annicotte JS, et al. Synergy between LRH-1 and beta-catenin induces G1 cyclin-mediated cell proliferation. Mol Cell. 2004;15:499-509.

[52] van Den Berg DL, Zhang W, Yates A, et al. Estrogen-related receptor beta interacts with Oct4 to positively regulate Nanog gene expression. Mol Cell Biol. 2008;28:5986-5995.
<figure><img src="image_46.png" /></figure>

[53] Zhang P, Andrianakos R, Yang Y, et al. Kruppel-like factor 4 (Klf4) prevents embryonic stem (ES) cell differentiation by regulating Nanog gene expression. *J Biol Chem.* 2010;285:9180-9189.

[54] Feng B, Jiang J, Kraus P, et al. Reprogramming of fibroblasts into induced pluripotent stem cells with orphan nuclear receptor Esrrb. *Nat Cell Biol.* 2009;11:197-203.

[55] Matsumoto Y, Maller JL. A centrosomal localization signal in cyclin E required for Cdk2-independent S phase entry. *Science*. 2004;306:885-888.

[56] Moroy T, Geisen C. Cyclin E. *Int J Biochem Cell Biol.* 2004;36:1424-1439.

[57] Jirmanova L, Afanassieff M, Gobert-Gosse S, et al. Differential contributions of ERK and PI3-kinase to the regulation of cyclin D1 expression and to the control of the G1/S transition in mouse embryonic stem cells. *Oncogene*. 2002;21:5515-5528.

[58] Savatier P, Huang S, Szekely L, et al. Contrasting patterns of retinoblastoma protein expression in mouse embryonic stem cells and embryonic fibroblasts. *Oncogene*. 1994;9:809-818.

[59] van der Laan S, Golfetto E, Vanacker JM, et al. Cell cycle-dependent expression of Dub3, Nanog and the p160 family of nuclear receptor coactivators (NCoAs) in mouse embryonic stem cells. *PLoS ONE*. 2014;9: e93663.

[60] Guo G, Smith A. A genome-wide screen in EpiSCs identifies Nr5a nuclear receptors as potent inducers of ground state pluripotency. *Development*. 2010;137:3185-3192.

[61] Rodriguez Varela MS, Mucci S, Videla Richardson GA, et al. Regulation of cyclin E1 expression in human pluripotent stem cells and derived neural progeny. *Cell Cycle*. 2018;7:1721-1744.