Polycomb-Like 3 Promotes Polycomb Repressive Complex 2 Binding to CpG Islands and Embryonic Stem Cell Self-Renewal

Polycomb repressive complex 2 (PRC2) trimethylates lysine 27 of histone H3 (H3K27me3) to regulate gene expression during diverse biological transitions in development, embryonic stem cell (ESC) differentiation, and cancer. Here, we show that Polycomb-like 3 (Pcl3) is a component of PRC2 that promotes ESC self-renewal. Using mass spectrometry, we identified Pcl3 as a Suz12 binding partner and confirmed Pcl3 interactions with core PRC2 components by co-immunoprecipitation. Knockdown of Pcl3 in ESCs increases spontaneous differentiation, yet does not affect early differentiation decisions as assessed in teratomas and embryoid bodies, indicating that Pcl3 has a specific role in regulating ESC self-renewal. Consistent with Pcl3 promoting PRC2 function, decreasing Pcl3 levels reduces H3K27me3 levels while overexpressing Pcl3 increases H3K27me3 levels. Furthermore, chromatin immunoprecipitation and sequencing (ChIP-seq) reveal that Pcl3 co-localizes with PRC2 core component, Suz12, and depletion of Pcl3 decreases Suz12 binding at over 60% of PRC2 targets. Mutation of conserved residues within the Pcl3 Tudor domain, a domain implicated in recognizing methylated histones, compromises H3K27me3 formation, suggesting that the Tudor domain of Pcl3 is essential for function. We also show that Pcl3 and its paralog, Pcl2, exist in different PRC2 complexes but bind many of the same PRC2 targets, particularly CpG islands regulated by Pcl3. Thus, Pcl3 is a component of PRC2 critical for ESC self-renewal, histone methylation, and recruitment of PRC2 to a subset of its genomic sites.

Whereas H3K27me3 is associated with repressed genes, H3K4me3 marks active genes. ESCs and a number of adult stem cells, however, contain a unique chromatin signature, termed bivalency, that is comprised of both H3K27me3 and H3K4me3 marks [8][9][10][11][12][13][14][15]. Many bivalent domains are at CpG islands, domains of DNA with elevated GC content that display low levels of DNA methylation. CpG islands are commonly found at vertebrate promoters and are associated with 70% of annotated genes including most housekeeping genes and many developmentally regulated genes [16][17][18]. CpG-rich domains commonly display H3K4me3, but GC-rich sequences also promote H3K27me3, creating opposing marks within the same domain [19,20]. By occupying CpG islands and marking them as bivalent domains in ESCs, PRC2 may keep the associated genes repressed but poised for rapid activation upon differentiation [11]. How PRC2 is recruited to CpG islands is not known.
Here, we show that mouse Pcl3 interacts with the core components of PRC2 and promotes complex function. By depleting Pcl3 in ESCs, we demonstrate that Pcl3 contributes to ESC self-renewal, but not differentiation, of the three germ layers. Using ChIP-seq of Pcl3 shRNA-treated cells, we show that Pcl3 knockdown causes decreased H3K27me3 and Suz12 binding to the genome, indicating that Pcl3 regulates PRC2 binding at diverse target genes. Furthermore, Pcl3 localizes with Suz12 at a subset of PRC2 targets, including genes and microRNAs associated with differentiation and development. We also show that several Pcl3 Tudor domain residues are necessary for H3K27me3. Finally, we identify two GC-rich binding motifs that are enriched at Pcl3-dependent PRC2 targets, indicating that Pcl3 promotes PRC2 binding at CpG islands. Taken together, these results reveal that Pcl3 is an important regulator of PRC2 at a subset of target genes.

Pcl3 is a component of PRC2
To identify PRC2 binding partners that could contribute to its function, we used the recently developed Floxin system to create a tandem affinity purification (TAP) tagged allele of Suz12 ( Figure  S1A) [54]. In brief, we reverted a Suz12 gene trap (Suz12 Gt/+ ) allele generated in a mouse ESC line to produce an allele that reexpresses Suz12 but that contains a loxP targeting site (Suz12 Rev/+ ). Via a modified Floxin shuttle vector, we inserted an exon encoding amino acids 277-741 of Suz12 fused to a carboxy-terminal 66His-36Flag TAP tag. The resultant allele (Suz12 Suz12TAP/+ ) expressed the full-length TAP-tagged Suz12 from the endogenous locus ( Figure 1A and Figure S1B). We measured protein and mRNA levels of Suz12 in all ESC lines by immunoblot and quantitative reverse transcription PCR (qRT-PCR) ( Figure 1A and Figure  S1C). As expected, Suz12 Gt/+ cells displayed reduced Suz12 expression, whereas Suz12 Rev/+ cells displayed levels restored to wild type amounts ( Figure 1A and Figure S1C). Suz12 Suz12TAP/+ cells displayed moderately increased mRNA and protein levels of Suz12 compared to wild type ( Figure 1A and Figure S1C).
We depleted Pcl3 using multiple Pcl3 shRNA targeting vectors ( Figure 2A and Figure S2A). Upon culturing multiple clones of Pcl3 knockdown ESCs, we observed an increased percentage of cells that were larger, less dense and displayed morphologies consistent with differentiation ( Figure 2B). These differentiated morphologies suggested that Pcl3 may play a role in ESC maintenance or self-renewal.
To assess whether Pcl3 contributes to ESC self-renewal, we examined ESC markers including Oct4, Nanog, and alkaline phosphatase. Consistent with the morphological changes, Oct4 and Nanog expression and protein levels were decreased in Pcl3 knockdown cells compared to control ESCs ( Figure 2C-2D and

Author Summary
Embryonic development requires coordinated changes in gene expression for the differentiation of specific cell types. Regulated changes in gene expression are also important for maintaining tissue homeostasis and preventing cancer. Histone modifications contribute to the control of gene expression by affecting chromatin structure and the recruitment of regulatory proteins. Polycomb repressive complex 2 (PRC2) catalyzes the methylation of a lysine residue on histone H3, an early step in gene repression. By investigating how PRC2 is recruited to genes, we have found that Polycomb-like 3 (Pcl3), a protein upregulated in diverse cancers, is a component of PRC2 that promotes its binding and function at target genes. Consistent with roles for PRC2 in regulating stem cell behaviors, Pcl3 is important for embryonic stem cell self-renewal. Thus, Pcl3 is a critical regulator of gene repression and stem cell self-renewal that acts by controlling PRC2 binding to target genes. Figure S2B). We also assessed scramble and Pcl3 shRNA ESCs for alkaline phosphatase activity. Alkaline phosphatase staining was slightly reduced in Pcl3 shRNA-treated ESCs ( Figure S2C). To quantitate this observation, we used a more sensitive colorimetric assay, which confirmed that Pcl3 knockdown reduces ESC alkaline phosphatase activity ( Figure 2E). To confirm that these phenotypes were specifically due to Pcl3 knockdown, we overexpressed a TAPtagged form of Pcl3 in wild type cells ( Figure S2D-S2E). Overexpression of Pcl3 in ESCs resulted in increased levels of Oct4 and Nanog, further indicating that Pcl3 levels correlate with ESC gene expression ( Figure 2C).
To determine whether decreased expression of Oct4 and Nanog affects self-renewal, we assessed Pcl3 knockdown ESCs for their ability to generate colonies and found that Pcl3 shRNA ESCs formed significantly fewer colonies than control cells ( Figure 2F and Figure S2F). These data suggest that Pcl3 promotes ESC selfrenewal. To substantiate this finding, we assayed the ability of wild type and Pcl3-overexpressing cells to self-renew and form colonies. We challenged wild type and Pcl3-overexpressing cells by growing them in media containing reduced LIF. Pcl3-overexpressing ESCs were able to form colonies modestly but significantly better than wild type cells, further indicating that Pcl3 enhances ESC selfrenewal ( Figure 2G).
In addition to its roles in self-renewal, PRC2 is critical for ESC differentiation and embryonic development [22,24,29,30]. Besides ESCs, Pcl3 is expressed in a number of differentiated tissues from Protein levels of Suz12 and Suz12-TAP were measured in wild type, Suz12 Gt/+ , Suz12 Rev/+ , and Suz12 Suz12TAP/+ cell lines by immunoblot. (B) Proteins detected by mass spectrometry that specifically co-purified with Suz12-TAP, their symbol, unique hits, and percent coverage. (C) Pcl3-V5 binds to Suz12-TAP. Suz12 Suz12TAP/+ ESCs were transfected with empty vector, Pcl3-V5, or Mks1-V5 (control), and lysates were immunoprecipitated with FlagM2 and probed with anti-V5. (D) Pcl3-V5 binds Suz12, Ezh2, and Eed. Lysates from ESCs transfected with empty vector, Pcl3-V5, and Mks1-V5 (control) vectors were subjected to immunoprecipitation with anti-V5. Samples were then probed with anti-Suz12, anti-Ezh3, and anti-Eed. All westerns and co-immunoprecipitations were performed three times. doi:10.1371/journal.pgen.1002576.g001 all three germ layers including bone, spleen, and prostate [41]. Pcl3 is not well expressed in the adult nervous system, but it expressed in the head during development [41]. To ascertain whether Pcl3 is essential for ESC differentiation, we performed teratoma and in vitro differentiation assays and assessed the presence of all three germ layers. Teratomas from scramble and Pcl3 shRNA-expressing cells were weighed and histologically examined. While Pcl3 knockdown was maintained, we did not observe differences in the size of the teratomas nor in the formation of endodermal, ectodermal, or mesodermal derivatives (data not shown and Figure 2H and Figure S2G). Moreover, immunofluorescent staining revealed the presence of neural tissue (NeuN and Tuj1), muscle (Actin), and skin basal cells (K14) in both Pcl3 expressing and Pcl3 knockdown teratomas ( Figure S2H).
In vitro differentiation of scramble and Pcl3 knockdown ESCs to embryoid bodies resulted in a similar finding. Pcl3 shRNA-treated EBs maintained Pcl3 knockdown and expressed markers of neuroectoderm (Nestin), mesoderm (T-brachyury), and endoderm (Hnf4) at least as well as control EBs, indicating that depletion of Pcl3 does not abrogate the ability of ESCs to differentiate to cell types of all three germ layers ( Figure S2I). Thus, Pcl3 enhances ESC self-renewal, but is not critical for ESC differentiation to the three germ layers.
To confirm that Pcl3 shRNA-associated phenotypes were due to decreased Pcl3 levels, we created a construct encoding a TAPtagged Pcl3 not recognized by the Pcl3 shRNAs. Pcl3-TAP was detected by immunoblot and immunofluorescence and was able to bind core PRC2 components as assessed by co-immunoprecipitation, indicating that Pcl3-TAP functions similarly to wild type Pcl3 ( Figure 3E-3F and Figure S3E). Pcl3 expression in Pcl3-TAPexpressing ESCs was nearly 60% of endogenous levels as measured by qRT-PCR ( Figure 3G). Re-expression of Pcl3-TAP increased H3K27me3 levels in Pcl3 knockdown ESCs ( Figure 3H). These data indicate that Pcl3 promotes the formation of H3K27me3 by PRC2.

Pcl3 promotes PRC2 binding and function at a subset of PRC2 targets
To elucidate which genes depend upon Pcl3 for H3K27me3 formation and to investigate whether decreased H3K27me3 is a result of Pcl3-dependent PRC2 binding, we performed ChIP-seq for H3K27me3 and Suz12 in Suz12 Suz12TAP/+ cells expressing either Pcl3 or control shRNA (GSE28325). Reads were aligned to the mouse reference genome (NCBI Build 37, MM9) using bowtie, and only uniquely aligned reads were considered for further analysis ( Figure S4A-S4B) [62].
Suz12 ChIP-seq was performed using the Flag tag of Suz12-TAP. To validate this approach, we compared ChIP-seq of Suz12-TAP using anti-Flag antibodies with ChIP-seq using anti-Suz12 antibodies in Suz12 Suz12TAP/+ cells and found they were equivalent. We, thus, used FlagM2 to detect Suz12 binding in all the following experiments. ChIP-seq with FlagM2 antibodies detected 95% of previously reported Suz12 binding sites, indicating that this approach accurately measured Suz12-TAP localization in the genome ( Figure S4C) [63,64]. In addition, we identified a number of novel Suz12 binding sites, typified by less significant Suz12 ChIP-seq peaks ( Figure S4C). These additional sites may reflect weak Suz12 binding or may be an artifact of increased Suz12 levels.
Comparison of H3K27me3 ChIP-seq in scramble and Pcl3-shRNA treated ESCs revealed a widespread decrease in H3K27me3 levels in Pcl3 shRNA ESCs ( Figure 4A). The decrease in H3K27me3 as assessed by ChIP-seq was similar to, but not as  profound as was detected biochemically. This may be attributable to the inability of ChIP-seq to detect H3K27me3 in regions of highly repetitive sequence, areas that are largely transcriptionally silent. At approximately 85% of sites showing decreased H3K27me3, Pcl3 knockdown cells also displayed reduced Suz12 binding compared to scramble controls, thus correlating decreased Suz12 occupation with decreased PRC2 function ( Figure 4A). Indeed, Pcl3 knockdown significantly reduced Suz12 binding at approximately 65% of PRC2 targets ( Figure 4A, 4C). Sites with decreased Suz12 binding outnumbered those with reduced H3K27me3, indicating that Pcl3 knockdown preferentially affected Suz12 ( Figure 4A).
The regions with the most decreased Suz12 occupation following Pcl3 knockdown were among the most significant Suz12 ChIP-seq peaks in control cells ( Figure S4D). In contrast, less significant Suz12 peaks were more likely to be unaffected by Pcl3 knockdown ( Figure S4D). The finding that the role of Pcl3 is most pronounced on the most significant Suz12 binding sites is consistent with a role for Pcl3 in mediating the recruitment of additional Suz12 to regions of low level Suz12 binding. Alternatively, changes in Suz12 binding may be more difficult to detect at these less significant peaks.
Knockdown of Pcl3 diminished Suz12 binding to chromatin particularly in areas of high gene density ( Figure 4B). Accordingly, chromosome 11, the chromosome with the highest gene density, showed the most profound decrease in Suz12 binding and H3K27me3 upon Pcl3 knockdown ( Figure 4A and Figure S4E). Consistent with PRC2 binding at repressed genes, sites with decreased Suz12 binding in Pcl3 shRNA-treated ESCs were inversely correlated with sites bound by the activating transcription factors E2f1, c-Myc, Zfx, Klf4, and Ctcf, ( Figure S4F and data not shown, Fisher test p-value = 9.5610 253 , 3.7610 215 , 1.2610 213 , 4.0610 26 , 3.6610 26 , respectively) [63]. This is in agreement with previous reports demonstrating nearly mutual exclusion of Suz12 binding with regions bound by this group of transcription factors [63].
Pcl3 co-localizes with Suz12 at many PRC2 targets to promote H3K27me3 As Pcl3 binds core components of PRC2 and promotes Suz12 localization to its target genes, we investigated whether Pcl3 colocalizes with Suz12 at target loci. To assess the distribution of Pcl3 binding to chromatin in ESCs, we performed ChIP-seq of Pcl3-TAP. Nearly all of the regions bound by Pcl3-TAP overlapped with Suz12 targets, suggesting that Pcl3 co-localizes at target genes as part of PRC2 ( Figure 5A and S5A-S5B, Fisher test p-value = 1.5610 2306 ). Quantification revealed that Pcl3-TAP bound nearly half of PRC2 targets, indicating that Pcl3 associates with a subset of PRC2 complexes ( Figure 5A). To substantiate the co-localization analysis, we assessed which Suz12 sites overlapped with Pcl3. Pcl3 and Suz12 co-localized predominantly at the most significant Suz12 binding sites ( Figure 5B). Suz12 sites at which Pcl3 did not co-localize were typified by less significant Suz12 binding ( Figure 5B). In addition, 80% of Pcl3 and Suz12 sites overlapped with CpG islands compared to 65% of sites bound only by Suz12, indicating that Pcl3 and Suz12 preferentially bind at CpG islands (p-value = 2.2610 216 ) [13].
To establish whether Pcl3 may directly contribute to PRC2 binding and function, we compared regions of Pcl3 binding to areas of Pcl3-dependent Suz12 chromatin occupation and H3K27me3 formation. Upon Pcl3 knockdown, 84% of Pcl3 binding regions showed decreased Suz12 binding, whereas less than half of regions not bound by Pcl3 showed decreased Suz12 binding ( Figure 5C-5E). In addition, Pcl3 binding regions showed a near two-fold enrichment for decreased H3K27me3 upon Pcl3 knockdown, as compared to regions not bound by Pcl3 ( Figure 5E). ChIP-qRT-PCR for Pcl3 confirmed that Pcl3 bound to many of the same genes that exhibited decreased Suz12 binding and H3K27me3 upon Pcl3 knockdown ( Figure 4F and Figure 5F). These data indicate that genes bound by Pcl3 are nearly twice as likely to have reduced Suz12 binding and H3K27me3 upon Pcl3 knockdown. Thus, Pcl3 not only co-localizes with Suz12 at many PRC2 targets, but also promotes Suz12 binding and PRC2 function.

Pcl3 regulates gene expression at a subset of PRC2 targets
Inhibition of PRC2 core components and complete abrogation of H3K27me3 dramatically alter gene expression, whereas inhibition of accessory proteins such as Jarid2 result in more modest changes [22,29,30,40]. To test whether Pcl3 affects gene expression, multiple scramble and Pcl3 shRNA clones were analyzed by microarray. Nearly 130 genes displayed altered expression upon Pcl3 knockdown. Unexpectedly, given that PRC2 functions to repress gene expression, nearly three times as many genes were down-regulated upon Pcl3 knockdown as were upregulated (Table S1). Notably, we identified more than half of the upregulated genes as PRC2 targets, which correlated well with depleted Suz12 binding and H3K27me3 ( Figure S6B). We confirmed the microarray data by measuring gene expression by qRT-PCR for several genes ( Figure 6SA).
Given that Pcl3 knockdown increases the proportion of differentiated cells, we tested whether the presence of differentiated cells masked changes in ESC gene expression by pre-plating twice to remove differentiated cells. Microarray analysis revealed that approximately 120 genes showed differential expression in the Pcl3 shRNA cells compared to control after removal of differentiated cells (Table S2). qRT-PCR assessment of a subset of the affected genes confirmed the microarray results ( Figure  S6B). Of these Pcl3-regulated genes, almost two thirds were up-regulated upon Pcl3 depletion. While these results demonstrate that Pcl3 participates in the repression of genes in ESCs, only a few of these genes also displayed decreased Suz12 binding and H3K27me3 levels (i.e., Hand1, Acta1, and Pmp22) (Table S2, Figure 4E-4F and Figure S6B) which may be attributable to the necessity of culturing Pcl3 shRNA ESCs for an extended time before subjecting them to microarray analysis.

Pcl3 regulates Suz12, but not complex stability
Loss of core PRC2 components can destabilize the complex [25,29]. To assess whether loss of Pcl3 affects complex stability, we tested whether the complex could form in the absence of Pcl3. All components, Suz12, Ezh2, and Eed, were found by coimmunoprecipitation to bind each other in both scramble and Pcl3 shRNA-expressing cells, suggesting that Pcl3 is not required for assembly or stabilization of the core complex ( Figure 6A).
To determine whether depleting Pcl3 influences levels of PRC2 components, we assessed the quantity of PRC2 core components in Pcl3 shRNA-treated ESCs. Interestingly, suppression of Pcl3 significantly lowered Suz12 mRNA and protein levels without affecting levels of Eed or Ezh2, suggesting that Pcl3 regulates Suz12 ( Figure 6B-6C). Overexpressing Pcl3 caused a concomitant increase in Suz12 mRNA and protein levels ( Figure 6D-6E). As Suz12 is not a known PRC2 target and Pcl3 is not enriched at the Suz12 locus, effects on Suz12 expression by Pcl3 are likely indirect ( Figure S6D) [7,75,76].

The Tudor domain of Pcl3 is critical for its function
Polycomb-like proteins contain three domains, an aminoterminal Tudor domain and two carboxy-terminal PHD fingers ( Figure 7A). Although the carboxy PHD finger of Pcl2 promotes PRC2 recruitment, the Tudor domain and amino PHD finger of Pcl2 are not required for its regulation of PRC2 [41,77]. In contrast, the Tudor domain and carboxy-terminal PHD finger of human Pcl3 are important for binding Ezh2 and self-association [53].
NMR studies indicate that the Drosophila Polycomb-like Tudor domain lacks a complete aromatic cage, needed to bind methylated lysines or arginines on histones [78]. Interestingly, a comparison of Drosophila Pcl Tudor domain sequence to those of mammalian Pcl Tudor domains suggests that mammalian Polycomb-like proteins may be able to form complete aromatic cages [78]. In particular, the conserved residues W48,Y54, F72, D74, and S76 are implicated in histone binding ( Figure 7A-7B) [78]. To ascertain if these Tudor domain residues are required for Pcl3 function, we mutated two sets (W48;Y54 and F72;D74), as well as a control set of residues (N75;Y78) within Pcl3-TAP ( Figure 7A-7B) [78]. These mutant forms of Pcl3 were expressed in Pcl3 knockdown ESCs to see if they could rescue H3K27me3 levels as well as wild type Pcl3-TAP ( Figure 7C). Similar to wild type Pcl3-TAP, Pcl3-TAP N75S;Y78S was able to promote H3K27me3 formation ( Figure 7D). In contrast, Pcl3-TAP W48A;Y54A and Pcl3-TAP F72S;D74S did not restore H3K27me3 formation ( Figure 7D). Thus, W48;Y54 and F72;D74, two pairs of residues implicated in histone binding, are necessary for Pcl3 function, while the control N75;Y78 pair is dispensable.
Pcl3 and Pcl2 are part of distinct PRC2 complexes with overlapping targets While Pcl1 is minimally transcribed in ESCs, Pcl2 is highly expressed and binds a subset of PRC2 target genes [41]. To assess whether depletion of Pcl3 affects Pcl1 or Pcl2, we measured Pcl1 and Pcl2 expression in scramble and Pcl3 knockdown ESCs. Pcl1 was expressed at extremely low levels but showed no difference upon Pcl3 knockdown ( Figure 8A). Likewise, we found that Pcl2 levels were similar between scramble and Pcl3 shRNA-treated ESCs ( Figure 8A).
Although Pcl2 and Pcl3 are both expressed in ESCs, it is unclear if they participate in same complex or regulate the same genes. To assess whether Pcl3 and Pcl2 incorporate into the same complex or exist in distinct PRC2 complexes, we tested Pcl3 and Pcl2 association by co-immunoprecipitation. Flag immunoprecipitation confirmed that Pcl2 and Pcl3 both bind Suz12 ( Figure 8B). However, Pcl3-TAP did not immunoprecipitate Pcl2, suggesting that Pcl2 and Pcl3 do not associate within the same PRC2 complex ( Figure 8B). To substantiate these data, we performed the reciprocal immunoprecipitation and found that Pcl2 immunoprecipitated Suz12-TAP but not Pcl3-TAP. These findings indicate that Pcl2 and Pcl3 both interact with the core PRC2 complex, but not with each other, suggesting that Pcl2 and Pcl3 participate in distinct PRC2 complexes.
To elucidate whether Pcl3 and Pcl2 bind distinct or overlapping sets of targets, we compared our Pcl3 ChIP-seq data with previously generated Pcl2 ChIP-seq data and found that Pcl3 and Pcl2 binding overlaps at 60-70% of sites [41]. This extensive overlap raises the possibility that many genes may be regulated by both Pcl3 and Pcl2.
Like Pcl3, Pcl2 can promote PRC2 activity at some loci and in some cell types [41,47,77,79]. To test whether Pcl2 can compensate for Pcl3 function, we analyzed whether Suz12 binding at Pcl2 targets was affected by Pcl3 depletion. Of sites bound by both Suz12 and Pcl2, 86% showed decreased Suz12 binding upon Pcl3 depletion ( Figure 8C-8D), suggesting that Pcl2 cannot compensate for Pcl3 function at Pcl2 sites. Unlike Pcl3, Pcl2 can either inhibit or promote PRC2 function [41,47,77,79]. To further investigate the functional relationship between Pcl2 and Pcl3, we assessed Suz12 binding and H3K27me3 at regions bound by both Pcl2 and Pcl3. If Pcl2 and Pcl3 act redundantly to promote Suz12 binding, sites bound by both would not show significant decreases in Suz12 binding upon Pcl3 knockdown as Pcl2 would be functional at these sites. Conversely, if Pcl2 hinders PRC2 binding while Pcl3 promotes PRC2 binding, Suz12 binding at shared Pcl2 and Pcl3 sites would be more dramatically reduced upon Pcl3 knockdown, as Pcl3 would no longer be able to counterbalance Pcl2 function. Among regions bound by Suz12, Pcl2, and Pcl3, 95% showed Suz12 depletion upon Pcl3 knockdown ( Figure 8C-8D). Moreover, Pcl3dependent Suz12 sites often co-localized with Pcl2 and Pcl3 binding, whereas Pcl3-independent Suz12 sites tended to be more distant from Pcl2 and Pcl3 binding sites ( Figure 8E). In addition, 88% of Pcl2 and Pcl3 shared sites showed diminished H3K27me3. Thus, these data suggest that not only does Pcl2 not compensate for Pcl3, but also that Pcl2 may hinder Suz12 binding and H3K27me3 at Pcl2 and Pcl3 shared binding sites.

Pcl2 and Pcl3 co-localize with Suz12 at CpG islands
Polycomb repressive elements (PREs), binding sites that recruit Polycomb proteins, have been well-described in Drosophila, but few have been found in mammals [80][81][82][83][84]. Nevertheless, mammalian PRC2 is known to preferentially bind CpG islands [13,19,85]. Similarly, we found that 86% of Pcl2 and Pcl3 colocalization sites were found in CpG islands [13,41,86]. To assess whether distinct DNA sequence motifs required for Pcl2 and Pcl3 recruitment exist, we searched for sequence features that can discriminate DNA sequences in genomic regions bound by both Pcl2 and Pcl3. We scanned 500 base pair regions surrounding Pcl3 ChIP-seq peaks for prevalent DNA sequences and identified two enriched 10-and 14-mer motifs ( Figure 9A). Motif 1 consisted of three consecutive G-C dinucleotides with two nucleotides on either side containing modest GC enrichment. Motif 2 contained two G-C-rich regions separated by seven nucleotides with little to no GC enrichment. Thus, motifs commonly found in Pcl2 and Pcl3 binding regions contain short GC sequences surrounded by dissimilar sequence.
To assess whether these motifs were predictive of Pcl2 and Pcl3 binding sites, Pcl3-dependent Suz12 binding sites or CpG islands, we constructed a position specific scoring matrix (PSSM) for each motif. The occurrence of motifs 1 and 2 correlated more tightly with Pcl3-dependent Suz12 binding sites that co-localize with both Pcl2 and Pcl3 binding sites than with Suz12 sites that are not associated with Pcl2 and Pcl3 binding, nor depend on Pcl3 activity ( Figure 9B-9C). In addition to the maximum PSSM scores, we also computed the CpG density of each site and used these three sequence features to build a support vector machine classifier to help predict whether specific motifs would co-localize with Pcl2 and Pcl3 binding and areas of Suz12 depletion following Pcl3 knockdown ( Figure 9D). The classifier had cross-validation accuracy of 75%, indicating that CpG density, Motif 1, and Motif 2 together account for the majority of Pcl3-dependent PRC2 binding. Even though the vast majority of Pcl2 and Pcl3 binding sites were found in CpG islands, they seemed to avoid regions with the highest CpG density, concordant with the moderate GC content of the motif sequences ( Figure 9B-9D). These data indicate that Pcl2 and Pcl3 may utilize specific binding motifs to recruit PRC2 to form H3K27me3 at CpG islands.

Discussion
In an effort to understand the mechanisms that regulate PRC2, we have discovered that a PRC2 interacting protein, Polycomblike 3, promotes ESC self-renewal and mediates PRC2 binding at a number of targets genes. In ESCs, Pcl3 promotes the expression of multiple markers of pluripotency. Consistent with a functional role for Pcl3 in PRC2, Pcl3 also promotes H3K27me3 formation. ChIP-seq analyses revealed that diminished Pcl3 decreased Suz12 binding, particularly in regions where Pcl3 co-localizes with Suz12, suggesting that one way that Pcl3 promotes PRC2 activity is by mediating PRC2 binding to chromatin. Pcl3 likely resides in a PRC2 complex distinct from that which incorporates its paralog, Pcl2, but both PRC2 complexes have overlapping targets predominantly at CpG-rich regions. Despite their overlap and homology, Pcl2 cannot compensate for Pcl3 and may inhibit PRC2 at shared genes, as regions of Pcl2 and Pcl3 co-binding show the most significant reduction in Suz12 binding upon Pcl3 knockdown.
Pcl3 likely has a direct role in promoting PRC2 binding at select targets, as it 1) biochemically interacts with PRC2, 2) binds to genomic regions overlapping with Suz12, and 3) is predominantly required for genomic Suz12 binding at Pcl3 binding sites. This cooccurrence and dependence of PRC2 on Pcl3 suggests that the biochemical interaction of Pcl3 with PRC2 underlies its role in PRC2 binding. In addition to this direct role for Pcl3 in PRC2 binding to chromatin, the reduction in Suz12 binding upon Pcl3 knockdown may also be partially attributable to decreased Suz12 levels. Suz12 is not known to be regulated by PRC2 and Pcl3 does not bind the Suz12 locus, indicating that this regulation of Suz12 by Pcl3 is likely to be indirect.
Many bivalent genes are involved in regulating ESC differentiation [8][9][10][11][12][13][14][15]. However, we observed that Pcl3 promotes ESC self-renewal, but is not required for differentiation into cell types of all germ layers. As we have found that Pcl3 promotes H3K27me3 at a subset of bivalent genes, it is possible that Pcl3 predominantly represses bivalent genes involved not in cell fate decisions in differentiating cells, but in genes that promote differentiation itself. Thus, Pcl3 depletion would derepress genes that promote differentiation, leading to specific defects in self-renewal. Conversely, Pcl3 overexpression may lead to increased self-renewal by repressing genes involved in differentiation that inhibit the expression of self-renewal genes such as Nanog.
Consistent with this idea, we found that inhibiting Pcl3 function increased H3K4me3 levels and expression at some but not all bivalent genes. Pcl3-dependent H3K4me3 may regulate be the genes that promote differentiation upon Pcl3 inhibition. Thus, diminishment of Pcl3 does not derepress all PRC2 targets, but instead may derepress a specific subset that restrain ESC differentiation.
PRC2 is recruited to chromatin by multiple mechanisms including interactions with CpG islands, non-coding RNAs, and histone binding proteins [36,38,39,[87][88][89]. Our analysis of Suz12 and Pcl3 binding sites revealed that they associate with CpGenriched regions, particularly two CpG-rich DNA sequence  motifs. These motifs contained GC rich regions surrounded by areas with minimal GC content, indicating that Pcl proteins do not associate with the most GC-rich regions, but rather with regions of moderate GC content and containing motifs 1 or 2. Our support vector machine classifier supported these findings by showing that Pcl2 and Pcl3 binding sites that overlap with Pcl3-dependent Suz12 bound sites localize in areas with moderate CpG density. The finding that Pcl proteins associate with CpG islands is consistent with prior observations that PRC2 binds to CpG islands and demonstrates that Pcl3 may participate in PRC2 recruitment to CpG islands [13,19,85].
Pcl proteins contain multiple domains by which they can associate with chromatin, including a Tudor domain and two PHD fingers implicated in recognizing methylated or unmethylated histone lysines or arginines [90][91][92]. Recent work has implicated several residues of the Pcl Tudor domain in histone binding [78]. Mutating two sets of these Pcl3 Tudor domain residues revealed that they are essential for Pcl3-dependent PRC2 activity. As these mutant forms were able to incorporate into PRC2, they may be acting as dominant negatives and poisoning the complex or as loss-of-function mutants by creating lessfunctional complexes. Thus, these mutants may perturb PRC2 function through inhibition of complex recruitment to histones or by decreasing PRC2 methyltransferase activity.
Two other paralogs of Pcl exist in mammals. Pcl1 is expressed minimally in ESCs, and Pcl2 has been implicated in either promoting or inhibiting PRC2 in a gene-and cell type-dependent manner [41]. In ESCs, Pcl2 restrains PRC2 activity globally, and in MEFs, Pcl2 inhibits PRC2 recruitment to certain loci [41,47,77]. In an effort to understand the relationship between these paralogs, we found that Pcl2 and Pcl3 do not associate and likely participate in separate PRC2 complexes, which both bind a subset of PRC2 targets. If Pcl2 and Pcl3 both promoted PRC2 binding and activity at these targets, we might expect Pcl2 to compensate in Suz12 binding and H3K27me3 upon loss of Pcl3. Instead, we observed the greatest dependence of Suz12 binding on Pcl3 at regions also bound by Pcl2, suggesting that Pcl2 may inhibit PRC2 binding at these genes ( Figure 9E).
It is unlikely that Pcl2 requires Pcl3 for activity as knockdown of Pcl2 and Pcl3 in ESCs results in opposite phenotypes, further suggesting that these two paralogs have distinct and opposing functions. For example, Pcl2 destabilizes Ezh2, whereas Pcl3 enhances Suz12 protein levels without affecting Ezh2 [41]. Moreover, Pcl2 and Pcl3 have opposite effects on ESC selfrenewal and differentiation. Pcl2 inhibits expression of Oct4, Sox2, and Nanog and promotes ESC differentiation [41]. In contrast, Pcl3 promotes Oct4 and Nanog levels and ESC self-renewal. Finally, whereas the Tudor domain of Pcl2 is dispensable for its function, the Pcl3 Tudor domain is essential for mediating PRC2 activity [41].
Despite its global requirement to restrain H3K27me3 levels in ESCs, Pcl2 can promote H3K27me3 and PRC2 binding at certain sites [41,47]. Pcl2 and Pcl3 overlap at approximately 65% of sites, indicating that each binds separately at a third of their sites. While our data indicates that Pcl3 and Pcl2 function oppositely at Pcl2 and Pcl3 targets, one possibility is that Pcl2 promotes PRC2 activity at regions devoid of Pcl3 ( Figure 9E) [41].
Based on these data, we propose that Pcl2 and Pcl3 have opposing functions on PRC2 binding and activity at sites they coregulate ( Figure 9E). At sites that Pcl2 and Pcl3 individually regulate, each may both promote PRC2 binding. Upon loss of Pcl3, sites dependent upon Pcl3 for promoting PRC2 function lose H3K27me3, and this loss may be exacerbated by negative regulation of PRC2 by Pcl2 at sites that Pcl2 and Pcl3 regulate together ( Figure 9E). Thus, Pcl2 and Pcl3 may function antagonistically at shared sites, but promote H3K27me3 at sites regulated by a single Pcl ( Figure 9E). This model predicts that the minority of sites at which Pcl2 promotes H3K27me3 levels will be regions not bound by Pcl3. The model does not resolve how Pcl2 could promote PRC2 activity at regions where Pcl3 does not bind, but could inhibit PRC2 activity at shared sites upon Pcl3 knockdown. One possibility is that regions of Pcl2 and Pcl3 cooccupancy are characterized by an inhibitory PRC2 activity independent of Pcl2.
Both bivalency and gene repression at CpG promoters are dependent on PRC2 activity. In ESCs, CpG promoters may be poised for activation, reflected by their ability to bind RNAPII even when inactive [93]. Pcl3-mediated deployment of PRC2 at CpG islands may be one mechanism by which these poised CpG promoters are kept in check in ESCs. When cell type-specific promoters are being repressed during development, Pcl3 may promote PRC2 binding at many of these same CpG islands in ways that recruit PRC1 to ensure more permanent transcriptional silencing.
Understanding how Pcl proteins regulate PRC2 at specific genes may also be critical to elucidating gene misregulation during cancer. Many cancers show a global silencing at gene promoters, particularly at CpG islands [94,95]. Ezh2 is one of the most commonly misregulated genes in cancer, and recent work has shown that Pcl3 is misregulated in diverse cancers as well [34,96]. Similar to our findings in ESCs, Pcl3 upregulation in cancer cells may promote PRC2 recruitment, inhibiting the expression of genes that inhibit self-renewal, and leading to inappropriate cell proliferation. Thus, studying the role of Pcl3 in PRC2 function in ESCs may provide insight into how Pcl3 up-regulation returns cancer cells to transcriptional states and self-renewal capacities similar to those of ESCs.

Tissue culture
Mouse E14 embryonic stem cells were grown as described previously [54,97]. Bay Genomics gene trap line XG122 was used for the establishment of Suz12 Suz12TAP/+ ESCs [54].

Lentiviral infection
ESCs were trypsinized and 1 ml of cells at 10 5 /ml were plated into one well of a 6-well plate along with 15 ml of concentrated lentivirus (Open Biosystems RMM4534 and Sigma SHC002). The plate was shaken every ten minutes for one hour and then incubated overnight at 37uC. Media was changed the next day and selection was started (2 mg/ml puromycin, 50 mg/ml zeocin). Cells were selected for six days, clones were picked, and immunoprecipitation. (D) Pcl3-TAP W48A;Y48A and Pcl3-TAP F72S;D74S do not support H3K27me3 whereas Pcl3-TAP N75S;Y78S does. Immunoblot displaying H3K27me3 levels in histones purified from cells expressing Pcl3 shRNA, Pcl3 shRNA with Pcl3-TAP, and Pcl3 shRNA with mutant Pcl3-TAP. Histone H3 and a-tubulin were used as loading controls. Each lane represents a different clone. All westerns and immunoprecipitations were performed 3-4 times with at least two clones from each mutant. doi:10.1371/journal.pgen.1002576.g007 knockdown was assessed. Five different Pcl3 shRNA lentiviral constructs were tested and all gave similar results.

Transfection
Cells were transfected using Lipofectamine 2000 (Invitrogen 11668-019). Briefly 0.5 mg of DNA was combined with 1 ml of Lipofectamine in 50 ml of Optimem and incubated for 20 min. ESCs were trypsinized and 2610 5 were plated into a well of a 6well along with Lipofectamine mixture. Cells were incubated for 18 hrs followed by media replacement. For siRNAs, 300 ng siRNA and 3 ml Lipofectamine in 300 ml of Optimem was used for the Lipofectamine mixture.

In vitro prepared siRNA pools
Libraries of siRNAs were made as previously described [56,57].

Purification and mass spectrometry of Suz12
Suz12-TAP and bound proteins were purified and subjected to mass spectrometry as previously described [79].

Histone purification
Histones were purified from ESCs as previously described [98].

Immunoprecipitation and Western immunoblot
Cells were treated as previously described [99]. For Flag immunoprecipitation, 40 ml of FlagM2 resin (Sigma A2220) was added to 800 mg protein lysate overnight. All experiments were performed three or more times. Pcl3 knockdown experiments contained at least two scramble and Pcl3 shRNA clones.

Antibodies
The list of antibodies used for immunoblotting and immunofluorescence is included in Text S1.

Immunofluorescence and staining
Cells were plated on coverslips coated with poly-lysine and Matrigel (BD Biosciences) and stained as previously described [99].

Expression analysis
Expression was assayed as described previously [99]. Expression levels are representative of three or more experiments using 2-6 different clones and assayed in quadruplet. Statistical significance was determined using an unpaired student's t-test. Primers for qRT-PCR and ChIP is included in Text S1.

DNA mutagenesis
Point mutations were created using the Agilent site-directed mutagenesis kit (Quikchange II XL Site-Directed Mutagenesis, Agilent Technologies, 200522) and Quikchange primer design application.

Alignment
Polycomb-like paralog alignment was performed using CLC sequence viewer (CLC Bio).

Alkaline phosphatase activity assay
Alkaline phosphatase activity was determined by staining and using StemTag Alkaline Phosphatase Activity Kit (Cell Biolabs CBA-301). ESCs were plated with ES media containing LIF and allowed to grow to 80% confluency. For staining, ESCs were fixed in 0.5% glutaraldehyde and washed three times in HBS followed by three washes in AP buffer (100 mM Tris pH 9.5, 100 mM NaCl, 50 mM MgCl 2 ). NBT/BCIP substrate was added and samples were incubated at 37u for 10 minutes. Reaction was stopped with 50 mM EDTA pH5.1 in PBS, and pictures were taken. For the StemTag colorimetric assay, cells were lysed, incubated with assay substrate, and assayed for absorbance. Experiment was performed in duplicate three independent times.

Colony formation assay
ESCs were trypsinized, counted, and plated at 100 cells/well of a 6-well. After 8 days, colonies were stained with methylene blue and quantitated. For wild type and Pcl3-overexpressing cells, ESCs were grown in 5% of normal LIF concentration. Experiment was performed in duplicate four independent times.

Teratoma formation
ESCs were trypsinized, counted, and prepared to a concentration of 6.7610 6 cells/ml. 300 ml (2610 6 ) cells were injected subcutaneously into severely combined immunodeficient (SCID) mice and allowed to form teratomas for 2-3 weeks. Tumors were removed and divided for paraffin sections, frozen sections, and RNA. Teratomas were formed from two clones each of scramble and Pcl3 shRNA ESCs.

Microarray
Sample preparation, labeling, and array hybridizations were performed according to standard protocols from the UCSF Shared Microarray Core Facilities and Agilent Technologies (http:// www.arrays.ucsf.edu and http://www.agilent.com). See Text S1 for more detailed description. The microarray was performed in 2-6 clones from scramble and Pcl3 shRNA ESCs.

ChIP
ChIP-seq experiments were done as previously described with 500 mg of chromatin and 5 mg of antibody [100]. ChIP-qRT-PCR experiments were performed at least two times with 10 7 ESCs and 10 mg of antibody. Antibodies used for immunoprecipitation Figure 8. Pcl3 and Pcl2 participate in separate PRC2 complexes but overlap at many PRC2 binding sites. (A) Pcl2 and Pcl1 mRNA levels are not significantly different between scramble and Pcl3 shRNA cells. Graph represents average expression from three different experiments in 2-6 different clones and assayed in quadruplet. Error bars indicate standard deviation. (B) Co-immunoprecipitation showing that immunoprecipitated Suz12-TAP associates with Pcl2 whereas immunoprecipitated Pcl3-TAP does not. Pcl3-TAP does bind Suz12 as seen in the lower blot. a-tubulin was used as a loading control. Reciprocal co-immunoprecipitations showing a Pcl2-specific association with Suz12-TAP but not an association between Pcl2 and Pcl3-TAP. b-actin was used as a loading control. Each blot is representative of three different experiments. (C) Graph depicts chromosomewise distributions of Pcl2 co-localization with sites depleted of Suz12 following Pcl3 knockdown (blue). Regions containing only Suz12 depletion and only Pcl2 binding are indicated in cyan and orange respectively. (D) Decreased Suz12 binding upon Pcl3 knockdown occurs at 86% of targets bound by Suz12 and Pcl2 (Fisher test p-value,p,10 2300 ). Of sites bound by Pcl2 and Pcl3, 95% show decreased Suz12 binding upon Pcl3 knockdown. (E) Pcl2 and Pcl3 binding sites co-localize with regions of Suz12 depletion upon Pcl3 knockdown (top), whereas Pcl2 and Pcl3 binding sites reside far from unaffected Suz12 binding sites (bottom). The figure shows the joint probability density of the nearest Pcl2 and Pcl3 locations relative to Suz12. Darker color represents higher probability density. Bivariate normal kernel density estimates were obtained by using the smoothScatter function in R. The distance coordinates were transformed to base pairs sign (Pcl22Suz12) log 10 |Suz122Pcl2|, and similarly for Pcl3. doi:10.1371/journal.pgen.1002576.g008

ChIP-seq analysis
The data consisted of short reads from single-end 30-nucleotide sequencing. The short reads were aligned to the reference mouse genome (version MM9) by using bowtie with the parameters -n 2y -best -m 1 [62]. The summary of alignment results is provided in Figure S4A. Paired samples were normalized by equalizing the number of reads in the background: partial sums of order statistics of binned read counts in paired samples were computed, and the ratios of the partial sums were plotted as a function of the quantile cutoff in the sum ( Figure S4B). The point where the curve deviates significantly from linearity was found by fitting a linear regression using half of the data between the 10th and 60th quantile and by finding the first point that deviates from the linear regression line by more than 3 standard deviations of the estimated error. The ratio of partial sums at this critical point was used to scale the read density in one of the samples. The significance of the difference in read counts was assessed by using the Skellam distribution as a null model. We used 400 bp running windows to scan the genome. A p-value cutoff of 10 27 was used to call significant windows, because this cutoff roughly corresponds to an adjusted p-value of 0.05 based on the method of Poisson clumping heuristic for approximating the probability of locally correlated rare events. The relation between the Skellam and adjusted p-values is given by p Skellam = 2log(12p adjusted ) E[C]/L, where L is the genome size, and E[C] the expected correlation length of Skellam p-values across adjacent windows. Overlapping significant windows were concatenated, and we report the new p-values recomputed in the joined windows. More detailed descriptions are included in Text S1.

Gene density analysis
The MM9 mouse genome was partitioned into disjoint 1 Mbp bins. The log fold-change of normalized counts in paired ChIP-seq data with Pcl3 knockdown or scramble was computed in each bin. Gene density was computed by counting the number of RefSeq genes in each 1 Mbp bin. Pearson correlation was computed between gene density and log fold-change across bins. In Figure 4B, we grouped the bins into units of 10; e.g., 10 denotes bins with gene density greater than or equal 0 and less than 10.

Suz12 binding sites in miRNA promoters
To find miRNAs targeted by Suz12, we used a prior annotation of miRNA promoters [64]. The differential Suz12 binding levels shown in Figure 4D-4E were computed within Suz12 ChIP-seq peaks overlapping with the annotated promoters.

Motif analysis
We used a greedy algorithm called LeitMotif which models the background nucleotide distribution with the 9 th order Markov dependence [86]. Motif lengths ranging between 6 and 15 were used for scanning. Top scoring motifs with length 10 and 14 were chosen for further analysis, because other motifs were either contained in or contained those two motifs. Sequences obtained from LeitMotif were aligned to produce the logos in Figure 9A and to generate the PSSM matrices used in Figure 9B-9D. A custom C code was used to scan the Suz12 sites with those PSSM matrices, assuming independent bases in the background. Support vector machine implemented in the e1071 R-package was used to classify the Suz12 binding sites with 10-fold cross validation.

Accession number
The ChIP-seq data are available from Gene Expression Omnibus under access number GSE28325.  tively, in scramble and Pcl3 shRNA expressing EBs. Experiment was performed with 3-6 clones. All immunoblots were performed 2-4 times, and a-tubulin was used as a loading control. All stainings were performed 2-3 times in 2-4 clones or teratoma samples. Expression analysis represents 3-4 experiments assayed in quadruplicate. Error bars indicate standard deviation and asterisks indicate statistical significance of p,0.005. shRNA cells were computed, and their ratios are plotted for the Suz12 FlagM2 ChIP and H3K27me3 ChIP. The theoretical ratio of partial sums is almost linear when two samples are identically distributed, as shown by the red line. The blue vertical bar marks the quantile at which the ratio begins to deviate from linearity, and it effectively separates the background bins from ChIP-enriched bins. The ratio at this quantile was used to scale the Pcl3 ChIP-seq counts. (C) Boxplot comparing Suz12 binding sites to Marson et al. [64]. The most significant Suz12 binding sites identified in our dataset overlap with Marson et al. while the less significant binding sites make up the majority of the remaining sites not identified in Marson et al. ChIP score = 2log p-value for Suz12 ChIP-seq/ Input. (D) Graph of 2log p-values indicates that the most significant Suz12 ChIP-seq binding sites are more likely to decrease following Pcl3 knockdown, while sites unaffected by Pcl3 knockdown are most often less significant Suz12 ChIP-seq binding sites. (E) Boxplot of log fold-changes in ChIP-seq read density within Suz12 binding sites. Suz12 binding and H3K27me3 in Pcl3 knockdown ESCs were decreased on all chromosomes. Chromosome 11 was the most significantly depleted (Pair-wise Wilcoxon rank sum test pvalue,1.3610 215 for Suz12 binding sites; p-value,5.2610 214 for H3K27me3). (F) Sites with decreased Suz12 binding upon Pcl3 knockdown tend to be devoid of E2f1 [63] (Fisher test pvalue = 9.5610 253 ). (G) Suz12 ChIP-qRT-PCR indicating that at some sites Pcl3 depletion does not affect Suz12 binding. (H) ChIP-qRT-PCR for H3K4me3 in scramble and Pcl3 shRNA ESCs. Error bars indicate standard deviation, and asterisks indicate statistical significance of p,0.04. ChIP-qRT-PCRs was performed 2-3 times and assayed in quadruplicate. (TIF) Figure S5 Pcl3 co-localizes with Suz12. (A) Aligning Suz12 and Pcl3 ChIP-seq reads at Pcl3 peak centers shows genome-wide co-localization of Pcl3 with Suz12. (B) To estimate the percentage overlap between Suz12 and Pcl3 binding sites, a sensitivity analysis was performed by varying the Skellam distribution p-value cutoff for calling peaks, ranging between 10 27 to 10 212 . The boxplot shows the percentage of Pcl3 peaks at each p-value cutoff found to be overlapping with Suz12 binding sites that pass the p-value cutoffs 10 27 , 10 28 , 10 29 , 10 210 , 10 211 , and 10 212 . (TIF) Figure S6 Pcl3 misregulates a subset of genes. (A, C) Relative expression of genes that either increased or decreased upon Pcl3 knockdown as measured by qRT-PCR and microarray. Cells were collected following selection for Pcl3 depletion, approximately 3-4 weeks. Microarray and qRT-PCR were performed with 2-6 clones each for control and Pcl3 knockdown ESCs. qRT-PCR analysis represents 2-3 different experiments performed with multiple clones and assayed in quadruplicate. For both graphs, error bars indicate standard deviation and data represents statistical significance of p,0.05. (A) Graph represents expression levels in ESCs that have not been pre-plated to remove differentiated cells. (B) Genes that show depleted Suz12 binding following Pcl3 knockdown and that are misregulated by microarray analysis. Cells in (C) were pre-plated to remove differentiated cells and then expression was assessed. (D) Binding profile of Suz12 and Pcl3 at the Suz12 locus. Turquoise puncta are background and not statistically significant. (TIF) Text S1 Supplemental methods. Detailed description of the microarray and ChIP-sequencing methods as well as a list of the antibodies and primers used. (DOC)