Mapping of Variable DNA Methylation Across Multiple Cell Types Defines a Dynamic Regulatory Landscape of the Human Genome

Active DNA demethylation is upstream of rod-photoreceptor fate determination and required for retinal development

Retinal cell fate specification from multipotent retinal progenitors is governed by dynamic changes in chromatin structure and gene expression. Methylation at cytosines in DNA (5mC) is actively regulated for proper control of gene expression and chromatin architecture. Numerous genes display active DNA demethylation across retinal development; a process that requires oxidation of 5mC to 5-hydroxymethylcytosine (5hmC) and is controlled by the ten-eleven translocation methylcytosine dioxygenase (TET) enzymes. Using an allelic series of conditional TET enzyme mutants, we determine that DNA demethylation is required upstream of NRL and NR2E3 expression for the establishment of rod-photoreceptor fate. Using histological, behavioral, transcriptomic, and base-pair resolution DNA methylation analyses, we establish that inhibition of active DNA demethylation results in global changes in gene expression and methylation patterns that prevent photoreceptor precursors from adopting a rod-photoreceptor fate, instead producing a retina in which all photoreceptors specify as cones. Our results establish the TET enzymes and DNA demethylation as critical regulators of retinal development and cell fate specification, elucidating a novel mechanism required for the specification of rod-photoreceptors.

Cross-tissue patterns of DNA hypomethylation reveal genetically distinct histories of cell development

BACKGROUND

Establishment of DNA methylation (DNAme) patterns is essential for balanced multi-lineage cellular differentiation, but exactly how these patterns drive cellular phenotypes is unclear. While > 80% of CpG sites are stably methylated, tens of thousands of discrete CpG loci form hypomethylated regions (HMRs). Because they lack DNAme, HMRs are considered transcriptionally permissive, but not all HMRs actively regulate genes. Unlike promoter HMRs, a subset of non-coding HMRs is cell type-specific and enriched for tissue-specific gene regulatory functions. Our data further argues not only that HMR establishment is an important step in enforcing cell identity, but also that cross-cell type and spatial HMR patterns are functionally informative of gene regulation.

RESULTS

To understand the significance of non-coding HMRs, we systematically dissected HMR patterns across diverse human cell types and developmental timepoints, including embryonic, fetal, and adult tissues. Unsupervised clustering of 126,104 distinct HMRs revealed that levels of HMR specificity reflects a developmental hierarchy supported by enrichment of stage-specific transcription factors and gene ontologies. Using a pseudo-time course of development from embryonic stem cells to adult stem and mature hematopoietic cells, we find that most HMRs observed in differentiated cells (~ 60%) are established at early developmental stages and accumulate as development progresses. HMRs that arise during differentiation frequently (~ 35%) establish near existing HMRs (≤ 6 kb away), leading to the formation of HMR clusters associated with stronger enhancer activity. Using SNP-based partitioned heritability from GWAS summary statistics across diverse traits and clinical lab values, we discovered that genetic contribution to trait heritability is enriched within HMRs. Moreover, the contribution of heritability to cell-relevant traits increases with both increasing HMR specificity and HMR clustering, supporting the role of distinct HMR subsets in regulating normal cell function.

CONCLUSIONS

Our results demonstrate that the entire HMR repertoire within a cell-type, rather than just the cell type-specific HMRs, stores information that is key to understanding and predicting cellular phenotypes. Ultimately, these data provide novel insights into how DNA hypo-methylation provides genetically distinct historical records of a cell's journey through development, highlighting HMRs as functionally distinct from other epigenomic annotations.

Profiling the immune epigenome across global cattle breeds

BACKGROUND

Understanding the variation between well and poorly adapted cattle breeds to local environments and pathogens is essential for breeding cattle with improved climate and disease-resistant phenotypes. Although considerable progress has been made towards identifying genetic differences between breeds, variation at the epigenetic and chromatin levels remains poorly characterized. Here, we generate, sequence and analyse over 150 libraries at base-pair resolution to explore the dynamics of DNA methylation and chromatin accessibility of the bovine immune system across three distinct cattle lineages.

RESULTS

We find extensive epigenetic divergence between the taurine and indicine cattle breeds across immune cell types, which is linked to the levels of local DNA sequence divergence between the two cattle sub-species. The unique cell type profiles enable the deconvolution of complex cellular mixtures using digital cytometry approaches. Finally, we show distinct sub-categories of CpG islands based on their chromatin and methylation profiles that discriminate between classes of distal and gene proximal islands linked to discrete transcriptional states.

CONCLUSIONS

Our study provides a comprehensive resource of DNA methylation, chromatin accessibility and RNA expression profiles of three diverse cattle populations. The findings have important implications, from understanding how genetic editing across breeds, and consequently regulatory backgrounds, may have distinct impacts to designing effective cattle epigenome-wide association studies in non-European breeds.

Epigenomic signatures reveal mechanistic clues and predictive markers for autism spectrum disorder

Autism spectrum disorder (ASD) comprises a heterogeneous group of neurodevelopmental outcomes in children with a commonality in deficits in social communication and language combined with repetitive behaviors and interests. The etiology of ASD is heterogeneous, as several hundred genes have been implicated as well as multiple in utero environmental exposures. Over the past two decades, epigenetic investigations, including DNA methylation, have emerged as a novel way to capture the complex interface of multivariate ASD etiologies. More recently, epigenome-wide association studies using human brain and surrogate accessible tissues have revealed some convergent genes that are epigenetically altered in ASD, many of which overlap with known genetic risk factors. Unlike transcriptomes, epigenomic signatures defined by DNA methylation from surrogate tissues such as placenta and cord blood can reflect past differences in fetal brain gene transcription, transcription factor binding, and chromatin. For example, the discovery of NHIP (neuronal hypoxia inducible, placenta associated) through an epigenome-wide association in placenta, identified a common genetic risk for ASD that was modified by prenatal vitamin use. While epigenomic signatures are distinct between different genetic syndromic causes of ASD, bivalent chromatin and some convergent gene pathways are consistently epigenetically altered in both syndromic and idiopathic ASD, as well as some environmental exposures. Together, these epigenomic signatures hold promising clues towards improved early prediction and prevention of ASD as well genes and gene pathways to target for pharmacological interventions. Future advancements in single cell and multi-omic technologies, machine learning, as well as non-invasive screening of epigenomic signatures during pregnancy or newborn periods are expected to continue to impact the translatability of the recent discoveries in epigenomics to precision public health.

epialleleR: an R/Bioconductor package for sensitive allele-specific methylation analysis in NGS data

Low-level mosaic epimutations within the BRCA1 gene promoter occur in 5-8% of healthy individuals and are associated with a significantly elevated risk of breast and ovarian cancer. Similar events may also affect other tumor suppressor genes, potentially being a significant contributor to cancer burden. While this opens a new area for translational research, detection of low-level mosaic epigenetic events requires highly sensitive and robust methodology for methylation analysis. We here present epialleleR, a computational framework for sensitive detection, quantification, and visualization of mosaic epimutations in methylation sequencing data. Analyzing simulated and real data sets, we provide in-depth assessments of epialleleR performance and show that linkage to epihaplotype data is necessary to detect low-level methylation events. The epialleleR is freely available at https://github.com/BBCG/epialleleR and https://bioconductor.org/packages/epialleleR/ as an open-source R/Bioconductor package.

DNMT3B overexpression downregulates genes with CpG islands, common motifs, and transcription factor binding sites that interact with DNMT3B

DNA methylation is a key epigenetic modification to regulate gene expression in mammalian cells. Abnormal DNA methylation in gene promoters is common across human cancer types. DNMT3B is the main de novo methyltransferase enhanced in several primary tumors. How de novo methylation is established in genes related to cancer is poorly understood. CpG islands (CGIs), common sequences, and transcription factors (TFs) that interact with DNMT3B have been associated with abnormal de novo methylation. We initially identified cis elements associated with DNA methylation to investigate the contribution of DNMT3B overexpression to the deregulation of its possible target genes in an epithelial cell model. In a set of downregulated genes (n = 146) from HaCaT cells with DNMT3B overexpression, we found CGI, common sequences, and TFs Binding Sites that interact with DNMT3B (we called them P-down-3B). PPL1, VAV3, IRF1, and BRAF are P-down-3B genes that are downregulated and increased their methylation in DNMT3B presence. Together these findings suggest that methylated promoters aberrantly have some cis elements that could conduce de novo methylation by DNMT3B.

Genetic and in utero environmental contributions to DNA methylation variation in placenta

Genetic and prenatal environmental factors shape fetal development and cardiometabolic health in later life. A key target of genetic and prenatal environmental factors is the epigenome of the placenta, an organ that is implicated in fetal growth and diseases in later life. This study had two aims: (1) to identify and functionally characterize placental variably methylated regions (VMRs), which are regions in the epigenome with high inter-individual methylation variability; and (2) to investigate the contributions of fetal genetic loci and 12 prenatal environmental factors (maternal cardiometabolic-,psychosocial-, demographic- and obstetric-related) on methylation at each VMR. Akaike's information criterion was used to select the best model out of four models [prenatal environment only, genotype only, additive effect of genotype and prenatal environment (G + E), and their interaction effect (G × E)]. We identified 5850 VMRs in placenta. Methylation at 70% of VMRs was best explained by G × E, followed by genotype only (17.7%), and G + E (12.3%). Prenatal environment alone best explained only 0.03% of VMRs. We observed that 95.4% of G × E models and 93.9% of G + E models included maternal age, parity, delivery mode, maternal depression or gestational weight gain. VMR methylation sites and their regulatory genetic variants were enriched (P < 0.05) for genomic regions that have known links with regulatory functions and complex traits. This study provided a genome-wide catalog of VMRs in placenta and highlighted that variation in placental DNA methylation at loci with regulatory and trait relevance is best elucidated by integrating genetic and prenatal environmental factors, and rarely by environmental factors alone.

Disease-related blood-based differential methylation in cystic fibrosis and its representation in lung cancer revealed a regulatory locus in PKP3 in lung epithelial cells

Cystic fibrosis (CF) is a monogenic disease, characterized by massive chronic lung inflammation. The observed variability in clinical phenotypes in monozygotic CF twins is likely associated with the extent of inflammation. This study sought to investigate inflammation-related aberrant DNA methylation in CF twins and to determine to what extent acquired methylation changes may be associated with lung cancer.

Blood-based genome-wide DNA methylation analysis was performed to compare the DNA methylomes of monozygotic twins, from the European CF Twin and Sibling Study with various degrees of disease severity. Putatively inflammation-related and differentially methylated positions were selected from a large lung cancer case-control study and investigated in blood by targeted bisulphite next-generation-sequencing. An inflammation-related locus located in the Plakophilin-3 (PKP3) gene was functionally analysed regarding promoter and enhancer activity in presence and absence of methylation using luciferase reporter assays.

We confirmed in a unique cohort that monozygotic twins, even if clinically discordant, have only minor differences in global DNA methylation patterns and blood cell composition. Further, we determined the most differentially methylated positions, a high proportion of which are blood cell-type-specific, whereas others may be acquired and thus have potential relevance in the context of inflammation as lung cancer risk factors. We identified a sequence in the gene body of PKP3 which is hypermethylated in blood from CF twins with severe phenotype and highly variably methylated in lung cancer patients and controls, independent of known clinical parameters, and showed that this region exhibits methylation-dependent promoter activity in lung epithelial cells.

i4mC-Deep: An Intelligent Predictor of N4-Methylcytosine Sites Using a Deep Learning Approach with Chemical Properties

DNA is subject to epigenetic modification by the molecule N4-methylcytosine (4mC). N4-methylcytosine plays a crucial role in DNA repair and replication, protects host DNA from degradation, and regulates DNA expression. However, though current experimental techniques can identify 4mC sites, such techniques are expensive and laborious. Therefore, computational tools that can predict 4mC sites would be very useful for understanding the biological mechanism of this vital type of DNA modification. Conventional machine-learning-based methods rely on hand-crafted features, but the new method saves time and computational cost by making use of learned features instead. In this study, we propose i4mC-Deep, an intelligent predictor based on a convolutional neural network (CNN) that predicts 4mC modification sites in DNA samples. The CNN is capable of automatically extracting important features from input samples during training. Nucleotide chemical properties and nucleotide density, which together represent a DNA sequence, act as CNN input data. The outcome of the proposed method outperforms several state-of-the-art predictors. When i4mC-Deep was used to analyze G. subterruneus DNA, the accuracy of the results was improved by 3.9% and MCC increased by 10.5% compared to a conventional predictor.

DNA methylation landscape of 16 canine somatic tissues by methylation-sensitive restriction enzyme-based next generation sequencing

DNA methylation plays important functions in gene expression regulation that is involved in individual development and various diseases. DNA methylation has been well studied in human and model organisms, but only limited data exist in companion animals like dog. Using methylation-sensitive restriction enzyme-based next generation sequencing (Canine DREAM), we obtained canine DNA methylation maps of 16 somatic tissues from two dogs. In total, we evaluated 130,861 CpG sites. The majority of CpG sites were either highly methylated (> 70%, 52.5-64.6% of all CpG sites analyzed) or unmethylated (< 30%, 22.5-28.0% of all CpG sites analyzed) which are methylation patterns similar to other species. The overall methylation status of CpG sites across the 32 methylomes were remarkably similar. However, the tissue types were clearly defined by principle component analysis and hierarchical clustering analysis with DNA methylome. We found 6416 CpG sites located closely at promoter region of genes and inverse correlation between DNA methylation and gene expression of these genes. Our study provides basic dataset for DNA methylation profiles in dogs.

Combined effects of genotype and childhood adversity shape variability of DNA methylation across age

Lasting effects of adversity, such as exposure to childhood adversity (CA) on disease risk, may be embedded via epigenetic mechanisms but findings from human studies investigating the main effects of such exposure on epigenetic measures, including DNA methylation (DNAm), are inconsistent. Studies in perinatal tissues indicate that variability of DNAm at birth is best explained by the joint effects of genotype and prenatal environment. Here, we extend these analyses to postnatal stressors. We investigated the contribution of CA, cis genotype (G), and their additive (G + CA) and interactive (G × CA) effects to DNAm variability in blood or saliva from five independent cohorts with a total sample size of 1074 ranging in age from childhood to late adulthood. Of these, 541 were exposed to CA, which was assessed retrospectively using self-reports or verified through social services and registries. For the majority of sites (over 50%) in the adult cohorts, variability in DNAm was best explained by G + CA or G × CA but almost never by CA alone. Across ages and tissues, 1672 DNAm sites showed consistency of the best model in all five cohorts, with G × CA interactions explaining most variance. The consistent G × CA sites mapped to genes enriched in brain-specific transcripts and Gene Ontology terms related to development and synaptic function. Interaction of CA with genotypes showed the strongest contribution to DNAm variability, with stable effects across cohorts in functionally relevant genes. This underscores the importance of including genotype in studies investigating the impact of environmental factors on epigenetic marks.

Maternal lipodome across pregnancy is associated with the neonatal DNA methylome

Aim: To classify the association between the maternal lipidome and DNA methylation in cord blood leukocytes. Materials & methods: Untargeted lipidomics was performed on first trimester maternal plasma (M1) and delivery maternal plasma (M3) in 100 mothers from the Michigan Mother-Infant Pairs cohort. Cord blood leukocyte DNA methylation was profiled using the Infinium EPIC bead array and empirical Bayes modeling identified differential DNA methylation related to maternal lipid groups. Results: M3-saturated lysophosphatidylcholine was associated with 45 differentially methylated loci and M3-saturated lysophosphatidylethanolamine was associated with 18 differentially methylated loci. Biological pathways enriched among differentially methylated loci by M3 saturated lysophosphatidylcholines were related to cell proliferation and growth. Conclusion: The maternal lipidome may be influential in establishing the infant epigenome.

DNA Methylation in Nonalcoholic Fatty Liver Disease

Nonalcoholic fatty liver disease (NAFLD) is a widespread hepatic disorder in the United States and other Westernized countries. Nonalcoholic steatohepatitis (NASH), an advanced stage of NAFLD, can progress to end-stage liver disease, including cirrhosis and liver cancer. Poor understanding of mechanisms underlying NAFLD progression from simple steatosis to NASH has limited the development of effective therapies and biomarkers. An accumulating body of studies has suggested the importance of DNA methylation, which plays pivotal roles in NAFLD pathogenesis. DNA methylation signatures that can affect gene expression are influenced by environmental and lifestyle experiences such as diet, obesity, and physical activity and are reversible. Hence, DNA methylation signatures and modifiers in NAFLD may provide the basis for developing biomarkers indicating the onset and progression of NAFLD and therapeutics for NAFLD. Herein, we review an update on the recent findings in DNA methylation signatures and their roles in the pathogenesis of NAFLD and broaden people's perspectives on potential DNA methylation-related treatments and biomarkers for NAFLD.

Peritoneal fibrosis and epigenetic modulation

Peritoneal dialysis (PD) is an effective treatment for patients with end-stage renal disease. However, peritoneal fibrosis (PF) is a common complication that ultimately leads to ultrafiltration failure and discontinuation of PD after long-term PD therapy. There is currently no effective therapy to prevent or delay this pathologic process. Recent studies have reported epigenetic modifications involved in PF, and accumulating evidence suggests that epigenetic therapies may have the potential to prevent and treat PF clinically. The major epigenetic modifications in PF include DNA methylation, histone modification, and noncoding RNAs. The mechanisms of epigenetic regulation in PF are complex, predominantly involving modification of signaling molecules, transcriptional factors, and genes. This review will describe the mechanisms of epigenetic modulation in PF and discuss the possibility of targeting them to prevent and treat this complication.

Glucose-induced, duration-dependent genome-wide DNA methylation changes in human endothelial cells

DNA methylation, a critical epigenetic mechanism, plays an important role in governing gene expressions during biological processes such as aging, which is well known to be accelerated in hyperglycemia (diabetes). In the present study, we investigated the effects of glucose on whole genome DNA methylation in small [human retinal microvascular endothelial cells (HRECs)] and large [human umbilical vein endothelial cells (HUVECs)] vessel endothelial cell (EC) lines exposed to basal or high glucose-containing media for variable lengths of time. Using the Infinium EPIC array, we obtained 773,133 CpG sites (probes) for analysis. Unsupervised clustering of the top 5% probes identified four distinct clusters within EC groups, with significant methylation differences attributed to EC types and the duration of cell culture rather than glucose stimuli alone. When comparing the ECs incubated for 2 days versus 7 days, hierarchical clustering analyses [methylation change >10% and false discovery rate (FDR) <0.05] identified 17,354 and 128 differentially methylated CpGs for HUVECs and HRECs, respectively. Predominant DNA hypermethylation was associated with the length of culture and was enriched for gene enhancer elements and regions surrounding CpG shores and shelves. We identified 88 differentially methylated regions (DMRs) for HUVECs and 8 DMRs for HRECs (all FDR <0.05). Pathway enrichment analyses of DMRs highlighted involvement of regulators of embryonic development (i.e., HOX genes) and cellular differentiation [transforming growth factor-β (TGF-β) family members]. Collectively, our findings suggest that DNA methylation is a complex process that involves tightly coordinated, cell-specific mechanisms. Such changes in methylation overlap genes critical for cellular differentiation and embryonic development.

Identification of DNA motifs that regulate DNA methylation

DNA methylation is an important epigenetic mark but how its locus-specificity is decided in relation to DNA sequence is not fully understood. Here, we have analyzed 34 diverse whole-genome bisulfite sequencing datasets in human and identified 313 motifs, including 92 and 221 associated with methylation (methylation motifs, MMs) and unmethylation (unmethylation motifs, UMs), respectively. The functionality of these motifs is supported by multiple lines of evidence. First, the methylation levels at the MM and UM motifs are respectively higher and lower than the genomic background. Second, these motifs are enriched at the binding sites of methylation modifying enzymes including DNMT3A and TET1, indicating their possible roles of recruiting these enzymes. Third, these motifs significantly overlap with "somatic QTLs" (quantitative trait loci) of methylation and expression. Fourth, disruption of these motifs by mutation is associated with significantly altered methylation level of the CpGs in the neighbor regions. Furthermore, these motifs together with somatic mutations are predictive of cancer subtypes and patient survival. We revealed some of these motifs were also associated with histone modifications, suggesting a possible interplay between the two types of epigenetic modifications. We also found some motifs form feed forward loops to contribute to DNA methylation dynamics.

Integrated analysis of environmental and genetic influences on cord blood DNA methylation in new-borns

Epigenetic processes, including DNA methylation (DNAm), are among the mechanisms allowing integration of genetic and environmental factors to shape cellular function. While many studies have investigated either environmental or genetic contributions to DNAm, few have assessed their integrated effects. Here we examine the relative contributions of prenatal environmental factors and genotype on DNA methylation in neonatal blood at variably methylated regions (VMRs) in 4 independent cohorts (overall n = 2365). We use Akaike’s information criterion to test which factors best explain variability of methylation in the cohort-specific VMRs: several prenatal environmental factors (E), genotypes in cis (G), or their additive (G + E) or interaction (GxE) effects. Genetic and environmental factors in combination best explain DNAm at the majority of VMRs. The CpGs best explained by either G, G + E or GxE are functionally distinct. The enrichment of genetic variants from GxE models in GWAS for complex disorders supports their importance for disease risk.

Environmental influences during prenatal development may have implications for health and disease later in life. Here, Czamara et al. assess DNA methylation in cord blood from new-born under various models including environmental and genetic effects individually and their additive or interaction effects.

Exploring the genetic basis of human population differences in DNA methylation and their causal impact on immune gene regulation

BACKGROUND

DNA methylation is influenced by both environmental and genetic factors and is increasingly thought to affect variation in complex traits and diseases. Yet, the extent of ancestry-related differences in DNA methylation, their genetic determinants, and their respective causal impact on immune gene regulation remain elusive.

RESULTS

We report extensive population differences in DNA methylation between 156 individuals of African and European descent, detected in primary monocytes that are used as a model of a major innate immunity cell type. Most of these differences (~ 70%) are driven by DNA sequence variants nearby CpG sites, which account for ~ 60% of the variance in DNA methylation. We also identify several master regulators of DNA methylation variation in trans, including a regulatory hub nearby the transcription factor-encoding CTCF gene, which contributes markedly to ancestry-related differences in DNA methylation. Furthermore, we establish that variation in DNA methylation is associated with varying gene expression levels following mostly, but not exclusively, a canonical model of negative associations, particularly in enhancer regions. Specifically, we find that DNA methylation highly correlates with transcriptional activity of 811 and 230 genes, at the basal state and upon immune stimulation, respectively. Finally, using a Bayesian approach, we estimate causal mediation effects of DNA methylation on gene expression in ~ 20% of the studied cases, indicating that DNA methylation can play an active role in immune gene regulation.

CONCLUSION

Using a system-level approach, our study reveals substantial ancestry-related differences in DNA methylation and provides evidence for their causal impact on immune gene regulation.

Enhancer DNA methylation in acute myeloid leukemia and myelodysplastic syndromes

DNA methylation (CpG methylation) exerts an important role in normal differentiation and proliferation of hematopoietic stem cells and their differentiated progeny, while it has also the ability to regulate myeloid versus lymphoid fate. Mutations of the epigenetic machinery are observed in hematological malignancies including acute myeloid leukemia (AML) and myelodysplastic syndromes (MDS) resulting in hyper- or hypo-methylation affecting several different pathways. Enhancers are cis-regulatory elements which promote transcription activation and are characterized by histone marks including H3K27ac and H3K4me1/2. These gene subunits are target gene expression 'fine-tuners', are differentially used during the hematopoietic differentiation, and, in contrast to promoters, are not shared by the different hematopoietic cell types. Although the interaction between gene promoters and DNA methylation has extensively been studied, much less is known about the interplay between enhancers and DNA methylation. In hematopoiesis, DNA methylation at enhancers has the potential to discriminate between fetal and adult erythropoiesis, and also is a regulatory mechanism in granulopoiesis through repression of neutrophil-specific enhancers in progenitor cells during maturation. The interplay between DNA methylation at enhancers is disrupted in AML and MDS and mainly hyper-methylation at enhancers raising early during myeloid lineage commitment is acquired during malignant transformation. Interactions between mutated epigenetic drivers and other oncogenic mutations also affect enhancers' activity with final result, myeloid differentiation block. In this review, we have assembled recent data regarding DNA methylation and enhancers' activity in normal and mainly myeloid malignancies.

Increased methylation at an unexplored glucocorticoid responsive element within exon 1D of NR3C1 gene is related to anxious-depressive disorders and decreased hippocampal connectivity

Among the major psychiatric disorders, anxious-depressive disorders stand out as one of the more prevalent and more frequently associated with hypothalamic-pituitary-adrenal (HPA) axis abnormalities. Methylation at the exon 1_F of the glucocorticoid receptor gene NR3C1 has been associated with both early stress exposure and risk for developing a psychiatric disorder; however, other NR3C1 promoter regions have been underexplored. Exon 1_D emerges as a suggestive new target in stress-related disorders epigenetically sensitive to early adversity. After assessment of 48 monozygotic twin pairs (n=96 subjects) informative for lifetime history of anxious-depressive disorders, they were classified as concordant, discordant or healthy in function of whether both, one or neither twin in each pair had a lifetime diagnosis of anxious-depressive disorders. DNA for epigenetic analysis was extracted from peripheral blood. Exon 1_F and exon 1_D CpG-specific methylation was analysed by means of pyrosequencing technology. Functional magnetic resonance imaging was available for 54 subjects (n=27 twin pairs). Exon 1_D CpG-specific methylation within a glucocorticoid responsive element (GRE) was correlated with familial burden of anxious-depressive disorders (r=0.35, z=2.26, p=0.02). Right hippocampal connectivity was significantly associated with CpG-specific GRE methylation (β=-2.33, t=-2.85, p=0.01). Exon 1_F was uniformly hypomethylated across all subgroups of the present sample. GRE hypermethylation at exon 1_D of the NR3C1 gene in monozygotic twins concordant for anxious-depressive disorders suggests this region plays a role in increasing vulnerability to psychosocial stress, partly mediated by altered hippocampal connectivity.

Obligatory and facilitative allelic variation in the DNA methylome within common disease-associated loci

Integrating epigenetic data with genome-wide association study (GWAS) results can reveal disease mechanisms. The genome sequence itself also shapes the epigenome, with CpG density and transcription factor binding sites (TFBSs) strongly encoding the DNA methylome. Therefore, genetic polymorphism impacts on the observed epigenome. Furthermore, large genetic variants alter epigenetic signal dosage. Here, we identify DNA methylation variability between GWAS-SNP risk and non-risk haplotypes. In three subsets comprising 3128 MeDIP-seq peripheral-blood DNA methylomes, we find 7173 consistent and functionally enriched Differentially Methylated Regions. 36.8% can be attributed to common non-SNP genetic variants. CpG-SNPs, as well as facilitative TFBS-motifs, are also enriched. Highlighting their functional potential, CpG-SNPs strongly associate with allele-specific DNase-I hypersensitivity sites. Our results demonstrate strong DNA methylation allelic differences driven by obligatory or facilitative genetic effects, with potential direct or regional disease-related repercussions. These allelic variations require disentangling from pure tissue-specific modifications, may influence array studies, and imply underestimated population variability in current reference epigenomes.

Genomic polymorphisms affect the epigenome, which in turn influences how epigenome- and genome-wide analysis are interpreted. Here, the authors characterise allelic differences in DNA methylation driven by obligatory or facilitative genetic effects, which may affect disease-related loci.

Cerebral Organoids Recapitulate Epigenomic Signatures of the Human Fetal Brain

Organoids derived from human pluripotent stem cells recapitulate the early three-dimensional organization of the human brain, but whether they establish the epigenomic and transcriptional programs essential for brain development is unknown. We compared epigenomic and regulatory features in cerebral organoids and human fetal brain, using genome-wide, base resolution DNA methylome and transcriptome sequencing. Transcriptomic dynamics in organoids faithfully modeled gene expression trajectories in early-to-mid human fetal brains. We found that early non-CG methylation accumulation at super-enhancers in both fetal brain and organoids marks forthcoming transcriptional repression in the fully developed brain. Demethylated regions (74% of 35,627) identified during organoid differentiation overlapped with fetal brain regulatory elements. Interestingly, pericentromeric repeats showed widespread demethylation in multiple types of in vitro human neural differentiation models but not in fetal brain. Our study reveals that organoids recapitulate many epigenomic features of mid-fetal human brain and also identified novel non-CG methylation signatures of brain development.

Defining baseline epigenetic landscapes in the rat liver

Aim

Characterization of the hepatic epigenome following exposure to chemicals and therapeutic drugs provides novel insights into toxicological and pharmacological mechanisms, however appreciation of genome-wide inter- and intra-strain baseline epigenetic variation, particularly in under-characterized species such as the rat is limited.

Material & methods

To enhance the utility of epigenomic endpoints safety assessment, we map both DNA modifications (5-methyl-cytosine and 5-hydroxymethyl-cytosine) and enhancer related chromatin marks (H3K4me1 and H3K27ac) across multiple male and female rat livers for two important outbred laboratory rat strains (Sprague–Dawley and Wistar).

Results & conclusion

Integration of DNA modification, enhancer chromatin marks and gene expression profiles reveals clear gender-specific chromatin states at genes which exhibit gender-specific transcription. Taken together this work provides a valuable baseline liver epigenome resource for rat strains that are commonly used in chemical and pharmaceutical safety assessment.

Integration of DNA methylation and gene transcription across nineteen cell types reveals cell type-specific and genomic region-dependent regulatory patterns

Despite numerous studies done on understanding the role of DNA methylation, limited work has focused on systems integration of cell type-specific interplay between DNA methylation and gene transcription. Through a genome-wide analysis of DNA methylation across 19 cell types with T-47D as reference, we identified 106,252 cell type-specific differentially-methylated CpGs categorized into 7,537 differentially (46.6% hyper- and 53.4% hypo-) methylated regions. We found 44% promoter regions and 75% CpG islands were T-47D cell type-specific methylated. Pyrosequencing experiments validated the cell type-specific methylation across three benchmark cell lines. Interestingly, these DMRs overlapped with 1,145 known tumor suppressor genes. We then developed a Bayesian Gaussian Regression model to measure the relationship among DNA methylation, genomic segment distribution, differential gene expression and tumor suppressor gene status. The model uncovered that 3′UTR methylation has much less impact on transcriptional activity than other regions. Integration of DNA methylation and 82 transcription factor binding information across the 19 cell types suggested diverse interplay patterns between the two regulators. Our integrative analysis reveals cell type-specific and genomic region-dependent regulatory patterns and provides a perspective for integrating hundreds of various omics-seq data together.

DNA methyltransferase 3b regulates articular cartilage homeostasis by altering metabolism

Osteoarthritis (OA) is the most common form of arthritis worldwide. It is a complex disease affecting the whole joint but is generally characterized by progressive degradation of articular cartilage. Recent genome-wide association screens have implicated distinct DNA methylation signatures in OA patients. We show that the de novo DNA methyltransferase (Dnmt) 3b, but not Dnmt3a, is present in healthy murine and human articular chondrocytes and its expression decreases in OA mouse models and in chondrocytes from human OA patients. Targeted deletion of Dnmt3b in murine articular chondrocytes results in an early-onset and progressive postnatal OA-like pathology. RNA-Seq and methylC-Seq analyses of Dnmt3b loss-of-function chondrocytes show that cellular metabolic processes are affected. Specifically, TCA metabolites and mitochondrial respiration are elevated. Importantly, a chondroprotective effect was found following Dnmt3b gain of function in murine articular chondrocytes in vitro and in vivo. This study shows that Dnmt3b plays a significant role in regulating postnatal articular cartilage homeostasis. Cellular pathways regulated by Dnmt3b in chondrocytes may provide novel targets for therapeutic approaches to treat OA.

Epigenetics and Liver Fibrosis

Liver fibrosis arises because prolonged injury combined with excessive scar deposition within hepatic parenchyma arising from overactive wound healing response mediated by activated myofibroblasts. Fibrosis is the common end point for any type of chronic liver injury including alcoholic liver disease, nonalcoholic fatty liver disease, viral hepatitis, and cholestatic liver diseases. Although genetic influences are important, it is epigenetic mechanisms that have been shown to orchestrate many aspects of fibrogenesis in the liver. New discoveries in the field are leading toward the development of epigenetic biomarkers and targeted therapies. This review considers epigenetic mechanisms as well as recent advances in epigenetic programming in the context of hepatic fibrosis.

NGSmethDB 2017: enhanced methylomes and differential methylation

The 2017 update of NGSmethDB stores whole genome methylomes generated from short-read data sets obtained by bisulfite sequencing (WGBS) technology. To generate high-quality methylomes, stringent quality controls were integrated with third-part software, adding also a two-step mapping process to exploit the advantages of the new genome assembly models. The samples were all profiled under constant parameter settings, thus enabling comparative downstream analyses. Besides a significant increase in the number of samples, NGSmethDB now includes two additional data-types, which are a valuable resource for the discovery of methylation epigenetic biomarkers: (i) differentially methylated single-cytosines; and (ii) methylation segments (i.e. genome regions of homogeneous methylation). The NGSmethDB back-end is now based on MongoDB, a NoSQL hierarchical database using JSON-formatted documents and dynamic schemas, thus accelerating sample comparative analyses. Besides conventional database dumps, track hubs were implemented, which improved database access, visualization in genome browsers and comparative analyses to third-part annotations. In addition, the database can be also accessed through a RESTful API. Lastly, a Python client and a multiplatform virtual machine allow for program-driven access from user desktop. This way, private methylation data can be compared to NGSmethDB without the need to upload them to public servers. Database website: http://bioinfo2.ugr.es/NGSmethDB.

The application of genome-wide 5-hydroxymethylcytosine studies in cancer research

Early detection and characterization of molecular events associated with tumorgenesis remain high priorities. Genome-wide epigenetic assays are promising diagnostic tools, as aberrant epigenetic events are frequent and often cancer specific. The deposition and analysis of multiple patient-derived cancer epigenomic profiles contributes to our appreciation of the underlying biology; aiding the detection of novel identifiers for cancer subtypes. Modifying enzymes and co-factors regulating these epigenetic marks are frequently mutated in cancers, and as epigenetic modifications themselves are reversible, this makes their study very attractive with respect to pharmaceutical intervention. Here we focus on the novel modified base, 5-hydoxymethylcytosine, and discuss how genome-wide 5-hydoxymethylcytosine profiling expedites our molecular understanding of cancer, serves as a lineage tracer, classifies the mode of action of potentially carcinogenic agents and clarifies the roles of potential novel cancer drug targets; thus assisting the development of new diagnostic/prognostic tools.

Comparative Epigenomic Profiling of the DNA Methylome in Mouse and Zebrafish Uncovers High Interspecies Divergence

The DNA methylation landscape is dynamically patterned during development and distinct methylation patterns distinguish healthy from diseased cells. However, whether tissue-specific methylation patterns are conserved across species is not known. We used comparative methylome analysis of base-resolution DNA methylation profiles from the liver and brain of mouse and zebrafish generated by reduced representation bisulfite sequencing to identify the conserved and divergent aspects of the methylome in these commonly used vertebrate model organisms. On average, 24% of CpGs are methylated in mouse livers and the pattern of methylation was highly concordant among four male mice from two different strains. The same level of methylation (24.2%) was identified in mouse brain. In striking contrast, zebrafish had 63 and 70% of CpG methylation in the liver and brain, respectively. This is attributed, in part, to the higher percentage of the zebrafish genome occupied by transposable elements (52% vs. 45% in mice). Thus, the species identity was more significant in determining methylome patterning than was the similarity in organ function. Conserved features of the methylome across tissues and species was the exclusion of methylation from promoters and from CpG islands near transcription start sites, and the clustering of methylated CpGs in gene bodies and intragenic regions. These data suggest that DNA methylation reflects species-specific genome structure, and supports the notion that DNA methylation in non-promoter regions may contribute to genome evolution.

DNA cytosine hydroxymethylation levels are distinct among non-overlapping classes of peripheral blood leukocytes

BACKGROUND

Peripheral blood leukocytes are the most commonly used surrogates to study epigenome-induced risk and epigenomic response to disease-related stress. We considered the hypothesis that the various classes of peripheral leukocytes differentially regulate the synthesis of 5-methylcytosine (5mCG) and its removal via Ten-Eleven Translocation (TET) dioxygenase catalyzed hydroxymethylation to 5-hydroxymethylcytosine (5hmCG), reflecting their responsiveness to environment. Although it is known that reductions in TET1 and/or TET2 activity lead to the over-proliferation of various leukocyte precursors in bone marrow and in development of chronic myelomonocytic leukemia and myeloproliferative neoplasms, the role of 5mCG hydroxymethylation in peripheral blood is less well studied.

RESULTS

We developed simplified protocols to rapidly and reiteratively isolate non-overlapping leukocyte populations from a single small sample of fresh or frozen whole blood. Among peripheral leukocyte types we found extreme variation in the levels of transcripts encoding proteins involved in cytosine methylation (DNMT1, 3A, 3B), the turnover of 5mC by demethylation (TET1, 2, 3), and DNA repair (GADD45A, B, G) and in the global and gene-region-specific levels of DNA 5hmCG (CD4+ T cells≫CD14+ monocytes>CD16+ neutrophils>CD19+ B cells>CD56+ NK cells>Siglec8+ eosinophils>CD8+ T cells).

CONCLUSIONS

Our data taken together suggest a potential hierarchy of responsiveness among classes of leukocytes with CD4+, CD8+ T cells and CD14+ monocytes being the most distinctly poised for a rapid methylome response to physiological stress and disease.