Genome-wide Methyl-Seq analysis of blood-brain targets of glucocorticoid exposure

Chronic exposure to glucocorticoids (GCs) can lead to psychiatric complications through epigenetic mechanisms such as DNA methylation (DNAm). We sought to determine whether epigenetic changes in a peripheral tissue can serve as a surrogate for those in a relatively inaccessible tissue such as the brain. DNA extracted from the hippocampus and blood of mice treated with GCs or vehicle solution was assayed using a genome-wide DNAm platform (Methyl-Seq) to identify differentially methylated regions (DMRs) induced by GC treatment. We observed that ～70% of the DMRs in both tissues lost methylation following GC treatment. Of the 3,095 DMRs that mapped to the same genes in both tissues, 1,853 DMRs underwent DNAm changes in the same direction. Interestingly, only 209 DMRs (<7%) overlapped in genomic coordinates between the 2 tissues, suggesting tissue-specific differences in GC-targeted loci. Pathway analysis showed that the DMR-associated genes were members of pathways involved in metabolism, immune function, and neurodevelopment. Also, changes in cell type composition of blood and brain were examined by fluorescence-activated cell sorting. Separation of the cortex into neuronal and non-neuronal fractions and the leukocytes into T-cells, B-cells, and neutrophils showed that GC-induced methylation changes primarily occurred in neurons and T-cells, with the blood tissue also undergoing a shift in the proportion of constituent cell types while the proportion of neurons and glia in the brain remained stable. From the current pilot study, we found that despite tissue-specific epigenetic changes and cellular heterogeneity, blood can serve as a surrogate for GC-induced changes in the brain.     

Aberrant signature methylome by DNMT1 hot spot mutation in hereditary sensory and autonomic neuropathy 1E

DNA methyltransferase 1 (DNMT1) is essential for DNA methylation, gene regulation and chromatin stability. We previously discovered DNMT1 mutations cause hereditary sensory and autonomic neuropathy type 1 with dementia and hearing loss (HSAN1E; OMIM 614116). HSAN1E is the first adult-onset neurodegenerative disorder caused by a defect in a methyltransferase gene. HSAN1E patients appear clinically normal until young adulthood, then begin developing the characteristic symptoms involving central and peripheral nervous systems. Some HSAN1E patients also develop narcolepsy and it has recently been suggested that HSAN1E is allelic to autosomal dominant cerebellar ataxia, deafness, with narcolepsy (ADCA-DN; OMIM 604121), which is also caused by mutations in DNMT1. A hotspot mutation Y495C within the targeting sequence domain of DNMT1 has been identified among HSAN1E patients. The mutant DNMT1 protein shows premature degradation and reduced DNA methyltransferase activity. Herein, we investigate genome-wide DNA methylation at single-base resolution through whole-genome bisulfite sequencing of germline DNA in 3 pairs of HSAN1E patients and their gender- and age-matched siblings. Over 1 billion 75-bp single-end reads were generated for each sample. In the 3 affected siblings, overall methylation loss was consistently found in all chromosomes with X and 18 being most affected. Paired sample analysis identified 564,218 differentially methylated CpG sites (DMCs; P < 0.05), of which 300?134 were intergenic and 264?084 genic CpGs. Hypomethylation was predominant in both genic and intergenic regions, including promoters, exons, most CpG islands, L1, L2, Alu, and satellite repeats and simple repeat sequences. In some CpG islands, hypermethylated CpGs outnumbered hypomethylated CpGs. In 201 imprinted genes, there were more DMCs than in non-imprinted genes and most were hypomethylated. Differentially methylated region (DMR) analysis identified 5649 hypomethylated and 1872 hypermethylated regions. Importantly, pathway analysis revealed 1693 genes associated with the identified DMRs were highly associated in diverse neurological disorders and NAD+/NADH metabolism pathways is implicated in the pathogenesis. Our results provide novel insights into the epigenetic mechanism of neurodegeneration arising from a hotspot DNMT1 mutation and reveal pathways potentially important in a broad category of neurological and psychological disorders.     

Aberrant signature methylome by DNMT1 hot spot mutation in hereditary sensory and autonomic neuropathy 1E

DNA methyltransferase 1 (DNMT1) is essential for DNA methylation, gene regulation and chromatin stability. We previously discovered DNMT1 mutations cause hereditary sensory and autonomic neuropathy type 1 with dementia and hearing loss (HSAN1E; OMIM 614116). HSAN1E is the first adult-onset neurodegenerative disorder caused by a defect in a methyltransferase gene. HSAN1E patients appear clinically normal until young adulthood, then begin developing the characteristic symptoms involving central and peripheral nervous systems. Some HSAN1E patients also develop narcolepsy and it has recently been suggested that HSAN1E is allelic to autosomal dominant cerebellar ataxia, deafness, with narcolepsy (ADCA-DN; OMIM 604121), which is also caused by mutations in DNMT1. A hotspot mutation Y495C within the targeting sequence domain of DNMT1 has been identified among HSAN1E patients. The mutant DNMT1 protein shows premature degradation and reduced DNA methyltransferase activity. Herein, we investigate genome-wide DNA methylation at single-base resolution through whole-genome bisulfite sequencing of germline DNA in 3 pairs of HSAN1E patients and their gender- and age-matched siblings. Over 1 billion 75-bp single-end reads were generated for each sample. In the 3 affected siblings, overall methylation loss was consistently found in all chromosomes with X and 18 being most affected. Paired sample analysis identified 564,218 differentially methylated CpG sites (DMCs; P < 0.05), of which 300 134 were intergenic and 264 084 genic CpGs. Hypomethylation was predominant in both genic and intergenic regions, including promoters, exons, most CpG islands, L1, L2, Alu, and satellite repeats and simple repeat sequences. In some CpG islands, hypermethylated CpGs outnumbered hypomethylated CpGs. In 201 imprinted genes, there were more DMCs than in non-imprinted genes and most were hypomethylated. Differentially methylated region (DMR) analysis identified 5649 hypomethylated and 1872 hypermethylated regions. Importantly, pathway analysis revealed 1693 genes associated with the identified DMRs were highly associated in diverse neurological disorders and NAD+/NADH metabolism pathways is implicated in the pathogenesis. Our results provide novel insights into the epigenetic mechanism of neurodegeneration arising from a hotspot DNMT1 mutation and reveal pathways potentially important in a broad category of neurological and psychological disorders.     

Comparative genome-wide DNA methylation analysis of colorectal tumor and matched normal tissues

Aberrant DNA methylation often occurs in colorectal cancer (CRC). In our study we applied a genome-wide DNA methylation analysis approach, MethylCap-seq, to map the differentially methylated regions (DMRs) in 24 tumors and matched normal colon samples. In total, 2687 frequently hypermethylated and 468 frequently hypomethylated regions were identified, which include potential biomarkers for CRC diagnosis. Hypermethylation in the tumor samples was enriched at CpG islands and gene promoters, while hypomethylation was distributed throughout the genome. Using epigenetic data from human embryonic stem cells, we show that frequently hypermethylated regions coincide with bivalent loci in human embryonic stem cells. DNA methylation is commonly thought to lead to gene silencing; however, integration of publically available gene expression data indicates that 75% of the frequently hypermethylated genes were most likely already lowly or not expressed in normal tissue. Collectively, our study provides genome-wide DNA methylation maps of CRC, comprehensive lists of DMRs, and gives insights into the role of aberrant DNA methylation in CRC formation.     

Comparative genome-wide DNA methylation analysis of colorectal tumor and matched normal tissues

Aberrant DNA methylation often occurs in colorectal cancer (CRC). In our study we applied a genome-wide DNA methylation analysis approach, MethylCap-seq, to map the differentially methylated regions (DMRs) in 24 tumors and matched normal colon samples. In total, 2687 frequently hypermethylated and 468 frequently hypomethylated regions were identified, which include potential biomarkers for CRC diagnosis. Hypermethylation in the tumor samples was enriched at CpG islands and gene promoters, while hypomethylation was distributed throughout the genome. Using epigenetic data from human embryonic stem cells, we show that frequently hypermethylated regions coincide with bivalent loci in human embryonic stem cells. DNA methylation is commonly thought to lead to gene silencing; however, integration of publically available gene expression data indicates that 75% of the frequently hypermethylated genes were most likely already lowly or not expressed in normal tissue. Collectively, our study provides genome-wide DNA methylation maps of CRC, comprehensive lists of DMRs, and gives insights into the role of aberrant DNA methylation in CRC formation.     

Adaptation of the CHARM DNA methylation platform for the rat genome reveals novel brain region-specific differences

Comprehensive High-throughput Arrays for Relative Methylation (CHARM) was recently developed as an experimental platform and analytic approach to assess DNA methylation (DNAm) at a genome-wide level. Its initial implementation was for human and mouse. We adapted it for rat and sought to examine DNAm differences across tissues and brain regions in this model organism. We extracted DNA from liver, spleen, and three brain regions: cortex, hippocampus, and hypothalamus from adult Sprague Dawley rats. DNA was digested with McrBC, and the resulting methyl-depleted fraction was hybridized to the rat CHARM array along with a mock-treated fraction. Differentially methylated regions (DMRs) between tissue types were detected using normalized methylation log-ratios. In validating 24 of the most significant DMRs by bisulfite pyrosequencing, we detected large mean differences in DNAm, ranging from 33-59%, among the most significant DMRs in the across-tissue comparisons. The comparable figures for the hippocampus vs. hypothalamus DMRs were 14-40%, for the cortex vs. hippocampus DMRs, 12-29%, and for the cortex vs. hypothalamus DMRs, 5-35%, with a correlation of r(2) = 0.92 between the methylation differences in 24 DMRs predicted by CHARM and those validated by bisulfite pyrosequencing. Our adaptation of the CHARM array for the rat genome yielded highly robust results that demonstrate the value of this method in detecting substantial DNAm differences between tissues and across different brain regions. This platform should prove valuable in future studies aimed at examining DNAm differences in particular brain regions of rats exposed to environmental stimuli with potential epigenetic consequences.     

Exploring genome-wide DNA methylation profiles altered in hepatocellular carcinoma using Infinium HumanMethylation 450 BeadChips

Hepatocellular carcinoma (HCC) incidence has increased in the US and also has one of the fastest growing death rates of any cancer. The purpose of the current study was to discover novel genome-wide aberrant DNA methylation patterns in HCC tumors that are predominantly HCV-related. Infinium HumanMethylation 450K BeadChip arrays were used to examine genome-wide DNA methylation profiles in 66 pairs of HCC tumor and adjacent non-tumor tissues. After Bonferroni adjustment, a total of 130,512 CpG sites significantly differed in methylation level in tumor compared with non-tumor tissues, with 28,017 CpG sites hypermethylated and 102,495 hypomethylated in tumor tissues. Absolute tumor/non-tumor methylation differences ≥ 20% were found in 24.9% of the hypermethylated and 43.1% of the hypomethylated CpG sites; almost 10,000 CpG sites have ≥ 30% DNA methylation differences. Most (60.1%) significantly hypermethylated CpG sites are located in CpG islands, with 21.6% in CpG shores and 3.6% in shelves. In contrast, only a small proportion (8.2%) of significantly hypomethylated CpG sites are situated in islands, while most are found in open sea (60.2%), shore (17.3%) or shelf (14.3%) regions. A total of 2,568 significant CpG sites (2,441 hypermethylated and 127 hypomethylated) covering 589 genes are located within 684 differentially methylated regions defined as regions with at least two significant CpG sites displaying > 20% methylation differences in the same direction within 250-bp. The top 500 significant CpG sites can significantly distinguish HCC tumor from adjacent tissues with one misclassification. Within adjacent non-tumor tissues, we also identified 75 CpG sites significantly associated with gender, 228 with HCV infection, 17,207 with cirrhosis, and 56 with both HCV infection and cirrhosis after multiple comparisons adjustment. Aberrant DNA methylation profiles across the genome were identified in tumor tissues from US HCC cases that are predominantly related to HCV infection. These results demonstrate the significance of aberrant DNA methylation in HCC tumorigenesis.     

Adaptation of the CHARM DNA methylation platform for the rat genome reveals novel brain region-specific differences

Comprehensive High-throughput Arrays for Relative Methylation (CHARM) was recently developed as an experimental platform and analytic approach to assess DNA methylation (DNAm) at a genome-wide level. Its initial implementation was for human and mouse. We adapted it for rat and sought to examine DNAm differences across tissues and brain regions in this model organism. We extracted DNA from liver, spleen, and three brain regions: cortex, hippocampus, and hypothalamus from adult Sprague Dawley rats. DNA was digested with McrBC, and the resulting methyl-depleted fraction was hybridized to the rat CHARM array along with a mock-treated fraction. Differentially methylated regions (DMRs) between tissue types were detected using normalized methylation log-ratios. In validating 24 of the most significant DMRs by bisulfite pyrosequencing, we detected large mean differences in DNAm, ranging from 33-59%, among the most significant DMRs in the across-tissue comparisons. The comparable figures for the hippocampus vs. hypothalamus DMRs were 14-40%, for the cortex vs. hippocampus DMRs, 12-29%, and for the cortex vs. hypothalamus DMRs, 5-35%, with a correlation of r² = 0.92 between the methylation differences in 24 DMRs predicted by CHARM and those validated by bisulfite pyrosequencing. Our adaptation of the CHARM array for the rat genome yielded highly robust results that demonstrate the value of this method in detecting substantial DNAm differences between tissues and across different brain regions. This platform should prove valuable in future studies aimed at examining DNAm differences in particular brain regions of rats exposed to environmental stimuli with potential epigenetic consequences.     

Existence and possible roles of independent non-CpG methylation in the mammalian brain

Methylated non-CpGs (mCpHs) in mammalian cells yield weak enrichment signals and colocalize with methylated CpGs (mCpGs), thus have been considered byproducts of hyperactive methyltransferases. However, mCpHs are cell type-specific and associated with epigenetic regulation, although their dependency on mCpGs remains to be elucidated. In this study, we demonstrated that mCpHs colocalize with mCpGs in pluripotent stem cells, but not in brain cells. In addition, profiling genome-wide methylation patterns using a hidden Markov model revealed abundant genomic regions in which CpGs and CpHs are differentially methylated in brain. These regions were frequently located in putative enhancers, and mCpHs within the enhancers increased in correlation with brain age. The enhancers with hypermethylated CpHs were associated with genes functionally enriched in immune responses, and some of the genes were related to neuroinflammation and degeneration. This study provides insight into the roles of non-CpG methylation as an epigenetic code in the mammalian brain genome.     

Genome-Wide Analysis of DNA Methylation before-and after Exercise in the Thoroughbred Horse with MeDIP-Seq

Athletic performance is an important criteria used for the selection of superior horses. However, little is known about exercise-related epigenetic processes in the horse. DNA methylation is a key mechanism for regulating gene expression in response to environmental changes. We carried out comparative genomic analysis of genome-wide DNA methylation profiles in the blood samples of two different thoroughbred horses before and after exercise by methylated-DNA immunoprecipitation sequencing (MeDIP-Seq). Differentially methylated regions (DMRs) in the pre-and post-exercise blood samples of superior and inferior horses were identified. Exercise altered the methylation patterns. After 30 min of exercise, 596 genes were hypomethylated and 715 genes were hypermethylated in the superior horse, whereas in the inferior horse, 868 genes were hypomethylated and 794 genes were hypermethylated. These genes were analyzed based on gene ontology (GO) annotations and the exercise-related pathway patterns in the two horses were compared. After exercise, gene regions related to cell division and adhesion were hypermethylated in the superior horse, whereas regions related to cell signaling and transport were hypermethylated in the inferior horse. Analysis of the distribution of methylated CpG islands confirmed the hypomethylation in the gene-body methylation regions after exercise. The methylation patterns of transposable elements also changed after exercise. Long interspersed nuclear elements (LINEs) showed abundance of DMRs. Collectively, our results serve as a basis to study exercise-based reprogramming of epigenetic traits.     

Epigenetic Basis of Regeneration: Analysis of Genomic DNA Methylation Profiles in the MRL/MpJ Mouse

Epigenetic regulation plays essential role in cell differentiation and dedifferentiation, which are the intrinsic processes involved in regeneration. To investigate the epigenetic basis of regeneration capacity, we choose DNA methylation as one of the most important epigenetic mechanisms and the MRL/MpJ mouse as a model of mammalian regeneration known to exhibit enhanced regeneration response in different organs. We report the comparative analysis of genomic DNA methylation profiles of the MRL/MpJ and the control C57BL/6J mouse. Methylated DNA immunoprecipitation followed by microarray analysis using the Nimblegen ‘3 × 720 K CpG Island Plus RefSeq Promoter’ platform was applied in order to carry out genome-wide DNA methylation profiling covering 20 404 promoter regions. We identified hundreds of hypo- and hypermethylated genes and CpG islands in the heart, liver, and spleen, and 37 of them in the three tissues. Decreased inter-tissue diversification and the shift of DNA methylation balance upstream the genes distinguish the genomic methylation patterns of the MRL/MpJ mouse from the C57BL/6J. Homeobox genes and a number of other genes involved in embryonic morphogenesis are significantly overrepresented among the genes hypomethylated in the MRL/MpJ mouse. These findings indicate that epigenetic patterning might be a likely molecular basis of regeneration capability in the MRL/MpJ mouse.     

Genome-wide DNA methylation scan in major depressive disorder

While genome-wide association studies are ongoing to identify sequence variation influencing susceptibility to major depressive disorder (MDD), epigenetic marks, such as DNA methylation, which can be influenced by environment, might also play a role. Here we present the first genome-wide DNA methylation (DNAm) scan in MDD. We compared 39 postmortem frontal cortex MDD samples to 26 controls. DNA was hybridized to our Comprehensive High-throughput Arrays for Relative Methylation (CHARM) platform, covering 3.5 million CpGs. CHARM identified 224 candidate regions with DNAm differences >10%. These regions are highly enriched for neuronal growth and development genes. Ten of 17 regions for which validation was attempted showed true DNAm differences; the greatest were in PRIMA1, with 12-15% increased DNAm in MDD (p = 0.0002-0.0003), and a concomitant decrease in gene expression. These results must be considered pilot data, however, as we could only test replication in a small number of additional brain samples (n = 16), which showed no significant difference in PRIMA1. Because PRIMA1 anchors acetylcholinesterase in neuronal membranes, decreased expression could result in decreased enzyme function and increased cholinergic transmission, consistent with a role in MDD. We observed decreased immunoreactivity for acetylcholinesterase in MDD brain with increased PRIMA1 DNAm, non-significant at p = 0.08.While we cannot draw firm conclusions about PRIMA1 DNAm in MDD, the involvement of neuronal development genes across the set showing differential methylation suggests a role for epigenetics in the illness. Further studies using limbic system brain regions might shed additional light on this role.     

Principal component analysis- and tensor decomposition-based unsupervised feature extraction to select more suitable differentially methylated cytosines: Optimization of standard deviation versus state-of-the-art methods

In contrast to RNA-seq analysis, which has various standard methods, no standard methods for identifying differentially methylated cytosines (DMCs) exist. To identify DMCs, we tested principal component analysis and tensor decomposition-based unsupervised feature extraction with optimized standard deviation, which has been shown to be effective for differentially expressed gene (DEG) identification. The proposed method outperformed certain conventional methods, including those that assume beta-binomial distribution for methylation as the proposed method does not require this, especially when applied to methylation profiles measured using high throughput sequencing. DMCs identified by the proposed method also significantly overlapped with various functional sites, including known differentially methylated regions, enhancers, and DNase I hypersensitive sites. The proposed method was applied to data sets retrieved from The Cancer Genome Atlas to identify DMCs using American Joint Committee on Cancer staging system edition labels. This suggests that the proposed method is a promising standard method for identifying DMCs.     

Genome-scale analysis of DNA methylation in colorectal cancer using Infinium HumanMethylation450 BeadChips

Illumina’s Infinium HumanMethylation450 BeadChip arrays were used to examine genome-wide DNA methylation profiles in 22 sample pairs from colorectal cancer (CRC) and adjacent tissues and 19 colon tissue samples from cancer-free donors. We show that the methylation profiles of tumors and healthy tissue samples can be clearly distinguished from one another and that the main source of methylation variability is associated with disease status. We used different statistical approaches to evaluate the methylation data. In general, at the CpG-site level, we found that common CRC-specific methylation patterns consist of at least 15,667 CpG sites that were significantly different from either adjacent healthy tissue or tissue from cancer-free subjects. Of these sites, 10,342 were hypermethylated in CRC, and 5,325 were hypomethylated. Hypermethylated sites were common in the maximum number of sample pairs and were mostly located in CpG islands, where they were significantly enriched for differentially methylated regions known to be cancer-specific. In contrast, hypomethylated sites were mostly located in CpG shores and were generally sample-specific. Despite the considerable variability in methylation data, we selected a panel of 14 highly robust candidates showing methylation marks in genes SND1, ADHFE1, OPLAH, TLX2, C1orf70, ZFP64, NR5A2, and COL4A. This set was successfully cross-validated using methylation data from 209 CRC samples and 38 healthy tissue samples from The Cancer Genome Atlas consortium (AUC = 0.981 [95% CI: 0.9677–0.9939], sensitivity = 100% and specificity = 82%). In summary, this study reports a large number of loci with novel differential methylation statuses, some of which may serve as candidate markers for diagnostic purposes.     

Short‐term calorie restriction ameliorates genomewide,age‐related alterations in DNA methylation

DNA methylation plays major roles in many biological processes, including aging, carcinogenesis, and development. Analyses of DNA methylation using next‐generation sequencing offer a new way to profile and compare methylomes across the genome in the context of aging. We explored genomewide DNA methylation and the effects of short‐term calorie restriction (CR) on the methylome of aged rat kidney. Whole‐genome methylation of kidney in young (6 months old), old (25 months old), and OCR (old with 4‐week, short‐term CR) rats was analyzed by methylated DNA immunoprecipitation and next‐generation sequencing (MeDIP‐Seq). CpG islands and repetitive regions were hypomethylated, but 5′‐UTR, exon, and 3′‐UTR hypermethylated in old and OCR rats. The methylation in the promoter and intron regions was decreased in old rats, but increased in OCR rats. Pathway enrichment analysis showed that the hypermethylated promoters in old rats were associated with degenerative phenotypes such as cancer and diabetes. The hypomethylated promoters in old rats related significantly to the chemokine signaling pathway. However, the pathways significantly enriched in old rats were not observed from the differentially methylated promoters in OCR rats. Thus, these findings suggest that short‐term CR could partially ameliorate age‐related methylation changes in promoters in old rats. From the epigenomic data, we propose that the hypermethylation found in the promoter regions of disease‐related genes during aging may indicate increases in susceptibility to age‐related diseases. Therefore, the CR‐induced epigenetic changes that ameliorate age‐dependent aberrant methylation may be important to CR's health‐ and life‐prolonging effects.     

Genome-wide DNA methylation analysis in precursor B-cells

DNA methylation is responsible for regulating gene expression and cellular differentiation and for maintaining genomic stability during normal human development. Furthermore, it plays a significant role in the regulation of hematopoiesis. In order to elucidate the influence of DNA methylation during B-cell development, genome-wide DNA methylation status of pro-B, pre-BI, pre-BII, and naïve-B-cells isolated from human umbilical cord blood was determined using the methylated CpG island recovery assay followed by next generation sequencing. On average, 182–200 million sequences were generated for each precursor B-cell subset in 10 biological replicates. An overall decrease in methylation was observed during the transition from pro-B to pre-BI, whereas no differential methylation was observed in the pre-BI to pre-BII transition or in the pre-BII to naïve B-cell transition. Most of the methylated regions were located within intergenic and intronic regions not present in a CpG island context. Putative novel enhancers were identified in these regions that were differentially methylated between pro-B and pre-BI cells. The genome-wide methylation profiles are publically available and may be used to gain a better understanding of the involvement of atypical DNA methylation in the pathogenesis of malignancies associated with precursor B-cells.