SUV39 SET domains mediate crosstalk of heterochromatic histone marks
Abstract
The SUV39 class of methyltransferase enzymes deposits histone H3 lysine 9 di- and trimethylation (H3K9me2/3), the hallmark of constitutive heterochromatin. How these enzymes are regulated to mark specific genomic regions as heterochromatic is poorly understood. Clr4 is the sole H3K9me2/3 methyltransferase in the fission yeast Schizosaccharomyces pombe, and recent evidence suggests that ubiquitination of lysine 14 on histone H3 (H3K14ub) plays a key role in H3K9 methylation. However, the molecular mechanism of this regulation and its role in heterochromatin formation remain to be determined. Our structure-function approach shows that the H3K14ub substrate binds specifically and tightly to the catalytic domain of Clr4, and thereby stimulates the enzyme by over 250-fold. Mutations that disrupt this mechanism lead to a loss of H3K9me2/3 and abolish heterochromatin silencing similar to clr4 deletion. Comparison with mammalian SET domain proteins suggests that the Clr4 SET domain harbors a conserved sensor for H3K14ub, which mediates licensing of heterochromatin formation.
Introduction
Histone methylation plays a critical role in regulating and organizing genomes to maintain genome integrity and establish appropriate gene regulation programs (Greer and Shi, 2012; Janssen et al., 2018). Histone H3 lysine 9 di- and trimethylation (H3K9me2/3) are highly conserved hallmarks of heterochromatin, where they provide a platform for the recruitment of chromatin effector complexes. H3K9me2/3 marks are dysregulated in various cancers, in neurological diseases and viral latency, and are targeted for therapeutic approaches (Kaniskan et al., 2018).
The SUV39 clade of the SET-domain family of protein lysine methyltransferases represents the enzyme class depositing H3K9me2/3 (Dillon et al., 2005). They share a SET domain as a catalytic core, which builds on a series of curved β-sheets that are sandwiched between distinct pre- and post-SET domains, providing structural elements based on zinc-finger motifs. The SUV39 enzymatic activity facilitates the transfer of a methyl group from S-adenosyl-L-methionine (SAM) to the ε-amino group of lysine 9 on histone H3 with high specificity, thereby converting the cofactor to S-adenosyl-L-homocysteine (SAH).
The mechanisms regulating SUV39 proteins to specifically deposit H3K9me2/3 in heterochromatic regions are not well understood. The fission yeast Schizosaccharomyces pombe heterochromatin system, closely related to the metazoan systems, with its sole SUV39 protein Clr4, serves as a paradigm for studying the role of this enzyme class in heterochromatin establishment and maintenance. Ectopic tethering experiments with Clr4 have established that H3K9me2/3 can support stable epigenetic transmission of transcriptional states (Audergon et al., 2015; Ragunathan et al., 2015). Clr4 has further served as a valuable tool to understand its animal orthologs, which share a very similar domain architecture consisting of an N-terminal chromodomain and a C-terminal catalytic SET domain (Figure 1A). The chromodomain binds the H3K9me2/3 mark and is critical for spreading the H3K9me2/3 mark across genomic domains (Melcher et al., 2000; Müller et al., 2016; Zhang et al., 2008). Recent evidence also suggests that an autoinhibitory loop in the post-SET domain regulates spreading of H3K9me2/3 by Clr4 (Iglesias et al., 2018).
Clr4 is a key part of S. pombe heterochromatin, which is found at pericentromeres, telomeres, and the mating type locus and depends on the tightly interconnected pathways of the nuclear RNA interference machinery and an array of chromatin modifiers (Grewal, 2010). Small RNAs are generated from a limited set of repeat sequences in heterochromatic regions and guide the RNAi machinery in the form of the RNA-induced initiation of transcriptional gene silencing complex (RITS) to nascent heterochromatic transcripts (Verdel et al., 2004; Bühler et al., 2007; Schalch et al., 2011; Shimada et al., 2016). RITS recruits Clr4 and its associated Cullin4-RING ubiquitin ligase complex, the CLRC complex (Zhang et al., 2008; Bayne et al., 2010), thereby driving the deposition of H3K9me2/3 marks, which provides a platform for binding of HP1 proteins that in turn recruit further transcriptional gene silencing complexes (Motamedi et al., 2008; Fischer et al., 2009; Leopold et al., 2019). Heterochromatin formation is further strongly dependent on deacetylation of H3K14 (Alper et al., 2013; Bjerling et al., 2002; Buscaino et al., 2013; Yamada et al., 2005), with H3K14 acetylation being strongly associated with active gene expression in euchromatin (Wang et al., 2012).
Even though it has long been recognized that rik1 plays an essential role in the S. pombe heterochromatin system (Ekwall and Ruusala, 1994; Allshire et al., 1995; Nakayama et al., 2001; Zhang et al., 2008), its molecular mechanism has remained elusive. Rik1 is a paralog of the DNA repair protein DDB1 and forms the central scaffold of the CLRC complex that includes the cullin protein Cul4, the WD-40 β-propeller proteins Dos1/Raf1, the RING finger protein Pip1, as well as the replication focus targeting sequence (RFTS) protein Dos2/Raf2 and Clr4 (Hong et al., 2005; Horn et al., 2005; Jia et al., 2005; Li et al., 2005; Thon et al., 2005). CLRC features the hallmarks of an intact CRL4-type E3 ligase, which uses the subunit Dos1/Raf1 as a substrate adapter (Buscaino et al., 2012; Horn et al., 2005; Jia et al., 2005; Kuscu et al., 2014). Recent work has revealed that the preferred substrate for ubiquitylation by CLRC is lysine 14 on histone H3, yielding H3K14ub (Oya et al., 2019), and that H3K14ub controls the activity of Clr4. However, how H3K14ub regulates Clr4 and its role in heterochromatin formation remain to be determined.
We combined biochemical and structural methods to decipher the molecular mechanism governing the stimulation of Clr4 by H3K14ub. These experiments identified a ubiquitin binding region (UBR) in the catalytic domain of Clr4, which mediates H3K14ub binding, and mutations that disrupt this interaction interface lose heterochromatin function. This reveals a critical regulatory mechanism that uses the SET domain of Clr4 for control of heterochromatin formation by an epigenetic crosstalk.
Results
The catalytic domain of Clr4 senses the presence of H3K14-linked ubiquitin
To understand the enzymatic reaction underlying the stimulation of Clr4 by H3K14ub, we set out to determine the domains of Clr4 that mediate this effect. Clr4’s domain architecture comprises an N-terminal chromodomain, which is connected by a linker region to the catalytic lysine methyltransferase (KMT) domain (Figure 1A). The KMT domain contains the N-terminal, pre-SET, SET, and post-SET regions and its structure has been determined (Min et al., 2002; Iglesias et al., 2018). To determine the role of the chromodomain and linker region, we recombinantly expressed both the full-length Clr4 and the isolated KMT domain and quantified their methyltransferase activity using two different assays: (1) incorporation of tritium-labeled methyl groups from a 3H-SAM donor into histone proteins and (2) fluorescence-based measurement of SAH generation (Cisbio EPIgeneous Methyltransferase Assay kit). As substrates we used unmodified H3 peptide (1–21) and H3 peptide ubiquitinated on K14 (H3K14ub) or K18 (H3K18ub) generated by native chemical ligation procedures. These branched peptides are identical to a native ubiquitin linkage except for a glycine to alanine change in the C-terminal residue of ubiquitin that is covalently linked to K14 (Figure 1—figure supplement 1A).
Figure 1B shows fluorographs of methyltransferase assays for Clr4 full-length and the KMT domain, which both manifest a similar degree of strong stimulation by H3K14ub (Figure 1C, Figure 1—figure supplement 1E–H). To get a sense for the steric requirements involved in stimulation by H3K14ub, we tested H3K18ub, the closest downstream lysine on the H3 tail. H3K18ub is a significantly poorer substrate than H3K14ub, demonstrating that the position of the ubiquitin modification on the H3 tail is critical for optimal Clr4 stimulation and that covalent linkage generating high local ubiquitin concentration is not sufficient. Quantification of the stimulation by the different substrates further showed that H3K14ub elicits a 40–50-fold boost in activity at a substrate concentration of 20 μM (Figure 1C). These experiments confirm the observation by Oya et al., 2019 that the H3K14ub substrate triggers a dramatic and specific increase in the methyltransferase activity of Clr4. However, in contrast to the previous study, we observe that the KMT domain is sufficient to mediate this regulatory mechanism. Since the primary effect of H3K14ub traces to the catalytic domain, we decided to further investigate the stimulation mechanism using the isolated KMT domain.
To fully characterize the H3K14ub-mediated stimulation, we measured enzyme kinetics of Clr4KMT on H3K14ub vs. the unmodified peptide (Figure 1D). An approximately hundred-fold difference for the Michaelis–Menten (KM) constants was measured with 0.33 ± 0.06 µM for H3K14ub and 28.1 ± 8.0 µM for H3. The turnover number (kcat) was about three times higher in the presence of H3K14ub (0.81 ± 0.03 min–1) compared with unmodified H3 (0.28 ± 0.02 min–1). This leads to an increase in overall enzymatic efficiency measured by the specificity constant ksp = kcat/KM of 250-fold (Table 1). Comparison of the kinetic parameters between H3K14ub and H3 substrate indicates that the presence of ubiquitin on lysine 14 leads to a tighter enzyme-substrate complex and to conformational changes in the active site that increase the rate of the methyltransferase reaction.
To determine whether H3K14ub uses an allosteric site for ubiquitin on Clr4, we challenged the methyltransferase reaction with increasing amounts of free ubiquitin. While we observed no significant increase in activity for unmodified H3, we observed a drop in the activity for H3K14ub at high concentrations of free ubiquitin (Figure 1E). This experiment failed to produce evidence of an allosteric site for free ubiquitin on Clr4, and we conclude that the stimulation of kcat is likely to depend on an induced-fit mechanism triggered by binding of H3K14ub to the Clr4KMT domain.
H3K14 ubiquitin mark affects the structural dynamics of Clr4
To understand the structural basis for the regulation of Clr4 by H3K14ub, we performed hydrogen/deuterium exchange coupled to mass spectrometry (HDX-MS) analysis on the free Clr4KMT domain and on Clr4KMT in complex with H3- and H3K14ub peptides. HDX-MS measures protein dynamics based on the rate of exchange of protein amide protons with the solvent (Kochert et al., 2018). Changes in HDX rates upon complex formation identify regions of the protein that are affected by the formation of the complex. Comparing the dynamics for Clr4KMT alone and in the presence of excess H3K14ub or unmodified H3 peptides, we observed a strong and unique reduction of the HDX rate for residues 243–261 of Clr4 in the presence of H3K14ub (Figure 2A and B, Figure 2—figure supplement 1, Supplementary files 2-4). A further region between residues 291–305 showed significant amide protection by both H3 and H3K14ub substrates, the latter peptide showing a more intense protection (Figure 2B, bottom graph). These results suggest that the H3 peptide interacts with residues 291–305, while the ubiquitin moiety binds to a region involving residues 243–261. The H3K14ub binding region identified by HDX-MS maps to the NT domain just before it transitions into the pre-SET domain, and we will refer to this region as UBR (Figure 2A). The UBR forms a ridge along the ‘back’ of Clr4, opposite to the active site pockets where cofactor and the substrate peptide bind (Figure 2C).
We also compared HDX-MS rates of isolated H3K14ub peptide with the H3K14ub-Clr4 complex and found that amino acids 25–43 on ubiquitin were protected from exchange upon interaction with Clr4 (Figure 2D and E, Figure 2—figure supplement 1, Supplementary file 3). These residues map to the α-helix of ubiquitin, indicating that this surface interacts with Clr4.
UBR mutants lose affinity for H3K14ub
To determine the functional importance of the UBR, we designed mutations that specifically disrupt the Clr4-H3K14ub interaction. While substitutions of residues 243–251 yielded unstable protein, substitution of residues 253–256 (sequence DPNF to GGSG, referred to as Clr4-GS253) resulted in a stable protein. Based on the Clr4 structure and sequence conservation, we further chose to mutate three phenylalanines that intersect orthogonally with the GS253 mutations on the surface of Clr4 (F256A, F310A, and F427A, referred to as Clr4-3FA).
We used isothermal titration calorimetry (ITC) to determine the effect of the mutations on the affinity of the Clr4-H3K14ub interaction. Consistent with the low KM observed previously, Clr4 binds to the H3K14ub peptide with a dissociation constant (Kd) of 80 ± 6 nM in a reaction that is dominated by enthalpy (Figure 3A). In contrast, Clr4-GS253 and Clr4-3FA showed complete loss of binding to H3K14ub under the same conditions (Figure 3B and C). We also attempted to determine the Kd for unmodified H3 peptides, but were unable to observe binding. The ITC data confirm that H3K14ub binds with high affinity to the KMT domain of Clr4 and that the UBR mutants disrupt the Clr4-H3K14ub enzyme-substrate complex.
The Clr4-3FA mutant showed very stable biochemical behavior and crystallized readily. To determine the effect of the 3FA mutations on the folding of Clr4, we solved the structure by X-ray crystallography. These efforts resulted in a 2.46 Å structure of Clr4-3FA (Table 2). The packing of the molecules in these crystals is similar to the packing observed in the structure of the autoinhibited Clr4, despite a difference in space group (PDBID:6BOX) (Iglesias et al., 2018). When comparing with previous X-ray structures of Clr4, we observed no significant difference in global or local protein folding between the wild-type Clr4 KMT domains and the Clr4-3FA mutant (RMSD = 0.37 Å) (Figure 3D and E, Figure 3—figure supplement 1). This indicates that the Clr4-3FA mutations do not induce significant structural changes.
To establish the impact of the mutants on the methyltransferase activity, we performed enzymatic assays using 3H-SAM and gel-based read-out. These assays show that the mutants lose activity on H3K14ub when compared to wild-type Clr4KMT (Figure 3F). To establish that the UBR mutants specifically target H3K14ub-mediated stimulation of Clr4 but not activity on unmodified peptide, we determined Michaelis–Menten kinetics for the Clr4-GS253 and Clr4-3FA mutants using a luminescence-based assay (Figure 3G, Table 1). Comparison of the specificity constants ksp shows that the efficiency of the UBR mutants on H3K14ub drops approximately 27-fold for Clr4-3FA and 48-fold for Clr4-GS253 compared to the wild-type. This drop in activity is caused by a corresponding drop in KM, consistent with the ITC results. The difference between Clr4-3FA and Clr4-GS253 is not statistically significant, but is consistent with small differences observed in the radioactive gel-based and filter-binding assays (Figure 3F, Figure 3—figure supplement 1C).
In contrast to the enzymatic activity on H3K14ub, the activity on unmodified H3 peptides is only mildly affected resulting in a ksp, that is, 91% of wild-type for Clr4-3FA and 78% for Clr4-GS253 (Figure 3G, Table 1). We conclude that the enzymatic and affinity measurements establish that both UBR mutants target the Clr4-H3K14ub interaction with high specificity.
Clr4 loss-of-function mutants affect heterochromatin
In S. pombe, heterochromatin formation depends strongly on the H3K9me2/3 methyltransferase activity of Clr4. Clr4 knock-out and mutant strains are defective for transcriptional gene silencing and display hyperacetylated heterochromatic regions (Nakayama et al., 2001; Bannister et al., 2001; Gerace et al., 2010; Iglesias et al., 2018). The UBR mutants, therefore, provide a unique opportunity to investigate the functional importance of the H3K14ub-mediated stimulation for the S. pombe heterochromatin system. We introduced the Clr4-GS253 and Clr4-3FA mutations at the endogenous clr4 locus and crossed them into a dual reporter (imr1L::ura4+/otr1R::ade6+) background for evaluating heterochromatin silencing using comparative growth assays (Ekwall et al., 1997). Strains with a functional silencing machinery are able to repress transcripts from these two centromeric reporters, and wild-type cells therefore show restricted growth on medium lacking uracil and turn red on medium containing low concentrations of adenine. Elevated levels of ura4, however, render the strains sensitive to growth on media containing 5-fluoroorotic acid (FOA), which is converted into the toxic product fluorodeoxyuridine when ura4 is expressed at elevated levels. The clr4-GS253 mutant showed increased growth on medium lacking uracil, no growth on FOA plates, and white colonies on low adenine, comparable to the clr4Δ mutant, which indicates a severe loss of transcriptional gene silencing (Figure 4A). The clr4-3FA mutant also showed a growth phenotype very similar to clr4Δ. However, the slightly pinkish color on low Ade plates suggests that the silencing defect is less severe than in clr4-GS253 or clr4Δ. The clr4 mutants are expressed at normal levels, and co-IP experiments with Rik1 show that they remain associated with the CLRC complex (Figure 4B and C). To get a quantitative measure of the loss of transcriptional gene silencing, we analyzed endogenous heterochromatic transcripts at centromeric dg/dh repeats and the subtelomeric tlh1 locus by RT-qPCR (Figure 4D). The levels of these transcripts were greatly elevated in the clr4-GS253 mutant, similar to clr4Δ. In contrast, the 3FA mutant showed a more nuanced silencing defect with high transcript levels for centromeric dh and telomeric transcripts, and a modest ~10-fold increase for centromeric dg. These findings are consistent with the weaker silencing defect of clr4-3FA observed in the growth assays and with the weaker loss of enzymatic stimulation in vitro. To test if the loss of gene silencing is associated with loss of H3K9 methylation, we performed chromatin immunoprecipitation (ChIP) against H3K9me1, H3K9me2, and H3K9me3 (Figure 4E and F). Both H3K9me2 and H3K9me3 were completely abolished in clr4-GS253 and clr4-3FA strains at centromeric dg/dh repeats, being indistinguishable from clr4Δ. Currently, we cannot explain why the dg transcripts in the 3FA mutant are only slightly elevated while completely losing H3K9me2/3. clr4-GS253 and clr4Δ also showed elevated levels of H3K9me1 when compared to wild-type clr4+, while this was not observed for clr4-3FA. The elevated H3K9me1 marks support the hypothesis that another unidentified methyltransferase could be depositing this mark (Jih et al., 2017). We further found that in agreement with the increased transcript levels RNA polymerase II occupancy at dh repeats was strongly increased.
In summary, we observe defects in transcriptional gene silencing and heterochromatin formation for the UBR mutants that are very similar to clr4Δ. Furthermore, the severity of the phenotype correlates with the degree of loss of enzymatic function on H3K14ub in vitro, consistent with Clr4’s KMT domain mediating the crosstalk between H3K14ub and H3K9me2/3 as an essential step in heterochromatin formation and maintenance.
H3K14ub stimulation is conserved in mammalian SUV39H2
To investigate if stimulation of H3K9 methylation by H3K14ub is a conserved process, we purified the KMT domains of both human G9a and SUV39H2, as well as plant SUVH4/KRYPTONITE and performed methyltransferase assays to investigate the substrate preference of these enzymes by radiometric and fluorescence-based assays (Figure 5A, Figure 5—figure supplement 1B and C). SUV39H1 and SUV39H2 are the closest human homologs of Clr4 and are tightly associated with constitutive heterochromatin. In contrast, G9a regulates H3K9me2/3 deposition in euchromatin and directly regulates gene expression. SUVH4 is one of 10 SUVH genes in Arabidopsis and is involved in the maintenance of DNA methylation. Of our candidate enzymes, SUV39H2, consistently showed stimulation by H3K14ub, while we could not detect stimulation for G9a or SUVH4. These results suggest that the H3K14ub-mediated stimulation is conserved in a subset of SUV39 family proteins in higher eukaryotes, and that H3K14ub is potentially implicated in the regulation of H3K9 methylation in mammalian heterochromatin.
Discussion
Histone methylation and acetylation are key post-translational modifications involved in the regulation of chromatin accessibility and transcription. We show here that the poorly understood histone modification H3K14ub leads to the strong stimulation of the methyltransferase activity of the fission yeast SUV39 family enzyme Clr4 through a sensing mechanism in its catalytic domain. The H3K14ub modification has, therefore, the biochemical potential to direct H3K9 methylation by Clr4 with high specificity. These results support a recently proposed regulatory role for H3K14ub in heterochromatin formation (Oya et al., 2019). Based on these and our findings, we propose that H3K14ub serves as a licensing mechanism for the deposition of H3K9me2/3 marks and heterochromatin formation (Figure 5B). In this model, acetylation of H3K14, which is strongly associated with euchromatin, assumes a much more specific role than normally associated with histone acetylation: it serves to prevent ubiquitination of H3K14 similar to a protecting group in organic chemical synthesis. The H3K14ac-specific HDACs Clr3 and Sir2 therefore indirectly control H3K9me2/3 by regulating the availability of H3K14 for ubiquitination by the CLRC ubiquitin ligase complex. Consistent with this model, mutations of H3K14 lead to loss of H3K9 methylation (Mellone et al., 2003), which is phenocopied by double deletion of Clr3 and Sir2 (Alper et al., 2013). The genetic data further suggests that ubiquitination is essential for the establishment of heterochromatin since mutants of the CLRC complex are completely devoid of heterochromatin (Horn et al., 2005; Jia et al., 2005; Li et al., 2005; Nakayama et al., 2001; Thon et al., 2005). In the clr4-GS253 and clr4-3FA mutants, we observe a severe loss of heterochromatin, which validates that recognition of H3K14ub by Clr4 is critically important for the establishment and maintenance of heterochromatin. In a parallel study, Shan et al. have shown that H3K14ub and the Clr4 UBR are critical factors for sequestration of Clr4 on H3K9M mutants, which mimic oncogenic histone lysine-to-methionine mutants (Shan et al., 2021).
We have consistently observed H3K14ub-mediated stimulation in reactions with the isolated KMT domain of Clr4, and we have shown that the N-terminus of Clr4 is not essential for recognition of H3K14ub. Our observations are difficult to reconcile with the results reported in Oya et al., 2019, where an intact N-terminus was required for stimulation. We speculate that this discrepancy is caused by the use of different protein expression constructs. While Oya et al. use GST-tagged proteins, we have removed the tags for all constructs. An influence on biochemical activity by the GST-tag has been observed previously for SUV39H2 (Schuhmacher et al., 2015), and it is therefore conceivable that a GST-tag also influences the behavior of Clr4. Shan et al. have confirmed independently that H3K14ub stimulates the isolated KMT domain as observed here (Shan et al., 2021).
The HDX-MS and mutagenesis results provide a valuable set of constraints that we employed to predict the atomic structure of the Clr4-H3K14ub complex by protein-protein docking using HADDOCK (Dominguez et al., 2003). The resulting model (Figure 5C) shows that the H3K14ub peptide can comfortably fit the H3 peptide into the binding groove and simultaneously cover Clr4 residues 253–256 mutated in Clr4-GS253 and F427 mutated in Clr4-3FA with its K14-linked ubiquitin moiety. The model suggests that F310, which is also mutated in Clr4-3FA, is unlikely to be involved directly in the interaction with ubiquitin as the N-terminal residues 253–256 shield F310 from interacting with ubiquitin (Figure 5D).
Comparing this structural model with the best understood regulatory mechanisms of other SET domain proteins reveals interesting parallels and suggests a mechanism for how H3K14ub might stimulate the enzymatic activity of Clr4. EZH2, the catalytic subunit of the polycomb repressive complex PRC2, which methylates H3K27, is stimulated through conformational changes triggered by binding of its own H3K27me product to the neighboring subunit EED in the PRC2 complex (Jiao and Liu, 2015; Justin et al., 2016; Brooun et al., 2016). Intriguingly, the UBR that we identified by HDX-MS corresponds structurally to the ‘SET activation loop’ (SAL) observed in PRC2 structures, which stabilizes the SET-I domain, a region inside the SET domain that provides one side of the peptide binding groove (Qian and Zhou, 2006; Wu et al., 2010). In a further example, the mixed lineage leukemia complexes with catalytic SET domain proteins MLL1 or MLL3 are activated by structural rearrangement of the regulatory subunits RBBP5 and ASH2L. RBBP5 thereby assumes a similar position as the SAL in PRC2 and stabilizes the SET-I domain in coordination with ASH2L (Li et al., 2016). In these experiments, apo MLL3 and the MLL3-ASH2L-RBBP5 complex showed little structural differences in the positioning of SET-I, but complex formation significantly reduced the conformational dynamics of the SET-I domain. The structural similarities between these regulatory mechanisms identify the SET-I domain as a potential mediator of the H3K14ub stimulation. In our model, the SET-I domain is sandwiched between the H3K14-linked ubiquitin moiety and the H3 peptide bound in the active site (Figure 5C and D). Therefore, binding of H3K14ub to Clr4 might not only contribute to high substrate affinity but also to an optimally organized active site.
Bioinformatic analysis of the UBR sequence using Hidden Markov Models suggests that the UBR sequence is well conserved in the Ascomycetes clade of fungi, which includes the Neurospora crassa Dim-5 protein for example (Figure 5—figure supplement 1D and E). Comparing the Clr4 motif with motifs obtained using homologous sequences from human SUV39H2, human G9a, and Arabidopsis SUVH4 shows that SUV39H2’s motif is similar to Clr4, while G9a and SUVH4 diverge significantly. This is consistent with our observation that H3K14ub can stimulate SUV39H2, but not G9A or SUVH4. The function of H3K14ub in human cells is poorly understood, but it has been detected by mass spectrometry (Harrison et al., 2016). In combination with the conservation of the H3K14ub mediated stimulation in SUV39H2, this suggests that the crosstalk observed in fission yeast is conserved in mammals. Further studies will be required to determine the prevalence of the crosstalk between H3K9me2/3 and H3K14ub, and the role it plays in the biology of eukaryotic organisms.
Materials and methods
Generation of S. pombe strains
Request a detailed protocolS. pombe strains were grown and manipulated using standard techniques unless differently stated. The list of strains used in this study can be found in the Key resources table. Strains were crossed and analyzed either by tetrad dissection or random spore analysis. For tetrad dissection, only spores from complete tetrads were further analyzed. In the absence of any selection markers, the genotype was analyzed by colony PCR. The strain clr4Δ::kanMX was generated by transforming a wild-type strain with a DNA fragment containing the kanMX cassette plus 80 bp of sequence found upstream and downstream of clr4 ORF. The transformants were selected on YES plates containing 100 µg/ml of G418, and correct integration was confirmed by PCR and sequencing. The strain 3xFLAG-clr4 was generated by CRISPR mutagenesis of a wild-type strain using as homologous repair template a fragment of DNA containing the 3xFLAG tag plus 500 bp of homology flanking the cleavage site and a silent point mutation disrupting the PAM site. The 3xFLAG-clr4-GS253 strain was generated with the same approach by transforming the 3xFLAG-clr4 strain with DNA containing the desired mutations and homologies. The sgRNA sequences used for CRISPR are indicated in Supplementary file 1C. The clr4Δ::rpl42-natMX strain was generated by transformation of a rpl42.sP56Q host strain with a DNA fragment containing the rpl42-natMX cassette plus 500 bp of homology for sequences flanking clr4 ORF. The integrants were selected by rounds of positive selection on rich media containing 100 µg/ml nourseothricin, and followed by negative selection on plates containing 100 µg/ml cycloheximide. The integrations were tested by PCR and sequencing. To generate the 3xFLAG-clr4-3FA mutant strain, the clr4Δ::rpl42-natMX strain was transformed with a DNA fragment containing the desired mutations and 500 bp homology arms outside the rpl42-natMX cassette. The strain was selected with rounds of positive and negative selection on YES plates containing cycloheximide and nourseothricin, respectively. The strain rik1-13xmyc-kanR was generated by transformation of a wild-type strain with a 13xmyc-kanR cassette with homology for the 3'end of rik1. The strain was selected on plates containing G418. All transformants were analyzed by PCR, sequenced, and the strains were outcrossed at least once. For plate-based silencing assays, reporter strains were generated by crossing the strains of interest with strains containing the marker genes ura4+ and ade6+ integrated at the centromere I (imr1L::ura4+, otr1R::ade6+).
Silencing assays in S. pombe
Request a detailed protocolStrains were grown overnight at 30°C in YES to a density of 5 × 106 cells/ml. Cells were harvested by centrifugation and resuspended in water. 10-fold serial dilutions were spotted onto PMG, PMG-Ura, and PMG +1 mg/ml 5-FOA, with 1 × 104 cells in the highest density spots.
RT-qPCR
Request a detailed protocolRT-qPCR was performed as previously described (Leopold et al., 2019). Briefly, strains were grown in YES to a density of 5 × 106 cells/ml, harvested by centrifugation, washed with water, and frozen at –70°C until further usage. RNA was purified by Trizol extraction followed by phenol/chloroform extraction and ethanol precipitation. Contaminant DNA was removed by DNaseI treatment and phenol/chloroform extraction followed by ethanol precipitation. cDNA was generated with Evo-Script Universal cDNA Master kit (Roche) and qPCR performed using SYBR Green I Master kit on LightCycler 480 instrument (Roche). qPCR primers used in this study are listed in Supplementary file 1A, and act1+ was used for normalization. Data were analyzed using the ΔΔCt method (Yuan et al., 2006).
Chromatin immunoprecipitation (ChIP)
Request a detailed protocolChIP was performed as previously described with some modifications (Leopold et al., 2019). Briefly, cells for ChIP were grown in 50 ml of YES medium to a density of 2 × 107 cells/ml. For RNA Polymerase II, H3K9me1, H3K9me2, and H3K9me3, cells were fixed in 1% formaldehyde for 30 min. For 3xFLAG-Clr4, a dual-crosslinking approach was employed as previously described (Tian et al., 2012). Cells were incubated at 18°C for 2 hr, resuspended in 5 ml PBS, and crosslinked at room temperature (RT) with 1.5 mM ethylene glycol bis-succinimidyl succinate (EGS, Thermo Scientific). After 30 min of incubation, 1% formaldehyde was added and cells were crosslinked for further 30 min. The residual formaldehyde was quenched with 125 mM glycine for 5 min. Cells were washed and stored at –70°C until further usage.
Pellets were resuspended in ChIP Buffer (50 mM HEPES/KOH pH 7.6, 150 mM NaCl, 1 mM EDTA, 1% Triton X-100, 1 mM PMSF, cOmplete EDTA free (Roche)), and lysed by bead-beating. Chromatin was enriched by centrifugation, resuspended in lysis buffer, and sonicated for 15 min (30 s/30 s ON/OFF) in a Bioruptor Pico. Cell debris was removed by centrifugation. For each IP, 50 μl of sheared soluble chromatin was diluted with 450 μl ChIP Buffer, mixed with 1 μg of antibody and incubated for 2 hr, followed by 45 min incubation with 30 µl of Dynabeads Protein A or G. For H3K9me3, 0.75 μg of anti-H3K9me3 (Diagenode C15500003) was pre-incubated with 30 µl of Dynabeads MyOne Streptavidin C1, followed by blocking with 5 μM biotin, according to the manufacturer’s instructions. The loaded beads were then added to the diluted chromatin and incubated for 2 hr.
The bead/protein complex was washed three times with ChIP Buffer, once with high salt buffer (50 mM HEPES/KOH pH 7.6, 500 mM NaCl, 1 mM EDTA), once with LiCl buffer (5 mM Tris-Cl pH 8, 250 mM LiCl, 0.5% Triton X-100, 0.05 % Tween 20), and once with TE (10 mM Tris-Cl pH 8, 1 mM EDTA). The protein-DNA complex was eluted in elution buffer (50 mM Tris-Cl pH 8, 10 mM EDTA, 1% SDS) at 65°C for 15 min, and the crosslinking reversed overnight at 65°C. The samples were then treated with proteinase K, and DNA was purified by phenol-chloroform. qPCR was performed using the primers given in Supplementary file 1A; act1+ was used as internal control.
Co-immunoprecipitation from S. pombe cells
Request a detailed protocolCells were grown overnight at 30°C to a density of 1.0 x 107 cells/ml, washed once with PBS, and the pellet frozen at –70°C until further usage. Pellets were resuspended in co-IP buffer (50 mM HEPES/KOH pH 7.5, 300 mM NaCl, 1 mM EDTA, 5 mM CHAPS, 10 mM DTT) supplemented with cOmplete EDTA free (Roche). The cells were mixed with the same volume of 0.5 mm glass beads and ruptured by bead beating in a Fastprep24 (MP Biomedical). Cell debris was removed by centrifugation and the supernatant incubated with Myc-Trap magnetic beads (ChromTek) for 1 hr at 4°C. The unbound fraction was removed, and the beads were washed five times with co-IP buffer. The bound proteins were eluted by boiling in PGLB (125 mM Tris-HCl, 4% SDS, 50% glycerol, 0.2% Orange G) supplemented with 100 mM DTT. The proteins of interest were analyzed by Western blot.
Western blotting
Request a detailed protocolSamples were run on 10% Bis-Tris gels and transferred to nitrocellulose membrane (Bio-Rad). Proteins were detected by Western blotting using antibodies against either tag or protein of interest, followed by incubation with secondary antibodies labeled with the DyLight fluorophores, and scanning with the Odyssey Imaging System (LI-COR).
Protein expression in Escherichia coli
Request a detailed protocolAll Clr4 protein constructs were cloned as N-terminal fusions to a 6xHIS-SUMO-tag in pSMT3 using Gibson cloning (Gibson et al., 2009; Mossessova and Lima, 2000). All constructs were expressed in the Rosetta2(DE3) strain. Cells were grown at 37°C to log-phase in 2xTY media supplemented with 2% ethanol, cooled down on ice, and induced with 0.4 mM IPTG at 16°C for 18 hr. Pellets were resuspended in NiNTA wash buffer (50 mM Tris-HCl pH 8.0, 300 mM NaCl, 20 mM imidazole pH 7.5, 1 mM β-mercaptoethanol) supplemented with protease inhibitor cocktail (PIC) (2 μM pepstatin, 2 µM leupeptin, 1 mM PMSF, 1 mM Benzamidine, 1 µg/ml aprotinin) and frozen at –80°C until further usage.
Protein purification from E. coli
Request a detailed protocolFor purification, frozen cells were quickly thawed at 37°C, supplemented with PIC, and ruptured by using either sonication or the Emulsiflex homogenizer. Cell debris was removed by centrifugation. The cleared lysate was filtered through a 0.45 µm filter and loaded on a HisTrap HP column (GE Healthcare). The full-length Clr4 was enriched from the cleared lysate by fractionation with 40% NH4SO4. The pellet containing Clr4 was then resuspended in NiNTA wash buffer (50 mM Tris-HCl pH 8.0, 300 mM NaCl, 20 mM imidazole pH 7.5, 1 mM β-mercaptoethanol), filtered through a 0.45 µm filter, and loaded on a HisTrap HP column. After washing with NiNTA wash buffer, the proteins were eluted by either gradient or step elution in elution buffer (50 mM Tris-HCl pH 8.0, 300 mM NaCl, 300 mM imidazole pH 7.5, 1 mM β-mercaptoethanol). All proteins containing a 6xHIS-SUMO-tag were treated with Ulp1 while being dialyzed overnight in cleavage buffer (5 mM Tris-HCl pH 8.0, 150 mM NaCl, 2 mM β-mercaptoethanol). SUV39H2 and G9a were treated with thrombin to remove the 6xHIS-tag and dialyzed into cleavage buffer. The cleaved tags and Ulp1 were removed by negative purification on a HisTrap HP column. All proteins were further purified by SEC (Superdex75 or Superdex200, GE Healthcare) in gel filtration buffer (5 mM Tris-HCl pH 8.0, 150 mM NaCl, 5 mM DTT). Peak fractions were concentrated as needed in Amicon Ultra spin concentrators. For methyltransferase assays, the concentrated proteins were mixed 1:1 with glycerol, frozen in liquid nitrogen, and stored at –80°C until further usage. Clr4FL, SUV39H2, and SUVH4 were frozen without glycerol.
Materials for peptide synthesis
Request a detailed protocolAll solvents and reagents were purchased from commercial sources and used without further purification. All amino acid derivatives, 2-chlorotrityl chloride (2 Cl Trt) resin and 2-(7-Aza-1H-benzotriazole-1-yl)–1,1,3,3-tetramethyluronium hexafluorophosphate (HATU) were purchased from Novabiochem, Merck (Darmstadt, Germany). N,N-Dimethylformamide (DMF), N,N-diisopropylethylamine (DIEA), and piperidine were from Acros Organics (Geel, Belgium). O-(6-Chlorobenzotriazol-1-yl)-N,N,N′,N′-tetramethyluronium hexafluorophosphate (HCTU) was from Carl Roth GmbH (Karsruhe, Netherlands). Hydrazine monohydrate was purchased from Alfa Aesar (Heysham, UK), acetonitrile (ACN) from Avantor Performance Materials (USA). Tentagel S RAM resin, hydroxybenzotriazole (HOBt), tris(2-carboxyethyl)phosphine (TCEP), diethylether, phenylsilane, tetrakis(triphenylphosphine)palladium(0), trifluoroacetic acid (TFA), dichloromethane (DCM), triisopropylsilane (TIS), L-glutathione reduced (GSH), sodium diethyldithiocarbamate trihydrate, and methyl thioglycolate (MTG) were from Sigma-Aldrich (Taufkirchen, Germany). 2,2'-Azobis[2-(2-imidazolin-2-yl)propane]dihydrochloride (VA-044) was purchased from Wako Pure Chemical Industries, Ltd. (Osaka, Japan). All other commonly used chemical reagents and buffer components were from Applichem (Darmstadt, Germany) and Fisher Scientific (Reinach, Switzerland).
Instrumentation for peptide synthesis and purification
Request a detailed protocolReaction vessels for manual peptide synthesis as well as the automated Tribute peptide synthesizer were from Protein Technologies Inc. Analytical RP-HPLC analysis was performed on an Agilent 1260 series instrument using an analytical Agilent Zorbax C18 column (column dimensions: 150 mm × 4.6 mm, 5 μm particle size) at a flow rate of 1 ml/min. All RP-HPLC analyses were done with 0.1% (v/v) TFA in H2O (RP-HPLC solvent A) and 90% ACN and 0.1% (v/v) TFA in H2O (RP-HPLC solvent B) as mobile phases. Typically, a gradient from 0% to 70% solvent A to solvent B over 30 min was used for analytical RP-HPLC analyses unless otherwise stated. Purification of proteins on a semi-preparative scale were performed on an Agilent 1260 series instrument using a semi-preparative Agilent Zorbax C18 column (column dimensions: 250 mm × 9.4 mm, 5 μm particle size) at a flow rate of 4 ml/min. Lyophilization was carried out with a Telstar LyoQuest freeze dryer. Electrospray ionization mass spectrometric (ESI-MS) analysis was conducted on a Shimadzu MS2020 single quadrupole instrument connected to a Nexera UHPLC system. Mass spectra were acquired by electrospray ionization in positive ion mode in the mass range of 200–2000 m/z.
Preparation of Fmoc-Aa-hydrazine-Cl-trityl-resin
Request a detailed protocolPreparation of Fmoc-Arg(Pbf)-NHNH-Cl-Trt-resin was performed as reported previously (Stavropoulos et al., 1996). In general, 0.50 g 2-Cl-Trt-resin (substitution: 1.63 mmol/g, 1.00 eq., 0.82 mmol) were swollen in 3.00 ml DMF for 30 min at RT. Subsequently, the resin was allowed to cool to 0°C and 1.00 ml of a solution containing DIEA (3.00 eq., 2.45 mmol, 427 μl) and hydrazine monohydrate (2.00 eq., 1.64 mmol, 80 μl) in DMF was added dropwise. The reaction mixture was stirred 1 hr at RT. 100 μl methanol (MeOH) was added, the resin was stirred 10 min at RT, transferred to a reaction vessel for manual peptide synthesis, and washed thoroughly with DMF.
Due to the low stability of the hydrazine resin, the first amino acid Fmoc-Arg(Pbf)-OH was coupled manually by standard Fmoc chemistry (Atherton and Sheppard, 1989). The loading of the 2-Cl-Trt-resin with hydrazine was assumed to be 50% (1.00 eq., 0.41 mmol). Fmoc-Arg(Pbf)-OH (5.00 eq., 2.05 mmol) was activated with 3.90 ml of a 0.50 M HATU solution (4.76 eq., 1.95 mmol) in DMF followed by 2 min incubation at RT. Then, 714 μl DIEA (10.00 eq., 4.10 mmol) were added and the reaction mixture was incubated another 1 min at RT. The activated amino acid was added to the resin, incubated 30 min at RT, and washed with DMF. To ensure high yield, the coupling procedure was repeated. Finally, the resin was washed with DMF, DCM, and MeOH and dried under vacuum. Resin substitution was determined by treating a defined amount of resin with 20% piperidine in DMF for 30 min at RT, followed by spectrophotometric quantification of released dibenzofulven-related species (resin loading: 0.51 mmol/g).
Automated solid phase peptide synthesis (SPPS)
Request a detailed protocolThe sequences of the two peptides synthesized by Fmoc-SPPS are summarized in Supplementary file 1B. The synthetic procedures and analytical data are presented in the following. General protocol: the peptides were synthesized by the Tribute peptide synthesizer (PTI) on the previously prepared Fmoc-Arg(Pbf)-hydrazine-Cl-Trt-resin (peptide P1) or on Tentagel S RAM resin (peptide P2) yielding peptides with C-terminal hydrazide (peptide P1) or amide (peptide P2). The syntheses were performed on 0.1 mmol scale using Fmoc chemistry. The following base-resistant groups were employed to protect amino acid side chains: Arg(Pbf), Thr(tBu), Lys(Boc), Gln(Trt), Asn(Trt), Ser(tBu). Of note, Lys14 (peptide P1) and Lys18 (peptide P2) were protected with Alloc protecting group. To maximize synthesis yield, pseudoproline dipeptide building blocks were used and amino acids were double coupled where necessary.
Briefly, the N-terminal Fmoc-group was deprotected with 20% (v/v) piperidine in DMF. Activation of amino acid (5.00 eq., 0.50 mmol) was achieved by addition of HCTU (4.76 eq. 0.48 mmol) and DIEA (10.00 eq., 1.00 mmol). The coupling step was performed by adding the activated amino acid to the resin, followed by 30 min incubation at RT. When the full-length peptides were assembled, the peptidyl-resin was washed with DMF, DCM, and MeOH and dried under vacuum. Importantly, the N-terminal Fmoc protecting group was not removed at this stage to allow further manipulation of the peptide.
Alloc deprotection
Request a detailed protocolThe peptidyl-resin was swollen for 30 min in DCM. Alloc deprotection of Lys14 (peptide P1) or Lys18 (peptide P2) was initiated by addition of 1.00 ml of dry DCM and PhSiH3 (24 eq., 24 mmol), followed by Pd(PPh3)4 (0.25 eq., 0.025 mmol) in 3.00 ml dry DCM. The peptidyl-resin was incubated for 30 min at RT and washed with DCM. The deprotection reaction with PhSiH3 and Pd(PPh3)4 was repeated two more times. The resin was thoroughly washed with DCM followed by washing with 0.5% (v/v) DIEA in DMF; 0.5% (w/v) sodium-diethyldithiocarbamate in DMF; 50% (v/v) DCM in DMF; 0.5% (w/v) HOBt in DMF and extensively washed with DMF.
Manual coupling of Boc-Cys(Trt)-OH
Request a detailed protocolThe cysteine residue was coupled manually to the ε-amino group of Lys14 (peptide P1) or Lys18 (peptide P2), yielding peptides P1′ and P2′ respectively (Figure 1—figure supplement 1A). Boc-Cys(Trt)-OH (5 eq., 0.50 mmol) was activated by addition of 0.95 ml of a 0.50 M HATU solution (4.76 eq., 0.48 mmol) in DMF, followed by 2 min incubation at RT. 172 μl DIEA (10 eq., 1 mmol) were added and the reaction mixture was incubated another 1 min at RT. The activated amino acid was added to the peptidyl-resin, incubated 30 min at RT, and washed with DMF. To ensure high yield, the coupling procedure was repeated. Finally, the resin was washed thoroughly with DMF.
N-terminal Fmoc deprotection and cleavage from the resin
Request a detailed protocolFmoc deprotection was achieved by treating the peptidyl-resin with 10.00 ml 20% (v/v) piperidine in DMF for 5 min. The deprotection step was repeated.
The peptides were cleaved from the resin using either 95% (v/v) TFA, 2.5% (v/v) TIS, and 2.5% (v/v) H2O. The crude peptides were precipitated by addition of ice-cold diethyl ether, recovered by centrifugation, dissolved in 50% (v/v) acetonitrile in H2O, flash-frozen, and lyophilized.
Preparation of ubiquitin-MESNa
Request a detailed protocolUbiquitin (G76C) was expressed in Rosetta2(DE3) cells as a C-terminal fusion to the N-terminal half of the Npu split intein, followed by a hexahistidine tag (Kilic et al., 2018). Bacterial cells were resuspended in lysis buffer (50 mM Na-phosphate buffer pH 7.0, 300 mM NaCl, 40 mM imidazole pH 7.0) and ruptured by sonication. Ubiquitin-NpuN-6xHIS was purified by Ni-affinity chromatography and eluted with 300 mM imidazole. The ubiquitin thioester was generated by intein derivatization with 100 mM sodium 2-mercaptoethanesulfonate (MESNa). The reaction was allowed to proceed for 16 hr at RT. Imidazole was removed by diafiltration and ubiquitin-MESNa was further purified by negative Ni-affinity chromatography to remove the NpuN-6xHIS tag. Finally, ubiquitin-MESNa was diafiltrated against 0.1% TFA, lyophilized, and stored at –20°C.
One-pot ligation and desulfurization
Request a detailed protocolUbiquitin-MESNa (1.00 eq., 1.16 μmol) was dissolved in 115.80 μl ligation buffer (6.00 M GmdHCl, 0.20 M sodium phosphate, pH 7, degassed) to a final concentration of 10.00 mM and 0.77 μl MTG (7.50 eq., 8.70 μmol) were added to the solution. The thiol-thioester exchange reaction was allowed to proceed for 20 min at RT. Ligation between ubiquitin-MTG thioester and the cysteine residue coupled to Lys14 (peptide P1) or Lys18 (peptide P2) was initiated by addition of peptide P1′ or peptide P2′ (1.30 eq., 1.51 μmol), followed by addition of 5.80 μl of 0.50 M TCEP solution in ligation buffer. The ligation mixture was incubated at 25°C for 16 hr. The progress of the reaction was monitored by RP-HPLC and ESI-MS analysis.
When the ligation was complete, radical desulfurization of the cysteine at the ligation site was performed in the same reaction tube without prior purification of the ligation product. TCEP desulfurization buffer (0.50 M TCEP, 6.00 M GdmHCl, 0.20 M phosphate, pH 7) was added to a final TCEP concentration of 0.25 M. The desulfurization reaction was initiated by addition of VA-044 and GSH to a final concentration of 30 mM and 40 mM, respectively. The reaction mixture was incubated at 42° C for 6 hr, and the progress of the reaction was monitored by RP-HPLC and ESI-MS analysis, yielding ubiquitin adducts P1′′ and P2′′ (Figure 1—figure supplement 1).
Purification of H3K14ub (P1′′) and H3K18ub (P2′′)
Request a detailed protocolThe ubiquitin adduct peptides P1′′ and P2′′ were purified by semi-preparative RP-HPLC on a linear gradient from 25% to 55% solvent B over 40 min. Pure fractions were pooled, lyophilized, and analyzed by analytical RP-HPLC and ESI-MS (Figure 1—figure supplement 1B and C).
Tritium-based histone methyltransferase assay
Request a detailed protocolThe HMT assays were performed in a total volume of 10 µl by mixing the substrates with the indicated constructs in reaction buffer (25 mM sodium phosphate pH 8.0, 100 mM NaCl, 10% v/v glycerol, 2.9 µM (0.5 µCi) S‐[Methyl‐3 H]‐adenosyl‐L-methionine [1 mCi/mL, Perkin Elmer]) and incubated at 30°C. Measurements were performed with 20 µM of substrate, and 20 nM or 200 nM of each enzyme. Assay-specific concentrations and incubation times are indicated in the text and captions. For gel-based assays, the reactions were stopped by addition of 5 µl of PGLB, 100 mM DTT and the unreacted SAM was removed by SDS-PAGE. The gels were stained with Coomassie Brilliant Blue, imaged, and incubated for 30 min in EN3HANCE (Perkin Elmer) prior to drying. The methylated histones were detected by fluorography.
For filter binding assays, the reactions were stopped by addition of 5 µl 75% acetic acid, and 10 µl of each mixture was spotted onto phosphocellulose filter paper disks. Each filter was washed for 5 min in 2 ml 100 mM sodium bicarbonate pH 9.0 for a total of three times. The activity on each disk was quantified in a liquid scintillation counter after 1 hr incubation in 3 ml ULTIMA GOLD F cocktail.
TR-FRET-based methyltransferase assay
Request a detailed protocolThe assays were performed using the EPIgeneous Methyltransferase kit (Cisbio). The reactions were performed by mixing the Clr4 constructs (30 nM) with the indicated substrates (1 nM –-200 µM for Figure 1D, 2 µM in Figure 1E, and 20 µM in Figure 5—figure supplement 1B) in Mtase Buffer (25 mM sodium phosphate pH 8.0, 100 mM NaCl, 10% v/v glycerol) supplemented with 10 µM SAM and incubated at 30°C for 15 min, unless differently stated. The reactions were stopped and measured using a TECAN SPARK plate reader. Detection was performed at 620 and 665 nm wavelength. The 665/620 signal ratio was used to extrapolate the concentration of SAH from a standard curve.
Bioluminescence-based methyltransferase assay
Request a detailed protocolMethyltransferase assays were performed with wild-type and mutant Clr4KMT domains on H3K14ub (using 20 nM Clr4) and unmodified H3 (using 1 µM Clr4) in Mtase Buffer supplemented with SAM (10 µM for H3K14ub, 100 µM for H3). Reactions were stopped by addition of TFA to 0.1% and the production of SAH was measured using the MTase-Glo Methyltransferase Assay kit (Promega) according to the manufacturer’s instructions using a Hidex Sense microplate reader.
Hydrogen deuterium exchange mass spectrometry (HDX-MS)
Request a detailed protocolHDX reactions were done in 50 µl volumes with a final protein concentration of 3 µM of enzyme. The dynamics of Clr4 was compared in three different conditions: (i) soluble Clr4, (ii) complex of Clr4 bound to ubiquitinated H3 peptide, and (iii) Clr4 in the presence of an excess of H3 peptide. Dynamics of ubiquitinated H3 peptide was investigated by comparing its D2O incorporation alone or when in complex with Clr4. Protein or protein complexes were prepared in 5 µl volume, and deuterium exchange was initiated by addition of 45 µl of deuterated buffer (10 mM Tris-HCl pH 7.4/150 mM NaCl/10 mM DTT). For the condition of Clr4 with an excess of H3 peptide, the deuterated buffer was supplemented with H3 peptide (from lyophilized powder). Exchange was carried out at RT for four timepoints (3 s, 30 s, 300 s, 3000 s) and terminated by the addition of 20 µl ice-cold quench buffer (3 M guanidine – HCl/100 mM NaH2PO4 pH 2.5/1% formic acid). All experiments were repeated in duplicates. Samples were immediately frozen in liquid nitrogen and stored at –80°C for up to 2 weeks.
Protein samples were thawed and injected in a UPLC system immersed in ice. The protein was digested via two immobilized pepsin columns (Thermo #23131), and peptides were collected onto a VanGuard precolumn trap (Waters). The trap was subsequently eluted and peptides separated with a C18, 300 Å, 1.7 µm particle size Fortis Bio column 100 × 2.1 mm over a gradient of 8–30%B over 18 min at 90 µl/min (Buffer A: 0.1% formic acid; buffer B: 100% acetonitrile/0.1% formic acid). Mass spectra were acquired on an Orbitrap Velos Pro (Thermo), for ions from 400 to 2200 m/z using an electrospray ionization source operated at 160°C, 5 kV of ion spray voltage. Peptides were identified by data-dependent acquisition after MS/MS, and data were analyzed by Mascot. Deuterium incorporation levels were quantified using HD examiner software (Sierra Analytics), and the quality of every peptide was checked manually. Results are presented as a percentage of maximal deuteration, with a theoretical maximal deuteration level of 86.3%. All experimental details and data of percentage deuterium incorporation for all peptides can be found in Supplementary files 2–4. Differences in exchange in a peptide were considered significant if they met all three of the following criteria: (1) for Clr4: ≥7% change in exchange, ≥ 0.6 Da difference in exchange, and a p-value<0.05 using a two-tailed Student’s t-test; and (2) for ubiquitin: ≥9% change in exchange, ≥1.3 Da difference in exchange, and a p-value<0.05 using a two-tailed Student’s t-test.
Isothermal titration calorimetry
Request a detailed protocolITC experiments were performed at 23°C using a MicroCal VP-ITC calorimeter (Malvern Panalytical). All proteins and peptides used were dialyzed overnight against ITC buffer (25 mM sodium phosphate pH 8.0, 100 mM NaCl, 10% v/v glycerol, 0.5 mM TCEP) prior to experiments. 10 μl of the indicated Clr4 constructs at 50 µM were injected in 180 s time intervals in the cell containing the H3K14ub peptide at 5 μM. After subtracting heat enthalpies for titrations of the respective proteins into buffer, the ITC data were analyzed using the Origin software provided by the manufacturer.
Molecular docking
Request a detailed protocolMolecular docking was performed using Clr4KMT192-490 (PDBID:1MVH), ubiquitin (PDBID:1UBQ), and a peptide comprising residues 7–17 of histone H3. The input peptide was derived by extension of the H3 peptide bound to N. crassa DIM-5 (PDBID:1PEG). Docking was performed with a local version of HADDOCK 2.0 using the interaction interface identified by HDX-MS as restraint.
Crystallization of Clr4-3FA
Request a detailed protocolClr4-3FA was purified as described above. Clr4-3FA at 20 mg/ml was mixed with H3K14ub peptide and S-adenosyl-homocysteine to a ratio of 1:1.2:10 (protein:peptide:SAH) and diluted to a concentration of 10 mg/ml. Protein was diluted 1:1 with reservoir solution (71.4 mM MES, 28.6 mM imidazole, 20% PEG 10 K, 30 mM magnesium acetate, 6.6% v/v MPD, 6.6% v/v PEG 1000, 6.6% v/v PEG 3350) and incubated at 18°C in sitting-drop vapor diffusion setup. Crystals were directly frozen in liquid nitrogen.
Crystallographic data collection and structure determination
Request a detailed protocolX-ray diffraction data were collected at the Diamond Light Source beamline I04 (Didcot, Oxford, UK) at a wavelength of 0.97950 Å. The datasets were analyzed using XDS and the CCP4 suite (Kabsch, 2010; Winn et al., 2011). Data were corrected for anisotropy using the STARANISO server (http://staraniso.globalphasing.org). Molecular replacement was performed using MOLREP and the coordinates from PDBID 6BOX (chain A) as search model (Vagin and Teplyakov, 2010). Refinement was performed using REFMAC, and Coot was used to build the Clr4-3FA model (Emsley et al., 2010; Murshudov et al., 2011). Structural alignments and figures were generated in PyMOL Molecular Graphics System, version 1.8, Schrödinger, LLC. Data collection and refinement statistics are presented in Table 2.
Accession codes
Request a detailed protocolCoordinates and structure factors have been deposited in the Protein Data Bank under accession code 6Z2A.
Data availability
Coordinates and structure factors have been deposited in the Protein Data Bank under accession code 6Z2A.
-
RCSB Protein Data BankID 6Z2A. Structure of Clr4 mutant - F256A/F310A/F427A bound to SAH.
References
-
BookSolid Phase Peptide Synthesis: A Practical Approach, Practical Approach SeriesOxford, England ; New York: IRL Press.
-
Functional divergence between histone deacetylases in fission yeast by distinct cellular localization and in vivo specificityMolecular and Cellular Biology 22:2170–2181.https://doi.org/10.1128/MCB.22.7.2170-2181.2002
-
HADDOCK: a protein-protein docking approach based on biochemical or biophysical informationJournal of the American Chemical Society 125:1731–1737.https://doi.org/10.1021/ja026939x
-
Features and development of CootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
Histone methylation: A dynamic mark in health, disease and inheritanceNature Reviews. Genetics 13:343–357.https://doi.org/10.1038/nrg3173
-
RNAi-dependent formation of heterochromatin and its diverse functionsCurrent Opinion in Genetics & Development 20:134–141.https://doi.org/10.1016/j.gde.2010.02.003
-
Heterochromatin: Guardian of the GenomeAnnual Review of Cell and Developmental Biology 34:265–288.https://doi.org/10.1146/annurev-cellbio-100617-062653
-
XDSActa Crystallographica. Section D, Biological Crystallography 66:125–132.https://doi.org/10.1107/S0907444909047337
-
Inhibitors of protein methyltransferases and demethylasesChemical Reviews 118:989–1068.https://doi.org/10.1021/acs.chemrev.6b00801
-
BookHydrogen-deuterium exchange mass spectrometry to study protein complexes inIn: Marsh JA, editors. Editor. Protein Complex Assembly. New York, NY: Springer. pp. 153–171.https://doi.org/10.1007/978-1-4939-7759-8_10
-
Structure of the SET domain histone lysine methyltransferase Clr4Nature Structural & Molecular Biology 9:828–832.https://doi.org/10.1038/nsb860
-
A two-state activation mechanism controls the histone methyltransferase Suv39h1Nature Chemical Biology 12:188–193.https://doi.org/10.1038/nchembio.2008
-
Refmac5 for the refinement of macromolecular crystal structuresActa Crystallographica. Section D, Biological Crystallography 67:355–367.https://doi.org/10.1107/S0907444911001314
-
HMMER web server: 2018 updateNucleic Acids Research 46:W200–W204.https://doi.org/10.1093/nar/gky448
-
SET domain protein lysine methyltransferases: Structure, specificity and catalysisCellular and Molecular Life Sciences 63:2755–2763.https://doi.org/10.1007/s00018-006-6274-5
-
The Chp1-Tas3 core is a multifunctional platform critical for gene silencing by RITSNature Structural & Molecular Biology 18:1351–1357.https://doi.org/10.1038/nsmb.2151
-
Activity and specificity of the human SUV39H2 protein lysine methyltransferaseBiochimica et Biophysica Acta (BBA) - Gene Regulatory Mechanisms 1849:55–63.https://doi.org/10.1016/j.bbagrm.2014.11.005
-
Two-step cross-linking for analysis of protein-chromatin interactionsMethods in Molecular Biology 809:7.https://doi.org/10.1007/978-1-61779-376-9_7
-
Molecular replacement with MOLREPActa Crystallographica. Section D, Biological Crystallography 66:22–25.https://doi.org/10.1107/S0907444909042589
-
Histone H3 lysine 14 acetylation is required for activation of a DNA damage checkpoint in fission yeastThe Journal of Biological Chemistry 287:4386–4393.https://doi.org/10.1074/jbc.M111.329417
-
Overview of the CCP4 suite and current developmentsActa Crystallographica. Section D, Biological Crystallography 67:235–242.https://doi.org/10.1107/S0907444910045749
-
Statistical analysis of real-time PCR dataBMC Bioinformatics 7:85.https://doi.org/10.1186/1471-2105-7-85
-
Roles of the Clr4 methyltransferase complex in nucleation, spreading and maintenance of heterochromatinNature Structural & Molecular Biology 15:381–388.https://doi.org/10.1038/nsmb.1406
Article and author information
Author details
Funding
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (PP00P3_139137)
- Thomas Schalch
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (PP00P3_163760_1)
- Thomas Schalch
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (PP00P3_172904)
- Thomas Schalch
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (PZ00P3_148269)
- Oscar Vadas
Biotechnology and Biological Sciences Research Council (BB/S018549)
- Thomas Schalch
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Janet Partridge, Robin Allshire, and Shiv Grewal for providing strains. We thank Steven Jacobsen, Chris Lima and Cheryl Arrowsmith for plasmids. We thank John Challiss, Jaswir Basran, and Raj Patel for support with enzymatic assays. We would like to thank Diamond Light Source for beamtime (proposal mx19880), and the staff of beamlines I04 for assistance with crystal testing and data collection. BF was supported by EPFL. OV was supported by a Swiss National Science Foundation (SNSF) Ambizione Fellowship (PZ00P3_148269). TS was supported by a SNSF Professorship (PP00P3_139137, PP00P3_163760_1, PP00P3_172904), a Biotechnology and Biological Sciences Research Council Project Grant (BB/S018549/1) and contributions by Fondation Ernst et Lucie Schmidheiny, Fonds Constantin Topali and Société Académique de Genève.
Copyright
© 2021, Stirpe et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,839
- views
-
- 311
- downloads
-
- 22
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Chromosomes and Gene Expression
- Developmental Biology
About 70% of human cleavage stage embryos show chromosomal mosaicism, falling to 20% in blastocysts. Chromosomally mosaic human blastocysts can implant and lead to healthy new-borns with normal karyotypes. Studies in mouse embryos and human gastruloids showed that aneuploid cells are eliminated from the epiblast by p53-mediated apoptosis while being tolerated in the trophectoderm. These observations suggest a selective loss of aneuploid cells from human embryos, but the underlying mechanisms are not yet fully understood. Here, we investigated the cellular consequences of aneuploidy in a total of 125 human blastocysts. RNA-sequencing of trophectoderm cells showed activated p53 pathway and apoptosis proportionate to the level of chromosomal imbalance. Immunostaining corroborated that aneuploidy triggers proteotoxic stress, autophagy, p53-signaling, and apoptosis independent from DNA damage. Total cell numbers were lower in aneuploid embryos, due to a decline both in trophectoderm and in epiblast/primitive endoderm cell numbers. While lower cell numbers in trophectoderm may be attributed to apoptosis, aneuploidy impaired the second lineage segregation, particularly primitive endoderm formation. This might be reinforced by retention of NANOG. Our findings might explain why fully aneuploid embryos fail to further develop and we hypothesize that the same mechanisms lead to the removal of aneuploid cells from mosaic embryos.
-
- Chromosomes and Gene Expression
- Developmental Biology
Transcription often occurs in bursts as gene promoters switch stochastically between active and inactive states. Enhancers can dictate transcriptional activity in animal development through the modulation of burst frequency, duration, or amplitude. Previous studies observed that different enhancers can achieve a wide range of transcriptional outputs through the same strategies of bursting control. For example, in Berrocal et al., 2020, we showed that despite responding to different transcription factors, all even-skipped enhancers increase transcription by upregulating burst frequency and amplitude while burst duration remains largely constant. These shared bursting strategies suggest that a unified molecular mechanism constraints how enhancers modulate transcriptional output. Alternatively, different enhancers could have converged on the same bursting control strategy because of natural selection favoring one of these particular strategies. To distinguish between these two scenarios, we compared transcriptional bursting between endogenous and ectopic gene expression patterns. Because enhancers act under different regulatory inputs in ectopic patterns, dissimilar bursting control strategies between endogenous and ectopic patterns would suggest that enhancers adapted their bursting strategies to their trans-regulatory environment. Here, we generated ectopic even-skipped transcription patterns in fruit fly embryos and discovered that bursting strategies remain consistent in endogenous and ectopic even-skipped expression. These results provide evidence for a unified molecular mechanism shaping even-skipped bursting strategies and serve as a starting point to uncover the realm of strategies employed by other enhancers.