Genome-wide Distribution Of 5hmC Research Articles

Epigenetic reprogramming by TET enzymes impacts co-transcriptional R-loops.

DNA oxidation by ten-eleven translocation (TET) family enzymes is essential for epigenetic reprogramming. The conversion of 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC) initiates developmental and cell-type-specific transcriptional programs through mechanisms that include changes in the chromatin structure. Here, we show that the presence of 5hmC in the transcribed gene promotes the annealing of the nascent RNA to the template DNA strand, leading to the formation of an R-loop. Depletion of TET enzymes reduced global R-loops in the absence of gene expression changes, whereas CRISPR-mediated tethering of TET to an active gene promoted the formation of R-loops. The genome-wide distribution of 5hmC and R-loops shows a positive correlation in mouse and human stem cells and overlap in half of all active genes. Moreover, R-loop resolution leads to differential expression of a subset of genes that are involved in crucial events during stem cell proliferation. Altogether, our data reveal that epigenetic reprogramming via TET activity promotes co-transcriptional R-loop formation, disclosing new mechanisms of gene expression regulation.

Genome-wide distribution of 5hmC in the dental pulp of mouse molars and incisors

The dental pulp is critical for the production of odontoblasts to create reparative dentin. In recent years, dental pulp has become a promising source of mesenchymal stem cells that are capable of differentiating into multiple cell types. To elucidate the transcriptional control mechanisms specifying the early phases of odontoblast differentiation, we analysed the DNA demethylation pattern associated with 5-hydroxymethylcytosine (5hmC) in the primary murine dental pulp. 5hmC plays an important role in chromatin accessibility and transcriptional control by modelling a dynamic equilibrium between DNA methylation and demethylation. Our research revealed 5hmC enrichment along genes and non-coding regulatory regions associated with specific developmental pathways in the genome of mouse incisor and molar dental pulp. Although the overall distribution of 5hmC is similar, the intensity and location of the 5hmC peaks significantly differs between the incisor and molar pulp genome, indicating cell type-specific epigenetic variations. Our study suggests that the differential DNA demethylation pattern could account for the distinct regulatory mechanisms underlying the tooth-specific ontogenetic programs.

Genome-wide 5-hydroxymethylcytosine (5hmC) emerges at early stage of in vitro differentiation of a putative hepatocyte progenitor

A basic question linked to differential patterns of gene expression is how cells reach different fates despite using the same DNA template. Since 5-hydroxymethylcytosine (5hmC) emerged as an intermediate metabolite in active DNA demethylation, there have been increasing efforts to elucidate its function as a stable modification of the genome, including a role in establishing such tissue-specific patterns of expression. Recently we described TET1-mediated enrichment of 5hmC on the promoter region of the master regulator of hepatocyte identity, HNF4A, which precedes differentiation of liver adult progenitor cells in vitro. Here, we studied the genome-wide distribution of 5hmC at early in vitro differentiation of human hepatocyte-like cells. We found a global increase in 5hmC as well as a drop in 5-methylcytosine after one week of in vitro differentiation from bipotent progenitors, at a time when the liver transcript program is already established. 5hmC was overall higher at the bodies of overexpressed genes. Furthermore, by modifying the metabolic environment, an adenosine derivative prevents 5hmC enrichment and impairs the acquisition of hepatic identity markers. These results suggest that 5hmC could be a marker of cell identity, as well as a useful biomarker in conditions associated with cell de-differentiation such as liver malignancies.

TET2 Loss Dysregulates the Behavior of Mesenchymal Stem Cells and Increases Their Ability to Promote TET2−/−-Driven Myeloid Malignancy Progression

Background: TET2 is a methylcytosine dioxygenase that regulates cytosine hydroxymethylation. Although there are extensive data implicating a pivotal role of TET2 in hematopoietic stem/progenitor cell (HSPCs), the importance of TET2 in bone marrow mesenchymal stromal cells (BMSCs) remains unknown. Aim: Tet2 loss may dysregulate the integrity of bone marrow niche, by which affects the malignant progression. Methods: Generation and maintenance of Tet2 conditional knockout MiceAnalysis of 5-hmC and 5-mC levels using dot blotFlow cytometry analysis, cell sorting, and hematopoietic progenitor cell (HPC) assayMSC Culture and long-term coculture with HSPCsReal-time PCR and RAN-Seq AnalysisMeDIP-seq to profile the genome-wide distribution of 5hmC. Results: Here, we show that while both Tet2 f/f ;Mx1Cre (conditional Tet2-inactivation in hematopoietic cells) and Tet2 f/f ;MxfCre (germ line Tet2-inactivation) developed myeloid malignancies in mice, Tet2 f/f ;MxfCre mice had a significantly shortened survival compared with Tet2 f/f ;Mx1Cre mice. Interestingly, Tet2 −/− recipient mice exhibited a higher incident of myeloid malignancies and a significantly reduced survival rate compare with WT recipient mice. These data indicate that Tet2 −/− bone marrow niche might promote the progression of myeloid malignancies in Tet2 −/− mice. Strikingly, deletion of Tet2 in mesenchymal stem cells (MSCs) using Prx1-cre is associated with a significantly accelerated malignancy progression and shortened survival, suggesting that MSCs are the cell components in Tet2 −/− mice play a role in the initiation/progression of Tet2 −/−-driven myeloid malignancies. Furthermore, Tet2 −/− MSCs displayed a significantly increased self-renewal, proliferating and differentiation capability as assayed by the frequency of CFU-F and commitment toward osteoblasts. In addition, Tet2 −/− but not WT MSCs exhibited a significantly increased supportive capacity to Tet2 −/− HSC/HPC proliferation. RNA-sequencing analysis revealed that Tet2 −/− MSCs exhibited a distinct gene expression profiles with 468 dysregulated genes as compared with WT MSCs. Furthermore, the number of 5-hmC peaks were significantly decreased in Tet2 −/− MSCs compared with WT MSCs based on whole genomic 5-hmC profiling. The majority of TET2-dependent 5hmC modifications in MSCs are located within genes. We then examined TET2 gene expression in MSCs derived from human myeloproliferative neoplasms (MPN) patients and healthy individuals and found that TET2 and 5-hmC was moderately down-regulated in MPN MSCs as compared with healthy controls. Conclusion: These results highlight the critical role of TET2 in the maintenance of BMSC functions and osteoblast differentiation, and provide evidence that dysregulation of epigenetic modifier in BMSCs contributes to the progression of myeloid malignancies.

Abstract 5322: Normal breast tissue 5-hydroxymethylcytosine is enriched at DNA regulatory elements of breast cancer cells

Abstract DNA methylation is an established modulator of gene expression in normal mammalian cells. Emerging evidence suggests that the oxidized 5-methylcytosine (5mC) product, 5-hydroxymethylcytosine (5hmC), has distinct biological functions and its dysregulation may contribute to cancer pathogenesis. While 5hmC levels in normal breast tissue have been shown to be elevated relative to several other tissue types, we lack a complete knowledge of the genomic distribution of 5hmC in breast tissue, as well as its potential biological functions. Here, we determined the genome-wide distribution of 5hmC and 5mC using methylation array data from 18 normal breast tissue samples. To disambiguate 5hmC from 5mC, we performed parallel processing of bisulfite and oxidative-bisulfite treated DNA, followed by statistical estimation of the 5mC and 5hmC proportions. Based on the distribution of 5hmC, we identified genomic loci with consistently elevated levels of 5hmC across normal breast tissues. In addition, we showed 5hmC is enriched within transcriptionally active chromatin in breast myoepithelial cells, marked by H3K4 trimethylation. Conversely, 5hmC was depleted within active promoters (H3K9ac-marked DNA) and transcriptionally inactive chromatin (H3K9me3 and H3K27me3-marked DNA). Importantly, the loci with high 5hmC in normal breast tissue were enriched among regulatory elements, including transcription factor binding sites, in transformed normal-like breast MCF10A cells, as well as in multiple breast cancer cell lines. These findings provide a genome-wide distribution of 5hmC in normal breast tissue. Furthermore, these findings identify DNA-regulatory regions in breast cancer cell lines that show statistical enrichment for those sites with elevated 5hmC in normal breast tissue, suggesting deregulation of 5hmC within these elements may contribute to breast carcinogenesis. Citation Format: Owen Michael Wilkins, Kevin C. Johnson, E. Andres Houseman, Jessica E. King, Carmen J. Marsit, Brock C. Christensen. Normal breast tissue 5-hydroxymethylcytosine is enriched at DNA regulatory elements of breast cancer cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5322.

Epigenetic Optical Mapping of 5-Hydroxymethylcytosine in Nanochannel Arrays.

The epigenetic mark 5-hydroxymethylcytosine (5-hmC) is a distinct product of active DNA demethylation that is linked to gene regulation, development, and disease. In particular, 5-hmC levels dramatically decline in many cancers, potentially serving as an epigenetic biomarker. The noise associated with next-generation 5-hmC sequencing hinders reliable analysis of low 5-hmC containing tissues such as blood and malignant tumors. Additionally, genome-wide 5-hmC profiles generated by short-read sequencing are limited in providing long-range epigenetic information relevant to highly variable genomic regions, such as the 3.7 Mbp disease-related Human Leukocyte Antigen (HLA) region. We present a long-read, highly sensitive single-molecule mapping technology that generates hybrid genetic/epigenetic profiles of native chromosomal DNA. The genome-wide distribution of 5-hmC in human peripheral blood cells correlates well with 5-hmC DNA immunoprecipitation (hMeDIP) sequencing. However, the long single-molecule read-length of 100 kbp to 1 Mbp produces 5-hmC profiles across variable genomic regions that failed to show up in the sequencing data. In addition, optical 5-hmC mapping shows a strong correlation between the 5-hmC density in gene bodies and the corresponding level of gene expression. The single-molecule concept provides information on the distribution and coexistence of 5-hmC signals at multiple genomic loci on the same genomic DNA molecule, revealing long-range correlations and cell-to-cell epigenetic variation.

5-Hydroxymethylcytosine alterations in the human postmortem brains of autism spectrum disorder.

Autism spectrum disorders (ASDs) include a group of syndromes characterized by impaired language, social and communication skills, in addition to restrictive behaviors or stereotypes. However, with a prevalence of 1.5% in developed countries and high comorbidity rates, no clear underlying mechanism that unifies the heterogeneous phenotypes of ASD exists. 5-hydroxymethylcytosine (5hmC) is highly enriched in the brain and recognized as an essential epigenetic mark in developmental and brain disorders. To explore the role of 5hmC in ASD, we used the genomic DNA isolated from the postmortem cerebellum of both ASD patients and age-matched controls to profile genome-wide distribution of 5hmC. We identified 797 age-dependent differentially hydroxymethylated regions (DhMRs) in the young group (age ≤ 18), while no significant DhMR was identified in the groups over 18 years of age. Pathway and disease association analyses demonstrated that the intragenic DhMRs were in the genes involved in cell-cell communication and neurological disorders. Also, we saw significant 5hmC changes in the larger group of psychiatric genes. Interestingly, we found that the predicted cis functions of non-coding intergenic DhMRs strikingly associate with ASD and intellectual disorders. A significant fraction of intergenic DhMRs overlapped with topologically associating domains. These results together suggest that 5hmC alteration is associated with ASD, particularly in the early development stage, and could contribute to the pathogenesis of ASD.

Profiling of oxBS-450K 5-hydroxymethylcytosine in human placenta and brain reveals enrichment at imprinted loci

ABSTRACTDNA methylation (5-methylcytosine, 5 mC) is involved in many cellular processes and is an epigenetic mechanism primarily associated with transcriptional repression. The recent discovery that 5 mC can be oxidized to 5-hydromethylcytosine (5hmC) by TET proteins has revealed the “sixth base” of DNA and provides additional complexity to what was originally thought to be a stable repressive mark. However, our knowledge of the genome-wide distribution of 5hmC in different tissues is currently limited. Here, we sought to define loci enriched for 5hmC in the placenta genome by combining oxidative bisulphite (oxBS) treatment with high-density Illumina Infinium HumanMethylation450 methylation arrays and to compare our results with those obtained in brain. Despite identifying over 17,000 high-confidence CpG sites with consistent 5hmC enrichment, the distribution of this modification in placenta is relatively sparse when compared to cerebellum and frontal cortex. Supported by validation using allelic T4 β-glucosyltransferase assays we identify 5hmC at numerous imprinted loci, often overlapping regions associated with parent-of-origin allelic 5 mC in both placenta and brain samples. Furthermore, we observe tissue-specific monoallelic enrichment of 5hmC overlapping large clusters of imprinted snoRNAs-miRNAs processed from long noncoding RNAs (lncRNAs) within the DLK1-DIO3 cluster on chromosome 14 and SNRPN-UBE3A domain on chromosome 15. Enrichment is observed solely on the transcribed alleles suggesting 5hmC is positively associated with transcription at these loci. Our study provides an extensive description of the 5hmC/5 mC landscape in placenta with our data available at www.humanimprints.net, which represents the most comprehensive resource for exploring the epigenetic profiles associated with human imprinted genes.

Environmental enrichment modulates 5-hydroxymethylcytosine dynamics in hippocampus

Gene–environment interactions mediated at the epigenetic level may provide an initial step in delivering an appropriate response to environmental changes. 5-Hydroxymethylcytosine (5hmC), a DNA base derived from 5-methylcytosine (5mC), accounts for ~40% of modified cytosine in the brain and has been implicated in DNA methylation-related plasticity. To identify the role of 5hmC in gene–environment interactions, we exposed both young (6-week-old) and aged (18-month-old) mice to both an enriched environment and a standard environment. Exposure to EE significantly improves learning and memory in aged mice and reduces 5hmC abundance in mouse hippocampus. Furthermore, we mapped the genome-wide distribution of 5hmC and found that the alteration of 5hmC modification occurred mainly at gene bodies. In particular, genes involved in axon guidance are enriched among the genes with altered 5hmC modification. These results together suggest that environmental enrichment could modulate the dynamics of 5hmC in hippocampus, which could potentially contribute to improved learning and memory in aged animals.

Genome-wide antagonism between 5-hydroxymethylcytosine and DNA methylation in the adult mouse brain.

Mounting evidence points to critical roles for DNA modifications, including 5-methylcytosine (5mC) and its oxidized forms, in the development, plasticity and disorders of the mammalian nervous system. The novel DNA base 5-hydroxymethylcytosine (5hmC) is known to be capable of initiating passive or active DNA demethylation, but whether and how extensively 5hmC functions in shaping the post-mitotic neuronal DNA methylome is unclear. Here we report the genome-wide distribution of 5hmC in dentate granule neurons from adult mouse hippocampus in vivo. 5hmC in the neuronal genome is highly enriched in gene bodies, especially in exons, and correlates with gene expression. Direct genome-wide comparison of 5hmC distribution between embryonic stem cells and neurons reveals extensive differences, reflecting the functional disparity between these two cell types. Importantly, integrative analysis of 5hmC, overall DNA methylation and gene expression profiles of dentate granule neurons in vivo reveals the genome-wide antagonism between these two states of cytosine modifications, supporting a role for 5hmC in shaping the neuronal DNA methylome by promoting active DNA demethylation.

Large Conserved Domains Of Low DNA Methylation Maintained By 5-Hydroxymethycytosine and Dnmt3a

Identification of recurrent leukemia-associated mutations in genes encoding regulators of DNA methylation such as DNMT3A and TET2 have underscored the critical importance of DNA methylation in maintenance of normal physiology. To gain insight into how DNA methylation exerts the central role, we sought to determine the genome-wide pattern of DNA methylation in the normal precursors of leukemia cells: the hematopoietic stem cell (HSC), and investigate the factors that affect alterations in DNA methylation and gene expression.We performed whole genome bisulfite sequencing (WGBS) on purified murine HSCs achieving a total of 1,121M reads, resulting in a combined average of 40X coverage. Using Hidden Markov Model we identified 32,325 under-methylated regions (UMRs) with average proportion of methylation ≤ 10% and by inspecting the UMR size distribution, we discovered exceptionally large “methylation Canyons” which span highly conserved domains frequently containing transcription factors and are quite distinct from CpG islands and shores. Methylation Canyons are a distinct genomic feature that is stable, albeit with subtle differences, across cell-types and species. Canyon-associated genes showed a striking pattern of enrichment for genes involved in transcriptional regulation (318 genes, P=6.2 x 10-123), as well as genes containing a homeobox domain (111 genes, P=3.9 x 10-85). We compared Canyons with TF binding sites as identified from more than 150 ChIP-seq data sets across a variety of blood lineages (>10)19 and found that TF binding peaks for 10 HSC pluripotency TFs are significantly enriched in entirety of Canyons compared with their surrounding regions. Low DNA methylation is usually associated with active gene expression. However, half of Canyon genes associated with H3K27me3 showed low or no expression regardless of their H3K4me3 association while H3K4me3-only Canyon genes were highly expressed.Because DNMT3A is mutated in a high frequency of human leukemias24, we examined the impact of loss of Dnmt3a on Canyon size. Upon knockout of Dnmt3a, the edges of the Canyons are hotspots of differential methylation while regions inside of Canyon are relatively resistant. The methylation loss in Dnmt3a KO HSCs led Canyon edge erosion, Canyon size expansion and addition of 861 new Canyons for a total of 1787 Canyons. Canyons marked with H3K4me3 only were most likely to expand after Dnmt3a KO and the canyons marked only with H3K27me3 or with both marks were more likely to contract. This suggests Dnmt3a specifically is acting to restrain Canyon size where active histone marks (and active transcription) are already present. WGBS cannot distinguish between 5mC and 5hmC, so we determined the genome-wide distribution of 5hmC in WT and Dnmt3a KO HSCs using the cytosine-5-methylenesulphonate (CMS)-Seq method in which sodium bisulfate treatment convert 5hmC to CMS; CMS-containing DNA fragments are then immunoprecipitated using a CMS specific antiserum. Strikingly, 5hmC peaks were enriched specifically at the borders of Canyons. In particular, expanding Canyons, typically associated with highest H3K4me3 marking, were highly enriched at the edges for the 5hmC signal suggesting a model in which Tet proteins and Dnmt3a act concomitantly on Canyon borders opposing each other in alternately effacing and restoring methylation at the edges, particularly at sites of active chromatin marks. Using Oncomine data, we tested whether Canyon-associated genes were likely to be associated with hematologic malignancy development and found Canyon genes were highly enriched in seven signatures of genes over-expressed in Leukemia patients compared to normal bone marrow; in contrast, four sets of control genes were not similarly enriched. Further using TCGA data, we found that expressed canyon genes are significantly enriched for differentially expressed genes between patients with and without DNMT3A mutation (p value<0.05) Overall, 76 expressed canyon genes, including multiple HOX genes, are significantly changed in patients with DNMT3A mutation (p=0.0031). Methylation Canyons, the novel epigenetic landscape we describe may provide a mechanism for the regulation of hematopoiesis and may contribute to leukemia development. Disclosures:No relevant conflicts of interest to declare.

Genome-wide mapping of 5-hydroxymethylcytosine in embryonic stem cells

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.

Bioorthogonal Labeling of 5-Hydroxymethylcytosine in Genomic DNA and Diazirine-Based DNA Photo-Cross-Linking Probes

DNA is not merely a combination of four genetic codes, namely A, T, C, and G. It also contains minor modifications that play crucial roles throughout biology. For example, the fifth DNA base, 5-methylcytosine (5-mC), which accounts for ∼1% of all the nucleotides in mammalian genomic DNA, is a vital epigenetic mark. It impacts a broad range of biological functions, from development to cancer. Recently, an oxidized form of 5-methylcytosine, 5-hydroxymethylcytosine (5-hmC), was found to constitute the sixth base in the mammalian genome; it was believed to be another crucial epigenetic mark. Unfortunately, further study of this newly discovered DNA base modification has been hampered by inadequate detection and sequencing methods, because current techniques fail to differentiate 5-hmC from 5-mC. The immediate challenge, therefore, is to develop robust methods for ascertaining the positions of 5-hmC within the mammalian genome. In this Account, we describe our development of the first bioorthogonal, selective labeling of 5-hmC to specifically address this challenge. We utilize β-glucosyltransferase (βGT) to transfer an azide-modified glucose onto 5-hmC in genomic DNA. The azide moiety enables further bioorthogonal click chemistry to install a biotin group, which allows for detection, affinity enrichment, and, most importantly, deep sequencing of the 5-hmC-containing DNA. With this highly effective and selective method, we revealed the first genome-wide distribution of 5-hmC in the mouse genome and began to shed further light on the biology of 5-hmC. The strategy lays the foundation for developing high-throughput, single-base-resolution sequencing methods for 5-hmC in mammalian genomes in the future. DNA and RNA are not static inside cells. They interact with protein and other DNA and RNA in fundamental biological processes such as replication, transcription, translation, and DNA and RNA modification and repair. The ability to investigate these interactions will also be enhanced by developing and utilizing bioorthogonal probes. We have chosen the photoreactive diazirine photophore as a bioorthogonal moiety to develop nucleic acid probes. The small size and unique photo-cross-linking activity of diazirine enabled us to develop a series of novel cross-linking probes to streamline the study of protein-nucleic acid and nucleic acid-nucleic acid interactions. In the second half of this Account, we highlight a few examples of these probes.

Selective chemical labeling reveals the genome-wide distribution of 5-hydroxymethylcytosine

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.