|Home | About | Journals | Submit | Contact Us | Français|
Embryonic stem (ES) cells are derived from blastocyst stage embryos and are believed to be functionally equivalent to the inner cell mass, which lacks the ability to produce all extraembryonic tissues. Here we report the identification of a rare transient cell population within mouse ES and induced pluripotent stem (iPS) cell cultures that express high levels of transcripts found in two-cell (2C) embryos in which the blastomeres are totipotent. We genetically tagged these 2C-like ES cells and show that they lack the ICM pluripotency proteins Oct4, Sox2, and Nanog and have acquired the ability to contribute to both embryonic and extraembryonic tissues. We show that nearly all ES cells cycle in and out of this privileged state, which we find is partially controlled by histone modifying enzymes. Transcriptome sequencing and bioinformatic analyses revealed that a significant number of 2C-transcripts are initiated from long terminal repeats derived from murine endogenous retroviruses, suggesting this foreign sequence has helped to drive cell fate regulation in placental mammals.
The zygote and its daughter cells are totipotent because they are able to develop into all embryonic and extraembryonic cell types1,2. The progeny of these first two daughter cells become progressively more fate restricted as they activate distinct patterns of gene expression that first direct them toward one of three broad lineages: Oct4/Sox2/Nanog+ epiblast cells that give rise to the embryo, Gata4/6+ primitive endoderm cells that contribute to extraembryonic membranes that encase the embryo, and Cdx2+ trophectoderm cells that form a large part of the placenta3. These early cell fate decisions represent a major and relatively recent advance in mammalian evolution in which the placenta and extraembryonic tissues that support the intrauterine nourishment of the fetus allow development to progress further before birth. The epigenetic landscape of the zygote changes dramatically during the first cell divisions. Shortly after fertilization the oocyte maternal transcripts are replaced with newly synthesized RNAs generated by activating transcription of the zygotic genome4–6. The unique transcriptional profile of the zygote and its daughter cells defines a brief period when the cells are totipotent.
Murine ES cells are isolated from the inner cell mass (ICM) of blastocysts that have already become a separate lineage from the trophectoderm7,8. ICM-derived ES cells are regarded as pluripotent because they have the capacity to generate tissues of the fetus but are extremely inefficient at colonizing the extraembryonic tissues9. The rare contribution of ES cells to extraembryonic tissues could be explained by contamination of ES cultures with trophectoderm or primitive endoderm-committed cells, or occur because rare ES cells have acquired the ability to produce extra-embryonic tissues in addition to embryonic tissues. This latter possibility is intriguing because recent evidence shows that ES cultures are a heterogeneous mixture of metastable cells with fluctuating expression of genes such as Zscan4, Stella, Nanog, Sox17 and Gata6, which could account for special attributes of individual cells10–14.
A large number of retrotransposons are expressed when the zygotic genome is first transcribed, including the endogenous retroviruses (ERVs), LINE-1 elements, and the non-autonomous SINE elements15. At the 2C stage, MuERV-L/MERVL retrovirus-like elements are transiently de-repressed and produce 3% of the transcribed mRNAs15–17. Following the 2C stage, MERVL-retroelement expression is silenced18,19. We discovered that this regulated pattern of MERVL expression overlapped with greater than one hundred 2C-specific genes that have co-opted regulatory elements from these foreign retroviruses to initiate their transcription. We exploited the regulated activity of these 2C virus-derived promoters to label cells and found that both ES and iPS cultures contain a small but relatively constant fraction of cells that have entered into the 2C-transcriptional state. Purification of these 2C-like cells reveals that they have unique developmental characteristics and efficiently produce progeny for extraembryonic and embryonic lineages.
To identify zygotically activated genes we performed deep RNA-sequencing on mouse oocytes and 2C stage embryos. A comparison of the transcripts in these cells identified a large number of genes and retrotransposons that became expressed in the 2C embryo as well as numerous transcripts that were downregulated (Fig. 1a, Supplementary Table 1). The most highly activated repeat was the MERVL-family of retroviruses and their corresponding long terminal repeat (LTR) promoters (Mt2_mm) which were activated greater than 300-fold (Supplementary Table 1). Sequence alignments revealed that greater than 25% of the nearly 700 copies of MERVL elements were activated, and that 307 genes generated chimeric transcripts with junctions to MERVL elements (Fig.1a, Supplementary Table 2), including 10 that were previously described15. Of the 626 chimeric transcripts generated, >90% were 5'LTR- exon fusions that generated open reading frames (ORFs), suggesting that these LTRs had become functional promoters for protein coding genes (Fig. 1b, Supplementary Fig. 1a). The most significantly enriched Gene Ontology categories representing these chimeric proteins were regulation of transcription, ion binding, translation, nucleotide binding, and mRNA transport (Fig. 1c). Two notable transcription factors that utilized alternate MERVL-LTR promoters were Gata4 and Tead4, factors important for the specification of primitive endoderm and trophectoderm, respectively20–22.
Since greater than 300 out of nearly 700 copies of the MERVL endogenous retroviruses still encode Gag viral protein, we stained 2C and early blastocyst embryos to confirm that viral Gag was expressed and developmentally regulated. We found that 2C embryos express Gag but lack the pluripotency marker Oct4, whereas blastula cells lack Gag but express Oct4 (Fig. 1d–e). Thus, MERVL activity is developmentally regulated and these retroviral promoters have been co-opted by many cellular genes to impose tight control over their expression.
Next we asked whether it was possible to use the regulatory sequences from MERVL elements to label 2C cells. We cloned the MERVL 5'LTR, primer binding site, and a portion of the Gag gene upstream of a tandem-tomato fluorescent reporter (2Ctomato). We injected fertilized eggs with the 2Ctomato construct and monitored the expression of tomato during culture in vitro. Tomato expression was highest in arrested zygotes and 2C embryos and became downregulated at the morula stage (Fig. 1f, Supplementary Movie 1). Interestingly, when we introduced the 2Ctomato construct into ES cells and selected for clonal stable integrants, we found several colonies that contained 1–5 cells that were strongly labeled with tomato amongst cells lacking expression of the reporter (Fig. 1g). Importantly, we also found that rare ES cells expressed MERVL mRNA and Gag protein, and that these overlapped with 2Ctomato+ cells (Fig. 1g, Supplementary Fig. 1b–c). The correspondence between the 2Ctomato reporter and the expression of MERVL was further confirmed by immunoblotting, and EM imaging of viral epsilon particles encoded by MERVL within the endoplasmic reticulum of tomato+ cells but not tomato− cells (Supplementary Fig. 1d, e). Thus MERVL expression is restricted in vivo to 1-4C embryos and is reactivated within a small subpopulation of ES cells derived from blastocysts.
To characterize the unexpected 2Ctomato-labeled cells within ES cultures we sorted tomato+ and tomato− cells and performed microarray and mRNA-sequencing analyses (Fig. 1h, Supplementary Tables 3–5). Tomato+ cells expressed 55 fold higher levels of MERVL transcripts than tomato− cells but the vast majority of other retrotransposons were unaffected (Supplementary Table 3). Strikingly, tomato+ cells had 165 transcripts activated more than 4-fold and zero genes repressed more than 4-fold compared with tomato− cells (Fig. 1h, Supplementary Table 4, Supplementary Fig. 2a–f). Amongst the genes that were highly enriched in tomato+ cells, several were previously shown to be restricted to the 2-4C stage of development, including Zscan4, Tcstv1/3, Eif1a, Tho4 (Gm4340) Tdpoz, and Zfp35223–25. In total, 525 genes that were enriched in 2Ctomato+ cells were also genes activated at the 2C stage, including 52 genes that generated chimeric transcripts linked to MERVL elements (Supplementary Tables 6–7).
A hallmark of the ICM and ES cells is their expression of Oct4, Sox2, and Nanog, whereas totipotent 2C embryos do not express Oct4 (Fig. 1d–e). We found that 2Ctomato+ cells within ES cultures also lacked Oct4, Sox2, and Nanog (Fig. 1i, Supplementary Fig. 1f). The reduction in Oct4, Sox2, and Nanog protein-labeling occurred despite changes in their mRNA levels, suggesting the regulation is occurring post-transcriptionally (Fig. 1h, Supplementary Fig. 2g). In summary, 2Ctomato labels a subset of ES cells that share transcriptional and proteomic features of 2C embryos and display strikingly different patterns of pluripotency markers from the majority of ES cells in culture.
We considered the possibility that the expression of the 2Ctomato reporter and MERVL-Gag protein in sporadic cells within ES cultures might arise from contamination with trophectoderm or primitive endoderm. To exclude this possibility we examined induced pluripotent cells (iPS) derived from mouse fibroblasts since they should not be contaminated with cells from blastocyst embryos. Similar to ES cells, we found that sporadic iPS cells express the MERVL-Gag protein and lack Oct4 (Fig. 1i). Thus, the heterogeneity within ES cultures is a property that is shared with iPS cultures and is unlikely to arise from a cell contaminant.
Next we examined whether the 2Ctomato+ cells represent a stable cell population or whether ES cells transition into and out of this 2C-like state. We utilized a Cre-LoxP fate mapping strategy to indelibly mark cells that had expressed 2C genes (Supplementary Fig. 3a–c). We generated a transgenic mouse line using the MERVL regulatory elements driving expression of a tamoxifen inducible Cre recombinase (2CERT2CreERT2, Supplementary Fig. 3a). These mice were then mated with Cre responsive reporter lines (ROSALSL-tomato and ROSALSL-LacZ, Supplementary Fig. 3b). ES cell lines were derived from double positive transgenic blastocysts (Supplementary Fig. 3c). After addition of 4 hydroxytamoxifen (4OHT) to the ES cultures we detected nuclear Cre expression in MERVL-Gag+ cells (Supplementary Fig. 3d). When ES cultures were grown for 2–6 days with 4OHT we found a steady increase in the percentage of reporter positive cells (Fig. 2a). Remarkably, over extended passages nearly every ES cell activated the reporter (Fig. 2b), demonstrating this transient state is regularly entered by ES cells.
To monitor the kinetics of the interconversion between 2Ctomato+ and− cells we performed flow cytometry to collect tomato+ and tomato− cells. When these purified subpopulations were cultured we found that tomato+ cells produced tomato− cells and vice versa (Fig. 2c). Within 24 hours nearly 50% of tomato+ cells convert to tomato−, independent of the starting percentages of the two cell populations (Fig. 2c and data not shown). Under hypoxic conditions (5% O2), the percentage of cells expressing the 2Ctomato reporter was decreased, which could be reversed by shifting the cultures back to 20% O2 (Fig. 2d). We also found that growing cells for 48 hours “ground state” media conditions (2i-media26) reduced but did not eliminate the presence of tomato+ cells relative to media containing knockout serum replacement (KOSR), suggesting extrinsic and intrinsic mechanisms regulate the MERVL/2C gene network (Fig. 2e).
After activation of the zygotic genome in mouse development, histone deacetylation and histone H1 synthesis leads to the formation of repressive chromatin which is thought to limit the broad pattern of transcription present in 2C embryos27,28. Using indirect immunofluorescence, we found that tomato+ cells had significantly higher levels of active histone marks including methylation of histone 3 lysine 4 (H3 K4) and acetylation of H3 and H4 (Supplementary Fig. 4a), a finding confirmed using immunoblot analysis of sorted cell populations (Fig. 3a). This type of chromatin mirrors that found in two cell embryos28. Next we tested whether tomato+ cells had different levels of DNA methylation compared to non-labeled ES cells. We found that the MERVL sequences were hypomethylated in tomato+ cells compared to tomato− cells. In contrast to the MERVL sequences, IAP retroviruses were highly methylated in both tomato+ and tomato− cells suggesting the altered pattern of methylation was not uniform across the genome (Supplementary Fig. 4b). In sum these data suggest that as ES cells (re)enter into the 2C state, their chromatin and DNA is altered to favor transcription in a way that mirrors the two-cell embryo.
We previously demonstrated that MERVL and 2C-specific genes were increased in the absence of the histone lysine specific demethylase Lsd1/Kdm1a29. To test whether other proviral co-repressors and histone-modifying enzymes also influence these genes we profiled the transcriptome of ES cells with homozygous mutations in the KRAB box associated transcriptional repressor Kap1 and the histone H3 K9 methyltranferase G9a29–31. We found that MERVL and a number of 2C genes were significantly upregulated in Kdm1a, Kap1, and G9a mutant ES cells (Fig. 3b, Supplementary Fig. 5b–c, Supplementary Table 7). These findings were confirmed using in situ hybridization and immunofluorescence microscopy (Fig. 3c, Supplementary Fig. 5d). Treatment of 2Ctomato ES lines with the histone deacetylase inhibitor TSA also increased the number of tomato+ cells four-fold (Fig. 3d). To better understand how 2C gene regulation is controlled when chromatin repressors are acutely eliminated we utilized a stably integrated 2Ctomato ES line that is homozygous for a floxed allele of Kdm1a and contains a Cre-ERT transgene that can be activated with 4OHT. Within 24 hours of deleting Kdm1a we found a ten-fold increase in tomato+ cells that was steadily maintained (Fig. 3e). In addition, FACS purified tomato− cells more rapidly became tomato+ in the absence of Kdm1a (Fig. 3f), and stayed in this state longer (Supplementary Movies 2–3). These findings suggest that Kdm1a, Kap1, G9a, and histone deacetylases all contribute to the repression of 2C genes in ES cells and that they function by altering the equilibrium between the 2Ctomato+ and− states.
Since 2C-like cells within ES cultures express high levels of 2C-restricted genes found in totipotent blastomeres and reduced levels of pluripotency associated proteins, we reasoned that this subpopulation of ES cells might have distinct functional characteristics. We tested whether 2Ctomato+ cells have acquired the ability to produce extraembryonic tissues, a characteristic that ES cells lack. We used FACS to collect tomato+ and tomato− cells from a 2Ctomato ES line and injected four cells into morula stage embryos. The tomato− cells contributed exclusively to the ICM of all 5 chimeric blastocysts analyzed (Fig. 4a). In contrast, the tomato+ cells contributed to the trophectoderm (in 4 of 5 chimeric embryos) in addition to the ICM (in 3 of 5 chimeric embryos) (Fig. 4a). To track the fate of the 2Ctomato+ ES cells later in development, we injected blastocysts with tomato+ or tomato− cells that were pre-infected with a lentivirus encoding GFP from a constitutively active Ef1a promoter (Ef1aGFP). Tomato−/GFP+ cells contributed exclusively to embryonic tissues, whereas tomato+/GFP+ cells contributed to embryonic endoderm, ectoderm, mesoderm, the germ lineage as well as the yolk sac, and placenta (Fig. 4b–c, Supplementary Fig. 6a–b). The extraembryonic contribution of the tomato+/GFP+ cells included trophoblast giant cells of the placenta (Fig. 4c). Thus, the developmental potential of 2Ctomato+ cells includes embryonic plus extraembryonic tissues in contrast to the majority of ES cells in culture, which are restricted to generating only embryonic cell types.
We next examined whether Kdm1a mutant ES lines, that contain higher frequencies of 2Ctomato+ cells also had increased potency in mouse chimera assays. As expected, Kdm1a heterozygous ES cells contributed exclusively to embryonic tissues (in 5/5 chimeric embryos) but never to extraembryonic tissues (Fig. 4d). In contrast, Kdm1a homozygous mutant ES cells generated both embryonic (in 4/6 chimeric embryos) and extraembryonic tissues (in 5/6 chimeric embryos) (Fig. 4d). To confirm the increased potential of Kdm1a mutant ES cells, we used a competition chimera assay. We co-injected a 1:1 mixture of Kdm1a Fl/Fl and KO/KO ES cells into five wild type blastocysts. PCR was then used to detect the appearance of Kdm1a Fl/Fl or KO/KO cells in dissected tissues. We detected Kdm1a Fl/Fl ES cells in the embryonic tissues and amnion, but not the yolk sac or placenta (Fig. 4e). In contrast, Kdm1a mutant ES cells contributed to embryonic tissues, the amnion, yolk sac, and placental tissues, including trophoblast giant cells and PGCs (Fig. 4e–f). Thus, the artificial activation of 2C genes achieved by removing Kdm1a is associated with expanded fate potential.
We have shown that 2C:tomato+ cells within ES cultures have increased potency, but it is unclear whether entrance into this state is essential for their long-term pluripotency. To test this possibility, we performed serial depletion of 2C-like ES cells by genetic ablation with diptheria toxin (DTA). We generated ES lines by crossing 2CERT2CreERT2 mice (Supplementary Fig. 3a) with a Cre responsive DTA allele (ROSALSL-DTA) and treated the ES line with 4OHT for 20 passages (Supplementary Fig. 7a). We found that these 2C-depleted ES cultures were still capable of generating high contribution chimeras (Supplementary Fig. 7b), although their differentiation was biased toward mesoderm and ectoderm lineages in vitro (Supplementary Fig. 7c). These data suggest that occasional entry into the 2C-like state might help to preserve the broad embryonic fate potential of ES cells.
In mammalian development, the zygote and its daughter cells progress from totipotent cells capable of generating an entire mouse to more lineage restricted inner and outer cells of the morula capable of generating embryonic or extraembryonic lineages, respectively. A key transcriptional feature of the totipotent cells is the onset of zygote genome activation in which the embryo switches from a maternal to a zygotic transcriptome. To mark cells at this early stage of embryonic development, we generated a reporter with the regulatory elements from the endogenous retrovirus MERVL, which is highly restricted to the zygote/2C stage. Surprisingly, we found that rare ES and iPS cells expressed the reporter. When we characterized these cells, we found that they lacked expression of the pluripotency proteins Oct4, Sox2, and Nanog. Instead, these rare cells expressed a large number of genes restricted to the 2C stage, and most importantly, were capable of forming both embryonic and extraembryonic lineages (Fig. 5a–b). Our studies identify a rare 2C-like cell in ES cultures that has expanded fate potential.
Although it is unclear how MERVL and other 2C genes regulate potency, several lines of evidence indicate that 2C-like cells are required for the health and maintenance of ES cultures. First, we found that nearly all cells enter into the 2C-like state over increasing passage. Second, when we depleted 2C-like ES cells from cell cultures we found that their differentiation characteristics were altered to generate more ectoderm and mesoderm derivatives. Third, functional studies of the Zscan4 gene, found adjacent to a full-length MERVL element and highly enriched within 2Ctomato+ cells, have shown that it is required for the maintenance of telomeres within ES cultures1414. Another important question that remains is whether the selection of these special ES cells can be utilized for practical purposes, such as reprogramming somatic nuclei. This idea is supported by the finding that 2C genes are not properly activated in cloned embryos, and that reprogramming efficiency is enhanced by inhibition of histone deacetylases and Kdm1a, which repress the 2C-state32,33,34. Thus overexpression of one or multiple MERVL-linked 2C genes or inhibition of other 2C-gene repressors may be useful strategies to facilitate reprogramming. This possibility is supported by the recent finding that forced Zscan4 expression in fibroblasts enhances their iPS reprogramming efficiency35.
Transposable elements are a major driving force of evolution. Our findings support the notion that the co-option of retrotransposable elements by cellular genes can serve as an evolutionary mechanism for coordinately linking the expression of many genes15,29,36. Transposon sequences have recently been shown to play a critical role in rewiring gene regulatory networks in ES cells and in the endometrium that contributed to the evolution of pregnancy in mammals37,38. It has also been speculated that ERVs played a role in the evolution of the placenta by providing fusogenic envelope genes adapted for formation of the syncytiotrophoblasts39. We suggest that ERVL elements, that are found in all placental mammals40, may have played an equally important gene regulatory role in early mammalian development, by contributing to the specification of cell types and leading to formation of placental tissues.
2Ctomato was created by digesting the MERVL LTR-Gag clone #929 with MluI and HindIII, resulting in MERVL LTR 1-730, and was ligated into pcDNA3 hygro tdtomato with CMV promoter removed. To generate 2Ctomato ES cells, Kdm1a Fl/Fl; Cre:ERT ES cells were transfected with 2Ctomato using Lipofectamine 2000 (Invitrogen) and selected with 150ug/ml hygromycin for 7 days. Colonies containing tomato+ cells were then picked and expanded. 2CERT2-Cre-ERT2 was generated by replacing tdtomato with an ERT2-Cre-ERT2 insert using EcoR1 and Not1 sites. DNA was linearized with Mlu1 and AvrII sites before injection into embryos to generate transgenic mice. The resulting mice were mated with ROSALSL-tomato mice (JAX 007905), ROSALSL-DTA mice (JAX 010527), or ROSALSL-LacZ mice (Gift of Anderson lab), and ES lines were derived using standard procedures. Kdm1a GT/GT, Kdm1a Fl/Fl, KAP1 ES3 Cre, and G9A TT2 ES cells were described previously29–31. RNA-Seq from oocytes and 2C embryos was performed by lysing litters of embryos (5–10 embryos) in Prelude Direct lysis buffer (Nugen) and amplifying RNA using the Ovation RNA-Seq system (Nugen) before library construction using the Tru-Seq RNA sample prep kit (Illumina). Microarray, QRT-PCR, immunostaining, and chimeric mouse injections were performed as described29. All animal experiments we performed in accordance with the Salk Institute's IACUC guidelines.
For RNA-Seq analysis of early stage embryos, three independent litters of superovulated oocytes or naturally fertilized superovulated oocytes were collected and lysed directly in 2 ul of Prelude Direct Lysis buffer (Nugen). RNA was then subject to amplification using the Ovation RNA-Seq System (Nugen). Amplified cDNA was fragmented using Covaris, and single end (oocytes) or paired end (2C embryos) libraries were then constructed using the mRNA-Seq sample prep kit (Illumina) or Tru-Seq RNA library construction kit (Illumina) starting with end repair. Sequencing was performed on either Illumina Genome Analyzer (oocytes) or Hi-Seq (2C embryos). 72 bp single-end reads (oocytes) or 100 bp paired-end reads (2C embryos) were aligned to the mouse genome using Bowtie, allowing up to three mismatches per alignment and up to 20 alignments per read, filtering out any read aligning in more than 20 locations. To compare the Oocyte data to the 2C data, read lengths were cut down to 72 bp (from the 3'end). The Oocyte data had an average of 33 million alignments per sample while the 2C data had an average of 49 million alignments per sample. Read counts were quantified using a custom gene reference based on UCSC's Knowngene reference. At each gene locus, all isoforms belonging to a single gene were fused into one transcript containing all exons from each isoform. Counts aligning in multiple locations were counted as a fraction of their total number of alignments at each location. Differential expression testing was performed with DESeq41. Genes with adjusted p-values less than 0.05 were marked as significant. Chimeric transcripts were identified utilizing the spliced alignment data produced by Tophat. Tophat identifies exons based on alignment pileup and it follows by aligning previously unaligned reads across potential splice junctions. We split the junction information into two lists, left and right side of each junction, and compared to both the UCSC known gene database for mm9 and to the RepeatMasker database, also from UCSC's database. Only junctions that hit an exon of a known model on one end and a repeat element on the other were retained. GO analysis was performed using the David Bioinformatic Resource (http://david.abcc.ncifcrf.gov/)42.
For RNA-Seq of 2Ctomato+ and− cells, Kdm1a KO ES cells, Kap1 KO ES cells, and G9a KO ES cells, sample libraries were prepared from 500-5ug of total RNA using the mRNA-Seq sample prep kit (Illumina) or Tru-Seq RNA library construction kit. Library samples were amplified on flow cells using cluster generation kit (Illumina) and then sequenced using consecutive 36 cycle sequencing kit on the Genome Analyzer (Illumina) or 100bp paired end reads on the Hi-Seq (Illumina). Raw sequence data was then aligned to the mouse genome using the short read aligner Bowtie and the default setting (2 mismatches per 25 bp and up to 40 genomic alignments) (http://bowtie-bio.sourceforge.net/index.shtml). RPKM values were also determined by Bowtie. For repetitive sequences, we aligned sequencing reads to the Repbase database using Bowtie (http://www.girinst.org/repbase/index.html). Differential expression was determined using DESeq as described above. To compare gene expression in G9a, Kdm1a, and Kap1 mutant ES cells, we identified upregulated genes by combining previously identified upregulated genes(29–31) and our own DESeq analysis.
The derivation and culture of Kdm1a GT/GT, Fl/Fl, and Fl/Fl, Cre-ERT ES cells were described previously29. The 2Ctomato construct was created by digesting the MERVL LTR-Gag clone #9 in pGL3 basic with MluI and HindIII, resulting in MERVL LTR 1-730, and was ligated into pcDNA3 hygro tdtomato digested with MluI and HindIII (to remove CMV promoter). To generate 2Ctomato ES cells, Kdm1a Fl/Fl; Cre:ERT2 ES cells were transfected with 2Ctomato using Lipofectamine 2000 (Invitrogen) and selected with 150ug/ml hygromycin for 7 days. Colonies containing tomato positive cells were then picked and expanded. 2Ctomato ES cells were also derived from a transgenic mouse generated by pronuclear injection of the 2Ctomato ES line. 2CERT2-Cre-ERT2 was generated by replacing tdtomato with an ERT2-Cre-ERT2 insert using EcoR1 and Not1 sites. DNA was linearized with Mlu1 and AvrII sites before injection into embryos to generate transgenic mice. The resulting mice were mated with ROSALSL-tomato mice (JAX 007905), ROSALSL-DTA mice (JAX 010527), or ROSALSL-LacZ mice (Gift of Anderson lab) and ES lines were derived using standard procedures.
KAP1 ES3Cre and G9A TT2 ES cells were described previously30,31. To recombine and delete the Kdm1a, KAP1, and G9A floxed alleles, cells were treated with 1uM 4OHT for 24 hours. Cells were then harvested at a minimum of 48 hours later to allow for loss of residual protein. To activate the 2CERT2-Cre-ERT2 transgene, cells were maintained in 1uM 4OHT and fed daily. 2Ctomato positive and negative cells were counted using a FACScan and sorted using a FACSDiVA. For differentiation assays, ES cells were grown in suspension in the absence of Lif as described29.
ES cells and iPS cells were plated on gelatinized glass coverslips on PMEFs. Cells were fixed with 4% PFA for 10 minutes, followed by washing with PBS-T (0.05% tween). Cells were then blocked in PBS-T containing 3% BSA for ten minutes and stained with primary antibody for 1 hour at room temperature. Antibodies used: mouse anti KAP1 (Abcam), 1:1000; mouse anti OCT3/4, Santa Cruz sc-5279, rabbit anti MERVL-GAG, gift of Heidmann lab, 1:2000; rat anti E-Cadherin, Abcam ab11512, 1:500, rabbit anti Pan Acetylated histone H3, Upstate #06-599, 1:1000; rabbit anti Pan acetyl H4, Upstate #06-598, 1:1000; and rabbit anti H3 DiMeK4, clone AW30, Abcam, 1:1000. After washing 3 times for 10 minutes with PBST, cells were stained with secondary antibody (1:1000 anti mouse, rat, or rabbit IgG Alexa fluor 488, 555, or 647) for one hour at room temperature and washed again 3 times with PBS-T. Coverslips were stained with DAPI in PBS for 5 minutes before inverting onto slides in mounting medium. Cells were then imaged using either an Olympus FV1000 confocal microscope and 60× oil objective, or a Zeiss Axioskop 2 epifluorescence microscope and 40× objective. Quantification of histone stains were performed with Fluoview. Preimplantation embryos were stained as described with minor modifications. Embryos were fixed in 4% PFA for 30 minutes and permeabilized in 0.25% Triton for 20 minutes, prior to blocking in 10% FBS for 1 hour in 0.1% Triton-PBS. Primary antibodies were incubated overnight at 4 degrees C in blocking buffer. Subsequent washes and secondary antibody incubations were at room temperature in 0.1% Triton-PBS.
A MERVL probe was generated by PCR from mouse ES cDNA using the forward primer 5' ccatccctgtcattgctca 3' and reverse primer 5' ccttttccaccccttgatt 3' and cloned into the PCR2.1 TOPO vector. A DIG labeled probe was prepared using in vitro transcription with the T7 polymerase. ES samples were fixed in 4% PFA, digested for two minutes with proteinase K, washed with PBS, acetylated, and hybridized with denatured probe overnight at 68 degrees. After washing with 5× SSC and 0.2 × SSC, DIG labeled probe was visualized using an anti DIG antibody coupled to alkaline phosphatase.
Whole cell extracts were prepared by pelleting ES cells at 200 × g and resuspending in 1:5 volume of 1% NP40 lysis buffer containing 10mM Tris, 150mM NaCl, and 1× protease inhibitors. To solubilize histones, extracts were also sonicated using a bioruptor on the high setting for 10 minutes. 10–50ug of total protein in LDS sample buffer (Invitrogen) was then loaded onto a 4–12% NuPage gel (Invitrogen), electrophoresed at 200V for 60 minutes, and transferred to nitrocellulose membranes at 30V for 90 minutes. Membranes were blocked in PBS-T containing 5% nonfat dry milk. Primary antibodies were incubated overnight at 4 degrees. Antibodies utilized: rabbit anti GAPDH, Santa Cruz sc25-778, 1:1000, rabbit and MERVL-GAG, gift of Heidmann lab, 1:1000, anti Pan AcH3, Upstate #06-599, 1:1000; anti Pan AcH4, Upstate #06-598, 1:1000; anti H3 DiMeK4 clone AW30, Abcam, 1:500; anti H4, Novus ab10158, 1:1000; anti H3, Novus, NB 500–171, 1:500. After washing extensively with PBS-T, secondary antibodies (anti rabbit or mouse HRP conjugate, 1:10,000 dilution) were incubated for one hour at room temperature. After washing extensively with PBS-T and water, blots were developed using ECL plus detection system (Amersham).
ES cells were lysed in tail lysis buffer (0.1M Tris pH8.5, 5mM EDTA, 0.2% SDS, 0.2M NaCl,) containing proteinase K (Roche) for 1 hour at 55 degrees C, followed by treatment with DNase free RNase for 30 minutes at 37 degrees C. DNA was then sonicated briefly and purified using Qiagen PCR purification columns. Bisulfite conversion of genomic DNA was carried out using the Epitect Bisulfite Kit (Qiagen). Bisulfite converted DNA was then PCR amplified using Accuprime Taq polymerase (Invitrogen) followed by TOPO TA cloning (Invitrogen). At least 10 individual clones per primer pair were sequenced (Eton Bio). Primer sequences were described previously29.
For QRT-PCR Analysis, first strand cDNA was generated from up to 5ug total RNA using Superscript III (Invitrogen) and polydT or random hexamer priming. QPCR was performed using SYBR green master mix (Applied Biosystems) in 96 well dishes in triplicate and repeated with at least two biological replicates with similar results. Standard curves were generated for each primer pair (described previously29) and expression levels were plotted relative to Gapdh (in arbitrary units).
Total RNA was prepared from 2C::tomato+ and− cells using RNEasy kits (Qiagen). Labeling of 100ng of total RNA was performed using the Whole Transcript (WT) Sense Target Labeling Assay kit (Affymetrix) before hybridization to Genechip Mouse Gene 1.0 ST Arrays. Probeset normalization and summarization were prepared using Robust Multichip Analysis (RMA) in Expression Console (Affymetrix).
ES cells were injected into either E2.5 or E3.5 C57Bl/6J embryos and cultured in vitro or implanted into pseudopregnant females. For PCR assays dissected tissues were placed in lysis buffer (1% SDS, 150mM NaCl, 10mM TriS pH8.0, 1mM EDTA pH 8.0) containing proteinase K overnight at 55 degrees C. DNA was then isolated by phenol chloroform extraction and ethanol precipitation, followed by PCR analysis with primers designed to amplify the Betageo cassette or the wild type Kdm1a foxed allele. For embryo imaging, chimeric mice were harvested between E9.5 and E12.5 and fixed with 4% PFA for two hours, washed extensively in PBS overnight, incubated in 30% sucrose for 4 hours, and frozen on dry ice in OCT. Cryosections were then taken and stained with DAPI before imaging.
We would like to thank members of the Pfaff laboratory for discussion; Shane Andrews, David Chambers, Yelena Dayn, Tsung-Chang Sung, James Fitzpatrick, Matthew Joens, Yury Sigal, Daniel Gibbs, and Ling Ouyang for technical assistance, Yoichi Shinkai and David Gilbert for G9a mutant ES cells, and T. Heidmann for MERVL Gag antibodies. This research was supported by the National Institute of Neurological Disorders and Stroke (R37NS037116) and the Marshall Heritage Foundation. T.S.M and W.D.G. were supported by CIRM and S.L.P. is an investigator of the Howard Hughes Medical Institute.
Supplementary Information is linked to the online version of the paper at www.nature.com/nature.
Author Contributions T.S.M. designed and performed all experiments with assistance from W.D.G., S.D., D.B. and K.L. under the supervision of S.L.P.. D.T. generated KAP1 ES3Cre cells and H.R. and D.T provided mRNA-Seq data from these cells. A.F. and O.S. generated and provided iPS lines and lentivirus constructs. T.S.M., W.D.G. and S.L.P. wrote the manuscript.
Microarray and RNA-Seq files were submitted to NCBI's GEO database at accession GSE33923. Reprints and permissions information is available at www.nature.com/reprints. The author's declare no competing financial interests.