Allelic reprogramming of the histone modification H3K4me3 in early mammalian development (original) (raw)

Accession codes

Primary accessions

Gene Expression Omnibus

Data deposits

All data have been deposited to GEO with the accession number GSE71434.

References

  1. Kouzarides, T. Chromatin modifications and their function. Cell 128, 693–705 (2007)
    Article CAS PubMed Google Scholar
  2. Ruthenburg, A. J., Allis, C. D. & Wysocka, J. Methylation of lysine 4 on histone H3: intricacy of writing and reading a single epigenetic mark. Mol. Cell 25, 15–30 (2007)
    Article CAS PubMed Google Scholar
  3. Peng, X. et al. TELP, a sensitive and versatile library construction method for next-generation sequencing. Nucleic Acids Res . 43, e35 (2015)
    Article PubMed Google Scholar
  4. Hamatani, T., Carter, M. G., Sharov, A. A. & Ko, M. S. Dynamics of global gene expression changes during mouse preimplantation development. Dev. Cell 6, 117–131 (2004)
    Article CAS PubMed Google Scholar
  5. Picelli, S. et al. Full-length RNA-seq from single cells using Smart-seq2. Nat. Protocols 9, 171–181 (2014)
    Article CAS PubMed Google Scholar
  6. Lee, M. T., Bonneau, A. R. & Giraldez, A. J. Zygotic genome activation during the maternal-to-zygotic transition. Annu. Rev. Cell Dev. Biol. 30, 581–613 (2014)
    Article CAS PubMed PubMed Central Google Scholar
  7. Peaston, A. E. et al. Retrotransposons regulate host genes in mouse oocytes and preimplantation embryos. Dev. Cell 7, 597–606 (2004)
    Article CAS PubMed Google Scholar
  8. Smith, Z. D. & Meissner, A. DNA methylation: roles in mammalian development. Nat. Rev. Genet. 14, 204–220 (2013)
    Article CAS PubMed Google Scholar
  9. Wu, H. & Zhang, Y. Reversing DNA methylation: mechanisms, genomics, and biological functions. Cell 156, 45–68 (2014)
    Article CAS PubMed PubMed Central Google Scholar
  10. Erkek, S. et al. Molecular determinants of nucleosome retention at CpG-rich sequences in mouse spermatozoa. Nat. Struct. Mol. Biol. 20, 868–875 (2013)
    Article CAS PubMed Google Scholar
  11. Xie, W. et al. Base-resolution analyses of sequence and parent-of-origin dependent DNA methylation in the mouse genome. Cell 148, 816–831 (2012)
    Article CAS PubMed PubMed Central Google Scholar
  12. Lepikhov, K. & Walter, J. Differential dynamics of histone H3 methylation at positions K4 and K9 in the mouse zygote. BMC Dev. Biol. 4, 12 (2004)
    Article PubMed PubMed Central Google Scholar
  13. Bernstein, B. E. et al. A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell 125, 315–326 (2006)
    Article CAS PubMed Google Scholar
  14. Wang, L. et al. Programming and inheritance of parental DNA methylomes in mammals. Cell 157, 979–991 (2014)
    Article CAS PubMed PubMed Central Google Scholar
  15. Tomizawa, S., Nowacka-Woszuk, J. & Kelsey, G. DNA methylation establishment during oocyte growth: mechanisms and significance. Int. J. Dev. Biol. 56, 867–875 (2012)
    Article CAS PubMed Google Scholar
  16. Lister, R. et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature 462, 315–322 (2009)
    Article ADS CAS PubMed PubMed Central Google Scholar
  17. Kageyama, S. et al. Alterations in epigenetic modifications during oocyte growth in mice. Reproduction 133, 85–94 (2007)
    Article CAS PubMed Google Scholar
  18. Andreu-Vieyra, C. V. et al. MLL2 is required in oocytes for bulk histone 3 lysine 4 trimethylation and transcriptional silencing. PLoS Biol . 8, e1000453 (2010)
    Article PubMed PubMed Central Google Scholar
  19. Shen, Y. et al. A map of the cis-regulatory sequences in the mouse genome. Nature 488, 116–120 (2012)
    Article ADS CAS PubMed PubMed Central Google Scholar
  20. Hon, G. C. et al. Epigenetic memory at embryonic enhancers identified in DNA methylation maps from adult mouse tissues. Nat. Genet. 45, 1198–1206 (2013)
    Article CAS PubMed PubMed Central Google Scholar
  21. Solter, D. & Knowles, B. B. Immunosurgery of mouse blastocyst. Proc. Natl Acad. Sci. USA 72, 5099–5102 (1975)
    Article ADS CAS PubMed PubMed Central Google Scholar
  22. Boroviak, T., Loos, R., Bertone, P., Smith, A. & Nichols, J. The ability of inner-cell-mass cells to self-renew as embryonic stem cells is acquired following epiblast specification. Nat. Cell Biol. 16, 516–528 (2014)
    Article CAS PubMed PubMed Central Google Scholar
  23. Demeestere, I. et al. Effect of preantral follicle isolation technique on in-vitro follicular growth, oocyte maturation and embryo development in mice. Hum. Reprod. 17, 2152–2159 (2002)
    Article CAS PubMed Google Scholar
  24. Su, Y. Q. et al. MARF1 regulates essential oogenic processes in mice. Science 335, 1496–1499 (2012)
    Article ADS CAS PubMed PubMed Central Google Scholar
  25. Sharif, B. et al. The chromosome passenger complex is required for fidelity of chromosome transmission and cytokinesis in meiosis of mouse oocytes. J. Cell Sci. 123, 4292–4300 (2010)
    Article CAS PubMed PubMed Central Google Scholar
  26. Bouniol-Baly, C. et al. Differential transcriptional activity associated with chromatin configuration in fully grown mouse germinal vesicle oocytes. Biol. Reprod. 60, 580–587 (1999)
    Article CAS PubMed Google Scholar
  27. Whitlock, K. E. & Westerfield, M. The olfactory placodes of the zebrafish form by convergence of cellular fields at the edge of the neural plate. Development 127, 3645–3653 (2000)
    Article CAS PubMed Google Scholar
  28. Gilfillan, G. D. et al. Limitations and possibilities of low cell number ChIP–seq. BMC Genomics 13, 645 (2012)
    Article CAS PubMed PubMed Central Google Scholar
  29. Hisano, M. et al. Genome-wide chromatin analysis in mature mouse and human spermatozoa. Nat. Protocols 8, 2449–2470 (2013)
    Article CAS PubMed Google Scholar
  30. Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protocols 7, 562–578 (2012)
    Article CAS PubMed Google Scholar
  31. Jurka, J. et al. Repbase Update, a database of eukaryotic repetitive elements. Cytogenet. Genome Res. 110, 462–467 (2005)
    Article CAS PubMed Google Scholar
  32. 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
  33. Zhang, Y. et al. Model-based analysis of ChIP–seq (MACS). Genome Biol . 9, R137 (2008)
    Article PubMed PubMed Central Google Scholar
  34. Feng, J., Liu, T., Qin, B., Zhang, Y. & Liu, X. S. Identifying ChIP–seq enrichment using MACS. Nat. Protocols 7, 1728–1740 (2012)
    Article CAS PubMed Google Scholar
  35. Yue, F. et al. A comparative encyclopedia of DNA elements in the mouse genome. Nature 515, 355–364 (2014)
    Article CAS PubMed PubMed Central Google Scholar
  36. Deng, Q., Ramsköld, D., Reinius, B. & Sandberg, R. Single-cell RNA-seq reveals dynamic, random monoallelic gene expression in mammalian cells. Science 343, 193–196 (2014)
    Article ADS CAS PubMed Google Scholar
  37. Kobayashi, H. et al. High-resolution DNA methylome analysis of primordial germ cells identifies gender-specific reprogramming in mice. Genome Res . 23, 616–627 (2013)
    Article CAS PubMed PubMed Central Google Scholar
  38. Zhang, Y. et al. Canonical nucleosome organization at promoters forms during genome activation. Genome Res . 24, 260–266 (2014)
    Article PubMed PubMed Central Google Scholar
  39. Bogdanovic, O. et al. Dynamics of enhancer chromatin signatures mark the transition from pluripotency to cell specification during embryogenesis. Genome Res . 22, 2043–2053 (2012)
    Article CAS PubMed PubMed Central Google Scholar
  40. Potok, M. E., Nix, D. A., Parnell, T. J. & Cairns, B. R. Reprogramming the maternal zebrafish genome after fertilization to match the paternal methylation pattern. Cell 153, 759–772 (2013)
    Article CAS PubMed PubMed Central Google Scholar

Download references

Acknowledgements

We appreciate help from members of the Xie laboratory for comments during preparation of the manuscript. We thank J. Zhao and G. Qin for helping with pronuclei isolation. We are grateful to the animal core facility, the sequencing core facility, and biocomputing facility at Tsinghua University. This work is supported by the funding provided by the National Basic Research Program of China 2016YFC0900301 (W.X.), 2015CB856201 (W.X.), 2012CB966701 (J.N.), the National Natural Science Foundation of China 31422031 (W.X.), 91519326 (W.X.), 31171381 (J.N.), the Beijing Natural Science Foundation grant 5152014 (J.N.), Tsinghua University Initiative Scientific Research Program 20161080043 (W.X.), the funding from the THU-PKU Center for Life Sciences (W.X., J.N.), and the Youth Thousand Scholar Program of China (W.X.).

Author information

Author notes

  1. Bingjie Zhang, Hui Zheng, Bo Huang and Wenzhi Li: These authors contributed equally to this work.

Authors and Affiliations

  1. Center for Stem Cell Biology and Regenerative Medicine, MOE Key Laboratory of Bioinformatics, THU-PKU Center for Life Sciences, School of Life Sciences, Tsinghua University, Beijing, 100084, China
    Bingjie Zhang, Hui Zheng, Bo Huang, Yunlong Xiang, Yu Zhang, Qianhua Xu, Yuanyuan Li, Qiujun Wang, Jing Ma, Qiangzong Yin & Wei Xie
  2. PKU-THU Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, 100871, China
    Bo Huang
  3. Center for Stem Cell Biology and Regenerative Medicine, School of Medicine, Tsinghua University, Beijing, 100084, China
    Wenzhi Li, Jia Ming, Bo Chen, Kehkooi Kee & Jie Na
  4. Singapore Institute for Clinical Sciences, Agency for Science, Technology and Research (A*STAR), Singapore, 117609, Singapore
    Xu Peng & Feng Xu
  5. State Key Laboratory of Biomembrane and Membrane Engineering, School of Life Sciences, Tsinghua University, Beijing, 100084, China
    Xiaotong Wu & Anming Meng
  6. Clinical and Translational Research Center of Shanghai First Maternity and Infant Hospital, School of Life Sciences and Technology, Tongji University, Shanghai, 200092, China
    Wenqiang Liu, Xiaochen Kou, Yanhong Zhao, Wenteng He, Chong Li & Shaorong Gao
  7. Institute of Molecular and Cell Biology, A*STAR, Singapore, 138673, Singapore
    Feng Xu

Authors

  1. Bingjie Zhang
    You can also search for this author inPubMed Google Scholar
  2. Hui Zheng
    You can also search for this author inPubMed Google Scholar
  3. Bo Huang
    You can also search for this author inPubMed Google Scholar
  4. Wenzhi Li
    You can also search for this author inPubMed Google Scholar
  5. Yunlong Xiang
    You can also search for this author inPubMed Google Scholar
  6. Xu Peng
    You can also search for this author inPubMed Google Scholar
  7. Jia Ming
    You can also search for this author inPubMed Google Scholar
  8. Xiaotong Wu
    You can also search for this author inPubMed Google Scholar
  9. Yu Zhang
    You can also search for this author inPubMed Google Scholar
  10. Qianhua Xu
    You can also search for this author inPubMed Google Scholar
  11. Wenqiang Liu
    You can also search for this author inPubMed Google Scholar
  12. Xiaochen Kou
    You can also search for this author inPubMed Google Scholar
  13. Yanhong Zhao
    You can also search for this author inPubMed Google Scholar
  14. Wenteng He
    You can also search for this author inPubMed Google Scholar
  15. Chong Li
    You can also search for this author inPubMed Google Scholar
  16. Bo Chen
    You can also search for this author inPubMed Google Scholar
  17. Yuanyuan Li
    You can also search for this author inPubMed Google Scholar
  18. Qiujun Wang
    You can also search for this author inPubMed Google Scholar
  19. Jing Ma
    You can also search for this author inPubMed Google Scholar
  20. Qiangzong Yin
    You can also search for this author inPubMed Google Scholar
  21. Kehkooi Kee
    You can also search for this author inPubMed Google Scholar
  22. Anming Meng
    You can also search for this author inPubMed Google Scholar
  23. Shaorong Gao
    You can also search for this author inPubMed Google Scholar
  24. Feng Xu
    You can also search for this author inPubMed Google Scholar
  25. Jie Na
    You can also search for this author inPubMed Google Scholar
  26. Wei Xie
    You can also search for this author inPubMed Google Scholar

Contributions

B.Z. developed STAR ChIP–seq and conducted ChIP–seq experiments; B.H., Y.X. and J. Ma collected mouse oocytes and early embryos. B.H. performed RNA-seq; H.Z. performed data analysis; W.L., J. Ming and B.Z. conducted KDM5B experiment. W.L., X.K., Y. Zhao, W.H., C.L. and B.H. collected mouse pronuclei. X.P., B.H. and F.X. advised the development of STAR ChIP–seq. Y. Zhang and Q.Y. performed oocyte DNA methylome profiling. Y.L. and Q.W. conducted high-throughput sequencing; X.W. collected zebrafish oocytes; K.K., A.M., S.G., J.N. and W.X. supervised the project or related experiments. H.Z. and W.X. wrote the manuscript.

Corresponding authors

Correspondence toJie Na or Wei Xie.

Ethics declarations

Competing interests

B.Z. and W.X. are co-inventors on a filed patent application for STAR ChIP–seq.

Additional information

Reviewer Information Nature thanks R. Schultz and the other anonymous reviewer(s) for their contribution to the peer review of this work.

Extended data figures and tables

Extended Data Figure 1 Validation of STAR ChIP–seq in mES cells.

a, A brief schematic of STAR ChIP–seq. Based on a previously described ChIP method28, micrococcal nuclease (MNase) is used to fragment native chromatin. Nucleosomes with certain histone modifications are immunoprecipitated using specific antibodies. Proteinase K is added for protein digestion and is heat-inactivated. Shrimp alkaline phosphatase (rSAP) is directly added to the reaction for 3′ DNA end repair, which is required for the poly-C addition in the subsequent TELP library preparation. After heat-inactivation of rSAP, the resulting DNA is subjected to the TELP library preparation as previously described3 without purification. Briefly, poly-C tailing is conducted on denatured single-strand DNA using dCTP and terminal deoxynucleotidyl transferase (TDT). Biotin-labelled anchor primer (indicated by ‘B’) containing poly-G is used for second strand DNA extension. The products are captured by magnetic streptavidin beads. After an adaptor ligation to the opposite end of poly-C, the double-stranded DNA is amplified by primers (P1 and P2) containing Illumina adaptor sequences. The resulting library is ready to be sequenced. b, A snapshot of the UCSC browser view shows enrichment of STAR ChIP–seq for histone modifications H3K4me3, H3K27ac and H3K27me3 using various numbers of mES cells. Data sets from mouse ENCODE are also shown for comparison35. c, The Pearson correlation coefficients showing the comparison between STAR ChIP–seq of H3K4me3 using various numbers of mES cells and conventional ChIP–seq of H3K4me3 from either this study using 1 × 106 mES cells or data from ENCODE35. d, A bar chart showing the percentages of conventional mESC H3K4me3 total/promoter peaks recaptured by STAR ChIP using various numbers of mES cells. e, Box plots showing the H3K4me3 enrichment determined by conventional ChIP–seq at peaks that are either recaptured or missed by STAR ChIP–seq.

Extended Data Figure 2 Validation of H3K4me3 STAR ChIP–seq and RNA-seq in mouse gametes and early embryos.

a, Representative bright-field microscopic images showing the pre-implantation embryos used for STAR ChIP–seq with time points of sample collection indicated. b, The Spearman correlation between the replicates (n = 2) of RNA-seq samples and between RNA-seq in this study and Deng et al.36 at stages when available. c, The Pearson correlation between the replicates (n = 2) of H3K4me3 for mouse pre-implantation embryos at each stage. d, Scatter plots showing the H3K4me3 enrichment levels (2-kb window) between replicates across the genome for mouse pre-implantation embryos at each stage. e, Scatter plot showing the H3K4me3 enrichment levels (2-kb window) in MII oocytes between data sets generated using an in-house antibody or a commercial antibody (Millipore 04-745). f, A snapshot of the UCSC browser view shows H3K4me3 STAR ChIP–seq results using two different H3K4me3 antibodies (in-house or commercial/Millipore 04-745) in mES cells and MII oocytes. Annotated genes and repetitive elements are also shown. g, Venn diagram showing the overlap between mES cell H3K4me3 (conventional ChIP–seq) promoter peaks and ICM H3K4me3 promoter peaks. h, Snapshots of the UCSC browser views showing the H3K4me3 and expression signals near Foxa1 (left), Atp1b1 (middle) and Tbx3 (right).

Extended Data Figure 3 Characterization of promoter and distal H3K4me3 in mouse gametes and early embryos.

a, Scatter plots showing the relationship between the expression levels of all genes and their promoter H3K4me3 enrichment (_Z_-score-normalized) for each development stage. Mouse ES cell data from ENCODE are used as control. Spearman correlation coefficients are also shown. b, Moving average plots (window = 100 genes, moving step = 1 gene) showing the relationship between the expression levels of non-maternal genes and their promoter H3K4me3 enrichment (_Z_-score-normalized) for each stage after major ZGA. c, The genomic distributions of H3K4me3 peaks are shown for early embryos and mES cells. TSS, transcription start site (±2.5 kb); TES, transcription end site (±2.5 kb). d, A bar chart showing the percentages of all mappable reads that can be aligned to RepBase for mouse embryos of various stages. The error bars denote the standard deviations of two biological replicates of RNA-seq for each stage. e, Enrichment of repeats (log ratio of observed/random) in H3K4me3 promoter and distal peaks compared to that in random peaks at each stage. f, The enrichment (log ratio of observed/random) of top repeat subfamily members in distal H3K4me3 peaks in mouse embryos of various stages and mES cells. g, The snapshots of distal H3K4me3 signals near various types of repeats.

Extended Data Figure 4 Reprogramming of H3K4me3 on the paternal genome.

a, A snapshot of UCSC genome browser shows the enrichment of H3K4me3 in sperm from a published data set10 and this study. b, Snapshots showing the allelic H3K4me3 signals and RNA-seq signals near the imprinted loci Snrpn and Impact. Please note the maternal RNAs in zygotes, early two-cell and late two-cell embryos are presumably inherited from oocytes. The paternal zygotic transcripts appear from the late two-cell stage onward for both genes. c, Bar chart showing the numbers of H3K4me3 peaks identified in sperm (in regions covered by SNPs) and the paternal alleles of early embryos. d, Hierarchical clustering of sperm, early embryos and somatic cortex based on the paternal allele H3K4me3 enrichment. For sperm H3K4me3 data, only regions covered by SNPs were included for analysis. e, A bar chart showing the percentages of raw reads assigned to the maternal or paternal allele for each stage. f, Scatter plots showing the H3K4me3 enrichment levels between SNP-assigned parental alleles from zygotes and pronuclei for both the paternal (left) and maternal genomes (right).

Extended Data Figure 5 Non-canonical H3K4me3 (ncH3K4me3) in oocytes and early embryos

a, A snapshot of the UCSC browser view shows the enrichment of H3K4me3 in MII oocytes, isolated maternal pronuclei, PN5 zygotes (maternal allele), and histone H3 (maternal allele) in PN5 zygotes. b, Scatter plots showing the comparisons between the enrichment of histone H3 (maternal) in PN5 zygotes and H3K4me3 from either MII oocytes (left) or the maternal allele of PN5 zygotes (right). c, Hierarchical clustering of early embryo stages (maternal) based on global H3K4me3 enrichment. d, A bar chart showing the expression levels (log2[FPKM + 1]) of several marker genes for granulosa cells and oocytes, as well as a house-keeping gene (Actb), in oocyte samples collected at different stages. ICMs are included as a non-oocyte control. e, Heat maps showing the expression (FPKM) of all genes and the H3K4me3 enrichment (normalized RPKM) at their promoters (TSS ± 2.5 kb) for oocytes and early embryos. f, Hierarchical clustering of stages of oocytes and early embryos based on global H3K4me3 enrichment. g, A snapshot of UCSC browser shows the H3K4me3 enrichment in oocyte of various stages as well as DNA methylation in MII oocyte. A PMD with ncH3K4me3 enrichment is shaded.

Extended Data Figure 6 Non-canonical distal H3K4me3 in oocytes correlates with DNA hypomethylation.

a, Box plots showing the CpG density of H3K4me3 promoter and distal peaks in early embryos, mES cells, and the cortex. b, Pie charts showing the percentage of the genome covered by PMDs (left) and the percentage of distal H3K4me3 marked regions covered by PMDs (right) in MII oocytes. c, A global view of H3K4me3 promoter and distal peak density of chr12 in MII oocytes (top). Gene densities and PMD densities are also shown (bottom). Three gene desert regions are shaded. d, DNA methylation levels around active genes are shown for those in embryonic day 16.5 (E16.5) female PGCs37, day 14 oocytes, and MII oocytes14.

Extended Data Figure 7 The transition from ncH3K4me3 to canonical H3K4me3 requires zygotic transcription but not DNA replication.

a, Box plots of gene expression levels (log2[FPKM + 1]) show the effects of α-amanitin and aphidicolin on genes activated during ZGA (left). Genes expressed at all stages (middle) and genes repressed at all stages (right) are also shown. b, Hierarchical clustering of two alleles of early embryos including those treated with α-amanitin or aphidicolin based on the global H3K4me3 enrichment. c, A snapshot of the UCSC browser view shows the global view of H3K4me3 enrichment for maternal and paternal alleles in various types of embryos including those treated with α-amanitin or aphidicolin. d, Bar charts showing normalized H3K4me3 and BrUTP fluorescence intensity of individual surrounded nucleolus (SN)-stage oocytes. Each pair of bars for H3K4me3 and BrUTP represents an individual oocyte. Data from three independent experiments were pooled. GV, germinal vesicle.

Extended Data Figure 8 ncH3K4me3 is not found in zebrafish.

a, A snapshot of the UCSC genome browser shows the global view of H3K4me3 enrichment in oocytes (this study), the dome stage embryos38, 48 hours post fertilization (hpf) embryos39 and DNA methylation level40 for various zebrafish developmental stages. b, The genomic distributions are shown for the H3K4me3 peaks in zebrafish oocytes, embryos at the dome stage38 and 48 hpf embryos39. c, Box plots showing the CpG densities of H3K4me3 promoter peaks in zebrafish oocytes, embryos at the dome stage38 and embryos at 48 hpf39. A random set of peaks that match the lengths of H3K4me3 peaks on the same chromosomes were used as a control in the analysis. d, Violin plots showing the distributions of DNA methylation levels from gametes to somatic tissues in zebrafish (left panel) or mice (right panel). e, Box plot showing the lengths of distal PMDs in zebrafish and mouse oocytes, sperms, and somatic tissues. f, Bar chart shows the percentages of the genome covered by distal PMDs in oocytes, sperms and somatic tissues in zebrafish or mouse.

Supplementary information

PowerPoint slides

Rights and permissions

About this article

Cite this article

Zhang, B., Zheng, H., Huang, B. et al. Allelic reprogramming of the histone modification H3K4me3 in early mammalian development.Nature 537, 553–557 (2016). https://doi.org/10.1038/nature19361

Download citation