Understanding the transcriptional regulation of pluripotent cells is of fundamental interest and will greatly inform efforts aimed at directing differentiation of embryonic stem (ES) cells or reprogramming somatic cells. We first analyzed the transcriptional profiles of mouse ES cells and primordial germ cells and identified genes upregulated in pluripotent cells both in vitro and in vivo. These genes are enriched for roles in transcription, chromatin remodeling, cell cycle, and DNA repair. We developed a novel computational algorithm, CompMoby, which combines analyses of sequences both aligned and non-aligned between different genomes with a probabilistic segmentation model to systematically predict short DNA motifs that regulate gene expression. CompMoby was used to identify conserved overrepresented motifs in genes upregulated in pluripotent cells. We show that the motifs are preferentially active in undifferentiated mouse ES and embryonic germ cells in a sequence-specific manner, and that they can act as enhancers in the context of an endogenous promoter. Importantly, the activity of the motifs is conserved in human ES cells. We further show that the transcription factor NF-Y specifically binds to one of the motifs, is differentially expressed during ES cell differentiation, and is required for ES cell proliferation. This study provides novel insights into the transcriptional regulatory networks of pluripotent cells. Our results suggest that this systematic approach can be broadly applied to understanding transcriptional networks in mammalian species.

Dimethylation of histone H3 lysine 9 (H3K9m2) and trimethylation of histone H3 lysine 27 (H3K27m3) are two hallmarks of transcriptional repression in many organisms. In Arabidopsis thaliana, H3K27m3 is targeted by Polycomb Group (PcG) proteins and is associated with silent protein-coding genes, while H3K9m2 is correlated with DNA methylation and is associated with transposons and repetitive sequences. Recently, ectopic genic DNA methylation in the CHG context (where H is any base except G) has been observed in globally DNA hypomethylated mutants such as met1, but neither the nature of the hypermethylated loci nor the biological significance of this epigenetic phenomenon have been investigated. Here, we generated high-resolution, genome-wide maps of both H3K9m2 and H3K27m3 in wild-type and met1 plants, which we integrated with transcriptional data, to explore the relationships between these two marks. We found that ectopic H3K9m2 observed in met1 can be due to defects in IBM1-mediated H3K9m2 demethylation at some sites, but most importantly targets H3K27m3-marked genes, suggesting an interplay between these two silencing marks. Furthermore, H3K9m2/DNA-hypermethylation at these PcG targets in met1 is coupled with a decrease in H3K27m3 marks, whereas CG/H3K9m2 hypomethylated transposons become ectopically H3K27m3 hypermethylated. Our results bear interesting similarities with cancer cells, which show global losses of DNA methylation but ectopic hypermethylation of genes previously marked by H3K27m3.

Gibbon species have accumulated an unusually high number of chromosomal changes since diverging from the common hominoid ancestor 15-18 million years ago. The cause of this increased rate of chromosomal rearrangements is not known, nor is it known if genome architecture has a role. To address this question, we analyzed sequences spanning 57 breaks of synteny between northern white-cheeked gibbons (Nomascus l. leucogenys) and humans. We find that the breakpoint regions are enriched in segmental duplications and repeats, with Alu elements being the most abundant. Alus located near the gibbon breakpoints (<150 bp) have a higher CpG content than other Alus. Bisulphite allelic sequencing reveals that these gibbon Alus have a lower average density of methylated cytosine that their human orthologues. The finding of higher CpG content and lower average CpG methylation suggests that the gibbon Alu elements are epigenetically distinct from their human orthologues. The association between undermethylation and chromosomal rearrangement in gibbons suggests a correlation between epigenetic state and structural genome variation in evolution.

The Drosophila MSL complex mediates dosage compensation by increasing transcription of the single X chromosome in males approximately two-fold. This is accomplished through recognition of the X chromosome and subsequent acetylation of histone H4K16 on X-linked genes. Initial binding to the X is thought to occur at "entry sites" that contain a consensus sequence motif ("MSL recognition element" or MRE). However, this motif is only ∼2 fold enriched on X, and only a fraction of the motifs on X are initially targeted. Here we ask whether chromatin context could distinguish between utilized and non-utilized copies of the motif, by comparing their relative enrichment for histone modifications and chromosomal proteins mapped in the modENCODE project. Through a comparative analysis of the chromatin features in male S2 cells (which contain MSL complex) and female Kc cells (which lack the complex), we find that the presence of active chromatin modifications, together with an elevated local GC content in the surrounding sequences, has strong predictive value for functional MSL entry sites, independent of MSL binding. We tested these sites for function in Kc cells by RNAi knockdown of Sxl, resulting in induction of MSL complex. We show that ectopic MSL expression in Kc cells leads to H4K16 acetylation around these sites and a relative increase in X chromosome transcription. Collectively, our results support a model in which a pre-existing active chromatin environment, coincident with H3K36me3, contributes to MSL entry site selection. The consequences of MSL targeting of the male X chromosome include increase in nucleosome lability, enrichment for H4K16 acetylation and JIL-1 kinase, and depletion of linker histone H1 on active X-linked genes. Our analysis can serve as a model for identifying chromatin and local sequence features that may contribute to selection of functional protein binding sites in the genome.

Facioscapulohumeral dystrophy (FSHD) is an autosomal dominant muscular dystrophy in which no mutation of pathogenic gene(s) has been identified. Instead, the disease is, in most cases, genetically linked to a contraction in the number of 3.3 kb D4Z4 repeats on chromosome 4q. How contraction of the 4qter D4Z4 repeats causes muscular dystrophy is not understood. In addition, a smaller group of FSHD cases are not associated with D4Z4 repeat contraction (termed "phenotypic" FSHD), and their etiology remains undefined. We carried out chromatin immunoprecipitation analysis using D4Z4-specific PCR primers to examine the D4Z4 chromatin structure in normal and patient cells as well as in small interfering RNA (siRNA)-treated cells. We found that SUV39H1-mediated H3K9 trimethylation at D4Z4 seen in normal cells is lost in FSHD. Furthermore, the loss of this histone modification occurs not only at the contracted 4q D4Z4 allele, but also at the genetically intact D4Z4 alleles on both chromosomes 4q and 10q, providing the first evidence that the genetic change (contraction) of one 4qD4Z4 allele spreads its effect to other genomic regions. Importantly, this epigenetic change was also observed in the phenotypic FSHD cases with no D4Z4 contraction, but not in other types of muscular dystrophies tested. We found that HP1gamma and cohesin are co-recruited to D4Z4 in an H3K9me3-dependent and cell type-specific manner, which is disrupted in FSHD. The results indicate that cohesin plays an active role in HP1 recruitment and is involved in cell type-specific D4Z4 chromatin regulation. Taken together, we identified the loss of both histone H3K9 trimethylation and HP1gamma/cohesin binding at D4Z4 to be a faithful marker for the FSHD phenotype. Based on these results, we propose a new model in which the epigenetic change initiated at 4q D4Z4 spreads its effect to other genomic regions, which compromises muscle-specific gene regulation leading to FSHD pathogenesis.

DNA methylation patterns are altered in numerous diseases and often correlate with clinically relevant information such as disease subtypes, prognosis and drug response. With suitable assays and after validation in large cohorts, such associations can be exploited for clinical diagnostics and personalized treatment decisions. Here we describe the results of a community-wide benchmarking study comparing the performance of all widely used methods for DNA methylation analysis that are compatible with routine clinical use. We shipped 32 reference samples to 18 laboratories in seven different countries. Researchers in those laboratories collectively contributed 21 locus-specific assays for an average of 27 predefined genomic regions, as well as six global assays. We evaluated assay sensitivity on low-input samples and assessed the assays' ability to discriminate between cell types. Good agreement was observed across all tested methods, with amplicon bisulfite sequencing and bisulfite pyrosequencing showing the best all-round performance. Our technology comparison can inform the selection, optimization and use of DNA methylation assays in large-scale validation studies, biomarker development and clinical diagnostics.