Genome-wide screen for differential DNA methylation associated with neural cell differentiation in mouse

Cellular differentiation involves widespread epigenetic reprogramming, including modulation of DNA methylation patterns. Using Differential Methylation Hybridization (DMH) in combination with a custom DMH array containing 51,243 features covering more than 16,000 murine genes, we carried out a genome-wide screen for cell- and tissue-specific differentially methylated regions (tDMRs) in undifferentiated embryonic stem cells (ESCs), in in-vitro induced neural stem cells (NSCs) and 8 differentiated embryonic and adult tissues. Unsupervised clustering of the generated data showed distinct cell- and tissue-specific DNA methylation profiles, revealing 202 significant tDMRs (p<0.005) between ESCs and NSCs and a further 380 tDMRs (p<0.05) between NSCs/ESCs and embryonic brain tissue. We validated these tDMRs using direct bisulfite sequencing (DBS) and methylated DNA immunoprecipitation on chip (MeDIP-chip). Gene ontology (GO) analysis of the genes associated with these tDMRs showed significant (absolute Z score>1.96) enrichment for genes involved in neural differentiation, including, for example, Jag1 and Tcf4. Our results provide robust evidence for the relevance of DNA methylation in early neural development and identify novel marker candidates for neural cell differentiation.     

Sequential changes in genome-wide DNA methylation status during adipocyte differentiation

DNA methylation is an epigenetic mark on the mammalian genome. There are numerous tissue-dependent and differentially methylated regions (T-DMRs) in the unique sequences distributed throughout the genome. To determine the epigenetic changes during adipocyte differentiation, we investigated the sequential changes in DNA methylation status of 3T3-L1 cells at the growing, confluent, postconfluent and mature adipocyte cell stages. Treatment of 3T3-L1 cells with 5-aza-2′-deoxycytidine inhibited differentiation in a stage-dependent manner, supporting the idea that formation of accurate DNA methylation profile, consisting of methylated and unmethylated T-DMRs, may be involved in differentiation. Analysis by methylation-sensitive quantitative real-time PCR of the 65 known T-DMRs which contain NotI sites detected 8 methylations that changed during differentiation, and the changes in the patterns of these methylations were diverse, confirming that the differentiation process involves epigenetic alteration at the T-DMRs. Intriguingly, the dynamics of the methylation change vary depending on the T-DMRs and differentiation stages. Restriction landmark genomic scanning detected 32 novel T-DMRs, demonstrating that differentiation of 3T3-L1 cells involves genome-wide epigenetic changes by temporal methylation/demethylation, in addition to maintenance of a static methylated/demethylated state, and both depend on differentiation stage.     

Methylation profiling identified novel differentially methylated markers including OPCML and FLRT2 in prostate cancer

To develop new methods to distinguish indolent from aggressive prostate cancers (PCa), we utilized comprehensive high-throughput array-based relative methylation (CHARM) assay to identify differentially methylated regions (DMRs) throughout the genome, including both CpG island (CGI) and non-CGI regions in PCa patients based on Gleason grade. Initially, 26 samples, including 8 each of low [Gleason score (GS) 6] and high (GS ≥7) grade PCa samples and 10 matched normal prostate tissues, were analyzed as a discovery cohort. We identified 3,567 DMRs between normal and cancer tissues, and 913 DMRs distinguishing low from high-grade cancers. Most of these DMRs were located at CGI shores. The top 5 candidate DMRs from the low vs. high Gleason comparison, including OPCML, ELAVL2, EXT1, IRX5, and FLRT2, were validated by pyrosequencing using the discovery cohort. OPCML and FLRT2 were further validated in an independent cohort consisting of 20 low-Gleason and 33 high-Gleason tissues. We then compared patients with biochemical recurrence (n=70) vs. those without (n=86) in a third cohort, and they showed no difference in methylation at these DMR loci. When GS 3+4 cases and GS 4+3 cases were compared, OPCML-DMR methylation showed a trend of lower methylation in the recurrence group (n=30) than in the no-recurrence (n=52) group. We conclude that whole-genome methylation profiling with CHARM revealed distinct patterns of differential DNA methylation between normal prostate and PCa tissues, as well as between different risk groups of PCa as defined by Gleason scores. A panel of selected DMRs may serve as novel surrogate biomarkers for Gleason score in PCa.     

Analysis of tissue-specific differentially methylated regions (TDMs) in humans

Alterations in DNA methylation have been implicated in mammalian development. Hence, the identification of tissue-specific differentially methylated regions (TDMs) is indispensable for understanding its role. Using restriction landmark genomic scanning of six mouse tissues, 150 putative TDMs were identified and 14 were further analyzed. The DNA sequences of the 14 mouse TDMs are analyzed in this study. Six of the human homologous regions show TDMs to both mouse and human and genes in five of these regions have conserved tissue-specific expression: preferential expression in testis. A TDM, DDX4, is further analyzed in nine testis tissues. An increase in methylation of the promoter region is significantly associated with a marked reduction of the gene expression and defects in spermatogenesis, suggesting that hypomethylation of the DDX4 promoter region regulates DDX4 gene expression in spermatogenic cells. Our results indicate that some genomic regions with tissue-specific methylation and expression are conserved between mouse and human and suggest that DNA methylation may have an important role in regulating differentiation and tissue-/cell-specific gene expression of some genes.     

Comprehensive DNA methylation analysis of human peripheral blood leukocytes and lymphoblastoid cell lines

DNA methylation is involved in development and in human diseases. Genomic DNA derived from lymphoblastoid cell lines (LCLs) is commonly used to study DNA methylation. There are potential confounding factors regarding the use of LCL-derived DNA, however, such as Epstein-Barr (EB) viral infection and artifacts induced during cell culture. Recently, several groups compared the DNA methylation status of peripheral blood leukocytes (PBLs) and LCLs and concluded that the DNA methylation profiles between them might be consistent. To confirm and extend theses results, we performed a comprehensive DNA methylation analysis using both PBLs and LCLs derived from the same individuals. Using the luminometric methylation assay, we revealed that the global DNA methylation level was different between PBLs and LCLs. Furthermore, the direction of change was not consistent. Comparisons of genome-wide DNA methylation patterns of promoter regions revealed that methylation profiles were largely conserved between PBLs and LCLs. A preliminary analysis in a small number of samples suggested that the methylation status of an LCL may be better correlated with PBLs from the same individual than with LCLs from other individuals. Expectedly, DNA methylation in promoter regions overlapping with CpG islands was associated with gene silencing in both PBLs and LCLs. With regard to methylation differences, we found that hypermethylation was more predominant than hypomethylation in LCLs compared with PBLs. These findings suggest that LCLs should be used for DNA methylation studies with caution as the methylation patterns of promoter regions in LCLs are not always the same as those in PBLs.     

Tissue specific differentially methylated regions (TDMR): Changes in DNA methylation during development

Tissue specific differentially methylated regions (TDMRs) were identified and localized in the mouse genome using second generation virtual RLGS (vRLGS). Sequenom MassARRAY quantitative methylation analysis was used to confirm and determine the fine structure of tissue specific differences in DNA methylation. TDMRs have a broad distribution of locations to intragenic and intergenic regions including both CpG islands, and non-CpG islands regions. Somewhat surprising, there is a strong bias for TDMR location in non-promoter intragenic regions. Although some TDMRs are within or close to repeat sequences, overall they are less frequently associated with repetitive elements than expected from a random distribution. Many TDMRs are methylated at early developmental stages, but unmethylated later, suggesting active or passive demethylation, or expansions of populations of cells with unmethylated TDMRs. This is notable during postnatal testis differentiation where many testis specific TDMRs become progressively "demethylated". These results suggest that methylation changes during development are dynamic, involve demethylation and methylation, and may occur at late stages of embryonic development or even postnatally.     

Bacterial DNA containing methylated CpG motifs retains immunostimulatory activity in synergy with modified lipopolysaccharides

We previously described the immunostimulatory activity of CIA07, a combination of bacterial DNA fragments and modified LPS, and demonstrated that CIA07 has antitumor activity in a mouse bladder cancer model. In this study, we investigated whether methylation of the CpG motifs on the bacterial DNA fragments affects the immunostimulatory potential of CIA07. E. coli DNA fragments were methylated with CpG methylase, and then combined with modified LPS for experiments. Our results revealed that methylated CIA07 (mCIA07) and unmethylated CIA07 were equally active in inducing cytokine secretion from human whole blood cells. In addition, both methylated DNA fragments and mCIA07 retained the ability to activate expression and nuclear translocation of NF-kappaB in RAW 264.7 cells. Finally, methylated DNA fragments and mCIA07 exhibited an antitumor activity comparable to those of their unmethylated counterparts in our mouse bladder cancer model. These data demonstrate that CpG methylation of E. coli DNA does not abrogate the immunostimulatory activity of DNA fragments or CIA07, suggesting that the synergistic activity by bacterial DNA in combination with LPS may be independent of the methylation status of CpG motifs.     

Inter-individual variation of DNA methylation and its implications for large-scale epigenome mapping

Genomic DNA methylation profiles exhibit substantial variation within the human population, with important functional implications for gene regulation. So far little is known about the characteristics and determinants of DNA methylation variation among healthy individuals. We performed bioinformatic analysis of high-resolution methylation profiles from multiple individuals, uncovering complex patterns of inter-individual variation that are strongly correlated with the local DNA sequence. CpG-rich regions exhibit low and relatively similar levels of DNA methylation in all individuals, but the sequential order of the (few) methylated among the (many) unmethylated CpGs differs randomly across individuals. In contrast, CpG-poor regions exhibit substantially elevated levels of inter-individual variation, but also significant conservation of specific DNA methylation patterns between unrelated individuals. This observation has important implications for experimental analysis of DNA methylation, e.g. in the context of epigenome projects. First, DNA methylation mapping at single-CpG resolution is expected to uncover informative DNA methylation patterns for the CpG-poor bulk of the human genome. Second, for CpG-rich regions it will be sufficient to measure average methylation levels rather than assaying every single CpG. We substantiate these conclusions by an in silico benchmarking study of six widely used methods for DNA methylation mapping. Based on our findings, we propose a cost-optimized two-track strategy for mammalian methylome projects.     

Buccals are likely to be a more informative surrogate tissue than blood for epigenome-wide association studies

There is increasing evidence that interindividual epigenetic variation is an etiological factor in common human diseases. Such epigenetic variation could be genetic or non-genetic in origin, and epigenome-wide association studies (EWASs) are underway for a wide variety of diseases/phenotypes. However, performing an EWAS is associated with a range of issues not typically encountered in genome-wide association studies (GWASs), such as the tissue to be analyzed. In many EWASs, it is not possible to analyze the target tissue in large numbers of live humans, and consequently surrogate tissues are employed, most commonly blood. But there is as yet no evidence demonstrating that blood is more informative than buccal cells, the other easily accessible tissue. To assess the potential of buccal cells for use in EWASs, we performed a comprehensive analysis of a buccal cell methylome using whole-genome bisulfite sequencing. Strikingly, a buccal vs. blood comparison reveals > 6X as many hypomethylated regions in buccal. These tissue-specific differentially methylated regions (tDMRs) are strongly enriched for DNaseI hotspots. Almost 75% of these tDMRs are not captured by commonly used DNA methylome profiling platforms such as Reduced Representational Bisulfite Sequencing and the Illumina Infinium HumanMethylation450 BeadChip, and they also display distinct genomic properties. Buccal hypo-tDMRs show a statistically significant enrichment near SNPs associated to disease identified through GWASs. Finally, we find that, compared with blood, buccal hypo-tDMRs show significantly greater overlap with hypomethylated regions in other tissues. We propose that for non-blood based diseases/phenotypes, buccal will be a more informative tissue for EWASs.     

ABCC9, NKAPL,and TMEM132C are potential diagnostic and prognostic markers in triple‐negative breast cancer

Triple‐negative breast cancer (TNBC) is a highly heterogeneous disease. The aim of this study is to identify the diagnostic and poor prognostic signatures in TNBC by exploring the aberrant DNA methylation and gene expression. Differential expression and methylation analysis of the TNBC and paracancer samples from The Cancer Genome Atlas were performed. Gene set enrichment and protein–protein interaction (PPI) network analysis was used to explore the mechanisms of TNBC. Methylation‐gene expression correlation analysis was performed, and multivariate Cox analysis and receiver operating characteristics analysis were used to further screen the hub genes for TNBC. We identified 1,525 differentially expressed genes and 150 differentially methylated genes between TNBC and paracancer samples. About 96.64% of the methylation sites were located on the CpG island. A total of 17 Gene Ontology biological process terms and 18 signal pathways were significantly enriched. GNG4, GNG11, PENK, MAOA, and AOX1 were identified as the core genes of the PPI network. Methylation‐expression correlations revealed that ABCC9 (cg06951626), NKAPL (cg18675097, cg01031101, and cg17384889), and TMEM132C (cg03530754) showed promise as diagnostic and prognostic markers in TNBC. ABCC9 (cg06951626), NKAPL (cg18675097, cg01031101, and cg17384889), and TMEM132C (cg03530754) were potential diagnostic and prognostic markers in TNBC.     

Comprehensive DNA methylation analysis of human peripheral blood leukocytes and lymphoblastoid cell lines

DNA methylation is involved in development and in human diseases. Genomic DNA derived from lymphoblastoid cell lines (LCLs) is commonly used to study DNA methylation. There are potential confounding factors regarding the use of LCL-derived DNA, however, such as Epstein-Barr (EB) viral infection and artifacts induced during cell culture. Recently, several groups compared the DNA methylation status of peripheral blood leukocytes (PBLs) and LCLs and concluded that the DNA methylation profiles between them might be consistent. To confirm and extend theses results, we performed a comprehensive DNA methylation analysis using both PBLs and LCLs derived from the same individuals. Using the luminometric methylation assay, we revealed that the global DNA methylation level was different between PBLs and LCLs. Furthermore, the direction of change was not consistent. Comparisons of genome-wide DNA methylation patterns of promoter regions revealed that methylation profiles were largely conserved between PBLs and LCLs. A preliminary analysis in a small number of samples suggested that the methylation status of an LCL may be better correlated with PBLs from the same individual than with LCLs from other individuals. Expectedly, DNA methylation in promoter regions overlapping with CpG islands was associated with gene silencing in both PBLs and LCLs. With regard to methylation differences, we found that hypermethylation was more predominant than hypomethylation in LCLs compared with PBLs. These findings suggest that LCLs should be used for DNA methylation studies with caution as the methylation patterns of promoter regions in LCLs are not always the same as those in PBLs.Key words: DNA methylation, lymphoblastoid cell lines, peripheral blood leukocytes, LUMA, promoter tiling array, gene expression     

CpG_MPs: identification of CpG methylation patterns of genomic regions from high-throughput bisulfite sequencing data

High-throughput bisulfite sequencing is widely used to measure cytosine methylation at single-base resolution in eukaryotes. It permits systems-level analysis of genomic methylation patterns associated with gene expression and chromatin structure. However, methods for large-scale identification of methylation patterns from bisulfite sequencing are lacking. We developed a comprehensive tool, CpG_MPs, for identification and analysis of the methylation patterns of genomic regions from bisulfite sequencing data. CpG_MPs first normalizes bisulfite sequencing reads into methylation level of CpGs. Then it identifies unmethylated and methylated regions using the methylation status of neighboring CpGs by hotspot extension algorithm without knowledge of pre-defined regions. Furthermore, the conservatively and differentially methylated regions across paired or multiple samples (cells or tissues) are identified by combining a combinatorial algorithm with Shannon entropy. CpG_MPs identified large amounts of genomic regions with different methylation patterns across five human bisulfite sequencing data during cellular differentiation. Different sequence features and significantly cell-specific methylation patterns were observed. These potentially functional regions form candidate regions for functional analysis of DNA methylation during cellular differentiation. CpG_MPs is the first user-friendly tool for identifying methylation patterns of genomic regions from bisulfite sequencing data, permitting further investigation of the biological functions of genome-scale methylation patterns.     

Genome-Wide Analysis of DNA Methylation Dynamics during Early Human Development

DNA methylation is globally reprogrammed during mammalian preimplantation development, which is critical for normal development. Recent reduced representation bisulfite sequencing (RRBS) studies suggest that the methylome dynamics are essentially conserved between human and mouse early embryos. RRBS is known to cover 5–10% of all genomic CpGs, favoring those contained within CpG-rich regions. To obtain an unbiased and more complete representation of the methylome during early human development, we performed whole genome bisulfite sequencing of human gametes and blastocysts that covered>70% of all genomic CpGs. We found that the maternal genome was demethylated to a much lesser extent in human blastocysts than in mouse blastocysts, which could contribute to an increased number of imprinted differentially methylated regions in the human genome. Global demethylation of the paternal genome was confirmed, but SINE-VNTR-Alu elements and some other tandem repeat-containing regions were found to be specifically protected from this global demethylation. Furthermore, centromeric satellite repeats were hypermethylated in human oocytes but not in mouse oocytes, which might be explained by differential expression of de novo DNA methyltransferases. These data highlight both conserved and species-specific regulation of DNA methylation during early mammalian development. Our work provides further information critical for understanding the epigenetic processes underlying differentiation and pluripotency during early human development.