Dynamic changes in the human methylome during differentiation

Human epigenome data reveal increased CpG methylation in alternatively spliced sites and putative exonic splicing enhancers

The role of gene body methylation, which represents a major part of methylation in DNA, remains mostly unknown. Evidence based on the CpG distribution associates its presence with nucleosome positioning and alternative splicing. Recently, it was also shown that cytosine methylation influences splicing. However, to date, there is no methylation-based data on the association of methylation with alternative splicing and the distribution in exonic splicing enhancers (ESEs). We presently report that, based on the computational analysis of the Human Epigenome Project data, CpG hypermethylation (>80%) is frequent in alternatively spliced sites (particularly in noncanonical) but not in alternate promoters. The methylation frequency increases in sequences containing multiple putative ESEs. However, significant differences in the extent of methylation are observed among different ESEs. Specifically, moderate levels of methylation, ranging from 20% to 80%, are frequent in SRp55-binding elements, which are associated with response to extracellular conditions, but not in SF2/ASF, primarily responsible for alternative splicing, or in CpG islands. Finally, methylation is more frequent in the presence of AT repeats and CpGs separated by 10 nucleotides and lower in adjacent CpGs, probably indicating its dependence on helical formations and on the presence of nucleosome positioning-related sequences. In conclusion, our results show the regulation of methylation in ESEs and support its involvement in alternative splicing.

High resolution detection and analysis of CpG dinucleotides methylation using MBD-Seq technology

Methyl-CpG binding domain protein sequencing (MBD-seq) is widely used to survey DNA methylation patterns. However, the optimal experimental parameters for MBD-seq remain unclear and the data analysis remains challenging. In this study, we generated high depth MBD-seq data in MCF-7 cell and developed a bi-asymmetric-Laplace model (BALM) to perform data analysis. We found that optimal efficiency of MBD-seq experiments was achieved by sequencing ∼100 million unique mapped tags from a combination of 500 mM and 1000 mM salt concentration elution in MCF-7 cells. Clonal bisulfite sequencing results showed that the methylation status of each CpG dinucleotides in the tested regions was accurately detected with high resolution using the proposed model. These results demonstrated the combination of MBD-seq and BALM could serve as a useful tool to investigate DNA methylome due to its low cost, high specificity, efficiency and resolution.

Dissecting DNA hypermethylation in cancer

There is compelling evidence to support the importance of DNA methylation alterations in cancer development. Both losses and gains of DNA methylation are observed, thought to contribute pathophysiologically by inactivating tumor suppressor genes, inducing chromosomal instability and ectopically activating gene expression. Lesser known are the causes of aberrant DNA methylation. Recent studies have pointed out that intrinsic gene susceptibility to DNA methylation, environmental factors and gene function all have an intertwined participation in this process. Overall, these data support a deterministic rather than a stochastic mechanism for de novo DNA methylation in cancer. In this review article, we discuss the technologies available to study DNA methylation and the endogenous and exogenous factors that influence the onset of de novo methylation in cancer.

A comparative analysis of DNA methylation across human embryonic stem cell lines

Background

We performed a comparative analysis of the genome-wide DNA methylation profiles from three human embryonic stem cell (HESC) lines. It had previously been shown that HESC lines had significantly higher non-CG methylation than differentiated cells, and we therefore asked whether these sites were conserved across cell lines.

Results

We find that heavily methylated non-CG sites are strongly conserved, especially when found within the motif TACAG. They are enriched in splice sites and are more methylated than other non-CG sites in genes. We next studied the relationship between allele-specific expression and allele-specific methylation. By combining bisulfite sequencing and whole transcriptome shotgun sequencing (RNA-seq) data we identified 1,020 genes that show allele-specific expression, and 14% of CG sites genome-wide have allele-specific methylation. Finally, we asked whether the methylation state of transcription factor binding sites affects the binding of transcription factors. We identified variations in methylation levels at binding sites and found that for several transcription factors the correlation between the methylation at binding sites and gene expression is generally stronger than in the neighboring sequences.

Conclusions

These results suggest a possible but as yet unknown functional role for the highly methylated conserved non-CG sites in the regulation of HESCs. We also identified a novel set of genes that are likely transcriptionally regulated by methylation in an allele-specific manner. The analysis of transcription factor binding sites suggests that the methylation state of cis-regulatory elements impacts the ability of factors to bind and regulate transcription.

Regulation of mammalian DNA methyltransferases: a route to new mechanisms

DNA methyltransferases (DNMTs) establish and maintain DNA methylation patterns at specific regions of the genome, thereby contributing to gene regulation. It is becoming evident that an intricate web of pathways target DNMTs to these genomic regions. Here, we review the understanding of these regulatory mechanisms and provide an overview of the new findings, emphasizing the emerging scenario in which several levels of regulation are coordinated to control DNMTs. The mechanisms involved include the dynamic interplay between interdependent post-translational modifications that regulate DNMTs, post-transcriptional regulation by miRNAs and the emerging role of non-coding RNA in targeting mammalian DNMTs. The analysis of these mechanisms is imperative to the understanding of the role of DNA methylation in regulating gene expression during development and in disease.

ICF syndrome mutations cause a broad spectrum of biochemical defects in DNMT3B-mediated de novo DNA methylation

The DNMT3B de novo DNA methyltransferase (DNMT) plays a major role in establishing DNA methylation patterns in early mammalian development, but its catalytic mechanism remains poorly characterized. Here, we provide a comprehensive biochemical analysis of human DNMT3B function through the characterization of a series of site-directed DNMT3B variants associated with immunodeficiency, centromere instability, and facial anomalies (ICF) syndrome. Our data reveal several novel and important aspects of DNMT3B function. First, DNMT3B, unlike DNMT3A, requires a DNA cofactor in order to stably bind to S-adenosyl-l-methionine (SAM), suggesting that it proceeds according to an ordered catalytic scheme. Second, ICF mutations cause a broad spectrum of biochemical defects in DNMT3B function, including defects in homo-oligomerization, SAM binding, SAM utilization, and DNA binding. Third, all tested ICF mutations, including the A766P and R840Q variants, result in altered catalytic properties without interfering with DNMT3L-mediated stimulation; this indicates that DNMT3L is not involved in the pathogenesis of ICF syndrome. Finally, our study reveals a novel level of coupling between substrate binding, oligomerization, and catalysis that is likely conserved within the DNMT3 family of enzymes.

Adaptive significance of natural variations in maternal care in rats: a translational perspective

A wealth of data from the last fifty years documents the potency of early life experiences including maternal care on developing offspring. A majority of this research has focused on the developing stress axis and stress-sensitive behaviors in hopes of identifying factors impacting resilience and risk-sensitivity. The power of early life experience to shape later development is profound and has the potential to increase fitness of individuals for their environments. Current findings in a rat maternal care paradigm highlight the complex and dynamic relation between early experiences and a variety of outcomes. In this review we propose adaptive hypotheses for alternate maternal strategies and resulting offspring phenotypes, and suggest means of distinguishing between these hypotheses. We also provide evidence underscoring the critical role of context in interpreting the adaptive significance of early experiences. If our goal is to identify risk-factors relevant to humans, we must better explore the role of the social and physical environment in our basic animal models.

BiQ Analyzer HT: locus-specific analysis of DNA methylation by high-throughput bisulfite sequencing

Bisulfite sequencing is a widely used method for measuring DNA methylation in eukaryotic genomes. The assay provides single-base pair resolution and, given sufficient sequencing depth, its quantitative accuracy is excellent. High-throughput sequencing of bisulfite-converted DNA can be applied either genome wide or targeted to a defined set of genomic loci (e.g. using locus-specific PCR primers or DNA capture probes). Here, we describe BiQ Analyzer HT (http://biq-analyzer-ht.bioinf.mpi-inf.mpg.de/), a user-friendly software tool that supports locus-specific analysis and visualization of high-throughput bisulfite sequencing data. The software facilitates the shift from time-consuming clonal bisulfite sequencing to the more quantitative and cost-efficient use of high-throughput sequencing for studying locus-specific DNA methylation patterns. In addition, it is useful for locus-specific visualization of genome-wide bisulfite sequencing data.

Genome-wide survey reveals dynamic widespread tissue-specific changes in DNA methylation during development

BACKGROUND

Changes in DNA methylation in the mammalian genome during development are frequent events and play major roles regulating gene expression and other developmental processes. It is necessary to identify these events so that we may understand how these changes affect normal development and how aberrant changes may impact disease.

RESULTS

In this study Methylated DNA ImmunoPrecipitation (MeDIP) was used in conjunction with a NimbleGen promoter plus CpG island (CpGi) array to identify Tissue and Developmental Stage specific Differentially Methylated DNA Regions (T-DMRs and DS-DMRs) on a genome-wide basis. Four tissues (brain, heart, liver, and testis) from C57BL/6J mice were analyzed at three developmental stages (15 day embryo, E15; new born, NB; 12 week adult, AD). Almost 5,000 adult T-DMRs and 10,000 DS-DMRs were identified. Surprisingly, almost all DS-DMRs were tissue specific (i.e. methylated in at least one tissue and unmethylated in one or more tissues). In addition our results indicate that many DS-DMRs are methylated at early development stages (E15 and NB) but are unmethylated in adult. There is a very strong bias for testis specific methylation in non-CpGi promoter regions (94%). Although the majority of T-DMRs and DS-DMRs tended to be in non-CpGi promoter regions, a relatively large number were also located in CpGi in promoter, intragenic and intergenic regions (>15% of the 15,979 CpGi on the array).

CONCLUSIONS

Our data suggests the vast majority of unique sequence DNA methylation has tissue specificity, that demethylation has a prominent role in tissue differentiation, and that DNA methylation has regulatory roles in alternative promoter selection and in non-promoter regions. Overall, our studies indicate changes in DNA methylation during development are a dynamic, widespread, and tissue-specific process involving both DNA methylation and demethylation.

DNA methylation in neurodegenerative disorders: a missing link between genome and environment?

The risk of developing neurodegenerative disorders such as Alzheimer's disease or Parkinson's disease is influenced by genetic and environmental factors. Environmental events occurring during development or later in life can be related to disease susceptibility. One way by which the environment may exert its effect is through epigenetic modifications, which might affect the functioning of genes. These include nucleosome positioning, post-translational histone modifications, and DNA methylation. In this review we will focus in the potential role of DNA methylation in neurodegenerative disorders and in the approaches to explore such epigenetic changes. Advances in deciphering the role of epigenetic modifications in phenotype are being uncovered for a variety of diseases, including cancer, autoimmune, neurodevelopmental and cognitive disorders. Epigenetic modifications are now being also associated with cardiovascular and metabolic traits, and they are expected to be especially involved in learning and memory processes, as well as in neurodegenerative disease. The study of the role of methylation and other epigenetic modifications in disease development will provide new insights in the etiopathogenesis of neurodegenerative disorders, and should hopefully shape new avenues in the development of therapeutic strategies.

A user's guide to the encyclopedia of DNA elements (ENCODE)

The mission of the Encyclopedia of DNA Elements (ENCODE) Project is to enable the scientific and medical communities to interpret the human genome sequence and apply it to understand human biology and improve health. The ENCODE Consortium is integrating multiple technologies and approaches in a collective effort to discover and define the functional elements encoded in the human genome, including genes, transcripts, and transcriptional regulatory regions, together with their attendant chromatin states and DNA methylation patterns. In the process, standards to ensure high-quality data have been implemented, and novel algorithms have been developed to facilitate analysis. Data and derived results are made available through a freely accessible database. Here we provide an overview of the project and the resources it is generating and illustrate the application of ENCODE data to interpret the human genome.

Emerging patterns of epigenomic variation

Fuelled by new sequencing technologies, epigenome mapping projects are revealing epigenomic variation at all levels of biological complexity, from species to cells. Comparisons of methylation profiles among species reveal evolutionary conservation of gene body methylation patterns, pointing to the fundamental role of epigenomes in gene regulation. At the human population level, epigenomic changes provide footprints of the effects of genomic variants within the vast nonprotein-coding fraction of the genome, and comparisons of the epigenomes of parents and their offspring point to quantitative epigenomic parent-of-origin effects confounding classical Mendelian genetics. At the organismal level, comparisons of epigenomes from diverse cell types provide insights into cellular differentiation. Finally, comparisons of epigenomes from monozygotic twins help dissect genetic and environmental influences on human phenotypes and longitudinal comparisons reveal aging-associated epigenomic drift. The development of new bioinformatic frameworks for comparative epigenome analysis is putting epigenome maps within the reach of researchers across a wide spectrum of biological disciplines.

Exposing the DNA methylome iceberg

DNA methylation was the first epigenetic modification discovered. Until recently, comprehensive coverage of the composition and distribution of methylated cytosines across the genome was lacking. Technological advances, however, are providing methylation maps that can reveal the genomic distribution of DNA methylation in different cell states or phenotypes. The emerging picture includes extensive gene body methylation that is highly conserved in eukaryotes, the presence of DNA methylation in previously unappreciated sequence contexts, and the discovery of another modified DNA base, 5-hydroxymethylcytosine. These new data point to the role of DNA methylation both in gene silencing and gene activation; reconciliation of these seemingly contradictory roles will be essential to fully unravel the biological function of DNA methylation in eukaryotes. Here we review how these recently exposed features of the DNA methylome are challenging previously held dogmas in the field.

Epigenetic modifications in plants: an evolutionary perspective

Plant genomes are modified by an array of epigenetic marks that help regulate plant growth and reproduction. Although plants share many epigenetic features with animals and fungi, some epigenetic marks are unique to plants. In different organisms, the same epigenetic mark can play different roles and/or similar functions can be carried out by different epigenetic marks. Furthermore, while the enzymatic systems responsible for generating or eliminating epigenetic marks are often conserved, there are also cases where they are quite divergent between plants and other organisms. DNA methylation and methylation of histone tails on the lysine 4, 9, and 27 positions are among the best characterized epigenetic marks in both plants and animals. Recent studies have greatly enhanced our knowledge about the pattern of these marks in various genomes and provided insights into how they are established and maintained and how they function. This review focuses on the conservation and divergence of the pathways that mediate these four types of epigenetic marks.

Accurate genome-scale percentage DNA methylation estimates from microarray data

DNA methylation is a key regulator of gene function in a multitude of both normal and abnormal biological processes, but tools to elucidate its roles on a genome-wide scale are still in their infancy. Methylation sensitive restriction enzymes and microarrays provide a potential high-throughput, low-cost platform to allow methylation profiling. However, accurate absolute methylation estimates have been elusive due to systematic errors and unwanted variability. Previous microarray preprocessing procedures, mostly developed for expression arrays, fail to adequately normalize methylation-related data since they rely on key assumptions that are violated in the case of DNA methylation. We develop a normalization strategy tailored to DNA methylation data and an empirical Bayes percentage methylation estimator that together yield accurate absolute methylation estimates that can be compared across samples. We illustrate the method on data generated to detect methylation differences between tissues and between normal and tumor colon samples.

DNA methylation and demethylation in mammals

Cell type-specific DNA methylation patterns are established during mammalian development and maintained in adult somatic cells. Understanding how these patterns of 5-methylcytosine are established and maintained requires the elucidation of mechanisms for both DNA methylation and demethylation. The enzymes involved in the de novo methylation of DNA and the maintenance of the resulting methylation patterns have been fairly well characterized. However, important remaining challenges are to understand how DNA methylation systems function in vivo and in the context of chromatin. In addition, the enzymes and mechanisms for demethylation remain to be elucidated. There is still no consensus as to how active enzymatic demethylation is achieved in mammalian cells, but recent studies implicate base excision repair for genome-wide DNA demethylation in germ cells and early embryos.

Genome-wide promoter DNA methylation dynamics of human hematopoietic progenitor cells during differentiation and aging

DNA methylation plays an important role in the self-renewal of hematopoietic stem cells and in the commitment to the lymphoid or myeloid lineages. Using purified CD34⁺ hematopoietic progenitor cells and differentiated myeloid cell populations from the same human samples, we obtained detailed methylation profiles at distinct stages of hematopoiesis. We identified a defined set of differentiation-related genes that are methylated in CD34⁺ hematopoietic progenitor cells but show pronounced DNA hypomethylation in monocytes and in granulocytes. In addition, by comparing hematopoietic progenitor cells from umbilical cord blood to hematopoietic progenitor cells from peripheral blood of adult donors we were also able to analyze age-related methylation changes in CD34⁺ cells. Interestingly, the methylation changes observed in older progenitor cells showed a bimodal pattern with hypomethylation of differentiation-associated genes and de novo methylation events resembling epigenetic mutations. Our results thus provide detailed insight into the methylation dynamics during differentiation and suggest that epigenetic changes contribute to hematopoietic progenitor cell aging.

Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling

Genome-wide mapping of 5-methylcytosine is of broad interest to many fields of biology and medicine. A variety of methods have been developed, and several have recently been advanced to genome-wide scale using arrays and next-generation sequencing approaches. We have previously reported reduced representation bisulfite sequencing (RRBS), a bisulfite-based protocol that enriches CG-rich parts of the genome, thereby reducing the amount of sequencing required while capturing the majority of promoters and other relevant genomic regions. The approach provides single-nucleotide resolution, is highly sensitive and provides quantitative DNA methylation measurements. This protocol should enable any standard molecular biology laboratory to generate RRBS libraries of high quality. Briefly, purified genomic DNA is digested by the methylation-insensitive restriction enzyme MspI to generate short fragments that contain CpG dinucleotides at the ends. After end-repair, A-tailing and ligation to methylated Illumina adapters, the CpG-rich DNA fragments (40-220 bp) are size selected, subjected to bisulfite conversion, PCR amplified and end sequenced on an Illumina Genome Analyzer. Note that alignment and analysis of RRBS sequencing reads are not covered in this protocol. The extremely low input requirements (10-300 ng), the applicability of the protocol to formalin-fixed and paraffin-embedded samples, and the technique's single-nucleotide resolution extends RRBS to a wide range of biological and clinical samples and research applications. The entire process of RRBS library construction takes ∼9 d.

Methylation-mediated deamination of 5-methylcytosine appears to give rise to mutations causing human inherited disease in CpNpG trinucleotides, as well as in CpG dinucleotides

The cytosine-guanine (CpG) dinucleotide has long been known to be a hotspot for pathological mutation in the human genome. This hypermutability is related to its role as the major site of cytosine methylation with the attendant risk of spontaneous deamination of 5-methylcytosine (5mC) to yield thymine. Cytosine methylation, however, also occurs in the context of CpNpG sites in the human genome, an unsurprising finding since the intrinsic symmetry of CpNpG renders it capable of supporting a semi-conservative model of replication of the methylation pattern. Recently, it has become clear that significant DNA methylation occurs in a CpHpG context (where H = A, C or T) in a variety of human somatic tissues. If we assume that CpHpG methylation also occurs in the germline, and that 5mC deamination can occur within a CpHpG context, then we might surmise that methylated CpHpG sites could also constitute mutation hotspots causing human genetic disease. To test this postulate, 54,625 missense and nonsense mutations from 2,113 genes causing inherited disease were retrieved from the Human Gene Mutation Database (http://www.hgmd.org). Some 18.2 per cent of these pathological lesions were found to be C → T and G → A transitions located in CpG dinucleotides (compatible with a model of methylation-mediated deamination of 5mC), an approximately ten-fold higher proportion than would have been expected by chance alone. The corresponding proportion for the CpHpG trinucleotide was 9.9 per cent, an approximately two-fold higher proportion than would have been expected by chance. We therefore estimate that ∼5 per cent of missense/nonsense mutations causing human inherited disease may be attributable to methylation-mediated deamination of 5mC within a CpHpG context.

Epigenomic communication systems in humans and honey bees: from molecules to behavior

A 2010 Nature editorial entitled "Time for the Epigenome" trumpets the appearance of the International Human Epigenome Consortium and likens it to Biology's equivalent of the Large Hadron Collider. It strongly endorses the viewpoint that selective modifications of "marks" on DNA and histones constitute the crucial codes of life, a proposition which is hotly contested (Ptashne et al., in 2010). This proposition reflects the current mindset that DNA and histone modifications are the prime movers in gene regulation during evolution. This claim is perplexing, since the well characterized organisms, Drosophila melanogaster and Caenorhabditis elegans, lack methylated DNA "marks" and the DNA methytransferase enzymology. Despite their complete absence, D. melanogaster nevertheless has extensive gene regulatory networks which drive sophisticated development, gastrulation, migration of germ cells and yield a nervous system with significant neural attributes. In stark contrast, the honey bee Apis mellifera deploys its human-type DNA methyltransferase enzymology to "mark" its DNA and it too has sophisticated development. What roles therefore is DNA methylation playing in different animals? The honey bee brings a fresh perspective to this question. Its combinatorial chemistry of pheromones, tergal and cuticular exudates provide an exquisite communication system between thousands of individuals. The development of queen and worker is strictly controlled by differential feeding of royal jelly and their adult behaviors are accompanied by epigenomic changes. Their interfaces with different "environments" are extensive, allowing an evaluation of the roles of epigenomes in behavior in a natural environment, in the space of a few weeks, and at requisite levels of experimental rigor.

Epigenetic and phenotypic profile of fibroblasts derived from induced pluripotent stem cells

Human induced pluripotent stem (hiPS) cells offer a novel source of patient-specific cells for regenerative medicine. However, the biological potential of iPS-derived cells and their similarities to cells differentiated from human embryonic stem (hES) cells remain unclear. We derived fibroblast-like cells from two hiPS cell lines and show that their phenotypic properties and patterns of DNA methylation were similar to that of mature fibroblasts and to fibroblasts derived from hES cells. iPS-derived fibroblasts (iPDK) and their hES-derived counterparts (EDK) showed similar cell morphology throughout differentiation, and patterns of gene expression and cell surface markers were characteristic of mature fibroblasts. Array-based methylation analysis was performed for EDK, iPDK and their parental hES and iPS cell lines, and hierarchical clustering revealed that EDK and iPDK had closely-related methylation profiles. DNA methylation analysis of promoter regions associated with extracellular matrix (ECM)-production (COL1A1) by iPS- and hESC-derived fibroblasts and fibroblast lineage commitment (PDGFRβ), revealed promoter demethylation linked to their expression, and patterns of transcription and methylation of genes related to the functional properties of mature stromal cells were seen in both hiPS- and hES-derived fibroblasts. iPDK cells also showed functional properties analogous to those of hES-derived and mature fibroblasts, as seen by their capacity to direct the morphogenesis of engineered human skin equivalents. Characterization of the functional behavior of ES- and iPS-derived fibroblasts in engineered 3D tissues demonstrates the utility of this tissue platform to predict the capacity of iPS-derived cells before their therapeutic application.

QDMR: a quantitative method for identification of differentially methylated regions by entropy

DNA methylation plays critical roles in transcriptional regulation and chromatin remodeling. Differentially methylated regions (DMRs) have important implications for development, aging and diseases. Therefore, genome-wide mapping of DMRs across various temporal and spatial methylomes is important in revealing the impact of epigenetic modifications on heritable phenotypic variation. We present a quantitative approach, quantitative differentially methylated regions (QDMRs), to quantify methylation difference and identify DMRs from genome-wide methylation profiles by adapting Shannon entropy. QDMR was applied to synthetic methylation patterns and methylation profiles detected by methylated DNA immunoprecipitation microarray (MeDIP-chip) in human tissues/cells. This approach can give a reasonable quantitative measure of methylation difference across multiple samples. Then DMR threshold was determined from methylation probability model. Using this threshold, QDMR identified 10 651 tissue DMRs which are related to the genes enriched for cell differentiation, including 4740 DMRs not identified by the method developed by Rakyan et al. QDMR can also measure the sample specificity of each DMR. Finally, the application to methylation profiles detected by reduced representation bisulphite sequencing (RRBS) in mouse showed the platform-free and species-free nature of QDMR. This approach provides an effective tool for the high-throughput identification of potential functional regions involved in epigenetic regulation.

A novel approach identifies new differentially methylated regions (DMRs) associated with imprinted genes

Imprinted genes are critical for normal human growth and neurodevelopment. They are characterized by differentially methylated regions (DMRs) of DNA that confer parent of origin-specific transcription. We developed a new strategy to identify imprinted gene-associated DMRs. Using genome-wide methylation profiling of sodium bisulfite modified DNA from normal human tissues of biparental origin, candidate DMRs were identified by selecting CpGs with methylation levels consistent with putative allelic differential methylation. In parallel, the methylation profiles of tissues of uniparental origin, i.e., paternally-derived androgenetic complete hydatidiform moles (AnCHMs), and maternally-derived mature cystic ovarian teratoma (MCT), were examined and then used to identify CpGs with parent of origin-specific DNA methylation. With this approach, we found known DMRs associated with imprinted genomic regions as well as new DMRs for known imprinted genes, NAP1L5 and ZNF597, and novel candidate imprinted genes. The paternally methylated DMR for one candidate, AXL, a receptor tyrosine kinase, was also validated in experiments with mouse embryos that demonstrated Axl was expressed preferentially from the maternal allele in a DNA methylation-dependent manner.

Hotspots of aberrant epigenomic reprogramming in human induced pluripotent stem cells

Induced pluripotent stem cells (iPSCs) offer immense potential for regenerative medicine and studies of disease and development. Somatic cell reprogramming involves epigenomic reconfiguration, conferring iPSCs with characteristics similar to embryonic stem (ES) cells. However, it remains unknown how complete the reestablishment of ES-cell-like DNA methylation patterns is throughout the genome. Here we report the first whole-genome profiles of DNA methylation at single-base resolution in five human iPSC lines, along with methylomes of ES cells, somatic cells, and differentiated iPSCs and ES cells. iPSCs show significant reprogramming variability, including somatic memory and aberrant reprogramming of DNA methylation. iPSCs share megabase-scale differentially methylated regions proximal to centromeres and telomeres that display incomplete reprogramming of non-CG methylation, and differences in CG methylation and histone modifications. Lastly, differentiation of iPSCs into trophoblast cells revealed that errors in reprogramming CG methylation are transmitted at a high frequency, providing an iPSC reprogramming signature that is maintained after differentiation.

Different binding properties and function of CXXC zinc finger domains in Dnmt1 and Tet1

Several mammalian proteins involved in chromatin and DNA modification contain CXXC zinc finger domains. We compared the structure and function of the CXXC domains in the DNA methyltransferase Dnmt1 and the methylcytosine dioxygenase Tet1. Sequence alignment showed that both CXXC domains have a very similar framework but differ in the central tip region. Based on the known structure of a similar MLL1 domain we developed homology models and designed expression constructs for the isolated CXXC domains of Dnmt1 and Tet1 accordingly. We show that the CXXC domain of Tet1 has no DNA binding activity and is dispensable for catalytic activity in vivo. In contrast, the CXXC domain of Dnmt1 selectively binds DNA substrates containing unmethylated CpG sites. Surprisingly, a Dnmt1 mutant construct lacking the CXXC domain formed covalent complexes with cytosine bases both in vitro and in vivo and rescued DNA methylation patterns in dnmt1^−/− embryonic stem cells (ESCs) just as efficiently as wild type Dnmt1. Interestingly, neither wild type nor ΔCXXC Dnmt1 re-methylated imprinted CpG sites of the H19a promoter in dnmt1^−/− ESCs, arguing against a role of the CXXC domain in restraining Dnmt1 methyltransferase activity on unmethylated CpG sites.

Conservation and divergence of DNA methylation in eukaryotes: new insights from single base-resolution DNA methylomes

DNA methylation is one of the most important heritable epigenetic modifications of the genome and is involved in the regulation of many cellular processes. Aberrant DNA methylation has been frequently reported to influence gene expression and subsequently cause various human diseases, including cancer. Recent rapid advances in next-generation sequencing technologies have enabled investigators to profile genome methylation patterns at single-base resolution. Remarkably, more than 20 eukaryotic methylomes have been generated thus far, with a majority published since November 2009. Analysis of this vast amount of data has dramatically enriched our knowledge of biological function, conservation and divergence of DNA methylation in eukaryotes. Even so, many specific functions of DNA methylation and their underlying regulatory systems still remain unknown to us. Here, we briefly introduce current approaches for DNA methylation profiling and then systematically review the features of whole genome DNA methylation patterns in eight animals, six plants and five fungi. Our systematic comparison provides new insights into the conservation and divergence of DNA methylation in eukaryotes and their regulation of gene expression. This work aims to summarize the current state of available methylome data and features informatively.

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.

Detection and removal of biases in the analysis of next-generation sequencing reads

Since the emergence of next-generation sequencing (NGS) technologies, great effort has been put into the development of tools for analysis of the short reads. In parallel, knowledge is increasing regarding biases inherent in these technologies. Here we discuss four different biases we encountered while analyzing various Illumina datasets. These biases are due to both biological and statistical effects that in particular affect comparisons between different genomic regions. Specifically, we encountered biases pertaining to the distributions of nucleotides across sequencing cycles, to mappability, to contamination of pre-mRNA with mRNA, and to non-uniform hydrolysis of RNA. Most of these biases are not specific to one analyzed dataset, but are present across a variety of datasets and within a variety of genomic contexts. Importantly, some of these biases correlated in a highly significant manner with biological features, including transcript length, gene expression levels, conservation levels, and exon-intron architecture, misleadingly increasing the credibility of results due to them. We also demonstrate the relevance of these biases in the context of analyzing an NGS dataset mapping transcriptionally engaged RNA polymerase II (RNAPII) in the context of exon-intron architecture, and show that elimination of these biases is crucial for avoiding erroneous interpretation of the data. Collectively, our results highlight several important pitfalls, challenges and approaches in the analysis of NGS reads.

DNA methylation restricts spontaneous multi-lineage differentiation of mesenchymal progenitor cells, but is stable during growth factor-induced terminal differentiation

The progressive restriction of differentiation potential from pluripotent embryonic stem cells, via multipotent progenitor cells to terminally differentiated, mature somatic cells, involves step-wise changes in transcription patterns that are tightly controlled by the coordinated action of key transcription factors and changes in epigenetic modifications. While previous studies have demonstrated tissue-specific differences in DNA methylation patterns that might function in lineage restriction, it is unclear at what exact developmental stage these differences arise. Here, we have studied whether terminal, multi-lineage differentiation of C2C12 myoblasts is accompanied by lineage-specific changes in DNA methylation patterns. Using bisulfite sequencing and genome-wide methylated DNA- and chromatin immunoprecipitation-on-chip techniques we show that in these cells, in general, myogenic genes are enriched for RNA polymerase II and hypomethylated, whereas osteogenic genes show lower polymerase occupancy and are hypermethylated. Removal of DNA methylation marks by 5-azacytidine (5AC) treatment alters the myogenic lineage commitment of these cells and induces spontaneous osteogenic and adipogenic differentiation. This is accompanied by upregulation of key lineage-specific transcription factors. We subsequently analyzed genome-wide changes in DNA methylation and polymerase II occupancy during BMP2-induced osteogenesis. Our data indicate that BMP2 is able to induce the transcriptional program underlying osteogenesis without changing the methylation status of the genome. We conclude that DNA methylation primes C2C12 cells for myogenesis and prevents spontaneous osteogenesis, but still permits induction of the osteogenic transcriptional program upon BMP2 stimulation. Based on these results, we propose that cell type-specific DNA methylation patterns are established prior to terminal differentiation of adult progenitor cells. This article is part of a Special Issue entitled: 11th European Symposium on Calcium.

Epigenetic modifications and human disease

Epigenetics is one of the most rapidly expanding fields in biology. The recent characterization of a human DNA methylome at single nucleotide resolution, the discovery of the CpG island shores, the finding of new histone variants and modifications, and the unveiling of genome-wide nucleosome positioning maps highlight the accelerating speed of discovery over the past two years. Increasing interest in epigenetics has been accompanied by technological breakthroughs that now make it possible to undertake large-scale epigenomic studies. These allow the mapping of epigenetic marks, such as DNA methylation, histone modifications and nucleosome positioning, which are critical for regulating gene and noncoding RNA expression. In turn, we are learning how aberrant placement of these epigenetic marks and mutations in the epigenetic machinery is involved in disease. Thus, a comprehensive understanding of epigenetic mechanisms, their interactions and alterations in health and disease, has become a priority in biomedical research.

Epigenetic modifications in pluripotent and differentiated cells

Epigenetic modifications constitute a complex regulatory layer on top of the genome sequence. Pluripotent and differentiated cells provide a powerful system for investigating how the epigenetic code influences cellular fate. High-throughput sequencing of these cell types has yielded DNA methylation maps at single-nucleotide resolution and many genome-wide chromatin maps. In parallel to epigenome mapping efforts, remarkable progress has been made in our ability to manipulate cell states; ectopic expression of transcription factors has been shown to override developmentally established epigenetic marks and to enable routine generation of induced pluripotent stem (iPS) cells. Despite these advances, many fundamental questions remain. The roles of epigenetic marks and, in particular, of epigenetic modifiers in development and in disease states are not well understood. Although iPS cells appear molecularly and functionally similar to embryonic stem cells, more genome-wide studies are needed to define the extent and functions of epigenetic remodeling during reprogramming.

Methylation at global LINE-1 repeats in human blood are affected by gender but not by age or natural hormone cycles

Previously, we reported on inter-individual and gender specific variations of LINE-1 methylation in healthy individuals. In this study, we investigated whether this variability could be influenced by age or sex hormones in humans. To this end, we studied LINE-1 methylation in vivo in blood-derived DNA from individuals aged 18 to 64 years and from young healthy females at various hormone levels during the menstrual cycle. Our results show that no significant association with age was observed. However, the previously reported increase of LINE-1 methylation in males was reconfirmed. In females, although no correlation between LINE-1 or Alu methylation and hormone levels was observed, a significant stable individual specific level of methylation was noted. In vitro results largely confirmed these findings, as neither estrogen nor dihydrotestosterone affected LINE-1 or Alu methylation in Hek293T, HUVEC, or MDA-kb2 cell lines. In contrast, a decrease in methylation was observed in estrogen-treated T47-Kbluc cell lines strongly expressing estrogen receptor. The very low expression of estrogen receptor in blood cells could explain the observed insensitivity of methylation at LINE-1 to natural hormonal variations in females. In conclusion, neither natural cycle of hormones nor age has a detectable effect on the LINE-1 methylation in peripheral blood cells, while gender remains an important factor.

The DNMT3 family of mammalian de novo DNA methyltransferases

The deposition of DNA methylation at promoters of transposons, X-linked genes, imprinted genes, and other lineage-specific genes is clearly associated with long-term transcriptional silencing. Thus, DNA methylation represents a key layer of epigenetic information in mammals that is required for embryonic development, germline differentiation, and, as shown more recently, for the function and maturation of neuronal tissues. The DNMT3A, DNMT3B, and DNMT3L proteins are primarily responsible for the establishment of genomic DNA methylation patterns and, as such, play an important role in human developmental, reproductive, and mental health. Progress in our understanding of this important protein family has been rapid in recent years and has been accompanied by stunning developments in the analysis of the human DNA methylome in multiple cell types. This review focuses on recent developments in the characterization of the DNMT3 family of DNA methyltransferases at the biochemical, structural, and functional levels. Interconnections between the DNA-based and histone-based layers of epigenetic information are particularly highlighted, as it is now clear that de novo methylation occurs chiefly in the context of nucleosomal templates.

SIRPH: an HPLC-based SNuPE for quantitative methylation measurement at specific CpG sites

Genome-wide sequence-specific methylation analysis has become a readily available and affordable procedure in epigenetic laboratories. Most of these procedures are based on immunoprecipitation, microarrays, or next generation deep bisulfite sequencing. However, most of these protocols are far from being quantitative. Moreover, abnormal or specific methylation patterns must always be further validated by quantitative sequence-specific methylation analysis. In this chapter, we describe a detailed and simplified protocol (using one universal HPLC gradient for all separations as well as one SNuPE annealing temperature for all primers) for the previously published SIRPH analysis, which is based on the single nucleotide primer extension combined with high-performance liquid chromatography. This method is highly accurate, reproducible, quantitative, and suitable for analysis of as little as 50 ng of PCR product derived from limited starting materials.

Profiling epigenetic alterations in disease

Nowadays, epigenetics is one of the fastest growing research areas in biomedicine. Studies have demonstrated that changes in the epigenome are not only common in cancer, but are also involved in the pathogenesis of noncancerous diseases like immunological, cardiovascular, developmental and neurological/psychiatric disorders. At the same time, during the last years, a technological revolution has taken place in the field of epigenomics, which is defined as the study of epigenetic changes throughout the whole genome. Microarray technologies and more recently, the development of next generation sequencing devices are now providing researchers with tools to draw high-resolution maps of DNA methylation and histone modifications in normal tissues and diseases. This chapter will review the currently available high-throughput techniques for studying the epigenome and their applications for characterizing human diseases.

Genome-wide distribution of DNA methylation at single-nucleotide resolution

DNA methylation, a well-known epigenetic modification in mammalian genomes, is important for development and health. Dysregulation of DNA methylation can cause abnormal gene regulation, leading to anomalous development and diseases. Until recently, the ability to understand the functions and dynamics of DNA methylation was limited by the availability of technologies for comprehensively characterizing methylation on a genome-wide scale. Rapid advances in high-throughput approaches (particularly next-generation sequencing), coupled with molecular techniques, have enabled unbiased genome-wide profiling of DNA modifications at single-base resolution and helped to elucidate their impact on gene regulation. Here, we discuss the development of genomic approaches to decipher the global methylome at single-base resolution, the challenges faced, and the emerging new insights. Our ability to decipher this important epigenetic modification and how it impacts gene expression will provide a framework for understanding numerous disease mechanisms, and suggest means to treat or prevent them in the future.

Evidence for a DNA-Based Mechanism of Intron-Mediated Enhancement

Many introns significantly increase gene expression through a process termed intron-mediated enhancement (IME). Introns exist in the transcribed DNA and the nascent RNA, and could affect expression from either location. To determine which is more relevant to IME, hybrid introns were constructed that contain sequences from stimulating Arabidopsis thaliana introns either in their normal orientation or as the reverse complement. Both ends of each intron are from the non-stimulatory COR15a intron in their normal orientation to allow splicing. The inversions create major alterations to the sequence of the transcribed RNA with relatively minor changes to the DNA structure. Introns containing portions of either the UBQ10 or ATPK1 intron increased expression to a similar degree regardless of orientation. Also, computational predictions of IME improve when both intron strands are considered. These findings are more consistent with models of IME that act at the level of DNA rather than RNA.

DNA methylation of genes in adipose tissue

Body fat distribution plays an important role in determining metabolic health. Whereas central obesity is closely associated with the development of CVD and type 2 diabetes, lower body fat appears to be protective and is paradoxically associated with improved metabolic and cardiovascular profiles. Physiological studies have demonstrated that fatty acid handling differs between white adipose tissue depots, with lower body white adipose tissue acting as a more efficient site for long-term lipid storage. The regulatory mechanisms governing these regional differences in function remain to be elucidated. Although the local microenvironment is likely to be a contributing factor, recent findings point towards the tissues being intrinsically distinct at the level of the adipocyte precursor cells (pre-adipocytes). The multi-potent pre-adipocytes are capable of generating cells of the mesenchymal lineage, including adipocytes. Regional differences in the adipogenic and replicative potential of these cells, as well as metabolic and biochemical activity, have been reported. Intriguingly, the genetic and metabolic characteristics of these cells can be retained through multiple generations when the cells are cultured in vitro. The rapidly emerging field of epigenetics may hold the key for explaining regional differences in white adipose tissue gene expression and function. Epigenetics describes the regulation of gene expression that occurs independently of changes in DNA sequence, for instance, DNA methylation or histone protein modification. In this review, we will discuss the contribution of DNA methylation to the determination of cells of adipogenic fate as well as the role DNA methylation may play during adipocyte terminal differentiation.

Targets and dynamics of promoter DNA methylation during early mouse development

DNA methylation is extensively reprogrammed during the early phases of mammalian development, yet genomic targets of this process are largely unknown. We optimized methylated DNA immunoprecipitation for low numbers of cells and profiled DNA methylation during early development of the mouse embryonic lineage in vivo. We observed a major epigenetic switch during implantation at the transition from the blastocyst to the postimplantation epiblast. During this period, DNA methylation is primarily targeted to repress the germline expression program. DNA methylation in the epiblast is also targeted to promoters of lineage-specific genes such as hematopoietic genes, which are subsequently demethylated during terminal differentiation. De novo methylation during early embryogenesis is catalyzed by Dnmt3b, and absence of DNA methylation leads to ectopic gene activation in the embryo. Finally, we identify nonimprinted genes that inherit promoter DNA methylation from parental gametes, suggesting that escape of post-fertilization DNA methylation reprogramming is prevalent in the mouse genome.

Epigenetic architecture and miRNA: reciprocal regulators

Deregulation of epigenetic and miRNA pathways are emerging as key events in carcinogenesis. miRNA genes can be epigenetically regulated and miRNAs can themselves repress key enzymes that drive epigenetic remodeling. Epigenetic and miRNA functions are thus tightly interconnected and crucial for maintaining correct local and global genomic architecture as well as gene-expression patterns, yet the underlying molecular mechanisms and their widespread effects remain poorly understood. Owing to the tissue specificity, versatility and relative stability of miRNAs, these small ncRNAs are considered especially promising in clinical applications, and their biogenesis and function is subject of active research. In this article, the current status of epigenetic miRNA regulation is summarized and future therapeutic prospects in the field are discussed with a focus on cancer.

The epigenetic regulator Cfp1

Numerous epigenetic modifications have been identified and correlated with transcriptionally active euchromatin or repressed heterochromatin and many enzymes responsible for the addition and removal of these marks have been characterized. However, less is known regarding how these enzymes are regulated and targeted to appropriate genomic locations. Mammalian CXXC finger protein 1 is an epigenetic regulator that was originally identified as a protein that binds specifically to any DNA sequence containing an unmethylated CpG dinucleotide. Mouse embryos lacking CXXC finger protein 1 die prior to gastrulation, and embryonic stem cells lacking CXXC finger protein 1 are viable but are unable to achieve cellular differentiation and lineage commitment. CXXC finger protein 1 is a regulator of both cytosine and histone methylation. It physically interacts with DNA methyltransferase 1 and facilitates maintenance cytosine methylation. Rescue studies reveal that CXXC finger protein 1 contains redundant functional domains that are sufficient to support cellular differentiation and proper levels of cytosine methylation. CXXC finger protein 1 is also a component of the Setd1 histone H3-Lys4 methyltransferase complexes and functions to target these enzymes to unmethylated CpG islands. Depletion of CXXC finger protein 1 leads to loss of histone H3-Lys4 tri-methylation at CpG islands and inappropriate drifting of this euchromatin mark into areas of hetero-chromatin. Thus, one function of CXXC finger protein 1 is to serve as an effector protein that interprets cytosine methylation patterns and facilitates crosstalk with histone-modifying enzymes.

ExprEssence--revealing the essence of differential experimental data in the context of an interaction/regulation net-work

Background

Experimentalists are overwhelmed by high-throughput data and there is an urgent need to condense information into simple hypotheses. For example, large amounts of microarray and deep sequencing data are becoming available, describing a variety of experimental conditions such as gene knockout and knockdown, the effect of interventions, and the differences between tissues and cell lines.

Results

To address this challenge, we developed a method, implemented as a Cytoscape plugin called ExprEssence. As input we take a network of interaction, stimulation and/or inhibition links between genes/proteins, and differential data, such as gene expression data, tracking an intervention or development in time. We condense the network, highlighting those links across which the largest changes can be observed. Highlighting is based on a simple formula inspired by the law of mass action. We can interactively modify the threshold for highlighting and instantaneously visualize results. We applied ExprEssence to three scenarios describing kidney podocyte biology, pluripotency and ageing: 1) We identify putative processes involved in podocyte (de-)differentiation and validate one prediction experimentally. 2) We predict and validate the expression level of a transcription factor involved in pluripotency. 3) Finally, we generate plausible hypotheses on the role of apoptosis, cell cycle deregulation and DNA repair in ageing data obtained from the hippocampus.

Conclusion

Reducing the size of gene/protein networks to the few links affected by large changes allows to screen for putative mechanistic relationships among the genes/proteins that are involved in adaptation to different experimental conditions, yielding important hypotheses, insights and suggestions for new experiments. We note that we do not focus on the identification of 'active subnetworks'. Instead we focus on the identification of single links (which may or may not form subnetworks), and these single links are much easier to validate experimentally than submodules. ExprEssence is available at http://sourceforge.net/projects/expressence/.

Structure and function of mammalian DNA methyltransferases

DNA methylation plays an important role in epigenetic signalling, having an impact on gene regulation, chromatin structure, development and disease. Here, we review the structures and functions of the mammalian DNA methyltransferases Dnmt1, Dnmt3a and Dnmt3b, including their domain structures, catalytic mechanisms, localisation, regulation, post-translational modifications and interaction with chromatin and other proteins, summarising data obtained in genetic, cell biology and enzymatic studies. We focus on the question of how the molecular and enzymatic properties of these enzymes are connected to the dynamics of DNA methylation patterns and to the roles the enzymes play in the processes of de novo and maintenance DNA methylation. Recent enzymatic and genome-wide methylome data have led to a new model of genomic DNA methylation patterns based on the preservation of average levels of DNA methylation in certain regions, rather than the methylation states of individual CG sites.

How genome-wide approaches can be used to unravel the remaining secrets of the imprintome

Genomic imprinting is the differential expression of genes according to their transmitting parent and is achieved by labelling of the two alleles with different epigenetic marks. The majority of described imprinted genes are present in clusters with coordinate regulation. Multiple mechanisms are known to regulate this differential expression, including repression of one allele by the action of cis-acting macro non-coding RNAs, insulator elements, allele specific histone modifications and DNA methylation. A hallmark of all imprinted regions described so far is the presence of one or more differentially methylated regions (DMRs). A DMR is a nucleotide sequence rich in CpG dinucleotides that is specifically methylated on one parental chromosome and unmethylated on the allele derived from the other parent. This parent-specific differential methylation may be imparted during spermatogenesis or oogenesis (as is the case for gametic DMRs) or may be acquired during embryogenesis (somatic DMRs). This review will describe the advantages and disadvantages of some of the techniques that can be used to compare epigenetic marks between parental chromosomes and to understand how these marks affect the 3D interactions and monoallelic expression at imprinted loci. Recent advances in sequencing technologies, in particular, provide exciting new opportunities to study imprinting. These analyses are likely to lead to the full characterization of the 'imprintome', which includes uncovering the totality of imprinted genes within a genome, their epigenetic landscape and unique features that render them resistant to epigenetic reprogramming in the early embryo.

DNA methylation in cell differentiation and reprogramming: an emerging systematic view

Embryonic stem cells have the unique ability to indefinitely self-renew and differentiate into any cell type found in the adult body. Differentiated cells can, in turn, be reprogrammed to embryonic stem-like induced pluripotent stem cells, providing exciting opportunities for achieving patient-specific stem cell therapy while circumventing immunological obstacles and ethical controversies. Since both differentiation and reprogramming are governed by major changes in the epigenome, current directions in the field aim to uncover the epigenetic signals that give pluripotent cells their unique properties. DNA methylation is one of the major epigenetic factors that regulates gene expression in mammals and is essential for establishing cellular identity. Recent analyses of pluripotent and somatic cell methylomes have provided important insights into the extensive role of DNA methylation during cell-fate commitment and reprogramming. In this article, the recent progress of differentiation and reprogramming research illuminated by high-throughput studies is discussed in the context of DNA methylation.

The DNA methylome of human peripheral blood mononuclear cells

Analysis across the genome of patterns of DNA methylation reveals a rich landscape of allele-specific epigenetic modification and consequent effects on allele-specific gene expression.

DNA methylation plays an important role in biological processes in human health and disease. Recent technological advances allow unbiased whole-genome DNA methylation (methylome) analysis to be carried out on human cells. Using whole-genome bisulfite sequencing at 24.7-fold coverage (12.3-fold per strand), we report a comprehensive (92.62%) methylome and analysis of the unique sequences in human peripheral blood mononuclear cells (PBMC) from the same Asian individual whose genome was deciphered in the YH project. PBMC constitute an important source for clinical blood tests world-wide. We found that 68.4% of CpG sites and <0.2% of non-CpG sites were methylated, demonstrating that non-CpG cytosine methylation is minor in human PBMC. Analysis of the PBMC methylome revealed a rich epigenomic landscape for 20 distinct genomic features, including regulatory, protein-coding, non-coding, RNA-coding, and repeat sequences. Integration of our methylome data with the YH genome sequence enabled a first comprehensive assessment of allele-specific methylation (ASM) between the two haploid methylomes of any individual and allowed the identification of 599 haploid differentially methylated regions (hDMRs) covering 287 genes. Of these, 76 genes had hDMRs within 2 kb of their transcriptional start sites of which >80% displayed allele-specific expression (ASE). These data demonstrate that ASM is a recurrent phenomenon and is highly correlated with ASE in human PBMCs. Together with recently reported similar studies, our study provides a comprehensive resource for future epigenomic research and confirms new sequencing technology as a paradigm for large-scale epigenomics studies.

Epigenetic modifications such as addition of methyl groups to cytosine in DNA play a role in regulating gene expression. To better understand these processes, knowledge of the methylation status of all cytosine bases in the genome (the methylome) is required. DNA methylation can differ between the two gene copies (alleles) in each cell. Such allele-specific methylation (ASM) can be due to parental origin of the alleles (imprinting), X chromosome inactivation in females, and other as yet unknown mechanisms. This may significantly alter the expression profile arising from different allele combinations in different individuals. Using advanced sequencing technology, we have determined the methylome of human peripheral blood mononuclear cells (PBMC). Importantly, the PBMC were obtained from the same male Han Chinese individual whose complete genome had previously been determined. This allowed us, for the first time, to study genome-wide differences in ASM. Our analysis shows that ASM in PBMC is higher than can be accounted for by regions known to undergo parent-of-origin imprinting and frequently (>80%) correlates with allele-specific expression (ASE) of the corresponding gene. In addition, our data reveal a rich landscape of epigenomic variation for 20 genomic features, including regulatory, coding, and non-coding sequences, and provide a valuable resource for future studies. Our work further establishes whole-genome sequencing as an efficient method for methylome analysis.