DNA hypermethylation and DNA hypomethylation is present at different loci in chronic kidney disease

Genetic risk factors for chronic kidney disease (CKD) are being identified through international collaborations. By comparison, epigenetic risk factors for CKD have only recently been considered using population-based approaches. DNA methylation is a major epigenetic modification that is associated with complex diseases, so we investigated methylome-wide loci for association with CKD. A total of 485,577 unique features were evaluated in 255 individuals with CKD (cases) and 152 individuals without evidence of renal disease (controls). Following stringent quality control, raw data were quantile normalized and β values calculated to reflect the methylation status at each site. The difference in methylation status was evaluated between cases and controls with resultant P values adjusted for multiple testing. Genes with significantly increased and decreased levels of DNA methylation were considered for biological relevance by functional enrichment analysis using KEGG pathways in Partek Genomics Suite. Twenty-three genes, where more than one CpG per loci was identified with P_adjusted < 10⁻⁸, demonstrated significant methylation changes associated with CKD and additional support for these associated loci was sought from published literature. Strong biological candidates for CKD that showed statistically significant differential methylation include CUX1, ELMO1, FKBP5, INHBA-AS1, PTPRN2, and PRKAG2 genes; several genes are differentially methylated in kidney tissue and RNA-seq supports a functional role for differential methylation in ELMO1 and PRKAG2 genes. This study reports the largest, most comprehensive, genome-wide quantitative evaluation of DNA methylation for association with CKD. Evidence confirming methylation sites influence development of CKD would stimulate research to identify epigenetic therapies that might be clinically useful for CKD.     

Distinctive patterns of DNA methylation associated with Parkinson disease

Parkinson disease (PD) is a multifactorial neurodegenerative disorder with high incidence in the elderly, where environmental and genetic factors are involved in etiology. In addition, epigenetic mechanisms, including deregulation of DNA methylation have been recently associated to PD. As accurate diagnosis cannot be achieved pre-mortem, identification of early pathological changes is crucial to enable therapeutic interventions before major neuropathological damage occurs. Here we investigated genome-wide DNA methylation in brain and blood samples from PD patients and observed a distinctive pattern of methylation involving many genes previously associated to PD, therefore supporting the role of epigenetic alterations as a molecular mechanism in neurodegeneration. Importantly, we identified concordant methylation alterations in brain and blood, suggesting that blood might hold promise as a surrogate for brain tissue to detect DNA methylation in PD and as a source for biomarker discovery.     

DNA hypermethylation and DNA hypomethylation is present at different loci in chronic kidney disease

Genetic risk factors for chronic kidney disease (CKD) are being identified through international collaborations. By comparison, epigenetic risk factors for CKD have only recently been considered using population-based approaches. DNA methylation is a major epigenetic modification that is associated with complex diseases, so we investigated methylome-wide loci for association with CKD. A total of 485,577 unique features were evaluated in 255 individuals with CKD (cases) and 152 individuals without evidence of renal disease (controls). Following stringent quality control, raw data were quantile normalized and β values calculated to reflect the methylation status at each site. The difference in methylation status was evaluated between cases and controls with resultant P values adjusted for multiple testing. Genes with significantly increased and decreased levels of DNA methylation were considered for biological relevance by functional enrichment analysis using KEGG pathways in Partek Genomics Suite. Twenty-three genes, where more than one CpG per loci was identified with P_adjusted < 10^?8, demonstrated significant methylation changes associated with CKD and additional support for these associated loci was sought from published literature. Strong biological candidates for CKD that showed statistically significant differential methylation include CUX1, ELMO1, FKBP5, INHBA-AS1, PTPRN2, and PRKAG2 genes; several genes are differentially methylated in kidney tissue and RNA-seq supports a functional role for differential methylation in ELMO1 and PRKAG2 genes. This study reports the largest, most comprehensive, genome-wide quantitative evaluation of DNA methylation for association with CKD. Evidence confirming methylation sites influence development of CKD would stimulate research to identify epigenetic therapies that might be clinically useful for CKD.     

Genome-wide DNA methylation profiling in rheumatoid arthritis identifies disease-associated methylation changes that are distinct to individual T- and B-lymphocyte populations

Changes to the DNA methylome have been described in patients with rheumatoid arthritis (RA). In previous work, we reported genome-wide methylation differences in T-lymphocyte and B-lymphocyte populations from healthy individuals. Now, using HumanMethylation450 BeadChips to interrogate genome-wide DNA methylation, we have determined disease-associated methylation changes in blood-derived T- and B-lymphocyte populations from 12 female patients with seropositive established RA, relative to 12 matched healthy individuals. Array data were analyzed using NIMBL software and bisulfite pyrosequencing was used to validate array candidates. Genome-wide DNA methylation, determined by analysis of LINE-1 sequences, revealed higher methylation in B-lymphocytes compared with T-lymphocytes (P ≤ 0.01), which is consistent with our findings in healthy individuals. Moreover, loci-specific methylation differences that distinguished T-lymphocytes from B-lymphocytes in healthy individuals were also apparent in RA patients. However, disease-associated methylation differences were also identified in RA. In these cases, we identified 509 and 252 CpGs in RA-derived T- and B-lymphocytes, respectively, that showed significant changes in methylation compared with their cognate healthy counterparts. Moreover, this included a restricted set of 32 CpGs in T-lymphocytes and 20 CpGs in B-lymphocytes (representing 15 and 10 genes, respectively, and including two, MGMT and CCS, that were common to both cell types) that displayed more substantial changes in methylation. These changes, apparent as hyper- or hypo-methylation, were independently confirmed by pyrosequencing analysis. Validation by pyrosequencing also revealed additional sites in some candidate genes that also displayed altered methylation in RA. In this first study of genome-wide DNA methylation in individual T- and B-lymphocyte populations in RA patients, we report disease-associated methylation changes that are distinct to each cell type and which support a role for discrete epigenetic regulation in this disease.     

Integrative Analysis of DNA Methylation and Gene Expression Data Identifies EPAS1 as a Key Regulator of COPD

Chronic Obstructive Pulmonary Disease (COPD) is a complex disease. Genetic, epigenetic, and environmental factors are known to contribute to COPD risk and disease progression. Therefore we developed a systematic approach to identify key regulators of COPD that integrates genome-wide DNA methylation, gene expression, and phenotype data in lung tissue from COPD and control samples. Our integrative analysis identified 126 key regulators of COPD. We identified EPAS1 as the only key regulator whose downstream genes significantly overlapped with multiple genes sets associated with COPD disease severity. EPAS1 is distinct in comparison with other key regulators in terms of methylation profile and downstream target genes. Genes predicted to be regulated by EPAS1 were enriched for biological processes including signaling, cell communications, and system development. We confirmed that EPAS1 protein levels are lower in human COPD lung tissue compared to non-disease controls and that Epas1 gene expression is reduced in mice chronically exposed to cigarette smoke. As EPAS1 downstream genes were significantly enriched for hypoxia responsive genes in endothelial cells, we tested EPAS1 function in human endothelial cells. EPAS1 knockdown by siRNA in endothelial cells impacted genes that significantly overlapped with EPAS1 downstream genes in lung tissue including hypoxia responsive genes, and genes associated with emphysema severity. Our first integrative analysis of genome-wide DNA methylation and gene expression profiles illustrates that not only does DNA methylation play a ‘causal’ role in the molecular pathophysiology of COPD, but it can be leveraged to directly identify novel key mediators of this pathophysiology.     

Mycoplasma CG- and GATC-specific DNA methyltransferases selectively and efficiently methylate the host genome and alter the epigenetic landscape in human cells

Aberrant DNA methylation is frequently observed in disease, including many cancer types, yet the underlying mechanisms remain unclear. Because germline and somatic mutations in the genes that are responsible for DNA methylation are infrequent in malignancies, additional mechanisms must be considered. Mycoplasmas spp., including Mycoplasma hyorhinis, efficiently colonize human cells and may serve as a vehicle for delivery of enzymatically active microbial proteins into the intracellular milieu. Here, we performed, for the first time, genome-wide and individual gene mapping of methylation marks generated by the M. hyorhinis CG- and GATC-specific DNA cytosine methyltransferases (MTases) in human cells. Our results demonstrated that, upon expression in human cells, MTases readily translocated to the cell nucleus. In the nucleus, MTases selectively and efficiently methylated the host genome at the DNA sequence sites free from pre-existing endogenous methylation, including those in a variety of cancer-associated genes. We also established that mycoplasma is widespread in colorectal cancers, suggesting that either the infection contributed to malignancy onset or, alternatively, that tumors provide a favorable environment for mycoplasma growth. In the human genome, ∼11% of GATC sites overlap with CGs (e.g., CGAT^mCG); therefore, the methylated status of these sites can be perpetuated by human DNMT1. Based on these results, we now suggest that the GATC-specific methylation represents a novel type of infection-specific epigenetic mark that originates in human cells with a previous exposure to infection. Overall, our findings unveil an entirely new panorama of interactions between the human microbiome and epigenome with a potential impact in disease etiology.     

DNA methylation analysis of Homeobox genes implicates HOXB7 hypomethylation as risk factor for neural tube defects

Neural tube defects (NTDs) are common birth defects of complex etiology. Though family- and population-based studies have confirmed a genetic component, the responsible genes for NTDs are still largely unknown. Based on the hypothesis that folic acid prevents NTDs by stimulating methylation reactions, epigenetic factors, such as DNA methylation, are predicted to be involved in NTDs. Homeobox (HOX) genes play a role in spinal cord development and are tightly regulated in a spatiotemporal and collinear manner, partly by epigenetic modifications. We have quantified DNA methylation for the different HOX genes by subtracting values from a genome-wide methylation analysis using leukocyte DNA from 10 myelomeningocele (MMC) patients and 6 healthy controls. From the 1575 CpGs profiled for the 4 HOX clusters, 26 CpGs were differentially methylated (P-value < 0.05; β-difference > 0.05) between MMC patients and controls. Seventy-seven percent of these CpGs were located in the HOXA and HOXB clusters, with the most profound difference for 3 CpGs within the HOXB7 gene body. A validation case-control study including 83 MMC patients and 30 unrelated healthy controls confirmed a significant association between MMC and HOXB7 hypomethylation (-14.4%; 95% CI: 11.9–16.9%; P-value < 0.0001) independent of the MTHFR 667C>T genotype. Significant HOXB7 hypomethylation was also present in 12 unaffected siblings, each related to a MMC patient, suggestive of an epigenetic change induced by the mother. The inclusion of a neural tube formation model using zebrafish showed that Hoxb7a overexpression but not depletion resulted in deformed body axes with dysmorphic neural tube formation. Our results implicate HOXB7 hypomethylation as risk factor for NTDs and highlight the importance for future genome-wide DNA methylation analyses without preselecting candidate pathways.     

Genome-wide DNA methylation profiling in rheumatoid arthritis identifies disease-associated methylation changes that are distinct to individual T- and B-lymphocyte populations

Changes to the DNA methylome have been described in patients with rheumatoid arthritis (RA). In previous work, we reported genome-wide methylation differences in T-lymphocyte and B-lymphocyte populations from healthy individuals. Now, using HumanMethylation450 BeadChips to interrogate genome-wide DNA methylation, we have determined disease-associated methylation changes in blood-derived T- and B-lymphocyte populations from 12 female patients with seropositive established RA, relative to 12 matched healthy individuals. Array data were analyzed using NIMBL software and bisulfite pyrosequencing was used to validate array candidates. Genome-wide DNA methylation, determined by analysis of LINE-1 sequences, revealed higher methylation in B-lymphocytes compared with T-lymphocytes (P ≤ 0.01), which is consistent with our findings in healthy individuals. Moreover, loci-specific methylation differences that distinguished T-lymphocytes from B-lymphocytes in healthy individuals were also apparent in RA patients. However, disease-associated methylation differences were also identified in RA. In these cases, we identified 509 and 252 CpGs in RA-derived T- and B-lymphocytes, respectively, that showed significant changes in methylation compared with their cognate healthy counterparts. Moreover, this included a restricted set of 32 CpGs in T-lymphocytes and 20 CpGs in B-lymphocytes (representing 15 and 10 genes, respectively, and including two, MGMT and CCS, that were common to both cell types) that displayed more substantial changes in methylation. These changes, apparent as hyper- or hypo-methylation, were independently confirmed by pyrosequencing analysis. Validation by pyrosequencing also revealed additional sites in some candidate genes that also displayed altered methylation in RA. In this first study of genome-wide DNA methylation in individual T- and B-lymphocyte populations in RA patients, we report disease-associated methylation changes that are distinct to each cell type and which support a role for discrete epigenetic regulation in this disease.     

Methylome sequencing for fibrolamellar hepatocellular carcinoma depicts distinctive features

With the goal of studying epigenetic alterations in fibrolamellar hepatocellular carcinoma (FLC) and establish an associated DNA methylation signature, we analyzed LINE-1 methylation in a cohort of FLC and performed next-generation sequencing of DNA methylation in a training set of pure-FLCs and non-cirrhotic hepatocellular carcinomas (nc-HCC). DNA methylation was correlated with gene expression. Furthermore, we established and validated an epigenetic signature differentiating pure-FLC from other HCCs. LINE-1 methylation correlated with shorter recurrence-free survival and overall survival in resected pure-FLC patients. Unsupervised clustering using CG sites located in islands distinguished pure-FLC from nc-HCC. Major DNA methylation changes occurred outside promoters, mainly in gene bodies and intergenic regions located in the vicinity of liver developmental genes (i.e., SMARCA4 and RXRA). Partially methylated domains were more prone to DNA methylation changes. Furthermore, we identified several putative tumor suppressor genes (e.g., DLEU7) and oncogenes (e.g., DUSP4). While ∼70% of identified gene promoters gaining methylation were marked by bivalent histone marks (H3K4me3/H3K27me3) in embryonic stem cells, ∼70% of those losing methylation were marked by H3K4me3. Finally, we established a pure FLC DNA methylation signature and validated it in an independent dataset. Our analysis reveals a distinct epigenetic signature of pure FLC as compared to nc-HCC, with DNA methylation changes occurring in the vicinity of liver developmental genes. These data suggest new options for targeting FLC based on cancer epigenome aberrations.     

Allele-specific DNA methylation analyses associated with siRNAs in Arabidopsis hybrids

Accumulating evidence has suggested that epigenetic marks including DNA methylation,small RNA and histone modification may involve hybrid vigor in plants.However,knowledge about how epigenetic marks in hybrids regulate gene expression is still limited.Based on genome-wide DNA methylation landscapes of Arabidopsis thaliana Ler and C24 ecotypes and their reciprocal F1 hybrids which were obtained in our previous work,we analyzed allele-specific DNA methylation and distinguished cis-and trans-regulated DNA methylation in hybrids.Our study indicated that both cis-and trans-regulated DNA methylation played roles in hybrids,when cis-regulation played a major role in CG methylation and trans-regulation played major roles in CHG and CHH methylation.In addition,we observed correlations between trans-regulated DNA methylation and siRNA densities.Enriched siRNA regions were significantly concurrent with highly trans-regulated DNA methylation regions.Our results illustrated DNA methylation regulation patterns integrated with siRNAs in Arabidopsis hybrids,and shed light on understanding the mechanism of epigenetic reprogramming for hybrid vigor.     

MB-COMT promoter DNA methylation is associated with working-memory processing in schizophrenia patients and healthy controls

Many genetic studies report mixed results both for the associations between COMT polymorphisms and schizophrenia and for the effects of COMT variants on common intermediate phenotypes of the disorder. Reasons for this may include small genetic effect sizes and the modulation of environmental influences. To improve our understanding of the role of COMT in the disease etiology, we investigated the effect of DNA methylation in the MB-COMT promoter on neural activity in the dorsolateral prefrontal cortex during working memory processing as measured by fMRI - an intermediate phenotype for schizophrenia. Imaging and epigenetic data were measured in 102 healthy controls and 82 schizophrenia patients of the Mind Clinical Imaging Consortium (MCIC) study of schizophrenia. Neural activity during the Sternberg Item Recognition Paradigm was acquired with either a 3T Siemens Trio or 1.5T Siemens Sonata and analyzed using the FMRIB Software Library (FSL). DNA methylation measurements were derived from cryo-conserved blood samples. We found a positive association between MB-COMT promoter methylation and neural activity in the left dorsolateral prefrontal cortex in a model using a region-of-interest approach and could confirm this finding in a whole-brain model. This effect was independent of disease status. Analyzing the effect of MB-COMT promoter DNA methylation on a neuroimaging phenotype can provide further evidence for the importance of COMT and epigenetic risk mechanisms in schizophrenia. The latter may represent trans-regulatory or environmental risk factors that can be measured using brain-based intermediate phenotypes.     

Age-related epigenome-wide DNA methylation and hydroxymethylation in longitudinal mouse blood

DNA methylation at cytosine-phosphate-guanine (CpG) dinucleotides changes as a function of age in humans and animal models, a process that may contribute to chronic disease development. Recent studies have investigated the role of an oxidized form of DNA methylation – 5-hydroxymethylcytosine (5hmC) – in the epigenome, but its contribution to age-related DNA methylation remains unclear. We tested the hypothesis that 5hmC changes with age, but in a direction opposite to 5-methylcytosine (5mC), potentially playing a distinct role in aging. To characterize epigenetic aging, genome-wide 5mC and 5hmC were measured in longitudinal blood samples (2, 4, and 10 months of age) from isogenic mice using two sequencing methods – enhanced reduced representation bisulfite sequencing and hydroxymethylated DNA immunoprecipitation sequencing. Examining the epigenome by age, we identified 28,196 unique differentially methylated CpGs (DMCs) and 8,613 differentially hydroxymethylated regions (DHMRs). Mouse blood showed a general pattern of epigenome-wide hypermethylation and hypo-hydroxymethylation with age. Comparing age-related DMCs and DHMRs, 1,854 annotated genes showed both differential 5mC and 5hmC, including one gene – Nfic – at five CpGs in the same 250 bp chromosomal region. At this region, 5mC and 5hmC levels both decreased with age. Reflecting these age-related epigenetic changes, Nfic RNA expression in blood decreased with age, suggesting that age-related regulation of this gene may be driven by 5hmC, not canonical DNA methylation. Combined, our genome-wide results show age-related differential 5mC and 5hmC, as well as some evidence that changes in 5hmC may drive age-related DNA methylation and gene expression.