PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Cell Rep. Author manuscript; available in PMC 2013 July 3.
Published in final edited form as:
PMCID: PMC3700580
NIHMSID: NIHMS475487

Bivalent Chromatin Marks Developmental Regulatory Genes in the Mouse Embryonic Germline in Vivo

Introduction

Pluripotency is dependent on the maintenance of a proper epigenetic landscape (Gaspar-Maia et al., 2011; Orkin and Hochedlinger, 2011). Bivalent domains, which are defined by the paradoxical coexistence of a permissive histone mark (H3K4me3) and a repressive mark (H3K27me3), are thought to play an important role in pluripotency by keeping developmental genes in a silenced state poised for activation upon differentiation (Azuara et al., 2006; Bernstein et al., 2006). In support of this model, bivalent domains have been shown to be present preferentially at developmental regulatory genes in undifferentiated embryonic stem cells (ESCs) and several adult tissues in vivo including sperm, testis, the cerebellum, and the hematopoietic compartment (Mikkelsen et al., 2007; Cui et al., 2009, 2012; Hammoud et al., 2009).The moderate de-repression of developmental regulators in eed mutant ESCs, which have reduced levels of H3K27me3, further supports this model (Azuara et al., 2006; Boyer et al., 2006). However, it has recently been proposed that the concomitant presence of H3K4me3 and H3K27me3 at developmental genes in ESCs is an in vitro artifact resulting from suboptimal culture conditions (Hong et al., 2011; Marks et al., 2012). In addition, direct evidence for bivalency in the developing mammalian embryo is limited and conflicting due to technical difficulties associated with the low amounts of material available (Alder et al., 2010; Dahl et al., 2010; Rugg-Gunn et al., 2010). The universal nature of bivalency has been further questioned due to conflicting reports from non-mammalian species, with bivalent domains demonstrated to exist in zebrafish (Vastenhouw et al., 2010) but not detected in Xenopus or Drosophila embryos (Akkers et al., 2009; Schuettengruber et al., 2009). We therefore lack a clear understanding of whether bivalency exists in embryonic cells in vivo, and how it relates to pluripotency.

Of particular interest to pluripotency are primordial germ cells (PGCs), the embryonic precursors of the germline. PGCs have a transcriptional profile similar to ESCs (Qin et al., 2012), including the silencing of developmental regulators, and can give rise to pluripotent stem cells when cultured in vitro. However, unlike ESCs, PGCs are unipotent and unable to contribute to chimeras. Additionally, PGCs are thought to undergo a period of ‘epigenetic erasure’ or reprogramming via global removal of both DNA methylation and histone modifications that may prevent transmission of epimutations (Greer and Shi, 2012; Seisenberger et al., 2013). Due to the low number of PGCs present during development, nearly all evidence for this model is from immunofluorescence (IF)-based methodologies that lack gene level resolution and have shown conflicting results (Seki et al., 2007; Hajkova et al., 2008; Kagiwada et al., 2012). To characterize the chromatin state of PGCs, we developed a low cell number Chromatin Immunoprecipitation (ChIP) for the analysis of histone marks using less than 10,000 cells per IP without the need for carrier chromatin or pre-amplification. Using this technique, we performed ChIP-Seq and ChIP-qPCR to show that bivalent domains are present at developmental regulatory genes at multiple stages of PGC development.

Results

Low cell number ChIP protocol

We developed a low cell number ChIP protocol using a magnetic bead-based capture with key improvements over existing approaches (O’Neill et al., 2006; Acevedo et al., 2007; Dahl and Collas, 2008; Adli and Bernstein, 2011; Shankaranarayanan et al., 2011), including reduced handling of material for cross-linking, improvements to sonication reproducibility, and no requirement for carrier DNA or DNA pre-amplification (see Methods). We validated the protocol using a transgenic mouse ESC line expressing an N-terminal myc-tagged H3.3 histone, due to technical simplicity and the availability of genome-wide data (Goldberg et al., 2010). ChIP was performed for myc-H3.3 and a control IgG with decreasing amounts of chromatin starting from 2.5 million to 12,500 cells per IP (Figure S1A). qPCR was used to interrogate the enrichment status of 15 different loci representing 6 categories of genomic regions. Importantly, 11 of the 12 loci expected to be bound by H3.3 (Goldberg et al., 2010) were as least 10 fold more enriched for myc-H3.3 compared to IgG, demonstrating the specificity and low background of the assay. Additionally, the relative enrichment between loci remains stable independent of starting cell number, demonstrating the robustness and sensitivity of the assay (Figure S1B). We envision that this protocol will be broadly useful to characterize the chromatin state of rare stem and precursor cell populations in vivo.

Genome-wide analysis of H3K4me3 and H3K27me3 in PGCs

To investigate the presence of bivalent domains in PGCs, we used a transgenic Oct4:GFP mouse line to isolate PGCs (GFP+) and the surrounding somatic cells (GFP−) by FACS at different stages of development as previously reported (Yoshimizu et al., 1999). Using our low cell ChIP protocol, we performed ChIP-Seq for H3K4me3 and H3K27me3 in two biological replicates of E11.5 PGCs, obtaining a total of approximately 386 million reads and between 15 and 32 million uniquely mappable reads post filtering per ChIP library (Table S1). The two biological replicates in our ChIP-Seq libraries show very strong correlation, with Pearson correlation coefficients of 0.978 for H3K4me3 and 0.996 for H3K27me3 (Figure 1A, S2A, Table S1), indicative of the high reproducibility of the protocol. Furthermore, all IPs are enriched for signal over input as assessed using CHANCE (Table S1) (Diaz et al., 2012).

Figure 1
Bivalent genes in E11.5 PGCs are enriched for developmental regulators

We identified H3K4me3-only regions, H3K27me3-only regions, and bivalent domains in both E11.5 PGCs (this study) and ESCs (Mikkelsen et al., 2007) using ChromHMM (Ernst and Kellis, 2012). The overall distribution of these histone marks is similar between PGCs and ESCs, except for an increase in H3K27me3-marked DNA in PGCs (Figure S2B), which is in agreement with IF data (Seki et al., 2005, 2007; Kagiwada et al., 2012). Genes marked only with H3K27me3 in PGCs (Supplemental File 2) are involved in cell motility (Figure S2C), which may be a consequence of PGCs having just completed their migration to the gonads at E11.5.

Bivalent genes in PGCs are enriched for developmental regulators

We next focused on bivalent regions. Consistent with their association with gene promoters, bivalent domains largely occur at transcriptional start sites (TSSs) in both E11.5 PGCs and ESCs (Figure S2D). E11.5 PGCs and ESCs have a similar distribution and a highly significant overlap of genes marked as bivalent (p<2.2X10−16, Figure 1B, 1C, S2E). Similarly to what has been described for ESCs (Bernstein et al., 2006), the 2715 genes marked as bivalent in PGCs and ESCs are strongly enriched for regulators of early development (Figure 1D). We note that while this manuscript was under review, Ng et al. reported a ChIP-Seq dataset for several histone marks in PGCs, including H3K4me3 and H3K27me3 (Ng et al., 2013). However, the authors did not analyze bivalency in detail and the low signal for H3K4me3 in E11.5 PGCs in that study precludes a direct comparison with our work (Figure S2A). Therefore, we sought to validate our data using ChIP-qPCR in independent biological samples of PGCs. ChIP-qPCR at 12 genes representing 8 developmental pathways confirmed that all are bivalently marked in PGCs, with bivalency defined as enrichment above a background of 10% (Supplemental File 1) for both H3K4me3 and H3K27me3 and a maximum 2-fold difference in enrichment between the two marks (Figure 1E). These bivalent regulators include transcription factors that are important for the specification of the 3 somatic germ layers, such as Pax6, Nkx2.2, MyoD, brachyury (T), Gata6 and FoxA1, and Hox genes such as HoxA3 and HoxB9. This is notable because ESCs are pluripotent and expected to maintain these genes in a poised transcriptional state for activation upon differentiation, but PGCs are unipotent and will not express any of these developmental regulators until post fertilization, in the next generation. In fact, the Hox cluster has been shown to be silenced as part of the early specification of PGC precursors at E7.25 (Saitou et al., 2002). A comparison of the expression pattern of genes bivalent in PGCs and ESCs demonstrates that these genes are transcriptionally repressed in ESCs and PGCs, but subsets are active in select lineage-restricted cell types such as MEFs, bone marrow, or brain tissue (Figure 1F). Taken together, our genome-wide analysis reveals that the mouse embryonic germline in vivo is remarkably similar to cultured ESCs with regards to the presence of bivalent chromatin at silenced developmental regulatory genes.

PGCs are enriched for bivalency at developmental genes compared to the surrounding soma

We next compared the chromatin landscape of PGCs to the surrounding soma. Using ChIP-qPCR, we examined the 12 previously tested developmental regulators plus Pou5f1/Oct4 in ESCs and E11.5 soma (Figure 2A, 2B). 12/12 of these developmental regulators are bivalent in both E11.5 PGCs and ESCs. However, in the soma most of these genes have stronger enrichment for either H3K4me3 or H3K27me3 and only four genes are bivalent (Figure 2B, 2C). These data indicate that developmental regulators of all somatic lineages exist in a bivalent state in ESCs and PGCs but not in somatic cells in vivo (Figure 2D), suggesting that this distribution of bivalent domains is a distinguishing characteristic of pluripotency-associated cells both in vitro and in vivo.

Figure 2
PGCs are enriched for bivalent developmental regulators compared to surrounding somatic cells at E11.5

Developmental regulators remain bivalent in PGCs during the initiation of sexual differentiation

E11.5 PGCs are sexually indifferent and capable of giving rise to pluripotent stem cells when cultured in vitro. After E11.5, PGCs initiate the process of sexual differentiation and express sex-specific genes by E13.5. Furthermore, E13.5 PGCs can no longer give rise to pluripotent stem cells (Labosky et al., 1994). We therefore sought to investigate whether the process of sexual differentiation and concomitant loss of the ability to give rise to pluripotent stem cell cells affects the bivalent chromatin state of PGCs. Using low cell ChIP-qPCR, we examined 8 representative developmental genes through a time-course of PGC development from E11.5 to E13.5. In agreement with previously published data (Qin et al., 2012), somatic developmental genes, including HoxA11, HoxB9, Pax6, Nkx2.2, T, MyoD, Gata6, and FoxA1 are expressed at very low to undetectable levels from E11.5 through E13.5 (Figure 3). Interestingly, these developmental regulators remain bivalent in PGCs throughout this time-period (Figure 3). In contrast, germline-specific genes, such as Dazl, Dppa3, and Vasa are enriched for only H3K4me3, coinciding with their transcriptional activation (Figure S3). Taken together, these data indicate that somatic developmental regulators remain bivalent in the mouse embryonic germline from E11.5 through initiation of sexual differentiation at E13.5.

Figure 3
Developmental regulators remain bivalent and transcriptionally repressed during sexual differentiation in the germline

Orphan bivalent domains are silent in PGCs and transcribed in a tissue-specific manner

Surprisingly, we identified a large number of ‘orphan’ bivalent domains that cannot be associated with any annotated promoter or gene in the RefSeq database, and are distal to annotated genes by a median distance of 13kb (Figure 4A). From a total of 9132 bivalent domains in E11.5 PGCs, 2886 are orphan, and many are also bivalent in ESCs (Table S2). Bivalent domains are known to have a strong association with unmethylated CpG islands (CGIs) (Ku et al., 2008). We found that 1413 (49.0%, p-value < 1×10−5) of the orphan bivalent domains in PGCs map to experimentally defined unmethylated CGIs (Illingworth et al., 2010). We next investigated the expression status of the regions containing orphan bivalent domains. If the orphan bivalent domains were similar to bivalent genes (Figure 1F), they would be silent in PGCs and expressed in somatic tissues. Indeed, only 23 (0.80%) orphan bivalent domains show signs of transcription in an RNA-Seq dataset from E11.5 PGCs (Seisenberger et al., 2012). Analysis of CAGE tag clusters showed that 927 of the 1413 (65.6%, p-value = 5.78×10−168) orphan bivalent domains with CGIs are expressed in at least 1 of 22 tissue types (Figure 4B, Table S3) (Kawaji et al., 2009). Similarly, 461 of the 1473 (31.3%, p-value = 3.61×10−35) orphan bivalent domains without CGIs are expressed in at least 1 of 22 tissue types (Figure 4B, Table S3). Furthermore, both classes of orphan bivalent domains (with or without CGI) display highly tissue-specific expression (Figure 4C), more so than bivalent developmental regulators (p-value = 1.1×10−141, 2.3×10−72 respectively), with a preferential activation in brain regions (Figure 4E). Interestingly, we found bivalent domains in PGCs at 300 non-coding RefSeq genes and 179 putative non-coding RNAs inferred from chromatin state (Guttman et al., 2009) and expression data (Carninci et al., 2005) (Figure 4D, S2F). These data suggest that, similar to protein-coding developmental regulators, some non-coding RNAs are bivalent in the germline and poised for activation upon lineage commitment.

Figure 4
Orphan bivalent domains overlap CpG Islands and show tissue specific expression

Discussion

In this work we describe an optimized ChIP protocol suitable for low numbers of cells and its use to examine the genome-wide and gene-specific distribution of H3K4me3 and H3K27me3 in mouse PGCs. We report that developmental regulatory genes remain bivalent and transcriptionally silent in vivo in PGCs, but not in adjacent somatic cells, throughout E11.5-E13.5, in a manner highly similar to cultured ESCs. In addition to developmental genes, we identify some ~3000 orphan bivalent domains that are enriched for CGIs and are expressed in a tissue-specific manner. Some of the bivalent domains identified here correspond to non-coding RNAs that we speculate may have developmental functions, similar to bivalent genes, although this remains to be determined. The findings presented here represent strong evidence for the existence of bivalent domains in the embryonic germline, and raise a number of important questions.

Do bivalent promoters escape the epigenetic reprogramming reported to occur in PGCs from E8.5–E13.5? Our data support the suggestion that H3K27me3 may act as a potential mechanism to compensate for the loss of DNA methylation or H3K9me2/3 in mid-gestation PGCs (Sasaki and Matsui, 2008). It is also possible that the loss of certain histone marks observed by IF may occur primarily at abundant repetitive sequences of the genome and mostly spare unique genes, such as developmental regulators. The application of the low cell ChIP protocol reported here should help answer this question for other histone marks and time-points.

Does the germline continuously maintain the somatic program in a bivalent state until activation in the next generation? It is unclear why such a large set of important developmental regulators would remain bivalent in the developing germline, but not in somatic cells. The hypothesis that bivalency maintains developmental genes in a transcriptionally poised state is intuitive in pluripotent cells that activate these genes upon differentiation, such as ESCs (Azuara et al., 2006; Bernstein et al., 2006) or the epiblast (Rugg-Gunn et al., 2010). The presence of bivalent domains in PGCs may contribute to preventing expression of the somatic program in the germline. Our observations suggest that developmental regulators may be kept in a repressed but accessible state in the germline for activation post-fertilization, in the next generation. The transmission of bivalency through the germline could provide a substratum for epigenetic inheritance. Surprisingly, recent work indicates that the sperm genome maintains a residual level of nucleosomes, and that these are enriched for H3K4me3/H3K27me3 bivalent marks at developmental regulators (Hammoud et al., 2009). It will be interesting to determine whether bivalency is detected at other stages of germline development and in oocytes. Functional studies of regulators of bivalency in PGCs should provide important insights into these questions.

Experimental Procedures

Isolation of PGCs

Male B6 mice homozygous for a transgenic Oct4ΔPE:GFP reporter were crossed with Swiss-Webster females. The gonadal regions from multiple embryos were isolated and pooled prior to enzymatic dissociation for FACS.

Low Cell ChIP-qPCR

Cells were pooled into batches of ~50,000 cells and cross-linked in 1% formaldehyde (Sigma F8775-25ml) in PBS for 5min at room temperature and quenched with 125mM glycine. For each IP, 11µl of protein A Dynabeads (Life Technologies 1001D) were pre-incubated with 2.4µg of either anti H3K4me3 (Diagenode pAb003-050), H3K27me3 (Millipore 07-449), or IgG (Abcam ab46540) in ChIP lysis buffer. Cross-linked cells were lysed in Lo-Bind microfuge tubes (Eppendorf 022431021) followed by a 40 minute sonication in a BioRuptor sonicator (Diagenode UCD-200) set to high power, 7” ON, 15” OFF, changing ice/water slurry every 10 minutes.

Lysate was diluted in lysis buffer and cleared of debris. Cleared lysate from ~50,000 cells was divided into four equal aliquots, and one aliquot used for Input. The remaining three aliquots (equivalent to ~12,500 cells each) were used for IP of H3K4me3, H3K27me3 and IgG. Lysates were incubated with preformed bead/antibody overnight at 4C with mixing. The chromatin/bead/antibody complexes were washed sequentially with lysis buffer three times, DOC buffer once, and TE buffer once, followed by a transfer to a fresh PCR tube. Chromatin was eluted using elution buffer (1% SDS, 0.1M NaHCO3). IP and Input samples were treated with RNaseA followed by Proteinase K treatment. Cross-linking was reversed by incubating overnight at 65C while shaking. DNA purification was done using a QiaQuick PCR Purification Kit (Qiagen 28104). E11.5 ChIP data are representative of 3–7 biological replicates expect for Dnmt3b which was technically replicated. Additional PGC data represent a minimum of 2 biological replicates except for male E12.5 which were technically replicated. All primers are listed in Supplemental File 1. A full protocol is located in the supplemental experimental procedures.

Low Cell ChIP-Seq

ChIP material was obtained and processed as described above with the following modifications. Cells were cross-linked in 0.25% formaldehyde (Thermo Scientific 28906) in PBS for 10 minutes at room temperature prior to quenching with 125mM glycine. Cross-linked material was sonicated using a Covaris sonicator for 12’ at Duty 5%, Intensity 3, and Bursts 200. Two aliquots of E11.5 PGCs, consisting of 104,000 and 97,000 cells, were each divided equally to perform ChIP for H3K4me3, H3K27me3, and an input control. IP-DNA and Inputs were purified using a MinElute PCR Purification Kit (Qiagen 28004) and libraries generated using the ThruPLEX-FDPrep Kit (Rubicon Genomics R40048) with 20 cycles of amplification for IP-DNA and 15 cycles for Input DNA.

ChIP-Seq Analysis

Reads from each ChIP-Seq library were filtered to retain only unique sequences (Table S1). Reads were aligned to the mm9/NCBI build 37 mouse genome using bowtie and reads mapping only once were retained. Data for ESCs were handled in the same manner, obtained from GEO accession GSE12241(Mikkelsen et al., 2007). To call and compare bivalent domains between the PGC and ESC datasets, a 13-state segmentation of the genome was generated using ChromHMM (v1.06) (Ernst and Kellis, 2012). The four samples (PGC H3K4me3, PGC H3K27me3, ESC H3K4me3, ESC H3K27me3) were treated as distinct marks from a single cell type to allow a direct and unbiased comparison of the genome segmentation in PGCs and ESCs. K-means clustering was used to group the state emission parameters for each sample into 2 groups ("on" or "off"), and these were used to classify each of the 13 ChromHMM segmentation states as "bivalent", "H3K4me3 only", "H3K27me3 only", or "none" for each cell type. Data from the replicate PGC libraries were pooled together for this analysis. H3 ChIP-Seq data from the Mikkelsen dataset was used as control for the ESC samples.

All gene set analysis is based on the mm9 RefSeq gene set. To ensure a conservative segregation of bivalent domains associated with regulation of annotated genes from orphan bivalent domains, genes were considered bivalent if an H3K4me3/H3K27me3 enriched region fell in the range 0.2kb upstream of the transcription start site (TSS) through the transcription end site (TES) of any isoform of the gene. H3K4me3/H3K27me3 enriched regions occurring outside the range from 0.2kb upstream of the TSS through the TES of any isoform of a gene were called orphans. Genes were considered H3K4me3-only if no isoforms were bivalent and if any isoform had H3K4me3-only regions falling in the window from − 0.2 kb upstream to 0.2 kb downstream of the TSS. Similarly, genes were considered to be H3K27me3-only if they were neither bivalent nor H3K4me3-only and if any isoform had regions of H3K27me3-only occurring in the +/− 0.2kb window. Analysis of enrichment for Gene Ontology terms was done using DAVID (Huang et al., 2009). The Gene Expression Omnibus accession number for the ChIIP-Seq data reported in this paper is GSE46396.

CAGE Analysis

CAGE tag cluster locations in mm9 and tags per million (tpm) values for 22 tissue types were obtained from FANTOM 4 (Kawaji et al., 2009). To assess the statistical significance of the accumulation of CAGE tags, the probability pi of observing a CAGE tag within 200bp of each orphan bivalent domain was first estimated based on CAGE tag representation in the region Xi containing the bivalent domain extended 10kb in each direction. Let Yi be the indicator random variable that indicates whether the observed bivalent domain in Xi had a detectable CAGE tag cluster within 200 bp. Assuming that Yi are independent Bernoulli random variables with probability pi, the sum of Yi is approximately normal with mean Σipi and variance Σipi(1−pi), via the Lyapunov central limit theorem (see supplemental experimental procedures).

Tissue Specificity Calculation

For each CAGE tag cluster, tissue specificity score was computed based on the Jensen-Shannon divergence between the relative abundance of tpm values across the tissue types and the extreme distribution of being expressed in only one tissue type where the tag cluster has the greatest expression value (Cabili et al., 2011).

Highlights

  • A ChIP-Seq/qPCR protocol for the analysis of histone marks using <10,000 cells
  • E11.5 PGCs have bivalent domains at developmental regulators, remarkably similar to ESCs
  • Developmental genes remain bivalent and silent from E11.5 to E13.5 in PGCs
  • Bivalent domains distal to known genes mark tissue-specific, un-annotated promoters

Supplementary Material

01

02

03

04

Click here to view.(1019K, xlsx)

05

Acknowledgements

MS conceived of the project. MRS and MS directed the project with participation from JSS. MS and MRS developed the study design with participation from KB and KTE. MS developed the low-cell ChIP protocol, performed ChIP and expression experiments and analyzed them. MS, KB, and KTE isolated PGCs. KTE conducted all FACS. MS generated ChIP-Seq libraries. Sequencing was done at UC DAVIS Expression Analysis Core. CO analyzed the majority of genome-wide data with aid and supervision from JSS. MS and MRS wrote the manuscript. All authors read, helped write, and edit the manuscript. We thank Jehnna Ronan for advice in developing the ChIP assay, members of the Santos lab for discussions and Robert Blelloch, Diana Laird, Marco Conti and Barbara Panning for critical reading of the manuscript. K.B. was the recipient of an NSF pre-doctoral fellowship. K.T.E. was the recipient of a CIRM post-doctoral training grant (TG2-01153). This work was funded by U54HD055764 to J.S.S. and an NIH New Innovator Award (DP2OD004698) to M. R.-S.

Footnotes

Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

References

  • Acevedo L, Leonardo Iniguez a, Holster H, Zhang X, Green R, Farnham P. Genome-scale ChIP-chip analysis using 10,000 human cells. BioTechniques. 2007;43:791–797. [PMC free article] [PubMed]
  • Adli M, Bernstein BE. Whole-genome chromatin profiling from limited numbers of cells using nano-ChIP-seq. Nature Protocols. 2011;6:1656–1668. [PubMed]
  • Akkers RC, Van Heeringen SJ, Jacobi UG, Janssen-Megens EM, Françoijs K-J, Stunnenberg HG, Veenstra GJC. A hierarchy of H3K4me3 and H3K27me3 acquisition in spatial gene regulation in Xenopus embryos. Developmental Cell. 2009;17:425–434. [PMC free article] [PubMed]
  • Alder O, Lavial F, Helness A, Brookes E, Pinho S, Chandrashekran A, Arnaud P, Pombo A, O’Neill L, Azuara V. Ring1B and Suv39h1 delineate distinct chromatin states at bivalent genes during early mouse lineage commitment. Development (Cambridge, England) 2010;137:2483–2492. [PubMed]
  • Azuara V, Perry P, Sauer S, Spivakov M, Jørgensen HF, John RM, Gouti M, Casanova M, Warnes G, Merkenschlager M, et al. Chromatin signatures of pluripotent cell lines. Nature Cell Biology. 2006;8:532–538. [PubMed]
  • Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Fry B, Meissner A, Wernig M, Plath K, Jaenisch R, et al. A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006;125:315–326. [PubMed]
  • Boyer LA, Plath K, Zeitlinger J, Brambrink T, Medeiros LA, Lee TI, Levine SS, Wernig M, Tajonar A, Ray MK, et al. Polycomb complexes repress developmental regulators in murine embryonic stem cells. Nature. 2006;441:349–353. [PubMed]
  • Cabili MN, Trapnell C, Goff L, Koziol M, Tazon-Vega B, Regev A, Rinn JL. Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. Genes & Development. 2011;25:1915–1927. [PubMed]
  • Carninci P, Kasukawa T, Katayama S, Gough J, Frith MC, Maeda N, Oyama R, Ravasi T, Lenhard B, Wells C, et al. The transcriptional landscape of the mammalian genome. Science (New York, N.Y.) 2005;309:1559–1563. [PubMed]
  • Cui K, Zang C, Roh T-Y, Schones DE, Childs RW, Peng W, Zhao K. Chromatin signatures in multipotent human hematopoietic stem cells indicate the fate of bivalent genes during differentiation. Cell Stem Cell. 2009;4:80–93. [PMC free article] [PubMed]
  • Cui P, Liu W, Zhao Y, Lin Q, Zhang D, Ding F, Xin C, Zhang Z, Song S, Sun F, et al. Comparative Analyses of H3K4 and H3K27 Trimethylations Between the Mouse Cerebrum and Testis. Genomics, Proteomics & Bioinformatics / Beijing Genomics Institute. 2012;10:82–93. [PubMed]
  • Dahl JA, Collas P. A rapid micro chromatin immunoprecipitation assay (microChIP) Nature Protocols. 2008;3:1032–1045. [PubMed]
  • Dahl JA, Reiner AH, Klungland A, Wakayama T, Collas P. Histone H3 lysine 27 methylation asymmetry on developmentally-regulated promoters distinguish the first two lineages in mouse preimplantation embryos. PloS One. 2010;5:e9150. [PMC free article] [PubMed]
  • Diaz A, Nellore A, Song JS. CHANCE: comprehensive software for quality control and validation of ChIP-seq data. Genome Biology. 2012;13:R98. [PMC free article] [PubMed]
  • Ernst J, Kellis M. ChromHMM: automating chromatin-state discovery and characterization. Nature Methods. 2012;9:215–216. [PMC free article] [PubMed]
  • Gaspar-Maia A, Alajem A, Meshorer E, Ramalho-Santos M. Open chromatin in pluripotency and reprogramming. Nature Reviews. Molecular Cell Biology. 2011;12:36–47. [PMC free article] [PubMed]
  • Goldberg AD, Banaszynski LA, Noh KM, Lewis PW, Elsaesser SJ, Stadler S, Dewell S, Law M, Guo X, Li X, et al. Distinct factors control histone variant H3.3 localization at specific genomic regions. Cell. 2010;140:678–691. [PMC free article] [PubMed]
  • Greer EL, Shi Y. Histone methylation: a dynamic mark in health, disease and inheritance. Nature Reviews. Genetics. 2012;13:343–357. [PMC free article] [PubMed]
  • Guttman M, Amit I, Garber M, French C, Lin MF, Huarte M, Zuk O, Carey BW, Cassady JP, Moran N, et al. Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. 2009;458:223–227. [PMC free article] [PubMed]
  • Hajkova P, Ancelin K, Waldmann T, Lacoste N, Lange UC, Cesari F, Lee C, Almouzni G, Schneider R, Surani MA. Chromatin dynamics during epigenetic reprogramming in the mouse germ line. Nature. 2008;452:877–881. [PMC free article] [PubMed]
  • Hammoud SS, Nix Da, Zhang H, Purwar J, Carrell DT, Cairns BR. Distinctive chromatin in human sperm packages genes for embryo development. Nature. 2009;460:473–478. [PMC free article] [PubMed]
  • Hong S-H, Rampalli S, Lee JB, McNicol J, Collins T, Draper JS, Bhatia M. Cell fate potential of human pluripotent stem cells is encoded by histone modifications. Cell Stem Cell. 2011;9:24–36. [PubMed]
  • Huang DW, Sherman BT, Lempicki Ra. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Research. 2009;37:1–13. [PMC free article] [PubMed]
  • Illingworth RS, Gruenewald-Schneider U, Webb S, Kerr ARW, James KD, Turner DJ, Smith C, Harrison DJ, Andrews R, Bird AP. Orphan CpG islands identify numerous conserved promoters in the mammalian genome. PLoS Genetics. 2010;6 [PMC free article] [PubMed]
  • Kagiwada S, Kurimoto K, Hirota T, Yamaji M, Saitou M. Replication-coupled passive DNA demethylation for the erasure of genome imprints in mice. The EMBO Journal. 2012;32:340–353. [PubMed]
  • Kawaji H, Severin J, Lizio M, Waterhouse A, Katayama S, Irvine KM, Hume Da, Forrest ARR, Suzuki H, Carninci P, et al. The FANTOM web resource: from mammalian transcriptional landscape to its dynamic regulation. Genome Biology. 2009;10:R40. [PMC free article] [PubMed]
  • Ku M, Koche RP, Rheinbay E, Mendenhall EM, Endoh M, Mikkelsen TS, Presser A, Nusbaum C, Xie X, Chi AS, et al. Genomewide analysis of PRC1 and PRC2 occupancy identifies two classes of bivalent domains. PLoS Genet. 2008;4:e1000242. [PMC free article] [PubMed]
  • Labosky Pa, Barlow DP, Hogan BL. Mouse embryonic germ (EG) cell lines: transmission through the germline and differences in the methylation imprint of insulin-like growth factor 2 receptor (Igf2r) gene compared with embryonic stem (ES) cell lines. Development (Cambridge, England) 1994;120:3197–3204. [PubMed]
  • Marks H, Kalkan T, Menafra R, Denissov S, Jones K, Hofemeister H, Nichols J, Kranz A, Francis Stewart A, Smith A, et al. The Transcriptional and Epigenomic Foundations of Ground State Pluripotency. Cell. 2012;149:590–604. [PMC free article] [PubMed]
  • Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, Alvarez P, Brockman W, Kim T-KK, Koche RP, et al. Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 2007;448:553–560. [PMC free article] [PubMed]
  • Ng J-H, Kumar V, Muratani M, Kraus P, Yeo J-C, Yaw L-P, Xue K, Lufkin T, Prabhakar S, Ng H-H. In Vivo Epigenomic Profiling of Germ Cells Reveals Germ Cell Molecular Signatures. Developmental Cell. 2013;24:324–333. [PubMed]
  • O’Neill LP, VerMilyea MD, Turner BM. Epigenetic characterization of the early embryo with a chromatin immunoprecipitation protocol applicable to small cell populations. Nature Genetics. 2006;38:835–841. [PubMed]
  • Orkin SH, Hochedlinger K. Chromatin connections to pluripotency and cellular reprogramming. Cell. 2011;145:835–850. [PubMed]
  • Qin H, Blaschke K, Wei G, Ohi Y, Blouin L, Qi Z, Yu J, Yeh R-F, Hebrok M, Ramalho-Santos M. Transcriptional analysis of pluripotency reveals the Hippo pathway as a barrier to reprogramming. Human Molecular Genetics. 2012;21:2054–2067. [PMC free article] [PubMed]
  • Rugg-Gunn PJ, Cox BJ, Ralston A, Rossant J, Handling E. Distinct histone modifications in stem cell lines and tissue lineages from the early mouse embryo. Proceedings of the National Academy of Sciences of the United States of America. 2010;107:10783–10790. [PubMed]
  • Saitou M, Barton SC, Surani MA. A molecular programme for the specification of germ cell fate in mice. Nature. 2002;418:293–300. [PubMed]
  • Sasaki H, Matsui Y. Epigenetic events in mammalian germ-cell development: reprogramming and beyond. Nature Reviews. Genetics. 2008;9:129–140. [PubMed]
  • Schuettengruber B, Ganapathi M, Leblanc B, Portoso M, Jaschek R, Tolhuis B, Van Lohuizen M, Tanay A, Cavalli G. Functional anatomy of polycomb and trithorax chromatin landscapes in Drosophila embryos. PLoS Biology. 2009;7:e13. [PMC free article] [PubMed]
  • Seisenberger S, Andrews S, Krueger F, Arand J, Walter J, Santos F, Popp C, Thienpont B, Dean W, Reik W. The Dynamics of Genome-wide DNA Methylation Reprogramming in Mouse Primordial Germ Cells. Molecular Cell. 2012;95:1–14. [PMC free article] [PubMed]
  • Seisenberger S, Peat JR, Reik W. Conceptual links between DNA methylation reprogramming in the early embryo and primordial germ cells. Current Opinion in Cell Biology. 2013;25:1–8. [PubMed]
  • Seki Y, Hayashi K, Itoh K, Mizugaki M, Saitou M, Matsui Y. Extensive and orderly reprogramming of genome-wide chromatin modifications associated with specification and early development of germ cells in mice. Developmental Biology. 2005;278:440–458. [PubMed]
  • Seki Y, Yamaji M, Yabuta Y, Sano M, Shigeta M, Matsui Y, Saga Y, Tachibana M, Shinkai Y, Saitou M. Cellular dynamics associated with the genome-wide epigenetic reprogramming in migrating primordial germ cells in mice. Development (Cambridge, England) 2007;134:2627–2638. [PubMed]
  • Shankaranarayanan P, Walia M, Wang L, Li N, Trindade LM, Gronemeyer H. single-tube linear dna amplification ( Linda ) for robust chiP-seq. Nature Methods. 2011;8:565–568. [PubMed]
  • Vastenhouw NL, Zhang Y, Woods IG, Imam F, Regev A, Liu XS, Rinn J, Schier AF. Chromatin signature of embryonic pluripotency is established during genome activation. Nature. 2010;464:922–926. [PMC free article] [PubMed]
  • Yoshimizu T, Sugiyama N, De Felice M, Yeom YI, Ohbo K, Masuko K, Obinata M, Abe K, Schöler HR, Matsui Y. Germline-specific expression of the Oct-4/green fluorescent protein (GFP) transgene in mice. Development, Growth & Differentiation. 1999;41:675–684. [PubMed]