Spatiotemporal clustering of the epigenome reveals rules of dynamic gene regulation

Spatial organization of different epigenomic marks was used to infer functions of the epigenome. It remains unclear what can be learned from the temporal changes of the epigenome. Here, we developed a probabilistic model to cluster genomic sequences based on the similarity of temporal changes of multiple epigenomic marks during a cellular differentiation process. We differentiated mouse embryonic stem (ES) cells into mesendoderm cells. At three time points during this differentiation process, we used high-throughput sequencing to measure seven histone modifications and variants--H3K4me1/2/3, H3K27ac, H3K27me3, H3K36me3, and H2A.Z; two DNA modifications--5-mC and 5-hmC; and transcribed mRNAs and noncoding RNAs (ncRNAs). Genomic sequences were clustered based on the spatiotemporal epigenomic information. These clusters not only clearly distinguished gene bodies, promoters, and enhancers, but also were predictive of bidirectional promoters, miRNA promoters, and piRNAs. This suggests specific epigenomic patterns exist on piRNA genes much earlier than germ cell development. Temporal changes of H3K4me2, unmethylated CpG, and H2A.Z were predictive of 5-hmC changes, suggesting unmethylated CpG and H3K4me2 as potential upstream signals guiding TETs to specific sequences. Several rules on combinatorial epigenomic changes and their effects on mRNA expression and ncRNA expression were derived, including a simple rule governing the relationship between 5-hmC and gene expression levels. A Sox17 enhancer containing a FOXA2 binding site and a Foxa2 enhancer containing a SOX17 binding site were identified, suggesting a positive feedback loop between the two mesendoderm transcription factors. These data illustrate the power of using epigenome dynamics to investigate regulatory functions.

[1]  W. Reik,et al.  Dynamic regulation of 5-hydroxymethylcytosine in mouse ES cells and during differentiation , 2011, Nature.

[2]  Megan F. Cole,et al.  Connecting microRNA Genes to the Core Transcriptional Regulatory Circuitry of Embryonic Stem Cells , 2008, Cell.

[3]  Tetsuya S. Tanaka,et al.  Short-Term Serum-Free Culture Reveals that Inhibition of Gsk3β Induces the Tumor-Like Growth of Mouse Embryonic Stem Cells , 2011, PloS one.

[4]  J. Lewerenz,et al.  Short-term serum supplementation improves glucose-oxygen deprivation in primary cortical cultures grown under serum-free conditions. , 2003, Methods in cell science : an official journal of the Society for In Vitro Biology.

[5]  Dan Xie,et al.  Comparative Epigenomic Annotation of Regulatory DNA , 2012, Cell.

[6]  A. Ribeiro,et al.  Genome wide study of NF-Y type CCAAT boxes in unidirectional and bidirectional promoters in human and mouse. , 2011, Journal of theoretical biology.

[7]  A. Perkins,et al.  Dynamic transcription programs during ES cell differentiation towards mesoderm in serum versus serum-freeBMP4 culture , 2007, BMC Genomics.

[8]  Manolis Kellis,et al.  Discovery and Characterization of Chromatin States for Systematic Annotation of the Human Genome , 2011, RECOMB.

[9]  P. Jin,et al.  Selective chemical labeling reveals the genome-wide distribution of 5-hydroxymethylcytosine , 2011, Nature Biotechnology.

[10]  Lee E. Edsall,et al.  Distinct epigenomic landscapes of pluripotent and lineage-committed human cells. , 2010, Cell stem cell.

[11]  Nicolas E. Buchler,et al.  On schemes of combinatorial transcription logic , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[12]  S. Schreiber,et al.  Histone Variant H2A.Z Marks the 5′ Ends of Both Active and Inactive Genes in Euchromatin , 2006, Cell.

[13]  E. Lander,et al.  The Mammalian Epigenome , 2007, Cell.

[14]  Saurabh Sinha,et al.  A Biophysical Model for Analysis of Transcription Factor Interaction and Binding Site Arrangement from Genome-Wide Binding Data , 2009, PloS one.

[15]  Ernest Fraenkel,et al.  A Quantitative Model of Transcriptional Regulation Reveals the Influence of Binding Location on Expression , 2010, PLoS Comput. Biol..

[16]  Yoshiakira Kanai,et al.  Depletion of definitive gut endoderm in Sox17-null mutant mice. , 2002, Development.

[17]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[18]  J. Ahringer,et al.  Differential chromatin marking of introns and expressed exons by H3K36me3 , 2008, Nature Genetics.

[19]  Keji Zhao,et al.  Genome-wide analysis of 5-hydroxymethylcytosine distribution reveals its dual function in transcriptional regulation in mouse embryonic stem cells. , 2011, Genes & development.

[20]  Ryan A. Flynn,et al.  A unique chromatin signature uncovers early developmental enhancers in humans , 2011, Nature.

[21]  Manolis Kellis,et al.  Dynamics of the epigenetic landscape during erythroid differentiation after GATA1 restoration. , 2011, Genome research.

[22]  G. Hon,et al.  Base-Resolution Analysis of 5-Hydroxymethylcytosine in the Mammalian Genome , 2012, Cell.

[23]  G. Ast,et al.  Chromatin organization marks exon-intron structure , 2009, Nature Structural &Molecular Biology.

[24]  Allen D. Delaney,et al.  Conserved Role of Intragenic DNA Methylation in Regulating Alternative Promoters , 2010, Nature.

[25]  R. Young,et al.  Histone H3K27ac separates active from poised enhancers and predicts developmental state , 2010, Proceedings of the National Academy of Sciences.

[26]  S. Andrews,et al.  Dynamic stage-specific changes in imprinted differentially methylated regions during early mammalian development and prevalence of non-CpG methylation in oocytes , 2011, Development.

[27]  B. Williams,et al.  Mapping and quantifying mammalian transcriptomes by RNA-Seq , 2008, Nature Methods.

[28]  David R. Liu,et al.  Conversion of 5-Methylcytosine to 5- Hydroxymethylcytosine in Mammalian DNA by the MLL Partner TET1 , 2009 .

[29]  Ken W. Y. Cho,et al.  Gastrulation in the mouse: The role of the homeobox gene goosecoid , 1992, Cell.

[30]  S. Balasubramanian,et al.  Quantitative Sequencing of 5-Methylcytosine and 5-Hydroxymethylcytosine at Single-Base Resolution , 2012, Science.

[31]  Tsutomu Chiba,et al.  Induction and monitoring of definitive and visceral endoderm differentiation of mouse ES cells , 2005, Nature Biotechnology.

[32]  M. Equihua,et al.  Analysis of finite mixture of distributions: a statistical tool for biological classification problems , 1988, Comput. Appl. Biosci..

[33]  D. Lees-Murdock,et al.  DNA methyltransferase expression in the mouse germ line during periods of de novo methylation , 2005, Developmental dynamics : an official publication of the American Association of Anatomists.

[34]  Philipp Kapranov,et al.  Genome-wide mapping of 5-hydroxymethylcytosine in embryonic stem cells , 2011, Nature.

[35]  Johannes Söding,et al.  Different Binding Properties and Function of CXXC Zinc Finger Domains in Dnmt1 and Tet1 , 2011, PloS one.

[36]  Timothy J. Durham,et al.  "Systematic" , 1966, Comput. J..

[37]  H. Blau,et al.  DNA Demethylation Dynamics , 2011, Cell.

[38]  Julia A. Lasserre,et al.  Histone modification levels are predictive for gene expression , 2010, Proceedings of the National Academy of Sciences.

[39]  Yi Zhang,et al.  Role of Tet proteins in 5mC to 5hmC conversion, ES-cell self-renewal and inner cell mass specification , 2010, Nature.

[40]  Haifan Lin,et al.  The biogenesis and function of PIWI proteins and piRNAs: progress and prospect. , 2009, Annual review of cell and developmental biology.

[41]  Zhiping Weng,et al.  Transcription factor binding and modified histones in human bidirectional promoters. , 2007, Genome research.

[42]  Pengzhi Yu,et al.  Dynamic chromatin states in human ES cells reveal potential regulatory sequences and genes involved in pluripotency , 2011, Cell Research.

[43]  Michael A. Beer,et al.  Predicting Gene Expression from Sequence , 2004, Cell.

[44]  Suhua Feng,et al.  5-Hydroxymethylcytosine is associated with enhancers and gene bodies in human embryonic stem cells , 2011, Genome Biology.

[45]  R. Durbin,et al.  Biological sequence analysis: Background on probability , 1998 .

[46]  G. R. Wyatt,et al.  A New Pyrimidine Base from Bacteriophage Nucleic Acids , 1952, Nature.

[47]  Kit T. Rodolfa,et al.  Sox17 promotes differentiation in mouse embryonic stem cells by directly regulating extraembryonic gene expression and indirectly antagonizing self-renewal. , 2010, Genes & development.