Epigenetic research aims to understand heritable gene regulation that is not directly encoded in the DNA sequence. Epigenetic mechanisms such as DNA methylation and histone modifications modulate the packaging of the DNA in the nucleus and thereby influence gene expression. Patterns of epigenetic information are faithfully propagated over multiple cell divisions, which makes epigenetic regulation a key mechanism for cellular differentiation and cell fate decisions. In addition, incomplete erasure of epigenetic information can lead to complex patterns of non-Mendelian inheritance. Stochastic and environment-induced epigenetic defects are known to play a major role in cancer and ageing, and they may also contribute to mental disorders and autoimmune diseases. Recent technical advances such as ChIP-on-chip and ChIP-seq have started to convert epigenetic research into a high-throughput endeavor, to which bioinformatics is expected to make significant contributions. Here, we review pioneering computational studies that have contributed to epigenetic research. In addition, we give a brief introduction into epigenetics—targeted at bioinformaticians who are new to the field—and we outline future challenges in computational epigenetics.

A unique feature of the germ cell lineage is the generation of totipotency. A critical event in this context is DNA demethylation and the erasure of parental imprints in mouse primordial germ cells (PGCs) on embryonic day 11.5 (E11.5) after they enter into the developing gonads1, 2. Little is yet known about the mechanism involved, except that it is apparently an active process. We have examined the associated changes in the chromatin to gain further insights into this reprogramming event. Here we show that the chromatin changes occur in two steps. The first changes in nascent PGCs at E8.5 establish a distinctive chromatin signature that is reminiscent of pluripotency. Next, when PGCs are residing in the gonads, major changes occur in nuclear architecture accompanied by an extensive erasure of several histone modifications and exchange of histone variants. Furthermore, the histone chaperones HIRA and NAP-1 (NAP111), which are implicated in histone exchange, accumulate in PGC nuclei undergoing reprogramming. We therefore suggest that the mechanism of histone replacement is critical for these chromatin rearrangements to occur. The marked chromatin changes are intimately linked with genome-wide DNA demethylation. On the basis of the timing of the observed events, we propose that if DNA demethylation entails a DNA repair-based mechanism, the evident histone replacement would represent a repair-induced response event rather than being a prerequisite.

The epigenetic profile of germ cells, which is defined by modifications of DNA and chromatin, changes dynamically during their development. Many of the changes are associated with the acquisition of the capacity to support post-fertilization development. Our knowledge of this aspect has greatly increased— for example, insights into how the re-establishment of parental imprints is regulated. In addition, an emerging theme from recent studies is that epigenetic modifiers have key roles in germ-cell development itself — for example, epigenetics contributes to the gene-expression programme that is required for germ-cell development, regulation of meiosis and genomic integrity. Understanding epigenetic regulation in germ cells has implications for reproductive engineering technologies and human health.

To gain insight into the function of DNA methylation at cis-regulatory regions and its impact on gene expression, we measured methylation, RNA polymerase occupancy and histone modifications at 16,000 promoters in primary human somatic and germline cells. We find CpG-poor promoters hypermethylated in somatic cells, which does not preclude their activity. This methylation is present in male gametes and results in evolutionary loss of CpG dinucleotides, as measured by divergence between humans and primates. In contrast, strong CpG island promoters are mostly unmethylated, even when inactive. Weak CpG island promoters are distinct, as they are preferential targets for de novo methylation in somatic cells. Notably, most germline-specific genes are methylated in somatic cells, suggesting additional functional selection. These results show that promoter sequence and gene function are major predictors of promoter methylation states. Moreover, we observe that inactive unmethylated CpG island promoters show elevated levels of dimethylation of Lys4 of histone H3, suggesting that this chromatin mark may protect DNA from methylation.

We describe the results of a genome-wide analysis of human cells that suggests that most protein-coding genes, including most genes thought to be transcriptionally inactive, experience transcription initiation. We found that nucleosomes with H3K4me3 and H3K9,14Ac modifications, together with RNA polymerase II, occupy the promoters of most protein-coding genes in human embryonic stem cells. Only a subset of these genes produce detectable full-length transcripts and are occupied by nucleosomes with H3K36me3 modifications, a hallmark of elongation. The other genes experience transcription initiation but show no evidence of elongation, suggesting that they are predominantly regulated at postinitiation steps. Genes encoding most developmental regulators fall into this group. Our results also identify a class of genes that are excluded from experiencing transcription initiation, at which mechanisms that prevent initiation must predominate. These observations extend to differentiated cells, suggesting that transcription initiation at most genes is a general phenomenon in human cells.

Although histones can form nucleosomes on virtually any genomic sequence, DNA sequences show considerable variability in their binding affinity. We have used DNA sequences of Saccharomyces cerevisiae whose nucleosome binding affinities have been experimentally determined (Yuan et al. 2005) to train a support vector machine to identify the nucleosome formation potential of any given sequence of DNA. The DNA sequences whose nucleosome formation potential are most accurately predicted are those that contain strong nucleosome forming or inhibiting signals and are found within nucleosome length stretches of genomic DNA with continuous nucleosome formation or inhibition signals. We have accurately predicted the experimentally determined nucleosome positions across a well-characterized promoter region of S. cerevisiae and identified strong periodicity within 199 center-aligned mononucleosomes studied recently (Segal et al. 2006) despite there being no periodicity information used to train the support vector machine. Our analysis suggests that only a subset of nucleosomes are likely to be positioned by intrinsic sequence signals. This observation is consistent with the available experimental data and is inconsistent with the proposal of a nucleosome positioning code. Finally, we show that intrinsic nucleosome positioning signals are both more inhibitory and more variable in promoter regions than in open reading frames in S. cerevisiae.

We have developed a method for mapping unmethylated sites in the human genome based on the resistance of TspRI-digested ends to ExoIII nuclease degradation. Digestion with TspRI and methylation-sensitive restriction endonuclease HpaII, followed by ExoIII and single-strand DNA nuclease allowed removal of DNA fragments containing unmethylated HpaII sites. We then used array comparative genomic hybridization (CGH) to map the sequences depleted by these procedures in human genomes derived from five human tissues, a primary breast tumor, and two breast tumor cell lines. Analysis of methylation patterns of the normal tissue genomes indicates that the hypomethylated sites are enriched in the 5' end of widely expressed genes, including promoter, first exon, and first intron. In contrast, genomes of the MCF-7 and MDA-MB-231 cell lines show extensive hypomethylation in the intragenic and intergenic regions whereas the primary tumor exhibits a pattern between those of the normal tissue and the cell lines. A striking characteristic of tumor cell lines is the presence of megabase-sized hypomethylated zones. These hypomethylated zones are associated with large genes, fragile sites, evolutionary breakpoints, chromosomal rearrangement breakpoints, tumor suppressor genes, and with regions containing tissue-specific gene clusters or with gene-poor regions containing novel tissue-specific genes. Correlation with microarray analysis shows that genes with a hypomethylated sequence 2 kb up- or downstream of the transcription start site are highly expressed, whereas genes with extensive intragenic and 3' untranslated region (UTR) hypomethylation are silenced. The method described herein can be used for large-scale screening of changes in the methylation pattern in the genome of interest.

In the universe of science, two worlds have recently collided—those of RNA and chromatin. The intersection of these two fields has been impending, but evidence for such a meaningful collision has only recently become apparent. In this review, we discuss the implications for noncoding RNAs and the formation of specialized chromatin domains in various epigenetic processes as diverse as dosage compensation, RNA interference-mediated heterochromatin assembly and gene silencing, and programmed DNA elimination. While mechanistic details as to how the RNA and chromatin worlds connect remain unclear, intriguing parallels exist in the overall design and machinery used in model organisms from all eukaryotic kingdoms. The role of potential RNA-binding chromatin-associated proteins will be discussed as one possible link between RNA and chromatin.