The DNA methylome of human sperm is distinct from blood with little evidence for tissue-consistent obesity associations.

Epidemiological research suggests that paternal obesity may increase the risk of fathering small for gestational age offspring. Studies in non-human mammals indicate that such associations could be mediated by DNA methylation changes in spermatozoa that influence offspring development in utero. Human obesity is associated with differential DNA methylation in peripheral blood. It is unclear, however, whether this differential DNA methylation is reflected in spermatozoa. We profiled genome-wide DNA methylation using the Illumina MethylationEPIC array in a cross-sectional study of matched human blood and sperm from lean (discovery n = 47; replication n = 21) and obese (n = 22) males to analyse tissue covariation of DNA methylation, and identify obesity-associated methylomic signatures. We found that DNA methylation signatures of human blood and spermatozoa are highly discordant, and methylation levels are correlated at only a minority of CpG sites (~1%). At the majority of these sites, DNA methylation appears to be influenced by genetic variation. Obesity-associated DNA methylation in blood was not generally reflected in spermatozoa, and obesity was not associated with altered covariation patterns or accelerated epigenetic ageing in the two tissues. However, one cross-tissue obesity-specific hypermethylated site (cg19357369; chr4:2429884; P = 8.95 × 10-8; 2% DNA methylation difference) was identified, warranting replication and further investigation. When compared to a wide range of human somatic tissue samples (n = 5,917), spermatozoa displayed differential DNA methylation across pathways enriched in transcriptional regulation. Overall, human sperm displays a unique DNA methylation profile that is highly discordant to, and practically uncorrelated with, that of matched peripheral blood. We observed that obesity was only nominally associated with differential DNA methylation in sperm, and therefore suggest that spermatozoal DNA methylation is an unlikely mediator of intergenerational effects of metabolic traits.

Protection against de novo methylation is instrumental in maintaining parent-of-origin methylation inherited from the gametes.

Identifying loci with parental differences in DNA methylation is key to unraveling parent-of-origin phenotypes. By conducting a MeDIP-Seq screen in maternal-methylation free postimplantation mouse embryos (Dnmt3L-/+), we demonstrate that maternal-specific methylation exists very scarcely at midgestation. We reveal two forms of oocyte-specific methylation inheritance: limited to preimplantation, or with longer duration, i.e. maternally imprinted loci. Transient and imprinted maternal germline DMRs (gDMRs) are indistinguishable in gametes and preimplantation embryos, however, de novo methylation of paternal alleles at implantation delineates their fates and acts as a major leveling factor of parent-inherited differences. We characterize two new imprinted gDMRs, at the Cdh15 and AK008011 loci, with tissue-specific imprinting loss, again by paternal methylation gain. Protection against demethylation after fertilization has been emphasized as instrumental in maintaining parent-of-origin methylation inherited from the gametes. Here we provide evidence that protection against de novo methylation acts as an equal major pivot, at implantation and throughout life.

Early life diet conditions the molecular response to post-weaning protein restriction in the mouse.

BACKGROUND: Environmental influences fluctuate throughout the life course of an organism. It is therefore important to understand how the timing of exposure impacts molecular responses. Herein, we examine the responses of two key molecular markers of dietary stress, namely variant-specific methylation at ribosomal DNA (rDNA) and small RNA distribution, including tRNA fragments, in a mouse model of protein restriction (PR) with exposure at pre- and/or post-weaning. RESULTS: We first confirm that pre-weaning PR exposure modulates the methylation state of rDNA in a genotype-dependent manner, whereas post-weaning PR exposure has no such effect. Conversely, post-weaning PR induces a shift in small RNA distribution, but there is no effect in the pre-weaning PR model. Intriguingly, mice exposed to PR throughout their lives show neither of these two dietary stress markers, similar to controls. CONCLUSIONS: The results show that the timing of the insult affects the nature of the molecular response but also, critically, that 'matching' diet exposure either side of weaning eliminates the stress response at the level of rDNA methylation and small RNA in sperm.

Maternal gestational diabetes is associated with genome-wide DNA methylation variation in placenta and cord blood of exposed offspring.

Exposure of a developing foetus to maternal gestational diabetes (GDM) has been shown to programme future risk of diabetes and obesity. Epigenetic variation in foetal tissue may have a mechanistic role in metabolic disease programming through interaction of the pregnancy environment with gene function. We aimed to identify genome-wide DNA methylation variation in cord blood and placenta from offspring born to mothers with and without GDM. Pregnant women of South Asian origin were studied and foetal tissues sampled at term delivery. The Illumina HumanMethylation450 BeadChip was used to assay genome-wide DNA methylation in placenta and cord blood from 27 GDM exposed and 21 unexposed offspring. We identified 1485 cord blood and 1708 placenta methylation variable positions (MVPs) achieving genome-wide significance (adjusted P-value <0.05) with methylation differences of >5%. MVPs were disproportionately located within first exons. A bioinformatic co-methylation algorithm was used to detect consistent directionality of methylation in 1000 bp window around each MVP was observed at 74% of placenta and 59% of cord blood MVPs. KEGG pathway analysis showed enrichment of pathways involved in endocytosis, MAPK signalling and extracellular triggers to intracellular metabolic processes. Replication studies should integrate genomics and transcriptomics with longitudinal sampling to elucidate stability, determine causality for translation into biomarker and prevention studies.

Epigenome-wide association studies for common human diseases.

Despite the success of genome-wide association studies (GWASs) in identifying loci associated with common diseases, a substantial proportion of the causality remains unexplained. Recent advances in genomic technologies have placed us in a position to initiate large-scale studies of human disease-associated epigenetic variation, specifically variation in DNA methylation. Such epigenome-wide association studies (EWASs) present novel opportunities but also create new challenges that are not encountered in GWASs. We discuss EWAS design, cohort and sample selections, statistical significance and power, confounding factors and follow-up studies. We also discuss how integration of EWASs with GWASs can help to dissect complex GWAS haplotypes for functional analysis.

Sexually dimorphic gene expression emerges with embryonic genome activation and is dynamic throughout development.

BACKGROUND: As sex determines mammalian development, understanding the nature and developmental dynamics of the sexually dimorphic transcriptome is important. To explore this, we generated 76 genome-wide RNA-seq profiles from mouse eight-cell embryos, late gestation and adult livers, together with 4 ground-state pluripotent embryonic (ES) cell lines from which we generated both RNA-seq and multiple ChIP-seq profiles. We complemented this with previously published data to yield 5 snap-shots of pre-implantation development, late-gestation placenta and somatic tissue and multiple adult tissues for integrative analysis. RESULTS: We define a high-confidence sex-dimorphic signature of 69 genes in eight-cell embryos. Sex-chromosome-linked components of this signature are largely conserved throughout pre-implantation development and in ES cells, whilst the autosomal component is more dynamic. Sex-biased gene expression is reflected by enrichment for activating and repressive histone modifications. The eight-cell signature is largely non-overlapping with that defined from fetal liver, neither was it correlated with adult liver or other tissues analysed. The number of sex-dimorphic genes increases throughout development. We identified many more dimorphic genes in adult compared to fetal liver. However, approximately two thirds of the dimorphic genes identified in fetal liver were also dimorphic in adult liver. Sex-biased expression differences unique to adult liver were enriched for growth hormone-responsiveness. Sexually dimorphic gene expression in pre-implantation development is driven by sex-chromosome based transcription, whilst later development is characterised by sex dimorphic autosomal transcription. CONCLUSION: This systematic study identifies three distinct phases of sex dimorphism throughout mouse development, and has significant implications for understanding the developmental origins of sex-specific phenotypes and disease in mammals.

A functional methylome map of ulcerative colitis.

The etiology of inflammatory bowel diseases is only partially explained by the current genetic risk map. It is hypothesized that environmental factors modulate the epigenetic landscape and thus contribute to disease susceptibility, manifestation, and progression. To test this, we analyzed DNA methylation (DNAm), a fundamental mechanism of epigenetic long-term modulation of gene expression. We report a three-layer epigenome-wide association study (EWAS) using intestinal biopsies from 10 monozygotic twin pairs (n = 20 individuals) discordant for manifestation of ulcerative colitis (UC). Genome-wide expression scans were generated using Affymetrix UG 133 Plus 2.0 arrays (layer 1). Genome-wide DNAm scans were carried out using Illumina 27k Infinium Bead Arrays to identify methylation variable positions (MVPs, layer 2), and MeDIP-chip on Nimblegen custom 385k Tiling Arrays to identify differentially methylated regions (DMRs, layer 3). Identified MVPs and DMRs were validated in two independent patient populations by quantitative real-time PCR and bisulfite-pyrosequencing (n = 185). The EWAS identified 61 disease-associated loci harboring differential DNAm in cis of a differentially expressed transcript. All constitute novel candidate risk loci for UC not previously identified by GWAS. Among them are several that have been functionally implicated in inflammatory processes, e.g., complement factor CFI, the serine protease inhibitor SPINK4, and the adhesion molecule THY1 (also known as CD90). Our study design excludes nondisease inflammation as a cause of the identified changes in DNAm. This study represents the first replicated EWAS of UC integrated with transcriptional signatures in the affected tissue and demonstrates the power of EWAS to uncover unexplained disease risk and molecular events of disease manifestation.

Guthrie card methylomics identifies temporally stable epialleles that are present at birth in humans.

A major concern in common disease epigenomics is distinguishing causal from consequential epigenetic variation. One means of addressing this issue is to identify the temporal origins of epigenetic variants via longitudinal analyses. However, prospective birth-cohort studies are expensive and time consuming. Here, we report DNA methylomics of archived Guthrie cards for the retrospective longitudinal analyses of in-utero-derived DNA methylation variation. We first validate two methodologies for generating comprehensive DNA methylomes from Guthrie cards. Then, using an integrated epigenomic/genomic analysis of Guthrie cards and follow-up samplings, we identify interindividual DNA methylation variation that is present both at birth and 3 yr later. These findings suggest that disease-relevant epigenetic variation could be detected at birth, i.e., before overt clinical disease. Guthrie card methylomics offers a potentially powerful and cost-effective strategy for studying the dynamics of interindividual epigenomic variation in a range of common human diseases.

DNA methylation profiling of human chromosomes 6, 20 and 22.

DNA methylation is the most stable type of epigenetic modification modulating the transcriptional plasticity of mammalian genomes. Using bisulfite DNA sequencing, we report high-resolution methylation profiles of human chromosomes 6, 20 and 22, providing a resource of about 1.9 million CpG methylation values derived from 12 different tissues. Analysis of six annotation categories showed that evolutionarily conserved regions are the predominant sites for differential DNA methylation and that a core region surrounding the transcriptional start site is an informative surrogate for promoter methylation. We find that 17% of the 873 analyzed genes are differentially methylated in their 5' UTRs and that about one-third of the differentially methylated 5' UTRs are inversely correlated with transcription. Despite the fact that our study controlled for factors reported to affect DNA methylation such as sex and age, we did not find any significant attributable effects. Our data suggest DNA methylation to be ontogenetically more stable than previously thought.

Comparative methylomics reveals gene-body H3K36me3 in Drosophila predicts DNA methylation and CpG landscapes in other invertebrates.

In invertebrates that harbor functional DNA methylation enzymatic machinery, gene-bodies are the primary targets for CpG methylation. However, virtually all other aspects of invertebrate DNA methylation have remained a mystery until now. Here, using a comparative methylomics approach, we demonstrate that Nematostella vectensis, Ciona intestinalis, Apis mellifera, and Bombyx mori show two distinct populations of genes differentiated by gene-body CpG density. Genome-scale DNA methylation profiles for A. mellifera spermatozoa reveal CpG-poor genes are methylated in the germline, as predicted by the depletion of CpGs. We find an evolutionarily conserved distinction between CpG-poor and GpC-rich genes: The former are associated with basic biological processes, the latter with more specialized functions. This distinction is strikingly similar to that recently observed between euchromatin-associated genes in Drosophila that contain intragenic histone 3 lysine 36 trimethylation (H3K36me3) and those that do not, even though Drosophila does not display CpG density bimodality or methylation. We confirm that a significant number of CpG-poor genes in N. vectensis, C. intestinalis, A. mellifera, and B. mori are orthologs of H3K36me3-rich genes in Drosophila. We propose that over evolutionary time, gene-body H3K36me3 has influenced gene-body DNA methylation levels and, consequently, the gene-body CpG density bimodality characteristic of invertebrates that harbor CpG methylation.

Comparative methylome analysis of benign and malignant peripheral nerve sheath tumors.

Aberrant DNA methylation (DNAm) was first linked to cancer over 25 yr ago. Since then, many studies have associated hypermethylation of tumor suppressor genes and hypomethylation of oncogenes to the tumorigenic process. However, most of these studies have been limited to the analysis of promoters and CpG islands (CGIs). Recently, new technologies for whole-genome DNAm (methylome) analysis have been developed, enabling unbiased analysis of cancer methylomes. By using MeDIP-seq, we report a sequencing-based comparative methylome analysis of malignant peripheral nerve sheath tumors (MPNSTs), benign neurofibromas, and normal Schwann cells. Analysis of these methylomes revealed a complex landscape of DNAm alterations. In contrast to what has been reported for other tumor types, no significant global hypomethylation was observed in MPNSTs using methylome analysis by MeDIP-seq. However, a highly significant (P < 10(-100)) directional difference in DNAm was found in satellite repeats, suggesting these repeats to be the main target for hypomethylation in MPNSTs. Comparative analysis of the MPNST and Schwann cell methylomes identified 101,466 cancer-associated differentially methylated regions (cDMRs). Analysis showed these cDMRs to be significantly enriched for two satellite repeat types (SATR1 and ARLα) and suggests an association between aberrant DNAm of these sequences and transition from healthy cells to malignant disease. Significant enrichment of hypermethylated cDMRs in CGI shores (P < 10(-60)), non-CGI-associated promoters (P < 10(-4)) and hypomethylated cDMRs in SINE repeats (P < 10(-100)) was also identified. Integration of DNAm and gene expression data showed that the expression pattern of genes associated with CGI shore cDMRs was able to discriminate between disease phenotypes. This study establishes MeDIP-seq as an effective method to analyze cancer methylomes.

Identification of type 1 diabetes-associated DNA methylation variable positions that precede disease diagnosis.

Monozygotic (MZ) twin pair discordance for childhood-onset Type 1 Diabetes (T1D) is ∼50%, implicating roles for genetic and non-genetic factors in the aetiology of this complex autoimmune disease. Although significant progress has been made in elucidating the genetics of T1D in recent years, the non-genetic component has remained poorly defined. We hypothesized that epigenetic variation could underlie some of the non-genetic component of T1D aetiology and, thus, performed an epigenome-wide association study (EWAS) for this disease. We generated genome-wide DNA methylation profiles of purified CD14+ monocytes (an immune effector cell type relevant to T1D pathogenesis) from 15 T1D-discordant MZ twin pairs. This identified 132 different CpG sites at which the direction of the intra-MZ pair DNA methylation difference significantly correlated with the diabetic state, i.e. T1D-associated methylation variable positions (T1D-MVPs). We confirmed these T1D-MVPs display statistically significant intra-MZ pair DNA methylation differences in the expected direction in an independent set of T1D-discordant MZ pairs (P = 0.035). Then, to establish the temporal origins of the T1D-MVPs, we generated two further genome-wide datasets and established that, when compared with controls, T1D-MVPs are enriched in singletons both before (P = 0.001) and at (P = 0.015) disease diagnosis, and also in singletons positive for diabetes-associated autoantibodies but disease-free even after 12 years follow-up (P = 0.0023). Combined, these results suggest that T1D-MVPs arise very early in the etiological process that leads to overt T1D. Our EWAS of T1D represents an important contribution toward understanding the etiological role of epigenetic variation in type 1 diabetes, and it is also the first systematic analysis of the temporal origins of disease-associated epigenetic variation for any human complex disease.

Ribosomal DNA copy number variation associates with hematological profiles and renal function in the UK Biobank.

The phenotypic impact of genetic variation of repetitive features in the human genome is currently understudied. One such feature is the multi-copy 47S ribosomal DNA (rDNA) that codes for rRNA components of the ribosome. Here, we present an analysis of rDNA copy number (CN) variation in the UK Biobank (UKB). From the first release of UKB whole-genome sequencing (WGS) data, a discovery analysis in White British individuals reveals that rDNA CN associates with altered counts of specific blood cell subtypes, such as neutrophils, and with the estimated glomerular filtration rate, a marker of kidney function. Similar trends are observed in other ancestries. A range of analyses argue against reverse causality or common confounder effects, and all core results replicate in the second UKB WGS release. Our work demonstrates that rDNA CN is a genetic influence on trait variance in humans.

Ribosomal DNA copy number is associated with body mass in humans and other mammals.

Body mass results from a complex interplay between genetics and environment. Previous studies of the genetic contribution to body mass have excluded repetitive regions due to the technical limitations of platforms used for population scale studies. Here we apply genome-wide approaches, identifying an association between adult body mass and the copy number (CN) of 47S-ribosomal DNA (rDNA). rDNA codes for the 18 S, 5.8 S and 28 S ribosomal RNA (rRNA) components of the ribosome. In mammals, there are hundreds of copies of these genes. Inter-individual variation in the rDNA CN has not previously been associated with a mammalian phenotype. Here, we show that rDNA CN variation associates with post-pubertal growth rate in rats and body mass index in adult humans. rDNA CN is not associated with rRNA transcription rates in adult tissues, suggesting the mechanistic link occurs earlier in development. This aligns with the observation that the association emerges by early adulthood.

Marmal-aid--a database for Infinium HumanMethylation450.

BACKGROUND: DNA methylation is indispensible for normal human genome function. Currently there is an increasingly large number of DNA methylomic data being released in the public domain allowing for an opportunity to investigate the relationships between the DNA methylome, genome function, and human phenotypes. The Illumina450K is one of the most popular platforms for assessing DNA methylation with over 10,000 samples available in the public domain. However, accessing all this data requires downloading each individual experiment and due to inconsistent annotation, accessing the right data can be a challenge. DESCRIPTION: Here we introduce 'Marmal-aid', the first standardised database for DNA methylation (freely available at http://marmal-aid.org). In Marmal-aid, the majority of publicly available Illumina HumanMethylation450 data is incorporated into a single repository allowing for re-processing of data including normalisation and imputation of missing values. The database is accessible in two ways: (1) Using an R package to allow for incorporation into existing analysis pipelines which can then be easily queried to gain insight into the functionality of certain CpG sites. This is aimed at a bioinformatician with experience in R. (2) Using a graphical interface allowing general biologists to query a pre-defined set of tissues (currently 15) providing a reference database of the methylation state in these tissues for the 450,000 CpG sites profiled by the Illumina HumanMethylation450. CONCLUSION: Marmal-aid is the largest publicly available Illumina HumanMethylation450 methylation database combining Illumina HumanMethylation450 data from a number of sources into a single location with a single common annotation format. This allows for automated extraction using the R package and inclusion into existing analysis pipelines. Marmal-aid also provides a easy to use GUI to visualise methylation data in user defined genomic regions for various reference tissues.

Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domains.

There is a growing realization that some aging-associated phenotypes/diseases have an epigenetic basis. Here, we report the first genome-scale study of epigenomic dynamics during normal human aging. We identify aging-associated differentially methylated regions (aDMRs) in whole blood in a discovery cohort, and then replicate these aDMRs in sorted CD4(+) T-cells and CD14(+) monocytes in an independent cohort, suggesting that aDMRs occur in precursor haematopoietic cells. Further replication of the aDMRs in buccal cells, representing a tissue that originates from a different germ layer compared with blood, demonstrates that the aDMR signature is a multitissue phenomenon. Moreover, we demonstrate that aging-associated DNA hypermethylation occurs predominantly at bivalent chromatin domain promoters. This same category of promoters, associated with key developmental genes, is frequently hypermethylated in cancers and in vitro cell culture, pointing to a novel mechanistic link between aberrant hypermethylation in cancer, aging, and cell culture.

Transgenerational inheritance of non-genetically determined phenotypes.

Inheritance of non-genetic factors permits ancestral environmental history to inform the development of subsequent generations. This form of soft inheritance has been shown in mammals, yet the molecular underpinnings of this phenomenon are poorly understood. In the present article, we focus on gametic inheritance of non-genetic factors, utilizing examples of paternal transmission to explore the core issues that need to be addressed in order to gain greater insight into the molecular mechanisms. Three essential processes are identified: (i) how the environment affects the germline to establish an altered molecular milieu, (ii) the molecular nature of the inherited mark, and (iii) how this affects genome function in the developing embryo to elicit an alternative developmental outcome.

Genetic variation at mouse and human ribosomal DNA influences associated epigenetic states.

BACKGROUND: Ribosomal DNA (rDNA) displays substantial inter-individual genetic variation in human and mouse. A systematic analysis of how this variation impacts epigenetic states and expression of the rDNA has thus far not been performed. RESULTS: Using a combination of long- and short-read sequencing, we establish that 45S rDNA units in the C57BL/6J mouse strain exist as distinct genetic haplotypes that influence the epigenetic state and transcriptional output of any given unit. DNA methylation dynamics at these haplotypes are dichotomous and life-stage specific: at one haplotype, the DNA methylation state is sensitive to the in utero environment, but refractory to post-weaning influences, whereas other haplotypes entropically gain DNA methylation during aging only. On the other hand, individual rDNA units in human show limited evidence of genetic haplotypes, and hence little discernible correlation between genetic and epigenetic states. However, in both species, adjacent units show similar epigenetic profiles, and the overall epigenetic state at rDNA is strongly positively correlated with the total rDNA copy number. Analysis of different mouse inbred strains reveals that in some strains, such as 129S1/SvImJ, the rDNA copy number is only approximately 150 copies per diploid genome and DNA methylation levels are < 5%. CONCLUSIONS: Our work demonstrates that rDNA-associated genetic variation has a considerable influence on rDNA epigenetic state and consequently rRNA expression outcomes. In the future, it will be important to consider the impact of inter-individual rDNA (epi)genetic variation on mammalian phenotypes and diseases.

DNA methylation-mediated down-regulation of DNA methyltransferase-1 (DNMT1) is coincident with, but not essential for, global hypomethylation in human placenta.

The genome of extraembryonic tissue, such as the placenta, is hypomethylated relative to that in somatic tissues. However, the origin and role of this hypomethylation remains unclear. The DNA methyltransferases DNMT1, -3A, and -3B are the primary mediators of the establishment and maintenance of DNA methylation in mammals. In this study, we investigated promoter methylation-mediated epigenetic down-regulation of DNMT genes as a potential regulator of global methylation levels in placental tissue. Although DNMT3A and -3B promoters lack methylation in all somatic and extraembryonic tissues tested, we found specific hypermethylation of the maintenance DNA methyltransferase (DNMT1) gene and found hypomethylation of the DNMT3L gene in full term and first trimester placental tissues. Bisulfite DNA sequencing revealed monoallelic methylation of DNMT1, with no evidence of imprinting (parent of origin effect). In vitro reporter experiments confirmed that DNMT1 promoter methylation attenuates transcriptional activity in trophoblast cells. However, global hypomethylation in the absence of DNMT1 down-regulation is apparent in non-primate placentas and in vitro derived human cytotrophoblast stem cells, suggesting that DNMT1 down-regulation is not an absolute requirement for genomic hypomethylation in all instances. These data represent the first demonstration of methylation-mediated regulation of the DNMT1 gene in any system and demonstrate that the unique epigenome of the human placenta includes down-regulation of DNMT1 with concomitant hypomethylation of the DNMT3L gene. This strongly implicates epigenetic regulation of the DNMT gene family in the establishment of the unique epigenetic profile of extraembryonic tissue in humans.

The methylome: approaches for global DNA methylation profiling.

DNA methylation plays a critical role in genome function both in health and disease. Almost 60 years after the discovery of 5-methyl cytosine and approximately 25 years since the discovery that altered DNA methylation plays a role in disease, the first high-resolution DNA methylation profile (or methylome) of any genome--Arabidopsis thaliana--was determined. Although only approximately 20% of the typical size of mammalian genomes, this milestone demonstrated that the methylomes of the human and similarly large genomes are now within reach. Here, we review current and emerging technologies that hold promise to deliver the first mammalian methylome and to facilitate comprehensive profiling of essentially any cell type in the context of development, disease and the environment.