SUMOylated PRC1 controls histone H3.3 deposition and genome integrity of embryonic heterochromatin

Abstract Chromatin integrity is essential for cellular homeostasis. Polycomb group proteins modulate chromatin states and transcriptionally repress developmental genes to maintain cell identity. They also repress repetitive sequences such as major satellites and constitute an alternative state of pericentromeric constitutive heterochromatin at paternal chromosomes (pat‐PCH) in mouse pre‐implantation embryos. Remarkably, pat‐PCH contains the histone H3.3 variant, which is absent from canonical PCH at maternal chromosomes, which is marked by histone H3 lysine 9 trimethylation (H3K9me3), HP1, and ATRX proteins. Here, we show that SUMO2‐modified CBX2‐containing Polycomb Repressive Complex 1 (PRC1) recruits the H3.3‐specific chaperone DAXX to pat‐PCH, enabling H3.3 incorporation at these loci. Deficiency of Daxx or PRC1 components Ring1 and Rnf2 abrogates H3.3 incorporation, induces chromatin decompaction and breakage at PCH of exclusively paternal chromosomes, and causes their mis‐segregation. Complementation assays show that DAXX‐mediated H3.3 deposition is required for chromosome stability in early embryos. DAXX also regulates repression of PRC1 target genes during oogenesis and early embryogenesis. The study identifies a novel critical role for Polycomb in ensuring heterochromatin integrity and chromosome stability in mouse early development.


Introduction
Major parts of metazoan genomes are composed of repetitive sequences such as endogenous retroviral sequences (ERVs) and tandem repeat sequences flanking and underlying centromeres and at telomeres. These sequences are packaged in constitutive heterochromatin characterized by H3 lysine 9 di-/trimethylation (H3K9me2/3) and occupancy of HP1 proteins and are largely transcriptionally repressed (Allshire et al, 1995;Peters et al, 2001;Ebert et al, 2004). While such chromatin states are effectively perpetuated during DNA replication (Alabert et al, 2017;Mendiratta et al, 2019), they also undergo nucleosome turnover outside of S phase and incorporate nucleosomes containing the histone variant H3.3. H3.3 is expressed throughout the cell cycle and is mostly deposited at transcription start sites (TSS) and along gene bodies of transcriptionally active genes, and at enhancers by the HIRA/CABIN1/UBN1 complex (Goldberg et al, 2010;Szenker et al, 2011;Martire et al, 2019;Mendiratta et al, 2019). In contrast, the chaperone DAXX drives H3.3 incorporation at and facilitates transcription of pericentromeric and telomeric repeat regions (Drane et al, 2010;Lewis et al, 2010;Gauchier et al, 2019). DAXX is required for the structural organization of PCH (Rapkin et al, 2015), and its localization at PCH is mediated via its interaction with ATRX (alpha thalassemia/mental retardation syndrome X-linked), a chromatin remodeler recruited via binding to H3K9me3 and HP1 Eustermann et al, 2011;Iwase et al, 2011;Hoelper et al, 2017). At ERVs, H3.3 deposition is not required for transcriptional silencing, yet it increases DAXX occupancy levels by stabilizing the protein (Hoelper et al, 2017).
In embryonic stem cells (ESCs) and other somatic cells, conjugation with the small ubiquitin-like modification (SUMO) is a prominent component of heterochromatic components at PCH, telomeres, and ERVs, with DAXX, ATRX, KAP1, and SETDB1 being SUMOylated (Hendriks et al, 2014). Moreover, the H3K9me3 histone methyltransferase (HMT) SUV39H1 functions as an E3 SUMO1 ligase for HP1 proteins (Maison et al, 2016), thereby promoting recruitment of HP1a to canonical H3K9me3/HP1 PCH in mouse fibroblasts (Maison et al, 2011). Importantly, deficiency of the SUMO E2 conjugating enzyme UBC9 induces major chromosome condensation and segregation defects suggesting a major role of SUMOylation in maintaining the heterochromatic state of H3K9me3/HP1 at PCH and genome integrity (Nacerddine et al, 2005). Indeed, SUMO2 was recently shown to be required in human cells for the association of DAXX and for H3.3 deposition at centromeres, presumably via SUMOylation of CENP-B. CENP-B and DAXX deficiency impaired the heterochromatin state at and function of centromeres (Morozov et al, 2017).
Besides at actively expressed genes, nucleosomes also undergo turnover at TSS of genes repressed by Polycomb group proteins. PcG proteins bind and robustly modify the chromatin state of their target loci in order to perpetuate a specific transcriptional program over cellular divisions (Beisel & Paro, 2011). They assemble in two major types of complexes, Polycomb Repressive Complex 1 and complex 2 (PRC1 and PRC2), which exist in multiple flavors harboring different subunits and enzymatic activities. The core subunits of PRC1, RING1 and RNF2, are responsible for catalyzing ubiquitination of lysine 119 of histone H2A (H2AK119ub1), while the EZH1/2 subunits of PRC2 mediate H3K27me3 (Cao et al, 2002;de Napoles et al, 2004). ESC models with synthetic targets showed that canonical PRC1 (cPRC1) complexes containing CBX proteins mediate heritable gene silencing by sustaining a positive feedback mechanism for recruitment of PRC2 and cPRC1 complexes through H3K27me3 and H2AK119u1 (Moussa et al, 2019). In contrast, variant PRC1 complexes (vPRC1) do not contain CBX subunits and do not rely on preexisting H3K27me3 for chromatin targeting Cooper et al, 2014). Instead, vPRC1 complexes mediate repression of targets genes in ESCs by catalyzing H2AK119u1 (Fursova et al, 2019;Scelfo et al, 2019) and stimulate binding of PRC2 to their targets (Kalb et al, 2014;Cooper et al, 2016). Intriguingly, H3.3 and HIRA are required to sustain H3K27me3 levels at and repression of some developmental transcription factor genes in ESCs, possibly by modulating PRC2 recruitment to chromatin (Banaszynski et al, 2013). Moreover, during heart development, the repressive activity of cPRC1 is controlled by CBX4, which functions as a SUMO E3 ligase for several chromatin proteins and itself (Kang et al, 2010). Auto-SUMOylation of CBX4 is counteracted by the SUMO protease SENP2. Loss of Senp2 induces increased recruitment of cPCR1 to PcG target genes and their repression (Kang et al, 2010).
Genetic and chemical perturbation studies in ESCs have revealed plasticity in heterochromatin states at PCH, telomeres, and ERVs between H3K9me3/HP1, DNA methylation and Polycomb repressive pathways (Peters et al, 2003;Saksouk et al, 2014;Walter et al, 2016;Gauchier et al, 2019). In development, natural plasticity in the formation of repressive heterochromatic states has been observed at PCH in mouse pre-implantation embryos (Santos et al, 2005;Puschendorf et al, 2008;Albert & Peters, 2009). PCH on maternal chromosomes (mat-PCH) is defined by Suv39h2-dependent and maternally inherited H3K9me3, and by HP1b and ATRX (Santos et al, 2005;Puschendorf et al, 2008;De La Fuente et al, 2015). In contrast, chromatin of the paternal genome becomes extensively reprogrammed in early embryos. Following the exchange of most nucleosomes by highly basic nuclear proteins during sperm formation (Gill et al, 2012), the nucleosomal configuration is re-established along the paternal genome shortly after fertilization. Such reconstitution occurs prior to replication and is mediated by HIRAdriven H3.3 incorporation (Loppin et al, 2005;Lin et al, 2014). It is, however, unknown whether HIRA controls H3.3 deposition at PCH (Lin et al, 2014). We demonstrated previously that proteins of the cPRC1 and PRC2, and their associated histone modifications, constitute an alternative repressive state at PCH of paternal chromosomes, in response to the absence of paternally inherited H3K9me3 (Puschendorf et al, 2008;Tardat et al, 2015). Remarkably, the H3.3 variant is incorporated at PCH after sperm decondensation and remains further enriched during the S/G2 phases of zygotes at paternal, but not maternal, PCH in a replication-independent manner (Santenard et al, 2010;Akiyama et al, 2011), suggesting the existence of an alternative mechanism for H3.3 deposition, independent of the classical repressive marking such as H3K9me3 and HP1, normally targeting ATRX and thereafter DAXX.
In this study, we unravel the molecular mechanism underlying specific H3.3 deposition at pat-PCH in mouse zygotes. We show that asymmetric loading of DAXX at pat-PCH depends on cPRC1 containing CBX2. Germline deficiency for PRC1 subunits Ring1 and Rnf2 results in loss of binding of DAXX and H3.3 occupancy at pat-PCH. The two SUMO-interacting motifs (SIMs) of DAXX are required for its association with pat-PCH implying a role for SUMOylation in DAXX chromatin targeting to these loci. Accordingly, mutation of specific residues in CBX2, which impair its SUMOylation, prevent DAXX targeting to PCH. Finally, we demonstrate that loss of H3.3 at pat-PCH upon Daxx knockout induces chromatin decompaction and breakage at PCH of exclusively paternal chromosomes and causes their mis-segregation. We show that H3.3 deposition by DAXX is required for chromosome stability in early embryos. Thus, we identify a novel pathway and role for SUMOylation and Polycomb in ensuring chromatin integrity. Genome-wide transcriptional analysis shows that Daxx regulates repression of PRC1 target genes in oocytes and 2-cell embryos. Our data suggest a regulatory function of the novel CBX2/cPRC1 ? SUMO2 ? DAXX ? H3.3 pathway in PRC1-mediated gene silencing during mouse development.

Results
The histone variant H3.3 is incorporated into pat-PCH prior to the first round of DNA replication The paternal genome undergoes extensive chromatin remodeling shortly after fertilization, with the replacement of sperm-born protamines by maternally provided histones. The remodeling process occurs many hours before the first round of replication arguing for de novo nucleosome deposition onto the paternal DNA template. To monitor the timing of incorporation of histone proteins at pat-PCH in mouse zygotes, we microinjected mRNAs encoding for EGFP-tagged H3.2 and mCherry-tagged H3.3 proteins into metaphase II (M-II) oocytes prior to their activation by intracytoplasmic sperm injection (ICSI). We monitored the localization of the tagged histones by fluorescence spinning-disk live microscopy in fertilized embryos (Fig EV1A; Movies EV1-EV3). As reported previously (Akiyama et al, 2011), canonical H3.2 became incorporated in both pronuclei at the time of replication ( Fig 1A). The deposition of the H3.3 variant into the maternal genome started with similar kinetics. In contrast, H3.3 was incorporated into the paternal genome immediately after fertilization (Lin et al, 2014), including the DAPI-brightly stained PCH domain of decondensing sperm heads (Fig 1A). H3.3 was clearly detectable at pat-PCH surrounding pre-nucleolar bodies in G2-stage zygotes (Santenard et al, 2010). Three variant residues in the globular domain of histones H3.3 versus H3.1/H3.2 confer specificity for associations with DAXX and HIRA/UBN1 versus CAF1 (Ahmad & Henikoff, 2002;Goldberg et al, 2010;Elsasser et al, 2012;Ricketts et al, 2015). Importantly, expression of a H3.3 mutant protein, carrying either a triple amino acid substitution mimicking H3.1/H3.2 or only the single G90M substitution, abrogated H3.3 incorporation into PCH of the decondensing sperm in zygotes (Fig EV1B and C). These data underscore the specificity of the assay. While Hira is required for H3.3 deposition in the decondensing sperm (Lin et al, 2014), it remains open whether H3.3 deposition at PCH is controlled by HIRA or another H3.3 chaperone, such as DEK or DAXX (Drane et al, 2010;Ivanauskiene et al, 2014).

The H3.3 chaperone DAXX localizes to pat-PCH in mouse zygotes
In somatic cells, the H3.3 chaperone DAXX and the chromatin remodeler ATRX drive incorporation of H3.3 at various heterochromatin regions including PCH (Drane et al, 2010;Goldberg et al, 2010;Lewis et al, 2010). During oogenesis and pre-implantation embryogenesis, Daxx and Atrx but also other H3.3 chaperones like Hira and Dek are abundantly expressed (Fig EV1D, and Park et al, 2013). Exogenous expression of DAXX was shown to mediate H3.3 incorporation at pat-PCH in zygotes deficient for Stella (Arakawa et al, 2015). To evaluate a possible contribution of endogenous DAXX and ATRX to H3.3 deposition at pat-PCH, we investigated their localization in zygotes by immunofluorescence (IF). Following fertilization, both proteins co-localized at PCH regions of decondensing maternal chromosomes (Fig 1B), marked by H3K9me3 and HP1b (Puschendorf et al, 2008). Such localization argues for targeting of ATRX to H3K9me3-modified nucleosomes via its ATRX-DNMT3-DNMT3L (ADD) domain and to HP1 proteins via its "Leu-X-Val-X-Leu" motif (Eustermann et al, 2011;Iwase et al, 2011). Indeed, ATRX and particularly DAXX localization was reduced at PCH in fully grown germinal vesicle (GV) oocytes deficient for Cbx1 (Hp1b) (Fig EV1E and F) (Tardat et al, 2015) and fully abrogated in GV oocytes deficient for the Suv39h2 HMT lacking both H3K9me3 and HP1b at PCH (Fig EV1E and F) (Peters et al, 2001). Following cell cycle progression, we observed a rapid reduction in DAXX levels at mat-PCH despite the persistence of ATRX (Fig 1B). In contrast, DAXX already localized at pat-PCH in decondensing sperm, while ATRX was not detectable there. Live-imaging analysis of decondensing sperm confirmed DAXX's early association to pat-PCH ( Fig EV1G). Shortly thereafter, prior to replication, ATRX started to co-localize with DAXX at pat-PCH (Figs 1B and EV1H). Using a custom developed method to quantify IF signal specifically at PCH domains in zygotes (Appendix Fig S1), we confirmed that DAXX is strongly enriched at pat-PCH compared to mat-PCH despite ATRX's presence maternally ( Fig 1C). Interestingly, for unidentified reasons, DAXX and ATRX co-localize to pat-PCH in a dotted "pearls on a string" pattern, while the signals are more continuous along the nucleolar precursor bodies of mat-PCH ( Fig 1B).
To determine recruitment dependencies between ATRX and DAXX, we efficiently depleted Daxx by siRNA injection (Fig. 1D) and investigated ATRX localization in late-stage zygotes. While the ATRX signal at euchromatin and mat-PCH was unaffected, ATRX was specifically lost from pat-PCH (Fig 1D), demonstrating that ATRX association with pat-PCH requires DAXX. In contrast, maternal deficiency of the ATRX protein, resulting from conditional deletion of the Atrx gene during oogenesis, caused complete loss of DAXX protein from mat-PCH only, while DAXX levels at pat-PCH were unaffected ( Fig 1E). Thus, these genetic experiments demonstrate unambiguously that DAXX localization at pat-PCH is upstream of ATRX, while it is downstream at mat-PCH.

Canonical PRC1 associates genetically and biochemically with DAXX
We next investigated whether PRC2 or cPRC1 controls DAXX localization to pat-PCH. A previous study implicated lysine 27 of H3.3 (H3.3K27) in formation and function of PCH in early embryos (Santenard et al, 2010). PRC2 and H3K27me3, though, do not localize to pat-PCH before replication (Puschendorf et al, 2008). Moreover, maternal deficiency for Ezh1 and/or the Ezh2 paralog did not affect DAXX localization at pat-PCH (Fig 2A) excluding PRC2 and H3K27me3 as key determinants of DAXX recruitment and function in mouse zygotes.
In mouse zygotes, CBX2 controls cPRC1 targeting to pat-PCH (Tardat et al, 2015). To investigate the biochemical association of DAXX with CBX2-PRC1 components, we performed co-immunoprecipitation (co-IP) experiments with extracts prepared from HEK293 cells expressing DAXX MycHis and CBX2 3×Flag proteins. We found that CBX2 3×Flag co-immunoprecipitated with DAXX MycHis and endogenous RNF2, even in the presence of benzonase, excluding a role for RNA or chromatin in these interactions ( Fig 2D). CBX2 3×Flag also immunoprecipitated H3K27me3, likely reflecting the interaction with the chromodomain of CBX2 as demonstrated previously (Tardat et al, 2015). Together, these results reveal a novel association of cPRC1 with DAXX that controls targeting of DAXX to chromatin.

The cPRC1 subunit CBX2 targets DAXX to chromatin
To determine whether CBX2, as part of cPRC1, can target DAXX to other loci in the genome, we co-expressed DAXX MycHis with human full-length CBX2 (CBX2.1), N-terminally tagged with a GAL4 DNAbinding domain in HEK293 5xUAS cells which carry a single copy of the 5xUAS-luciferase reporter gene (Stielow et al, 2008
To investigate which domains of CBX2 are required for DAXX recruitment, we co-expressed DAXX MycHis with CBX2.2, a shorter, natural splice isoform of human CBX2 (Völkel et al, 2012), that contains the N-terminal chromo domain and the AT hook of CBX2.1, yet harbors a different C-terminus. CBX2.2 lacks most of the middle region and C-terminus of CBX2.1 including the Pc-Box, a sequence known to mediate interactions with other cPRC1 components ( Fig 3B) (Schoorlemmer et al, 1997). Though both GAL4DBD CBX2.1 and GAL4DBD CBX2.2 proteins were similarly expressed and showed robust enrichment at the 5x-UAS site, the truncated CBX2.2 failed to recruit RNF2 or DAXX, or to mediate H3.3 incorporation (Fig 3C and D). These data argue that DAXX ◀ Figure 1. Parent-of-origin-specific targeting pathways for DAXX and ATRX to PCH in mouse zygotes.
A Still images of time-lapse imaging of H3.3-mCherry and H3.2-EGFP proteins in wild-type zygotes (n = 17). B IF of DAXX and ATRX in zygotes at different pronuclear (PN) stages (PN0, n = 8; PN2, n = 6; PN5, n = 33). C Boxplot displaying signal intensity of DAXX at pat-and mat-PCH in mouse zygotes at PN5 stage, relative to mean signal at pat-PCH. In boxplots, the enter lines show the medians; box limits indicate the 25 th and 75 th percentiles as determined by R software; and whiskers extend 1.5 times the interquartile range from the 25 th and 75 th percentiles. Datapoints are indicated by circles. Experiments were replicated three times. ***P < 0.001, n = 28, paired t-test. D IF of DAXX and ATRX in PN5 zygotes upon siRNA-mediated knockdown of Daxx (n = 14) or with control siRNAs (n = 8). E IF of DAXX and ATRX in control (n = 15) and Atrx mÀz+ (n = 14) PN5 zygotes.
Data information: All scale bars, 10 lm. Source data are available online for this figure.

of 24
The EMBO Journal 39: e103697 | 2020 ª 2020 The Authors directly interacts with the middle to C-terminal parts of CBX2.1 or with other cPRC1 complex members, bound via its Pc-Box. To investigate these possibilities, we performed ChIP assays for CBX7, a CBX2 homolog sharing sequence similarity only in the CD and the Pc-Box ( Fig EV2A). As compared to GAL4DBD CBX2.1, we measured comparably low enrichment values for DAXX MycHis at the 5-UAS site when expressing either GAL4DBD CBX7 or GAL4DBD CBX2.2, though both proteins themselves were enriched at the site ( Fig EV2B). In co-IP experiments, as compared to CBX2, we observed only weak interactions between DAXX MycHis and CBX7 or with RYBP, a core component of the variant PRC1 complex (Garcia et al, 1999;Gao et al, 2012) (Fig EV2C). These results indicate that amino acids unique to CBX2.1, located between the AT hook and the Pc-Box, enable DAXX recruitment to chromatin.

CBX2 mediates DAXX targeting genome-wide
In HEK293 cells, endogenous CBX2 localizes to cPRC1 target genes (Gao et al, 2012). To address whether CBX2.1-mediated DAXX targeting could represent a general mechanism, we performed ChIPsequencing analyses in HEK293 5xUAS expressing both DAXX MycHis and GAL4DBD CBX2.1. We observed strong co-occupancy for CBX2.1 and DAXX at cPRC1 target genes like CCND2 and RUNX1 (Fig 3E) which we confirmed by ChIP-PCR ( Fig EV2D). Further global analysis revealed uniform co-occupancy of DAXX MycHis and GAL4DBD CBX2.1 at the TSS of many genes in a CBX2.1-dependent manner (Figs 3F-H and EV2E). In contrast, TSS occupancy levels of the shorter splice variant GAL4DBD CBX2.2 protein were lower compared to those of the full-length GAL4DBD CBX2.1 protein suggesting local stabilization of the CBX2.1 protein by interaction with other PRC1 components (Figs 3F and EV2F). Moreover, CBX2.2 failed to recruit DAXX to such targets ( Fig EV2G).
We also observed enrichments of the anti-MYC antibody in narrow peaks in the absence of exogenous CBX2.1 or CBX2.2 expression (red dots in Figs 3G and H,and EV2G and H). Further comparison to data from independent ChIP-sequencing studies for the transcription factor MYC shows that such TSSs are bound by MYC in different cell types, suggesting reactivity of the anti-MYC antibody to endogenously expressed MYC protein in HEK293 cells (ENCODE Project Consortium, 2012).
DAXX targeting to pat-PCH requires its SUMO-interacting motifs DAXX contains two conserved SIMs, enabling binding to SUMOylated proteins, such as PML ( Fig EV3A) (Lin et al, 2006;Santiago et al, 2009;Morozov et al, 2017). To study the role of SUMOylation in paternal heterochromatin formation, we performed IF for SUMO modifications in zygotes. While we failed to detect SUMO1 at PCH, we observed, in both early and late zygotes, enrichment for SUMO2/3 at pat-PCH but not mat-PCH. The SUMO2/3 signal colocalized with CBX2 and ATRX (Figs 4A and EV3B). To investigate the relevance of SUMOylated protein recognition in DAXX targeting to pat-PCH, we generated a DAXX mutant protein lacking both SIMs (DAXX DSIM1/2 ). As expected, this mutant did not localize to PML/ SUMO2 foci in mouse ES cells compared to the wild-type DAXX protein ( Fig EV3C). When expressed in HEK293 5xUAS cells, DAXX DSIM1/2 binding to the 5xUAS site was strongly decreased compared to wild-type DAXX, despite the presence of GAL4DBD CBX2.1 (Fig 4B). Since both wild-type and mutant DAXX were expressed at similar levels ( Fig EV3D), these results suggest that DAXX SIM domains are required for efficient chromatin targeting by CBX2-containing cPRC1. Enrichment of DAXX DSIM1/2 to canonical CBX2 target genes such as CCND2 and RUNX1 was also reduced, albeit to a lower extent possibly due to other interactions with endogenous and chromatin factors bound to these loci ( Fig EV3E).
To determine whether the SIMs of DAXX are required for targeting to pat-PCH in early embryo as well, we performed genetic complementation experiments in which we expressed exogenous DAXX or DAXX DSIM1/2 , translated from microinjected mRNAs, in newly generated Daxx conditionally deficient embryos, lacking maternal Daxx (Daxx mÀz+ ) expression ( Fig EV3F). While DAXX protein was not detectable in Daxx mÀz+ zygotes, cPRC1-catalyzed H2AK119ub1 levels at pat-PCH were unaltered showing that Daxx mÀz+ deficiency does not abrogate cPRC1 targeting and catalytic function (Fig EV3G and H). As expected, exogenously provided DAXX-EGFP localized at pat-PCH in Daxx mÀz+ zygotes in a manner comparable to that of endogenous protein in wild-type zygotes (Figs 1, 4C and EV1G). In contrast, DAXX DSIM1/2 -EGFP failed to localize to pat-PCH despite being properly expressed in zygotes ( Fig 4C and D), arguing that SUMOylation of a pat-PCH-associated protein underlies DAXX recruitment.
Source data are available online for this figure.

H P V P QK K AQ I L V A K P E L K D P H P V P QK K AQ I L V A K P E L K D P
146 147 154 CBX2 mouse CBX2.1 human

of 24
The EMBO Journal 39: e103697 | 2020 ª 2020 The Authors SUMO2-modification of CBX2 drives DAXX binding to pat-PCH CBX2 has been identified as a SUMO2-modified protein in human cells (Hendriks et al, 2014;Lamoliatte et al, 2014;Tammsalu et al, 2014). Sequence comparison shows that lysines K146, K153, and K362, which are SUMOylated in human CBX2, are conserved in the mouse protein ( Fig 4E). To study SUMOylation of mouse CBX2, we purified 6xHis-tagged SUMO-conjugated proteins from HEK293 cells co-expressing mouse CBX2 fused to EGFP. We show that mouse CBX2-EGFP is indeed conjugated to SUMO2, but not SUMO1 or SUMO3 ( Fig 4F). In contrast, CBX2 lacking the Pc-Box, required for the interaction of CBX2 with RNF2 and other cPRC1 components, failed to become SUMOylated (Fig EV3I), suggesting that incorporation of CBX2 in Polycomb complexes is required for its SUMOylation. Sequence analysis of mouse CBX2 using SUMOylation prediction software suggested additional putative SUMOylated sites ( Fig 4E) (Zhao et al, 2014). Thus, we tested combinations of lysine to arginine point mutations at several predicted SUMOylated lysines in mouse full-length CBX2. Mutation of at least four lysines (K146R, K147R, K154R, and K350R, referred to as CBX2 4KR ) was required to reduce significantly SUMO2 levels on CBX2 in vitro (Figs 4G and EV3I).
In wild-type zygotes, exogenously expressed CBX2 4KR -EGFP was bound to pat-PCH (Fig 4H and I) and along paternal mitotic chromosomes of the first cleavage division in a similar manner as CBX2-EGFP ( Fig EV3J) (Puschendorf et al, 2008), indicating that mutation of the four lysines did not impair CBX2 binding to paternal chromatin. In contrast, recruitment of endogenous DAXX to pat-PCH was impaired in these conditions arguing for a dominant negative effect of exogenous CBX2 4KR -EGFP over endogenous CBX2 in zygote embryos (Fig 4H and I). Together, these data demonstrate that (i) mouse fulllength CBX2 is a genuine SUMO2 substrate, (ii) SUMOylation of CBX2 requires its interaction with RNF2 and possibly other cPRC1 members, and (iii) SUMO2 on CBX2 mediates targeting of DAXX to pat-PCH.
The H3K9me3/HP1 pathway prevents DAXX/CBX2-PRC1 targeting at maternal PCH When compared to pat-PCH, DAXX is only weakly associated with mat-PCH despite robust localization of ATRX ( Fig 1B). This is surprising given the efficient recruitment of DAXX by ATRX to H3K9me3-modified chromatin at satellite and telomeric sequences in somatic cells Hoelper et al, 2017). Previously, we demonstrated that maternal deficiency of Suv39h2 or Hp1b, causing a loss of H3K9me3 and/or of HP1b localization at mat-PCH, resulted in the recruitment of the CBX2-containing canonical PRC1 complex to mat-PCH (Puschendorf et al, 2008;Tardat et al, 2015). Here, we show that DAXX becomes abundantly localized to mat-PCH in late zygote embryos of both mouse mutants as well, mimicking the dotted "pearls on a string" pattern as seen at pat-PCH in wild-type zygotes, downstream of CBX2 recruitment (Fig 5A and B). Similarly, ATRX co-localized with DAXX at mat-PCH in the "pearls on a string" configuration in both mutants, suggesting DAXX-dependent recruitment. Next, we tested the impact of the CBX2 F12A -EGFP mutant protein on DAXX localization (Tardat et al, 2015). The F12A mutation in the aromatic cage of the CD of CBX2 enables the neighboring AT hook to efficiently bind to AT-rich major satellite sequences even in the presence of H3K9me3 and HP1b, thus overcoming the steric hinderance imposed by the CD of CBX2 (Tardat et al, 2015). Expression of CBX2 F12A -EGFP in wild-type zygotes induced a strong enrichment of DAXX at mat-PCH, also in a dotted pattern, despite the presence of HP1b ( Fig 5C). Therefore, these results indicate that the CBX2/ cPRC1 ? SUMO2 ? DAXX ? ATRX recruitment pathway even functions in a H3K9me3/HP1b chromatin environment. It further indicates that DAXX recruitment via the H3K9me3/ Hp1b ? ATRX ? DAXX pathway is inefficient in early embryos. These data argue for different modes of interactions between the ATRX and DAXX proteins and/or possible regulation thereof in the two pathways.
DAXX and CBX2-PRC1 ensure stability of major satellite repeats of pat-PCH Daxx has been implicated in controlling expression of genes and repetitive sequences and in nuclear organization in ESCs and differentiated cells (Drane et al, 2010;Rapkin et al, 2015;Hoelper et al, 2017). To address the role of DAXX in early embryogenesis, we generated embryos lacking either maternal (Daxx m-z+ )or ◀ Figure 4. SUMOylation of CBX2 mediates DAXX targeting to pat-PCH.
Right: Absolute intensities of EGFP in same samples. In boxplots, the enter lines show the medians; box limits indicate the 25 th and 75 th percentiles as determined by R software; and whiskers extend 1.5 times the interquartile range from the 25 th and 75 th percentiles. Datapoints are indicated by circles. Experiments were replicated two times. ***P < 0.001, t-test.
Data information: Scale bars, 10 lm. Source data are available online for this figure.
ª 2020 The Authors The EMBO Journal 39: e103697 | 2020 maternal and zygotic expression of Daxx (Daxx mÀzÀ ), by crossing conditional Daxx F/F ; Zp3-cre females with heterozygous Daxx +/À males ( Fig EV3F) (Michod et al, 2012). Daxx mÀzÀ deficiency resulted in major impairment of pre-implantation development, starting at the 8-cell stage and reducing the efficiency of blastocyst development by 50% (Fig 6A). The developmental impairment correlated well with major defects in chromosome segregation, cytokinesis, and occurrence of endo-reduplication in Daxx mÀzÀ morulae and blastocysts (Fig 6B and C). In contrast, Daxx mÀz+ embryos developed normally to the blastocyst stage. This finding is consistent with an enrichment of DAXX at pericentromeric chromocenters in Daxx mÀz+ 4-cell stage embryos (Fig EV4A), arguing for a rescue of the developmental phenotype by expression of the paternal Daxx allele. Nonetheless, we noticed micronuclei in both Daxx mÀzÀ and Daxx mÀz+ 4-cell embryos suggesting the occurrence of chromosomal segregation defects at earlier stages as well (Fig EV4B). To investigate the process and mechanism underlying micronuclei formation, we performed live-imaging experiments on zygotes expressing H2B-mCherry and Tubulin-EGFP. We measured delays in segregation of one or more chromosomes during the first cleavage division in 33% of Daxx mÀz+ embryos and 25% of Ring1 mÀz+ ; Rnf2 mÀz+ embryos lacking PRC1 expression (Fig 6D and E, Movie EV4). All lagging chromosomes or micronuclei were of paternal origin, as revealed by staining for 5-methylcytosine and 5-hydroxymethylcytosine as markers for parental identity (Inoue & Zhang, 2011) and by the absence of ATRX labeling at PCH in resulting micronuclei (Figs 6F and EV4B and C). To identify sequences subjected to instability, we probed by DNA fluorescence in situ hybridization (FISH) major and minor satellites demarcating pericentromeric and centromeric regions, respectively, during the metaphase-to-anaphase transition of the first cleavage division. All lagging chromosomes from anaphase Daxx mÀz+ embryos marked for major satellites sequences lacked signals for minor satellites indicating breakage of proximal ends (Fig 6G) Figure 5. CBX2 and the Suv39h/Hp1 pathway control DAXX targeting to pat-PCH in mouse zygotes.
Data information: Panels (A and C): Paternal and maternal pronuclei of a given condition originate from same zygotes. All scale bars 10 lm. Source data are available online for this figure.
10 of 24 The EMBO Journal 39: e103697 | 2020 ª 2020 The Authors     these embryos showed indeed that major satellite sequences undergo stretching and breakage in the absence of DAXX or cPRC1 (Fig EV4D, Movie EV5). Consistently, we detected some mini-chromosomes positive for FISH signals of major and minor satellites ( Fig 6G). We also observed that lagging chromosomes and micronuclei from Daxx mÀz+ embryos at the 1-cell to 2-cell transition stage lacked CREST staining. Altogether, these data argue that DAXX functions to ensure PCH stability of paternal chromosomes in zygotes.

DAXX controls H3.3 deposition at pat-PCH in a PRC1-dependent manner
Since DAXX functions as an H3.3 chaperone, we investigated whether H3.3 incorporation was altered in Daxx-deficient embryos. Due to our inability to immuno-label endogenous H3.3 in zygotes, we microinjected mRNA encoding for H3.3-EGFP in both wild-type and maternally deficient Daxx mÀz+ embryos and assayed its incorporation at PCH localized around nucleolar precursor bodies in paternal pronuclei via live imaging. We observed a clear reduction of H3.3-EGFP fluorescence signal at such sites in embryos maternally deficient for Daxx or for Ring1 and Rnf2 (Fig 6H; Movies EV6-EV8). These data show that DAXX mediates H3.3 incorporation at pat-PCH in a PRC1-dependent manner.
DAXX-mediated H3.3 chaperone function safeguards nucleosome density and chromosomal stability of pat-PCH domains Since H3.3 is the predominant histone variant deposited at pat-PCH in early zygotes (Fig 1A) (Akiyama et al, 2011), we hypothesized that nucleosome density at pat-PCH would be impaired in Daxx mÀz+ embryos lacking local H3.3 incorporation. Indeed, we measured reduced levels of H3 and ATRX signals at pat-PCH in Daxx mÀz+ late zygote embryos compared to their wild-type counterparts (Fig 6I and J), arguing that DAXX-mediated H3.3 deposition is required to control nucleosome density at pat-PCH.
To assess for a direct role of H3.3 deposition in maintaining genome integrity, we generated the R244A substitution mutation in the histone-binding domain of mouse DAXX, which reduces its binding capacity for H3.3 (Elsasser et al, 2012). Relative to wild-type DAXX-EGFP, we co-immunoprecipitated reduced amounts of endogenous H3.3 with the DAXX R244A -EGFP mutant from HEK293 cell extracts ( Fig 6K). Next, we expressed DAXX-EGFP and DAXX R244A -EGFP at different levels in Daxx mÀz+ embryos from DNA replication onwards by microinjection of corresponding mRNAs and measured the impact on chromosome segregation. While injection of Tubulin-EGFP or low levels of Daxx-EGFP mRNA failed to rescue the segregation defects in the first cleavage division of Daxx mÀz+ embryos, higher levels of Daxx-EGFP were sufficient to do so. In contrast, high expression of DAXX R244A -EGFP, despite being recruited to pat-PCH normally, did not rescue the segregation defects (Figs 6L-N and EV4E). Quantification of total H3 IF signals revealed significantly higher levels at pat-PCH when expressing wild-type DAXX compared to the DAXX R244A mutant in Daxx mÀz+ embryos (Fig 6M and N). We therefore conclude that DAXX-mediated H3.3 incorporation into chromatin is essential for maintaining nucleosome density at pat-PCH and ensuring stability of the underlying major satellite sequences during early embryonic divisions.

DAXX deficiency abrogates forward and reverse expression of major satellites
As we reported previously, PRC1 represses nascent transcription of major satellite at pat-PCH in mouse zygotes (Puschendorf et al, 2008). To assess the role of DAXX in forward and reverse ◀ Figure 6. DAXX and PRC1 control nucleosome density and chromatin stability at pat-PCH and chromosome segregation in early mouse embryos.
A Developmental progression of Daxx m+z+ , Daxx mÀz+ , and Daxx mÀzÀ early embryos. B Examples of morula and blastocyst embryos with normal and large polyploid nuclei. C Percentages of embryos with normal or polyploid nuclei. D Still images of time-lapse imaging of the first cleavage division in wild-type and maternally deficient zygotes expressing H2B-mCherry. Time points represent minutes after prometaphase. Lagging chromosomes and micronuclei are indicated by yellow arrowheads. E Percentage of zygotes with lagging chromosomes during 1st cleavage division. Data represent mean AE SEM in n ≥ 3 replicate experiments. F IF of 5-methylcytosine and 5-hydroxymethylcytosine staining of wild-type (n = 42), Daxx mÀz+ (n = 14) and Ring1 mÀz+; Rnf2 mÀz+ (n = 6) zygotes with lagging chromosomes. Yellow arrowheads indicate lagging chromosomes. G DNA FISH for major (green) and minor (red) satellites on cleavage chromosomes of a Daxx mÀz+ zygote with a chromosome broken within the major satellite region. H H3.3-EGFP live-imaging signal in pat-and mat-PCH of control Daxx m+z+ (n = 18), Daxx mÀz+ (n = 8), and Ring1 mÀz+ ; Rnf2 mÀz+ (n = 12) zygotes. Scale bar, 10 lm. 12 of 24 The EMBO Journal 39: e103697 | 2020 ª 2020 The Authors transcription of major satellites (Probst et al, 2009;Casanova et al, 2013), we performed RNA FISH with strand-specific probes in late zygotes. Surprisingly, expression of both strands was majorly downregulated in Daxx mÀz+ embryos (Fig 7A and B), showing a crucial function for DAXX and/or H3.3 deposition in enabling major satellite transcription and/or transcript retention at pat-PCH. Moreover, these data argue for a dominant role of PRC1 and possibly H2AK119 ubiquitination over DAXX in transcriptional repression of major satellites repetitive sequences.
DAXX contributes to PRC1-mediated gene repression in oocytes and 2-cell embryos To determine whether DAXX regulates gene expression during oocyte growth and zygotic genome activation in early embryos, we performed RNA sequencing on Daxx deficient and control MII oocytes as well as Daxx mÀz+ and control 2-cell embryos ( Fig EV5A). Compared to control oocytes, Daxx transcripts levels were over 10-fold reduced in Daxx mutant oocytes, with > 65% of residual transcripts lacking the floxed exon 3, showing an efficient depletion of Daxx expression. In contrast, Daxx mRNA levels were only < 2-fold reduced in maternally deficient 2-cell embryos, with > 65% being full length, indicating early and potent transcriptional activation of the paternal allele in Daxx mÀz+ embryos (Fig EV5B), which relates to immuno-detectable DAXX protein from the four-cell stage onwards (Fig EV4A and B).
In Daxx mutant oocytes, 25 and four genes were at least twofold up-or downregulated, respectively (minlogFC = 1; FDR < 0.05; Figs 7C and EV5C, Table EV1). In Daxx mÀz+ 2-cell embryos, only six genes were upregulated, which may relate to clear paternal Daxx expression (minlogFC = 1; FDR < 0.05; Figs 7D and EV5B and C). Four of these genes were also upregulated in mutant oocytes, suggesting inheritance of some aberrantly expressed maternal mRNAs (Fig EV5C). To investigate whether aberrant gene expression could contribute to the chromosome instability observed in Daxx mÀz+ embryos, we performed gene ontology analyses for genes differentially expressed in oocytes and embryos as identified using slightly more relaxed significance criteria (minlogFC = 1; P < 0.01; Fig 7C and D; Table EV2). We did not observe any enrichment for biological process or molecular function terms related to genome instability, chromosome segregation failure, or cell cycle for up-or downregulated genes. These data argue against a potential indirect transcriptional effect and support the model that maternal DAXX maintains pat-PCH stability in embryos directly, by controlling local nucleosome density and chromatin configuration.
Intriguingly, GO terms related to development and differentiation were over-represented among differentially expressed genes (Table EV2), possibly hinting to a function of DAXX in PRC1mediated gene regulation. To investigate this notion further, we calculated the occupancy of H3K27me3 around transcriptional start sites of genes in wild-type oocytes, and classified genes as being lowly, mediumly, or highly trimethylated at H3K27 (Fig 7E)  . As anticipated, benchmarking this chromatin classification to the transcriptional response of genes to Ring1; Rnf2 double deficiency (dko) in GV oocytes revealed a clear positive relation (Fig 7F and G) (Posfai et al, 2012). While H3K27me3 highly marked genes were majorly upregulated in PRC1 dko oocytes, H3K27me3 lowly marked genes were not (Fig 7G). Importantly, we observed a significant, positive association between upregulated gene expression in Daxx-deficient oocytes and embryos and H3K27me3 levels in wild-type oocytes as well, indicating a regulatory role for DAXX in Polycomb-mediated gene repression (Fig 7H and I). Indeed, direct comparison of differential expression in Daxx-deficient versus control MII oocytes relative to Ring1; Rnf2 dko versus control GV oocytes showed significant positive associations (Fig 7J). Genes upregulated in Ring1; Rnf2 dko oocytes were generally upregulated in Daxx mutant oocytes ( Fig 7J) and vice versa ( Fig EV5D) and include various transcription and signaling factors such as Ascl1, Bex3, Dlx5, Dmrt2, Lhx6, Jun, Klf12, Lmo4, Myb, Notch1, Sox2, Sox6, Tbx2, Tbx18, Tcf19, Tfec, Tle3, and Zbtb7c. Together, these results indicate that DAXX contributes to gene repression by PRC1 complexes in oocytes, plausibly by interacting with CBX2 and controlling DAXX-mediated nucleosome deposition at PRC1-target genes during oocyte development.

Discussion
In this study, we identify and characterize a novel chromatin pathway required for maintaining genome integrity at major satellite tandem repeats in mouse early embryos. Our work elucidates a modular "plug and play" principle in chromatin biology in which SUMOylation functions as a molecular glue that drives interactions between components of the H3K9 methylation and Polycomb pathways that to date have been considered to function separately. We demonstrate that chromatin of major satellite sequences becomes modified by SUMO2/3 shortly after fertilization and throughout the first cell cycle, in a parent-of-origin-specific manner only at the paternal genome. Our biochemical and molecular genetic experiments show that CBX2, as part of canonical PRC1, is the major and functionally relevant substrate of SUMOylation at pat-PCH. CBX2 SUMOylation drives recruitment of DAXX via its SUMO-interacting domains (SIMs). DAXX enables H3.3 incorporation at pat-PCH marked by cPRC1, independently or possibly in conjunction with ATRX. We unambiguously show that H3.3 deposition is essential for maintaining nucleosomal density and DNA stability at pat-PCH. We also demonstrate that DAXX contributes to repression of genes silenced by PRC1 in oocytes and 2-cell embryos. These data suggest a modulatory function of the novel CBX2/ cPRC1 ? SUMO2 ? DAXX ? H3.3 deposition module in Polycomb gene silencing (Fig 7K).
The asymmetry in heterochromatic states at pat-PCH and mat-PCH within one cell is remarkable. While on the maternal genome a H3K9me3 state is inherited from the oocyte and maintained in the embryo, a canonical PRC1 state linked to DAXX-mediated H3.3 deposition is newly formed at pat-PCH due to the absence of H3K9me3 in sperm-derived chromatin (Puschendorf et al, 2008;Tardat et al, 2015). Here, cPRC1 is initially recruited via the binding of CBX2 to AT-rich major satellite sequences (Tardat et al, 2015). This interaction is further strengthened via the interaction of CBX2 with de novo PRC2-mediated H3K27me3 established during the S and G2 phases of the zygote cell cycle (Albert & Peters, 2009;Tardat et al, 2015). Though lysine 27 of H3.3 and its methylation was suggested to be directly required for heterochromatin formation in early embryos (Santenard et al, 2010), we show here that PRC2 and ª 2020 The Authors The EMBO Journal 39: e103697 | 2020   Figure 7.
14 of 24 The EMBO Journal 39: e103697 | 2020 ª 2020 The Authors therefore de novo H3K27me3 at pat-PCH are not critical to DAXX recruitment (and thus H3.3 deposition) to pat-PCH. PRC2 is also superfluous for PRC1 recruitment to pat-PCH (Puschendorf et al, 2008). The developmental defects induced by exogenous H3.3K27R (Santenard et al, 2010) likely relate to the interference of other H3.3K27 functions, such as regulating genome activation by lysine 27 acetylation (Dahl et al, 2016). It remains to be determined which factors control the target specificity of SUMOylation at paternal and maternal PCH in early embryos. The basic machinery for this is available as the genes encoding SUMO E1-activating enzymes Uba2 and Sae1, the E2conjugating enzyme Ube2i (encoding for UBC9), and three Sumo1, Sumo2, and Sumo3 and several sentrin/SUMO-specific proteases (Senp), removing SUMOylation, are expressed in oocytes and early embryos ( Fig EV1D) (Kang et al, 2010). The identity of the E3 ligase for CBX2 is unclear. It is highly unlikely to be the Polycomb protein Pc2 (CBX4), a paralog of CBX2 and E3 SUMO1 ligase for the transcriptional repressor CtBP and the kinase HIPK2 (Kagey et al, 2003), since Cbx4 is not expressed during oogenesis and in zygote stage embryos (Fig EV1D). In contrast, four E3 SUMO ligases belonging to the Pias gene family are all expressed at these stages ( Fig EV1D). Importantly, in mouse fibroblasts the H3K9me3 HMT SUV39H1 (but not its paralog SUV39H2), functioning as an E3 SUMO1 ligase for HP1 proteins, and the SUMO protease SENP7 regulate HP1a targeting to canonical H3K9me3/HP1 constitutive heterochromatin (Maison et al, 2011(Maison et al, , 2012(Maison et al, , 2016Romeo et al, 2015). In early embryos, however, we did not detect SUMO1 staining at mat-PCH indicating that this pathway is dysfunctional at this early stage of development. This may be due to that Suv39h2 is more prominently expressed than Suv39h1 in zygotes. Moreover, HP1a (CBX5), though expressed, is not localized to chromatin in zygotes (Tardat et al, 2015). Finally, Senp7 is not expressed in oocytes and zygotes, possibly explaining the absence of HP1a in early embryonic chromatin.
At pat-PCH, the order of DAXX and ATRX recruitment to chromatin is inversed relative to mat-PCH in embryos and to PCH in somatic cells. In mouse fibroblasts, ATRX is constitutively localized at PCH via its interactions with H3K9me3-/HP1-modified chromatin. The SIMs of DAXX drive the association of DAXX and ATRX at SUMOylated PML bodies in G1-, early S-, and G2-phase cells (Ishov et al, 2004;Shastrula et al, 2019). At the end of S phase, however, DAXX and ATRX are released from SUMOylated PML bodies and DAXX is transiently recruited to PCH, already marked by ATRX. The S phase restricted accumulation of DAXX at PCH is inhibited by roscovitine treatment, blocking CDK2-mediated phosphorylation of ATRX and possibly other proteins (Ishov et al, 2004). Biochemically, the N-terminal four-helix bundle domain (4HB) of DAXX interacts directly with a short evolutionary highly conserved fragment of ATRX (Hoelper et al, 2017). It will be important to study whether DAXX and ATRX at pat-PCH and mat-PCH in zygotes are differentially modified by post-translational modifications, explaining the parent-of-origin-specific targeting dependencies and cell cycle kinetics of PCH enrichments. Intriguingly, DAXX and ATRX co-localize at multiple foci in euchromatin of zygotes, particularly in maternal pronuclei (Fig 1). At such foci, the co-localization is dependent on Daxx but not Atrx suggesting that DAXX may define the chromatin state at other sequences in the genome similarly as at pat-PCH.
We demonstrate that DAXX recruitment to pat-PCH is required for maintaining nucleosome density at major satellite sequences. The absence of DAXX in Daxx-deficient or Ring1; Rnf2 dko embryos leads to chromatin unfolding and DNA breakage at major satellites when paternally inherited chromosomes segregate to opposite poles during the first cleavage divisions. This phenotype is suppressed by expressing DAXX from S phase onwards suggesting that the timing of function of the CBX2/cPRC1 ? SUMO2 ? DAXX ? H3.3 deposition pathway during late S phase and/or G2 is sufficient. This timing within the embryonic cell cycle is similar to the reported ◀ Figure 7. DAXX and PRC1 control nucleosome density and chromatin stability at PCH in early mouse embryos.
A RNA FISH of forward and reverse strand major satellite sequences in Daxx m+z+ and Daxx mÀz+ PN5 zygotes. Scale bars, 10 lm. B Box plots indicating the foci number (left) and total intensity (right) of forward (top) and reverse (bottom) transcription as measured by RNA FISH on 3D paternal pronuclei in Daxx m+z+ and Daxx mÀz+ PN5 zygotes. In boxplots, the center lines show the medians; box limits indicate the 25 th and 75 th percentiles as determined by R software; and whiskers extend 1.5 times the interquartile range from the 25 th and 75 th percentiles. Datapoints are indicated by circles. Experiments were replicated three times. ***P < 0.001; t-test. C, D MA-plots showing fold change (FC) (in log 2 ) in gene expression between indicated samples as a function of averaged expression across all samples. Dotted lines indicated log 2 (FC) of 1. E Histogram showing classification of genes into three groups according to distribution of ChIP enrichment (log 2 ) of H3K27me3 around TSS in MII oocytes . Lower hinge, central line, and upper hinge represent 25 th , 50 th (median), and 75 th percentiles, respectively. Upper/lower whiskers extend to the largest/smallest values no further than 1.5*IQR from the upper/lower hinge, where IQR is interquartile range or distance between 25 th and 75 th percentiles. Outliers are not displayed. The notches extend 1.5*IQR/√n where n is a number of genes in each H3K27me3 class. A number of genes (n) in "Low", "Medium", and "High" H3K27me3 class used for boxplots are 9,824, 2,311, and 2,647, respectively. Expression changes were estimated using two or three biological replicates for each condition. Statistical significance was estimated using two-sided t-tests. J Boxplot showing FC (log 2 ) in gene expression between Daxx-deficient and control MII oocytes for five groups of genes classified according to FC (log 2 ) in gene expression between Ring1 Rnf2-deficient and control GV oocytes. Boxplots were created similarly to (G-I), and the number of genes used for each boxplot is displayed below. Expression changes were estimated using three biological replicates for each condition. Statistical significance was estimated using two-sided ttests. K Schematic diagram showing molecular pathways controlling repression and chromatin integrity at paternal and maternal major satellite sequences in mouse early pre-implantation embryos.
Source data are available online for this figure.
ª 2020 The Authors The EMBO Journal 39: e103697 | 2020 relocation of DAXX to PCH in fibroblasts during late S phase (Ishov et al, 2004). Nonetheless, DAXX was also shown to deposit H3.3 nucleosomes onto major satellites prior to S phase in fibroblast (Drane et al, 2010). Thus, given its localization in pat-PCH of early zygotes (Fig 2), DAXX may also deposit H3.3 during sperm decondensation and pronuclear formation. Together, our data underscore the relevance of the H3.3-H4 chaperone function of DAXX for controlling embryonic chromatin integrity. In mouse early embryos, major satellites are predominantly expressed from the forward strand in paternal genomes, while the reverse strands are expressed in both genomes in a developmentally coordinated manner (Probst et al, 2010). In Daxx mÀz+ zygotes, transcripts from both strands were undetectable by RNA FISH suggesting that DAXX is generally required for transcription of major satellites (Fig 7). Indeed, reduced major satellite expression has also been observed in Daxx-deficient fibroblasts (Drane et al, 2010). In strong contrast, maternal deficiency of cPRC1 enhances major satellite expression pat-PCH (Puschendorf et al, 2008). Thus, though working in the same pathway, DAXX and cPRC1 serve opposing roles in balancing expression of the major satellite tandem repeats. Since mitotic stability of pat-PCH is affected in Daxx mÀz+ as in Ring1 mÀz+ ; Rnf2 mÀz+ embryos, it is likely that the derailment of transcription of major satellites does not cause heterochromatin destabilization at pat-PCH.
By exposing embryos to LNA-DNA gapmers from the zygote stage onwards, knockdown of reverse but not of forward strand major satellite transcripts has been shown to impair the development of regular and parthenote embryos beyond the two-cell stage (Probst et al, 2010;Casanova et al, 2013). Moreover, the LNA-DNA gapmer treatment of specifically reverse transcripts interferes with the clustering of major satellites into chromocenters normally occurring in 2-and 4-cell embryos (Probst et al, 2010;Casanova et al, 2013). These studies underscore the developmental importance of certain transcriptional and chromatin reprogramming events occurring at mat-PCH. In Daxx mÀz+ embryos, however, neither zygotic expression of major satellites at mat-PCH, nor the two-to four-cell transition or the formation of chromocenters was altered ( Fig EV4F). In contrast, maternally provided ATRX has been implicated in suppressing major satellite expression at mat-PCH (De La Fuente et al, 2015). Moreover, reduction of maternal Atrx induces mitotic recombination at centromeric minor satellite sequences and chromosome instability in mouse early embryos (De La Fuente et al, 2015). Unfortunately, it is unclear whether these instability phenotypes are restricted to the maternal genome or apply to the paternal genome as well. Even so, we did not observe defects in major satellite stability and segregation of maternal chromosomes in Daxx mÀz+ embryos, supporting the notion that the low DAXX levels and absence of H3.3 deposition at mat-PCH do not compromise H3K9me3/HP1 heterochromatin integrity. Our results also suggest that ATRX functions independently of DAXX at mat-PCH.
In 2004, Ishov et al reported that~10% of Daxx deficient embryonic fibroblasts contained two nuclei, a level threefold higher than in wild-type cells. We observed that lack of maternal and paternal Daxx expression severely impaired progression of embryos toward the blastocyst stage. Notably,~45% of Daxx mÀzÀ morula and blastocysts contained large polyploid nuclei compared to 18% of wild-type and Daxx mÀz+ embryos, pointing to major defects in chromosome segregation and/or cytokinesis in the absence of DAXX expression later in development (Fig 6). The more dramatic ploidy phenotype in later stages suggests a different molecular role for DAXX in maintaining genome integrity. Interestingly, the canonical H3K9me3/HP1 state as present at mat-PCH is de novo established at pat-PCH at the 8-cell stage and co-localizes with DAXX (Puschendorf et al, 2008). Likewise, we detected PML bodies from the 8-cell stage onwards as well, which co-localize with DAXX ( Fig EV4G). Interestingly, morpholino-mediated depletion of total H3.3 protein showed that H3.3 sustains proper chromosome segregation throughout pre-implantation development. Moreover, H3.3 maintains a decondensed chromatin state, in part by directing MOF-mediated H4K16 acetylation and counteracting H1 incorporation, and which is required for development beyond the morula stage (Lin et al, 2013). To understand the possibly multiple mechanisms underlying the developmental and cytokinesis defects observed in Daxx mÀzÀ embryos, it will thus be important to assess the impact of Daxx deficiency on general mitotic chromatin structure as well as on gene expression during late preimplantation development.
Our study in HEK293 cells shows that exogenous CBX2 can recruit DAXX to and drive H3.3 deposition at a heterologous gene target in a cPRC1-dependent manner (Fig 4). The recruitment is likely mediated via SUMOylation given that it depends on the SIM domains of DAXX. Our genome-wide analysis shows that exogenous CBX2 recruits DAXX to transcriptional start sites of many genes in HEK293 cells. To identify a potential role of this interaction in a physiologically relevant setting, we compared differential gene expression in Daxx-deficient versus Ring1; Rnf2 double deficient oocytes and 2-cell embryos where CBX2 is naturally highly expressed. These analyses revealed many commonly upregulated, but not downregulated genes (Fig 7). Moreover, showing further specificity of the correlation, we did not observe any correspondence in derailed gene expression in oocytes deficient for Daxx versus Setdb1, a well-known histone H3 lysine 9 methyltransferase and SUMOylated interaction partner of DAXX in mouse ESCs (Fig EV5E and F) (Hendriks et al, 2014;Eymery et al, 2016).
Cbx2 is one of five mammalian orthologues (Cbx2, 4, 6, 7, 8) of the Drosophila Polycomb Pc gene. In post-implantation embryos, it is expressed in most tissues and serves non-redundant gene regulatory functions (Schoorlemmer et al, 1997;Lau et al, 2017). For example, Cbx2 is required for proper axial patterning of the skeleton, for limb development (Core et al, 1997) as well as for male sex determination in mice and men (Katoh-Fukui et al, 1998;Biason-Lauber et al, 2009). In developing oocytes, Cbx2 is the only Pc orthologue that is highly expressed (Fig EV1D) (Puschendorf et al, 2008). Accordingly, a high level of maternal Cbx2 transcripts is contributed from the egg to the embryo (Puschendorf et al, 2008). The other orthologs become progressively activated during development toward the blastocyst stage (Fig EV1D). Our data demonstrate that DAXX specifically contributes to PRC1-mediated gene repression in oocytes and early embryos suggesting an unique role for CBX2 in cPRC1 function at the onset of embryogenesis.
CBX2 integrates many functions that are served by different domains of the protein (Fig EV5G). While the N-terminal CD confers binding to H3K27 and H3K9 trimethylated nucleosomes (Bernstein et al, 2006), the neighboring AT hook (ATH) controls binding to ATrich DNA (Tardat et al, 2015). The linker sequence between CD and ATH contributes to nucleosome binding (Tardat et al, 2015). Importantly, phosphorylation of the N-terminal serine-rich region (SRR) prevents binding of the adjacent ATH to major satellite doublestranded DNA and single-stranded RNA, while it promotes CD binding to H3K27me3-modified nucleosomes (Kawaguchi et al, 2017). Interestingly, the CD-ATH domain was shown to prevent binding to H3K9me3-modified chromatin containing HP1 protein, presumably via steric hindrance requiring the linker sequence and the ability of the CD to bind to methylated histones (Tardat et al, 2015). The CD-ATH domain also enables CBX2 binding to mitotic chromatin, independent of H3K27me3, PRC2, or other cPRC1 components (Fig EV5G) (Zhen et al, 2014). For immobilization along mitotic chromosomes, however, the C-terminal Pc-box is needed (Zhen et al, 2014) which drives interactions with cPRC1 components and mediates gene silencing (Schoorlemmer et al, 1997). The protein contains a large low complexity disordered region (LCDR) with several positively charged lysine residues which are required for biochemical chromatin compaction in vitro and axial patterning in vivo (Grau et al, 2011;Lau et al, 2017). Recently, CBX2 has been shown to undergo phase separation in vitro (Plys et al, 2019;Tatavosian et al, 2019). This phase separation is dependent on the positively charged lysine residues in the LCDR (Plys et al, 2019). Moreover, mutation of ATH, ATHL1 or 2, or of the SRR is sufficient to reduce or fully abrogate phase separation in vitro as well as the formation of CBX2containing foci in cells (Tatavosian et al, 2019). Interestingly, the major SUMOylation sites on CBX2 are located between the ATHL domains suggesting that SUMOylation could perturb charge-based intra-and intermolecular interactions, thereby modulating phase separation, Polycomb body formation in vivo (Schoorlemmer et al, 1997;Bel et al, 1998) and possibly long range interactions (Schoenfelder et al, 2015). Indeed, in fly, Pc-containing foci coagulate in larger aggregates in the absence of SUMO and demonstrate slower recovery of Pc-GFP signal upon FRAP treatment. In contrast, genetic abrogation of the SUMO peptidase Veloren (Velo) function disperses Pc foci and reduces Pc occupancy on known PREs in wing imaginal disk cells (Gonzalez et al, 2014). Thus, given the presence of CBX2 foci at pat-PCH and euchromatin in mouse pre-implantation embryos (Puschendorf et al, 2008;Terranova et al, 2008;Tardat et al, 2015), the level of CBX2 SUMOylation may modulate its protein function. On the one hand, it may control interactions with SIM containing proteins like DAXX. On the other hand, it may fine-tune the level of interactions between PRC1 complexes thereby determining gene repression and nuclear organization in early embryos, throughout development and in disease (Schoenfelder et al, 2015;Kundu et al, 2017;Lau et al, 2017;Sproll et al, 2018).
Oocytes and zygotes maternally deficient for Daxx were generated by crossing Daxx flox/flox ; Zp3-cre females with Daxx flox/flox males. Maternally and zygotically deficient embryos (Daxx mÀzÀ ) were generated by crossing Daxx flox/flox ; Zp3-cre females with Daxx flox/D males. Control embryos were generated by crossing Daxx flox/flox females with Daxx flox/flox males. Animal housing and experiments were performed in accordance with the Swiss animal protection laws and institutional guidelines.
Daxx MycHis was a gift from P. Leder (Addgene plasmid #1852; Ecsedy et al, 2003). Daxx-EGFP was generated by PCR amplification of the Daxx ORF from Daxx MycHis , cloning this fragment into pCDNA3.1-polyA (Yamagata et al, 2005) and subsequently cloning an EGFP ORF after the 3 0 end of the Daxx cDNA.
Expression vectors for 6xHis-tagged SUMO1, SUMO2, and SUMO3 were a kind gift from L. Linares (Institut de Recherche et Cancerologie de Montpellier, France).

siRNA microinjection and intracytoplasmic sperm injection
Fully grown GV oocytes were collected from ovaries of 8-week-old C57BL/6J female mice by puncturing ovarian follicles with needles (22G × 1 ¼"). The cumulus-oocytes complexes were transferred to M2 medium (Sigma) containing 2.5 mM milrinone (Sigma), and the cumulus cells were gently removed with a narrow-bore glass pipette. The oocytes were then transferred to M16 medium (Sigma) containing 2.5 mM milrinone. The GV oocytes were injected with 10 lM siRNAs specific for Daxx (5 0 -AUCGAACCUUGGCCAUGAA) or 10 lM control scrambled siRNA (5 0 -GAUCCAAGAAAUGGUUAGA) (all from Microsynth, Switzerland). Microinjection of siRNAs to GV oocytes was performed using an Olympus IX71 microscope equipped with micromanipulators (Eppendorf, NK2) and FemtoJet (Eppendorf). Two hours after injection, oocytes were incubated in milrinone-free M16 medium at 37°C, 5% CO 2 and 5% oxygen to induce in vitro meiotic maturation (IVM). Fourteen hours after IVM, matured MII oocytes containing the first polar body were selected for ICSI.
ICSI was carried out according to the piezo-assisted micromanipulation method as described previously (Yoshida & Perry, 2007). In brief, MII oocytes were placed into a 10 ll "micromanipulation droplet" of M2 medium containing 5 lg/ml of cytochalasin B (Sigma). Micromanipulation droplets were covered by mineral oil and placed on the stage of a microscope (Olympus X71) equipped with a micromanipulator (Narishige) and a PiezoXpert (Eppendorf). A sperm head was aspirated into the tip of an injection pipette, with a 10 lm internal diameter, and then deposited into the cytoplasm of a MII oocyte by applying a pulse of piezo to penetrate the zona pellucida and plasma membrane of the oocyte. ICSI-generated zygotes were cultured in M16 medium at 37°C, 5% CO 2 , and 5% oxygen. Following 4-5 h of culture, two pronuclei were visible in fertilized zygotes.
For ESC staining, transiently transfected ESCs were seeded after approximately 16 h of culture on poly-L-lysine-coated coverslips prior to fixation with 2% PFA in PBS. After the immunolabeling with the indicated antibodies, slides were mounted in Vectashield containing DAPI (Vector).
Image acquisition was performed with a Zeiss LSM 700 confocal microscope equipped with a 63 × Plan-Apochromat objective (N.A. 1.4). Z-stacks were acquired using a frame size of 512 × 512 pixels, a pixel depth of 16 bits, and 0.33 lm z-step, with sequential single track scanning using 639, 555, 488, and 405 nm wavelengths of lasers.
3D DNA fluorescence in situ hybridization 3D DNA FISH was performed according to previously described protocols with modifications (Maalouf et al, 2009;Liu et al, 2012b). Unless otherwise specified, steps were performed at room temperature. The embryos at anaphase stage were regularly monitored and handled using a stereo-microscope (Leica, MZ16) from 26 hp-hCG onwards. Following disappearance of the nuclear membrane, the embryos were separated to individual M16 droplets and cultured for an additional 2 h before fixation. The zona pellucida of embryos was removed by acidic tyrode (Sigma) prior to fixation with 4% PFA in PBS. The embryos were washed in PBS and deposited on glass slides to allow adherence. They were re-fixed with 4% PFA in PBS for 10 min and permeabilized for another 30 min in 0.5% Triton X-100 in PBS. Following two 5-min washes in 2× saline sodium citrate (SSC), pH 6.3, embryos were incubated in 2× SSC containing 200 lg/ml RNase (Microsynth) for 30 min at 37°C. After a rinse in 2× SSC, the slide was equilibrated in the hybridization buffer (50% formamide, SSC 2×, Denhardt 1×, 40 mM NaH2PO4, 10% dextran sulfate) overnight at RT. The DNA probes were denatured in hybridization buffer solution at 85°C for 10 min and incubated at 37°C for 30 min. Embryos were denatured in hybridization 18 of 24 The EMBO Journal 39: e103697 | 2020 ª 2020 The Authors buffer solution at 85°C for 10 min and incubated with denatured DNA probes for 24 h at 37°C in a humidified chamber for hybridization. After two times 5-min washes with 2× SSC at 42°C, slides were directly mounted in Vectashield medium containing DAPI (H-1200, Vector Laboratories). Probes for major and minor satellites were amplified from Daxx f/f genomic DNA by PCR (Table EV3) and labeled with Cy3 and Cy5, respectively, using the random priming procedure (Invitrogen, Ref 18095-011).

RNA FISH
Wash buffer A and hybridization buffer containing probes were freshly prepared before performing RNA FISH. For wash buffer A, 2 ml Stellaris RNA FISH Wash Buffer A (Biosearch Technologies Cat# SMF-WA1-60), 7 ml nuclease-free water, and 1 ml deionized formamide were mixed together with gentle vortexing. Hybridization buffer contained 450 ll Stellaris RNA FISH Hybridization Buffer (Biosearch Technologies Cat# SMF-HB1-10) and 50 ll deionized formamide. A mixture of mouse major satellite RNA probes labeled directly by Quasar 570 (forward) and Quasar 670 (reverse) (Biosearch Technologies) was added to 500 ll hybridization buffer (final working concentration of each probe was 125 nM). The zona pellucida of zygotes at G2 phase was removed by protease before fixation (Sigma, P5147-1 G, 5 mg/ml in M2 medium) for 10 min at 37°C prior to fixation. Embryos were fixed in 4% paraformaldehyde in PBS (pH 7.2) for 15 min at RT and permeabilized for another 5 min in 0.5% Triton X-100/ PBS (RNase-free). The embryos were washed in PBS twice and deposited on glass slides to allow adherence. Following two times 5-min washes in wash buffer A, 50 ll hybridization buffer containing probe was added to embryos and they were incubated in a humidified chamber overnight at 37°C. Following aspiration of hybridization buffer embryos were washed three times 10 min in wash buffer A at 37°C. After washing two times with wash buffer B (Biosearch Technologies Cat# SMF-WB1-20) at room temperature, embryos were directly mounted in Vectashield medium containing DAPI (H-1200, Vector Laboratories). Each individual transcription signal was detected and quantified in 3D by the ForciPicker algorithm from ImageJ software (Liu et al, 2012b).

mRNA microinjection and live-cell imaging
For microinjection, constructs were linearized and transcribed in vitro using the mMessage mMachine T7/T3 kit (Ambion). After purification using the RNeasy Mini Kit (Qiagen), the synthesized mRNA was diluted using nuclease-free water (Ambion). To image of 1st mitosis of embryos, we microinjected 3-5 pl in vitro-transcribed RNAs encoding H2B-mCherry (20 ng/ll) and Tubulin-EGFP (100 ng/lg) into zygotes as described above and cultured in incubator for at least 2 h to allow maturation of fluorescent protein. We started live confocal imaging around 26 hp-hCG. To image zygotes from PN0 to PN5 stage, we microinjected 3-5 picoliter mRNAs encoding H3.3-mCherry and H3.2-EGFP or DAXX-EGFP (20 ng/ll each) into MII oocytes. Following 3 h in vitro culture, MII oocytes were fertilized by ICSI and processed for live imaging immediately. For imaging, 10-20 embryos were cultured in 50 ll M16 medium without phenol red (Millipore) covered by mineral oil in a multi-well l-slide (Ibidi, 81506) designed for high end live-cell microscopy. We imaged embryonic development with a spinning-disk multipoint confocal (Olympus IX81) equipped with the Yokogawa CSU-X1 scanhead, 491-and 561-nm laser, and a 40× silicon immersion objective lens (Olympus, UplanSApo, N.A. 1.25). The Ibidi slides were enclosed in a temperature and humidity controlled environment set at 37°C and 5% CO 2 (Life Cell imaging, Basel). Fifteen z-axis confocal sections (every 2 lm) of 512 × 512 pixels images were acquired every 10 min. for monitoring the 1st mitosis and every 20 min to monitor development from PN0 to PN5. We reconstructed the images in 3D using ImageJ.

3D image processing and quantitative analysis
To quantify the level of fluorescence signals within distinct nuclear compartments such as pericentromeric heterochromatin (PCH) regions in mouse early embryos, we developed a new image processing approach that enables automated identification of compartments in paternal and maternal pronuclei and quantification of different fluorescence signals representing, e.g., various chromatin proteins and modifications in wild-type and genetically modified embryos. The method couples supervised learning-based segmentation to an automated iterative rule-based watershed merger that corrects for over-segmentation produced in classical watershed algorithms. The separation of parental pronuclei and identification of their parental origin (maternal or paternal) is based on measured features, such as pronuclear size and intensity distribution of the fluorescent DNA dye DAPI. In detail, raw 4D microscope data (xyz and color channel) are converted into a new 3D data stack with isotropic resolution along xyz using an ImageJ-based macro, working recursively along experiment folders. The macro also extracts and saves separately the DAPI channel for subsequent pixel classification. Next, supervised learning 3D pixel classification is used to generate prediction maps for the different nuclear compartments using several key features, such as the intensity of the DAPI fluorophore, the local DAPI gradients, and the structure tensor. Training is realized with the interactive learning and classification toolkit Ilastik (Sommer & Gerlich, 2013). Finally, each prediction map is morphologically filtered and binarized to define masks for separate nuclear compartments. Quantification of different fluorophore channels relative to the distinct nuclear compartments including PCH is realized by summing up the intensities over the corresponding masks. Such volume normalized read outs are referred to as signal intensities in boxplot figures, which were generated using the online "http://shiny.chemgrid.org/boxplotr/" tool.

Immunoprecipitation and immunoblotting
Cells were transiently transfected with the appropriated constructs and harvested 24 h later. Lysate preparation and subsequent immunoprecipitation were performed as described previously with minor changes (Tardat et al, 2015). Briefly, the lysates were cleared by centrifugation before pre-clearing with 50 ll of mouse IgG agarose beads (Sigma) for at least 1 h at 4°C with rotation. After centrifugation, 25 ll of mouse anti-Flag M2 agarose beads (Sigma) were added to the supernatant and allowed to precipitate overnight at 4°C with rotation. Beads were pelleted by centrifugation, washed three times, and resuspended in SDS buffer. For the analysis of immunoprecipitation, the proteins were separated by SDS-PAGE and detected with the following antibodies: anti-Flag (Sigma, Cat#F1804), anti-GFP (Roche, Cat#11814460001), anti-Histone H3.3 (Millipore, Cat#09-838), anti-H3K27me3 (Cell Signaling, Cat#9733), anti-Myc Tag (Abcam, Cat#Ab9132), and anti-RNF2 (Active Motif, Cat#39663). For immunoblots, anti-aTubulin (Sigma, Cat#T9026) was used as a loading control.

SUMOylation assay
HEK293 cells were co-transfected with expression vectors of either 6xHis-tagged SUMO1, SUMO2, or SUMO3 with pcDNA3.1-Cbx2-EGFP wild-type or mutant as indicated in the corresponding figures. After 24 h, cells were harvested followed by protein extraction under denaturing conditions (6 M guanidinium-HCl, 0.5% (v/v) Triton, 10 mM imidazole, 20 mM Tris at pH 7.5, 0.5 mM DTT, and 0.5 mM Iodoacetamide). To prevent the loss of SUMO conjugation, 10 lM N-Ethylmaleimide (NEM) was added to lysis and washing buffers. His-SUMO-conjugated proteins were purified by nickel chromatography (Qiagen) and separated by SDS-PAGE. Immunoblot was performed with the an anti-GFP antibody (Roche, Cat#11814460001).
Chromatin immunoprecipitation, ChIP-Sequencing, and quantitative PCR HEK293 5xUAS cells were transfected twice with Lipofectamine 3000 with the appropriate constructs and fixed 24 h later with 1% (v/v) formaldehyde (Sigma) for 10 min at room temperature before quenching of the cross-linker with 125 mM glycine. Preparation of chromatin was performed as described previously (Tardat et al, 2015). Chromatin was sheared with a Bioruptor apparatus (Diagenode) to obtain DNA fragments with size between 200 and 800 bp. The lysates were incubated overnight at 4°C with Dynabeads protein G (Invitrogen) coupled to the appropriate antibody: GAL4 DNAbinding domain (Millipore, Cat#06-262), Myc-tag (Abcam, Cat#Ab9132), RNF2 (Active Motif, Cat#39663), Histone variant H3.3 (Millipore, 09-838), and Histone H3 (Abcam, Cat#1791). Recovered DNA was purified using the MinElute PCR purification kit (Qiagen Cat#28006) following the manufacturer's instructions. For Sequencing, the libraries were prepared with the ChIP-Seq NEB Ultra Kit (New England Biolabs Cat# E6200) and sequenced with the Illumina Hiseq (50 cycles and single-end reads). Quantitative PCR analyses of immunoprecipitated DNA were performed in triplicate using a Fast SYBR Green Master Mix (Thermo Fisher Scientific Cat#4309155) on an ABI 7500 Fast Real-Time PCR System (Applied Biosystems). The primer pairs used are indicated in Table EV3.
Potential PCR duplicates were removed from libraries which were generated using SoLo NuGEN RNA-Seq protocol by merging BAM files for technical replicates and deduplicating resulting BAM files using UMIs and python script nudup.py from NuGEN Technologies (version 2.2).
Gene Ontology analysis was performed using R package topGO with parameters "method = weight01; statistic = fisher". GO annotation for genes was taken from Bioconductor annotation package org.Mm.eg.db.
Differential expression analysis of RNA-Seq data for Setdb1 KO GV oocytes Analysis was done as described previously (Eymery et al, 2016). Briefly, differentially expressed genes were identified by fitting a two-factor model of the form "cre + genotype", with cre corresponding to expression status of the Zp3 cre-recombinase ("expressed" or "not expressed"), and genotype to the Setdb1 genotype ("+/À" or "À/À"). Statistical significance of differential expression was estimated using log-likelihood tests, and Benjamini-Hochberg method was used to correct for multiple testing.

of 24
The EMBO Journal 39: e103697 | 2020 ª 2020 The Authors Reads for ChIP and input samples were counted for promoters (AE 2.5 kb around TSS) using QuasR package selecting only uniquely mapped reads, and enrichments were calculated by subtracting corresponding input log 2 (RPKM) from ChIP log 2 (RPKM). Classification of promoters into "Low", "Medium", and "High" classes with respect to H3K27me3 enrichments for each stage was done by fitting Gaussian mixture model with 3 components using mixtools R package (Benaglia et al, 2009). Promoters were classified as "Low H3K27me3" if posterior probability for belonging to the lowest Gaussian component was above 0.05, and the rest were classified as "Medium H3K27me3" or "High H3K27me3" choosing class with the highest posterior probability.
MycHis and GAL4DBD ChIP-Seq data processing, analysis, and visualization MycHis and GAL4DBD ChIP-seq reads were aligned to the genome (BSgenome.Hsapiens.UCSC.hg19) using the QuasR::qAlign function; only uniquely mapping reads were considered. ChIP occupancy levels were quantified using the QuasR::qCount function at all mappable TSS AE 1 kb. Counts were normalized by library size, using the number of reads mapping to the whole genome, excluding TSS regions (TSS AE 5 kb) as a normalizer. Enrichment values were obtained by subtracting corresponding input values from the same regions (in log 2 scale). Genome browser tracks were generated in R using the Gviz package.
To understand the anti-MYC enrichment at TSS regions in the absence of exogenous DAXX MycHis -protein expression, we analyzed ChIP-seq data of transcription factors in multiple human cell lines, available at the UCSC browser as part of the "ENCODE Transcription Factor Binding Tracks" ("ENC TF Binding") (ENCODE Project Consortium, 2012).

Quantification and statistical analysis
The number (n) of independent experimental replications, the definition of center and precisions measures are reported in the figure legends. Statistical analysis was performed using SigmaStat 3.5 statistical software. The two-tailed t-test or paired t-test was used for statistical analysis except when specified in the figure legends. For all statistical analyses, a value of P < 0.05 was statistically significant. For box plots, whiskers extend to datapoints that are < 1.5 × interquartile range away from the 1 st /3 rd quartile.

Data availability
The datasets produced in this study are available in the following databases: Expanded View for this article is available online.