Differential DNA methylation with age displays both common and dynamic features across human tissues that are influenced by CpG landscape

Background

DNA methylation is an epigenetic modification that changes with age in human tissues, although the mechanisms and specificity of this process are still poorly understood. We compared CpG methylation changes with age across 283 human blood, brain, kidney, and skeletal muscle samples using methylation arrays to identify tissue-specific age effects.

Results

We found age-associated CpGs (ageCGs) that are both tissue-specific and common across tissues. Tissue-specific ageCGs are frequently located outside CpG islands with decreased methylation, and common ageCGs show the opposite trend. AgeCGs are significantly associated with poorly expressed genes, but those with decreasing methylation are linked with higher tissue-specific expression levels compared with increasing methylation. Therefore, tissue-specific gene expression may protect against common age-dependent methylation. Distinguished from other tissues, skeletal muscle ageCGs are more associated with expression, enriched near genes related to myofiber contraction, and closer to muscle-specific CTCF binding sites. Kidney-specific ageCGs are more increasingly methylated compared to other tissues as measured by affiliation with kidney-specific expressed genes. Underlying chromatin features also mark common and tissue-specific age effects reflective of poised and active chromatin states, respectively. In contrast with decreasingly methylated ageCGs, increasingly methylated ageCGs are also generally further from CTCF binding sites and enriched within lamina associated domains.

Conclusions

Our data identified common and tissue-specific DNA methylation changes with age that are reflective of CpG landscape and suggests both common and unique alterations within human tissues. Our findings also indicate that a simple epigenetic drift model is insufficient to explain all age-related changes in DNA methylation.     

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.     

Longitudinal,genome-scale analysis of DNA methylation in twins from birth to 18 months of age reveals rapid epigenetic change in early life and pair-specific effects of discordance

Background

The extent to which development- and age-associated epigenetic changes are influenced by genetic, environmental and stochastic factors remains to be discovered. Twins provide an ideal model with which to investigate these influences but previous cross-sectional twin studies provide contradictory evidence of within-pair epigenetic drift over time. Longitudinal twin studies can potentially address this discrepancy.

Results

In a pilot, genome-scale study of DNA from buccal epithelium, a relatively homogeneous tissue, we show that one-third of the CpGs assayed show dynamic methylation between birth and 18 months. Although all classes of annotated genomic regions assessed show an increase in DNA methylation over time, probes located in intragenic regions, enhancers and low-density CpG promoters are significantly over-represented, while CpG islands and high-CpG density promoters are depleted among the most dynamic probes. Comparison of co-twins demonstrated that within-pair drift in DNA methylation in our cohort is specific to a subset of pairs, who show more differences at 18 months. The rest of the pairs show either minimal change in methylation discordance, or more similar, converging methylation profiles at 18 months. As with age-associated regions, sites that change in their level of within-pair discordance between birth and 18 months are enriched in genes involved in development, but the average magnitude of change is smaller than for longitudinal change.

Conclusions

Our findings suggest that DNA methylation in buccal epithelium is influenced by non-shared stochastic and environmental factors that could reflect a degree of epigenetic plasticity within an otherwise constrained developmental program.     

Reducing the risk of false discovery enabling identification of biologically significant genome-wide methylation status using the HumanMethylation450 array

Background

The Illumina HumanMethylation450 BeadChip (HM450K) measures the DNA methylation of 485,512 CpGs in the human genome. The technology relies on hybridization of genomic fragments to probes on the chip. However, certain genomic factors may compromise the ability to measure methylation using the array such as single nucleotide polymorphisms (SNPs), small insertions and deletions (INDELs), repetitive DNA, and regions with reduced genomic complexity. Currently, there is no clear method or pipeline for determining which of the probes on the HM450K bead array should be retained for subsequent analysis in light of these issues.

Results

We comprehensively assessed the effects of SNPs, INDELs, repeats and bisulfite induced reduced genomic complexity by comparing HM450K bead array results with whole genome bisulfite sequencing. We determined which CpG probes provided accurate or noisy signals. From this, we derived a set of high-quality probes that provide unadulterated measurements of DNA methylation.

Conclusions

Our method significantly reduces the risk of false discoveries when using the HM450K bead array, while maximising the power of the array to detect methylation status genome-wide. Additionally, we demonstrate the utility of our method through extraction of biologically relevant epigenetic changes in prostate cancer.

Electronic supplementary material

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

Meta-analysis of human methylation data for evidence of sex-specific autosomal patterns

Background

Several individual studies have suggested that autosomal CpG methylation differs by sex both in terms of individual CpG sites and global autosomal CpG methylation. However, these findings have been inconsistent and plagued by spurious associations due to the cross reactivity of CpG probes on commercial microarrays. We collectively analysed 76 published studies (n = 6,795) for sex-associated differences in both autosomal and sex chromosome CpG sites.

Results

Overall autosomal methylation profiles varied substantially by study, and we encountered substantial batch effects. We accounted for these by conducting random effects meta-analysis for individual autosomal CpG methylation associations. After excluding non-specific probes, we found 184 autosomal CpG sites differentially methylated by sex after correction for multiple testing. In line with previous studies, average beta differences were small. Many of the most significantly associated CpG probes were new. Of note was differential CpG methylation in the promoters of genes thought to be involved in spermatogenesis and male fertility, such as SLC9A2, SPESP1, CRISP2, and NUPL1. Pathway analysis revealed overrepresentation of genes differentially methylated by sex in several broad Gene Ontology biological processes, including RNA splicing and DNA repair.

Conclusions

This study represents a comprehensive analysis of sex-specific methylation patterns. We demonstrate the existence of sex-specific methylation profiles and report a large number of novel DNA methylation differences in autosomal CpG sites between sexes.

Electronic supplementary material

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

Genome-wide signatures of differential DNA methylation in pediatric acute lymphoblastic leukemia

Background

Although aberrant DNA methylation has been observed previously in acute lymphoblastic leukemia (ALL), the patterns of differential methylation have not been comprehensively determined in all subtypes of ALL on a genome-wide scale. The relationship between DNA methylation, cytogenetic background, drug resistance and relapse in ALL is poorly understood.

Results

We surveyed the DNA methylation levels of 435,941 CpG sites in samples from 764 children at diagnosis of ALL and from 27 children at relapse. This survey uncovered four characteristic methylation signatures. First, compared with control blood cells, the methylomes of ALL cells shared 9,406 predominantly hypermethylated CpG sites, independent of cytogenetic background. Second, each cytogenetic subtype of ALL displayed a unique set of hyper- and hypomethylated CpG sites. The CpG sites that constituted these two signatures differed in their functional genomic enrichment to regions with marks of active or repressed chromatin. Third, we identified subtype-specific differential methylation in promoter and enhancer regions that were strongly correlated with gene expression. Fourth, a set of 6,612 CpG sites was predominantly hypermethylated in ALL cells at relapse, compared with matched samples at diagnosis. Analysis of relapse-free survival identified CpG sites with subtype-specific differential methylation that divided the patients into different risk groups, depending on their methylation status.

Conclusions

Our results suggest an important biological role for DNA methylation in the differences between ALL subtypes and in their clinical outcome after treatment.     

An Epigenetic Signature in Peripheral Blood Predicts Active Ovarian Cancer

Background

Recent studies have shown that DNA methylation (DNAm) markers in peripheral blood may hold promise as diagnostic or early detection/risk markers for epithelial cancers. However, to date no study has evaluated the diagnostic and predictive potential of such markers in a large case control cohort and on a genome-wide basis.

Principal Findings

By performing genome-wide DNAm profiling of a large ovarian cancer case control cohort, we here demonstrate that active ovarian cancer has a significant impact on the DNAm pattern in peripheral blood. Specifically, by measuring the methylation levels of over 27,000 CpGs in blood cells from 148 healthy individuals and 113 age-matched pre-treatment ovarian cancer cases, we derive a DNAm signature that can predict the presence of active ovarian cancer in blind test sets with an AUC of 0.8 (95% CI (0.74–0.87)). We further validate our findings in another independent set of 122 post-treatment cases (AUC = 0.76 (0.72–0.81)). In addition, we provide evidence for a significant number of candidate risk or early detection markers for ovarian cancer. Furthermore, by comparing the pattern of methylation with gene expression data from major blood cell types, we here demonstrate that age and cancer elicit common changes in the composition of peripheral blood, with a myeloid skewing that increases with age and which is further aggravated in the presence of ovarian cancer. Finally, we show that most cancer and age associated methylation variability is found at CpGs located outside of CpG islands.

Significance

Our results underscore the potential of DNAm profiling in peripheral blood as a tool for detection or risk-prediction of epithelial cancers, and warrants further in-depth and higher CpG coverage studies to further elucidate this role.     

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.     

Adult monozygotic twins discordant for intra-uterine growth have indistinguishable genome-wide DNA methylation profiles

Background

Low birth weight is associated with an increased adult metabolic disease risk. It is widely discussed that poor intra-uterine conditions could induce long-lasting epigenetic modifications, leading to systemic changes in regulation of metabolic genes. To address this, we acquire genome-wide DNA methylation profiles from saliva DNA in a unique cohort of 17 monozygotic monochorionic female twins very discordant for birth weight. We examine if adverse prenatal growth conditions experienced by the smaller co-twins lead to long-lasting DNA methylation changes.

Results

Overall, co-twins show very similar genome-wide DNA methylation profiles. Since observed differences are almost exclusively caused by variable cellular composition, an original marker-based adjustment strategy was developed to eliminate such variation at affected CpGs. Among adjusted and unchanged CpGs 3,153 are differentially methylated between the heavy and light co-twins at nominal significance, of which 45 show sensible absolute mean β-value differences. Deep bisulfite sequencing of eight such loci reveals that differences remain in the range of technical variation, arguing against a reproducible biological effect. Analysis of methylation in repetitive elements using methylation-dependent primer extension assays also indicates no significant intra-pair differences.

Conclusions

Severe intra-uterine growth differences observed within these monozygotic twins are not associated with long-lasting DNA methylation differences in cells composing saliva, detectable with up-to-date technologies. Additionally, our results indicate that uneven cell type composition can lead to spurious results and should be addressed in epigenomic studies.