Global Analysis of DNA Methylation Variation in Adipose Tissue from Twins Reveals Links to Disease-Associated Variants in Distal Regulatory Elements

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h²_median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.     

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.     

Novel Epigenetic Changes Unveiled by Monozygotic Twins Discordant for Smoking Habits

Exposure to cigarette smoking affects the epigenome and could increase the risk of developing diseases such as cancer and cardiovascular disorders. Changes in DNA methylation associated with smoking may help to identify molecular pathways that contribute to disease etiology. Previous studies are not completely concordant in the identification of differentially methylated regions in the DNA of smokers. We performed an epigenome-wide DNA methylation study in a group of monozygotic (MZ) twins discordant for smoking habits to determine the effect of smoking on DNA methylation. As MZ twins are considered genetically identical, this model allowed us to identify smoking-related DNA methylation changes independent from genetic components. We investigated the whole blood genome-wide DNA methylation profiles in 20 MZ twin pairs discordant for smoking habits by using the Illumina HumanMethylation450 BeadChip. We identified 22 CpG sites that were differentially methylated between smoker and non-smoker MZ twins by intra-pair analysis. We confirmed eight loci already described by other groups, located in AHRR, F2RL3, MYOG1 genes, at 2q37.1 and 6p21.33 regions, and also identified several new loci. Moreover, pathway analysis showed an enrichment of genes involved in GTPase regulatory activity. Our study confirmed the evidence of smoking-related DNA methylation changes, emphasizing that well-designed MZ twin models can aid the discovery of novel DNA methylation signals, even in a limited sample population.     

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.     

Genome-wide DNA methylation studies suggest distinct DNA methylation patterns in pediatric embryonal and alveolar rhabdomyosarcomas

Rhabdomyosarcoma is the most common soft-tissue sarcoma in children. While cytogenetic abnormalities have been well characterized in this disease, aberrant epigenetic events such as DNA hypermethylation have not been described in genome-wide studies. We have analyzed the methylation status of 25,500 promoters in normal skeletal muscle, and in cell lines and tumor samples of embryonal and alveolar rhabdomyosarcoma from pediatric patients. We identified over 1,900 CpG islands that are hypermethylated in rhabdomyosarcomas relative to skeletal muscle. Genes involved in tissue development, differentiation, and oncogenesis such as DNAJA4, HES5, IRX1, BMP8A, GATA4, GATA6, ALX3, and P4HTM were hypermethylated in both RMS cell lines and primary samples, implicating aberrant DNA methylation in the pathogenesis of rhabdomyosarcoma. Furthermore, cluster analysis revealed embryonal and alveolar subtypes had distinct DNA methylation patterns, with the alveolar subtype being enriched in DNA hypermethylation of polycomb target genes. These results suggest that DNA methylation signatures may aid in the diagnosis and risk stratification of pediatric rhabdomyosarcoma and help identify new targets for therapy.     

GATC motifs may alter the conformation of DNA depending on sequence context and N6-adenine methylation status: possible implications for DNA-protein recognition

As part of our analysis of the role of a uniquely clustered set of dam methylation sites (the motif GATC) within the origin of DNA replication in Escherichia coli, we have studied the effect of GATCs in various methylation states on the intrinsic curvature of DNA. We have designed a set of DNA linkers and used commercially available linkers containing GATC motifs. The linkers were ligated and the electrophoretic mobility of the resulting multimers in different states of methylation was tested relative to reference fragments. We report that properly phased GATCs in certain sequence environments modulate DNA curvature and that these effects may be enhanced by N⁶-adenine methylation of the GATCs. These structural alterations may in turn affect DNA-protein interactions, especially those involving proteins that rely on both primary sequence and structure for recognition. We present an example, where introduction of a GATC within an integration host factor (IHF) binding site, which does not alter the consensus sequence, reduces the binding affinity of the protein for the modified site.     

Development of a Quantitative Methylation-Specific Polymerase Chain Reaction Method for Monitoring Beta Cell Death in Type 1 Diabetes

DNA methylation is a mechanism by which cells control gene expression, and cell-specific genes often exhibit unique patterns of DNA methylation. We previously reported that the mouse insulin-2 gene (Ins2) promoter has three potential methylation (CpG) sites, all of which are unmethylated in insulin-producing cells but methylated in other tissues. In this study we examined Ins2 exon 2 and found a similar tissue-specific methylation pattern. These methylation patterns can differentiate between DNA from insulin-producing beta cells and other tissues. We hypothesized that damaged beta cells release their DNA into circulation at the onset of type 1 diabetes mellitus (T1DM) and sought to develop a quantitative methylation-specific polymerase chain reaction (qMSP) assay for circulating beta cell DNA to monitor the loss of beta cells. Methylation-specific primers were designed to interrogate two or more CpG in the same assay. The cloned mouse Ins2 gene was methylated in vitro and used for development of the qMSP assay. We found the qMSP method to be sensitive and specific to differentiate between insulin-producing cells and other tissues with a detection limit of 10 copies in the presence of non-specific genomic DNA background. We also compared different methods for data analysis and found that the Relative Expression Ratio method is the most robust method since it incorporates both a reference value to normalize day-to-day variability as well as PCR reaction efficiencies to normalize between the methylation-specific and bisulfite-specific components of the calculations. The assay was applied in the streptozotocin-treated diabetic mouse model and detected a significant increase in circulating beta cell DNA before the rise in blood glucose level. These results demonstrate that this qMSP assay can be used for monitoring circulating DNA from insulin-producing cells, which will provide the basis for development of assays to detect beta cell destruction in early T1DM.     

Multilocus loss of DNA methylation in individuals with mutations in the histone H3 Lysine 4 Demethylase KDM5C

Background

A number of neurodevelopmental syndromes are caused by mutations in genes encoding proteins that normally function in epigenetic regulation. Identification of epigenetic alterations occurring in these disorders could shed light on molecular pathways relevant to neurodevelopment.

Results

Using a genome-wide approach, we identified genes with significant loss of DNA methylation in blood of males with intellectual disability and mutations in the X-linked KDM5C gene, encoding a histone H3 lysine 4 demethylase, in comparison to age/sex matched controls. Loss of DNA methylation in such individuals is consistent with known interactions between DNA methylation and H3 lysine 4 methylation. Further, loss of DNA methylation at the promoters of the three top candidate genes FBXL5, SCMH1, CACYBP was not observed in more than 900 population controls. We also found that DNA methylation at these three genes in blood correlated with dosage of KDM5C and its Y-linked homologue KDM5D. In addition, parallel sex-specific DNA methylation profiles in brain samples from control males and females were observed at FBXL5 and CACYBP.

Conclusions

We have, for the first time, identified epigenetic alterations in patient samples carrying a mutation in a gene involved in the regulation of histone modifications. These data support the concept that DNA methylation and H3 lysine 4 methylation are functionally interdependent. The data provide new insights into the molecular pathogenesis of intellectual disability. Further, our data suggest that some DNA methylation marks identified in blood can serve as biomarkers of epigenetic status in the brain.     

Epigenetic instability at imprinting control regions in a KrasG12D-induced T-cell neoplasm

Although aberrant DNA methylation within imprinted domains has been reported in a variety of neoplastic diseases, it remains largely uncharacterized in the context of carcinogenesis. In this study, we induced T-cell lymphoma in mice by employing a breeding scheme involving mouse strains, LSL-Kras^G12D and MMTV-Cre. We then systematically surveyed imprinted domains for DNA methylation changes during tumor progression using combined bisulfite restriction analysis and NGS-based bisulfite sequencing. We detected hyper- or hypo-methylation at the imprinting control regions (ICRs) of the Dlk1, Peg10, Peg3, Grb10, and Gnas domains. These DNA methylation changes at ICRs were more prevalent and consistent than those observed at the promoter regions of well-known tumor suppressors, such as Mgmt, Fhit, and Mlh1. Thus, the changes observed at these imprinted domains are the outcome of isolated incidents affecting DNA methylation settings. Within imprinted domains, DNA methylation changes tend to be restricted to ICRs as nearby somatic differentially methylated regions and promoter regions experience no change. Furthermore, detailed analyses revealed that small cis-regulatory elements within ICRs tend to be resistant to DNA methylation changes, suggesting potential protection by unknown trans-factors. Overall, this study demonstrates that DNA methylation changes at ICRs are dynamic during carcinogenesis and advocates that detection of aberrant DNA methylation at ICRs may serve as a biomarker to enhance diagnostic procedures.     

Aberrant DNA methylation is associated with disease progression, resistance to imatinib and shortened survival in chronic myelogenous leukemia

The epigenetic impact of DNA methylation in chronic myelogenous leukemia (CML) is not completely understood. To elucidate its role we analyzed 120 patients with CML for methylation of promoter-associated CpG islands of 10 genes. Five genes were identified by DNA methylation screening in the K562 cell line and 3 genes in patients with myeloproliferative neoplasms. The CDKN2B gene was selected for its frequent methylation in myeloid malignancies and ABL1 as the target of BCR-ABL translocation. Thirty patients were imatinib-naïve (mostly treated by interferon-alpha before the imatinib era), 30 were imatinib-responsive, 50 were imatinib-resistant, and 10 were imatinib-intolerant. We quantified DNA methylation by bisulfite pyrosequencing. The average number of methylated genes was 4.5 per patient in the chronic phase, increasing significantly to 6.2 in the accelerated and 6.4 in the blastic phase. Higher numbers of methylated genes were also observed in patients resistant or intolerant to imatinib. These patients also showed almost exclusive methylation of a putative transporter OSCP1. Abnormal methylation of a Src suppressor gene PDLIM4 was associated with shortened survival independently of CML stage and imatinib responsiveness. We conclude that aberrant DNA methylation is associated with CML progression and that DNA methylation could be a marker associated with imatinib resistance. Finally, DNA methylation of PDLIM4 may help identify a subset of CML patients that would benefit from treatment with Src/Abl inhibitors.     

The Identification of Specific Methylation Patterns across Different Cancers

Abnormal DNA methylation is known as playing an important role in the tumorgenesis. It is helpful for distinguishing the specificity of diagnosis and therapeutic targets for cancers based on characteristics of DNA methylation patterns across cancers. High throughput DNA methylation analysis provides the possibility to comprehensively filter the epigenetics diversity across various cancers. We integrated whole-genome methylation data detected in 798 samples from seven cancers. The hierarchical clustering revealed the existence of cancer-specific methylation pattern. Then we identified 331 differentially methylated genes across these cancers, most of which (266) were specifically differential methylation in unique cancer. A DNA methylation correlation network (DMCN) was built based on the methylation correlation between these genes. It was shown the hubs in the DMCN were inclined to cancer-specific genes in seven cancers. Further survival analysis using the part of genes in the DMCN revealed high-risk group and low-risk group were distinguished by seven biomarkers (PCDHB15, WBSCR17, IGF1, GYPC, CYGB, ACTG2, and PRRT1) in breast cancer and eight biomarkers (ZBTB32, OR51B4, CCL8, TMEFF2, SALL3, GPSM1, MAGEA8, and SALL1) in colon cancer, respectively. At last, a protein-protein interaction network was introduced to verify the biological function of differentially methylated genes. It was shown that MAP3K14, PTN, ACVR1 and HCK sharing different DNA methylation and gene expression across cancers were relatively high degree distribution in PPI network. The study suggested that not only the identified cancer-specific genes provided reference for individual treatment but also the relationship across cancers could be explained by differential DNA methylation.     

Array-based profiling of reference-independent methylation status (aPRIMES) identifies frequent promoter methylation and consecutive downregulation of ZIC2 in pediatric medulloblastoma

Existing microarray-based approaches for screening of DNA methylation are hampered by a number of shortcomings, such as the introduction of bias by DNA copy-number imbalances in the test genome and negligence of tissue-specific methylation patterns. We developed a method designated array-based profiling of reference-independent methylation status (aPRIMES) that allows the detection of direct methylation status rather than relative methylation. Array-PRIMES is based on the differential restriction and competitive hybridization of methylated and unmethylated DNA by methylation-specific and methylation-sensitive restriction enzymes, respectively. We demonstrate the accuracy of aPRIMES in detecting the methylation status of CpG islands for different states of methylation. Application of aPRIMES to the DNA from desmoplastic medulloblastomas of monozygotic twins showed strikingly similar methylation profiles. Additional analysis of 18 sporadic medulloblastomas revealed an overall correlation between highly methylated tumors and poor clinical outcome and identified ZIC2 as a frequently methylated gene in pediatric medulloblastoma.     

DNA Hypomethylation Affects Cancer-Related Biological Functions and Genes Relevant in Neuroblastoma Pathogenesis

Neuroblastoma (NB) pathogenesis has been reported to be closely associated with numerous genetic alterations. However, underlying DNA methylation patterns have not been extensively studied in this developmental malignancy. Here, we generated microarray-based DNA methylation profiles of primary neuroblastic tumors. Stringent supervised differential methylation analyses allowed us to identify epigenetic changes characteristic for NB tumors as well as for clinical and biological subtypes of NB. We observed that gene-specific loss of DNA methylation is more prevalent than promoter hypermethylation. Remarkably, such hypomethylation affected cancer-related biological functions and genes relevant to NB pathogenesis such as CCND1, SPRR3, BTC, EGF and FGF6. In particular, differential methylation in CCND1 affected mostly an evolutionary conserved functionally relevant 3′ untranslated region, suggesting that hypomethylation outside promoter regions may play a role in NB pathogenesis. Hypermethylation targeted genes involved in cell development and proliferation such as RASSF1A, POU2F2 or HOXD3, among others. The results derived from this study provide new candidate epigenetic biomarkers associated with NB as well as insights into the molecular pathogenesis of this tumor, which involves a marked gene-specific hypomethylation.     

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.     

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.     

Deregulated Expression of SRC,LYN and CKB Kinases by DNA Methylation and Its Potential Role in Gastric Cancer Invasiveness and Metastasis

Kinases are downstream modulators and effectors of several cellular signaling cascades and play key roles in the development of neoplastic disease. In this study, we aimed to evaluate SRC, LYN and CKB protein and mRNA expression, as well as their promoter methylation, in gastric cancer. We found elevated expression of SRC and LYN kinase mRNA and protein but decreased levels of CKB kinase, alterations that may have a role in the invasiveness and metastasis of gastric tumors. Expression of the three studied kinases was also associated with MYC oncogene expression, a possible biomarker for gastric cancer. To understand the mechanisms that regulate the expression of these genes, we evaluated the DNA promoter methylation of the three kinases. We found that reduced SRC and LYN methylation and increased CKB methylation was associated with gastric cancer. The reduced SRC and LYN methylation was associated with increased levels of mRNA and protein expression, suggesting that DNA methylation is involved in regulating the expression of these kinases. Conversely, reduced CKB methylation was observed in samples with reduced mRNA and protein expression, suggesting CKB expression was found to be only partly regulated by DNA methylation. Additionally, we found that alterations in the DNA methylation pattern of the three studied kinases were also associated with the gastric cancer onset, advanced gastric cancer, deeper tumor invasion and the presence of metastasis. Therefore, SRC, LYN and CKB expression or DNA methylation could be useful markers for predicting tumor progression and targeting in anti-cancer strategies.     

Differential DNA methylation identified in the blood and retina of AMD patients

Age-related macular degeneration (AMD) is a major cause of blindness in the western world. While genetic studies have linked both common and rare variants in genes involved in regulation of the complement system to increased risk of development of AMD, environmental factors, such as smoking and nutrition, can also significantly affect the risk of developing the disease and the rate of disease progression. Since epigenetics has been implicated in mediating, in part, the disease risk associated with some environmental factors, we investigated a possible epigenetic contribution to AMD. We performed genome-wide DNA methylation profiling of blood from AMD patients and controls. No differential methylation site reached genome-wide significance; however, when epigenetic changes in and around known GWAS-defined AMD risk loci were explored, we found small but significant DNA methylation differences in the blood of neovascular AMD patients near age-related maculopathy susceptibility 2 (ARMS2), a top-ranked GWAS locus preferentially associated with neovascular AMD. The methylation level of one of the CpG sites significantly correlated with the genotype of the risk SNP rs10490924, suggesting a possible epigenetic mechanism of risk. Integrating genome-wide DNA methylation analysis of retina samples with and without AMD together with blood samples, we further identified a consistent, replicable change in DNA methylation in the promoter region of protease serine 50 (PRSS50). These methylation changes may identify sites in novel genes that are susceptible to non-genetic factors known to contribute to AMD development and progression.