Decoding the dynamic DNA methylation and hydroxymethylation landscapes in endodermal lineage intermediates during pancreatic differentiation of hESC

Dynamic changes in DNA methylation and demethylation reprogram transcriptional outputs to instruct lineage specification during development. Here, we applied an integrative epigenomic approach to unveil DNA (hydroxy)methylation dynamics representing major endodermal lineage intermediates during pancreatic differentiation of human embryonic stem cells (hESCs). We found that 5-hydroxymethylcytosine (5hmC) marks genomic regions to be demethylated in the descendent lineage, thus reshaping the DNA methylation landscapes during pancreatic lineage progression. DNA hydroxymethylation is positively correlated with enhancer activities and chromatin accessibility, as well as the selective binding of lineage-specific pioneer transcription factors, during pancreatic differentiation. We further discovered enrichment of hydroxymethylated regions (termed ‘5hmC-rim’) at the boundaries of large hypomethylated functional genomic regions, including super-enhancer, DNA methylation canyon and broad-H3K4me3 peaks. We speculate that ‘5hmC-rim’ might safeguard low levels of cytosine methylation at these regions. Our comprehensive analysis highlights the importance of dynamic changes of epigenetic landscapes in driving pancreatic differentiation of hESC.

Related papers:

title: Functions of DNA methylation: islands, start sites, gene bodies and beyond
abstract: DNA methylation is frequently described as a 'silencing' epigenetic mark, and indeed this function of 5-methylcytosine was originally proposed in the 1970s. Now, thanks to improved genome-scale mapping of methylation, we can evaluate DNA methylation in different genomic contexts: transcriptional start sites with or without CpG islands, in gene bodies, at regulatory elements and at repeat sequences. The emerging picture is that the function of DNA methylation seems to vary with context, and the relationship between DNA methylation and transcription is more nuanced than we realized at first. Improving our understanding of the functions of DNA methylation is necessary for interpreting changes in this mark that are observed in diseases such as cancer.

title: DNA methylation patterns and epigenetic memory.
abstract: The character of a cell is defined by its constituent proteins, which are the result of specific patterns of gene expression. Crucial determinants of gene expression patterns are DNA-binding transcription factors that choose genes for transcriptional activation or repression by recognizing the sequence of DNA bases in their promoter regions. Interaction of these factors with their cognate sequences triggers a chain of events, often involving changes in the structure of chromatin, that leads to the assembly of an active transcription complex (e.g., Cosma et al. 1999). But the types of transcription factors present in a cell are not alone sufficient to define its spectrum of gene activity, as the transcriptional potential of a genome can become restricted in a stable manner during development. The constraints imposed by developmental history probably account for the very low efficiency of cloning animals from the nuclei of differentiated cells (Rideout et al. 2001; Wakayama and Yanagimachi 2001). A “transcription factors only” model would predict that the gene expression pattern of a differentiated nucleus would be completely reversible upon exposure to a new spectrum of factors. Although many aspects of expression can be reprogrammed in this way (Gurdon 1999), some marks of differentiation are evidently so stable that immersion in an alien cytoplasm cannot erase the memory. The genomic sequence of a differentiated cell is thought to be identical in most cases to that of the zygote from which it is descended (mammalian B and T cells being an obvious exception). This means that the marks of developmental history are unlikely to be caused by widespread somatic mutation. Processes less irrevocable than mutation fall under the umbrella term “epigenetic” mechanisms. A current definition of epigenetics is: “The study of mitotically and/or meiotically heritable changes in gene function that cannot be explained by changes in DNA sequence” (Russo et al. 1996). There are two epigenetic systems that affect animal development and fulfill the criterion of heritability: DNA methylation and the Polycomb-trithorax group (Pc-G/trx) protein complexes. (Histone modification has some attributes of an epigenetic process, but the issue of heritability has yet to be resolved.) This review concerns DNA methylation, focusing on the generation, inheritance, and biological significance of genomic methylation patterns in the development of mammals. Data will be discussed favoring the notion that DNA methylation may only affect genes that are already silenced by other mechanisms in the embryo. Embryonic transcription, on the other hand, may cause the exclusion of the DNA methylation machinery. The heritability of methylation states and the secondary nature of the decision to invite or exclude methylation support the idea that DNA methylation is adapted for a specific cellular memory function in development. Indeed, the possibility will be discussed that DNA methylation and Pc-G/trx may represent alternative systems of epigenetic memory that have been interchanged over evolutionary time. Animal DNA methylation has been the subject of several recent reviews (Bird and Wolffe 1999; Bestor 2000; Hsieh 2000; Costello and Plass 2001; ). For recent reviews of plant and fungal DNA methylation, see Finnegan et al. (2000), Martienssen and Colot (2001), and Matzke et al. (2001).

title: DNA methylation and cellular reprogramming.
abstract: The recent discovery that a small number of defined factors are sufficient to reprogram somatic cells into pluripotent stem cells has significantly expanded our knowledge of the plasticity of the epigenome. In this review we discuss some aspects of cell fate plasticity and epigenetic alterations, with emphasis on DNA methylation during cellular reprogramming. Recent data suggest that DNA methylation is a major barrier to induced pluripotent stem (iPS) cell reprogramming. The demethylating agent 5-azacytidine can enhance the efficiency of iPS cells generation and the putative DNA demethylase protein activation-induced cytidine deaminase (AID/AICDA) can erase DNA methylation at pluripotency gene promoters, thereby allowing cellular reprogramming. Elucidation of the epigenetic changes taking place during cellular reprogramming will enhance our understanding of stem cell biology and facilitate therapeutic applications.

title: DNA methylation: roles in mammalian development
abstract: DNA methylation is among the best studied epigenetic modifications and is essential to mammalian development. Although the methylation status of most CpG dinucleotides in the genome is stably propagated through mitosis, improvements to methods for measuring methylation have identified numerous regions in which it is dynamically regulated. In this Review, we discuss key concepts in the function of DNA methylation in mammals, stemming from more than two decades of research, including many recent studies that have elucidated when and where DNA methylation has a regulatory role in the genome. We include insights from early development, embryonic stem cells and adult lineages, particularly haematopoiesis, to highlight the general features of this modification as it participates in both global and localized epigenetic regulation.

title: Hotspots of aberrant epigenomic reprogramming in human induced pluripotent stem cells
abstract: Induced pluripotent stem cells (iPSCs) offer immense potential for regenerative medicine and studies of disease and development. Somatic cell reprogramming involves epigenomic reconfiguration, conferring iPSCs with characteristics similar to embryonic stem (ES) cells. However, it remains unknown how complete the reestablishment of ES-cell-like DNA methylation patterns is throughout the genome. Here we report the first whole-genome profiles of DNA methylation at single-base resolution in five human iPSC lines, along with methylomes of ES cells, somatic cells, and differentiated iPSCs and ES cells. iPSCs show significant reprogramming variability, including somatic memory and aberrant reprogramming of DNA methylation. iPSCs share megabase-scale differentially methylated regions proximal to centromeres and telomeres that display incomplete reprogramming of non-CG methylation, and differences in CG methylation and histone modifications. Lastly, differentiation of iPSCs into trophoblast cells revealed that errors in reprogramming CG methylation are transmitted at a high frequency, providing an iPSC reprogramming signature that is maintained after differentiation.

title: Two DNA methyltransferases from murine erythroleukemia cells: purification, sequence specificity, and mode of interaction with DNA.
abstract: Dye-ligand chromatography on Cibacron blue F3GA-agarose has been used to resolve two species of DNA (cytosine-5-)-methyltransferase from nuclear extracts of uninduced Friend murine erythroleukemia cells. Each species has been highly purified; the activities in the first and second peaks were associated with polypeptides of Mr 150,000 and 175,000, respectively. Analysis of substrate specificity with synthetic DNAs and restriction fragments of phi X174 replicative form DNA and pBR322 DNA showed that neither enzyme had dependence on the sequence context of CpG dinucleotides; poly(dG-dC) had the greatest methyl-accepting activity of any unmethylated DNA substrate tested. De novo methylation by both enzymes was inefficient relative to methylation of hemimethylated sites. Methyl-accepting activity was strongly dependent on DNA chain length. This observation suggests that binding to DNA, followed by one-dimensional diffusion of enzyme along the DNA molecule, is important in the mechanism by which DNA methyltransferase locates its recognition sites.

title: DNA methylation in mammals.
abstract: DNA methylation is one of the best characterized epigenetic modifications. In mammals it is involved in various biological processes including the silencing of transposable elements, regulation of gene expression, genomic imprinting, and X-chromosome inactivation. This article describes how DNA methylation serves as a cellular memory system and how it is dynamically regulated through the action of the DNA methyltransferase (DNMT) and ten eleven translocation (TET) enzymes. Its role in the regulation of gene expression, through its interplay with histone modifications, is also described, and its implication in human diseases discussed. The exciting areas of investigation that will likely become the focus of research in the coming years are outlined in the summary.

title: Conversion of 5-Methylcytosine to 5-Hydroxymethylcytosine in Mammalian DNA by MLL Partner TET1
abstract: Methylation Mediation Methylation of cytosine bases, 5-methylcytosine (5mC), in DNA plays an important regulatory role in mammalian genomes. Methylation patterns are often inherited across generations, but they can also be dynamic, suggesting that active DNA demethylation pathways exist. One such pathway, best characterized in plants, involves the removal of the 5mC base, and its replacement by C, via a DNA repair mechanism. Kriaucionis and Heintz (p. 929, published online 16 April) now show that, as well as 5mC in mammalian genomes, there are also significant amounts of 5-hydroxymethylcytosine (5hmC) in DNA of Purkinje neurons, which have large nuclei with apparently very little heterochromatin. Tahiliani et al. (p. 930, published online 16 April) find that the protein TET1 is capable of converting 5mC into 5hmC both in vitro and in vivo. 5-Hydroxymethylcytosine is also present in embryonic stem cells, and levels of 5hmC and TET1 show correlated variation during cell differentiation. Methylated C bases, an important epigenetic mark in genomic DNA, can be enzymically converted to 5-hydroxymethylcytosine. DNA cytosine methylation is crucial for retrotransposon silencing and mammalian development. In a computational search for enzymes that could modify 5-methylcytosine (5mC), we identified TET proteins as mammalian homologs of the trypanosome proteins JBP1 and JBP2, which have been proposed to oxidize the 5-methyl group of thymine. We show here that TET1, a fusion partner of the MLL gene in acute myeloid leukemia, is a 2-oxoglutarate (2OG)- and Fe(II)-dependent enzyme that catalyzes conversion of 5mC to 5-hydroxymethylcytosine (hmC) in cultured cells and in vitro. hmC is present in the genome of mouse embryonic stem cells, and hmC levels decrease upon RNA interference–mediated depletion of TET1. Thus, TET proteins have potential roles in epigenetic regulation through modification of 5mC to hmC.

title: Tet1 and Tet2 regulate 5-hydroxymethylcytosine production and cell lineage specification in mouse embryonic stem cells.
abstract: TET family enzymes convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) in DNA. Here, we show that Tet1 and Tet2 are Oct4-regulated enzymes that together sustain 5hmC in mouse embryonic stem cells (ESCs) and are induced concomitantly with 5hmC during reprogramming of fibroblasts to induced pluripotent stem cells. ESCs depleted of Tet1 by RNAi show diminished expression of the Nodal antagonist Lefty1 and display hyperactive Nodal signaling and skewed differentiation into the endoderm-mesoderm lineage in embryoid bodies in vitro. In Fgf4- and heparin-supplemented culture conditions, Tet1-depleted ESCs activate the trophoblast stem cell lineage determinant Elf5 and can colonize the placenta in midgestation embryo chimeras. Consistent with these findings, Tet1-depleted ESCs form aggressive hemorrhagic teratomas with increased endoderm, reduced neuroectoderm, and ectopic appearance of trophoblastic giant cells. Thus, 5hmC is an epigenetic modification associated with the pluripotent state, and Tet1 functions to regulate the lineage differentiation potential of ESCs.

title: Tet-Mediated Formation of 5-Carboxylcytosine and Its Excision by TDG in Mammalian DNA
abstract: Evidence for a possible route for DNA demethylation in animals is suggested. The prevalent DNA modification in higher organisms is the methylation of cytosine to 5-methylcytosine (5mC), which is partially converted to 5-hydroxymethylcytosine (5hmC) by the Tet (ten eleven translocation) family of dioxygenases. Despite their importance in epigenetic regulation, it is unclear how these cytosine modifications are reversed. Here, we demonstrate that 5mC and 5hmC in DNA are oxidized to 5-carboxylcytosine (5caC) by Tet dioxygenases in vitro and in cultured cells. 5caC is specifically recognized and excised by thymine-DNA glycosylase (TDG). Depletion of TDG in mouse embyronic stem cells leads to accumulation of 5caC to a readily detectable level. These data suggest that oxidation of 5mC by Tet proteins followed by TDG-mediated base excision of 5caC constitutes a pathway for active DNA demethylation.

title: Tet Proteins Can Convert 5-Methylcytosine to 5-Formylcytosine and 5-Carboxylcytosine
abstract: Evidence for a possible route for DNA demethylation in animals is suggested. 5-methylcytosine (5mC) in DNA plays an important role in gene expression, genomic imprinting, and suppression of transposable elements. 5mC can be converted to 5-hydroxymethylcytosine (5hmC) by the Tet (ten eleven translocation) proteins. Here, we show that, in addition to 5hmC, the Tet proteins can generate 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC) from 5mC in an enzymatic activity–dependent manner. Furthermore, we reveal the presence of 5fC and 5caC in genomic DNA of mouse embryonic stem cells and mouse organs. The genomic content of 5hmC, 5fC, and 5caC can be increased or reduced through overexpression or depletion of Tet proteins. Thus, we identify two previously unknown cytosine derivatives in genomic DNA as the products of Tet proteins. Our study raises the possibility that DNA demethylation may occur through Tet-catalyzed oxidation followed by decarboxylation.

title: Recognition and potential mechanisms for replication and erasure of cytosine hydroxymethylation
abstract: Cytosine residues in mammalian DNA occur in at least three forms, cytosine (C), 5-methylcytosine (M; 5mC) and 5-hydroxymethylcytosine (H; 5hmC). During semi-conservative DNA replication, hemi-methylated (M/C) and hemi-hydroxymethylated (H/C) CpG dinucleotides are transiently generated, where only the parental strand is modified and the daughter strand contains native cytosine. Here, we explore the role of DNA methyltransferases (DNMT) and ten eleven translocation (Tet) proteins in perpetuating these states after replication, and the molecular basis of their recognition by methyl-CpG-binding domain (MBD) proteins. Using recombinant proteins and modified double-stranded deoxyoligonucleotides, we show that DNMT1 prefers a hemi-methylated (M/C) substrate (by a factor of >60) over hemi-hydroxymethylated (H/C) and unmodified (C/C) sites, whereas both DNMT3A and DNMT3B have approximately equal activity on all three substrates (C/C, M/C and H/C). Binding of MBD proteins to methylated DNA inhibited Tet1 activity, suggesting that MBD binding may also play a role in regulating the levels of 5hmC. All five MBD proteins generally have reduced binding affinity for 5hmC relative to 5mC in the fully modified context (H/M versus M/M), though their relative abilities to distinguish the two varied considerably. We further show that the deamination product of 5hmC could be excised by thymine DNA glycosylase and MBD4 glycosylases regardless of context.

title: Replication-Dependent Loss of 5-Hydroxymethylcytosine in Mouse Preimplantation Embryos
abstract: The oxidation product of methylated cytosine is passively lost from DNA in the zygote as cell division progresses. Although global erasure of DNA methylation has been observed in zygotes and primordial germ cells, the responsible enzyme(s) have been elusive. The demonstration that members of the Tet (ten eleven translocation) family of proteins are capable of catalyzing conversion of 5-methylcytosine (5mC) of DNA to 5-hydroxymethylcytosine (5hmC) raises the possibility that Tet proteins may participate in this process. Indeed, recent studies have implicated the involvement of Tet3 in the conversion of 5mC to 5hmC in zygotes. This result, combined with the demonstration that Tet proteins can further oxidize 5hmC to 5-carboxylcytosine followed by excision by thymine-DNA glycosylase, raises the possibility that active demethylation may take place in a process that involves Tet3-mediated oxidation followed by base excision repair. We demonstrated by immunostaining of mitotic chromosome spreads of preimplantation embryos that the 5hmC associated with the paternal genome in zygotes is gradually lost during preimplantation development. Our study suggests that, although the conversion of 5mC to 5hmC in zygotes is an enzyme-catalyzed process, loss of 5hmC during preimplantation appears to be a DNA replication–dependent passive process.

title: Effects of Tet-induced oxidation products of 5-methylcytosine on Dnmt1- and DNMT3a-mediated cytosine methylation.
abstract: We investigated systematically the effects of Tet-induced oxidation products of 5-methylcytosine on Dnmt1- and DNMT3a-mediated cytosine methylation in synthetic duplex DNA. We found that the replacement of 5-methylcytosine at a CpG site with a 5-hydroxymethylcytosine, 5-formylcytosine, 5-carboxylcytosine or 5-hydroxymethyluracil resulted in altered methylation of cytosine at both the opposite and the neighboring CpG sites. Our results provided important new knowledge about the implications of the 5-methylcytosine oxidation products in maintenance cytosine methylation.

title: Active DNA demethylation at enhancers during the vertebrate phylotypic period
abstract: The vertebrate body plan and organs are shaped during a conserved embryonic phase called the phylotypic stage. However, the mechanisms that guide the epigenome through this transition and their evolutionary conservation remain elusive. Here we report widespread DNA demethylation of enhancers during the phylotypic period in zebrafish, Xenopus tropicalis and mouse. These enhancers are linked to developmental genes that display coordinated transcriptional and epigenomic changes in the diverse vertebrates during embryogenesis. Binding of Tet proteins to (hydroxy)methylated DNA and enrichment of 5-hydroxymethylcytosine in these regions implicated active DNA demethylation in this process. Furthermore, loss of function of Tet1, Tet2 and Tet3 in zebrafish reduced chromatin accessibility and increased methylation levels specifically at these enhancers, indicative of DNA methylation being an upstream regulator of phylotypic enhancer function. Overall, our study highlights a regulatory module associated with the most conserved phase of vertebrate embryogenesis and suggests an ancient developmental role for Tet dioxygenases.

title: Epigenetic memory at embryonic enhancers identified in DNA methylation maps from adult mouse tissues
abstract: Mammalian development requires cytosine methylation, a heritable epigenetic mark of cellular memory believed to maintain a cell's unique gene expression pattern. However, it remains unclear how dynamic DNA methylation relates to cell type–specific gene expression and animal development. Here, by mapping base-resolution methylomes in 17 adult mouse tissues at shallow coverage, we identify 302,864 tissue-specific differentially methylated regions (tsDMRs) and estimate that >6.7% of the mouse genome is variably methylated. Supporting a prominent role for DNA methylation in gene regulation, most tsDMRs occur at distal cis-regulatory elements. Unexpectedly, some tsDMRs mark enhancers that are dormant in adult tissues but active in embryonic development. These 'vestigial' enhancers are hypomethylated and lack active histone modifications in adult tissues but nevertheless exhibit activity during embryonic development. Our results provide new insights into the role of DNA methylation at tissue-specific enhancers and suggest that epigenetic memory of embryonic development may be retained in adult tissues.

title: TET-mediated DNA demethylation controls gastrulation by regulating Lefty–Nodal signalling
abstract: Mammalian genomes undergo epigenetic modifications, including cytosine methylation by DNA methyltransferases (DNMTs). Oxidation of 5-methylcytosine by the Ten-eleven translocation (TET) family of dioxygenases can lead to demethylation. Although cytosine methylation has key roles in several processes such as genomic imprinting and X-chromosome inactivation, the functional significance of cytosine methylation and demethylation in mouse embryogenesis remains to be fully determined. Here we show that inactivation of all three Tet genes in mice leads to gastrulation phenotypes, including primitive streak patterning defects in association with impaired maturation of axial mesoderm and failed specification of paraxial mesoderm, mimicking phenotypes in embryos with gain-of-function Nodal signalling. Introduction of a single mutant allele of Nodal in the Tet mutant background partially restored patterning, suggesting that hyperactive Nodal signalling contributes to the gastrulation failure of Tet mutants. Increased Nodal signalling is probably due to diminished expression of the Lefty1 and Lefty2 genes, which encode inhibitors of Nodal signalling. Moreover, reduction in Lefty gene expression is linked to elevated DNA methylation, as both Lefty–Nodal signalling and normal morphogenesis are largely restored in Tet-deficient embryos when the Dnmt3a and Dnmt3b genes are disrupted. Additionally, a point mutation in Tet that specifically abolishes the dioxygenase activity causes similar morphological and molecular abnormalities as the null mutation. Taken together, our results show that TET-mediated oxidation of 5-methylcytosine modulates Lefty–Nodal signalling by promoting demethylation in opposition to methylation by DNMT3A and DNMT3B. These findings reveal a fundamental epigenetic mechanism featuring dynamic DNA methylation and demethylation crucial to regulation of key signalling pathways in early body plan formation.

title: An epigenomic road map for endoderm development.
abstract: While studies of organ development have traditionally relied on model organisms, recent advances in embryonic stem cell (ESC) culture allow investigation of organogenesis in human cells. Wang et al. (2015) employ this system to map the dynamic enhancer landscape during ESC differentiation to the endoderm derivatives pancreas and liver.

title: Enrichment of human embryonic stem cell‐derived NKX6.1‐expressing pancreatic progenitor cells accelerates the maturation of insulin‐secreting cells in vivo
abstract: Human embryonic stem cells (hESCs) are considered a potential alternative to cadaveric islets as a source of transplantable cells for treating patients with diabetes. We previously described a differentiation protocol to generate pancreatic progenitor cells from hESCs, composed of mainly pancreatic endoderm (PDX1/NKX6.1‐positive), endocrine precursors (NKX2.2/synaptophysin‐positive, hormone/NKX6.1‐negative), and polyhormonal cells (insulin/glucagon‐positive, NKX6.1‐negative). However, the relative contributions of NKX6.1‐negative versus NKX6.1‐positive cell fractions to the maturation of functional β‐cells remained unclear. To address this question, we generated two distinct pancreatic progenitor cell populations using modified differentiation protocols. Prior to transplant, both populations contained a high proportion of PDX1‐expressing cells (∼85%–90%) but were distinguished by their relatively high (∼80%) or low (∼25%) expression of NKX6.1. NKX6.1‐high and NKX6.1‐low progenitor populations were transplanted subcutaneously within macroencapsulation devices into diabetic mice. Mice transplanted with NKX6.1‐low cells remained hyperglycemic throughout the 5‐month post‐transplant period whereas diabetes was reversed in NKX6.1‐high recipients within 3 months. Fasting human C‐peptide levels were similar between groups throughout the study, but only NKX6.1‐high grafts displayed robust meal‐, glucose‐ and arginine‐responsive insulin secretion as early as 3 months post‐transplant. NKX6.1‐low recipients displayed elevated fasting glucagon levels. Theracyte devices from both groups contained almost exclusively pancreatic endocrine tissue, but NKX6.1‐high grafts contained a greater proportion of insulin‐positive and somatostatin‐positive cells, whereas NKX6.1‐low grafts contained mainly glucagon‐expressing cells. Insulin‐positive cells in NKX6.1‐high, but not NKX6.1‐low grafts expressed nuclear MAFA. Collectively, this study demonstrates that a pancreatic endoderm‐enriched population can mature into highly functional β‐cells with only a minor contribution from the endocrine subpopulation. Stem Cells 2013;31:2432–2442

title: Dynamic chromatin remodeling mediated by polycomb proteins orchestrates pancreatic differentiation of human embryonic stem cells.
abstract: Embryonic development is characterized by dynamic changes in gene expression, yet the role of chromatin remodeling in these cellular transitions remains elusive. To address this question, we profiled the transcriptome and select chromatin modifications at defined stages during pancreatic endocrine differentiation of human embryonic stem cells. We identify removal of Polycomb group (PcG)-mediated repression on stage-specific genes as a key mechanism for the induction of developmental regulators. Furthermore, we discover that silencing of transitory genes during lineage progression associates with reinstatement of PcG-dependent repression. Significantly, in vivo- but not in vitro-differentiated endocrine cells exhibit close similarity to primary human islets in regard to transcriptome and chromatin structure. We further demonstrate that endocrine cells produced in vitro do not fully eliminate PcG-mediated repression on endocrine-specific genes, probably contributing to their malfunction. These studies reveal dynamic chromatin remodeling during developmental lineage progression and identify possible strategies for improving cell differentiation in culture.

title: Epigenetic priming of enhancers predicts developmental competence of hESC-derived endodermal lineage intermediates.
abstract: Embryonic development relies on the capacity of progenitor cells to appropriately respond to inductive cues, a cellular property known as developmental competence. Here, we report that epigenetic priming of enhancers signifies developmental competence during endodermal lineage diversification. Chromatin mapping during pancreatic and hepatic differentiation of human embryonic stem cells revealed the en masse acquisition of a poised chromatin state at enhancers specific to endoderm-derived cell lineages in gut tube intermediates. Experimentally, the acquisition of this poised enhancer state predicts the ability of endodermal intermediates to respond to inductive signals. Furthermore, these enhancers are first recognized by the pioneer transcription factors FOXA1 and FOXA2 when competence is acquired, while subsequent recruitment of lineage-inductive transcription factors, such as PDX1, leads to enhancer and target gene activation. Together, our results identify the acquisition of a poised chromatin state at enhancers as a mechanism by which progenitor cells acquire developmental competence.

title: Human DNA methylomes at base resolution show widespread epigenomic differences
abstract: DNA cytosine methylation is a central epigenetic modification that has essential roles in cellular processes including genome regulation, development and disease. Here we present the first genome-wide, single-base-resolution maps of methylated cytosines in a mammalian genome, from both human embryonic stem cells and fetal fibroblasts, along with comparative analysis of messenger RNA and small RNA components of the transcriptome, several histone modifications, and sites of DNA–protein interaction for several key regulatory factors. Widespread differences were identified in the composition and patterning of cytosine methylation between the two genomes. Nearly one-quarter of all methylation identified in embryonic stem cells was in a non-CG context, suggesting that embryonic stem cells may use different methylation mechanisms to affect gene regulation. Methylation in non-CG contexts showed enrichment in gene bodies and depletion in protein binding sites and enhancers. Non-CG methylation disappeared upon induced differentiation of the embryonic stem cells, and was restored in induced pluripotent stem cells. We identified hundreds of differentially methylated regions proximal to genes involved in pluripotency and differentiation, and widespread reduced methylation levels in fibroblasts associated with lower transcriptional activity. These reference epigenomes provide a foundation for future studies exploring this key epigenetic modification in human disease and development.

title: Genome-wide mapping of 5-hydroxymethylcytosine in embryonic stem cells
abstract: 5-hydroxymethylcytosine (5hmC) is a modified base present at low levels in diverse cell types in mammals. 5hmC is generated by the TET family of Fe(II) and 2-oxoglutarate-dependent enzymes through oxidation of 5-methylcytosine (5mC). 5hmC and TET proteins have been implicated in stem cell biology and cancer, but information on the genome-wide distribution of 5hmC is limited. Here we describe two novel and specific approaches to profile the genomic localization of 5hmC. The first approach, termed GLIB (glucosylation, periodate oxidation, biotinylation) uses a combination of enzymatic and chemical steps to isolate DNA fragments containing as few as a single 5hmC. The second approach involves conversion of 5hmC to cytosine 5-methylenesulphonate (CMS) by treatment of genomic DNA with sodium bisulphite, followed by immunoprecipitation of CMS-containing DNA with a specific antiserum to CMS. High-throughput sequencing of 5hmC-containing DNA from mouse embryonic stem (ES) cells showed strong enrichment within exons and near transcriptional start sites. 5hmC was especially enriched at the start sites of genes whose promoters bear dual histone 3 lysine 27 trimethylation (H3K27me3) and histone 3 lysine 4 trimethylation (H3K4me3) marks. Our results indicate that 5hmC has a probable role in transcriptional regulation, and suggest a model in which 5hmC contributes to the ‘poised’ chromatin signature found at developmentally-regulated genes in ES cells.

title: The anti-CMS technique for genome-wide mapping of 5-hydroxymethylcytosine
abstract: 5-hydroxymethylcytosine (5hmC) is a recently discovered base in the mammalian genome, produced upon oxidation of 5-methylcytosine (5mC) in a process catalyzed by TET proteins. The biological functions of 5hmC and further oxidation products of 5mC are under intense investigation, as they are likely intermediates in DNA demethylation pathways. Here we describe a novel protocol to profile 5hmC at a genome-wide scale. This approach is based on sodium bisulfite–mediated conversion of 5hmC to cytosine-5-methylenesulfonate (CMS); CMS-containing DNA fragments are then immunoprecipitated using a CMS-specific antiserum. The anti-CMS technique is highly specific with a low background, and is much less dependent on 5hmC density than anti-5hmC immunoprecipitation (IP). Moreover, it does not enrich for CA and CT repeats, as noted for 5hmC DNA IP using antibodies to 5hmC. The anti-CMS protocol takes 3 d to complete.

title: Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position
abstract: We describe an assay for transposase-accessible chromatin using sequencing (ATAC-seq), based on direct in vitro transposition of sequencing adaptors into native chromatin, as a rapid and sensitive method for integrative epigenomic analysis. ATAC-seq captures open chromatin sites using a simple two-step protocol with 500–50,000 cells and reveals the interplay between genomic locations of open chromatin, DNA-binding proteins, individual nucleosomes and chromatin compaction at nucleotide resolution. We discovered classes of DNA-binding factors that strictly avoided, could tolerate or tended to overlap with nucleosomes. Using ATAC-seq maps of human CD4+ T cells from a proband obtained on consecutive days, we demonstrated the feasibility of analyzing an individual's epigenome on a timescale compatible with clinical decision-making.

title: Model-based Analysis of ChIP-Seq (MACS)
abstract: We present Model-based Analysis of ChIP-Seq data, MACS, which analyzes data generated by short read sequencers such as Solexa's Genome Analyzer. MACS empirically models the shift size of ChIP-Seq tags, and uses it to improve the spatial resolution of predicted binding sites. MACS also uses a dynamic Poisson distribution to effectively capture local biases in the genome, allowing for more robust predictions. MACS compares favorably to existing ChIP-Seq peak-finding algorithms, and is freely available.

title: The Behaviour of 5-Hydroxymethylcytosine in Bisulfite Sequencing
abstract: We recently showed that enzymes of the TET family convert 5-mC to 5-hydroxymethylcytosine (5-hmC) in DNA. 5-hmC is present at high levels in embryonic stem cells and Purkinje neurons. The methylation status of cytosines is typically assessed by reaction with sodium bisulfite followed by PCR amplification. Reaction with sodium bisulfite promotes cytosine deamination, whereas 5-methylcytosine (5-mC) reacts poorly with bisulfite and is resistant to deamination. Since 5-hmC reacts with bisulfite to yield cytosine 5-methylenesulfonate (CMS), we asked how DNA containing 5-hmC behaves in bisulfite sequencing.  We used synthetic oligonucleotides with different distributions of cytosine as templates for generation of DNAs containing C, 5-mC and 5-hmC. The resulting DNAs were subjected in parallel to bisulfite treatment, followed by exposure to conditions promoting cytosine deamination. The extent of conversion of 5-hmC to CMS was estimated to be 99.7%. Sequencing of PCR products showed that neither 5-mC nor 5-hmC undergo C-to-T transitions after bisulfite treatment, confirming that these two modified cytosine species are indistinguishable by the bisulfite technique. DNA in which CMS constituted a large fraction of all bases (28/201) was much less efficiently amplified than DNA in which those bases were 5-mC or uracil (the latter produced by cytosine deamination). Using a series of primer extension experiments, we traced the inefficient amplification of CMS-containing DNA to stalling of Taq polymerase at sites of CMS modification, especially when two CMS bases were either adjacent to one another or separated by 1–2 nucleotides.  We have confirmed that the widely used bisulfite sequencing technique does not distinguish between 5-mC and 5-hmC. Moreover, we show that CMS, the product of bisulfite conversion of 5-hmC, tends to stall DNA polymerases during PCR, suggesting that densely hydroxymethylated regions of DNA may be underrepresented in quantitative methylation analyses.

title: MOABS: model based analysis of bisulfite sequencing data
abstract: Bisulfite sequencing (BS-seq) is the gold standard for studying genome-wide DNA methylation. We developed MOABS to increase the speed, accuracy, statistical power and biological relevance of BS-seq data analysis. MOABS detects differential methylation with 10-fold coverage at single-CpG resolution based on a Beta-Binomial hierarchical model and is capable of processing two billion reads in 24 CPU hours. Here, using simulated and real BS-seq data, we demonstrate that MOABS outperforms other leading algorithms, such as Fisher’s exact test and BSmooth. Furthermore, MOABS analysis can be easily extended to differential 5hmC analysis using RRBS and oxBS-seq. MOABS is available at http://code.google.com/p/moabs/.

title: Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities.
abstract: Genome-scale studies have revealed extensive, cell type-specific colocalization of transcription factors, but the mechanisms underlying this phenomenon remain poorly understood. Here, we demonstrate in macrophages and B cells that collaborative interactions of the common factor PU.1 with small sets of macrophage- or B cell lineage-determining transcription factors establish cell-specific binding sites that are associated with the majority of promoter-distal H3K4me1-marked genomic regions. PU.1 binding initiates nucleosome remodeling, followed by H3K4 monomethylation at large numbers of genomic regions associated with both broadly and specifically expressed genes. These locations serve as beacons for additional factors, exemplified by liver X receptors, which drive both cell-specific gene expression and signal-dependent responses. Together with analyses of transcription factor binding and H3K4me1 patterns in other cell types, these studies suggest that simple combinations of lineage-determining transcription factors can specify the genomic sites ultimately responsible for both cell identity and cell type-specific responses to diverse signaling inputs.

title: Dissecting the dynamic changes of 5-hydroxymethylcytosine in T-cell development and differentiation
abstract: Hydroxymethylcytosine (5hmC) is an epigenetic DNA modification produced through the enzymatic activity of TET proteins. Here we present the first genome-wide mapping of 5hmC in T cells during sequential steps of lineage commitment in the thymus and the periphery (thymic DP, CD4 SP, and CD8 SP cells; peripheral naive CD8 and CD4 T cells; and in vitro-differentiated Th1 and Th2 cells). We show that 5hmC is enriched at gene bodies and cell type-specific enhancers, that its levels in the gene body correlate strongly with gene expression and histone modifications, and that its levels change dynamically during the course of T-cell development and differentiation. Our analysis will facilitate increased understanding of the role of 5hmC in T-cell development and differentiation. The discovery of Ten Eleven Translocation proteins, enzymes that oxidize 5-methylcytosine (5mC) in DNA, has revealed novel mechanisms for the regulation of DNA methylation. We have mapped 5-hydroxymethylcytosine (5hmC) at different stages of T-cell development in the thymus and T-cell differentiation in the periphery. We show that 5hmC is enriched in the gene body of highly expressed genes at all developmental stages and that its presence correlates positively with gene expression. Further emphasizing the connection with gene expression, we find that 5hmC is enriched in active thymus-specific enhancers and that genes encoding key transcriptional regulators display high intragenic 5hmC levels in precursor cells at those developmental stages where they exert a positive effect. Our data constitute a valuable resource that will facilitate detailed analysis of the role of 5hmC in T-cell development and differentiation.

title: Tet2 and Tet3 cooperate with B-lineage transcription factors to regulate DNA modification and chromatin accessibility
abstract: Ten-eleven translocation (TET) enzymes oxidize 5-methylcytosine, facilitating DNA demethylation and generating new epigenetic marks. Here we show that concomitant loss of Tet2 and Tet3 in mice at early B cell stage blocked the pro- to pre-B cell transition in the bone marrow, decreased Irf4 expression and impaired the germline transcription and rearrangement of the Igκ locus. Tet2/3-deficient pro-B cells showed increased CpG methylation at the Igκ 3’ and distal enhancers that was mimicked by depletion of E2A or PU.1, as well as a global decrease in chromatin accessibility at enhancers. Importantly, re-expression of the Tet2 catalytic domain in Tet2/3-deficient B cells resulted in demethylation of the Igκ enhancers and restored their chromatin accessibility. Our data suggest that TET proteins and lineage-specific transcription factors cooperate to influence chromatin accessibility and Igκ enhancer function by modulating the modification status of DNA.

title: DNMT3A and TET2 compete and cooperate to repress lineage-specific transcription factors in hematopoietic stem cells
abstract: Mutations in the epigenetic modifiers DNMT3A and TET2 non-randomly co-occur in lymphoma and leukemia despite their epistasis in the methylation–hydroxymethylation pathway. Using Dnmt3a and Tet2 double-knockout mice in which the development of malignancy is accelerated, we show that the double-knockout methylome reflects regions of independent, competitive and cooperative activity. Expression of lineage-specific transcription factors, including the erythroid regulators Klf1 and Epor, is upregulated in double-knockout hematopoietic stem cells (HSCs). DNMT3A and TET2 both repress Klf1, suggesting a model of cooperative inhibition by epigenetic modifiers. These data demonstrate a dual role for TET2 in promoting and inhibiting HSC differentiation, the loss of which, along with DNMT3A, obstructs differentiation, leading to transformation.

title: deepTools2: a next generation web server for deep-sequencing data analysis
abstract: We present an update to our Galaxy-based web server for processing and visualizing deeply sequenced data. Its core tool set, deepTools, allows users to perform complete bioinformatic workflows ranging from quality controls and normalizations of aligned reads to integrative analyses, including clustering and visualization approaches. Since we first described our deepTools Galaxy server in 2014, we have implemented new solutions for many requests from the community and our users. Here, we introduce significant enhancements and new tools to further improve data visualization and interpretation. deepTools continue to be open to all users and freely available as a web service at deeptools.ie-freiburg.mpg.de. The new deepTools2 suite can be easily deployed within any Galaxy framework via the toolshed repository, and we also provide source code for command line usage under Linux and Mac OS X. A public and documented API for access to deepTools functionality is also available.

title: An Integrated Encyclopedia of DNA Elements in the Human Genome
abstract: The human genome encodes the blueprint of life, but the function of the vast majority of its nearly three billion bases is unknown. The Encyclopedia of DNA Elements (ENCODE) project has systematically mapped regions of transcription, transcription factor association, chromatin structure and histone modification. These data enabled us to assign biochemical functions for 80% of the genome, in particular outside of the well-studied protein-coding regions. Many discovered candidate regulatory elements are physically associated with one another and with expressed genes, providing new insights into the mechanisms of gene regulation. The newly identified elements also show a statistical correspondence to sequence variants linked to human disease, and can thereby guide interpretation of this variation. Overall, the project provides new insights into the organization and regulation of our genes and genome, and is an expansive resource of functional annotations for biomedical research.

title: Dynamic regulation of 5-hydroxymethylcytosine in mouse ES cells and during differentiation
abstract: Methylation at the 5′ position of cytosine in DNA has important roles in genome function and is dynamically reprogrammed during early embryonic and germ cell development. The mammalian genome also contains 5-hydroxymethylcytosine (5hmC), which seems to be generated by oxidation of 5-methylcytosine (5mC) by the TET family of enzymes that are highly expressed in embryonic stem (ES) cells. Here we use antibodies against 5hmC and 5mC together with high throughput sequencing to determine genome-wide patterns of methylation and hydroxymethylation in mouse wild-type and mutant ES cells and differentiating embryoid bodies. We find that 5hmC is mostly associated with euchromatin and that whereas 5mC is under-represented at gene promoters and CpG islands, 5hmC is enriched and is associated with increased transcriptional levels. Most, if not all, 5hmC in the genome depends on pre-existing 5mC and the balance between these two modifications is different between genomic regions. Knockdown of Tet1 and Tet2 causes downregulation of a group of genes that includes pluripotency-related genes (including Esrrb, Prdm14, Dppa3, Klf2, Tcl1 and Zfp42) and a concomitant increase in methylation of their promoters, together with an increased propensity of ES cells for extraembryonic lineage differentiation. Declining levels of TETs during differentiation are associated with decreased hydroxymethylation levels at the promoters of ES cell-specific genes together with increased methylation and gene silencing. We propose that the balance between hydroxymethylation and methylation in the genome is inextricably linked with the balance between pluripotency and lineage commitment.

title: TET1 and hydroxymethylcytosine in transcription and DNA methylation fidelity
abstract: Enzymes catalysing the methylation of the 5-position of cytosine (mC) have essential roles in regulating gene expression and maintaining cellular identity. Recently, TET1 was found to hydroxylate the methyl group of mC, converting it to 5-hydroxymethyl cytosine (hmC). Here we show that TET1 binds throughout the genome of embryonic stem cells, with the majority of binding sites located at transcription start sites (TSSs) of CpG-rich promoters and within genes. The hmC modification is found in gene bodies and in contrast to mC is also enriched at CpG-rich TSSs. We provide evidence further that TET1 has a role in transcriptional repression. TET1 binds a significant proportion of Polycomb group target genes. Furthermore, TET1 associates and colocalizes with the SIN3A co-repressor complex. We propose that TET1 fine-tunes transcription, opposes aberrant DNA methylation at CpG-rich sequences and thereby contributes to the regulation of DNA methylation fidelity.

title: Selective chemical labeling reveals the genome-wide distribution of 5-hydroxymethylcytosine
abstract: In contrast to 5-methylcytosine (5-mC), which has been studied extensively, little is known about 5-hydroxymethylcytosine (5-hmC), a recently identified epigenetic modification present in substantial amounts in certain mammalian cell types. Here we present a method for determining the genome-wide distribution of 5-hmC. We use the T4 bacteriophage β-glucosyltransferase to transfer an engineered glucose moiety containing an azide group onto the hydroxyl group of 5-hmC. The azide group can be chemically modified with biotin for detection, affinity enrichment and sequencing of 5-hmC–containing DNA fragments in mammalian genomes. Using this method, we demonstrate that 5-hmC is present in human cell lines beyond those previously recognized. We also find a gene expression level–dependent enrichment of intragenic 5-hmC in mouse cerebellum and an age-dependent acquisition of this modification in specific gene bodies linked to neurodegenerative disorders.

title: Impaired hydroxylation of 5-methylcytosine in myeloid cancers with mutant TET2
abstract: TET2 is a close relative of TET1, an enzyme that converts 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) in DNA. The gene encoding TET2 resides at chromosome 4q24, in a region showing recurrent microdeletions and copy-neutral loss of heterozygosity (CN-LOH) in patients with diverse myeloid malignancies. Somatic TET2 mutations are frequently observed in myelodysplastic syndromes (MDS), myeloproliferative neoplasms (MPN), MDS/MPN overlap syndromes including chronic myelomonocytic leukaemia (CMML), acute myeloid leukaemias (AML) and secondary AML (sAML). We show here that TET2 mutations associated with myeloid malignancies compromise catalytic activity. Bone marrow samples from patients with TET2 mutations displayed uniformly low levels of 5hmC in genomic DNA compared to bone marrow samples from healthy controls. Moreover, small hairpin RNA (shRNA)-mediated depletion of Tet2 in mouse haematopoietic precursors skewed their differentiation towards monocyte/macrophage lineages in culture. There was no significant difference in DNA methylation between bone marrow samples from patients with high 5hmC versus healthy controls, but samples from patients with low 5hmC showed hypomethylation relative to controls at the majority of differentially methylated CpG sites. Our results demonstrate that Tet2 is important for normal myelopoiesis, and suggest that disruption of TET2 enzymatic activity favours myeloid tumorigenesis. Measurement of 5hmC levels in myeloid malignancies may prove valuable as a diagnostic and prognostic tool, to tailor therapies and assess responses to anticancer drugs.

title: TET enzymes, TDG and the dynamics of DNA demethylation
abstract: DNA methylation has a profound impact on genome stability, transcription and development. Although enzymes that catalyse DNA methylation have been well characterized, those that are involved in methyl group removal have remained elusive, until recently. The transformative discovery that ten-eleven translocation (TET) family enzymes can oxidize 5-methylcytosine has greatly advanced our understanding of DNA demethylation. 5-Hydroxymethylcytosine is a key nexus in demethylation that can either be passively depleted through DNA replication or actively reverted to cytosine through iterative oxidation and thymine DNA glycosylase (TDG)-mediated base excision repair. Methylation, oxidation and repair now offer a model for a complete cycle of dynamic cytosine modification, with mounting evidence for its significance in the biological processes known to involve active demethylation.

title: A unique regulatory phase of DNA methylation in the early mammalian embryo
abstract: DNA methylation is highly dynamic during mammalian embryogenesis. It is broadly accepted that the paternal genome is actively depleted of 5-methylcytosine at fertilization, followed by passive loss that reaches a minimum at the blastocyst stage. However, this model is based on limited data, and so far no base-resolution maps exist to support and refine it. Here we generate genome-scale DNA methylation maps in mouse gametes and from the zygote through post-implantation. We find that the oocyte already exhibits global hypomethylation, particularly at specific families of long interspersed element 1 and long terminal repeat retroelements, which are disparately methylated between gametes and have lower methylation values in the zygote than in sperm. Surprisingly, the oocyte contributes a unique set of differentially methylated regions (DMRs)—including many CpG island promoters—that are maintained in the early embryo but are lost upon specification and absent from somatic cells. In contrast, sperm-contributed DMRs are largely intergenic and become hypermethylated after the blastocyst stage. Our data provide a genome-scale, base-resolution timeline of DNA methylation in the pre-specified embryo, when this epigenetic modification is most dynamic, before returning to the canonical somatic pattern.

title: GREAT improves functional interpretation of cis-regulatory regions
abstract: We developed the Genomic Regions Enrichment of Annotations Tool (GREAT) to analyze the functional significance of cis-regulatory regions identified by localized measurements of DNA binding events across an entire genome. Whereas previous methods took into account only binding proximal to genes, GREAT is able to properly incorporate distal binding sites and control for false positives using a binomial test over the input genomic regions. GREAT incorporates annotations from 20 ontologies and is available as a web application. Applying GREAT to data sets from chromatin immunoprecipitation coupled with massively parallel sequencing (ChIP-seq) of multiple transcription-associated factors, including SRF, NRSF, GABP, Stat3 and p300 in different developmental contexts, we recover many functions of these factors that are missed by existing gene-based tools, and we generate testable hypotheses. The utility of GREAT is not limited to ChIP-seq, as it could also be applied to open chromatin, localized epigenomic markers and similar functional data sets, as well as comparative genomics sets.

title: Reversing DNA Methylation: Mechanisms, Genomics, and Biological Functions
abstract: Methylation of cytosines in the mammalian genome represents a key epigenetic modification and is dynamically regulated during development. Compelling evidence now suggests that dynamic regulation of DNA methylation is mainly achieved through a cyclic enzymatic cascade comprised of cytosine methylation, iterative oxidation of methyl group by TET dioxygenases, and restoration of unmodified cytosines by either replication-dependent dilution or DNA glycosylase-initiated base excision repair. In this review, we discuss the mechanism and function of DNA demethylation in mammalian genomes, focusing particularly on how developmental modulation of the cytosine-modifying pathway is coupled to active reversal of DNA methylation in diverse biological processes.

title: Quantitative Sequencing of 5-Methylcytosine and 5-Hydroxymethylcytosine at Single-Base Resolution
abstract: Distinguishing Epigenetic Marks Methylation of the cytosine base in eukaryotic DNA (5mC) is an important epigenetic mark involved in gene silencing and genome stability. Methylated cytosine can be enzymatically oxidized to 5-hydroxymethylcytosine (5hmC), which may function as a distinct epigenetic mark—possibly involved in pluripotency—and it may also be an intermediate in active DNA demethylation. To be able to detect 5hmC genome-wide and at single-base resolution, Booth et al. (p. 934, published online 26 April) developed a 5hmC sequencing chemistry that selectively oxidizes 5hmC to 5-formylcytosine and then to uracil while leaving 5mC unchanged. Using this method, mouse embryonic stem cell genomic DNA was sequenced to reveal that 5hmC is found enriched at intragenic CpG islands and long interspersed nuclear element–1 retrotransposons. A sequencing method can discriminate epigenetically modified cytosine nucleotides within embryonic stem cell DNA. 5-Methylcytosine can be converted to 5-hydroxymethylcytosine (5hmC) in mammalian DNA by the ten-eleven translocation (TET) enzymes. We introduce oxidative bisulfite sequencing (oxBS-Seq), the first method for quantitative mapping of 5hmC in genomic DNA at single-nucleotide resolution. Selective chemical oxidation of 5hmC to 5-formylcytosine (5fC) enables bisulfite conversion of 5fC to uracil. We demonstrate the utility of oxBS-Seq to map and quantify 5hmC at CpG islands (CGIs) in mouse embryonic stem (ES) cells and identify 800 5hmC-containing CGIs that have on average 3.3% hydroxymethylation. High levels of 5hmC were found in CGIs associated with transcriptional regulators and in long interspersed nuclear elements, suggesting that these regions might undergo epigenetic reprogramming in ES cells. Our results open new questions on 5hmC dynamics and sequence-specific targeting by TETs.

title: Oxidative bisulfite sequencing of 5-methylcytosine and 5-hydroxymethylcytosine
abstract: To uncover the function of and interplay between the mammalian cytosine modifications 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), new techniques and advances in current technology are needed. To this end, we have developed oxidative bisulfite sequencing (oxBS-seq), which can quantitatively locate 5mC and 5hmC marks at single-base resolution in genomic DNA. In bisulfite sequencing (BS-seq), both 5mC and 5hmC are read as cytosines and thus cannot be discriminated; however, in oxBS-seq, specific oxidation of 5hmC to 5-formylcytosine (5fC) and conversion of the newly formed 5fC to uracil (under bisulfite conditions) means that 5hmC can be discriminated from 5mC. A positive readout of actual 5mC is gained from a single oxBS-seq run, and 5hmC levels are inferred by comparison with a BS-seq run. Here we describe an optimized second-generation protocol that can be completed in 2 d.

title: The bHLH protein PTF1-p48 is essential for the formation of the exocrine and the correct spatial organization of the endocrine pancreas.
abstract: We have generated a mouse bearing a null allele of the gene encoding basic helix-loop-helix (bHLH) protein p48, the cell-specific DNA-binding subunit of hetero-oligomeric transcription factor PTF1 that directs the expression of genes in the exocrine pancreas. The null mutation, which establishes a lethal condition shortly after birth, leads to a complete absence of exocrine pancreatic tissue and its specific products, indicating that p48 is required for differentiation and/or proliferation of the exocrine cell lineage. p48 is so far the only developmental regulator known to be required exclusively for committing cells to an exocrine fate. The hormone secreting cells of all four endocrine lineages are present in the mesentery that normally harbors the pancreatic organ until day 16 of gestation. Toward the end of embryonic life, cells expressing endocrine functions are no longer detected at their original location but are now found to colonize the spleen, where they persist in a functional state until postnatal death of the organism occurs. These findings suggest that the presence of the exocrine pancreas is required for the correct spatial assembly of the endocrine pancreas and that, in its absence, endocrine cells are directed by default to the spleen, a site that, in some reptiles, harbors part of this particular cellular compartment.

title: The role of the transcriptional regulator Ptf1a in converting intestinal to pancreatic progenitors
abstract: Pancreas development begins with the formation of buds at specific sites in the embryonic foregut endoderm. We used recombination-based lineage tracing in vivo to show that Ptf1a (also known as PTF1-p48) is expressed at these early stages in the progenitors of pancreatic ducts, exocrine and endocrine cells, rather than being an exocrine-specific gene as previously described. Moreover, inactivation of Ptf1a switches the character of pancreatic progenitors such that their progeny proliferate in and adopt the normal fates of duodenal epithelium, including its stem-cell compartment. Consistent with the proposal that Ptf1a supports the specification of precursors of all three pancreatic cell types, transgene-based expression of Pdx1, a gene essential to pancreas formation, from Ptf1a cis-regulatory sequences restores pancreas tissue to Pdx1-null mice that otherwise lack mature exocrine and endocrine cells because of an early arrest in organogenesis. These experiments provide evidence that Ptf1a expression is specifically connected to the acquisition of pancreatic fate by undifferentiated foregut endoderm.

title: Transcription factors as readers and effectors of DNA methylation
abstract: Recent technological advances have made it possible to decode DNA methylomes at single-base-pair resolution under various physiological conditions. Many aberrant or differentially methylated sites have been discovered, but the mechanisms by which changes in DNA methylation lead to observed phenotypes, such as cancer, remain elusive. The classical view of methylation-mediated protein–DNA interactions is that only proteins with a methyl-CpG binding domain (MBD) can interact with methylated DNA. However, evidence is emerging to suggest that transcription factors lacking a MBD can also interact with methylated DNA. The identification of these proteins and the elucidation of their characteristics and the biological consequences of methylation-dependent transcription factor–DNA interactions are important stepping stones towards a mechanistic understanding of methylation-mediated biological processes, which have crucial implications for human development and disease.

title: CUX1: a modulator of tumour aggressiveness in pancreatic neuroendocrine neoplasms.
abstract: Pancreatic neuroendocrine neoplasms (PNENs) constitute a rare tumour entity, and prognosis and treatment options depend on tumour-mediating hallmarks such as angiogenesis, proliferation rate and resistance to apoptosis. The molecular pathways that determine the malignant phenotype are still insufficiently understood and this has limited the use of effective combination therapies in the past. In this study, we aimed to characterise the effect of the oncogenic transcription factor Cut homeobox 1 (CUX1) on proliferation, resistance to apoptosis and angiogenesis in murine and human PNENs. The expression and function of CUX1 were analysed using knockdown and overexpression strategies in Ins-1 and Bon-1 cells, xenograft models and a genetically engineered mouse model of insulinoma (RIP1Tag2). Regulation of angiogenesis was assessed using RNA profiling and functional tube-formation assays in HMEC-1 cells. Finally, CUX1 expression was assessed in a tissue microarray of 59 human insulinomas and correlated with clinicopathological data. CUX1 expression was upregulated during tumour progression in a time- and stage-dependent manner in the RIP1Tag2 model, and associated with pro-invasive and metastatic features of human insulinomas. Endogenous and recombinant CUX1 expression increased tumour cell proliferation, tumour growth, resistance to apoptosis, and angiogenesis in vitro and in vivo. Mechanistically, the pro-angiogenic effect of CUX1 was mediated via upregulation of effectors such as HIF1α and MMP9. CUX1 mediates an invasive pro-angiogenic phenotype and is associated with malignant behaviour in human insulinomas.

title: Sox9 coordinates a transcriptional network in pancreatic progenitor cells
abstract: During pancreas development, both the exocrine and endocrine lineages differentiate from a common pool of progenitor cells with similarities to mature pancreatic duct cells. A small set of transcription factors, including Tcf2, Onecut1, and Foxa2, has been identified in these pancreatic progenitor cells. The Sry/HMG box transcription factor Sox9 is also expressed in the early pancreatic epithelium and is required for normal pancreatic exocrine and endocrine development in humans. In this study, we found Sox9 in mice specifically expressed with the other progenitor transcription factors in both pancreatic progenitor cells and duct cells in the adult pancreas. Sox9 directly bound to all three genes in vitro and in intact cells, and regulated their expression. In turn, both Foxa2 and Tcf2 regulated Sox9 expression, demonstrating feedback circuits between these genes. Furthermore, Sox9 activated the expression of the proendocrine factor Neurogenin3, which also depends on the other members of the progenitor transcription network. These studies indicate that Sox9 plays a dual role in pancreatic progenitor cells: both maintaining a stable transcriptional network and supporting the programs by which these cells differentiate into distinct lineages.

title: Large conserved domains of low DNA methylation maintained by Dnmt3a
abstract: Gains and losses in DNA methylation are prominent features of mammalian cell types. To gain insight into the mechanisms that promote shifts in DNA methylation and contribute to changes in cell fate, including malignant transformation, we performed genome-wide mapping of 5-methylcytosine and 5-hydroxymethylcytosine in purified mouse hematopoietic stem cells. We discovered extended regions of low methylation (canyons) that span conserved domains frequently containing transcription factors and are distinct from CpG islands and shores. About half of the genes in these methylation canyons are coated with repressive histone marks, whereas the remainder are covered by activating histone marks and are highly expressed in hematopoietic stem cells (HSCs). Canyon borders are demarked by 5-hydroxymethylcytosine and become eroded in the absence of DNA methyltransferase 3a (Dnmt3a). Genes dysregulated in human leukemias are enriched for canyon-associated genes. The new epigenetic landscape we describe may provide a mechanism for the regulation of hematopoiesis and may contribute to leukemia development.

title: Broad H3K4me3 is associated with increased transcription elongation and enhancer activity at tumor-suppressor genes
abstract: Tumor suppressors are mostly defined by inactivating mutations in tumors, yet little is known about their epigenetic features in normal cells. Through integrative analysis of 1,134 genome-wide epigenetic profiles, mutations from >8,200 tumor-normal pairs and our experimental data from clinical samples, we discovered broad peaks for trimethylation of histone H3 at lysine 4 (H3K4me3; wider than 4 kb) as the first epigenetic signature for tumor suppressors in normal cells. Broad H3K4me3 is associated with increased transcription elongation and enhancer activity, which together lead to exceptionally high gene expression, and is distinct from other broad epigenetic features, such as super-enhancers. Genes with broad H3K4me3 peaks conserved across normal cells may represent pan-cancer tumor suppressors, such as TP53 and PTEN, whereas genes with cell type–specific broad H3K4me3 peaks may represent cell identity genes and cell type–specific tumor suppressors. Furthermore, widespread shortening of broad H3K4me3 peaks in cancers is associated with repression of tumor suppressors. Thus, the broad H3K4me3 epigenetic signature provides mutation-independent information for the discovery and characterization of new tumor suppressors.

title: Distinct features of H3K4me3 and H3K27me3 chromatin domains in pre-implantation embryos
abstract: Histone modifications have critical roles in regulating the expression of developmental genes during embryo development in mammals. However, genome-wide analyses of histone modifications in pre-implantation embryos have been impeded by the scarcity of the required materials. Here, by using a small-scale chromatin immunoprecipitation followed by sequencing (ChIP–seq) method, we map the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 trimethylation (H3K27me3), which are associated with gene activation and repression, respectively, in mouse pre-implantation embryos. We find that the re-establishment of H3K4me3, especially on promoter regions, occurs much more rapidly than that of H3K27me3 following fertilization, which is consistent with the major wave of zygotic genome activation at the two-cell stage. Furthermore, H3K4me3 and H3K27me3 possess distinct features of sequence preference and dynamics in pre-implantation embryos. Although H3K4me3 modifications occur consistently at transcription start sites, the breadth of the H3K4me3 domain is a highly dynamic feature. Notably, the broad H3K4me3 domain (wider than 5 kb) is associated with higher transcription activity and cell identity not only in pre-implantation development but also in the process of deriving embryonic stem cells from the inner cell mass and trophoblast stem cells from the trophectoderm. Compared to embryonic stem cells, we found that the bivalency (that is, co-occurrence of H3K4me3 and H3K27me3) in early embryos is relatively infrequent and unstable. Taken together, our results provide a genome-wide map of H3K4me3 and H3K27me3 modifications in pre-implantation embryos, facilitating further exploration of the mechanism for epigenetic regulation in early embryos.

title: Rfx6 Directs Islet Formation and Insulin Production in Mice and Humans
abstract: Insulin from the β-cells of the pancreatic islets of Langerhans controls energy homeostasis in vertebrates, and its deficiency causes diabetes mellitus. During embryonic development, the transcription factor neurogenin 3 (Neurog3) initiates the differentiation of the β-cells and other islet cell types from pancreatic endoderm, but the genetic program that subsequently completes this differentiation remains incompletely understood. Here we show that the transcription factor Rfx6 directs islet cell differentiation downstream of Neurog3. Mice lacking Rfx6 failed to generate any of the normal islet cell types except for pancreatic-polypeptide-producing cells. In human infants with a similar autosomal recessive syndrome of neonatal diabetes, genetic mapping and subsequent sequencing identified mutations in the human RFX6 gene. These studies demonstrate a unique position for Rfx6 in the hierarchy of factors that coordinate pancreatic islet development in both mice and humans. Rfx6 could prove useful in efforts to generate β-cells for patients with diabetes.

title: Establishment and Maintenance of Genomic Methylation Patterns in Mouse Embryonic Stem Cells by Dnmt3a and Dnmt3b
abstract: We have previously shown that the DNA methyltransferases Dnmt3a and Dnmt3b carry out de novo methylation of the mouse genome during early postimplantation development and of maternally imprinted genes in the oocyte. In the present study, we demonstrate that Dnmt3a and Dnmt3b are also essential for the stable inheritance, or “maintenance,” of DNA methylation patterns. Inactivation of both Dnmt3a and Dnmt3b in embryonic stem (ES) cells results in progressive loss of methylation in various repeats and single-copy genes. Interestingly, introduction of the Dnmt3a, Dnmt3a2, and Dnmt3b1 isoforms back into highly demethylated mutant ES cells restores genomic methylation patterns; these isoforms appear to have both common and distinct DNA targets, but they all fail to restore the maternal methylation imprints. In contrast, overexpression of Dnmt1 and Dnmt3b3 failed to restore DNA methylation patterns due to their inability to catalyze de novo methylation in vivo. We also show that hypermethylation of genomic DNA by Dnmt3a and Dnmt3b is necessary for ES cells to form teratomas in nude mice. These results indicate that genomic methylation patterns are determined partly through differential expression of different Dnmt3a and Dnmt3b isoforms.

title: Severe Global DNA Hypomethylation Blocks Differentiation and Induces Histone Hyperacetylation in Embryonic Stem Cells
abstract: It has been reported that DNA methyltransferase 1-deficient (Dnmt1−/−) embryonic stem (ES) cells are hypomethylated (20% CpG methylation) and die through apoptosis when induced to differentiate. Here, we show that Dnmt[3a−/−,3b−/−] ES cells with just 0.6% of their CpG dinucleotides behave differently: the majority of cells within the culture are partially or completely blocked in their ability to initiate differentiation, remaining viable while retaining the stem cell characteristics of alkaline phosphatase and Oct4 expression. Restoration of DNA methylation levels rescues these defects. Severely hypomethylated Dnmt[3a−/−,3b−/−] ES cells have increased histone acetylation levels, and those cells that can differentiate aberrantly express extraembryonic markers of differentiation. Dnmt[3a−/−,3b−/−] ES cells with >10% CpG methylation are able to terminally differentiate, whereas Dnmt1−/− ES cells with 20% of the CpG methylated cannot differentiate. This demonstrates that successful terminal differentiation is not dependent simply on adequate methylation levels. There is an absolute requirement that the methylation be delivered by the maintenance enzyme Dnmt1.

title: Tet2 facilitates the derepression of myeloid target genes during CEBPα-induced transdifferentiation of pre-B cells.
abstract: The methylcytosine hydroxylase Tet2 has been implicated in hematopoietic differentiation and the formation of myeloid malignancies when mutated. An ideal system to study the role of Tet2 in myelopoeisis is CEBPα-induced transdifferentiation of pre-B cells into macrophages. Here we found that CEBPα binds to upstream regions of Tet2 and that the gene becomes activated. Tet2 knockdowns impaired the upregulation of macrophage markers as well as phagocytic capacity, suggesting that the enzyme is required for both early and late stage myeloid differentiation. A slightly weaker effect was seen in primary cells with a Tet2 ablation. Expression arrays of transdifferentiating cells with Tet2 knockdowns permitted the identification of a small subset of myeloid genes whose upregulation was blunted. Activation of these target genes was accompanied by rapid increases of promoter hydroxy-methylation. Our observations indicate that Tet2 helps CEBPα rapidly derepress myeloid genes during the conversion of pre-B cells into macrophages.

title: Genome-wide Analysis Identifies Bcl6-Controlled Regulatory Networks during T Follicular Helper Cell Differentiation.
abstract: T follicular helper (Tfh) cell is a unique T cell subset specialized in promoting humoral immunity. B-cell lymphoma 6 protein (Bcl6) has been identified as an obligatory transcription factor in Tfh cells; however, the molecular mechanism underlying Bcl6 function remains largely unknown. Here, we defined Bcl6 target genes in Tfh cells by analyzing genome-wide Bcl6 occupancy together with transcriptome profiling. With consensus sequences being different from those in Th9, B cells, and macrophages, Bcl6 binding in Tfh cell was closely associated with a decrease in 5-hydroxymethylcytosine (5hmC). Importantly, Bcl6 promoted Tfh cell differentiation through antagonizing IL-7R (CD127)/signal transducer and activator of transcription (STAT) 5 axis; deletion of the Bcl6 gene in T cells resulted in enhanced IL-7R-STAT5 signaling and substantial expansion of CD127(hi) non-Tfh cells. Thus, our study systemically examines Bcl6-controlled regulatory networks and provides important insights into Bcl6's biological functions in Tfh cells.

title: Targeted DNA demethylation and activation of endogenous genes using programmable TALE-TET1 fusion proteins
abstract: Genome-wide studies have defined cell type–specific patterns of DNA methylation that are important for regulating gene expression in both normal development and disease. However, determining the functional significance of specific methylation events remains challenging, owing to the lack of methods for removing such modifications in a targeted manner. Here we describe an approach for efficient targeted demethylation of specific CpGs in human cells using fusions of engineered transcription activator–like effector (TALE) repeat arrays and the TET1 hydroxylase catalytic domain. Using these TALE-TET1 fusions, we demonstrate that modification of critical methylated promoter CpG positions can lead to substantial increases in the expression of endogenous human genes. Our results delineate a strategy for understanding the functional significance of specific CpG methylation marks in the context of endogenous gene loci and validate programmable DNA demethylation reagents with potential utility for research and therapeutic applications.

title: Engineered
Split-TET2 Enzyme for Inducible Epigenetic Remodeling
abstract: The Ten-eleven translocation (TET) family of 5-methylcytosine (5mC) dioxygenases catalyze the conversion of 5mC into 5-hydroxymethylcytosine (5hmC) and further oxidized species to promote active DNA demethylation. Here we engineered a split-TET2 enzyme to enable temporal control of 5mC oxidation and subsequent remodeling of epigenetic states in mammalian cells. We further demonstrate the use of this chemically inducible system to dissect the correlation between DNA hydroxymethylation and chromatin accessibility in the mammalian genome. This chemical-inducible epigenome remodeling tool will find broad use in interrogating cellular systems without altering the genetic code, as well as in probing the epigenotype-phenotype relations in various biological systems.

title: Dynamic regulation of Pdx1 enhancers by Foxa1 and Foxa2 is essential for pancreas development.
abstract: The onset of pancreas development in the foregut endoderm is marked by activation of the homeobox gene Pdx1 (IPF1). Pdx1 is essential for the expansion of the pancreatic primordium and the development of endocrine islets. The control of Pdx1 expression has been only partially elucidated. We demonstrate here that the winged-helix transcription factors Foxa1 and Foxa2 co-occupy multiple regulatory domains in the Pdx1 gene. Compound conditional ablation of both Foxa1 and Foxa2 in the pancreatic primordium results in complete loss of Pdx1 expression and severe pancreatic hypoplasia. Mutant mice exhibit hyperglycemia with severely disrupted acinar and islet development, and die shortly after birth. Assessment of developmental markers in the mutant pancreas revealed a failure in the expansion of the pancreatic anlage, a blockage of exocrine and endocrine cell differentiation, and an arrest at the primitive duct stage. Comparing their relative developmental activity, we find that Foxa2 is the major regulator in promoting pancreas development and cell differentiation. Using chromatin immunoprecipitations (ChIP) and ChIP sequencing (ChIPSeq) of fetal pancreas and islet chromatin, we demonstrate that Foxa1 and Foxa2 predominantly occupy a distal enhancer at -6.4 kb relative to the transcriptional start site in the Pdx1 gene. In addition, occupancy of the well-characterized proximal Pdx1 enhancer by Foxa1 and Foxa2 is developmental stage-dependent. Thus, the regulation of Pdx1 expression by Foxa1 and Foxa2 is a key early event controlling the expansion and differentiation of the pancreatic primordia.

title: Nucleation of DNA repair factors by FOXA1 links DNA demethylation to transcriptional pioneering
abstract: FOXA1 functions in epigenetic reprogramming and is described as a 'pioneer factor'. However, exactly how FOXA1 achieves these remarkable biological functions is not fully understood. Here we report that FOXA1 associates with DNA repair complexes and is required for genomic targeting of DNA polymerase β (POLB) in human cells. Genome-wide DNA methylomes demonstrate that the FOXA1 DNA repair complex is functionally linked to DNA demethylation in a lineage-specific fashion. Depletion of FOXA1 results in localized reestablishment of methylation in a large portion of FOXA1-bound regions, and the regions with the most consistent hypermethylation exhibit the greatest loss of POLB and are represented by active promoters and enhancers. Consistently, overexpression of FOXA1 commits its binding sites to active DNA demethylation in a POLB-dependent manner. Finally, FOXA1-associated DNA demethylation is tightly coupled with estrogen receptor genomic targeting and estrogen responsiveness. Together, these results link FOXA1-associated DNA demethylation to transcriptional pioneering by FOXA1.

title: Insulator dysfunction and oncogene activation in IDH mutant gliomas
abstract: Gain-of-function IDH mutations are initiating events that define major clinical and prognostic classes of gliomas. Mutant IDH protein produces a new onco-metabolite, 2-hydroxyglutarate, which interferes with iron-dependent hydroxylases, including the TET family of 5′-methylcytosine hydroxylases. TET enzymes catalyse a key step in the removal of DNA methylation. IDH mutant gliomas thus manifest a CpG island methylator phenotype (G-CIMP), although the functional importance of this altered epigenetic state remains unclear. Here we show that human IDH mutant gliomas exhibit hypermethylation at cohesin and CCCTC-binding factor (CTCF)-binding sites, compromising binding of this methylation-sensitive insulator protein. Reduced CTCF binding is associated with loss of insulation between topological domains and aberrant gene activation. We specifically demonstrate that loss of CTCF at a domain boundary permits a constitutive enhancer to interact aberrantly with the receptor tyrosine kinase gene PDGFRA, a prominent glioma oncogene. Treatment of IDH mutant gliomaspheres with a demethylating agent partially restores insulator function and downregulates PDGFRA. Conversely, CRISPR-mediated disruption of the CTCF motif in IDH wild-type gliomaspheres upregulates PDGFRA and increases proliferation. Our study suggests that IDH mutations promote gliomagenesis by disrupting chromosomal topology and allowing aberrant regulatory interactions that induce oncogene expression.

title: Tumor development is associated with decrease of TET gene expression and 5-methylcytosine hydroxylation
abstract: The TET (ten–eleven translocation) family of α-ketoglutarate (α-KG)-dependent dioxygenases catalyzes the sequential oxidation of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine and 5-carboxylcytosine, leading to eventual DNA demethylation. The TET2 gene is a bona fide tumor suppressor frequently mutated in leukemia, and TET enzyme activity is inhibited in IDH1/2-mutated tumors by the oncometabolite 2-hydroxyglutarate, an antagonist of α-KG, linking 5mC oxidation to cancer development. We report here that the levels of 5hmC are dramatically reduced in human breast, liver, lung, pancreatic and prostate cancers when compared with the matched surrounding normal tissues. Associated with the 5hmC decrease is the substantial reduction of the expression of all three TET genes, revealing a possible mechanism for the reduced 5hmC in cancer cells. The decrease of 5hmC was also observed during tumor development in different genetically engineered mouse models. Together, our results identify 5hmC as a biomarker whose decrease is broadly and tightly associated with tumor development.