The polycomb group protein Yaf2 regulates the pluripotency of embryonic stem cells in a phosphorylation-dependent manner

The polycomb group (PcG) proteins are key epigenetic regulators in stem cell maintenance. PcG proteins have been thought to act through one of two polycomb repressive complexes (PRCs), but more recent biochemical analyses have challenged this model in the identification of noncanonical PRC1 (nc-PRC1) complexes characterized by the presence of Rybp or Yaf2 in place of the canonical Chromobox proteins. However, the biological significance of these nc-PRC1s and the potential mechanisms by which they mediate gene repression are largely unknown. Here, we explore the functional consequences of Yaf2 disruption on stem cell regulation. We show that deletion of Yaf2 results in compromised proliferation and abnormal differentiation of mouse embryonic stem cells (mESCs). Genome-wide profiling indicates Yaf2 functions primarily as a transcriptional repressor, particularly impacting genes associated with ectoderm cell fate in a manner distinct from Rybp. We confirm that Yaf2 assembles into a noncanonical PRC complex, with deletion analysis identifying the region encompassing amino acid residues 102–150 as required for this assembly. Furthermore, we identified serine 166 as a Yaf2 phosphorylation site, and we demonstrate that mutation of this site to alanine (S166A) compromises Ring1B-mediated H2A monoubiquitination and in turn its ability to repress target gene expression. We therefore propose that Yaf2 and its phosphorylation status serve as dual regulators to maintain the pluripotent state in mESCs.

The polycomb group (PcG) proteins are key epigenetic regulators in stem cell maintenance. PcG proteins have been thought to act through one of two polycomb repressive complexes (PRCs), but more recent biochemical analyses have challenged this model in the identification of noncanonical PRC1 (nc-PRC1) complexes characterized by the presence of Rybp or Yaf2 in place of the canonical Chromobox proteins. However, the biological significance of these nc-PRC1s and the potential mechanisms by which they mediate gene repression are largely unknown. Here, we explore the functional consequences of Yaf2 disruption on stem cell regulation. We show that deletion of Yaf2 results in compromised proliferation and abnormal differentiation of mouse embryonic stem cells (mESCs). Genomewide profiling indicates Yaf2 functions primarily as a transcriptional repressor, particularly impacting genes associated with ectoderm cell fate in a manner distinct from Rybp. We confirm that Yaf2 assembles into a noncanonical PRC complex, with deletion analysis identifying the region encompassing amino acid residues 102-150 as required for this assembly. Furthermore, we identified serine 166 as a Yaf2 phosphorylation site, and we demonstrate that mutation of this site to alanine (S166A) compromises Ring1B-mediated H2A monoubiquitination and in turn its ability to repress target gene expression. We therefore propose that Yaf2 and its phosphorylation status serve as dual regulators to maintain the pluripotent state in mESCs.
How the PcG proteins contribute to gene silencing or activation in diverse physiological states is of interest. The Ring-Pcgf module constitutes a minimal core on which the different PRC1s are assembled and deposits H2AK119ub1, a histone modification that impedes RNA polymerase II elongation (2,15), through the E3 protein ligase subunits Ring1A/B. Indeed, our recent studies suggest that of the Pcgf family members, only Pcgf3/5 are essential for PRC1-mediated deposition of H2AK119ub1 in mESC (16), supporting a conclusion that some Pcgf proteins confer differing E3 ligase activities on PRC1. Furthermore, the E3 ligase activity of the Pcgf4/Ring1B complex has been reported to be modulated by its auxiliary subunits RYBP, PHC2, or CBX2/8 (6). Interestingly, it has been shown that Mel-18 phosphorylation is required for targeting the Mel-18/Ring1B complex to chromatin (17). In contrast, phosphorylation of Ring1B by CK2 impairs its E3 ubiquitin-ligase activity (18). Therefore, post-translational modifications such as phosphorylation add another layer of regulation to E3 ligase activities on PRC1. Although several PcG proteins of the PRC1 complexes have been found to be phosphorylated (17,19,20), the defined contribution of this modification to Polycomb chromatin domain function in vivo has remained largely unknown.
In this study, we provide evidence for a phosphorylationbased mechanism that affects Yaf2-mediated pluripotency programs in mESCs. Yaf2 protein was first identified in a yeast two-hybrid screen for Yin-Yang-1 (YY-1) transcriptional factor (21) and was further shown to modify the transcriptional activity of E4TF1/hGABP and Myc-mediated transactivation and transformation (22,23). Genetic studies in zebrafish have implicated Yaf2 as a survival factor during early development and organogenesis (24). Although Yaf2 has recently been identified as a component of noncanonical PRC1 complexes (6), its defined contribution and relevance to Polycomb function in vivo still await elucidation. Using the CRISPR/Cas9 technique, we found that deletion of Yaf2 in mESCs resulted in defective proliferation and differentiation, revealing critical roles for Yaf2 in the maintenance of pluripotency in mESCs. Furthermore, Yaf2 is phosphorylated at serine 166, and this modification mediates its transcriptional repression and is essential for fully catalyzing monoubiquitination of histone H2A in vivo. In summary, these findings therefore provide new insights into the functional and mechanistic role of Yaf2 in regulation of pluripotency.

Yaf2 has an indispensable function in mESC self-renewal
To investigate potential functions of Yaf2, a characteristic subunit of noncanonical PRC1 complexes in mESCs, several Yaf2-deficient mESC clones (referred to as Yaf2 Ϫ/Ϫ ) were established by utilizing CRISPR/Cas9 technique (Fig. S1) (25). The Yaf2 knockout mESCs were viable and exhibited typical mESC morphological features when grown on a feeder layer of mouse embryonic fibroblasts (MEFs), indistinguishable from those formed by WT counterparts (Fig. 1A). They did, however, display a stable and mild proliferation defect and moderate reduction in colony-forming ability relative to their controls (Fig. 1, B and C). The knockout mESCs accumulated in the G 0/1 phase of the cell cycle, increasing from ϳ14% of WT in G 0/1 to 20% of the Yaf2 knockout cells (Fig. 1D). Remarkably, lentivirus-based re-introduction of Yaf2 back into Yaf2 Ϫ/Ϫ mESCs rescued the population doubling time and G 0/1 phenotypes. The smaller size of Yaf2-deficient mESC colonies was not due to increased apoptosis as we found no significant difference in the percentage of mESCs undergoing apoptosis between the WT and Yaf2 Ϫ/Ϫ when assessed by annexin-V staining (Fig.  1E). Thus, the impaired growth of Yaf2 Ϫ/Ϫ mESCs was likely due to their altered cell cycle profile. Notably, Yaf2 knockout mESCs continued to express high levels of undifferentiated state-specific markers such as Oct4, Nanog, Sox2, and alkaline phosphatase (Fig. 1, F and G). Taken together, these results indicate that Yaf2 is required to retain proper self-renewal capacity of mESCs.

Yaf2 controls differentiation programs in mESCs
To test the potential of Yaf2 Ϫ/Ϫ mESCs to undergo in vitro differentiation, mutant and control mESCs were tested for their capacity to form embryoid bodies (EBs). As shown in Fig. 2A, Yaf2 Ϫ/Ϫ mESCs did form EBs, but they were reduced in size compared with those of the controls during differentiation (Fig.  2B). RT-PCR analysis of these EBs, which were kept in differentiation medium for 3, 7, and 12 days, revealed up-regulation of genes associated with each of the three germ layers (26), including Gata6 and Foxa2 for endoderm; Fgf5, Nestin, and Sox11 for ectoderm; and Brachyury and Flk1 for mesoderm (Fig. 2C). Further analysis indicated that these changes were accompanied by diminished expression of mESC pluripotency markers, including Oct4 and Nanog. Importantly, throughout Yaf2 Ϫ/Ϫ EB differentiation and even in undifferentiated Yaf2 Ϫ/Ϫ mESCs, the expression of Fgf5, Nestin, and Sox11 was markedly increased with respect to control EBs, whereas higher levels of endoderm and mesoderm markers were only observed at days 7 and 12 of Yaf2 Ϫ/Ϫ EBs, indicating that Yaf2 mainly prevents the differentiation of mESCs into ectoderm (Fig. 2C). Strikingly, lentiviral expression of FLAG-tagged Yaf2 was able to rescue EB size and the pattern of expression of these germ layer markers to levels similar to those in control cells (Fig. 2, A-C). To further investigate the differentiation capacity of Yaf2 Ϫ/Ϫ mESCs, WT and mutant mESCs were subcutaneously A, representative phase images of EBs derived from mESCs of the indicated genotypes at different time points following differentiation. Images were taken at ϫ4 magnification at days 3, 7, and 12. B, bar graphs indicate mean diameter of 20 random EBs of the indicated genotypes. C, real-time qPCR was used to assess the expression levels of pluripotency and lineage-specific markers in mRNA samples derived from mESCs and differentiated EBs of the indicated genotypes at days 3, 7, and 12. Bar graphs represent the mean of three independent biological repeats. Error bars indicate Ϯ S.D. Student's t test analysis.

Yaf2 in embryonic stem cells
injected into the flanks of 6-week-old immune-deficient mice to form teratomas (Fig. S2). Immunohistochemistry showed that various cell types and tissues derived from all three germ layers were present in WT mESC-derived teratomas. However, teratomas from Yaf2 Ϫ/Ϫ mESCs were significantly smaller, and we noted a larger proportion of cells with an ectodermal origin, as evidenced in an immunostaining using antibody specific to glial fibrillary acidic protein (ectoderm marker), when compared with WT teratomas. Notably, Yaf2 deletion has no effect on endoderm and mesoderm specification, as demonstrated by TROMA1 and SMA1, respectively, staining (Fig. S2). These results indicate that Yaf2 is necessary for the differentiation or maintenance of the ectodermal lineage.

Genome-wide profiling of gene expression in Yaf2 ؊/؊ mESCs
To identify the downstream target genes of Yaf2, we used RNA-seq to compare the transcriptome profiling of WT and Yaf2 Ϫ/Ϫ mESCs. 351 genes were up-regulated, and 146 genes ( Fig. 3A and Table S3) were down-regulated based on a foldchange cutoff of 2 and a p value less than 0.05, indicating that Yaf2 functions predominantly as a transcriptional repressor in mESCs (Fig. 3, A and B). Gene ontology (GO) analysis indicated that the up-regulated genes in Yaf2 Ϫ/Ϫ are statistically associated with ectoderm differentiation genes (e.g. central nervous system development and epithelial cell differentiation) (Fig.  3C), implying a particular role for Yaf2 in ectoderm cell fate. Consistently, of the top 31 up-regulated genes with Ͼ10-fold increase in expression, 10 have been reported to be essential or are expected to have a function in ectoderm development (Fig.  3D). To determine the potential cross-talk between Yaf2 and its homolog Rybp, Rybp null mESCs were also generated, and their transcriptional status was analyzed by RNA-seq (Fig. S3). There were 904 differentially expressed genes between the WT and Rybp Ϫ/Ϫ mESCs, with 609 up-and 295 down-regulated genes A, volcano plots depicting significant differentially expressed genes in Yaf2 Ϫ/Ϫ versus WT mESCs. Up-regulated genes were labeled in red and down-regulated genes were labeled in green if they had a log2 fold change of Ͼ1 or less than Ϫ1, respectively. The number of differentially expressed genes they represented is indicated. B, RNA-seq heat map of the 497 expressed transcripts with Ͼ2-fold expression differences in WT and Yaf2 Ϫ/Ϫ mESCs. Red indicates high expression, and blue indicates low expression. C, Gene Ontology analysis of biological functions of Yaf2 target genes. D, fold changes in the expression of the top 31 genes, each of which displayed more than 10-fold up-regulation in Yaf2 Ϫ/Ϫ compared with WT. Ectoderm-specific genes are highlighted in red.

Yaf2 in embryonic stem cells
in Rybp Ϫ/Ϫ mESCs (Table S4). Of the 48 up-regulated genes with Ͼ10-fold increase in expression, 16 were strongly associated with germline development (meiosis and spermatogenesis), consistent with a role of Rybp in PRC1.6 activity (26,27). Further analysis revealed no significant overlap between the genes regulated by Yaf2 and Rybp (Fig. S3), indicating Yaf2 and Rybp serve distinct functions in mESCs. Together, these results strongly suggest that the repression of developmental genes by Yaf2 may contribute to the maintenance of mESC identity.

Yaf2 assembles a noncanonical PRC1 in mESCs
To investigate the molecular background of Yaf2 function, we characterized the yet-unknown physical interactome of Yaf2 in mESCs. Yaf2 Ϫ/Ϫ mESCs were stably transfected with a C-terminal FLAG-tagged form of Yaf2, and nuclear extracts were subjected to immunoprecipitation with anti-FLAG antibody, and then the immunoprecipitates were analyzed by mass spectrometric analysis to determine associated proteins. Fig.  4A shows a summary of Yaf2 interactome, where we focused on known canonical and noncanonical PRC1 complexes. Yaf2 copurified many proteins that reside in noncanonical PRC1 complexes (PRC1.1, PRC1.3/5, and PRC1.6). In FLAG-Yaf2 immunoprecipitates, Rybp was not detected, in line with previous reports showing that Yaf2 and Rybp are mutually exclusive in PRC1 complexes. Interestingly, we found that Yaf2 specifically co-purified with Pcgf2 and Pcgf4, but not with canonical components, such as polyhomeotic homolog proteins (Phcs), chromobox proteins (Cbxs), and Sex comb on midleg proteins (Scms), which is consistent with the finding that Cbx proteins and Yaf2 are mutually exclusive within the PRC1 complex context, as both compete for the same surface on the C-terminal domain of Ring1B (28). These interactions were confirmed using specific antibodies with deficient mESCs rescued with Yaf2-FLAG (Fig. 4, C and D). Of note, we did not detect Yaf2 interactions with Rybp by Western blotting, consistent with the MS results.
To define the domains in Yaf2 that mediate its association with other members of the complex, we constructed a series of Additionally, the ability of Yaf2 mutants to assemble into the PRC1 complexes is also shown. 5% of the cell lysates (Input) used in each assay was loaded as input. Rybp, Suz12, and Yy1 are shown as negative controls. Of note, the C-terminal region of 29 residues (151-179) is required for keeping the doublet of Yaf2 protein on SDS-polyacrylamide gels. E, quantitative RT-PCR detecting Fgf5, Csn3, and Gfra1 expression levels in Yaf2 Ϫ/Ϫ mESCs following re-expression of WT and Yaf2 mutants. Error bars indicate Ϯ S.D. Bar graphs represent the mean of three independent biological repeats. Student's t test analysis.

Yaf2 in embryonic stem cells
deletion mutants of Yaf2 and first tested their ability to interact with other subunits in the complex in mESCs by immunoprecipitation followed by immunoblotting. As summarized in Fig.  4B, a deletion mutant of Yaf2 lacking amino acids 102-150 (⌬102-150) failed to interact with other components of the complexes, consistent with its inability to rescue growth defects associated with Yaf2 deficiency in mESCs. The C-terminal deletion mutant (⌬151-179) and the RanBP2-type zinc finger domain deletion (⌬19 -48) and point mutants (4C-A, in which Cys-25, Cys-28, Cys-39, and Cys-42 sites were mutated from cysteine to alanine), however, retained their ability to incorporate into the complexes in vivo, although they failed to restore the proliferation defect observed in Yaf2 Ϫ/Ϫ mESCs (Fig. 4,  B-D). These studies thus define a novel domain in Yaf2 that mediates its interaction with other components of the complexes. Additionally, these deletion mutants were also assessed for their ability to inhibit the expression of Yaf2 target genes in mESCs. Three Yaf2 deletion mutants, including ⌬19 -48, ⌬102-150, and ⌬151-179 (Fig. 4E), failed to restore target gene repression. In contrast, mutants deleted of amino acids 2-18 or 49 -101 (⌬2-18 and ⌬49 -101, respectively) showed a WT level of gene repression.

Critical role for Yaf2 phosphorylation in maintaining mESC identity
The WT Yaf2 protein migrated as a doublet on SDS-PAGE, and deletion of amino acids 151-179, but not other regions, caused the disappearance of the top band of the doublet (Fig. 4,  C and D). Furthermore, calf intestine phosphatase (CIP) treatment resulted in the same change in WT Yaf2 protein (Fig. 5A), indicating that phosphorylation within the region between amino acids 151 and 179 is causal to the motility shift. Strikingly, mutation of three consecutive serines (164 -166) within this region to alanines abolished the slowly migrating band (Fig.  5, B and C). To identify the exact amino acid(s) that is phosphorylated, we substituted these three serine residues with alanine individually to prevent phosphorylation. We found that only the S166A mutant migrated as a single band with mobility Whole-cell lysate (Input) was also subjected to Western blotting to determine the expression of Gsk3␣/␤ and Yaf2. H, quantitative RT-PCR analysis to determine mRNA expression of Yaf2 target genes in mESCs treated with the Gsk3 inhibitor CHIR (6 M) for 24 h. Expression levels of these genes were normalized and compared with those of the untreated control cells. Error bars indicate Ϯ S.D. Bar graphs represent the mean of three independent biological repeats. Student's t test analysis was used.

Yaf2 in embryonic stem cells
identical to the lower band of the doublet, whereas the S166D mutant, where the serine was mutated to aspartic acid to mimic phosphoserine, migrated as a single band with mobility identical to the slower migrating band of WT Yaf2 (Fig. 5C).
To evaluate the role of Yaf2 Ser-166 phosphorylation in regulating mESC self-renewal and differentiation, we stably introduced full-length or phosphorylation mutant Yaf2 into Yaf2 Ϫ/Ϫ mESCs. As shown in Fig. 5D, S166D fully restored WT levels of growth in mESCs. However, phosphodefective S166A abolished its ability to rescue the colony-growth defect associated with the deficiency of the Yaf2 gene in mESCs. We examined the ability of the S166A mESCs to undergo differentiation by method of EB and teratoma formation. These S166A mESCderived EB exhibited higher expression of the ectodermal marker Nestin, Fgf5, and Sox11 throughout EB differentiation (Fig. S4). Consistent with these observations, teratomas derived from S166A mESCs contained more advanced neuron-selected structures (Fig. S2). Furthermore, S166A also failed to suppress Yaf2-mediated target gene repression in undifferentiated mESCs (Fig. 5E). Thus, these results provide strong evidence that the Ser-166 can be phosphorylated and highlight the pivotal role of this residue for maintaining the function of Yaf2 protein in mESCs.
Close examination of the amino acid residues surrounding Ser-166 revealed the presence of a glycogen synthesis kinase 3 (Gsk3) consensus sequence (Ser/Thr-(Xaa-Xaa-Xaa)-Ser/Thr, with Xaa being any amino acid) (29). To gain further support that Gsk3 is responsible for Ser-166 phosphorylation, we used a most selective Gsk3 inhibitor, CHIR99021. To this end, mESCs expressing FLAG-tagged Yaf2 were treated in the presence or absence of CHIR99021. Results shown in Fig. 5F demonstrated that phosphorylation of Yaf2 was greatly inhibited by the presence of CHIR99021 as evidenced by Western blotting. Similar analysis also demonstrated that CHIR99021 treatment resulted in decreased endogenous phospho-Yaf2 levels in WT mESCs (Fig. 5F), indicating that Gsk3 at least partially contributes to Yaf2 phosphorylation. This notion was further supported by the co-immunoprecipitation experiment that detected an interaction between Yaf2 and Gsk3␣, but not Gsk3␤, in mESCs (Fig. 5G). In addition, treatment of mESCs with CHIR99021 also derepressed the activity of Yaf2-mediated target gene expression (Fig. 5H). Taken together, Gsk3 is likely responsible for the phosphorylation of Yaf2 at the Ser-166 site to control developmental gene expression.

Yaf2 phosphorylation status correlates with levels of H2AK119ub1 at specific promoters
To test whether Yaf2 can directly regulate target genes, we examined Yaf2 occupancy on the promoter of a panel of selected target genes, whose expressions were significantly derepressed upon Yaf2 deletion (Fig. 6A), with the use of the ChIP-qPCR assay. ChIP analysis revealed robust binding by Yaf2 to target promoters in WT mESCs, but the signal was dramatically reduced to background levels in Yaf2 Ϫ/Ϫ mESCs (Fig. 6B). As expected, the PRC1.6 (Tdrkh and Tex101) and canonical PRC1 (Emoes and Foxa2) target genes were not significantly enriched and served as a control in our ChIP experiments. As an additional control for specificity, we showed that an irrelevant antibody (IgG) failed to immunoprecipitate Yaf2 target promoter sequences in ChIP experiments (Fig. 6C). These strongly suggest that Yaf2 targets to its target promoters directly in mESCs. Similar ChIP experiments were also performed with Yaf2 Ϫ/Ϫ and WT mESCs to determine whether the absence of Yaf2 protein would have an impact on the recruitment of Ring1B, Rybp (components of PRC1), or Suz12 (PRC2 complex). Unexpectedly, ChIP data indicate that recruitment of these components was minimally affected by loss of Yaf2 (Fig. 6C and Fig. S5), which is inconsistent with the results obtained in HeLa cells (30).
Importantly, despite our observation that promoter occupancy by H3K27me3 was not affected by loss of Yaf2, we observed that H2AK119ub1 enrichment was reduced by 2-3fold in the Yaf2 Ϫ/Ϫ mESCs at these specific targets of Yaf2 ( Fig.  6C and Fig. S5). Furthermore, deletion of Yaf2 also resulted in moderate reduction in global H2AK119ub1 as measured by Western blotting (Fig. S5). To verify that the reduction of H2AK119ub1 levels at genes bound by Yaf2 in knockout cells is only due to the lack of Yaf2, we reexpressed both WT and mutant Yaf2 in the knockout mESCs using lentiviral expression vectors. WT Yaf2 displayed full rescue of H2AK119ub1 levels as well as target gene derepression. However, a Yaf2 phosphorylation-defective mutant S166A virtually failed to restore H2AK119ub1 levels at target promoters or target gene repression in null mESCs, although it displayed a binding affinity for target promoters similar to that of the WT (Fig. 6, A-C). This indicated that phosphorylation of Yaf2 at Ser-166 enhances Ring1B-mediated deposition of H2AK119ub1, thus contributing to target gene repression.

Discussion
The central conclusion from this study is the demonstration that Yaf2 contributes to the maintenance of mESC pluripotency in a manner that depends upon its phosphorylation state. Interestingly, we find that phosphorylation of Yaf2 is essential for Ring1B-mediated E3 ubiquitin ligase activity. Thus, this study suggests a mechanism for how signaling affects PcG-mediated chromatin-based epigenetic processes to ensure proper regulation of cell fate at the molecular level and how the enzymatic activity of PRC1 is regulated.
Recent studies show that the composition of PRC1 complexes is considerably far more diverse than previously recognized. In mammals, six major groups of PRC1 complexes, referred to as PRC1.1-1.6 based on available Pcgf subunits, and their potential roles in H2AK119ub1 have been characterized (6). Although distinct types of PRC1 are believed to catalyze H2AK119ub1 due to the presence of a common subunit Ring1B, they do not appear to have an equal effect on enzymatic activity. Indeed, we have shown that among Pcgf family members only, Pcgf3/5 are required for maintaining global H2AK119ub1 levels in mESCs (16). In vitro ubiquitylation assays demonstrated that Ring1A/B with either Pcgf2 or Pcgf4 can catalyze H2AK119ub1 (17,31), but incorporation of the auxiliary subunits Rybp, CBX2/8, or PHC2 into these minimal complexes has been reported to affect their E3 ubiquitin ligase activity (6). According to our findings, PRC1-mediated H2AK119ub1 in mESCs is significantly reduced after deletion Yaf2 in embryonic stem cells of Yaf2. This scenario is reminiscent of the reductions in H2AK119ub1 at PRC1 targets that occur after deletion of Rybp in mESCs (32), suggesting that Yaf2 exert a function similar to that of Rybp in mESCs.
Recent literature has revealed that the properties of self-renewal and pluripotency are governed by an intricate set of cellintrinsic programs in response to the extracellular signals, including developmental cues (33). These intracellular signaling cascades ultimately converge to the chromatin directly (i.e. via histone phosphorylation) or indirectly (via post-translational modifications of chromatin-modifying complexes) to enable gene regulation (33,34). Recent observations have shown that PcG proteins, as central players in above epigenetic programming events, are targeted for post-translational modi-fications such as phosphorylation (17,19,20,35). However, the functional relevance of these phosphorylation events has remained largely unknown. Here, we report Yaf2 was found to be phosphorylated at serine 166 by, at least partially, Gsk3. Intriguingly, Yaf2 Ser-166 phosphodefective mutants were still assembled into PRC1 complexes and did not significantly differ in its promoter tethering activity from WT Yaf2. However, this mutant failed to rescue H2AK119ub1 levels at Yaf2-specific target promoters and thus lost its ability to repress target gene expression in null mESCs. Collectively, these findings suggest Yaf2 is a phosphoprotein that may integrate diverse signals and directly converts them into chromatin modifications (i.e. H2AK119ub1) that control developmental gene expression and the regulation of mESC maintenance and lineage commitment Figure 6. Phosphorylation of Yaf2 promotes H2AK119 monoubiquitination and represses target gene expression. A, quantitative RT-PCR analysis to determine mRNA expression of target genes of Yaf2 (Fgf5, Csn3, Gfra1, Otx2, Wnt8a, and Lef1), PRC1.6 (Tex101 and Tdrkh), and canonical PRC1/2 (Eomes and Hoxb8) in mESCs of the indicated genotypes. Hprt1 was used as a control. Expression was normalized to actin. B, anti-FLAG ChIP-qPCR assays were done using Yaf2 Ϫ/Ϫ and Yaf2 Ϫ/Ϫ rescued with FLAG-Yaf2 mESCs. C, ChIP-qPCR was used to analyze the occupancy of H2AK119ub1 and Ring1B on the promoters of Yaf2 target genes in mESCs of the indicated genotypes. Error bars indicate Ϯ S.D. Bar graphs represent the mean of three independent biological repeats. p value: Student's t test analysis. D, model for Yaf2-mediated transcriptional repression. See under "Discussion" for details. (Fig. 6D). Our findings also implicate PcG-mediated epigenetic silencing as a dynamically controlled process.

Yaf2 in embryonic stem cells
Our data demonstrate that phosphorylation of Yaf2 on serine 166 is at least partially regulated through the Gsk3 pathway. Previous studies (35) revealed that phosphorylation of EZH2 on Ser-21 by the Akt/Gsk3 pathway inhibits its methyltransferase activity. We found that Gsk3-mediated phosphorylation of Yaf2 on serine 166 supports Ring1B-mediated H2AK119ub1 and is critical for maintenance of pluripotency in mESCs. Gsk3 signaling has not been fully explored in regulation of stem cell pluripotency, although a combination of small molecule inhibitors for the mitogen-activated protein kinase kinase (MEK) and Gsk3 (2i conditions) has been implicated in the maintenance of mESC pluripotency under minimal growth conditions (36,37). Therefore, our findings that this signaling directly modifies PRC1 activity warrant further investigation of this signaling cue in control of these processes.
It has been shown by transgenic studies that Yy1, the mammalian homologue of Drosophila PHO, can phenotypically correct PHO mutant flies (38) and is required for recruitment of PRC1/2 to certain Polycomb-responsive elements (PRE) in Drosophila (39,40). Because of its potential interaction with Yaf2, it has been proposed that Yaf2 might act as a bridge between Yy1 and a subset of PRC1 to target these complexes to specific DNA sequences (PREs) (41,42). Intriguingly, however, a stable association of Yy1 with PRC1 has never been detected in studies from us and other groups (Fig. 4) (6,26,43). It is possible that the interaction between Yy1 and PRC1 components may be transient. Thus, further experiments are needed to carefully examine the influence of Yy1 on the recruitment of PRC1 to chromatin, at least in mESCs.
Transcriptome analysis reveals that Yaf2 deficiency in mESCs resulted in selective derepression of ectoderm-specific genes (Fig. 3D). Moreover, Yaf2 Ϫ/Ϫ EBs expressed increased levels of ectodermal markers (Fig. 2C). Consistent with these observations, Yaf2 deletion skewed differentiation toward neural ectoderm lineages in a teratoma assay (Fig. S2). Importantly, ChIP analysis revealed direct binding by Yaf2 to target promoters in mESCs (Fig. 6), suggesting a potential requirement for Yaf2 in repression of polycomb-controlled developmental regulators. Given the evidence that Yaf2 regulates gene expression, particularly by repressing essential developmental genes involved in neural ectoderm differentiation, it is possible that subtle deregulation of these developmental regulators interferes with timely progression of neural ectoderm development. Further studies using mice with a conditional Yaf2 mutant in the nervous system will shed light about the in vivo contributions of Yaf2 during nervous system development.
The full-length cDNA of mouse Yaf2 (Q99LW6) was added to a C-terminal FLAG-tag (DYKDDDDK) sequence by PCR and inserted into pBS KS II vector. Deletions and point mutations of the Yaf2 cDNA were generated by QuikChange II sitedirected mutagenesis kit (Agilent). All the verified products were introduced into the lentiviral vector as described previously (26).

Generation of mutant mESCs
Knockout mESCs were generated as described previously (44). Briefly, constructs expressing a pair of sgRNAs were cotransfected along with a puromycin resistance plasmid into mESCs by Lipofectamine 2000 (Life Technologies, Inc.). After 24 h of transfection, mESCs were selected by puromycin for 2 days and then seeded on mitomycin-irradiated feeder cells on gelatinized tissue culture plates in mESC medium. Correctly targeted knockout mESC clones were identified by genomic PCR analysis. Subsequently, the null mESCs were confirmed by Western blotting.

Lentiviral supernatant production and infections
Lentiviral supernatants were prepared as described (26). The mESCs were infected by lentiviral supernatants with Polybrene (Sigma, final concentration of 8 g/ml). After 24 h, positive mESCs were selected by puromycin.

Colony formation, AP staining, and flow cytometry
WT or mutant mESCs (about 1000 cells) were dissociated by trypsin into single-cell suspension and plated onto 10-cm MEFcovered dishes and incubated for 7 days before imaging by microscopy.
AP staining of mESC colonies was done following the manufacturer's instructions (Yeasen). Cell cycle analysis has been described previously (46). Propidium iodide (PI) staining of mESCs was performed according to the manufacturer's instructions (Vazyme). Apoptosis was detected using an annexin V-FITC/PI detection kit (Yeasen) and followed by flow cytometry analysis on a FACS LSRF or tessa (BD Biosciences).

Embryoid body and teratoma formation and analysis
Embryoid body formation was performed as described previously (26). Briefly, trypsinized mESCs were resuspended in Embryoid body medium. A 30-l sample (300 -500 cells) was pipetted onto the 15-cm Petri-dish plate, which was then placed upside down, and the hanging drops were cultured for about 3 days. Then EBs were collected and cultured on a rotating shaker. Total RNA was collected (TRIzol, ThermoFisher Scientific) on days 3, 7, and 12.
Teratoma formation was performed as described previously (26). Briefly, mESCs were trypsinized, washed twice with PBS, and then injected subcutaneously into the flanks of 8-week-old immunodeficient mice. After 4 weeks, mature teratomas were fixed in 4% phosphate-buffered formaldehyde. Paraffin-embedded tissue was sliced and stained with hematoxylin and eosin (H&E) and then immunohistochemistry. The experimental animal facility has been accredited by the Association for Assessment and Accreditation of Laboratory Animal Care International (AAALAC), and the Institutional Animal Care and Use Committee (IACUC) of the Model Animal Research Center of Nanjing University approved all animal protocols used in this study.

Expression analysis
Total RNA was isolated from mESCs or EBs using TRIzol reagent (ThermoFisher Scientific) and was reverse-transcribed according to the manufacturer's protocol. Briefly, RT-PCRs were carried out with the HiScript TM First Strand cDNA synthesis kit (Vazyme). Quantitative RT-PCR was performed in triplicate from independent biological samples using Pow-erUp TM SYBR Green Master Mix (Life Technologies, Inc.) on a StepOne TM software version 2.3 (Applied Biosystems). Expression levels of each gene were normalized to actin. Primers for RT-PCR are listed in Table S1.

Whole-cell lysates, histone extraction, and Western blot analysis
Whole-cell lysates of mESCs were prepared in RIPA buffer as described (16). Proteins were separated via SDS-PAGE and analyzed by Western blotting. Histone extraction was performed as described (47). Briefly, mESCs were resuspended in 1 ml of Hypotonic lysis buffer and incubated for 30 min on a rotator at 4°C. The intact nuclei were spun by centrifuge at 10,000 ϫ g for 10 min at 4°C. The nuclear pellet was resuspended in 200 l of 0.2 N HCl. Histone modifications were analyzed by immunoblotting. List of the antibodies used in this study is shown in Table S2.

RNA-seq analysis
RNA-seq was performed as described previously (16) and supported by Annoroad Gene Tech. Total RNA was isolated using TRIzol reagent (ThermoFisher Scientific). Sequencing libraries were implemented by the Next Ultra RNA Library Prep Kit for Illumina (New England Biolabs), and then three independent biological samples were sequenced on an Illumina Hiseq4000 platform and generated PEI150 strategy. The expression levels of each gene were calculated by fragments per kilobase of transcript per million fragments mapped. Then differential gene expression analysis between WT and knockout mESCs was measured by DESeq (version 1.16). The list of deregulated genes in knockout is shown in Tables S3 and S4. In addition, RNA-seq data were deposited at the Gene Expression Omnibus under accession number GSE113830.

CHIR99021 treatment
Gsk3 inhibitor CHIR99021 was obtained from Sigma. mESC lines were treated with CHIR99021 (6 M) for 24 h.

Dephosphorylation assay
mESCs were resuspended in the dephosphorylation buffer (100 mM NaCl, 50 mM Tris-HCl, pH 7.9, 10 mM MgCl 2 , 1 mM DTT, pH 7.9 at 25°C) with protease inhibitors. Then, CIP (New England Biolabs) was added to the lysates and incubated up to 60 min at 37°C. The reaction was stopped by the addition of phosphatase inhibitors (10 mM sodium fluoride, 10 mM sodium vanadate, 10 mM potassium pyrophosphate, and 5 mM sodium phosphate). The dephosphorylated lysates were processed for SDS-PAGE.

Immunoprecipitation and mass spectrometry (MS)
Immunoprecipitation and MS were performed as described previously (16). The nuclear lysates of Yaf2 Ϫ/Ϫ mESCs rescued with FLAG-Yaf2 were used for immunoprecipitation by M2 beads. For MS identification, the products were run into the SDS-polyacrylamide gel about 1-2 cm, followed by Coomassie Blue staining. The gel slice, including all proteins, was cut and identified by LC-MS analysis (2D-NanoLC/TripleTOF5600).

ChIP
ChIP analysis was performed as described previously (26). The proteins and DNAs of mESCs were cross-linked in 1% fresh formaldehyde solution diluted from 37% formaldehyde (Sigma) followed by 0.125 M glycine to quench. Cells were lysed in ChIP buffer, and then sonication was performed using Bioruptor system (Diagenode) to produce DNA fragments of 200 -1000 bp. For immunoprecipitation, 1 volume of chromatin solution was mixed with 4 volumes of ChIP dilution buffer and appropriate antibodies. Immunocomplexes were captured by protein A/G-Sepharose overnight at 4°C. The beads were washed twice by low-salt buffer and high-salt buffer, respectively. Immunocomplexes were eluted from beads using eluent buffer, and then the DNAs were decross-linked and purified using DNA gel extraction kit (Axygen). Quantitative PCRs were performed using Pow-erUp TM SYBR Green Master Mix (ThermoFisher Scientific) on a StepOne TM software version 2.3 (Applied Biosystems). The enrichment was calculated as 2 Ϫ⌬Ct , where ⌬Ct ϭ Ct (ChIP) Ϫ Ct(input). Primers utilized for ChIP are listed in supporting Table  1.