Integrated analyses of DNA methylation and hydroxymethylation reveal tumor suppressive roles of ECM1, ATF5, and EOMES in human hepatocellular carcinoma

Background

Differences in 5-hydroxymethylcytosine, 5hmC, distributions may complicate previous observations of abnormal cytosine methylation statuses that are used for the identification of new tumor suppressor gene candidates that are relevant to human hepatocarcinogenesis. The simultaneous detection of 5-methylcytosine and 5-hydroxymethylcytosine is likely to stimulate the discovery of aberrantly methylated genes with increased accuracy in human hepatocellular carcinoma.

Results

Here, we performed ultra-performance liquid chromatography/tandem mass spectrometry and single-base high-throughput sequencing, Hydroxymethylation and Methylation Sensitive Tag sequencing, HMST-seq, to synchronously measure these two modifications in human hepatocellular carcinoma samples. After identification of differentially methylated and hydroxymethylated genes in human hepatocellular carcinoma, we integrate DNA copy-number alterations, as determined using array-based comparative genomic hybridization data, with gene expression to identify genes that are potentially silenced by promoter hypermethylation.

Conclusions

We report a high enrichment of genes with epigenetic aberrations in cancer signaling pathways. Six genes were selected as tumor suppressor gene candidates, among which, ECM1, ATF5 and EOMES are confirmed via siRNA experiments to have potential anti-cancer functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0533-9) contains supplementary material, which is available to authorized users.     

Ontogeny,distribution and potential roles of 5-hydroxymethylcytosine in human liver function

Background

Interindividual differences in liver functions such as protein synthesis, lipid and carbohydrate metabolism and drug metabolism are influenced by epigenetic factors. The role of the epigenetic machinery in such processes has, however, been barely investigated. 5-hydroxymethylcytosine (5hmC) is a recently re-discovered epigenetic DNA modification that plays an important role in the control of gene expression.

Results

In this study, we investigate 5hmC occurrence and genomic distribution in 8 fetal and 7 adult human liver samples in relation to ontogeny and function. LC-MS analysis shows that in the adult liver samples 5hmC comprises up to 1% of the total cytosine content, whereas in all fetal livers it is below 0.125%. Immunohistostaining of liver sections with a polyclonal anti-5hmC antibody shows that 5hmC is detected in most of the hepatocytes. Genome-wide mapping of the distribution of 5hmC in human liver samples by next-generation sequencing shows significant differences between fetal and adult livers. In adult livers, 5hmC occupancy is overrepresented in genes involved in active catabolic and metabolic processes, whereas 5hmC elements which are found in genes exclusively in fetal livers and disappear in the adult state, are more specific to pathways for differentiation and development.

Conclusions

Our findings suggest that 5-hydroxymethylcytosine plays an important role in the development and function of the human liver and might be an important determinant for development of liver diseases as well as of the interindividual differences in drug metabolism and toxicity.     

Comparative characterization of the PvuRts1I family of restriction enzymes and their application in mapping genomic 5-hydroxymethylcytosine

PvuRts1I is a modification-dependent restriction endonuclease that recognizes 5-hydroxymethylcytosine (5hmC) as well as 5-glucosylhydroxymethylcytosine (5ghmC) in double-stranded DNA. Using PvuRts1I as the founding member, we define a family of homologous proteins with similar DNA modification-dependent recognition properties. At the sequence level, these proteins share a few uniquely conserved features. We show that these enzymes introduce a double-stranded cleavage at the 3'-side away from the recognized modified cytosine. The distances between the cleavage sites and the modified cytosine are fixed within a narrow range, with the majority being 11-13 nt away in the top strand and 9-10 nt away in the bottom strand. The recognition sites of these enzymes generally require two cytosines on opposite strand around the cleavage sites, i.e. 5'-CN(11-13)↓N(9-10)G-3'/3'-GN(9-10)↓N(11-13)C-5', with at least one cytosine being modified for efficient cleavage. As one potential application for these enzymes is to provide useful tools for selectively mapping 5hmC sites, we have compared the relative selectivity of a few PvuRts1I family members towards different forms of modified cytosines. Our results show that the inherently different relative selectivity towards modified cytosines can have practical implications for their application. By using AbaSDFI, a PvuRts1I homolog with the highest relative selectivity towards 5ghmC, to analyze rat brain DNA, we show it is feasible to map genomic 5hmC sites close to base resolution. Our study offers unique tools for determining more accurate hydroxymethylomes in mammalian cells.     

Simple and accurate single base resolution analysis of 5-hydroxymethylcytosine by catalytic oxidative bisulfite sequencing using micelle incarcerated oxidants

Oxidation of 5-methylcytosine (5mC) is catalyzed by ten-eleven translocation (TET) enzymes to produce 5-hydroxymethylcytosine (5hmC) and following oxidative products. The oxidized nucleotides were shown to be the intermediates for DNA demethylation, as the nucleotides are removed by base excision repair system initiated by thymine DNA glycosylase. A simple and accurate method to determine initial oxidation product 5hmC at single base resolution in genomic DNA is necessary to understand demethylation mechanism. Recently, we have developed a new catalytic oxidation reaction using micelle-incarcerated oxidants to oxidize 5hmC to form 5-formylcytosine (5fC), and subsequent bisulfite sequencing can determine the positions of 5hmC in DNA. In the present study, we described the optimization of the catalytic oxidative bisulfite sequencing (coBS-seq), and its application to the analysis of 5hmC in genomic DNA at single base resolution in a quantitative manner. As the oxidation step showed quite low damage on genomic DNA, the method allows us to down scale the sample to be analyzed.     

A novel method for the efficient and selective identification of 5-hydroxymethylcytosine in genomic DNA

Recently, 5-hydroxymethylcytosine (5hmC) was identified in mammalian genomic DNA. The biological role of this modification remains unclear; however, identifying the genomic location of this modified base will assist in elucidating its function. We describe a method for the rapid and inexpensive identification of genomic regions containing 5hmC. This method involves the selective glucosylation of 5hmC residues by the β-glucosyltransferase from T4 bacteriophage creating β-glucosyl-5-hydroxymethylcytosine (β-glu-5hmC). The β-glu-5hmC modification provides a target that can be efficiently and selectively pulled down by J-binding protein 1 coupled to magnetic beads. DNA that is precipitated is suitable for analysis by quantitative PCR, microarray or sequencing. Furthermore, we demonstrate that the J-binding protein 1 pull down assay identifies 5hmC at the promoters of developmentally regulated genes in human embryonic stem cells. The method described here will allow for a greater understanding of the temporal and spatial effects that 5hmC may have on epigenetic regulation at the single gene level.     

TET蛋白：一个新的DNA修饰酶家族

DNA的胞嘧啶（C）5-甲基化是一种重要的表观修饰,它参与基因调节、基因组印记、X-染色体失活、重复序列抑制和癌症发生等过程. 5-甲基胞嘧啶（5mC）可被TET (ten-eleven translocation)蛋白家族进一步转化为5-羟甲基胞嘧啶（5hmC）,该过程是DNA去甲基化的1个必要阶段. 5hmC可在活性转录基因起始位点和Polycomb抑制基因启动子延伸区域富集.TET蛋白包括3个成员TET1、TET2和TET3,均属于α-酮戊二酸和Fe²⁺依赖的双加氧酶,其催化涉及氧化过程.小鼠Tet1在胚胎干细胞发育中拥有双重作用,即促进全能因子的转录,又参与发育调节因子的抑制.人TET蛋白的破坏与造血系统肿瘤相关,如在骨髓增生性疾病/肿瘤存在频繁的TET2基因突变.TET蛋白和5hmC的研究为DNA甲基化/去甲基化及其生物学功能提供了新的视点.     

Single-base resolution analysis of DNA epigenome via high-throughput sequencing

Epigenetic changes caused by DNA methylation and histone modifications play important roles in the regulation of various cellular processes and development. Recent discoveries of 5-methylcytosine (5mC) oxidation derivatives including 5-hydroxymethylcytosine (5hmC), 5-formylcytsine (5fC) and 5-carboxycytosine (5caC) in mammalian genome further expand our understanding of the epigenetic regulation. Analysis of DNA modification patterns relies increasingly on sequencing-based profiling methods. A number of different approaches have been established to map the DNA epigenomes with single-base resolution, as represented by the bisulfite-based methods, such as classical bisulfite sequencing (BS-seq), TAB-seq (TET-assisted bisulfite sequencing), oxBS-seq (oxidative bisulfite sequencing) and etc. These methods have been used to generate base-resolution maps of 5mC and its oxidation derivatives in genomic samples. The focus of this review will be to discuss the chemical methodologies that have been developed to detect the cytosine derivatives in the genomic DNA.     

Emerging roles of TET proteins and 5-Hydroxymethylcytosines in active DNA demethylation and beyond

Cytosine methylation is the major epigenetic modification of metazoan DNA. Although there is strong evidence that active DNA demethylation occurs in animal cells, the molecular details of this process are unknown. The recent discovery of the TET protein family (TET1–3) 5-methylcytosine hydroxylases has provided a new entry point to reveal the identity of the long-sought DNA demethylase. Here, we review the recent progress in understanding the function of TET proteins and 5-hydroxymethylcytosine (5hmC) through various biochemical and genomic approaches, the current evidence for a role of 5hmC as an early intermediate in active DNA demethylation and the potential functions of TET proteins and 5hmC beyond active DNA demethylation. We also discuss how future studies can extend our knowledge of this novel epigenetic modification.Key words: TET1, 5-hydroxymethylcytosine, active DNA demethylation, epigenetic, DNA methylation, hippocampus, electroconvulsive stimulation, Gadd45b, BER     

Lead exposure induces changes in 5-hydroxymethylcytosine clusters in CpG islands in human embryonic stem cells and umbilical cord blood

Prenatal exposure to neurotoxicants such as lead (Pb) may cause stable changes in the DNA methylation (5mC) profile of the fetal genome. However, few studies have examined its effect on the DNA de-methylation pathway, specifically the dynamic changes of the 5-hydroxymethylcytosine (5hmC) profile. Therefore, in this study, we investigate the relationship between Pb exposure and 5mC and 5hmC modifications during early development. To study the changes in the 5hmC profile, we use a novel modification of the Infinium™ HumanMethylation450 assay (Illumina, Inc.), which we named HMeDIP-450K assay, in an in vitro human embryonic stem cell model of Pb exposure. We model Pb exposure-associated 5hmC changes as clusters of correlated, adjacent CpG sites, which are co-responding to Pb. We further extend our study to look at Pb-dependent changes in high density 5hmC regions in umbilical cord blood DNA from 48 mother-infant pairs from the Early Life Exposure in Mexico to Environmental Toxicants (ELEMENT) cohort. For our study, we randomly selected umbilical cord blood from 24 male and 24 female children from the 1st and 4th quartiles of Pb levels. Our data show that Pb-associated changes in the 5hmC and 5mC profiles can be divided into sex-dependent and sex-independent categories. Interestingly, differential 5mC sites are better markers of Pb-associated sex-dependent changes compared to differential 5hmC sites. In this study we identified several 5hmC and 5mC genomic loci, which we believe might have some potential as early biomarkers of prenatal Pb exposure.     

Sensitive enzymatic quantification of 5-hydroxymethylcytosine in genomic DNA

The recent discovery of genomic 5-hydroxymethylcytosine (hmC) and mutations affecting the respective Tet hydroxylases in leukemia raises fundamental questions about this epigenetic modification. We present a sensitive method for fast quantification of genomic hmC based on specific transfer of radiolabeled glucose to hmC by a purified glucosyltransferase. We determined hmC levels in various adult tissues and differentiating embryonic stem cells and show a correlation with differential expression of tet genes.