The prognostic significance of whole blood global and specific DNA methylation levels in gastric adenocarcinoma

Background

Epigenetics, particularly DNA methylation, has recently been elucidated as important in gastric cancer (GC) initiation and progression. We investigated the clinical and prognostic importance of whole blood global and site-specific DNA methylation in GC.

Methods

Genomic DNA was extracted from the peripheral blood of 105 Omani GC patients at diagnosis. DNA methylation was quantified by pyrosequencing of global DNA and specific gene promoter regions at 5 CpG sites for CDH1, 7 CpG sites for p16, 4 CpG sites for p53, and 3 CpG sites for RUNX3. DNA methylation levels in patients were categorized into low, medium, and high tertiles. Associations between methylation level category and clinicopathological features were evaluated using χ² tests. Survival analyses were carried out using the Kaplan-Meier method and log rank test. A backward conditional Cox proportional hazards regression model was used to identify independent predictors of survival.

Results

Older GC patients had increased methylation levels at specific CpG sites within the CDH1, p53, and RUNX-3 promoters. Male gender was significantly associated with reduced global and increased site-specific DNA methylation levels in CDH1, p16, and p53 promoters. Global DNA low methylation level was associated with better survival on univariate analysis. Patients with high and medium methylation vs. low methylation levels across p16 promoter CpG sites, site 2 in particular, had better survival. Multivariate analysis showed that global DNA hypermethylation was a significant independent predictor of worse survival (hazard ratio (HR) = 2.0, 95% CI: 1.1–3.8; p = 0.02) and high methylation mean values across p16 promoter sites 1–7 were associated with better survival with HR of 0.3 (95% CI, 0.1–0.8; p = 0.02) respectively.

Conclusions

Analysis of global and site-specific DNA methylation in peripheral blood by pyrosequencing provides quantitative DNA methylation values that may serve as important prognostic indicators.     

LINE-1 methylation levels in leukocyte DNA and risk of renal cell cancer

Purpose

Leukocyte global DNA methylation levels are currently being considered as biomarkers of cancer susceptibility and have been associated with risk of several cancers. In this study, we aimed to examine the association between long interspersed nuclear elements (LINE-1) methylation levels, as a biomarker of global DNA methylation in blood cell DNA, and renal cell cancer risk.

Experimental Design

LINE-1 methylation of bisulfite-converted genomic DNA isolated from leukocytes was quantified by pyrosequencing measured in triplicate, and averaged across 4 CpG sites. A total of 328 RCC cases and 654 controls frequency-matched(2∶1) on age(±5years), sex and study center, from a large case-control study conducted in Central and Eastern Europe were evaluated.

Results

LINE-1 methylation levels were significantly higher in RCC cases with a median of 81.97% (interquartile range[IQR]: 80.84–83.47) compared to 81.67% (IQR: 80.35–83.03) among controls (p = 0.003, Wilcoxon). Compared to the lowest LINE-1 methylation quartile(Q1), the adjusted ORs for increasing methylation quartiles were as follows: OR(Q2) = 1.84(1.20−2.81), OR(Q3) = 1.72(1.11−2.65) and OR(Q4) = 2.06(1.34−3.17), with a p-trend = 0.004. The association was stronger among current smokers (p-trend<0.001) than former or never smokers (p-interaction = 0.03). To eliminate the possibility of selection bias among controls, the relationship between LINE-1 methylation and smoking was evaluated and confirmed in a case-only analysis, as well.

Conclusions

Higher levels of LINE-1 methylation appear to be positively associated with RCC risk, particularly among current smokers. Further investigations using both post- and pre-diagnostic genomic DNA is warranted to confirm findings and will be necessary to determine whether the observed differences occur prior to, or as a result of carcinogenesis.     

Lung Fibroblasts from Patients with Idiopathic Pulmonary Fibrosis Exhibit Genome-Wide Differences in DNA Methylation Compared to Fibroblasts from Nonfibrotic Lung

Excessive fibroproliferation is a central hallmark of idiopathic pulmonary fibrosis (IPF), a chronic, progressive disorder that results in impaired gas exchange and respiratory failure. Fibroblasts are the key effector cells in IPF, and aberrant expression of multiple genes contributes to their excessive fibroproliferative phenotype. DNA methylation changes are critical to the development of many diseases, but the DNA methylome of IPF fibroblasts has never been characterized. Here, we utilized the HumanMethylation 27 array, which assays the DNA methylation level of 27,568 CpG sites across the genome, to compare the DNA methylation patterns of IPF fibroblasts (n = 6) with those of nonfibrotic patient controls (n = 3) and commercially available normal lung fibroblast cell lines (n = 3). We found that multiple CpG sites across the genome are differentially methylated (as defined by P value less than 0.05 and fold change greater than 2) in IPF fibroblasts compared to fibroblasts from nonfibrotic controls. These methylation differences occurred both in genes recognized to be important in fibroproliferation and extracellular matrix generation, as well as in genes not previously recognized to participate in those processes (including organ morphogenesis and potassium ion channels). We used bisulfite sequencing to independently verify DNA methylation differences in 3 genes (CDKN2B, CARD10, and MGMT); these methylation changes corresponded with differences in gene expression at the mRNA and protein level. These differences in DNA methylation were stable throughout multiple cell passages. DNA methylation differences may thus help to explain a proportion of the differences in gene expression previously observed in studies of IPF fibroblasts. Moreover, significant variability in DNA methylation was observed among individual IPF cell lines, suggesting that differences in DNA methylation may contribute to fibroblast heterogeneity among patients with IPF. These results demonstrate that IPF fibroblasts exhibit global differences in DNA methylation that may contribute to the excessive fibroproliferation associated with this disease.     

A Genome-Wide Methylation Study on Essential Hypertension in Young African American Males

Objective

There is emerging evidence from animal studies suggesting a key role for methylation in the pathogenesis of essential hypertension. However, to date, very few studies have investigated the role of methylation in the development of human hypertension, and none has taken a genome-wide approach. Based on the recent studies that highlight the involvement of inflammation in the development of hypertension, we hypothesize that changes in DNA methylation of leukocytes are involved in the pathogenesis of hypertension.

Method & Results

We conducted a genome-wide methylation analysis on 8 hypertensive cases and 8 normotensive age-matched controls aged 14–23 years and performed validation of the most significant CpG sites in 2 genes in an independent sample of 36 hypertensive cases and 60 normotensive controls aged 14–30 years. Validation of the CpG sites in the SULF1 gene was further conducted in a second replication sample of 36 hypertensive cases and 34 controls aged 15.8–40 years. A CpG site in the SULF1 gene showed higher methylation levels in cases than in healthy controls in the genome-wide step (p = 6.2×10⁻⁵), which was confirmed in the validation step (p = 0.011) for subjects ≤30 years old but was not significant for subjects of all ages combined (p = 0.095).

Conclusion

The identification of a difference in a blood leukocyte DNA methylation site between hypertensive cases and normotensive controls suggests that changes in DNA methylation may play an important role in the pathogenesis of hypertension. The age dependency of the effect further suggests complexity of epigenetic regulation in this age-related disease.     

The Relationship of DNA Methylation with Age,Gender and Genotype in Twins and Healthy Controls

Cytosine-5 methylation within CpG dinucleotides is a potentially important mechanism of epigenetic influence on human traits and disease. In addition to influences of age and gender, genetic control of DNA methylation levels has recently been described. We used whole blood genomic DNA in a twin set (23 MZ twin-pairs and 23 DZ twin-pairs, N = 92) as well as healthy controls (N = 96) to investigate heritability and relationship with age and gender of selected DNA methylation profiles using readily commercially available GoldenGate bead array technology. Despite the inability to detect meaningful methylation differences in the majority of CpG loci due to tissue type and locus selection issues, we found replicable significant associations of DNA methylation with age and gender. We identified associations of genetically heritable single nucleotide polymorphisms with large differences in DNA methylation levels near the polymorphism (cis effects) as well as associations with much smaller differences in DNA methylation levels elsewhere in the human genome (trans effects). Our results demonstrate the feasibility of array-based approaches in studies of DNA methylation and highlight the vast differences between individual loci. The identification of CpG loci of which DNA methylation levels are under genetic control or are related to age or gender will facilitate further studies into the role of DNA methylation and disease.     

Identification of type 1 diabetes-associated DNA methylation variable positions that precede disease diagnosis

Monozygotic (MZ) twin pair discordance for childhood-onset Type 1 Diabetes (T1D) is ∼50%, implicating roles for genetic and non-genetic factors in the aetiology of this complex autoimmune disease. Although significant progress has been made in elucidating the genetics of T1D in recent years, the non-genetic component has remained poorly defined. We hypothesized that epigenetic variation could underlie some of the non-genetic component of T1D aetiology and, thus, performed an epigenome-wide association study (EWAS) for this disease. We generated genome-wide DNA methylation profiles of purified CD14⁺ monocytes (an immune effector cell type relevant to T1D pathogenesis) from 15 T1D–discordant MZ twin pairs. This identified 132 different CpG sites at which the direction of the intra-MZ pair DNA methylation difference significantly correlated with the diabetic state, i.e. T1D–associated methylation variable positions (T1D–MVPs). We confirmed these T1D–MVPs display statistically significant intra-MZ pair DNA methylation differences in the expected direction in an independent set of T1D–discordant MZ pairs (P = 0.035). Then, to establish the temporal origins of the T1D–MVPs, we generated two further genome-wide datasets and established that, when compared with controls, T1D–MVPs are enriched in singletons both before (P = 0.001) and at (P = 0.015) disease diagnosis, and also in singletons positive for diabetes-associated autoantibodies but disease-free even after 12 years follow-up (P = 0.0023). Combined, these results suggest that T1D–MVPs arise very early in the etiological process that leads to overt T1D. Our EWAS of T1D represents an important contribution toward understanding the etiological role of epigenetic variation in type 1 diabetes, and it is also the first systematic analysis of the temporal origins of disease-associated epigenetic variation for any human complex disease.     

Blood diagnostic biomarkers for major depressive disorder using multiplex DNA methylation profiles: discovery and validation

Aberrant DNA methylation in the blood of patients with major depressive disorder (MDD) has been reported in several previous studies. However, no comprehensive studies using medication-free subjects with MDD have been conducted. Furthermore, the majority of these previous studies has been limited to the analysis of the CpG sites in CpG islands (CGIs) in the gene promoter regions. The main aim of the present study is to identify DNA methylation markers that distinguish patients with MDD from non-psychiatric controls. Genome-wide DNA methylation profiling of peripheral leukocytes was conducted in two set of samples, a discovery set (20 medication-free patients with MDD and 19 controls) and a replication set (12 medication-free patients with MDD and 12 controls), using Infinium HumanMethylation450 BeadChips. Significant diagnostic differences in DNA methylation were observed at 363 CpG sites in the discovery set. All of these loci demonstrated lower DNA methylation in patients with MDD than in the controls, and most of them (85.7%) were located in the CGIs in the gene promoter regions. We were able to distinguish patients with MDD from the control subjects with high accuracy in the discriminant analysis using the top DNA methylation markers. We also validated these selected DNA methylation markers in the replication set. Our results indicate that multiplex DNA methylation markers may be useful for distinguishing patients with MDD from non-psychiatric controls.     

Fetal DNA Methylation Associates with Early Spontaneous Preterm Birth and Gestational Age

Spontaneous preterm birth (PTB, <37 weeks gestation) is a major public health concern, and children born preterm have a higher risk of morbidity and mortality throughout their lives. Recent studies suggest that fetal DNA methylation of several genes varies across a range of gestational ages (GA), but it is not yet clear if fetal epigenetic changes associate with PTB. The objective of this study is to interrogate methylation patterns across the genome in fetal leukocyte DNA from African Americans with early PTB (24^1/7–34^0/7 weeks; N = 22) or term births (39^0/7–40^6/7weeks; N = 28) and to evaluate the association of each CpG site with PTB and GA. DNA methylation was assessed across the genome with the HumanMethylation450 BeadChip. For each individual sample and CpG site, the proportion of DNA methylation was estimated. The associations between methylation and PTB or GA were evaluated by fitting a separate linear model for each CpG site, adjusting for relevant covariates. Overall, 29 CpG sites associated with PTB (FDR<.05; 5.7×10⁻¹⁰<p<2.9×10⁻⁶) independent of GA. Also, 9637 sites associated with GA (FDR<.05; 9.5×10⁻¹⁶<p<1.0×10⁻³), with 61.8% decreasing in methylation with shorter GA. GA-associated CpG sites were depleted in the CpG islands of their respective genes (p<2.2×10⁻¹⁶). Gene set enrichment analysis (GSEA) supported enrichment of GA-associated CpG sites in genes that play a role in embryonic development as well as the extracellular matrix. Additionally, this study replicated the association of several CpG sites associated with gestational age in other studies (CRHBP, PIK3CD and AVP). Dramatic differences in fetal DNA methylation are evident in fetuses born preterm versus at term, and the patterns established at birth may provide insight into the long-term consequences associated with PTB.     

Comprehensive analysis of MGMT promoter methylation: correlation with MGMT expression and clinical response in GBM

O⁶-methylguanine DNA-methyltransferase (MGMT) promoter methylation has been identified as a potential prognostic marker for glioblastoma patients. The relationship between the exact site of promoter methylation and its effect on gene silencing, and the patient''s subsequent response to therapy, is still being defined. The aim of this study was to comprehensively characterize cytosine-guanine (CpG) dinucleotide methylation across the entire MGMT promoter and to correlate individual CpG site methylation patterns to mRNA expression, protein expression, and progression-free survival. To best identify the specific MGMT promoter region most predictive of gene silencing and response to therapy, we determined the methylation status of all 97 CpG sites in the MGMT promoter in tumor samples from 70 GBM patients using quantitative bisulfite sequencing. We next identified the CpG site specific and regional methylation patterns most predictive of gene silencing and improved progression-free survival. Using this data, we propose a new classification scheme utilizing methylation data from across the entire promoter and show that an analysis based on this approach, which we call 3R classification, is predictive of progression-free survival (HR  = 5.23, 95% CI [2.089–13.097], p<0.0001). To adapt this approach to the clinical setting, we used a methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA) test based on the 3R classification and show that this test is both feasible in the clinical setting and predictive of progression free survival (HR  = 3.076, 95% CI [1.301–7.27], p = 0.007). We discuss the potential advantages of a test based on this promoter-wide analysis and compare it to the commonly used methylation-specific PCR test. Further prospective validation of these two methods in a large independent patient cohort will be needed to confirm the added value of promoter wide analysis of MGMT methylation in the clinical setting.     

Identification of New Differentially Methylated Genes That Have Potential Functional Consequences in Prostate Cancer

Many differentially methylated genes have been identified in prostate cancer (PCa), primarily using candidate gene-based assays. Recently, several global DNA methylation profiles have been reported in PCa, however, each of these has weaknesses in terms of ability to observe global DNA methylation alterations in PCa. We hypothesize that there remains unidentified aberrant DNA methylation in PCa, which may be identified using higher resolution assay methods. We used the newly developed Illumina HumanMethylation450 BeadChip in PCa (n = 19) and adjacent normal tissues (n = 4) and combined these with gene expression data for identifying new DNA methylation that may have functional consequences in PCa development and progression. We also confirmed our methylation results in an independent data set. Two aberrant DNA methylation genes were validated among an additional 56 PCa samples and 55 adjacent normal tissues. A total 28,735 CpG sites showed significant differences in DNA methylation (FDR adjusted P<0.05), defined as a mean methylation difference of at least 20% between PCa and normal samples. Furthermore, a total of 122 genes had more than one differentially methylated CpG site in their promoter region and a gene expression pattern that was inverse to the direction of change in DNA methylation (e.g. decreased expression with increased methylation, and vice-versa). Aberrant DNA methylation of two genes, AOX1 and SPON2, were confirmed via bisulfate sequencing, with most of the respective CpG sites showing significant differences between tumor samples and normal tissues. The AOX1 promoter region showed hypermethylation in 92.6% of 54 tested PCa samples in contrast to only three out of 53 tested normal tissues. This study used a new BeadChip combined with gene expression data in PCa to identify novel differentially methylated CpG sites located within genes. The newly identified differentially methylated genes may be used as biomarkers for PCa diagnosis.     

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.     

Integrated genetic and epigenetic analysis identifies haplotype-specific methylation in the FTO type 2 diabetes and obesity susceptibility locus

Recent multi-dimensional approaches to the study of complex disease have revealed powerful insights into how genetic and epigenetic factors may underlie their aetiopathogenesis. We examined genotype-epigenotype interactions in the context of Type 2 Diabetes (T2D), focussing on known regions of genomic susceptibility. We assayed DNA methylation in 60 females, stratified according to disease susceptibility haplotype using previously identified association loci. CpG methylation was assessed using methylated DNA immunoprecipitation on a targeted array (MeDIP-chip) and absolute methylation values were estimated using a Bayesian algorithm (BATMAN). Absolute methylation levels were quantified across LD blocks, and we identified increased DNA methylation on the FTO obesity susceptibility haplotype, tagged by the rs8050136 risk allele A (p = 9.40×10⁻⁴, permutation p = 1.0×10⁻³). Further analysis across the 46 kb LD block using sliding windows localised the most significant difference to be within a 7.7 kb region (p = 1.13×10⁻⁷). Sequence level analysis, followed by pyrosequencing validation, revealed that the methylation difference was driven by the co-ordinated phase of CpG-creating SNPs across the risk haplotype. This 7.7 kb region of haplotype-specific methylation (HSM), encapsulates a Highly Conserved Non-Coding Element (HCNE) that has previously been validated as a long-range enhancer, supported by the histone H3K4me1 enhancer signature. This study demonstrates that integration of Genome-Wide Association (GWA) SNP and epigenomic DNA methylation data can identify potential novel genotype-epigenotype interactions within disease-associated loci, thus providing a novel route to aid unravelling common complex diseases.     

Genetic Control of Individual Differences in Gene-Specific Methylation in Human Brain

We have observed extensive interindividual differences in DNA methylation of 8590 CpG sites of 6229 genes in 153 human adult cerebellum samples, enriched in CpG island “shores” and at further distances from CpG islands. To search for genetic factors that regulate this variation, we performed a genome-wide association study (GWAS) mapping of methylation quantitative trait loci (mQTLs) for the 8590 testable CpG sites. cis association refers to correlation of methylation with SNPs within 1 Mb of a CpG site. 736 CpG sites showed phenotype-wide significant cis association with 2878 SNPs (after permutation correction for all tested markers and methylation phenotypes). In trans analysis of methylation, which tests for distant regulation effects, associations of 12 CpG sites and 38 SNPs remained significant after phenotype-wide correction. To examine the functional effects of mQTLs, we analyzed 85 genes that were with genetically regulated methylation we observed and for which we had quality gene expression data. Ten genes showed SNP-methylation-expression three-way associations—the same SNP simultaneously showed significant association with both DNA methylation and gene expression, while DNA methylation was significantly correlated with gene expression. Thus, we demonstrated that DNA methylation is frequently a heritable continuous quantitatively variable trait in human brain. Unlike allele-specific methylation, genetic polymorphisms mark both cis- and trans-regulatory genetic sites at measurable distances from their CpG sites. Some of the genetically regulated DNA methylation is directly connected with genetically regulated gene expression variation.