Two are better than one: HPoxBS - hairpin oxidative bisulfite sequencing.

The controlled and stepwise oxidation of 5mC to 5hmC, 5fC and 5caC by Tet enzymes is influencing the chemical and biological properties of cytosine. Besides direct effects on gene regulation, oxidised forms influence the dynamics of demethylation and re-methylation processes. So far, no combined methods exist which allow to precisely determine the strand specific localisation of cytosine modifications along with their CpG symmetric distribution. Here we describe a comprehensive protocol combining conventional hairpin bisulfite with oxidative bisulfite sequencing (HPoxBS) to determine the strand specific distribution of 5mC and 5hmC at base resolution. We apply this method to analyse the contribution of local oxidative effects on DNA demethylation in mouse ES cells. Our method includes the HPoxBS workflow and subsequent data analysis using our developed software tools. Besides a precise estimation and display of strand specific 5mC and 5hmC levels at base resolution we apply the data to predict region specific activities of Dnmt and Tet enzymes. Our experimental and computational workflow provides a precise double strand display of 5mC and 5hmC modifications at single base resolution. Based on our data we predict region specific Tet and Dnmt enzyme efficiencies shaping the distinct locus levels and patterns of 5hmC and 5mC.

Genomic integrity of ground-state pluripotency.

Pluripotent cells appear to be in a transient state during early development. These cells have the capability to transition into embryonic stem cells (ESCs). It has been reported that mouse pluripotent cells cultivated in chemically defined media sustain the ground state of pluripotency. Because the epigenetic pattern of pluripotent cells reflects their environment, culture under different conditions causes epigenetic changes, which could lead to genomic instability. This study focused on the DNA methylation pattern of repetitive elements (REs) and their activation levels under two ground-state conditions and assessed the genomic integrity of ESCs. We measured the methylation and expression level of REs in different media. The results indicated that although the ground-state conditions show higher REs activity, they did not lead to DNA damage; therefore, the level of genomic instability is lower under the ground-state compared with the conventional condition. Our results indicated that when choosing an optimum condition, different features of the condition must be considered to have epigenetically and genomically stable stem cells.

H(O)TA: estimation of DNA methylation and hydroxylation levels and efficiencies from time course data.

MOTIVATION: Methylation and hydroxylation of cytosines to form 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) belong to the most important epigenetic modifications and their vital role in the regulation of gene expression has been widely recognized. Recent experimental techniques allow to infer methylation and hydroxylation levels at CpG dinucleotides but require a sophisticated statistical analysis to achieve accurate estimates. RESULTS: We present H(O)TA, a software tool based on a stochastic modeling approach, which simultaneously analyzes time course data from hairpin bisulfite sequencing and hairpin oxidative bisulfite sequencing. AVAILABILITY AND IMPLEMENTATION: : https://mosi.uni-saarland.de/HOTA. CONTACT: charalampos.kyriakopoulos@uni-saarland.de or verena.wolf@uni-saarland.de.

The Influence of Hydroxylation on Maintaining CpG Methylation Patterns: A Hidden Markov Model Approach.

DNA methylation and demethylation are opposing processes that when in balance create stable patterns of epigenetic memory. The control of DNA methylation pattern formation by replication dependent and independent demethylation processes has been suggested to be influenced by Tet mediated oxidation of 5mC. Several alternative mechanisms have been proposed suggesting that 5hmC influences either replication dependent maintenance of DNA methylation or replication independent processes of active demethylation. Using high resolution hairpin oxidative bisulfite sequencing data, we precisely determine the amount of 5mC and 5hmC and model the contribution of 5hmC to processes of demethylation in mouse ESCs. We develop an extended hidden Markov model capable of accurately describing the regional contribution of 5hmC to demethylation dynamics. Our analysis shows that 5hmC has a strong impact on replication dependent demethylation, mainly by impairing methylation maintenance.

Active demethylation of the Foxp3 locus leads to the generation of stable regulatory T cells within the thymus.

Stable expression of Foxp3 in regulatory T cells (Tregs) depends on DNA demethylation at the Treg-specific demethylated region (TSDR), a conserved, CpG-rich region within the Foxp3 locus. The TSDR is selectively demethylated in ex vivo Tregs purified from secondary lymphoid organs, but it is unclear at which stage of Treg development demethylation takes place. In this study, we show that commitment to a stable lineage occurred during early stages of murine thymic Treg development by engraving of lineage-specific epigenetic marks in parallel with establishment of a Treg-specific gene expression profile. TSDR demethylation was achieved through an active mechanism and involved enzymes of the ten-eleven-translocation family and hydroxylation of methylated cytosines, a modification that is implicated as an initiating step of mitosis-independent DNA demethylation pathways and has not yet been observed at specific loci during immune cell differentiation. Together, our results demonstrate that initiating TSDR demethylation during early stages of thymic Treg development commences stabilization of Foxp3 expression and guarantees full functionality and long-term lineage stability of Tregs.

A comprehensive approach for genome-wide efficiency profiling of DNA modifying enzymes.

A precise understanding of DNA methylation dynamics is of great importance for a variety of biological processes including cellular reprogramming and differentiation. To date, complex integration of multiple and distinct genome-wide datasets is required to realize this task. We present GwEEP (genome-wide epigenetic efficiency profiling) a versatile approach to infer dynamic efficiencies of DNA modifying enzymes. GwEEP relies on genome-wide hairpin datasets, which are translated by a hidden Markov model into quantitative enzyme efficiencies with reported confidence around the estimates. GwEEP predicts de novo and maintenance methylation efficiencies of Dnmts and furthermore the hydroxylation efficiency of Tets. Its design also allows capturing further oxidation processes given available data. We show that GwEEP predicts accurately the epigenetic changes of ESCs following a Serum-to-2i shift and applied to Tet TKO cells confirms the hypothesized mutual interference between Dnmts and Tets.

Dnmt1 has de novo activity targeted to transposable elements.

DNA methylation plays a critical role during development, particularly in repressing retrotransposons. The mammalian methylation landscape is dependent on the combined activities of the canonical maintenance enzyme Dnmt1 and the de novo Dnmts, 3a and 3b. Here, we demonstrate that Dnmt1 displays de novo methylation activity in vitro and in vivo with specific retrotransposon targeting. We used whole-genome bisulfite and long-read Nanopore sequencing in genetically engineered methylation-depleted mouse embryonic stem cells to provide an in-depth assessment and quantification of this activity. Utilizing additional knockout lines and molecular characterization, we show that the de novo methylation activity of Dnmt1 depends on Uhrf1, and its genomic recruitment overlaps with regions that enrich for Uhrf1, Trim28 and H3K9 trimethylation. Our data demonstrate that Dnmt1 can catalyze DNA methylation in both a de novo and maintenance context, especially at retrotransposons, where this mechanism may provide additional stability for long-term repression and epigenetic propagation throughout development.

Hidden Markov Modelling Reveals Neighborhood Dependence of Dnmt3a and 3b Activity.

DNA methylation is an epigenetic mark whose important role in development has been widely recognized. This epigenetic modification results in heritable information not encoded by the DNA sequence. The underlying mechanisms controlling DNA methylation are only partly understood. Several mechanistic models of enzyme activities responsible for DNA methylation have been proposed. Here, we extend existing Hidden Markov Models (HMMs) for DNA methylation by describing the occurrence of spatial methylation patterns over time and propose several models with different neighborhood dependences. Furthermore, we investigate correlations between the neighborhood dependence and other genomic information. We perform numerical analysis of the HMMs applied to comprehensive hairpin and non-hairpin bisulfite sequencing measurements and accurately predict wild-type data. We find evidence that the activities of Dnmt3a and Dnmt3b responsible for de novo methylation depend on 5' (left) but not on 3' (right) neighboring CpGs in a sequencing string.

Hairpin Bisulfite Sequencing: Synchronous Methylation Analysis on Complementary DNA Strands of Individual Chromosomes.

The accurate and quantitative detection of 5-methylcytosine is of great importance in the field of epigenetics. The method of choice is usually bisulfite sequencing because of the high resolution and the possibility to combine it with next generation sequencing. Nevertheless, also this method has its limitations. Following the bisulfite treatment DNA strands are no longer complementary such that in a subsequent PCR amplification the DNA methylation patterns information of only one of the two DNA strand is preserved. Several years ago Hairpin Bisulfite sequencing was developed as a method to obtain the pattern information on complementary DNA strands. The method requires fragmentation (usually by enzymatic cleavage) of genomic DNA followed by a covalent linking of both DNA strands through ligation of a short DNA hairpin oligonucleotide to both strands. The ligated covalently linked dsDNA products are then subjected to a conventional bisulfite treatment during which all unmodified cytosines are converted to uracils. During the treatment the DNA is denatured forming noncomplementary ssDNA circles. These circles serve as a template for a locus specific PCR to amplify chromosomal patterns of the region of interest. As a result one ends up with a linearized product, which contains the methylation information of both complementary DNA strands.

CD137+CD154- Expression As a Regulatory T Cell (Treg)-Specific Activation Signature for Identification and Sorting of Stable Human Tregs from In Vitro Expansion Cultures.

Regulatory T cells (Tregs) are an attractive therapeutic tool for several different immune pathologies. Therapeutic Treg application often requires prolonged in vitro culture to generate sufficient Treg numbers or to optimize their functionality, e.g., via genetic engineering of their antigen receptors. However, purity of clinical Treg expansion cultures is highly variable, and currently, it is impossible to identify and separate stable Tregs from contaminating effector T cells, either ex vivo or after prior expansion. This represents a major obstacle for quality assurance of expanded Tregs and raises significant safety concerns. Here, we describe a Treg activation signature that allows identification and sorting of epigenetically imprinted Tregs even after prolonged in vitro culture. We show that short-term reactivation resulted in expression of CD137 but not CD154 on stable FoxP3+ Tregs that displayed a demethylated Treg-specific demethylated region, high suppressive potential, and lack of inflammatory cytokine expression. We also applied this Treg activation signature for rapid testing of chimeric antigen receptor functionality in human Tregs and identified major differences in the signaling requirements regarding CD137 versus CD28 costimulation. Taken together, CD137+CD154- expression emerges as a universal Treg activation signature ex vivo and upon in vitro expansion allowing the identification and isolation of epigenetically stable antigen-activated Tregs and providing a means for their rapid functional testing in vitro.

Histone and DNA methylation control by H3 serine 10/threonine 11 phosphorylation in the mouse zygote.

BACKGROUND: In the mammalian zygote, epigenetic reprogramming is a tightly controlled process of coordinated alterations of histone and DNA modifications. The parental genomes of the zygote show distinct patterns of histone H3 variants and distinct patterns of DNA and histone modifications. The molecular mechanisms linking histone variant-specific modifications and DNA methylation reprogramming during the first cell cycle remain to be clarified. RESULTS: Here, we show that the degree and distribution of H3K9me2 and of DNA modifications (5mC/5hmC) are influenced by the phosphorylation status of H3S10 and H3T11. The overexpression of the mutated histone variants H3.1 and 3.2 at either serine 10 or threonine 11 causes a decrease in H3K9me2 and 5mC and a concomitant increase in 5hmC in the maternal genome. Bisulphite sequencing results indicate an increase in hemimethylated CpG positions following H3.1T10A overexpression suggesting an impact of H3S10 and H3T11 phosphorylation on DNA methylation maintenance. CONCLUSIONS: Our data suggest a crosstalk between the cell-cycle-dependent control of S10 and T11 phosphorylation of histone variants H3.1 and H3.2 and the maintenance of the heterochromatic mark H3K9me2. This histone H3 "phospho-methylation switch" also influences the oxidative control of DNA methylation in the mouse zygote.

MERVL/Zscan4 Network Activation Results in Transient Genome-wide DNA Demethylation of mESCs.

Mouse embryonic stem cells are dynamic and heterogeneous. For example, rare cells cycle through a state characterized by decondensed chromatin and expression of transcripts, including the Zscan4 cluster and MERVL endogenous retrovirus, which are usually restricted to preimplantation embryos. Here, we further characterize the dynamics and consequences of this transient cell state. Single-cell transcriptomics identified the earliest upregulated transcripts as cells enter the MERVL/Zscan4 state. The MERVL/Zscan4 transcriptional network was also upregulated during induced pluripotent stem cell reprogramming. Genome-wide DNA methylation and chromatin analyses revealed global DNA hypomethylation accompanying increased chromatin accessibility. This transient DNA demethylation was driven by a loss of DNA methyltransferase proteins in the cells and occurred genome-wide. While methylation levels were restored once cells exit this state, genomic imprints remained hypomethylated, demonstrating a potential global and enduring influence of endogenous retroviral activation on the epigenome.