The PRC2.1 Subcomplex Opposes G1 Progression through Regulation of CCND1 and CCND2

Progression through the G1 phase of the cell cycle is the most highly regulated step in cellular division. We employed a chemogenomics approach to discover novel cellular networks that regulate cell cycle progression. This approach uncovered functional clusters of genes that altered sensitivity of cells to inhibitors of the G1/S transition. Mutation of components of the Polycomb Repressor Complex 2 rescued growth inhibition caused by the CDK4/6 inhibitor palbociclib, but not to inhibitors of S phase or mitosis. In addition to its core catalytic subunits, mutation of the PRC2.1 accessory protein MTF2, but not the PRC2.2 protein JARID2, rendered cells resistant to palbociclib treatment. We found that PRC2.1 (MTF2), but not PRC2.2 (JARID2), was critical for promoting H3K27me3 deposition at CpG islands genome-wide and in promoters. This included the CpG islands in the promoter of the CDK4/6 cyclins CCND1 and CCND2, and loss of MTF2 lead to upregulation of both CCND1 and CCND2. Our results demonstrate a role for PRC2.1, but not PRC2.2, in promoting G1 progression.


Introduction
Cellular decisions to grow and divide are made by assessing the balance of activating and inhibitory inputs that govern the transition between cell cycle phases.Regulated progression through the cell cycle is crucial for normal cellular growth and organismal development [1][2][3] .Progression from G1 into S phase is the most highly regulated step of the cell cycle, as initiating DNA replication commits a cell to divide and is frequently mutationally activated in tumors.Cyclin Dependent Kinase 4 (CDK4) and the related CDK6 (henceforth referred to collectively as CDK4/6) play critical roles in promoting G1 progression through phosphorylation of the retinoblastoma protein (RB1).Phosphorylation relieves RB1-mediated transcriptional repression of E2F transcription factors, which are then competent to drive transcription of genes necessary for progression into S phase [3][4][5][6] .Because of their crucial role in regulating G1 progression, specific inhibitors targeting CDK4/6 have proven to be effective therapeutics.Palbociclib was the first FDA approved CDK4/6 inhibitor and highly efficacious in the treatment of HR+/HER2-breast cancers, followed by the structurally related molecules ribociclib and abemaciclib [7][8][9][10][11] .However, this classical model of G1 regulation has recently given way to a more complex model [12][13][14] , underscored by the complexity of genetic alterations that lead to resistance to treatment with CDK4/6 inhibitors 4,12,13,15 .Thus, while G1 progression has been the focus of intense study, our understanding of its regulation remains incomplete.
The Polycomb Repressive Complex 2 (PRC2) was initially identified in Drosophila as a developmental regulator that represses the expression of Hox genes 16 .The PRC2 complex is conserved in throughout eukaryotes [17][18][19][20][21] and catalyzes the mono-, di-and tri-methylation of Histone 3 Lysine 27 (referred to collectively as H3K27me3, the fully methylated form of H3K27) and thereby acts as a transcriptional repressor 20 .The core PRC2 complex is composed of a H3K27me3 "reader" EED, a scaffold protein SUZ12, and the catalytic subunit EZH2 (or the more poorly-expressed and less catalytically active paralog EZH1 22 ).This core complex is capable of catalyzing H3K27me3 deposition and chromatin association, but how PRC2 achieves full spatiotemporal regulation of chromatin localization and catalytic activity has been an area of active investigation.Recent studies have identified additional accessory factors that modify the localization and enzymatic activity of these core components 23 .The associated auxiliary factors define different PRC2 subcomplexes, which are called PRC2.1 and PRC2.2, based on the composition of the subunits associated with the core PRC2 complex (reviewed in 24,25 ).In addition to the core PRC2 subunits, PRC2 .1 consists of two modules, one module containing a Polycomb-like (PCL) protein PHF1, MTF2 or PHF19 and a second module of either PALI1/2 or EPOP.The more homogenous PRC2.2 always consists of the core PRC2 subunits in complex with both JARID2 and AEBP2.The role of these complexes in different cellular processes and contexts is debated.Despite the lack of an a clear analogous sequence to the Polycomb Response Elements which promotes PRC2 chromatin association in Drosophila 26,27 , the presence of a DNA-binding extended homology domain in each PCL protein has been proposed to recruit PRC2.1 to unmethylated CpG islands and establish H3K27me3 28,29 .In contrast, PRC2.2 localizes to sites utilizing pre-existing mono-ubiquitinated H2AK119 (H2AK119ub1), which is deposited by the PRC1 complex [30][31][32][33] , through a ubiquitin interaction motif contained within JARID2 31,32,34,35 .Regardless of their specific roles in the propagation of H3K27me3 histone marks, members of both PRC2.1 and PRC2.2 have been implicated as both positive and negative regulators of stem cell maintenance, differentiation and cancer, depending on the cellular context 25,[36][37][38][39][40][41][42] .All of the PRC2 core subunits (EZH2, SUZ12 and EED) have been shown to inhibit that transcription of both positive and negative regulators of G1/S progression, including the CDK4/6 protein inhibitor p16 [43][44][45][46][47] .However, the net result of these opposing effects on cell cycle progression, and the contribution of the individual subcomplexes to this regulation, rained unclear.
To identify novel regulators of cellular proliferation, we utilized a whole-genome chemogenetic approach to identify genes that sensitized or lent resistance to inhibitors of different cell cycle stages.We uncovered novel resistance mechanisms to three known inhibitors of cell cycle progression in the human haploid cell line HAP1.This approach revealed that mutations in mitochondrial function or the Polycomb complexes rescued slow growth in palbociclib.We could recapitulate these positive genetic interactions pharmacologically using small molecule inhibitors of either PRC2 activity or mitochondrial respiration.Loss of core PRC2 members or PCL subunits of the PRC2.1 subcomplex, particularly MTF2, resulted in resistance to palbociclib, while loss of PRC2.2-specific subunits resulted in sensitivity.Data from CUT&RUN and RNA sequencing experiments performed on clonal MTF2∆ and JARID2∆ knockout mutant cell lines suggest that PRC2.1 plays a more critical role in repressing gene expression when compared with PRC2.2 in HAP1 cells, particularly at promoters containing CpG islands.D-type cyclins are among the genes that are repressed by PRC2.1 and loss of MTF2 results in increased expression of both CCND1 and CCND2 through loss of H3K27me3 in their promoters.This increased expression resulted in an apparent increase in CDK4/6 kinase activity and S-phase entry of cells, driving resistance to CDK4/6 inhibition.Our results suggest that PRC2.1 plays a strong role in G1 progression.

Chemogenetic CRISPR-Cas9 Screen Utilizing Cell Cycle Inhibitors Identified Novel Players in the Cell Cycle
Recently, CRISPR-Cas9 knockout genetic screens have emerged as a powerful way in which to probe genetic interactions [48][49][50] , with the haploid human cell line HAP1serving a popular model for these studies [51][52][53][54][55][56] .To identify novel genes involved in cell cycle regulation, we carried out genome-wide CRISPR-Cas9 chemogenomic screens in HAP1 cells treated with each of three well-characterized inhibitors of cell cycle progression: palbociclib (a CDK4/6 and G1 progression inhibitor), colchicine (a microtubule polymerization and mitosis inhibitor) and camptothecin (a Topoisomerase I and S/G2 inhibitor).We used a concentration for each inhibitor that reduced cellular proliferation by 30-50% (Fig S1A ) (see also 51 ) and confirmed their effects on cell cycle progression (Fig. S1B).We then performed a CRISPR-Cas9 whole-genome screen for each of the three inhibitors (Fig. 1A) by introducing the Toronto Knockout Library 57 via lentiviral transduction into a HAP1 cell line constitutively expressing Cas9.Following puromycin selection for two days, cells were propagated in either DMSO (Mock) or in the presence of drug (Treated) for eighteen days.Following propagation, genomic DNA was extracted from the initial and final pools and subjected to deep sequencing, and gene-compound interactions were determined using the Orobas pipeline (source code -Supplemental File 1, chemogenomics screen results -Table S1 and S2).A gene was considered as being significantly enriched or de-enriched in a given condition if both the Loess-adjusted differential gene effect between mean Treated and Mock control was ±0.5 and the FDR<0.4.This approach resulted in the recovery of predicted compound-gene interactions demonstrating the robustness of both the screen and our analysis approach.For example, targeting of genes known to play roles in DNA Damage Repair (DDR) 58,59 , including RAD54L, MUS81 and sixteen proteins in the Fanconi Anemia pathway, strongly sensitized cells to camptothecin, which generates protein-DNA adducts (Fig. 1B and 1C).The molecular target of camptothecin, TOP1 60 , is the most resistant gene in the camptothecin screen, as are proteins involved in p53 transcriptional regulation, such as TP53BP1 and STAGA members TAF2, TAF4, TAF5, TAF11 and TAF13 (Fig. 1B).Similarly, colchicine sensitized cells to the mutation of genes encoding proteins involved in mitotic spindle assembly, nuclear division and cytoskeletal assembly, such as CLASP1, DLGAP5 and KNTC1 (Fig. 1D and 1E).Interestingly, inactivation of genes involved in the adaptive immune system, such as BIRC6, UBA5 and USP14, also resulted in sensitivity to colchicine.This observation is intriguing, as colchicine is used clinically as an immunomodulator in the treatment of gout 61 .CCNE1, CDK6, CDK2, CCND2 and CCND1, all of which are integral to promoting the G1/S phase transition, ranked as the 2 nd , 24 th , 27 th , 29 th and 46 th most important genes for palbociclib resistance, respectively (Fig. 1F and 1G).CCND1 and CCND2 bind either CDK4 or CDK6, the molecular targets of palbociclib, whereas CDK2 and CCNE1 form a related CDK kinase that promotes the G1/S transition.Similarly, cells with sgRNAs targeting RB1, whose phosphorylation by CDK4/6 is a critical step in G1 progression, displayed substantial resistance to palbociclib.The recovery of genes known to function in the relevant biological processes supports the strength of this dataset and bolstered our confidence to use the results obtained to identify novel chemicalgenetic interactions.

Chemogenomics Screen Uncovered Novel Genetic Interactions Involved in Response to Inhibitors of Cell Cycle Progression
To identify genes whose inactivation rendered cells sensitive or resistant to a specific cell cycle perturbation, we compared how the Orobas-calculated differential gene effect for a given targeted gene varied in each compound across our CRISPR-Cas9 screen.The majority of genes that conferred either sensitivity or resistance were specific to only one cell cycle inhibitor, with little overlap between the multiple conditions, suggesting that we identified genes that play roles in distinct biological processes (Fig. 2A and Fig S1C ).For example, genes encoding DNA repair proteins, mitotic spindle components, and CDK2/4/6 holoenzyme components were only required for proliferation in camptothecin, colchicine, and palbociclib, respectively.We found that only thirteen and twenty genes resulted in sensitivity or resistance, respectively, in every conditions tested and were deemed non-specific and excluded from any further analysis (see Table S2).
We next turned our attention to unexpected and novel compound-gene interactions.To probe these interactions, we analyzed genes that significantly altered response to our three compounds using the gene annotation and analysis portal Metascape and the protein-protein interaction network analysis STRING.In addition to DDR genes, Metascape and STRING analysis of the results of our camptothecin treatment revealed de-enrichment for sgRNAs targeting genes encoding members of the KICSTOR complex (KPTN, SZT2, ITFG2 and KICS2), which negatively regulates mTOR.In contrast, sgRNAs targeting of genes involved in RNA metabolism and chromatin organization increased resistance to this drug (Fig. 2B-2D).It has been suggested that mTOR is involved in attenuating the DDR response through phosphorylation of RNF168, leading to its degradation 62 , which could provide one explanation of the observed sensitivity.The loss of genes involved in chromatin structure and the metabolism of RNA conferred resistance to camptothecin is unexpected, given that both these processes have been implicated in DNA repair after damage 63,64 .Genes whose inactivation enhanced sensitivity to colchicine included those involved in the amino acid starvation response (DEPDC5 TSC1, SZT2 and NPRL2) and mRNA splicing (SNRPB2, SF3B2, PPIL1, RBM22 and DHX35), while mutation of genes that control vesicle trafficking (VPS16, VPS18, VPS29, VPS41, VPS51 and VPS52) or encode members of the Mediator complex (CCNC, CDK8 and MED26, MED1, MED7, MED12, MED18 and MED11) attenuated the cellular response to the drug.Unexpectedly, inactivation of genes encoding members of the TP53 signaling pathway (TFDP1 and HIPK2), SAGA H3 acetylation complex components (KAT2A, TRRAP, TADA3, TAF5L, TADA2B, SGF29 and TADA1), and the Fanconi Anemia complex (FANCA, FANCC, FANCE, FANCF, FANCG, FANCL and FAAP100), all implicated in DNA damage repair, resulted in resistance to colchicine (Fig. 2E-2G).Sensitivity to palbociclib was enhanced in cells expressing sgRNAs targeting H4 acetylation, positive regulators of Pol II transcription and regulators of the DNA Damage Response pathway (Fig. 3A and 3B), although this sensitivity was much weaker than that seen with DNA damaging agents.This observation is consistent with long-term treatment with palbociclib inducing DNA damage, as has been suggested by a number of recent publications 65,66 .Unexpectedly, Metascape analysis of our palbociclib chemogenetic screen revealed that sensitivity to palbociclib was enhanced when genes involved in chromatin organization were targeted (Fig. 3A).Inactivation of members of the SIN3 histone deacetylase (SIN3B, SINHCAF and ARID4B), the NuA4 histone acetyltransferase (ING3, DMAP1, MORF4L2, YEATS4 and VPS73), the STAGA histone acetyltransferase (KAT2A,TADA1, TADA2B, TAF5L, and SUPT20H) and the Mediator (MED13, MED25, MED10, MED15, TAF7, TAG13 and CCNC) complexes all resulted in palbociclib sensitivity.(Fig. 3B).Each of these protein complexes promotes gene expression, suggesting that palbociclib sensitivity might be a result of a reduction in the transcription of genes important for the G1/S transition.
Because mechanisms of clinical resistance to palbociclib is an area of active investigation, we turned our attention to focus on these pathways.Metascape analysis of genes whose loss conferred palbociclib resistance were highly enriched for splicing factors, oxidative phosphorylation and mitochondrial translation, in addition to chromatin modification (Fig. 3A).STRING analysis of the high-confidence, physical interactions of proteins important for palbociclib sensitivity revealed multiple highly-connected interaction networks (Fig. 3B).Strikingly, almost 25% (170 out of the 689) of the genes whose mutation conferred unique resistance to palbociclib have terms associated with mitochondrial respiratory chain complex assembly, ATP synthesis or mitochondrial gene expression.Specifically, we see many components implicated in assembly of Mitochondrial Respiratory Chain Complex I and IV, as well as core mitochondrial ribosome and mitochondrial translation initiation and termination (Fig. 3A-3C).To confirm this positive genetic interaction between mitochondrial homeostasis and resistance to palbociclib, and to dissect whether specific electron transport chain steps might be implicated in this resistance, we asked whether chemical inhibition of oxidative phosphorylation could rescue sensitivity to palbociclib.To target different stages of the oxidative phosphorylation, we employed Rotenone, TTFA and Oligomycin, which inhibit Complex I, Complex II and ATP synthase, respectively.Cells were grown in the presence of palbociclib alone or in combination with each drug for 48 hours and viability was determined by PrestoBlue assay.Cells exposed to Rotenone, TTFA and Oligomycin all showed positive, suppressive growth interactions with palbociclib (Fig. 3D).While alternative explanations could explain the observed novel chemical-genetic interactions we uncovered here, such as either changes in phenotypic lag rates due to alterations in protein stability or more general screen variability 67 , these results suggest a connection between mitochondrial gene function and CDK4/6 inhibitors.

Polycomb Repressive Complex Components Display Differing Responses to Palbociclib Treatment
Intriguingly, inactivation of EZH2, SUZ12 and EED, the three core members of the PRC2 complex, resulted in profound resistance to palbociclib, being the 8 th , 4 th and 3 rd strongest resistance hits out of the 18,053 genes examined when ranked by the score of differential effect (Fig. 4A, Table S1).Mutation of RBBP7, which associates with the core PRC2 complex 20 along with a number of histone deacetylases 68 , also desensitized cells to palbociclib, but to a more modest extent.Satisfyingly, inactivation of RING1, RNF2, and PCGF6, which are members of PRC1, also displayed resistance to palbociclib.The PRC1 complex contains a ubiquitin ligase that works in concert with PRC2 through H2AK119ub1 deposition, a histone mark that influences both PRC2 chromatin localization and catalytic activity 30,32 .As expected, PRC1 and PRC2 components identified in our palbociclib chemogenomics screen formed a highly interconnected STRING physical interaction network (Fig. 4B), indicating that loss of either H3K27me3 or H2AK119ub1 reduced sensitivity to this drug.In contrast, inactivation of genes encoding OGT, ASXL1 and HAT1, which are members of the H2AK119ub1 deubiquitinase complex that opposes PRC2-mediated gene repression 69 , resulted in sensitivity to palbociclib (Fig. 4A).Importantly, no component of any PRC1 or PRC2 core complex displayed significant resistance or sensitivity to camptothecin and colchicine in our chemogenomics screens (Fig. 4C), implicating PRC2 in the regulation of G1 specifically, and not to other phases of the cell cycle or the DNA damage response pathway.PR-DUB components ASXL1 and OGT did show resistance to camptothecin, consistent with their role in repressing the homologous recombination DNA repair pathway 70 .We sought to confirm the role of the core PRC2 complex in palbociclib resistance by treating cells with combinations of palbcociclib and the EZH2 inhibitor GSK126 using a quantitative Crystal Violet assay.After nine days of drug combination treatments, we found that cells treated with increasing doses of GSK126 withstood palbociclib-induced growth suppression (Fig. 4D), confirming that inactivation of the PRC2 core complex, either through genetic inactivation or chemical inhibition, resulted in resistance to palbociclib.
The PRC2 core binds to auxiliary proteins to create biochemically distinct subcomplexes, termed PRC2.1 and PRC2.2 23,24 .These alternative complexes are thought to modify the chromatin localization and enzymatic activity of PRC2, reenforcing existing H3K27me3 in certain contexts [71][72][73] , while initiating H3K27me3 deposition at new loci in others 33,39 .Mutation of the PRC2.1 complex members PHF1, MTF2, PHF19 and EPOP/C17orf96 all display significant resistance to palbociclib, with MTF2 being the strongest of these (Fig. 4A, 4C).Conversely, targeting the genes encoding the PRC2.2-specificaccessory proteins AEBP2 or JARID2 resulted in enhanced or neutral palbociclib sensitivity, respectively.These data suggest that PRC2.1 plays a previously uncharacterized role in promoting G1 progression, while PRC2.2 antagonizes it.To confirm the results from our palbociclib chemogenomics screen, we generated polyclonal knockout mutant pools of the individual core and accessory proteins of PRC2.We generated these populations by independently infecting three distinct sgRNAs targeting genes for each PRC2 complex member, or positive and negative control genes, in GFP-positive and doxycycline-inducible Cas9 cells and induced DNA cleavage for three days (henceforth referred to as pooled knockouts).We performed Western blots to confirm reduction in protein levels for the genes targeted by the sgRNAs used to generate our pooled knockouts (Fig. S2B).After confirming reduction in the targeted proteins, we carried out a competitive growth assay using these pooled knockouts (schematic in Fig. S2A or see 51 ).Briefly, GFP-positive pooled knockouts were mixed with GFP-negative wild-type cells and propagated in the presence or absence of palbociclib.The ratio of GFP-positive to GFP-negative cells was recorded every three days by flow cytometry for eighteen days.Pools containing sgRNAs against EZH2, SUZ12, EED and MTF2 all displayed resistance to palbociclib, similar to the level observed with sgRNAs targeting RB1, our positive control for palbociclib-induced growth defects (Fig. 4E).Conversely, compared with the dramatic reduction seen in pools transduced with CCNE1 sgRNAs, our positive control for enhanced palbociclib sensitivity CCNE1, sgRNAs targeting PRC2.2 components showed a slight, but statistically signification reduction in proliferation in palbociclib over the eighteen day assay (AEBP2: p-value = 0.002 and JARID2: p-value = 0.0148, unpaired two-tailed Student's t-test).Thus, we confirmed the results of our chemogenomics screen that MTF2-containing PRC2.1 inhibits G1 progression, while PRC2.2 does not, and could in fact promote it.
To further interrogate the role of PRC2.1 and PRC2.2 in the regulation of G1 progression, we generated SUZ12, MTF2 and JARID2-null monoclonal cell lines (which we will refer to as SUZ12∆, MTF2∆ and JARID2∆).We chose to mutate MTF2 to probe the function of PRC2.1, as it has been shown to be more highly expressed and the predominant PCL subunit associated with the PRC2 core complex in a variety of contexts 74,75 .Furthermore, we selected SUZ12 for inactivation out of the core PRC2 complex members, and not the catalytic subunit EZH2, because the presence of the EZH2 paralogue EZH1 might compensate for loss of EZH2 76 .Additionally, SUZ12 has a critical role in bridging accessory proteins with the catalytic core in all known PRC2 complexes 39 .MTF2∆ cells displayed resistance to palbociclib when compared with wild-type cells in a nine-day quantitative Crystal Violet assay (MTF2∆ IC 50 = 1.033µM , wild-type IC 50 = 0.3936µM) while JARID2∆ cells were slightly more sensitive than wild-type (JARID2∆ IC 50 = 0.2216µM) (Fig. 4F).In addition to showing sensitivity to palbociclib, MTF2∆ cells also displayed resistance to, ribociclib and abemaciclib, two CDK4/6 inhibitors that are structurally related to palbociclib, in a competitive growth assay (Fig. S2C).These results confirmed our screen results that mutation of MTF2 leads in CDK4/6 inhibitor resistance with verified clonal mutants.
Palbociclib exerts its chemotherapeutic effects by inducing a G1 arrest and senescence in tumor cells with a functional RB-E2F pathway [77][78][79] .However, a recent report demonstrates that palbociclib treatment induces both G1 arrest and apoptosis through the increase in DNA damage in cultured cells 80 , introducing the possibility that PRC2.1 could be altering regulators of the DDR pathway, resulting in the observed palbociclib resistance.To determine if inactivation of PRC2.1 or PRC2.2 altered cell cycle progression, we examined how wild-type, SUZ12∆, MTF2∆ and JARID2∆ cells responded to palbociclib-induced G1 arrest.To assess this directly, we performed a BrdU incorporation assay by growing each mutant for 24 hours in palbociclib, pulsed the cells with BrdU for one hour prior to harvest and then measured BrdU incorporation by flow cytometry.Each of the four cell lines had similar levels of BrdU incorporation in the absence of drug (Fig. 4G).However, MTF2∆ and SUZ12∆ mutants showed significantly more cells than wild-type in S-phase in the presence of palbociclib (Fig. 4G).JARID2∆ mutants showed fewer cells in S phase, but our assay was unable to establish that this reduction was statistically significant.To rule out the possibility that cellular viability was not compromised in our monoclonal knockout cell lines, we used Western blotting to monitor changes in PARP cleavage or increased BCL2L11/BIM expression, which both serve as apoptosis indicators 81 .There was no detectible basal increase in apoptosis in the monoclonal knockout mutant cell lines or when cells were treated with palbociclib for 48 hours, (Fig. S2D), supporting the conclusion that the resistance to palbociclib observed in the MTF2∆ and SUZ12∆ cells was due to the repressive role the MTF2-containing PRC2.1 complex plays in the canonical CDK4/6-RB1-E2F pathway.

PRC2.1 and PRC2.2 Mutants Display Altered H3K27me3 and Transcriptional Landscapes
To determine why the mutation of PRC2 subcomplex components altered the cellular response to palbociclib, we sought to see how H3K27me3 levels and gene expression changed in MTF2∆ and JARID2∆ cells.Western blotting of total H3K27me3 levels in three independently-generated clones indicated that there was no change in the bulk levels of H3K27me3 (Fig. 5A), suggesting that any change of phenotype observed in the mutants was due to a change in the localization of this mark and not due to an overall reduction in its abundance.This is in contrast to SUZ12∆ cells, which displayed a significant reduction in the H3K27me3 mark (unpaired Student's t-test, p-value = 0.0104).To probe the changes in transcription and H3K27me3 distribution genome-wide, we generated CUT&RUN libraries with an anti-H3K27me3 antibody and RNA-Seq libraries from total RNA isolated from our MTF2∆ and JARID2∆ cell lines, grown either in the presence or absence of palbociclib for 24 hours.Changes in H3K27me3 levels and mRNA expression were determined by comparing the enrichment of reads in the MTF2∆ and JARID2∆ libraries to the wild-type cell line (Tables S3 and S4).Because cancer cells have been known to adapt to palbociclib treatment partially through changes to histone marks, chromatin structure and gene expression [82][83][84][85] , we also investigated how both transcript levels and H3K27me3 distribution responded to treatment with palbociclib in our clonal knockout cell lines.Primary component analysis (PCA) of our called, reproducible H3K27me3 peaks and transcript abundance from our CUT&RUN and RNA-seq data, respectively, showed a high percentage of variance between each of the genotypes tested, along with good clustering of repeats of the same genotype and treatments (Fig. S3A), suggesting a shift in the epigenetic and transcriptional landscapes when either MTF2 or JARID2 are absent.PCA analysis of our RNA-Seq experiment revealed substantial shifts in variance between palbociclib-treated and Mock-treated samples for each genotype (Fig. S3A, bottom), suggesting that exposure to palbociclib resulted in changes in gene expression, consistent with previous reports 86,87 .However, the PCA of our CUT&RUN experiment did not reveal large differences in H3K27me3 distribution between palbociclib-treated and untreated samples (Fig. S3A -top).In line with this observation, when we analyzed the change in distribution of H3K27me3 peaks between palbociclib and Mock treated cells using DESeq2, we found no significant changes in the location of H3K27me3 reproducible peaks in the presence or absence of palbociclib (data not shown).This suggests that MTF2∆ mutants are not resistant to palbociclib because MTF2 is required for a transcriptional adaptation to the drug, but instead because MTF2 alters expression of genes important for G1/S progression, even in unperturbed cells.
Due to the known role of PRC2 in repressing gene expression, we next asked how H3K27me3 distribution changed in promoters of genes.We defined promoters as 4kb upstream and 1kb downstream of all annotated transcription start sites, and calculated the total number of reads within each of these regions.Our parameters led to ~61,000 genomic regions being designated as promoters.In addition to annotated protein coding genes, this included the promoters of non-coding transcribable units such as rRNA, miRNAs, lncRNAs, and pseudogenes.We observed a greater number of promoters with significantly decreased H3K27me3 (log 2 fold-change ±1, adjusted p-value < 0.1) in the MTF2∆ compared to JARID2∆ cell lines (5,808 vs 1,034 promoters, respectively).Of these, 5,149 promoters displayed MTF2-dependent H3K27me3, 392 were JARID2-dependent and 629 were co-dependent on MTF2 and JARID2 for wild-type levels of H3K27me3 (Fig. 5B and 5D).Consistent with the greater change in H3K27me3 signal at promoters in MTF2∆ cells, 733 vs 114 transcripts were significantly upregulated upon MTF2 vs JARID2 inactivation, respectively, with 666 transcripts that were exclusively MTF2-dependent, 47 transcripts that were exclusively JARID2-dependent, and 67 transcripts displaying co-dependence on both MTF2 and JARID2.These results indicate that the MTF2containing PRC2.1 complexes affect the deposition of H3K27me3 in the promoter regions of more genes than the JARID2-containing PRC2.2, and therefore, are more important for transcriptional repression in HAP1 cells.
Given the diverse regulatory roles of PRC2 in different biological contexts, and the limited information on PRC2.1 and PRC2.2 outside of stem cells, we were curious to see what classes of genes were being differentially regulated in the MTF2∆ and JARID2∆ cell lines.Only ~30-40% of the promoters with significantly changed levels of H3K27me3 were upstream of protein coding genes (Fig S3B ), while ~80-90% of the significantly differentially expressed transcripts encoded proteins (Fig. S3C).Therefore, we focused a Metascape analysis on the promoters and mRNAs of protein coding genes with differential H3K27me3 and transcript levels, respectively (Fig. 5C).Analysis of the promoters of genes with decreased H3K27me3 and increased transcript levels in both MTF2∆ and JARID2∆ lines were terms associated with embryonic morphogenesis, cell fate commitment and developmental growth, all processes previously been shown to be regulated, at least in part, by PRC2 88 .Intriguingly, terms for genes that specifically displayed decreased promoters H3K27me3 and upregulated mRNA in MTF2∆ cells included the pro-growth pathways cGMP and ERBB4 signaling.Conversely, terms for genes which displayed increased promoter H3K27me3 signal and decreased transcript levels in JARID2∆ cells contained pathways that could reduce cellular proliferation and viability, such as positive transcriptional regulation of RUNX1 and positive regulators of program cell death.We also saw terms that had opposite effects on H3K27me3 and transcript levels in MTF2∆ compared to JARID2∆ cells, such as secretion by the cell and regulation of cellular component biogenesis, which could potentially exacerbate palbociclib-induced growth defects 89,90 .Together, these data support the notion that MTF2 antagonizes cell growth in normal cellular conditions, while JARID2 promotes it.

PRC2.1 and PRC2.2 Mutants Display Differential H3K27me3 Modification in Promoters in Cell Cycle Related Genes with CpG Islands
MTF2-containing PRC2.1 have been previously shown to localize to chromatin using a winged helix in its extended homology domain that has affinity for CG-rich sequences 28,29 , whereas PRC2.2 localization is dependent on chromatin context, specifically H2AK119ub1 deposited by PRC1 [31][32][33] .To determine whether CpG island targeting by PRC2.1 could help explain the palbociclib resistance we observed in the absence of MTF2, we identified and plotted 1,877 peaks that overlapped with CpG islands in wild-type cells and had the greatest H3K27me3 signal in a 10kb window surrounding the CpG islands.We then plotted the H3K27me3 signal observed in the MTF2∆ and JARID2∆ cells for these same loci (Fig. 5E).We observed a complete loss of H3K27me3 signal intensity at CpG islands in the MTF2∆ mutants, but only a partial loss at these loci in JARID2∆ cells (Fig. 5E and Fig. S3D).When we expanded our findings genome-wide, we found a significant loss of H3K27me3 peaks at CpG islands in MTF2∆ cells (Fisher's exact test, odds ratio = 20.4,p-value 2.2x10 - 308 ), compared with JARID2∆, where this loss was much more modest (fisher's exact test, odds ratio = 9.8, pvalue = 6.5x10 -7 ).This result is consistent with the interpretation that the MTF2-containing PRC2.1 is required for all H3K27me3 deposition at CpG islands, whereas JARID2-containing PRC2.2 is only required to achieve full wild-type H3K27me3 levels at these sites.
CpG islands are a very common feature of mammalian promoters, with 50-70% human promoters estimated to contain at least one CpG island 91 .Since promoters are highly associated with CpG islands, we examined 2,000 promoters with the highest level of H3K27me3 signal intensity that overlapped with CpG islands in wild-type cells, then plotted the H3K27me3 signal intensity at those same loci in our mutant cell lines (Fig. 5F, left plots).Consistent with the result seen at CpG islands genome-wide, we observed a complete loss of high signal intensity in the MTF2∆ cells, but only a slight loss in JARID2∆ cells.When we averaged the H3K27me3 signal intensity over all 25,124 promoters that contain CpG islands, we observed a pattern of MTF2∆ cells having greatly decreased H3K27me3 levels in these regions, particularly surrounding the transcription start site (Fig. 5F, right plot).We also observe a partial loss of H3K27me3 at the 2,000 promoters that contain CpG islands with the greatest intensity of H3K27me3 in the JARID2∆ lines compared to wildtype (Fig 5F, left plots), in line with what was seen at CpG islands genome-wide.However, when we averaged the signal intensity over all 25,124 promoters that contain CpG islands in JARID2∆ cells, the distribution of H3K27me3 in the promoter regions looks similar to the wild-type distribution (Fig. 5F, right plot).Reactome and MSigDB analysis of the promoters of protein coding genes that overlapped with CpG islands showed strong enrichment for terms associated with cell cycle and E2F target genes (Fig. 5G) as well as enrichment binding E2F6 (Fig. S3E, p-value: 1.9x10 -91 ), which both regulates transcription of G1 progression genes 92 and is a well characterize component of Polycomb complexes 23,93 .These results suggest that MTF2 is required for H3K27me3 deposition at promoters containing CpG islands involved in cell cycle regulation and can explain why MTF2∆ cells display a greater change in gene expression than do JARID2∆ cell lines.

PRC2.1 Represses Expression of CCND1 and CCND2
Our CUT&RUN results suggest that MTF2-containing PRC2.1 impacts gene expression, at least in part, through deposition of H3K27me3 at promoters with CpG islands.Therefore, we hypothesized this PRC2 complex must be antagonizing G1 progression through repression of cell cycle-promoting genes.When inspecting the results of our CUT&RUN and RNA-Seq experiments, we found that the promoters of both CCND1 and CCND2 had lost H3K27me3 signal (Fig. 6A) and displayed strong transcriptional induction in MTF2∆ cells (Fig 6B).In fact, while CCND1 and CCND2 were both among the most up-regulated statisticallysignificant transcripts within the MTF2∆ cell line, their transcription and promoter H3K27 methylation were unaltered in JARID2∆ cells (Fig. S4A and S4B).Furthermore, CCND1 and CCND2 displayed both decreased H3K27me3 and increased mRNA levels in MTF2∆ cells (Fig. 6C), suggesting that the increase in these transcripts was due directly to a change in H3K27me3 in their promoters.Given that increased CCND1 levels is sufficient to drive increased CDK4/6 kinase activity, upregulation of these D-type cyclins is likely to be a significant contributor to the palbociclib resistance in MTF2∆ cells.DESeq2 analysis of H3K27me3 density in MTF2∆ cells displayed a statistically significant 4.3 and 2.7 log 2 fold-decrease in H3K27me3 signal in the promoter region of CCND1 and CCND2, respectively, when compared to wild-type H3K27me3 levels (Fig. 6D, Table S3), whereas changes in H3K27me3 levels in the CCND3 promoter were not statistically significant (Fig S4D, Table S3).Given our observation that H3K27me3 signal is lost at CpG islands in MTF2∆ cells, we inspected the D-type cyclin promoters for CpG islands.Indeed, the regions upstream of all three D-type cyclins contained CpG islands, but CCND1 and CCND2 had regions of GC density about seven times larger (7,460bp and 6,003bp, respectively) than CCND3 (996bp).Furthermore, the promoter of CCND1 contained about twice as many CpG repeats than did CCND2 (575 vs 379) and about six times as many CpG repeats as CCND3 (575 vs 95) (Fig. S4D).These results suggest that the levels of CCND1 and CCND2 mRNA transcripts, but not CCND3, were regulated by MTF2 in a CpG island-dependent manner.
We sought to confirm our observation that ablation of MTF2 resulted in increased levels of CCND1 and CCND2 protein.We generated pooled knockouts of MTF2, JARID2 and the core PRC2 components SUZ12, EZH2 and EED using three independent sgRNAs.In pooled knockouts of MTF2, EZH2, EED and SUZ12, we observed an increase in both CCND1 and CCND2 protein levels by Western blot, but not for CCND3 (Fig. 6E).Consistent with the results from our CUT&RUN and RNA-Seq datasets, we did not observe a significant change in either CCND1 or CCND2 levels in JARID2 pooled knockouts.We next examined mRNA and protein levels of the D-type cyclins in MTF2∆ and JARID2∆ clones by qRT-PCR and Western blotting, respectively.Again clonal knockouts of SUZ12∆ and MTF2∆, but not JARID2∆ lines, had increased mRNA (Fig. 6F) and protein levels (Fig. 6G and Fig. S6E) for both CCND1 and CCND2, but not CCND3.To determine whether other genes involved in the canonical CDK4/6-RB-E2F pathway were also altered, we examined mRNA and protein levels of known cell cycle regulators in our knockout cell lines.In contrast to CCND1 and CCND2, none of the E-type cyclins, CIP/KIP CDK inhibitors, RB1 or E2F proteins displayed significantly altered mRNA transcript abundance in our RNA-seq experiment in either MTF2∆ or JARID2∆ lines (Fig. S4F).To confirm that protein stability of these factors was not altered in our knockout lines, we also examined protein levels of a panel of known G1 regulators by Western blot (Fig. S4G).Similarly, we did not observe an increase in levels of any of the tested proteins, confirming that CCND1 and CCND2 were the only up-regulated canonical CDK4/6-RB-E2F pathway regulators in MTF2∆ cells.
While D-type cyclins are necessary to activate the kinase activity of CDK4/6, they have also been shown to play roles outside of the RB1-E2F pathway [87][88][89] .We sought to test if CDK4/6 activity was increased in MTF2∆.To do this, we determined the extent of RB1 phosphorylation at S807/S811, which are wellcharacterized CDK4/6 targeted residues.To do this, we titrated wild-type, MTF2∆, and JARID2∆ cells with increasing amounts of palbociclib and determined the levels of total RB1 and phosphorylated RB1 levels to calculate the ratio at each concentration.In each of our cell lines, higher concentrations of palbociclib resulted in decreased levels of phosphorylated RB1, as expected.However, compared to WT or JARID2∆ cells, MTF2∆ mutant cells maintained a higher ratio of phosphorylated to unphosphorylated RB1 at each concentration of palbociclib tested (Fig. 6H).This result suggests that the increased levels of CCND1 and CCND2 in MTF2∆ cells increases CDK4/6 kinase activity, driving cells into S-phase (Fig. 4G).In total, our results suggest that MTF2-contiaing PRC2.1 antagonizes G1 progression by repressing expression of the D-type cyclins CCND1 and CCND2.

Discussion
Regulated progression through cell cycle phases is critical to normal cellular function and viability, while disordered progression is the hallmark of many disease states.Although the cell cycle has been an area of active research for decades, our understanding of its regulation remains incomplete.Using a chemogenetic approach, we found that inactivation of members of PRC2.1, but not factors specific to PRC2.2, resulted in profound resistance to the CDK4/6 inhibitor palbociclib.Loss of PRC2.1 complex members led to upregulation of the D-type cyclins CCND1 and CCND2, resulting in increased RB1 phosphorylation and S-phase entry in palbociclib-treated cells.We propose that PRC2.1, but not PRC2.2,meditates H3K27me3 deposition in the promoters CCND1 and CCND2 through the recognition of the CpG islands.These results tie PRC2.1 directly to the regulation of G1 progression.
In the chemogenomic screens reported here, we recovered genes in a diverse array of biological pathways that resulted in sensitivity or resistance to well-characterized cell cycle inhibitors.In addition, we observed that inactivation of genes involved in mitochondrial homeostasis resulted in resistance to palbociclib.Small molecule inhibitors of EZH2 or the electron transport chain co-administered with palbociclib resulted in enhanced cell growth (Fig. 3D and 4D), supporting the observed chemical-genetic interaction seen in our screen.However, genes identified in genetic screens should be interpreted with caution.Reproducible, and sometimes robust interactions can sometimes result from complicated changes in doubling time or alterations to the physiologic state of the cell 61 .It was recently demonstrated that genes encoding members of the electron transport chain are over-represented in DepMap co-dependency data, due to the remarkable stability of these protein complexes, which results in phenotypic lag that can vary in different backgrounds 61 .While mitochondrial complex assembly factors as enriched in Metascape analysis of our camptothecin screen as well as in palbociclib (Fig. 2C), the enrichment was greater than 1600-fold more significant in palbociclib (Fig. 3A).Moreover, a number of reports have found increased oxygen consumption and ROS production due to greater number and size of mitochondria in cells treated palbociclib 82,83,90 .This is consistent with a direct effect of CDK4/6 activity on mitochondrial function.Thus, in the case of both the PRC2 and the mitochondrial gene cluster, our data and that of others suggest that these results represent a direct link between these pathways and CDK4/6 biology.
Recently, PRC2 subcomplex accessory proteins have been implicated in an increasing number of processes that define cellular identity, including stem cell maintenance, differentiation and cancer 35,36,91,92 .Despite the importance of controlled cellular division to each one of these processes, few reports have interrogated the roles of the different subcomplexes outside of stem cell model systems or specifically on their role in cell cycle regulation.Here, we show that cells that lose either MTF2 or SUZ12 continue to proliferate despite palbociclib blockade (Fig. 4E and 4G).These mutants show no apparent change in the proportion of cells undergoing apoptosis and display a greater proportion of cells entering S-phase in the presence of palbociclib, compared to wild-type or JARID2∆ cells.This increase is consistent with our findings that in MTF2∆ cell lines treated with palbociclib, a higher percentage of RB1 remains phosphorylated, while a similar increase is not seen in JARID2∆ cells.We surmise that the upregulation of CCND1 and CCND2 expression observed in cells lacking MTF2 results in increased CDK4/6 kinase activity that is sufficient to overcome palbociclib-mediated inhibition.Critically, we did not observe any significant changes in expression of other classic regulators of the CDK4/6-RB1-E2F pathway in either our CUT&RUN or RNA-Seq datasets (Fig. S4F  and S4G).While we cannot exclude the possibility that MTF2 inactivation alters the expression of other factors that influence G1 progression, we propose CCND1 and CCND2 represent major targets of PRC2.1 repression restraining G1 progression in HAP1 cells.Though a recent report demonstrated that low MTF2 expression leads to increased chemotherapeutic resistance in leukemia 93 and down regulation of MTF2 was correlated with poorer clinical outcomes in breast cancer 94 , more work is needed to determine whether D type cyclins are the critical PRC2.1 targets in tumors themselves.
Work over the past decade has implicated accessory proteins as critical for proper genomic localization of the PRC2 enzymatic core.However, reports differ on in what chromatin and cellular contexts these subcomplexes act.Data from both mouse and human ES cells has suggested that PRC2.1 and PRC2.2 have overlapping genomic occupancy 65,67 , and that either subcomplex alone is capable of maintaining pluripotency 65,67,95 .However, recent reports have found differing dependencies on these subcomplexes for proper distribution of H3K27me3 in cellular models of differentiation 27,33,[96][97][98] .For example, a recent study in a model of induced differentiation suggested that MTF2 is involved in the maintenance of repression of PRC2 genes, whereas JARID2 is important for de novo deposition of H3K27me3 critical for gene silencing through genes "pre-marked" with H2AK119ub1 33 .Conversely, PRC2.1 was shown to be required for the majority of H3K27me3 deposition during induced cell-fate transitions in mESCs, whereas PRC2.2 was not 27 .This study generated a triple knockout of all three PCL proteins (PHF1, MTF2 and PHF19), resulting in complete ablation of all PRC2.1 activity and did not probe the contribution of each accessory protein individually.Furthermore, MTF2 transcript levels are down regulated upon differentiation, whereas PHF1 and PHF19 levels increase 99 .These data suggest that the subunit composition of PRC2.1 changes during this process.In our experiments, MTF2 is the only PCL subunit important for D-type cyclin repression.These data are consistent with cell typespecific contributions of this class of proteins.
Using mutants of genes encoding subunits specific to either PRC2.1 or PRC2.2, we investigated the role of each subcomplex in cell cycle progression in HAP1 cells.In contrast to what has been demonstrated for ES cell lines where the two subcomplexes work synergically at the majority of sites 65,67 , we show that MTF2 is required for the majority of H3K27me3 deposition at CpG islands genome-wide and JARID2 was only partially required for H3K27me3 at these loci (Fig. 5E).Importantly, the presence of MTF2 is more critical than JARID2 for the accumulation of H3K27me3 directly upstream of annotated transcription start sites in CpG islands-containing promoters in HAP1 cells (Fig. 5F).Concordant with the patterns in H3K27me3 in promoters, we found that MTF2 loss resulted in a greater number of upregulated transcripts than JARID2 loss (Fig. 5B and Fig. S3C).Finally, MTF2∆ cells displayed a stronger correlation between genes with decreased promoter H3K27me3 levels and increased transcription than did JARID2∆ lines (Fig. 6C and Fig. S4C).However, we cannot exclude the possibility that AEBP2 plays a larger role in the activity of PRC2.2 than does JARID2 in these cells, as we identified AEBP2 as significantly, albeit modestly, increasing sensitivity to palbociclib in pooled knockout cells (Fig. 4A and Fig. 4E).As H3K27me3 peak distribution was altered in the JARID2∆ cell lines (Fig. S3A), loss of JARID2 could alter H3K27me3 sites distal to promoters to change chromosome architecture or enhancer-promoter interactions.Alternatively, genes upregulated by loss of either MTF2 or JARID2 which did not have a significant alteration in promoter H3K27 methylation could be indirect effects.A recent report found that while PRC2.2 activity was not required for establishment of H3K27me3 during differentiation, but was instead required for recruitment of a PRC1 complex required for higher level chromatin interactions 27 .Future studies will be necessary to fully understand the coordination between these complexes.
The efficacy of CDK4/6 inhibitors in the treatment of HR+/HER2-breast cancer demonstrates the success of applying basic knowledge of cell cycle regulation to the generation of clinically-relevant drugs.However, despite this success in the treatment of breast cancer, the efficacy of CDK4/6 inhibition is variable, with 10-20% of tumors primarily resistant and an additional 40% becoming resistant to these drugs within the first two years 100,101 .Moreover, CDK4/6 inhibitors are currently being explored for other tumor types, and these are each likely to have novel resistance mechanisms 92,102,103 .Thus, understanding perturbations in molecular pathways that can result in resistance to CDK4/6 inhibition could lead to improved patient responses and outcomes.In this study, we found that mutation of the PRC2.1 accessory protein MTF2 results in the development of resistance to palbociclib-induced proliferation reduction.Previously, EZH2, SUZ12, EED, MTF2 and JARID2 have all been suggested to not only act as oncogenes [104][105][106][107][108][109][110][111] , but also to have tumor suppressor activities 30,91,93,94,[112][113][114] , depending on the type of cancer.These observations not only underscore the context-dependent ramifications of mutation of these PRC2 complex members, but also may help inform the context in which CDK4/6 inhibitors are most efficacious.Clinical trials using CDK4/6 inhibitors in combination with other therapeutics are underway and the mutational status and expression levels of PRC2 subunits might serve as predictors of efficacy.
Requests for further information and reagents should be directed to and will be fulfilled by the Lead Contact, David Toczyski (dpt4darwin@gmail.com).

Cell Lines
Cas9 expressing HAP1 cells were cultured in IMDM supplemented with 4mM glutamine (Gibco), 10% Tetracycline-free FBS (Sigma-Aldrich) and either 1x Antibiotic, Antimycotic (Invitrogen) or 1% penicillinstreptomycin (Sigma-Aldrich).HAP1 cells stably expressed Cas9 were employed for the whole-genome screen, while for subsequent experiments, a HAP1 line harboring a doxycycline-inducible Cas9 was utilized.HEK293T cells used for the production of virus were cultured in DMEM supplemented with 2mM glutamine (Gibco), 10% Tetracycline-free FBS (Sigma-Aldrich) in 1x Antibiotic, Antimycotic (Invitrogen).Cells were detached from tissue culture dishes using 0.25% Trypsin (Gibco) and maintained at 37˚C, 5% CO2.Our laboratory conducts regular mycoplasma testing of cultured cells with the MycoAlert Mycoplasma Detection kit (Lonza), and no mycoplasma contamination of any cell line was detected during this study.

Genome-wide chemical screening
The lentiviral TKOv3 sgRNA library (Addgene #90294) was used to perform pooled genome-wide CRISPR knockout screens.The library contains 70,948 guides, targeting 18,053 protein-coding genes (4 guides/gene).Ninety million HAP1 cells stably expressing Cas9 were seeded into 15cm dishes and infected with TKOv3 lentivirus at a multiplicity of infection of roughly 0.3, such that every sgRNA is represented in approximately 200-300 cells after selection (>200-fold coverage).After 24 hours of infection, cells with successful viral integration were selected in 25mL IMDM medium containing 1µg/mL puromycin (Sigma-Aldrich).Selection took place for 48 hours.Following selection, cells were harvested, pooled, and split into 3 replicates of 15 million cells each to maintain >200-fold coverage of the sgRNA library (day 0).At day 3, each replicate was split such that every drug screen had a at least 15 million cells per replicate to maintain >200-fold coverage.The drug concentrations (IC 30 -IC 50 determined as described below) used in the genome-wide chemical screens were as follows: Palbociclib -0.7µM, Colchicine -9.2nM, Camptothecin -1nM.An increase in potency was observed for most drugs when used in the pooled screens, thus screening concentrations were adjusted to preserve IC 30 -IC 50 throughout each passage.Cells were subject to treatment with drug in 0.1% DMSO, or 0.1% DMSO alone.Drug-containing media was refreshed every 3 days, along with the passaging of cells and the collection of cell pellets.To preserve >200-fold coverage, 20 million cells were pelleted with every passage, from day 0 to day 18.
Genomic DNA extraction and sequencing library preparation were performed as described previously 53 .Briefly, genomic DNA from cell pellets were extracted using the Wizard Genomics DNA Purification Kit (Promega) and quantified using the Qubit dsDNA Broad Range Assay kit (Invitrogen).Sequencing libraries were prepared as described previously 122 .Briefly, two PCR amplification steps were performed to first enrich for the gRNA-regions in the genome and second, attach Illumina sequencing indices to the amplified regions.Sequencing libraries were prepared from 50µg of genomic DNA (200-fold library coverage) using the NEBNext Ultra II Q5 Polymerase (NEB).Primers used included Illumina TruSeq adapters with i5 and i7 indices.Barcoded libraries were gel-purified using the PureLink Quick Gel Extract kit (ThermoFisher) and sequenced on an Illumina HiSeq2500.

Drug concentrations for chemical screening
Drug dosing experiments were performed to determine screening concentrations.HAP1 cells stably expressing Cas9 were seeded at a density of 2.5 million cells per 15-cm dish.Cells were treated with 0.1% DMSO, or drug in 0.1% DMSO, 2 hours after seeding.Viable adherent cells were counted 2 days post-treatment on a Coulter counter, and inhibitory concentrations were determined.The following are ranges of drug concentrations used in the dosing experiments: Palbociclib: 1.5-10µM, Colchicine: 1.5-150nM, Camptothecin: 1-5nM.

Orobas Pipeline for Scoring Chemical Genetic interactions
The Orobas pipeline (version 0.5.0 -Supplemental File 1) was used to score chemical genetic interactions from the genome-wide CRISPR/Cas9 screen data.The process is summarized here, and the complete R code is provided as a supplemental file.sgRNAs were normalized to sequencing depth for each sample and the log 2 fold change (LFC) in sgRNA abundance was calculated for each condition relative to the corresponding T0 sample.Guides with fewer than 30 read counts in the T0 sample were filtered out from further analysis, and genes with fewer than 3 remaining guides post-filtering were also filtered out from scoring.Residual effects were computed for each gene by calculating the residual LFC between sgRNAs in treated vs. DMSO samples after averaging technical replicate LFCs.Residual effects were then M-A transformed and loess-normalized to account for potential skew and non-linearity present in the data, and per-gene effect sizes and FDRs were computed by applying the moderated t-test to normalized residual effects.Hits were called as genes with FDRs less than 0.4 and per-gene effect sizes greater than 0.5 or less than -0.5 (a complete list of effect sizes and FDRs is included as Table S1).
STRING Interaction Network Generation STRING networks were set to only display physical interactions scores that were returned with high confidence (0.7) and taken from text-mining, experiments and databases.
sgRNA Lentiviral Vector Cloning Oligos for sgRNA targets were designed to contain the 5' overhang CACCG-for the sense oligo 5' and for that antisense 3' over hang AAAC-and -C, respectively.10µM each of sense and antisense oligos (Integrated DNA Technologies) were mixed in 1x T4 DNA Ligase buffer and water to a total volume of 10µL.This mixture was heated to 95ºC for 5 mins, then oligos were annealed by decreasing the temperature at a rate of -0.1ºC/second till the mix reached 25ºC.Annealing reactions were diluted 1:10 with water and then 1µL was used to ligate into 100ng of BsmBI digested pLentiGuidePuro vector (Addgene #52963) in 1x T4 DNA Ligase Buffer.600 units of T4 DNA Ligase (NEB) and water to a total volume of 25µL.After incubating for 1hr at 37˚C, 2µL of the ligation reaction was transformed into BME pre-treated XL10-Gold cells (Agilent) per the manufacturer's instructions and plated on LB + 100µg/mL carbenicillin plates for selection.Plasmids recovered from single colonies were confirmed by Sanger sequencing.
Polyclonal and Monoclonal Knockout Generation Cas9 expressing HAP1 cells were transduced with pLentiGuidePuro vectors (Addgene #52963) expressing a single sgRNA (see Table S5 for sgRNA sequence).Lentiviral transduction was conducted at low MOI (~30%) following standard protocols.Integration of the sgRNA was selected with 1µg/mL puromycin for up to two days, followed by combined puromycin selection and Cas9 induction for three days with 1µg/mL doxycycline.This polyclonal pool of pLentiGuidePuro transduced cells was then used for "pooled" knockout experiments or used to generate monoclonal cell lines.Trypsinized, single cells were then sorted into individual wells in a 96-well plate using the Sony SH800 sorter (UCSF, LCA).Isolated single cell-derived colonies were screened for mutation by PCR, followed by Sanger sequencing of the purified PCR product and ICE analysis (Synthego) of the resulting chromatographs.Candidate clonal knockouts were then confirmed by Western blot.Only monoclonal lines that clearly displayed knockout alleles and had no protein product by Western blot were utilized further.
Competitive growth assays (GFP/BFP pooled knockouts and GFP+/monoclonal knockout pools) For pooled knockout competitive growth assays, HAP1 cells harboring an inducible Cas9 and expressing GFP and HAP1 cells expressing BFP (Hundley et al.) were mixed at a ratio of 1:4 GFP:BFP HAP1 cells into a single well, with three GFP/BFP cell mixtures for each gene targeted for inactivation.Mixtures were transduced at a low MOI with a pLentiGuidePuro vector expressing one sgRNA (three biological replicates per gene, sgRNAs in Table S5).After 24 hours of lentiviral transduction, pools of cells were selected with 1µg/mL puromycin for one day, followed by 1µg/mL puromycin and doxycycline for three days to select for sgRNA integration and to induce Cas9 expression.After three days of Cas9 induction, pools were split into media with or without palbociclib every three days, for eighteen days.The GFP/BFP ratio was monitored on the Attune NxT (Invitrogen) flow cytometer every three days.FlowJo v10 was used to determine the GFP/BFP ratio at each time point.The ratio of GFP to BFP was normalized to the day zero ratio (prior to splitting into palbociclib), and subsequently to the matched untreated ratio at each time point.
For HAP1 GFP+/-competitive growth assays, clonal MTF2∆ or JARID2∆ GFP-cells were mixed with HAP1 GFP + clones at 1:4 GFP-:GFP+ ratio, split into media with or without drug, and analyzed by flow cytometry as described above.
Western blotting Harvested cell pellets were lysed in 1xRIPA buffer supplemented with 1x EDTA-free cOmplete protease inhibitor (Roche) and 1x PhosphoSTOP phosphatase inhibitor (Roche) for 30 minutes on ice with two rounds of 15 second vortexing.Lysates were cleared at 21,000 xg for 10mins at 4˚C.Protein concentration was determined by BCA assay and BSA standard curve (Pierce), and samples were adjusted to 1µg/µL total protein with 1xRIPA and SDS-PAGE sample loading buffer was added (62.5 mM Tris-HCl (pH 6.8), 2.5% SDS, 0.1% bromophenol blue, 10% glycerol, 5% β-mercaptoethanol (v/v)).10µg of total protein was loaded per lane onto a 4-20% Criterion Tris-HCl Protein gel (Bio-Rad) and separated by electrophoresis at 150 V for 1hr.Proteins were transferred and immobilized onto a nitrocellulose membrane (GE Healthcare) by electrophoresis for 1h at 100V in standard transfer buffer containing 20% methanol.Membranes were blocked for an hour at room temperature and then then probed overnight in a 1:1000 dilution of 1˚ antibody (unless otherwise indicated) at 4˚C and in a 1:10,000 2˚ antibody at room temperature for 1hr at in the appropriate blocking buffer.Chemiluminescent and fluorescent signals were visualized with an Odyssey FC imager (LICOR).
Cell Cycle Analysis by Propidium Iodine 200,000 cells/well were plated in 6 well dishes at, as to be 10-20% confluent at the time of treatment.Cells were treated with inhibitors 24h after plating, then harvested 48h later by trypsinization, washed twice with cold 1x PBS, fixed by dropwise addition of ice-cold 70% ethanol, and incubated at 4˚C overnight.Fixed samples were washed twice with 1x PBS + 1% BSA prior to resuspension in a solution of 1x PBS, 1mg/mL RNase A and 50 µg/mL propidium iodide for 1 hour at 37°C.DNA content of at least 20,000 single cells was determined by Attune NxT flow cytometer (Invitrogen), and data were analyzed using FlowJo v10.
BrdU Incorporation Assay 250,000 cells/well were plated in 6-well dishes and grown for 24hrs prior to treatment.Cells were then treated with either DMSO (mock) or 1.5µM palbociclib for a total of 24hrs, with 10µM BrdU being added to the culture medium 1hr prior to harvesting.Cells were counted using the Countess automatic hemocytometer (Invitrogen) to ensure that only 1 million cells were stained.Cells were prepared for analysis using BD Pharmagen BrdU Flow Kits (BD Biosciences) according to the manufacturer's instructions.BrdU incorporation was determined for at least 20,000 single cells by Attune NxT flow cytometer (Invitrogen), and data were analyzed using FlowJo v10.
Quantitative Crystal Violet Proliferation Assay 1mL of a 1,000 cells/mL suspension were seeded into a per well in a 6-well plate containing 1mL IMDM supplemented with double the indicated concentration of palbociclib and GSK126 in technical triplicate.Cells were allowed to proliferate for nine days, with the media supplemented with the drug at the concentration indicated replaced every three days.After nine days, cells were washed once with 1x PBS, followed by staining and fixation in a 0.25% Crystal Violet, 20% methanol solution for 10mins at room temperature.Following staining, cells were washed six times with 1x PBS and lysed in a 100mM sodium citrate and 50% ethanol solution for 30mins at room temperature on an orbital shaker.Lysates were recovered and absorbance at 590nM was detected using a Synergy Neo2 Microplate Reader (BioTek).Proliferation at each concentration was determined relative to untreated wells.
PrestoBlue Proliferation Assay 45µL of a 50,000 cells/mL cell suspension was seeded into a 96-well plate containing 45µL of IMDM supplemented with the indicated concentration of palbociclib, antimycin A, TTFA or oligomycin in triplicate.After proliferation for 48hrs, 10µL of PrestoBlue (Invitrogen) was added to each well and incubated for 30mins at 37˚C.Conversion of PrestoBlue was determined by recording the fluorescence excitation at 560nM and emission at 590nM using a Synergy Neo2 Microplate Reader (BioTek).Proliferation at each concentration was determined relative to untreated wells.RNA-Extraction 150,000 cells were seeded into 6-well plates and allowed to grow overnight.Cells were treated with DMSO (Mock) or 1.5µM palbociclib for 24hrs prior to harvesting directly in TRIzol reagent (Invitrogen).After chloroform extraction, the aqueous phase was transferred to a fresh tube and 1 volume of 100% ethanol was added and mixed thoroughly before binding to an RNA Clean & Concentrator (Zymo).RNA was DNase I digested on-column (Zymo), purified according the manufacturer's instructions and eluted in ]nuclease-free water.To prepare RNA-Seq libraries, 2µg of total RNA was polyA, followed by Illimina adaptor ligation and paired-end sequencing on an Illumina HiSeq at a depth of at least 22 million reads per sample by Azenta.
First Strand cDNA Synthesis and qRT-PCR 2µg of total RNA was first heat denatured in the presence of dNTPs and oligo-dT at 65˚C for 5mins.RNase inhibitor and Tetro reverse transcriptase (Bioline) was then added to heat denatured total RNA and cDNA was synthesized at 45˚C for 1hr, followed by heat inactivation at 85˚C for 5mins.cDNA synthesis reactions were then diluted 1:5 and 2µL was added into qRT-PCR reaction mix, utilizing SensiFast Lo-ROX qRT-PCR Mastermix (Bioline) in both biological and technical triplicate.Reactions were carried out and analyzed using a QuantStudio5 machine (Applied Biosystems).See Table S5 for qRT-PCR primer sequences.CUT&RUN Library Preparation CUT&RUN libraries were generated by first lysing 300,000-500,000 cells in 500µL of Nuclei Extraction Buffer (20mM HEPES-KOH pH 7.9, 10mM KCl, 1mM MgCl 2 , 0.1% Triton X-100, 20% glycerol, and 1x protease inhibitor) for 10 minutes on ice.. Next, samples were spun down and washed twice with Nuclei Extraction Buffer before being resuspended in 500µL nuclei extraction buffer.10µL of Concanavalin A-coated beads (EpiCypher) previously washed in Wash Buffer (20mM HEPES-KOH pH 7.5, 150mM NaCl, 2mM EDTA, 0.5mM spermidine, and 1x protease inhibitor) and resuspended in Binding Buffer (20mM HEPES-KOH pH 7.5, 1mM CaCl 2 , and 1mM MnCl 2 ) were then added to the samples and incubated with rotation for 15 min at 4°C.Next, samples were washed once with Binding Buffer before being resuspended in 50 µL of Buffer 2 containing 0.1%BSA, 2µM EDTA and 0.5µL H3K27me3 1˚ antibody, followed by overnight incubation with rotation at 4°C.Following the incubation, samples were washed twice with Buffer 2 before being incubated in 50µL of Buffer 2 containing ~700ng/mL protein A-MNase fusion protein (Batch #6 from the Henikoff Lab) for 1 hour with rotation at 4°C.Samples were washed two more times and resuspended in 100 µL of Buffer 2 before starting the MNase digestion by adding CaCl 2 to a concentration of 2mM on ice for 30 minutes, after which the reaction was quenched with the addition of 100 µL 2X Stop Buffer (200mM NaCl, 20mM EDTA, 4mM EGTA, 50µg/mL RNase A, 40µg/mL GlycoBlue (Ambion), and 2pg/mL spike-in DNA) to inactivate the MNase.Samples were incubated for 30 min at 37°C and spun down for 5 minutes at 4°C to release DNA fragments.DNA was phenol:chloroform extracted and 200µL of the recovered aqueous phase was ethanol precipitated with 500µL ethanol, 20µL 3M NaOAc, 2µL GlycoBlue at -80˚C.Libraries were prepared using 2S Plus DNA Library Kit adapters (Swift Biosciences) and size-selected using SPRIselect beads (Beckman Coulter) before being amplified and sent for paired-end sequencing on the NovaSeq 6000 (150 bp reads).
CUT&RUN Processing and Analysis CUT&RUN paired reads were aligned to a reference human genome (hg38) by the bwa-mem algorithm.PCR duplicate reads were removed by Picardand peaks were called using macs2 with the broad flag and an FDR of 0.05.Bedtools intersect was used to identify reproducible peaks between biological replicates of each condition, and reproducible peaks from each condition were compiled into a list.Bedtools multicov was used to build a matrix with the number of reads from each dataset falling in each region in this list.This matrix was used for all 'genome wide' analyses.Bedtools multicov was also used to build a matrix with the number of reads from each dataset in a 5,000bp window around the transcription start site (4kb upstream, 1kb downstream) of all hg38 genes defined by gencode v41.The gencode v41annotation for CCND2 was originally incorrectly assigned to chr12:4,265,771-4,270,771 and reassigned using the Refseq coordinates chr12:4,269,762-4,274,762. Count matrices were analyzed with DESeq2 to compare changes in H3K27me3 deposition globally, and changes in H3K27me3 deposition in promoters.For heatmaps, deduplicated BAM files were converted to bigwigs and BED files and normalized reads per kilobase per million mapped read using deepTools bamCoverage.For genome wide analyses, H3K27me3 CUT&RUN signal in normalized bigwigs was measured using deepTools computeMatrix in 10kb regions centered around WT peaks overlapping with CpG islands.For promoter analyses, H3K27me3 CUT&RUN signal in normalized bigwigs was measured using deepTools computeMatrix in 5kb regions (4kb upstream, 1kb downstream) around transcription start sites for promoters overlapping with CpG islands.Promoters with the highest average H3K27me3 signal intensity in Wild-type replicate 1 and sorted in descending order.RNA-Seq Processing and Analysis RNA-seq paired reads were quantified using Salmon.Transcript-level abundance estimates from Salmon and gene-level count matrices were created using Tximport and analyzed using DESeq2.Paired reads were aligned using STAR to generate BAM files.BAM files were converted to BED files using bamCoverage and normalized using RPKMs and to effective genome size of hg38 (2,913,022,398) with a bin size of 10. (B): Volcano plots and 1nM camptothecin whole-genome screen results.The Orobas calculated "Differential Gene Effect" was plotted against the -log 10 (p-value) for this effect for each gene targeted in the screen, as calculated by the Orobas pipeline.Red dotted line indicates the established cut-off.Highlighted dots are genes with known roles in response to each treatment, with blue or yellow dots indicate genes that when inactivated resulted in sensitivity or resistance, respectively, to camptothecin.(C): Representative STRING analysis networks for protein complexes with known roles in pathways that we identified as sensitive in our 1nM camptothecin chemogenomics screen.Blue dots in the STRING network indicate genes that when inactivated resulted in sensitivity to camptothecin.(D): Same as in (B) but for 9.2nM colchicine whole-genome screen results.(E): Same as in (C) but for 9.2nM colchicine whole-genome screen results.(F): Same as in (B) but for 0.7µM palbociclib whole-genome screen results.(G): Same as in (C) but for 0.7µM palbociclib whole-genome screen results.(B): Volcano plot of novel genes identified in the 1nM camptothecin chemogenomics screen.Volcano plots are plotted as in (Fig. 1B), with highlighted dots representing novel genes identified in the camptothecin screen.(C): Dot plot of Metascape analysis of significant genes that sensitized or de-sensitized cells to 1nM camptothecin.The -log 10 (p-value) was plotted against each term.The Metascape determined enrichment for observed over expected genes associated with a given term is indicated by color of circle.The percentage of the genes above our established cut-off associated with a given term is indicated by the size of the circle.(D): STRING analysis of genes identified from our analysis of the 1nM camptothecin chemogenomics screen.(E), (F) and (G) Same as in (B), (C) and (D) except for 9.2nM colchicine chemogenomics screen.

Figure 1 :
Figure 1: Chemogenomics CRISPR-Cas9 Screen to Study Cell Cycle Progression (A): Schematic of whole-genome CRISPR-Cas9 screen.(B):Volcano plots and 1nM camptothecin whole-genome screen results.The Orobas calculated "Differential Gene Effect" was plotted against the -log 10 (p-value) for this effect for each gene targeted in the screen, as calculated by the Orobas pipeline.Red dotted line indicates the established cut-off.Highlighted dots are genes with known roles in response to each treatment, with blue or yellow dots indicate genes that when inactivated resulted in sensitivity or resistance, respectively, to camptothecin.(C): Representative STRING analysis networks for protein complexes with known roles in pathways that we identified as sensitive in our 1nM camptothecin chemogenomics screen.Blue dots in the STRING network indicate genes that when inactivated resulted in sensitivity to camptothecin.(D): Same as in (B) but for 9.2nM colchicine whole-genome screen results.(E): Same as in (C) but for 9.2nM colchicine whole-genome screen results.(F): Same as in (B) but for 0.7µM palbociclib whole-genome screen results.(G): Same as in (C) but for 0.7µM palbociclib whole-genome screen results.

Figure 2 :
Figure 2: Analysis of Camptothecin and Colchicine Chemogenomics Screen Reveals Novel Players in Cell Cycle Regulation (A): Dot plot comparison of the effect of gene mutation across three different screen conditions.Circle color indicates the strength of the positive or negative differential gene effect, circle size indicates its -log 10 (p-value) from our Orobas analysis.(B):Volcano plot of novel genes identified in the 1nM camptothecin chemogenomics screen.Volcano plots are plotted as in (Fig.1B), with highlighted dots representing novel genes identified in the camptothecin screen.(C): Dot plot of Metascape analysis of significant genes that sensitized or de-sensitized cells to 1nM camptothecin.The -log 10 (p-value) was plotted against each term.The Metascape determined enrichment for observed over expected genes associated with a given term is indicated by color of circle.The percentage of the genes above our established cut-off associated with a given term is indicated by the size of the circle.(D): STRING analysis of genes identified from our analysis of the 1nM camptothecin chemogenomics screen.(E), (F) and (G) Same as in (B), (C) and (D) except for 9.2nM colchicine chemogenomics screen.

Figure 3 :
Figure 3: Mutation of Mitochondria Genes Attenuates the Sensitivity to Palbociclib(A): Dot plot of Metascape analysis of significant genes in the 0.7µM palbociclib chemogenomics screen.The -log 10 (p-value) was plotted against each term.The Metascape determined enrichment for observed over expected genes associated with a given term is indicated by color of circle.The percentage of the genes above our established cut-off are associated with a given term is indicated by the size of the circle.(B): Volcano plot of genes identified from our analysis of the 0.7µM palbociclib chemogenomics screen.Volcano plots are plotted as in Fig.1D, with highlighted dots representing novel genes identified in the palbociclib screen.(C): STRING networks of novel protein complexes identified in palbociclib screen.Dots in the STRING network indicate genes that when inactivated resulted in sensitivity (blue, left) of resistance (yellow, right) to 0.7µM palbociclib.(D): Dose-response curve of palbociclib-induced growth inhibition with oxidative phosphorylation inhibitors.HAP1 cells were exposed to palbociclib with or without increasing concentrations of the oxidative phosphorylation inhibitors rotenone, TTFA and oligomycin.Cell proliferation in each treatment was determined after 48 hours by PrestoBlue assay.All data were normalized to the initial dose of each oxidative phosphorylation inhibitor in indicated concentration of palbociclib.Error bars represent the mean of three replicates ±StdDev.

Figure 4 :
Figure 4: Loss of Polycomb Repressive Complex Components Display Specific Resistance to Palbociclib (A): Volcano plots in Fig. 3B except with members of PR-DUB, Polycomb Repressive Complex 1 and 2 highlighted.(B): STRING analysis network of Polycomb Repressive Complex components.Yellow dots in STRING network indicate that inactivation of these genes conferred resistance to 0.7µM palbociclib in our chemogenomics screen.(C): Dot plot of comparison of the effect of PRC2 complex member gene mutation across three different screen conditions, as in Fig. 2B.(D): Dose-response curve of a nine-day quantitative Crystal Violet assay demonstrating rescue of palbociclibinduced growth inhibition with GSK126.All data were normalized to untreated cells and represents the mean of three technical replicates ±StdDev.(E):Results of competitive growth assay for each indicated time point, normalized to the initial GFP+/GFPratio of the pool.The performance of each sgRNA in 1.5µM palbociclib vs Mock is shown, after normalizing to the control sgRNA.Data represents the mean ±SEM of the GFP+/GFP-ratios of three independent sgRNAs.(F): Dose-response curve of nine-day quantitative Crystal Violet assay results demonstrating rescue of palbociclib-induced growth in MTF2∆ and JARID2∆ cells.Data represent mean Crystal Violet staining of cell culture wells of three independently isolated monoclonal knockout cell lines ±StdDev.