SS18 regulates pluripotent-somatic transition through phase separation

The transition from pluripotent to somatic states marks a critical event in mammalian development, but remains largely unresolved. Here we report the identification of SS18 as a regulator for pluripotent to somatic transition or PST by CRISPR-based whole genome screens. Mechanistically, SS18 forms microscopic condensates in nuclei through a C-terminal intrinsically disordered region (IDR) rich in tyrosine, which, once mutated, no longer form condensates nor rescue SS18−/− defect in PST. Yet, the IDR alone is not sufficient to rescue the defect even though it can form condensates indistinguishable from the wild type protein. We further show that its N-terminal 70aa is required for PST by interacting with the Brg/Brahma-associated factor (BAF) complex, and remains functional even swapped onto unrelated IDRs or even an artificial 24 tyrosine polypeptide. Finally, we show that SS18 mediates BAF assembly through phase separation to regulate PST. These studies suggest that SS18 plays a role in the pluripotent to somatic interface and undergoes liquid-liquid phase separation through a unique tyrosine-based mechanism.

A t the start of mammalian development, a fertilized egg is totipotent, then undergoes successive cell divisions to become a blastocyst from which naive pluripotent stem cells can be derived [1][2][3][4][5] . Then, blastocysts undergo implantation to mark the beginning of fetal development for both the fetus and the placenta. During or immediately after implantation, cells from the inner cell mass in a blastocyst begin to differentiate and eventually endow the three germ layers and germ cells. While this process occurs without any appreciable difficulty in vivo and naive ESCs undergo spontaneous differentiation once cultured under conditions without LIF or other supporting factors 2,6-8 , one may assume that little regulatory mechanism is required for or built into this process.
Given the precision nature during mammalian development, one would argue that each cell division is a fate decision process meticulously regulated. One such example is the pluripotent to somatic transition or PST during early embryonic development. Derived from inner cell mass or ICM of mouse blastocysts 1,5 , mESCs are in vitro copies of ICM, representing the pluripotent state quite faithfully by virtue of their ability to regenerate all fetal tissues when reintroduced to blastocysts even after prolong culture and storage in vitro 9 . In vitro, mESCs are capable of differentiating into cells for all lineages as well 3,10 . Tcf3 (a.k.a. Tcf7l1) has been reported to play a dominant role in pluripotency exit by mediating lineage commitment [11][12][13][14] . On the other hand, Tfe3, a master regulator of lysosomal biogenesis 15 , mediates pluripotency exit by targeting Esrrb through stage-specific subcellular localization 16 . In addition, pluripotency exit can be regulated at various levels, including but not limited to transcription factors 11,12,16 , signal transduction pathway 12,16,17 , epigenetics 18 , lincRNAs 19 . While these sophisticated methods and elegant models used for mESC differentiation in vitro and in vivo, mechanistic insights at the molecular and chromatin level are still lacking at the present.
The ATP-dependent chromatin remodeling complexes BAFs (Brg/Brahma-associated factors), are also mammalian switch/ sucrose nonfermentable (SWI-SNF) initially found in yeast, that consist of 15 subunits encoded by 29 genes, have been shown to play critical roles in self-renewal and pluripotency of ESCs, development 20 and human tumor 21 . Among these subunits, BRG1 is the enzyme core of BAF complexes and has been shown to be required for the maintenance of stem cells [22][23][24] and DPF2, as the canonical BAF specific subunit, is required for maintaining pluripotency and ESCs differentiation 25 . However, it remains to be determined if any of the BAFs are required for pluripotency exit.
CRISPR-Cas9 is an emerging technology that has been adopted for whole genome screen 26 . In this report, we describe a PST system based on inducible over-expression of cJUN in naive mESCs and its adaption for whole genome screening to identify regulators for PST. We identify SS18, a component of BAFs, as a regulator at the pluripotent-somatic interface.

Results
Genome-wide Cas9 screen identifies factors regulating PST. Mouse ESCs or mESCs are naïve pluripotent stem cells capable of self-renewal indefinitely while preserving the capacity to differentiate into any of the three somatic germ layer lineages, thus representing an ideal system to study pluripotent to somatic transition or PST. To identify regulators for PST, we took advantage of our early finding that cJUN behaves as a guardian for the somatic fate in MEFs and is capable of driving mESCs to differentiate rapidly 27 (Supplementary Fig. 1a-c). Using this inducible system, we further show that exposure to cJUN for up to 6 h results in about 50% of the starting mESCs capable of regaining naive pluripotency ( Supplementary Fig. 1d). However, exposure beyond 8 h leads to <15% of the cells regaining the naive state, suggesting that there is a window of sensitivity between 6-8 h that marks the point of no return or a checkpoint between pluripotent and somatic states ( Supplementary Fig. 1d, e). RNAseq analysis further reveals rapid down-and up-regulation of pluripotent genes and somatic ones between 4-8 h, respectively (Supplementary Fig. 1f, g), thus confirming the observed fate transition at the molecular level.
We hypothesize that cJUN accomplishes this rapid PST by cooperating with a cellular machinery critical to cell fate transitions. To gain insight into such machinery, we design a CRISPR based screening system as illustrated (Fig. 1a, Supplementary Fig. 1h). The rationale behind this design is our prediction that any colonies become resistant to cJUN induced PST must harbor a guide RNA targeting a critical component of this machinery (Fig. 1a). We first performed quality control tests to assure proper execution of the screens (Fig. 1b, c, Supplementary Fig. 1i). We then performed the screen as illustrated and identified multiple hits (Fig. 1d). Among the hits, SS18 scored high (Fig. 1d) and can be validated by knockout experiments with the corresponding sgRNA sequences (Fig. 1e) as measured by Nanog expression during cJUN induced PST. Interestingly, p53 appears to play a role in PST, in addition to Ss18 (Fig. 1e). We then performed a second screen and identified similar hits plus Brg1 (also known as Smarca4) ( Supplementary Fig. 1j). Since both SS18 and BRG1 are subunits of the BAF complex 28 , we decided to focus on SS18 for further analysis.
SS18 regulates cJUN-induced PST. Unlike p53 29 , Ss18 has not been reported to be involved in pluripotency nor mESC differentiation. To probe its likely function, we generated Ss18 −/− mESCs as illustrated in Fig. 2a. We obtained KO cell lines (Fig. 2b) that can be confirmed by western blotting for SS18 protein ( Supplementary  Fig. 2a). Ss18 −/− mESCs are morphologically identical to their wild type counterparts (Fig. 2b) and express similar naive genes under naive conditions (Fig. 2c, left panel). When induced to undergo PST, Ss18 −/− mESCs are resistant to cJUN induced PST with substantial cells maintaining Nanog expression at 8 h (Fig. 2c, right panel; Supplementary Fig. 2b). Over-expression of SS18 in Ss18 −/− cells can rescue the defect as measured by Nanog expression (Fig. 2d, Supplementary Fig. 2c) or morphology/GFP expression ( Fig. 2g). At 6 h post induction,~90% Ss18 −/− colonies remain pluripotent morphologically, while only 30% wild type ESCs can do so ( Fig. 2e-g, Supplementary Fig. 2d). Taken together, these results demonstrate that SS18 deficiency confers resistance to cJUN induced PST. We further validated the role of SS18 in other differentiation systems ( Supplementary Fig. 2e), such as spontaneous differentiation, naïve-primed differentiation and neuroectodermal differentiation, and show that the deficiency of SS18 delays all these differentiation processes as measured by Nanog/Esrrb/Klf4 expression ( Supplementary Fig. 2f) or morphology ( Supplementary  Fig. 2g).
SS18 forms condensates in nuclei. To probe the potential mechanisms associated with SS18 mediated PST, we carried out immunostaining experiments and show that endogenous SS18 forms punctate condensates during PST ( Supplementary Fig. 3a).
To help visualize SS18, we tagged it with EGFP and show that SS18-EGFP forms condensates readily when expressed in mESCs. We then analyzed the response of the condensates to 3% 1,6-hexanediol, an aliphatic alcohol which could disrupt weak hydrophobic interactions specifically 30 and show that it can effectively disperse SS18-GFP condensates in 10 s, and more than half of the condensates recover after 1,6-hexanediol removal in 10 s (Supplementary Fig. 3b, c). To characterize the condensates in more detail, we performed FRAP experiments, and show that SS18-EGFP condensates recover rapidly in~10 s after bleaching (Fig. 3a, Supplementary Fig. 3d). We further bleached the SS18-GFP condensates in the edge (Case1) or middle (Case 2) and observed a quick florescence recovery within 6 s in both conditions (Fig. 3b). When we measured the pixel intensity along the white line in case 2, we can show a clear recovery in kymography (Fig. 3c). We then quantified the mobility of SS18 protein by measuring the pixel intensity during fluorescence recovery inside or outside of eight bleached regions within condensates, and show that τ(inside) = 1.841 s, τ(outside) = 1.905 s; Plateau(inside) = 66.18%, Plateau(outside) = 79.91% (Fig. 3d)   of SS18 in the condensates. Moreover, the condensates appear to be able to divide or fuse quite rapidly as well (Fig. 3e). We further performed quantitative analysis of fusion and fission events for the condensates in cells ( Supplementary Fig. 3e) (Fig. 4a). We then constructed a series of mutants for the IDR ( Supplementary Fig. 4a), confirming their expression by western blot ( Supplementary  Fig. 4b), and show that its ability to rescue Ss18 −/− mESCs is dependent on the length of IDR left in the mutants (Supplementary Fig. 4c). Based on the rescue activities, we reason that IDR223-379 may be most critical (pink box in Fig. 4b).
This sub region appears to be enriched with Q, Y, P and G (Fig. 4b). We

Fig. 2 SS18
Regulates cJUN-induced PST. a Construct design for SS18 knockout in mESC. The cleavage sites of Cas9 were designed flanking exon1 and exon2 within Ss18 loci. b SS18KO mESCs are similar to WT ones in morphology. Three biological replicates. c Gene expression profiles by RNA-seq for WT and SS18 KO mESCs. SS18 knockout has little impact on most pluripotent genes but delays the c-JUN induced PST process significantly. d Rescue of the delayed PST as indicated by Nanog expression in SS18 KO mESCs with SS18 in lentiviral vector. Data are mean ± s.d., two-sided, unpaired t test; n = 3 independent experiments, **p < 0.01. e The experimental design to test the window of sensitivity for SS18 response to c-Jun-induced PST. f Results from (e) to show that OCT4-GFP positive colony numbers recovered for 48 h in 2i/LIF medium after c-Jun induction at the indicated time points for wild type and SS18 KO ESCs. Data are mean ± s.d., two-sided, unpaired t test; n = 3 independent experiments. **p < 0.01, ****p < 0.0001. Note no change for 4 h and the greatest change at 6 h for the window of sensitivity. g Representative images for the c-Jun-induced PST in WT, Ss18 −/− and Ss18 −/− overexpressing Ss18, respectively. Three biological replicates.
then constructed four mutants, Q-A, Y-A, P-A and G-A by substituting each of those residues with A ( Fig. 4b, Supplementary Fig. 4d). Among them, Y-A mutant lost its ability to form condensates, while the other three mutants appear to be retaining the ability to form condensates, compared to the WT SS18 or the truncated form SS18-IDR (Fig. 4c, d, Supplementary  Fig. 4e-g). Importantly, we show that this Y-A mutant fails to rescue Ss18 −/− mESCs in PST (Fig. 4c). We further investigated the co-localization of SS18 or mutants with BRG1 by immunostaining and pixel intensity analysis, and show that except Y-A mutant, all the other three mutants can colocalize with BRG1 in the same condensates (Fig. 4d, e). Surprisingly, we find that SS18-IDR (71-418) is sufficient to form condensates, but unable to colocalize with BRG1 (Fig. 4d, e), suggesting that the Nterminal residues are not required for condensates formation, but critical for recruiting BRG1. We further performed IP-MS experiments to investigate whether the SS18 Y-A mutant can pull down the components in BAF complex ( Supplementary  Fig. 4h) and show that Y-A mutant lost the capability to interact with BAF subunits (Fig. 4f). Taking together, these results suggest that SS18 contains a Y-type IDR required for BRG1 interaction and PST.    nor BRD4-IDR, can rescue condensates formation (Fig. 5a). Consistently, poly-Y, TAF15-and FUS-PLDs can functionally replace SS18-IDR in rescuing SS18 −/− defect during cJUN-based PST (Fig. 5b). We further performed IP-MS experiments to investigate whether these SS18 truncation/fusions can pull down the subunits in BAF complex ( Supplementary Fig. 5b), and show by pairwise comparison scatterplots that SS18 N-terminal, N-polyY and N-FUS(IDR) can pull down the subunits of BAF complex while SS18(IDR) could not (Fig. 5c). Notably, we show that N-BRD4 could pull down the subunits in BAF complex, but still fails to mediate PST (Fig. 5b, c), perhaps due to the type of condensates formed (Fig. 5a). We summarize the results into a table in Fig. 5d and show that SS18 mediates PST by forming BAF complex through its N-70 domain and condensates with a Cterminal IDR (Fig. 5e). These results further suggest that IDRs are interchangeable among the Y-type ones and can be designed artificially like the poly-Y.
To investigate condensates associated with SS18, we stained SS18KO mESCs expressing SS18-GFP with antibodies against endogenous BAF specific subunit DPF2 or PBAF specific subunit BRD7 and show that DPF2, not BRD7, is enriched with SS18 ( Fig. 6a, b, Supplementary Fig. 6b), which is also consistent with endogenous SS18 (Supplementary Fig. 6c, d, e). These results demonstrate that SS18 is associated exclusively with BAF, not PBAF, in mESCs.
To test the role of SS18 in BAF assembly, we immunoprecipitated the common subunit BRG1 in SS18 WT and SS18 KO mESCs and then analyzed the complexes by IP-MS to detect the relative abundance of each components. We found SS18 deficiency decreases the assembly of canonical BAF and ncBAF, but promotes that of PBAF (Fig. 6c). In addition, the full length not N terminal of SS18 can rescue BAF assembly in SS18 KO mESCs ( Supplementary Fig. 7a, b). Given the two unbiased screens we performed that identified both SS18 and Brg1 but no PBAF specific subunits, we speculate that PST relies on primarily BAF or/and ncBAF, not PBAF. We indeed can show that BAF ncBAF and PBAF serve different roles in PST by shRNA knockdown of BAF/ncBAF/PBAF specific subunits and shared ones (Fig. 6d, Supplementary Fig. 7c, Supplementary Data 2). The depletion of BAF specific components, DPF1/2, and ncBAF specific ones can delay this transition, while depleting PBAF specific subunit BRD7 promotes PST (Fig. 6d). Unlike Ss18, Brg1 knockdown in mESCs leads to defects in self-renewal and survival morphologically ( Supplementary Fig. 7d), consistent with a previous report [22][23][24] . We further performed RNA-Seq experiments to investigate the downstream targets of Brg1 in mESC, and show by GO analysis that the 572 genes affected by Brg1 knockdown are involved in cell-substrate adhesion, positive regulation of cell migration and epithelial cell development ( Supplementary Fig. 7e). However, BRD9 inhibitor BI-9564 has no effect in the process (Supplementary Fig. 7f). We knocked down the specific components of ncBAF and PBAF in SS18 KO mESCs and found the deficiency of ncBAF will delay PST further, on the contrary, the knockdown of PBAF will rescue that ( Supplementary Fig. 7g), which suggests that BAF and ncBAF play similar/redundant roles against PBAF in cJUN induced PST.
To further investigate the mechanism underlying the different effects between BAF and PBAF within PST, we performed ChIPseq of H3K27ac, cJUN, SS18, and PBAF specific subunits PBRM1 and BRD7 at different time points during PST ( Supplementary  Fig. 7h). All the H3K27ac occupancy sites ± 5 kb of 0 h and 8 h are divided into three categories (8 h down, permanent and 8 h up) according to the dynamic changes of intensity by fold change >2. Both BAF and PBAF complexes relocate to the same sites occupied by cJUN. We can estimate the relocation rate from the ratio of average intensity of up sites and down sites from the graph. Contrary to PBAF subunits, SS18 has higher relocation rate at 4 h and 8 h during PST ( Supplementary Fig. 7h), suggesting that BAF containing SS18, compared with PBAF, can facilitate cell fate transition more efficiently. Taken together, these data suggest that SS18 promotes BAF but inhibits PBAF assembly to mediate PST (Fig. 6e).

Discussion
Broadly speaking, the mammalian cell fate continuum can be divided into two categories: pluripotent and somatic. Pluripotent cells have the capacity to give rise to all somatic cells 9 . On the other hand, somatic cells can be reprogrammed to a pluripotent state [38][39][40] . The inter-conversion between somatic and pluripotent cells through differentiation and reprogramming suggests that there may be an interface between somatic and pluripotent states. We have proposed previously that cJUN serves as a guardian of somatic fate as opposed to factors such as Oct4 that specify the pluripotent fate 27 . In this report, we further define a molecular link at the pluripotent-somatic interface. Mechanistically, we propose that cJUN drives a rapid transition through this interface with the help of SS18, a well-known member of the BAF complex. We provided evidence that SS18 is a bipartite protein with an Nterminal 70aa domain specific for BAF interaction and a Cterminal IDR for forming condensate through a unique tyrosinebased mechanism.
It is of interest to note that the IDR of SS18 can be replaced by similar IDRs from FUS and TAF15 33 , but not BRD4 and MED1 32 (Fig. 5b). We noted that the IDRs of FUS and TAF15, like that of SS18, are enriched with tyrosine residues, while those of BRD4 and MED1 are not (Supplementary Fig. 5a). Based on these observations, we further reasoned that perhaps a polymer of tyrosine residues can also function as an IDR to form condensates. Indeed, when we replaced the IDR of SS18 with a 24tyrosine residue polypeptide, similar condensates can be also observed (Fig. 5A, N-Poly Y). This construct can also rescue full length SS18 function in PST assays (Fig. 5b). Based on these results, we would like to propose that the 24-residue polymer can serve as a model tool to further analyze liquid-liquid phase separation or condensate formation in vitro and in vivo. The simple composition as well as well-behaved pattern of condensates may help further illuminate principles and rules associated with phase separation. More generally, it may be possible to design additional polymers that behave similarly to other classes of IDRs such as those from BRD4 or MED1.
Chromosome translocation t(X;18) (p11.2; q11.2) has been suggested as a cause of synovial sarcoma 28 . Our dissection of its structure into two main domains, i.e., the N-terminal BAF interacting domain and the C-terminal IDR, may further encourage innovative strategies to design targeted therapies.
Gene knockout in mESCs. Knockout mouse ESCs lines were generated using the CRISPR/Cas9 system 41,42 . Guide RNAs were designed using crispr.mit.edu. and cloned downstream of the human U6 promoter of pX330 (Addgene). For SS18 KO in mouse ESCs, two pairs of Guide RNAs were designed for each gene to delete critical exons. For Ss18, the gene locus of exon1 to exon2 were deleted. Mouse ESCs clones carrying deletions in both alleles were identified by genotyping and western blot.
Immunoblotting. Cells were collected and lysed in lysis buffer supplemented with protease inhibitor cocktail (Roche) on ice for 15 min, and then cells were boiled in 100°C for 10 min. After centrifugation, the cell supernatants were subjected to SDS-PAGE and incubated with corresponding primary antibody and secondary antibodies. The following antibodies were used in the project: anti-SS18 (CST no. Immunofluorescence. Cells growing on coverslips were washed three times with PBS, then fixed with 4% PFA for 30 min, and subsequently penetrated and blocked with 0.1% Triton X-100 and 3% BSA for 30 min at room temperature. Then, the cells were incubated with primary antibody for half hour. After five washes in PBS, 1 h of incubation in secondary antibodies, cells were then incubated in DAPI (Sigma-Aldrich D9542) for 2 min. Then, the coverslips were mounted on the slides for observation on the confocal microscope (Zeiss LSM 900 with Airyscan detector). The following antibodies were used in this project: anti-SS18, (CST no. 21792.  The graph on the right shows pixel intensity across five condensates of SS18-EGFP (green) and DPF2 (red). Scale bars, 3 μm. SS18-EGFP, laser intensity 2.5%, detector gain 500 V. DPF2 immunofluorescence, laser intensity 3.5%, detector gain 750 V. Both SS18-EGFP and DPF2 immunofluorescence were processed by background subtraction. Three independent experiments. b Representative images of the immunofluorescence for endogenous PBAF specific subunit BRD7 in SS18KO mESCs with exogenous SS18-EGFP. The graph on the right shows pixel intensity across five condensates of SS18-EGFP and BRD7 (red). Note no co-localization between SS18 and BRD7. Scale bars, 3 μm. SS18-EGFP, laser intensity 2.5%, detector gain 500 V. BRD7 immunofluorescence, laser intensity 3.5%, detector gain 750 V. Both SS18-EGFP and BRD7 immunofluorescence were processed by background subtraction. Three independent experiments. c The histogram shows the relative precipitated abundance of indicated proteins to BRG1 as quantified by normalized LFQ intensity from IP-MS experiments by anti-BRG1 antibody performed on WT and SS18KO mESCs. IP-MS experiments were performed in triplicates and a two-sample t test was applied. Data are mean ± s.d., two-sided, unpaired t test; n = 3 independent experiments, *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. d The expression levels of Nanog by Jun-FlagTetON mESCs at 8 h during PST with knocking down of indicated components of BAF ncBAF or PBAF. G1 + G1l, knockdown Gltscr1 and Gltscr1l jointly. Data are mean ± s.d., two -sided, unpaired t test; n = 3 independent experiments, **p < 0.01, ***p < 0.001, ****p < 0.0001. e A working model for SS18-mediated PST. SS18 functions by interacting with BAF or ncBAF subunits through its N-terminal 70aa, and forming condensates through a tyrosine-rich IDR. Meanwhile, the full length of SS18, co-localizing with BAF-specific component and isolating PBAF-specific one, can promote the BAF complex assembly and impair that of PBAF. In WT cells, BAF ncBAF and PBAF reach a balance. In SS18KO cells, the balance shifts towards PBAF to cause a delay in PST.
1:1000). anti-BRG1, (CST no. 52251 1:1000). anti-DPF2, (Proteintech, no. 12111-1-AP-50 μl,1:100). anti-BRD7, (Proteintech, no.51009-2-AP-50 μl, 1:100). Genome-wide CRISPR-Cas9 screens. The lentiviral gRNA plasmid library and CRISPR/Cas9 expressing plasmid for genome-wide CRISPR-Cas9 screen were obtained from Addgene (#1000000053). Amplification of the library and preparation of the lentivirus were performed following the protocol provided by Addgene with electrocompetent E. coli (TaKaRa). HygroR-cJunTetONESCs stably expressing Cas9 were seeded in the 15-cm dishes (about 3 × 10 6 cells per dish) and a total of 2.5 × 10 7 cells were infected with the gRNA lentivirus library to keep the diversity of library. Meanwhile, the MOI of infection was <0.3. Puromycin (1 µg/ml GIBCO A1113802) was used to eliminate non-infected cells after 48 h. The infected cells were collected as the sample of input 1 and the cells of further three passages were collected as the sample of input 2. Next, 2 µg/ml doxycycline was used to trigger the expression of HygroR-cJun fusion protein and 10 µg/ml hygromycin was added in the culture medium to eliminate the cells without response to doxycycline. Dox was added only for the first 24 h at passage1 to drive PST but withdrew in the following three passages. The GFP positive cells were analyzed by FACS when passed in 24 h after Dox withdrawal. Almost all the cells would have transitioned from pluripotency to somatic state. The cells resisting PST or apoptosis were enriched for three passages in mESC medium, while the rest ones were eliminated gradually. Finally, the enriched cells were collected as the sample of output group.
Illumina sequencing of gRNAs and statistical analysis. The genomic DNA of the three samples (input 1, 2 and output group) were extracted and the Illumina sequencing of gRNAs were conducted as described previously 26 . The counts of unique sgRNA for a given sample were normalized as follows: normalized sgRNA counts ¼ countsofpersgRNA thesumofcountsofallsgRNA Enrichment and depletion of guides and genes were analysed using MAGeCK statistical package5 43 by comparing read counts between input 2 and input 1 or output group and input 2.
FACS analysis. For OCT4-GFP analysis in the process of screen, cells were suspended in DPBS supplemented with 2% FBS (FACS buffer) for direct detection. DAPI was used to exclude the nonviable cells. The cells were then analyzed on Accuri C6 flow cytometer (BD Biosciences). The GFP fluorescence intensity was detected in the FITC channel. Data analysis was done using FlowJo7.6.1. software (LLC).
ChIP-Seq. Briefly, cells were fixed with 1% formaldehyde for 10 min and then followed by the reaction with 0.125 M glycine. Cells were then lysed in ChIP-buffer A (50 mM HEPES-KOH, 140 mM NaCl, 1 mM EDTA (pH 8.0), 10% glycerol, 0.5% NP-40, 0.25% Triton X-100, 50 mM Tris-HCl (pH 8.0), and protease inhibitor cocktail) for 10 min at 4°C. Samples were centrifuged at 1400 g for 5 min at 4°C. Pellets were lysed in ChIP-buffer B (1% SDS, 50 mM Tris-HCl (pH 8.0), 10 mM EDTA and protease inhibitor cocktail) for 5 min at 4°C. The DNA was fragmented to 100-500 bp by sonication, and then centrifuged at 12,000 g for 2 min. The supernatant was diluted with ChIP IP buffer (0.01% SDS, 1% Triton X-100, 2 mM EDTA, 50 mM Tris-HCl (pH 8.0), 150 mM NaCl and protease inhibitor cocktail). Immunoprecipitation was performed with 1 µg anti-H3K27ac antibody (Abcam, ab4729) coupled to Dynabeads with proteinA/G overnight at 4°C. Beads were washed, eluted and reverse crosslinked. DNA was extracted by phenol/chloroform for sequencing. The ChIP DNA library was constructed following the Illumina ChIP-seq library generation protocol. Briefly, 5 ng ChIP DNA was blunt-ended, and then a dA tail was added. Illumina genomic adapters with index sequences were ligated to the DNA. The adapter-ligated DNA was amplified by PCR for 18 cycles, and the resulting DNA libraries were quantified and tested by qPCR with positive primers to assess the quality of the library. The ChIP DNA library for NextSeq 500 sequencing was constructed with VAHTS Turbo DNA Library Prep Kit for Illumina (Vazyme) according to manufacturer's instructions. The ChIP-seq for other factors in mESCs was constructed with NovoNGS CUT&Tag 2.0 High-Sensitivity Kit for Illumina (novoprotein N259-YH01) according to manufacturer's instructions. ChIP-seq data are mapped to the 10 mm mouse genome assembly using bowtie2, version 2.4.1.
SS18 endogenous immunoprecipitation and on-bead digestion. Whole cell extracts of mES cells with cJun overexpression were prepared using lysis buffer (50 mM Tris pH 8.0, 150 mM NaCl, 10% Glycerol, 0.5% NP40) with fresh added 1x Complete Protease inhibitors (Sigma, 1187358001). Cells were incubated for 2 h at 4°C on a rotation wheel. Soluble cell lysates were collected after maximum speed centrifugation at 4°C for 15 min. 1 mg of cell lysates were incubated with either SS18 antibody or matched IgG overnight at 4°C on a rotation wheel. Combined Protein A/G magnetic beads (Bio-rad, 1614833) were added for another 1.5 h. Beads were then washed three times with wash cell lysis buffer and one time with PBS. After completely removal of PBS, immunoprecipitated proteins were digested using on-bead digestion protocol as described before 17. Briefly, beads were incubated with 100 mL of elution buffer (2 M urea, 10 mM DTT and 100 mM Tris pH 8.5) for 20 min. Afterwards, iodoacetamide (Sigma, I1149) was added to a final concentration of 50 mM for a 10 min in the dark, following with 250 ng of trypsin (Promega, V5280) partially digestion for 2 hr. After incubation, the supernatant was collected in a separate tube. The beads were then incubated with 100 mL of elution buffer for another 5 min, and the supernatant was collected in the same tube. All these steps were performed at RT in a thermo shaker at 500 g. Combined elutes were digested with 100 ng of trypsin overnight at RT. Finally, tryptic peptides were acidified to pH < 2 by adding 10 mL of 10% TFA (Sigma, 1002641000) and desalted using C18 Stage tips (Sigma, 66883-U) prior to MS analyses. Each experiment was performed in technical triplicate.
Mass spectrometry analysis. Tryptic peptides were separated by AcclaimTM PepMapTM 100 C18 column (Thermo, 164941) using a 140 min of total data collection (100 min of 2-22%, 20 min 22-28% and 12 min of 28-36% gradient of B buffer (80% acetonitrile and 0.1% formic acid in H2O) for peptide separation, following with two steps washes: 2 min of 36-100% and 6 min of 100% B buffer) with an Easy-nLC 1200 connected online to a Fusion Lumos mass spectrometer (Thermo). Scans were collected in data-dependent top-speed mode with dynamic exclusion at 90 s. Raw data were analyzed using MaxQuant version 1.6.0.1 search against Mouse Fasta database, with label free quantification and match between runs functions enabled. The output protein group was analyzed and visualized using DEP package as described before 44 .
Statistical information. Data are presented as mean ± s.d. as indicated in the figure legends. Unpaired two-tailed student t test, The P value was calculated with the Prism 6 software. A P < 0.05 was considered as statistically, *P < 0.05, **P < 0.01, ***P < 0.001. No statistical method was used to predetermine sample size. The experiments were not randomized. The investigators were not blinded to allocation during experiment and outcome assessment.
Reporting summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability
The RNA-Seq, ChIP-seq data have been deposited in the Gene Expression Omnibus database under the accession code GSE135451. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (http://proteomecentral. proteomexchange.org) via the iProX partner repository 45 with the dataset identifier PXD026208. Source data are provided with this paper. The source data underlying Figs. 1e, 2d, f, 3c, d, 4c, 5b, 6c, d and Supplementary Figs. 1d, g, 2b, d, f, 3c-e, 4c, e-g, 6e, 7b, c, f, g are provided as a Source Data file. Supplementary Data 1 and 2 are primers used in this study. The Cas9 screen data and IP-MS data are provided in Supplementary Data 3 and 4, respectively. The authors declare that all data supporting the findings of this study are available within the article and its supplementary information files or from the corresponding author upon reasonable request. Source data are provided with this paper.