A systematic comparison of quantitative high-resolution DNA methylation analysis and methylation-specific PCR

Assessment of DNA methylation has become a critical factor for the identification, development and application of methylation based biomarkers. Here we describe a systematic comparison of a quantitative high-resolution mass spectrometry-based approach (MassARRAY), pyrosequencing and the broadly used methylation-specific PCR (MSP) technique analyzing clinically relevant epigenetically silenced genes in acute myeloid leukemia (AML). By MassARRAY and pyrosequencing, we identified significant DNA methylation differences at the ID4 gene promoter and in the 5′ region of members of the SFRP gene family in 62 AML patients compared with healthy controls. We found a good correlation between data obtained by MassARRAY and pyrosequencing (correlation coefficient R² = 0.88). MSP-based assessment of the identical samples showed less pronounced differences between AML patients and controls. By direct comparison of MSP-derived and MassARRAY-based methylation data as well as pyrosequencing, we could determine overestimation of DNA methylation data by MSP. We found sequence-context dependent highly variable cut-off values of quantitative DNA methylation values serving as discriminator for the two MSP methylation categories. Moreover, good agreements between quantitative methods and MSP could not be achieved for all investigated loci. Significant correlation of the quantitative assessment but not of MSP-derived methylation data with clinically important characteristics in our patient cohort demonstrated clinical relevance of quantitative DNA methylation assessment. Taken together, while MSP is still the most commonly applied technique for DNA methylation assessment, our data highlight advantages of quantitative approaches for precise characterization and reliable biomarker use of aberrant DNA methylation in primary patient samples, particularly.     

Constitutional promoter methylation and risk of familial melanoma

Constitutional epigenetic changes detected in blood or non-disease involving tissues have been associated with disease susceptibility. We measured promoter methylation of CDKN2A (p16 and p14ARF) and 13 melanoma-related genes using bisulfite pyrosequencing of blood DNA from 114 cases and 122 controls in 64 melanoma-prone families (26 segregating CDKN2A germline mutations). We also obtained gene expression data for these genes using microarrays from the same blood samples. We observed that CDKN2A epimutation is rare in melanoma families, and therefore is unlikely to cause major susceptibility in families without CDKN2A mutations. Although methylation levels for most gene promoters were very low (<5%), we observed a significantly reduced promoter methylation (odds ratio = 0.63, 95% confidence interval = 0.50, 0.80, P < 0.001) and increased expression (fold change = 1.27, P = 0.048) for TNFRSF10C in melanoma cases. Future research in large prospective studies using both normal and melanoma tissues is required to assess the significance of TNFRSF10C methylation and expression changes in melanoma susceptibility.     

Association between length of gestation and cervical DNA methylation of PTGER2 and LINE 1-HS

Worldwide, more than 1 in 10 infants is born prior to 37 weeks gestation. Preterm birth can lead to increased mortality risk and poor life-long health and neurodevelopmental outcomes. Whether environmental risk factors affect preterm birth through epigenetic phenomena is largely unstudied. We sought to determine whether preterm risk factors, such as smoke exposure and education, were associated with cervical DNA methylation in the prostaglandin E receptor 2 gene (PTGER2) and a repetitive element, long interspersed nuclear element-1 Homo sapiens-specific (LINE 1-HS). Second, we aimed to determine whether mid-pregnancy DNA methylation of these regions in cervical samples could predict the length of gestation. We obtained a cervical swab between 16–19 weeks gestation from 80 women participating in a Mexico City birth cohort, used pyrosequencing to analyze DNA methylation of PTGER2 and LINE 1-HS, and examined associations with maternal covariates. We used accelerated failure time models to analyze associations of DNA methylation with the length of gestation. DNA methylation of both sequences was associated with Pap smear inflammation. LINE 1-HS methylation was associated with smoke exposure, BMI and parity. In adjusted models, gestations were 3.3 days longer (95%CI 0.6, 6.0) for each interquartile range of PTGER2 DNA methylation. Higher LINE 1-HS methylation was associated with shorter gestations (-3.3 days, 95%CI -6.5, -0.2). In conclusion, cervical DNA methylation was associated with risk factors for preterm birth and the length of gestation.     

Identification of differentially methylated markers among cytogenetic risk groups of acute myeloid leukemia

Aberrant DNA methylation is known to occur in cancer, including hematological malignancies such as acute myeloid leukemia (AML). However, less is known about whether specific methylation profiles characterize specific subcategories of AML. We examined this issue by using comprehensive high-throughput array-based relative methylation analysis (CHARM) to compare methylation profiles among patients in different AML cytogenetic risk groups. We found distinct profiles in each group, with the high-risk group showing overall increased methylation compared with low- and mid-risk groups. The differentially methylated regions (DMRs) distinguishing cytogenetic risk groups of AML were enriched in the CpG island shores. Specific risk-group associated DMRs were located near genes previously known to play a role in AML or other malignancies, such as MN1, UHRF1, HOXB3, and HOXB4, as well as TRIM71, the function of which in cancer is not well characterized. These findings were verified by quantitative bisulfite pyrosequencing and by comparison with results available at the TCGA cancer genome browser. To explore the potential biological significance of the observed methylation changes, we correlated our findings with gene expression data available through the TCGA database. The results showed that decreased methylation at HOXB3 and HOXB4 was associated with increased gene expression of both HOXB genes specific to the mid-risk AML, while increased DNA methylation at DCC distinctive to the high-risk AML was associated with increased gene expression. Our results suggest that the differential impact of cytogenetic changes on AML prognosis may, in part, be mediated by changes in methylation.     

Correlation between BRAF mutation and promoter methylation of TIMP3, RARβ2 and RASSF1A in thyroid cancer

Our aim was to comprehensively analyze promoter hypermethylation of a panel of novel and known methylation markers for thyroid neoplasms and to establish their relationship with BRAF mutation and clinicopathologic parameters of thyroid cancer. A cohort of thyroid tumors, consisting of 44 cancers and 44 benign thyroid lesions, as well as 15 samples of adjacent normal thyroid tissue, was evaluated for BRAF mutation and promoter hypermethylation. Genes for quantitative methylation specific PCR (QMSP) were selected by a candidate gene approach. Twenty-two genes were tested: TSHR, RASSF1A, RARβ2, DAPK, hMLH1, ATM, S100, p16, CTNNB1, GSTP1, CALCA, TIMP3, TGFßR2, THBS1, MINT1, CTNNB1, MT1G, PAK3, NISCH, DCC, AIM1 and KIF1A. The PCR-based “mutector assay” was used to detect BRAF mutation. All p values reported are two sided. Considerable overlap was seen in the methylation markers among the different tissue groups. Significantly higher methylation frequency and level were observed for KIF1A and RARß2 in cancer samples compared with benign tumors. A negative correlation between BRAF mutation and RASSF1A methylation, and a positive correlation with RARß2 methylation were observed in accordance with previous results. In addition, positive correlation with TIMP3 and a marginal correlation with DCC methylation were observed. The present study constitutes a comprehensive promoter methylation profile of thyroid neoplasia and shows that results must be analyzed in a tissue-specific manner to identify clinically useful methylation markers. Integration of genetic and epigenetic changes in thyroid cancer will help identify relevant biologic pathways that drive its development.     

Non-referenced genome assembly from epigenomic short-read data

Current computational methods used to analyze changes in DNA methylation and chromatin modification rely on sequenced genomes. Here we describe a pipeline for the detection of these changes from short-read sequence data that does not require a reference genome. Open source software packages were used for sequence assembly, alignment, and measurement of differential enrichment. The method was evaluated by comparing results with reference-based results showing a strong correlation between chromatin modification and gene expression. We then used our de novo sequence assembly to build the DNA methylation profile for the non-referenced Psammomys obesus genome. The pipeline described uses open source software for fast annotation and visualization of unreferenced genomic regions from short-read data.     

Comparative analysis of methods to determine DNA methylation levels of a tumor-related microRNA gene

Quantifying levels of DNA methylation in tumors is a useful approach for the identification of potential tumor suppressors and to find biomarkers that can be used as prognostic or therapeutic indicators. In the current study, we compared three methods commonly used for quantifying DNA methylation—bisulfite pyrosequencing, quantitative methylation-specific PCR (Q-MSP), and MethyLight—by focusing on the CpG island of the gene encoding the microRNA-34b and microRNA-34c (miR-34b/c); aberrant regulation of this miR is associated with various human malignancies, including gastric cancer. Standard curve analysis using control DNA samples demonstrated the highest quantitative accuracy in Q-MSP analysis. We also carried out methylation analysis using gastric mucosa specimens obtained from gastric cancer patients. We found a high correlation between methylation levels determined by Q-MSP and those determined by MethyLight (R² = 0.952), whereas the results of bisulfite pyrosequencing and the other two methods were less well correlated (R² = 0.864 and R² = 0.804 for Q-MSP and MethyLight, respectively). This may reflect possible PCR bias in the pyrosequencing technique, which we show can be corrected for by applying a cubic approximate equation to the original data. Thus, although results obtained by the different DNA methylation analysis techniques are largely comparable, an appropriate correction may be necessary for stringent comparison.     

Synergistic chromatin repression of the tumor suppressor gene RARB in human prostate cancers

DNA methylation and polycomb proteins are well-known mediators of epigenetic silencing in mammalian cells. Usually described as mutually exclusive, this statement is today controversial and recent in vitro studies suggest the co-existence of both repressor systems. We addressed this issue in the study of Retinoic Acid Receptor β (RARβ), a tumor suppressor gene frequently silenced in prostate cancer. We found that the RARβ promoter is hypermethylated in all studied prostate tumors and methylation levels are positively correlated with H3K27me3 enrichments. Thus, by using bisulfite conversion and pyrosequencing of immunoprecipitated H3K27me3 chromatin, we demonstrated that DNA methylation and polycomb repression co-exist in vivo at this locus. We found this repressive association in 6/6 patient tumor samples of different Gleason score, suggesting a strong interplay of DNA methylation and EZH2 to silence RARβ during prostate tumorigenesis.     

Gene promoter methylation in colorectal cancer and healthy adjacent mucosa specimens: Correlation with physiological and pathological characteristics,and with biomarkers of one-carbon metabolism

We evaluated the promoter methylation levels of the APC, MGMT, hMLH1, RASSF1A and CDKN2A genes in 107 colorectal cancer (CRC) samples and 80 healthy adjacent tissues. We searched for correlation with both physical and pathological features, polymorphisms of folate metabolism pathway genes (MTHFR, MTRR, MTR, RFC1, TYMS, and DNMT3B), and data on circulating folate, vitamin B12 and homocysteine, which were available in a subgroup of the CRC patients. An increased number of methylated samples were found in CRC respect to adjacent healthy tissues, with the exception of APC, which was also frequently methylated in healthy colonic mucosa. Statistically significant associations were found between RASSF1A promoter methylation and tumor stage, and between hMLH1 promoter methylation and tumor location. Increasing age positively correlated with both hMLH1 and MGMT methylation levels in CRC tissues, and with APC methylation levels in the adjacent healthy mucosa. Concerning gender, females showed higher hMLH1 promoter methylation levels with respect to males. In CRC samples, the MTR 2756AG genotype correlated with higher methylation levels of RASSF1A, and the TYMS 1494 6bp ins/del polymorphism correlated with the methylation levels of both APC and hMLH1. In adjacent healthy tissues, MTR 2756AG and TYMS 1494 6bp del/del genotypes correlated with APC and MGMT promoter methylation, respectively. Low folate levels were associated with hMLH1 hypermethylation. Present results support the hypothesis that DNA methylation in CRC depends from both physiological and environmental factors, with one-carbon metabolism largely involved in this process.     

Genome-wide and single-base resolution DNA methylomes of the Pacific oyster Crassostrea gigas provide insight into the evolution of invertebrate CpG methylation

Background

Studies of DNA methylomes in a wide range of eukaryotes have revealed both conserved and divergent characteristics of DNA methylation among phylogenetic groups. However, data on invertebrates particularly molluscs are limited, which hinders our understanding of the evolution of DNA methylation in metazoa. The sequencing of the Pacific oyster Crassostrea gigas genome provides an opportunity for genome-wide profiling of DNA methylation in this model mollusc.

Results

Homologous searches against the C. gigas genome identified functional orthologs for key genes involved in DNA methylation: DNMT1, DNMT2, DNMT3, MBD2/3 and UHRF1. Whole-genome bisulfite sequencing (BS-seq) of the oyster’s mantle tissues revealed that more than 99% methylation modification was restricted to cytosines in CpG context and methylated CpGs accumulated in the bodies of genes that were moderately expressed. Young repeat elements were another major targets of CpG methylation in oysters. Comparison with other invertebrate methylomes suggested that the 5’-end bias of gene body methylation and the negative correlation between gene body methylation and gene length were the derived features probably limited to the insect lineage. Interestingly, phylostratigraphic analysis showed that CpG methylation preferentially targeted genes originating in the common ancestor of eukaryotes rather than the oldest genes originating in the common ancestor of cellular organisms.

Conclusions

Comparative analysis of the oyster DNA methylomes and that of other animal species revealed that the characteristics of DNA methylation were generally conserved during invertebrate evolution, while some unique features were derived in the insect lineage. The preference of methylation modification on genes originating in the eukaryotic ancestor rather than the oldest genes is unexpected, probably implying that the emergence of methylation regulation in these ''relatively young’ genes was critical for the origin and radiation of eukaryotes.

Electronic supplementary material

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

Major Sperm Protein Genes from Globodera rostochiensis

Three genes in the major sperm protein (MSP) gene family from the potato cyst nematode Globodera rostochiensis were cloned and sequenced. In contrast to the absence of introns in Caenorhabditis elegans MSP genes, these genes in G. rostochiensis contained a 57 nucleotide intron, with normal exon-intron boundaries, in the same relative location as the intron in Onchocerca volvulus. The MSP genes of G. rostochiensis had putative CAAT, TATA, and polyadenylation signals. The predicted G. rostochiensis MSP gene product is 126 amino acids long, one residue shorter than the products in the other species. The comparison of MSP amino acid sequences from four diverse nematode species suggests that O. volvulus, Ascaris suum, and C. elegans may be more closely related to each other than they are to G. rostochiensis.     

Prognostic value of three different methods of MGMT promoter methylation analysis in a prospective trial on newly diagnosed glioblastoma

Hypermethylation in the promoter region of the MGMT gene encoding the DNA repair protein O⁶-methylguanine-DNA methyltransferase is among the most important prognostic factors for patients with glioblastoma and predicts response to treatment with alkylating agents like temozolomide. Hence, the MGMT status is widely determined in most clinical trials and frequently requested in routine diagnostics of glioblastoma. Since various different techniques are available for MGMT promoter methylation analysis, a generally accepted consensus as to the most suitable diagnostic method remains an unmet need. Here, we assessed methylation-specific polymerase chain reaction (MSP) as a qualitative and semi-quantitative method, pyrosequencing (PSQ) as a quantitative method, and methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA) as a semi-quantitative method in a series of 35 formalin-fixed, paraffin-embedded glioblastoma tissues derived from patients treated in a prospective clinical phase II trial that tested up-front chemoradiotherapy with dose-intensified temozolomide (UKT-05). Our goal was to determine which of these three diagnostic methods provides the most accurate prediction of progression-free survival (PFS). The MGMT promoter methylation status was assessable by each method in almost all cases (n = 33/35 for MSP; n = 35/35 for PSQ; n = 34/35 for MS-MLPA). We were able to calculate significant cut-points for the continuous methylation signals at each CpG site analysed by PSQ (range, 11.5 to 44.9%) and at one CpG site assessed by MS-MLPA (3.6%) indicating that a dichotomisation of continuous methylation data as a prerequisite for comparative survival analyses is feasible. Our results show that, unlike MS-MLPA, MSP and PSQ provide a significant improvement of predicting PFS compared with established clinical prognostic factors alone (likelihood ratio tests: p<0.001). Conclusively, taking into consideration prognostic value, cost effectiveness and ease of use, we recommend pyrosequencing for analyses of MGMT promoter methylation in high-throughput settings and MSP for clinical routine diagnostics with low sample numbers.     

Assessing combined methylation-sensitive high resolution melting and pyrosequencing for the analysis of heterogeneous DNA methylation

Heterogeneous DNA methylation leads to difficulties in accurate detection and quantification of methylation. Methylation-sensitive high resolution melting (MS-HRM) is unique among regularly used methods for DNA methylation analysis in that heterogeneous methylation can be readily identified, although not quantified, by inspection of the melting curves. Bisulfite pyrosequencing has been used to estimate the level of heterogeneous methylation by quantifying methylation levels present at individual CpG dinucleotides. Sequentially combining the two methodologies using MS-HRM to screen the amplification products prior to bisulfite pyrosequencing would be advantageous. This would not only replace the quality control step using agarose gel analysis prior to the pyrosequencing step but would also provide important qualitative information in its own right. We chose to analyze DAPK1 as it is an important tumor suppressor gene frequently heterogeneously methylated in a number of malignancies, including chronic lymphocytic leukemia (CLL). A region of the DAPK1 promoter was analyzed in ten CLL samples by MS-HRM. By using a biotinylated primer, bisulfite pyrosequencing could be used to directly analyze the samples. MS-HRM revealed the presence of various extents of heterogeneous DAPK1 methylation in all CLL samples. Further analysis of the biotinylated MS-HRM products by bisulfite pyrosequencing provided quantitative information for each CpG dinucleotide analyzed, and confirmed the presence of heterogeneous DNA methylation. Whereas each method could be used individually, MS-HRM and bisulfite pyrosequencing provided complementary information for the assessment of heterogeneous methylation.Key words: MS-HRM, pyrosequencing, digital PCR, heterogeneous DNA methylation, DAPK1, chronic lymphocytic leukemia     

Genome-wide DNA methylation changes in skeletal muscle between young and middle-aged pigs

Background

Age-related physiological, biochemical and functional changes in mammalian skeletal muscle have been shown to begin at the mid-point of the lifespan. However, the underlying changes in DNA methylation that occur during this turning point of the muscle aging process have not been clarified. To explore age-related genomic methylation changes in skeletal muscle, we employed young (0.5 years old) and middle-aged (7 years old) pigs as models to survey genome-wide DNA methylation in the longissimus dorsi muscle using a methylated DNA immunoprecipitation sequencing approach.

Results

We observed a tendency toward a global loss of DNA methylation in the gene-body region of the skeletal muscle of the middle-aged pigs compared with the young group. We determined the genome-wide gene expression pattern in the longissimus dorsi muscle using microarray analysis and performed a correlation analysis using DMR (differentially methylated region)-mRNA pairs, and we found a significant negative correlation between the changes in methylation levels within gene bodies and gene expression. Furthermore, we identified numerous genes that show age-related methylation changes that are potentially involved in the aging process. The methylation status of these genes was confirmed using bisulfite sequencing PCR. The genes that exhibited a hypomethylated gene body in middle-aged pigs were over-represented in various proteolysis and protein catabolic processes, suggesting an important role for these genes in age-related muscle atrophy. In addition, genes associated with tumorigenesis exhibited aged-related differences in methylation and expression levels, suggesting an increased risk of disease associated with increased age.

Conclusions

This study provides a comprehensive analysis of genome-wide DNA methylation patterns in aging pig skeletal muscle. Our findings will serve as a valuable resource in aging studies, promoting the pig as a model organism for human aging research and accelerating the development of comparative animal models in aging research.

Electronic supplementary material

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

Prenatal stress decreases Bdnf expression and increases methylation of Bdnf exon IV in rats

There is ample evidence that exposure to stress during gestation increases the risk of the offspring to develop mood disorders. Brain-derived neurotrophic factor (Bdnf) plays a critical role during neuronal development and is therefore a prime candidate to modulate neuronal signaling in adult offspring of rat dams that were stressed during gestation. In the current study, we tested the hypothesis that alterations in Bdnf expression in prenatally stressed (PNS) offspring are mediated by changes in DNA methylation in exons IV and VI of the Bdnf gene. We observed decreased Bdnf expression in the amygdala and hippocampus of prenatally stressed rats both at weaning and in adulthood. This decrease in Bdnf expression was accompanied by increased DNA methylation in Bdnf exon IV in the amygdala and hippocampus, suggesting that PNS-induced reduction in Bdnf expression may, at least in part, be mediated by increased DNA methylation of Bdnf exon IV. Expression of DNA methyltransferases (Dnmt) 1 and 3a was increased in PNS rats in the amygdala and hippocampus. Our data suggest that PNS induces decreases in Bdnf expression that may at least in part be mediated by increased DNA methylation of Bdnf exon IV.