Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Targeting pericentric non-consecutive motifs for heterochromatin initiation

Abstract

Pericentric heterochromatin is a critical component of chromosomes marked by histone H3 K9 (H3K9) methylation1,2,3. However, what recruits H3K9-specific histone methyltransferases to pericentric regions in vertebrates remains unclear4, as does why pericentric regions in different species share the same H3K9 methylation mark despite lacking highly conserved DNA sequences2,5. Here we show that zinc-finger proteins ZNF512 and ZNF512B specifically localize at pericentric regions through direct DNA binding. Notably, both ZNF512 and ZNF512B are sufficient to initiate de novo heterochromatin formation at ectopically targeted repetitive regions and pericentric regions, as they directly recruit SUV39H1 and SUV39H2 (SUV39H) to catalyse H3K9 methylation. SUV39H2 makes a greater contribution to H3K9 trimethylation, whereas SUV39H1 seems to contribute more to silencing, probably owing to its preferential association with HP1 proteins. ZNF512 and ZNF512B from different species can specifically target pericentric regions of other vertebrates, because the atypical long linker residues between the zinc-fingers of ZNF512 and ZNF512B offer flexibility in recognition of non-consecutively organized three-nucleotide triplets targeted by each zinc-finger. This study addresses two long-standing questions: how constitutive heterochromatin is initiated and how seemingly variable pericentric sequences are targeted by the same set of conserved machinery in vertebrates.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: ZNF512 and ZNF512B specifically localize to pericentric regions.
Fig. 2: ZNF512 and ZNF512B are sufficient to initiate heterochromatin formation through interaction with SUV39H.
Fig. 3: Distinct properties of SUV39H1 and SUV39H2.
Fig. 4: Split zinc-fingers of ZNF512 and ZNF512B recognize non-consecutive TTC triplets in pericentric repeats.
Fig. 5: ZNF512 and ZNF512B facilitate SUV39H-dependent H3K9me3 formation at LINEs.

Similar content being viewed by others

Data availability

The raw sequencing data reported in this paper have been deposited to the Genome Sequence Archive of the National Genomics Data Center, China National Center for Bioinformation/Beijing Institute of Genomics, Chinese Academy of Sciences (GSA: CRA011890) and are publicly accessible at https://ngdc.cncb.ac.cn/gsa. The processed data have been deposited in BIG OMIX under accession number OMIX004553 and are publicly accessible at https://ngdc.cncb.ac.cn/omix. Accession codes of the public data in ENCODE used in this study are as follows: ZNF512 ChIP–seq in human K562 cells, ENCSR696URR; ZNF512B ChIP–seq in human MCF-7 cells, ENCSR761LRR. The expression atlas data for Pax3, Pax9, Foxd3, Zfp512 and Zfp512b used in Extended Data Fig. 1 were downloaded from https://www.ebi.ac.uk/gxa/home. We used mouse genome version mm10 and human genome versions hg38 and CHM13 in this study. Source data for Figs. 1b, 2d, 3b,c,e,f and Extended Data Figs. 2b, 4b and 8d are provided with this paper. The gel for Extended Data Fig. 3d and uncropped film scans for Fig. 4e,f and Extended Data Figs. 3c and  8j are presented in Supplementary Fig. 1.

Code availability

The software used to analyse these data is listed in the Methods and is all publicly available.

References

  1. Allshire, R. C. & Madhani, H. D. Ten principles of heterochromatin formation and function. Nat. Rev. Mol. Cell Biol. 19, 229–244 (2018).

    Article  CAS  PubMed  Google Scholar 

  2. Janssen, A., Colmenares, S. U. & Karpen, G. H. Heterochromatin: guardian of the genome. Annu. Rev. Cell Dev. Biol. 34, 265–288 (2018).

    Article  CAS  PubMed  Google Scholar 

  3. Grewal, S. I. S. & Jia, S. Heterochromatin revisited. Nat. Rev. Genet. 8, 35–46 (2007).

    Article  CAS  PubMed  Google Scholar 

  4. Becker, J. S., Nicetto, D. & Zaret, K. S. H3K9me3-dependent heterochromatin: barrier to cell fate changes. Trends Genet. 32, 29–41 (2016).

    Article  CAS  PubMed�� Google Scholar 

  5. Saksouk, N., Simboeck, E. & Dejardin, J. Constitutive heterochromatin formation and transcription in mammals. Epigenetics Chromatin 8, 3 (2015).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Heitz, E. Das heterochromatin der moose (Bornträger, 1928).

  7. Muller, H. J. Types of visible variations induced by x-rays in Drosophila. J. Genet. 22, 299–334 (1930).

    Article  Google Scholar 

  8. Rea, S. et al. Regulation of chromatin structure by site-specific histone H3 methyltransferases. Nature 406, 593–599 (2000).

    Article  ADS  CAS  PubMed  Google Scholar 

  9. Lachner, M., O’Carroll, D., Rea, S., Mechtler, K. & Jenuwein, T. Methylation of histone H3 lysine 9 creates a binding site for HP1 proteins. Nature 410, 116–120 (2001).

    Article  ADS  CAS  PubMed  Google Scholar 

  10. Nakayama, J., Rice, J. C., Strahl, B. D., Allis, C. D. & Grewal, S. I. S. Role of histone H3 lysine 9 methylation in epigenetic control of heterochromatin assembly. Science 292, 110–113 (2001).

    Article  ADS  CAS  PubMed  Google Scholar 

  11. Peters, A. H. F. M. et al. Loss of the Suv39h histone methyltransferases impairs mammalian heterochromatin and genome stability. Cell 107, 323–337 (2001).

    Article  CAS  PubMed  Google Scholar 

  12. Peters, A. H. F. M. et al. Partitioning and plasticity of repressive histone methylation states in mammalian chromatin. Mol. Cell 12, 1577–1589 (2003).

    Article  CAS  PubMed  Google Scholar 

  13. Bannister, A. J. et al. Selective recognition of methylated lysine 9 on histone H3 by the HP1 chromo domain. Nature 410, 120–124 (2001).

    Article  ADS  CAS  PubMed  Google Scholar 

  14. Volpe, T. A. et al. Regulation of heterochromatic silencing and histone H3 lysine-9 methylation by RNAi. Science 297, 1833–1837 (2002).

    Article  ADS  CAS  PubMed  Google Scholar 

  15. Verdel, A. et al. RNAi-mediated targeting of heterochromatin by the RITS complex. Science 303, 672–676 (2004).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  16. Hall, I. M. et al. Establishment and maintenance of a heterochromatin domain. Science 297, 2232–2237 (2002).

    Article  ADS  CAS  PubMed  Google Scholar 

  17. Jia, S. T., Noma, K. & Grewal, S. I. S. RNAi-independent heterochromatin nucleation by the stress-activated ATF/CREB family proteins. Science 304, 1971–1976 (2004).

    Article  ADS  CAS  PubMed  Google Scholar 

  18. Maison, C. et al. SUMOylation promotes de novo targeting of HP1alpha to pericentric heterochromatin. Nat. Genet. 43, 220–227 (2011).

    Article  CAS  PubMed  Google Scholar 

  19. Velazquez Camacho, O. et al. Major satellite repeat RNA stabilize heterochromatin retention of Suv39h enzymes by RNA-nucleosome association and RNA:DNA hybrid formation. eLife 6, e25293 (2017).

    Article  PubMed  PubMed Central  Google Scholar 

  20. Johnson, W. L. et al. RNA-dependent stabilization of SUV39H1 at constitutive heterochromatin. eLife 6, e25299 (2017).

    Article  PubMed  PubMed Central  Google Scholar 

  21. Shirai, A. et al. Impact of nucleic acid and methylated H3K9 binding activities of Suv39h1 on its heterochromatin assembly. eLife 6, e25317 (2017).

    Article  PubMed  PubMed Central  Google Scholar 

  22. Bulut-Karslioglu, A. et al. A transcription factor-based mechanism for mouse heterochromatin formation. Nat. Struct. Mol. Biol. 19, 1023–1030 (2012).

    Article  CAS  PubMed  Google Scholar 

  23. Puri, D. et al. Foxd3 controls heterochromatin-mediated repression of repeat elements and 2-cell state transcription. EMBO Rep. 22, e53180 (2021).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Gao, X. D. et al. C-BERST: defining subnuclear proteomic landscapes at genomic elements with dCas9-APEX2. Nat. Methods 15, 433–436 (2018).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Myers, S. A. et al. Discovery of proteins associated with a predefined genomic locus via dCas9-APEX-mediated proximity labeling. Nat. Methods 15, 437–439 (2018).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Qiu, W. Q. et al. Determination of local chromatin interactions using a combined CRISPR and peroxidase APEX2 system. Nucleic Acids Res. 47, e52 (2019).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Ong, S. E., Foster, L. J. & Mann, M. Mass spectrometric-based approaches in quantitative proteomics. Methods 29, 124–130 (2003).

    Article  CAS  PubMed  Google Scholar 

  28. Saksouk, N. et al. Redundant mechanisms to form silent chromatin at pericentromeric regions rely on BEND3 and DNA methylation. Mol. Cell 56, 580–594 (2014).

    Article  CAS  PubMed  Google Scholar 

  29. Schmidtmann, E., Anton, T., Rombaut, P., Herzog, F. & Leonhardt, H. Determination of local chromatin composition by CasID. Nucleus 7, 476–484 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Robinett, C. C. et al. In vivo localization of DNA sequences and visualization of large-scale chromatin organization using lac operator/repressor recognition. J. Cell Biol. 135, 1685–1700 (1996).

    Article  CAS  PubMed  Google Scholar 

  31. Wang, T. et al. Crystal structure of the human SUV39H1 chromodomain and its recognition of histone H3K9me2/3. PLoS ONE 7, e52977 (2012).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  32. Zhang, K., Mosch, K., Fischle, W. & Grewal, S. I. S. Roles of the Clr4 methyltransferase complex in nucleation, spreading and maintenance of heterochromatin. Nat. Struct. Mol. Biol. 15, 381–388 (2008).

    Article  CAS  PubMed  Google Scholar 

  33. McCarthy, R. L., Zhang, J. & Zaret, K. S. Diverse heterochromatin states restricting cell identity and reprogramming. Trends Biochem. Sci 48, 513–526 (2023).

    Article  CAS  PubMed  Google Scholar 

  34. Grewal, S. I. S. The molecular basis of heterochromatin assembly and epigenetic inheritance. Mol. Cell 83, 1767–1785 (2023).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Wang, C., Zhu, B. & Xiong, J. Recruitment and reinforcement: maintaining epigenetic silencing. Sci. China Life Sci. 61, 515–522 (2018).

    Article  ADS  PubMed  Google Scholar 

  36. Ragunathan, K., Jih, G. & Moazed, D. Epigenetic inheritance uncoupled from sequence-specific recruitment. Science 348, 1258699 (2015).

    Article  PubMed  Google Scholar 

  37. Kudithipudi, S., Schuhmacher, M. K., Kebede, A. F. & Jeltsch, A. The SUV39H1 protein lysine methyltransferase methylates chromatin proteins involved in heterochromatin formation and VDJ recombination. ACS Chem. Biol. 12, 958–968 (2017).

    Article  CAS  PubMed  Google Scholar 

  38. Schuhmacher, M. K., Kudithipudi, S., Kusevic, D., Weirich, S. & Jeltsch, A. Activity and specificity of the human SUV39H2 protein lysine methyltransferase. Biochim. Biophys. Acta 1849, 55–63 (2015).

    Article  CAS  PubMed  Google Scholar 

  39. Burton, A. et al. Heterochromatin establishment during early mammalian development is regulated by pericentromeric RNA and characterized by non-repressive H3K9me3. Nat. Cell Biol. 22, 767–778 (2020).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Yang, P., Wang, Y. X. & Macfarlan, T. S. The role of KRAB-ZFPs in transposable element repression and mammalian evolution. Trends Genet. 33, 871–881 (2017).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Persikov, A. V. & Singh, M. De novo prediction of DNA-binding specificities for Cys2His2 zinc finger proteins. Nucleic Acids Res. 42, 97–108 (2014).

    Article  CAS  PubMed  Google Scholar 

  42. Wolfe, S. A., Nekludova, L. & Pabo, C. O. DNA recognition by Cys2His2 zinc finger proteins. Annu. Rev. Biophys. Biomol. Struct. 29, 183–212 (2000).

    Article  CAS  PubMed  Google Scholar 

  43. Bulut-Karslioglu, A. et al. Suv39h-dependent H3K9me3 marks intact retrotransposons and silences LINE elements in mouse embryonic stem cells. Mol. Cell 55, 277–290 (2014).

    Article  CAS  PubMed  Google Scholar 

  44. Liu, N. et al. Selective silencing of euchromatic L1s revealed by genome-wide screens for L1 regulators. Nature 553, 228–232 (2018).

    Article  ADS  CAS  PubMed  Google Scholar 

  45. Seczynska, M. & Lehner, P. J. The sound of silence: mechanisms and implications of HUSH complex function. Trends Genet. 39, 251–267 (2023).

    Article  CAS  PubMed  Google Scholar 

  46. Schotta, G. et al. A silencing pathway to induce H3-K9 and H4-K20 trimethylation at constitutive heterochromatin. Gene Dev. 18, 1251–1262 (2004).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Kourmouli, N. et al. Heterochromatin and tri-methylated lysine 20 of histone H4 in animals. J. Cell Sci. 117, 2491–2501 (2004).

    Article  CAS  PubMed  Google Scholar 

  48. Lehnertz, B. et al. Suv39h-mediated histone H3 lysine 9 methylation directs DNA methylation to major satellite repeats at pericentric heterochromatin. Curr. Biol. 13, 1192–1200 (2003).

    Article  CAS  PubMed  Google Scholar 

  49. Maison, C. et al. Higher-order structure in pericentric heterochromatin involves a distinct pattern of histone modification and an RNA component. Nat. Genet. 30, 329–334 (2002).

    Article  PubMed  Google Scholar 

  50. Booth, L. N. & Brunet, A. The aging epigenome. Mol. Cell 62, 728–744 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Tachibana, M. et al. G9a histone methyltransferase plays a dominant role in euchromatic histone H3 lysine 9 methylation and is essential for early embryogenesis. Genes Dev. 16, 1779–1791 (2002).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Xiong, J. et al. Cooperative action between SALL4A and TET proteins in stepwise oxidation of 5-methylcytosine. Mol. Cell 64, 913–925 (2016).

    Article  CAS  PubMed  Google Scholar 

  53. Xu, M., Chen, S. & Zhu, B. Investigating the cell cycle-associated dynamics of histone modifications using quantitative mass spectrometry. Methods Enzymol. 512, 29–55 (2012).

    Article  CAS  PubMed  Google Scholar 

  54. Hung, V. et al. Spatially resolved proteomic mapping in living cells with the engineered peroxidase APEX2. Nat. Protoc. 11, 456–475 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Kim, D., Paggi, J. M., Park, C., Bennett, C. & Salzberg, S. L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915 (2019).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Liao, Y., Smyth, G. K. & Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930 (2014).

    Article  CAS  PubMed  Google Scholar 

  58. Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).

    Article  PubMed  PubMed Central  Google Scholar 

  59. Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).

    Article  PubMed  PubMed Central  Google Scholar 

  61. Ramírez, F. et al. deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res. 44, W160–W165 (2016).

    Article  PubMed  PubMed Central  Google Scholar 

  62. Zhang, Y. et al. Model-based analysis of ChIP-seq (MACS). Genome Biol. 9, R137 (2008).

    Article  PubMed  PubMed Central  Google Scholar 

  63. Amemiya, H. M., Kundaje, A. & Boyle, A. P. The ENCODE blacklist: identification of problematic regions of the genome. Sci. Rep. 9, 9354 (2019).

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  64. Heinz, S. et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell 38, 576–589 (2010).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Ross-Innes, C. S. et al. Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature 481, 389–393 (2012).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  66. Tamura, K., Stecher, G. & Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis version 11. Mol. Biol. Evol. 38, 3022–3027 (2021).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank G. Li and L.-L. Du for helpful discussion; F. Suo, Q. Zhang and H. Zhang for assistance with protein structure prediction using AlphaFold; J. Jia and S. Meng for assistance with FACS experiments; and J. Wang and M. Zhang from the Laboratory of Proteomics of the Institute of Biophysics for assistance with mass spectrometry. This work was supported by grants from the China Natural Science Foundation (32288102, 32130021 and 92153302), the National Key R&D Programme of China (2018YFE0203302), the Chinese Academy of Sciences (XDB39010100 and JZHKYPT-2021-05), and the New Cornerstone Science Laboratory. Z.Z. is supported by the Youth Innovation Promotion Association (2017133) of the Chinese Academy of Sciences.

Author information

Authors and Affiliations

Authors

Contributions

B.Z. designed and supervised the project. R.M. designed and performed most of the experiments. Y.Z. performed most of the bioinformatics analysis. J.Z. participated in performing experiments during revision. Y.F. performed experiments for Extended Data Fig. 3c. H.-T.W. helped to prepare Extended Data Figs. 6b7 and 9b. P.Z., Z.L. and Z.Z. assisted with data analysis. R.M. and B.Z. wrote the manuscript, and Y.Z. participated in writing the methods and figure legends. All the authors read and commented on the manuscript.

Corresponding author

Correspondence to Bing Zhu.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature thanks the anonymous reviewer(s) for their contribution to the peer review of this work.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 The spatiotemporal expression pattern and binding sites of pericentric sequence-specific DNA binding proteins.

(a) The expression atlas of Pax3, Pax9, Foxd3, Znf512, and Znf512b in 51 experiments. TPM: transcript per million. (b) upper: major satellite consensus sequence in mouse pericentric regions and three types of satellite repeats in human pericentric regions. The binding sites of PAX3 and PAX9 in major satellite consensus sequence were highlighted in orange, and the binding site of FOXD3 in major satellite consensus sequence was highlighted in blue. TTC triples in mouse major satellite consensus sequence and human satellite repeats were highlighted in bold. lower: binding motifs of PAX3, PAX9 and FOXD3.

Extended Data Fig. 2 Workflow and the overlap analysis to characterize the protein composition of pericentric heterochromatin.

(a) Schematic diagram of the approach to map specific genomic site-associated proteome based on proximity labeling. (b) Overlap between proteins identified at pericentric heterochromatin by dCas9-APEX2 and candidates identified by PICh. A full list of pericentric heterochromatin associated proteins was included in Source Data for Extended Data Fig. 2.

Source Data

Extended Data Fig. 3 ZNF512 and ZNF512B can interact with SUV39H1 and SUV39H2 in vivo and in vitro.

(a) upper: A schematic diagram showing the LacO-LacI targeting system for testing whether a factor can establish H3K9me3 modified heterochromatin de novo. lower: A schematic diagram showing the LacO-LacI targeting system for detecting protein-protein interaction. (b) Validation of the LacO-LacI reporter system for testing whether a factor can establish H3K9me3 modified heterochromatin de novo. A representative cross-section immunofluorescence image showing the relative localization of LacI-mCherry fused SUV39H1 or SUV39H2, H3K9me3, and DAPI in CHO A03-1 cells. Condensed chromatin was stained by DAPI. Scale bar, 5 μm. (c) Co-IP and western blot analyzing the interaction between ZNF512 or ZNF512B and SUV39H1 or SUV39H2 in vivo. The mobilities of protein markers were indicated on the left of each panel. The asterisk refers to the degradation band of proteins. (d) MBP-Flag-ZNF512 or MBP-Flag-ZNF512B and MBP-HA-SUV39H1 or MBP-HA-SUV39H2 or MBP-HA-GFP-APEX2 were pulled down by HA beads. Bound proteins were eluted by HA peptide, separated on SDS-PAGE gels, and stained with Coomassie Brilliant Blue. The mobilities of protein markers were indicated on the right of each panel. For panels b, c, and d, representative data from two independent biological repeats are shown.

Extended Data Fig. 4 ZNF512 and ZNF512B do not appear to be essential for maintaining the molecular signature of pericentric heterochromatin.

(a) Transcriptome browser tracks representing the knockout of Znf512, Znf512b, Suv39h1 and Suv39h2 in wild type cell, Znf512 single KO cell, Znf512b single KO cell, Znf512 and Znf512b double knockout cell, and QKO cell. (b) RT-qPCR showing the relative expression level of major satellite RNA in wild-type and different KO cells. The expression level was normalized to Gapdh. Bars indicate the mean; dots indicate individual replicate. (c) A representative cross-section image showing the localization of H3K9me3 in wild type, Znf512 single KO, Znf512b single KO, Znf512 and Znf512b double KO cells. Scale bar, 2 μm. For panel b, data are from two independent biological replicates. For panel c, representative data from four independent biological repeats are shown.

Source Data

Extended Data Fig. 5 The localization of H3K9me3 and mutated EGFP-SUV39H in QKO cells and the localization of endogenous SUV39H2 in H1 TKO cells.

(a) A representative cross-section image showing the localization of H3K9me3 in Znf512, Znf512b, Suv39h1 and Suv39h2 quadruple KO (QKO) cells. Scale bar, 5 μm. (b) The localization of EGFP-SUV39H1, EGFP-SUV39H2 wild type or mutated versions and H3K9me3 were analyzed in Znf512, Znf512b, Suv39h1, Suv39h2 quadruple KO (QKO) cell line. SUV39H1-W64A/Y67A, the aromatic cage mutant, was abbreviated as SUV39H1-2A, and SUV39H2-W139A/W142 A, was abbreviated as SUV39H2-2A. SUV39H1-H324K and SUV39H2-H398K were catalytically dead mutants. Scale bar, 5 μm. (c) A representative cross-section image showing the localization of endogenous SUV39H2 and DAPI in wild-type and different KO cells. Scale bar, 5 μm. For panel a, b, and c, representative data from two independent biological repeats are shown.

Extended Data Fig. 6 ZNF512 and ZNF512B utilize their zinc fingers to target the pericentric sequences.

(a) Schematic diagram of ZNF512 and ZNF512B protein sequences. Zinc finger domains and linkers’ lengths were labeled. (b) Multi-sequence alignments of zinc finger domains in mouse ZNF512 and ZNF512B. The zinc fingerprint residues at positions −1, 3, 6 of each C2H2 zinc finger were highlighted. (c) The binding motifs of mouse ZNF512 and ZNF512B were predicted by a polynomial SVM-based algorithm. (d) Representative cross-section fluorescence images showing the relative localization of EGFP-ZNF512 with one or two zinc fingerprints mutated and DAPI in Znf512, Znf512b, Suv39h1, Suv39h2 quadruple KO (QKO) cells. Zinc fingerprint residues mutated versions, in which one zinc finger of ZNF512 mutated to AAA, were abbreviated as ZF-3A. Zinc fingerprint residues mutated versions, which two zinc fingers of ZNF512 mutated to AAA, were abbreviated as ZF-6A. Scale bar, 5 μm. (e) Representative cross-section fluorescence images showing the relative localization of EGFP-ZNF512B with one or two zinc fingerprints mutated and DAPI in Znf512, Znf512b, Suv39h1, Suv39h2 quadruple KO (QKO) cells. Zinc fingerprint residues mutated versions, which one zinc finger of ZNF512B mutated to AAA, were abbreviated as ZF-3A. Zinc fingerprint residues mutated versions, which two zinc fingers of ZNF512B mutated to AAA, were abbreviated as ZF-6A. Scale bar, 5 μm. Data are representative of three (d) or one (e) independent biological experiments.

Extended Data Fig. 7 Evolution of the zinc finger domains and the length of linkers in ZNF512 and ZNF512B.

Phylogenetic tree showing the evolutionary relationship of zinc finger domains of ZNF512 and ZNF512B in humans, mice, chicken, mainland tiger snake, western clawed frog, zebrafish, sea lamprey, and Florida lancelet. Zinc fingerprints were highlighted in red, and predicted binding motifs were highlighted in blue. The length of linkers in ZNF512 and ZNF512B was also evolutionarily conserved among vertebrates.

Extended Data Fig. 8 The long linkers play an important role in ZNF512 and ZNF512B for targeting the pericentric sequences.

(a) Bar plots showing the ratio of the percentage of reads mapped on mouse major satellite repeats between biotin ChIPed samples and the corresponding input samples in indicated cell lines. (b) Bar plot showing the fold enrichment of TTC triplets in the consensus sequence of mouse major satellite repeats. (c) Box plots showing the fold enrichment of TTC triplets in human pericentric satellite III/II/I. Centre line, median; box, 25th and 75th percentiles; whiskers, 1.5 × the interquartile range (IQR). Numbers of satellites: n = 189 (satellite III), n = 61 (satellite II), n = 88 (satellite I). (d) Bar plots showing the ratio of the percentage of reads mapped on human pericentric satellite repeats between ZNF512 or ZNF512B ChIPed samples and the corresponding control samples in indicated cell lines. Bars indicate the mean; dots indicate individual replicate. (e, f) Predicted structures of wild-type or linkers-swapped ZNF512 or ZNF512B proteins using AlphaFold. Zinc finger domains were highlighted in red color. The fourth zinc finger domain of ZNF512B was colored green. (g) Representative cross-section fluorescence images showing the relative localization of EGFP-ZNF512 or EGFP-ZNF512B linker-swapped versions in QKO cells. Scale bar, 5 μm. (h) After overexpressing Flag-tagged APEX2, ZNF512 or ZNF512B wild type or linker-swapped version in Znf512,Znf512b,Suv39h2 triple KO (H2 TKO) cells, H3K9me3 modification at pericentric regions was visualized. Scale bar, 5 μm. (i) After overexpressing mCherry-tagged APEX2, ZNF512 or ZNF512B wild type or linker-swapped version in H2 TKO cells, endogenous SUV39H1 at pericentric regions was visualized. Scale bar, 5 μm. (j) DNA pull-down with biotin-major satellite probe or no probe followed by western blot analysis for MBP tag using purified recombinant MBP-Flag-tagged ZNF512 or ZNF512B. The mobilities of protein markers were indicated on the left of the panel. Data are representative of two (h, i and j) or one (g) independent biological experiments.

Source Data

Extended Data Fig. 9 H3K9me3 ChIP-seq signals around LINEs.

(a) Profile (top) and heatmap (bottom) showing H3K9me3 signals in WT and Suv39h1 and Suv39h2 double knockout mESCs around LINEs (n = 6,996). The difference in H3K9me3 between Suv39h1 and Suv39h2 double knockout and WT mESCs was also shown. The x-axis represented the distance from the start or end of LINE in kilobases (Kb). Heatmaps were sorted by H3K9me3 signal. The average signals of two replicates were shown. (b) Percentage of ZNF512-bound full-length LINEs per subfamily in mouse ESCs arranged from the youngest to the oldest subfamily. LINE ages obtained from previously published studies. (Myr) Million years.

Extended Data Fig. 10 Reproducibility of replicates in this study.

(a) Overall similarity of ZNF512 signals within ZNF512 peaks in WT mouse ES cells (n = 3,218). Pearson correlation coefficients were indicated. (b) Overall similarity of H3K9me3 signals within H3K9me3 peaks in WT mouse ES cells (n = 29,982). Pearson correlation coefficients were indicated.

Supplementary information

Supplementary Information

Supplementary Figs. 1 and 2.

Reporting Summary

Source data

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ma, R., Zhang, Y., Zhang, J. et al. Targeting pericentric non-consecutive motifs for heterochromatin initiation. Nature 631, 678–685 (2024). https://doi.org/10.1038/s41586-024-07640-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41586-024-07640-5

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing