Development of a Novel Output Value for Quantitative Assessment in Methylated DNA Immunoprecipitation-CpG Island Microarray Analysis

In DNA methylation microarray analysis, quantitative assessment of intermediate methylation levels in samples with various global methylation levels is still difficult. Here, specifically for methylated DNA immunoprecipitation-CpG island (CGI) microarray analysis, we developed a new output value. The signal log ratio reflected the global methylation levels, but had only moderate linear correlation (r = 0.72) with the fraction of DNA molecules immunoprecipitated. By multiplying the signal log ratio using a coefficient obtained from the probability value that took account of signals in neighbouring probes, its linearity was markedly improved (r = 0.94). The new output value, Me value, reflected the global methylation level, had a strong correlation also with the fraction of methylated CpG sites obtained by bisulphite sequencing (r = 0.88), and had an accuracy of 71.8 and 83.8% in detecting completely methylated and unmethylated CGIs. Analysis of gastric cancer cell lines using the Me value showed that methylation of CGIs in promoters and gene bodies was associated with low and high, respectively, gene expression. The degree of demethylation of promoter CGIs after 5-aza-2''-deoxycytidine treatment had no association with that of induction of gene expression. The Me value was considered to be useful for analysis of intermediate methylation levels of CGIs.     

MBD-isolated Genome Sequencing provides a high-throughput and comprehensive survey of DNA methylation in the human genome

DNA methylation is an epigenetic modification involved in both normal developmental processes and disease states through the modulation of gene expression and the maintenance of genomic organization. Conventional methods of DNA methylation analysis, such as bisulfite sequencing, methylation sensitive restriction enzyme digestion and array-based detection techniques, have major limitations that impede high-throughput genome-wide analysis. We describe a novel technique, MBD-isolated Genome Sequencing (MiGS), which combines precipitation of methylated DNA by recombinant methyl-CpG binding domain of MBD2 protein and sequencing of the isolated DNA by a massively parallel sequencer. We utilized MiGS to study three isogenic cancer cell lines with varying degrees of DNA methylation. We successfully detected previously known methylated regions in these cells and identified hundreds of novel methylated regions. This technique is highly specific and sensitive and can be applied to any biological settings to identify differentially methylated regions at the genomic scale.     

Influence of DNA methylation on positioning and DNA flexibility of nucleosomes with pericentric satellite DNA

DNA methylation occurs on CpG sites and is important to form pericentric heterochromatin domains. The satellite 2 sequence, containing seven CpG sites, is located in the pericentric region of human chromosome 1 and is highly methylated in normal cells. In contrast, the satellite 2 region is reportedly hypomethylated in cancer cells, suggesting that the methylation status may affect the chromatin structure around the pericentric regions in tumours. In this study, we mapped the nucleosome positioning on the satellite 2 sequence in vitro and found that DNA methylation modestly affects the distribution of the nucleosome positioning. The micrococcal nuclease assay revealed that the DNA end flexibility of the nucleosomes changes, depending on the DNA methylation status. However, the structures and thermal stabilities of the nucleosomes are unaffected by DNA methylation. These findings provide new information to understand how DNA methylation functions in regulating pericentric heterochromatin formation and maintenance in normal and malignant cells.     

Identification of a DNA Methylome Profile of Esophageal Squamous Cell Carcinoma and Potential Plasma Epigenetic Biomarkers for Early Diagnosis

DNA methylation is a critical epigenetic mechanism involved in key cellular processes. Its deregulation has been linked to many human cancers including esophageal squamous cell carcinoma (ESCC). This study was designed to explore the whole methylation status of ESCC and to identify potential plasma biomarkers for early diagnosis. We used Infinium Methylation 450k array to analyze ESCC tissues (n = 4), paired normal surrounding tissues (n = 4) and normal mucosa from healthy individuals (n = 4), and combined these with gene expression data from the GEO database. One hundred and sixty eight genes had differentially methylated CpG sites in their promoter region and a gene expression pattern inverse to the direction of change in DNA methylation. These genes were involved in several cancer-related pathways. Three genes were validated in additional 42 ESCC tissues and paired normal surrounding tissues. The methylation frequency of EPB41L3, GPX3, and COL14A1 were higher in tumor tissues than in normal surrounding tissues (P<0.017). The higher methylation frequency of EPB41l3 was correlated with large tumor size (P = 0.044) and advanced pT tumor stage (P = 0.001). The higher methylation frequency of GPX3 and COL14A1 were correlated with advanced pN tumor stage (P = 0.001 and P<0.001). The methylation of EPB41L3, GPX3, and COL14A1 genes were only found in ESCC patients'' plasma, but not in normal individuals upon testing 42 ESCC patients and 50 healthy individuals. Diagnostic sensitivity was increased when methylation of any of the 3 genes were counted (64.3% sensitivity and 100% specificity). These differentially methylated genes in plasma may be used as biomarkers for early diagnosis of ESCC.     

Genome-wide DNA methylation profiling of chronic lymphocytic leukemia allows identification of epigenetically repressed molecular pathways with clinical impact

We performed a genome-wide analysis of aberrant DNA methylation in chronic lymphocytic leukemia (CLL) using methylated CpG island amplification (MCA) coupled with a promoter microarray. We identified 280 potential targets of aberrant DNA methylation in CLL. These genes were located more frequently in chromosomes 19 (16%, p = 0.001), 16 (11%, p = 0.001), 17 (10%, p = 0.02) and 11 (9%, p = 0.02) and could be grouped in several functional networks. Methylation status was confirmed for 22 of these genes (SOX11, DLX1, FAM62C, SOX14, RSPO1, ADCY5, HAND2, SPOCK, MLL, ING1, PRIMA1, BCL11B, LTBP2, BNC1, NR2F2, SALL1, GALGT2, LHX1, DLX4, KLK10, TFAP2 and APP) in 78 CLL patients by pyrosequencing. As a proof of principle, we analyzed the expression of 2 genes, PRIMA1 and APP, in primary cells and of GALGT2, TFAP2C and PRIMA1 in leukemia cells. There was an inverse association between methylation and gene expression. This could be reversed by treatment with 5-aza-2′-deoxycytidine in cell lines. Treatment in a clinical trial with 5-azacitidine resulted in decreased methylation of LINE, DLX4 and SALL1 in the peripheral blood B-cells of patients with CLL. IgV_H mutational status or ZAP-70 expression were not associated with specific methylation profiles. By multivariate analysis, methylation of LINE and APP was associated with shorter overall survival (p = 0.045 and 0.0035, respectively). This study demonstrates that aberrant DNA methylation is common and has potential prognostic and therapeutic value in CLL.Key words: chronic lymphocytic leukemia, DNA methylation, MCA/promoter microarray, epigenetics     

A DNA methylation microarray-based study identifies ERG as a gene commonly methylated in prostate cancer

DNA methylation of promoter regions is a common event in prostate cancer, one of the most common cancers in men worldwide. Because prior reports demonstrating that DNA methylation is important in prostate cancer studied a limited number of genes, we systematically quantified the DNA methylation status of 1505 CpG dinucleotides for 807 genes in 78 paraffin-embedded prostate cancer samples and three normal prostate samples. The ERG gene, commonly repressed in prostate cells in the absence of an oncogenic fusion to the TMPRSS2 gene, was one of the most commonly methylated genes, occurring in 74% of prostate cancer specimens. In an independent group of patient samples, we confirmed that ERG DNA methylation was common, occurring in 57% of specimens, and cancer-specific. The ERG promoter is marked by repressive chromatin marks mediated by polycomb proteins in both normal prostate cells and prostate cancer cells, which may explain ERG’s predisposition to DNA methylation and the fact that tumors with ERG DNA methylation were more methylated, in general. These results demonstrate that bead arrays offer a high-throughput method to discover novel genes with promoter DNA methylation such as ERG, whose measurement may improve our ability to more accurately detect prostate cancer.     

HOP/OB1/NECC1 promoter DNA is frequently hypermethylated and involved in tumorigenic ability in esophageal squamous cell carcinoma

Promoter DNA hypermethylation with gene silencing is a common feature of human cancer, and cancer-prone methylation is believed to be a landmark of tumor suppressor genes (TSG). Identification of novel methylated genes would not only aid in the development of tumor markers but also elucidate the biological behavior of human cancers. We identified several epigenetically silenced candidate TSGs by pharmacologic unmasking of esophageal squamous cell carcinoma (ESCC) cell lines by demethylating agents (5-aza-2'-deoxycitidine and trichostatin A) combined with ESCC expression profiles using expression microarray. HOP/OB1/NECC1 was identified as an epigenetically silenced candidate TSG and further examined for (a) expression status, (b) methylation status, and (c) functional involvement in cancer cell lines. (a) The HOP gene encodes two putative promoters (promoters A and B) associated with two open reading frames (HOPalpha and HOPbeta, respectively), and HOPalpha and HOPbeta were both down-regulated in ESCC independently. (b) Promoter B harbors dense CpG islands, in which we found dense methylation in a cancer-prone manner (55% in tumor tissues by TaqMan methylation-specific PCR), whereas promoter A does not harbor CpG islands. HOPbeta silencing was associated with DNA methylation of promoter B in nine ESCC cell lines tested, and reactivated by optimal conditions of demethylating agents, whereas HOPalpha silencing was not reactivated by such treatments. Forced expression of HOP suppressed tumorigenesis in soft agar in four different squamous cell carcinoma cell lines. More convincingly, RNA interference knockdown of HOP in TE2 cells showed drastic restoration of the oncogenic phenotype. In conclusion, HOP is a putative TSG that harbors tumor inhibitory activity, and we for the first time showed that the final shutdown process of HOP expression is linked to promoter DNA hypermethylation under the double control of the discrete promoter regions in cancer.     

Impact of DNA demethylation of the G0S2 gene on the transcription of G0S2 in squamous lung cancer cell lines with or without nuclear receptor agonists

We recently identified that DNA methylation of the G0S2 gene was significantly more frequent in squamous lung cancer than in non-squamous lung cancer. However, the significance of G0S2 methylation levels on cancer cells is not yet known. We investigated the effect of G0S2 methylation levels on cell growth, mRNA expression, and chromatin structure using squamous lung cancer cell lines and normal human bronchial epithelial cells. DNA methylation and mRNA expression of G0S2 were inversely correlated, and in one of the squamous lung cancer cell lines, LC-1 sq, G0S2 was completely methylated and suppressed. Overexpression of G0S2 in LC-1 sq did not show growth arrest or apoptosis. The G0S2 gene has been reported to be a target gene of all-trans retinoic acid and peroxisome proliferator-activated receptor agonists. We treated LC-1 sq with 5-Aza-2′-deoxycytidine, Trichostatin A, all-trans retinoic acid, Wy 14643, or Pioglitazone either alone or in combination. Only 5-Aza-2′-deoxycytidine restored mRNA expression of G0S2. Chromatin immunoprecipitation revealed that histone H3 lysine 9 was methylated regardless of DNA methylation or mRNA expression. In summary, mRNA expression of G0S2 was regulated mainly by DNA methylation in squamous lung cancer cell lines. When the G0S2 gene was methylated, nuclear receptor agonists could not restore mRNA expression of G0S2 and did not show any additive effect on mRNA expression of G0S2 even after the treatment with 5-Aza-2′-deoxycytidine.     

Genome-wide DNA methylation profiles in noncancerous gastric mucosae with regard to Helicobacter pylori infection and the presence of gastric cancer

Background and Aims: To determine genome‐wide DNA methylation profiles induced by Helicobacter pylori (H. pylori) infection and to identify methylation markers in H. pylori‐induced gastric carcinogenesis. Methods: Gastric mucosae obtained from controls (n = 20) and patients with gastric cancer (n = 28) were included. A wide panel of CpG sites in cancer‐related genes (1505 CpG sites in 807 genes) was analyzed using Illumina bead array technology. Validation of the results of Illumina bead array technique was performed using methylation‐specific PCR method for four genes (MOS, DCC, CRK, and PTPN6). Results: The Illumina bead array showed that a total of 359 CpG sites (269 genes) were identified as differentially methylated by H. pylori infection (p < .0001). The correlation between methylation‐specific PCR and bead array analysis was significant (p < .0001, Spearman coefficient = 0.5054). Methylation profiles in noncancerous gastric mucosae of the patients with gastric cancer showed quite distinct patterns according to the presence or absence of the current H. pylori infection; however, 10 CpG sites were identified to be hypermethylated and three hypomethylated in association with the presence of gastric cancer regardless of H. pylori infection (p < .01). Conclusions: Genome‐wide methylation profiles showed a number of genes differentially methylated by H. pylori infection. Methylation profiles in noncancerous gastric mucosae from the patients with gastric cancer can be affected by H. pylori‐induced gastritis. Differentially methylated CpG sites in this study needs to be validated in a larger population using quantitative methylation‐specific PCR method.     

Hierarchical Clustering of Breast Cancer Methylomes Revealed Differentially Methylated and Expressed Breast Cancer Genes

Oncogenic transformation of normal cells often involves epigenetic alterations, including histone modification and DNA methylation. We conducted whole-genome bisulfite sequencing to determine the DNA methylomes of normal breast, fibroadenoma, invasive ductal carcinomas and MCF7. The emergence, disappearance, expansion and contraction of kilobase-sized hypomethylated regions (HMRs) and the hypomethylation of the megabase-sized partially methylated domains (PMDs) are the major forms of methylation changes observed in breast tumor samples. Hierarchical clustering of HMR revealed tumor-specific hypermethylated clusters and differential methylated enhancers specific to normal or breast cancer cell lines. Joint analysis of gene expression and DNA methylation data of normal breast and breast cancer cells identified differentially methylated and expressed genes associated with breast and/or ovarian cancers in cancer-specific HMR clusters. Furthermore, aberrant patterns of X-chromosome inactivation (XCI) was found in breast cancer cell lines as well as breast tumor samples in the TCGA BRCA (breast invasive carcinoma) dataset. They were characterized with differentially hypermethylated XIST promoter, reduced expression of XIST, and over-expression of hypomethylated X-linked genes. High expressions of these genes were significantly associated with lower survival rates in breast cancer patients. Comprehensive analysis of the normal and breast tumor methylomes suggests selective targeting of DNA methylation changes during breast cancer progression. The weak causal relationship between DNA methylation and gene expression observed in this study is evident of more complex role of DNA methylation in the regulation of gene expression in human epigenetics that deserves further investigation.     

Inheritance Patterns and Stability of DNA Methylation Variation in Maize Near-Isogenic Lines

DNA methylation is a chromatin modification that contributes to epigenetic regulation of gene expression. The inheritance patterns and trans-generational stability of 962 differentially methylated regions (DMRs) were assessed in a panel of 71 near-isogenic lines (NILs) derived from maize (Zea mays) inbred lines B73 and Mo17. The majority of DMRs exhibit inheritance patterns that would be expected for local (cis) inheritance of DNA methylation variation such that DNA methylation level was coupled to local genotype. There are few examples of DNA methylation that exhibit trans-acting control or paramutation-like patterns. The cis-inherited DMRs provide an opportunity to study the stability of inheritance for DNA methylation variation. There was very little evidence for alterations of DNA methylation levels at these DMRs during the generations of the NIL population development. DNA methylation level was associated with local genotypes in nearly all of the >30,000 potential cases of inheritance. The majority of the DMRs were not associated with small RNAs. Together, our results suggest that a significant portion of DNA methylation variation in maize exhibits locally (cis) inherited patterns, is highly stable, and does not require active programming by small RNAs for maintenance.DNA methylation may contribute to heritable epigenetic information in many eukaryotic genomes. In this study, we have documented the inheritance patterns and trans-generational stability for nearly 1000 DNA methylation variants in a segregating maize population. At most loci studied, the DNA methylation differences are locally inherited and are not influenced by the other allele or other genomic regions. The inheritance of DNA methylation levels across generations is quite robust with almost no examples of unstable inheritance, suggesting that DNA methylation differences can be quite stably inherited, even in segregating populations.     

Direct detection of methylation in genomic DNA

The identification of methylated sites on bacterial genomic DNA would be a useful tool to study the major roles of DNA methylation in prokaryotes: distinction of self and nonself DNA, direction of post-replicative mismatch repair, control of DNA replication and cell cycle, and regulation of gene expression. Three types of methylated nucleobases are known: N⁶-methyladenine, 5-methylcytosine and N⁴-methylcytosine. The aim of this study was to develop a method to detect all three types of DNA methylation in complete genomic DNA. It was previously shown that N⁶-methyladenine and 5-methylcytosine in plasmid and viral DNA can be detected by intersequence trace comparison of methylated and unmethylated DNA. We extended this method to include N⁴-methylcytosine detection in both in vitro and in vivo methylated DNA. Furthermore, application of intersequence trace comparison was extended to bacterial genomic DNA. Finally, we present evidence that intrasequence comparison suffices to detect methylated sites in genomic DNA. In conclusion, we present a method to detect all three natural types of DNA methylation in bacterial genomic DNA. This provides the possibility to define the complete methylome of any prokaryote.     

Aberrant CpG methylation of the TFAP2A gene constitutes a mechanism for loss of TFAP2A expression in human metastatic melanoma

Metastatic melanoma is a deadly treatment-resistant form of skin cancer whose global incidence is on the rise. During melanocyte transformation and melanoma progression the expression profile of many genes changes. Among these, a gene implicated in several steps of melanocyte development, TFAP2A, is frequently silenced; however, the molecular mechanism of TFAP2A silencing in human melanoma remains unknown. In this study, we measured TFAP2A mRNA expression in primary human melanocytes compared to 11 human melanoma samples by quantitative real-time RT-PCR. In addition, we assessed CpG DNA methylation of the TFAP2A promoter in these samples using bisulfite sequencing. Compared to primary melanocytes, which showed high TFAP2A mRNA expression and no promoter methylation, human melanoma samples showed decreased TFAP2A mRNA expression and increased promoter methylation. We further show that increased CpG methylation correlates with decreased TFAP2A mRNA expression. Using The Cancer Genome Atlas, we further identified TFAP2A as a gene displaying among the most decreased expression in stage 4 melanomas vs. non-stage 4 melanomas, and whose CpG methylation was frequently associated with lack of mRNA expression. Based on our data, we conclude that TFAP2A expression in human melanomas can be silenced by aberrant CpG methylation of the TFAP2A promoter. We have identified aberrant CpG DNA methylation as an epigenetic mark associated with TFAP2A silencing in human melanoma that could have significant implications for the therapy of human melanoma using epigenetic modifying drugs.     

A suite of DNA methylation markers that can detect most common human cancers

Cancer-specific DNA methylation from the tumor derived fraction of cell free DNA found in blood samples could be used for minimally invasive detection and monitoring of cancer. The knowledge of marker regions with cancer-specific DNA methylation is necessary to the success of such a process. We analyzed the largest cancer DNA methylation dataset available—TCGA Illumina HumanMethylation450 data with over 8,500 tumors—in order to find cancer-specific DNA methylation markers for most common human cancers. First, we identified differentially methylated regions for individual cancer types and those were further filtered against data from normal tissues to obtain marker regions with cancer-specific methylation, resulting in a total of 1,250 hypermethylated and 584 hypomethylated marker CpGs. From hypermethylated markers, optimal sets of six markers for each TCGA cancer type were chosen that could identify most tumors with high specificity and sensitivity [area under the curve (AUC): 0.969-1.000] and a universal 12 marker set that can detect tumors of all 33 TCGA cancer types (AUC >0.84). In addition to hundreds of new DNA methylation markers, our approach also identified markers that are in current clinical use, SEPT9 and GSTP1, indicating the validity of our approach and a significant potential utility for the newly discovered markers. The hypermethylated markers are linked to polycomb associated loci and a significant fraction of the discovered markers is within noncoding RNA genes; one of the best markers is MIR129-2. Future clinical testing of herein discovered markers will confirm new markers that will improve minimally invasive diagnosis and monitoring for multiple cancers.     

Deciphering the heterogeneity in DNA methylation patterns during stem cell differentiation and reprogramming

Background

Human induced pluripotent stem cells (iPSCs) have a wide range of applications throughout the fields of basic research, disease modeling and drug screening. Epigenetic instable iPSCs with aberrant DNA methylation may divide and differentiate into cancer cells. Unfortunately, little effort has been taken to compare the epigenetic variation in iPSCs with that in differentiated cells. Here, we developed an analytical procedure to decipher the DNA methylation heterogeneity of mixed cells and further exploited it to quantitatively assess the DNA methylation variation in the methylomes of adipose-derived stem cells (ADS), mature adipocytes differentiated from ADS cells (ADS-adipose) and iPSCs reprogrammed from ADS cells (ADS-iPSCs).

Results

We observed that the degree of DNA methylation variation varies across distinct genomic regions with promoter and 5’UTR regions exhibiting low methylation variation and Satellite showing high methylation variation. Compared with differentiated cells, ADS-iPSCs possess globally decreased methylation variation, in particular in repetitive elements. Interestingly, DNA methylation variation decreases in promoter regions during differentiation but increases during reprogramming. Methylation variation in promoter regions is negatively correlated with gene expression. In addition, genes showing a bipolar methylation pattern, with both completely methylated and completely unmethylated reads, are related to the carbohydrate metabolic process, cellular development, cellular growth, proliferation, etc.

Conclusions

This study delivers a way to detect cell-subset specific methylation genes in a mixed cell population and provides a better understanding of methylation dynamics during stem cell differentiation and reprogramming.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-978) contains supplementary material, which is available to authorized users.     

Leukocyte DNA as Surrogate for the Evaluation of Imprinted Loci Methylation in Mammary Tissue DNA

There is growing interest in identifying surrogate tissues to identify epimutations in cancer patients since primary target tissues are often difficult to obtain. Methylation patterns at imprinted loci are established during gametogenesis and post fertilization and their alterations have been associated with elevated risk of cancer. Methylation at several imprinted differentially methylated regions (GRB10 ICR, H19 ICR, KvDMR, SNRPN/SNURF ICR, IGF2 DMR0, and IGF2 DMR2) were analyzed in DNA from leukocytes and mammary tissue (normal, benign diseases, or malignant tumors) from 87 women with and without breast cancer (average age of cancer patients: 53; range: 31–77). Correlations between genomic variants and DNA methylation at the studied loci could not be assessed, making it impossible to exclude such effects. Methylation levels observed in leukocyte and mammary tissue DNA were close to the 50% expected for monoallellic methylation. While no correlation was observed between leukocyte and mammary tissue DNA methylation for most of the analyzed imprinted genes, Spearman''s correlations were statistically significant for IGF2 DMR0 and IGF2 DMR2, although absolute methylation levels differed. Leukocyte DNA methylation levels of selected imprinted genes may therefore serve as surrogate markers of DNA methylation in cancer tissue.