Hendrik Marks

The Transcriptional and Epigenomic Foundations of Ground State Pluripotency

Summary Mouse embryonic stem (ES) cells grown in serum exhibit greater heterogeneity in morphology and expression of pluripotency factors than ES cells cultured in defined medium with inhibitors of two kinases (Mek and GSK3), a condition known as “2i” postulated to establish a naive ground state. We show that the transcriptome and epigenome profiles of serum- and 2i-grown ES cells are distinct. 2i-treated cells exhibit lower expression of lineage-affiliated genes, reduced prevalence at promoters of the repressive histone modification H3K27me3, and fewer bivalent domains, which are thought to mark genes poised for either up- or downregulation. Nonetheless, serum- and 2i-grown ES cells have similar differentiation potential. Precocious transcription of developmental genes in 2i is restrained by RNA polymerase II promoter-proximal pausing. These findings suggest that transcriptional potentiation and a permissive chromatin context characterize the ground state and that exit from it may not require a metastable intermediate or multilineage priming.

Sequential ChIP-bisulfite sequencing enables direct genome-scale investigation of chromatin and DNA methylation cross-talk

Cross-talk between DNA methylation and histone modifications drives the establishment of composite epigenetic signatures and is traditionally studied using correlative rather than direct approaches. Here, we present sequential ChIP-bisulfite-sequencing (ChIP-BS-seq) as an approach to quantitatively assess DNA methylation patterns associated with chromatin modifications or chromatin-associated factors directly. A chromatin-immunoprecipitation (ChIP)-capturing step is used to obtain a restricted representation of the genome occupied by the epigenetic feature of interest, for which a single-base resolution DNA methylation map is then generated. When applied to H3 lysine 27 trimethylation (H3K27me3), we found that H3K27me3 and DNA methylation are compatible throughout most of the genome, except for CpG islands, where these two marks are mutually exclusive. Further ChIP-BS-seq-based analysis in Dnmt triple-knockout (TKO) embryonic stem cells revealed that total loss of CpG methylation is associated with alteration of H3K27me3 levels throughout the genome: H3K27me3 in localized peaks is decreased while broad local enrichments (BLOCs) of H3K27me3 are formed. At an even broader scale, these BLOCs correspond to regions of high DNA methylation in wild-type ES cells, suggesting that DNA methylation prevents H3K27me3 deposition locally and at a megabase scale. Our strategy provides a unique way of investigating global interdependencies between DNA methylation and other chromatin features.

Whole-Genome Bisulfite Sequencing of Two Distinct Interconvertible DNA Methylomes of Mouse Embryonic Stem Cells

The use of two kinase inhibitors (2i) enables derivation of mouse embryonic stem cells (ESCs) in the pluripotent ground state. Using whole-genome bisulfite sequencing (WGBS), we show that male 2i ESCs are globally hypomethylated compared to conventional ESCs maintained in serum. In serum, female ESCs are hypomethyated similarly to male ESCs in 2i, and DNA methylation is further reduced in 2i. Regions with elevated DNA methylation in 2i strongly correlate with the presence of H3K9me3 on endogenous retroviruses (ERVs) and imprinted loci. The methylome of male ESCs in serum parallels postimplantation blastocyst cells, while 2i stalls ESCs in a hypomethylated, ICM-like state. WGBS analysis during adaptation of 2i ESCs to serum suggests that deposition of DNA methylation is largely random, while loss of DNA methylation during reversion to 2i occurs passively, initiating at TET1 binding sites. Together, our analysis provides insight into DNA methylation dynamics in cultured ESCs paralleling early developmental processes.

A role for cohesin in T cell receptor rearrangement and thymocyte differentiation

Cohesin enables post-replicative DNA repair and chromosome segregation by holding sister chromatids together from the time of DNA replication in S phase until mitosis. There is growing evidence that cohesin also forms long-range chromosomal cis-interactions and may regulate gene expression in association with CTCF, mediator or tissue-specific transcription factors. Human cohesinopathies such as Cornelia de Lange syndrome are thought to result from impaired non-canonical cohesin functions, but a clear distinction between the cell-division-related and cell-division-independent functions of cohesion—as exemplified in Drosophila—has not been demonstrated in vertebrate systems. To address this, here we deleted the cohesin locus Rad21 in mouse thymocytes at a time in development when these cells stop cycling and rearrange their T-cell receptor (TCR) α locus (Tcra). Rad21-deficient thymocytes had a normal lifespan and retained the ability to differentiate, albeit with reduced efficiency. Loss of Rad21 led to defective chromatin architecture at the Tcra locus, where cohesion-binding sites flank the TEA promoter and the Eα enhancer, and demarcate Tcra from interspersed Tcrd elements and neighbouring housekeeping genes. Cohesin was required for long-range promoter–enhancer interactions, Tcra transcription, H3K4me3 histone modifications that recruit the recombination machinery and Tcra rearrangement. Provision of pre-rearranged TCR transgenes largely rescued thymocyte differentiation, demonstrating that among thousands of potential target genes across the genome, defective Tcra rearrangement was limiting for the differentiation of cohesin-deficient thymocytes. These findings firmly establish a cell-division-independent role for cohesin in Tcra locus rearrangement and provide a comprehensive account of the mechanisms by which cohesin enables cellular differentiation in a well-characterized mammalian system.

High-resolution analysis of epigenetic changes associated with X inactivation

Differentiation of female murine ES cells triggers silencing of one X chromosome through X-chromosome inactivation (XCI). Immunofluorescence studies showed that soon after Xist RNA coating the inactive X (Xi) undergoes many heterochromatic changes, including the acquisition of H3K27me3. However, the mechanisms that lead to the establishment of heterochromatin remain unclear. We first analyze chromatin changes by ChIP-chip, as well as RNA expression, around the X-inactivation center (Xic) in female and male ES cells, and their day 4 and 10 differentiated derivatives. A dynamic epigenetic landscape is observed within the Xic locus. Tsix repression is accompanied by deposition of H3K27me3 at its promoter during differentiation of both female and male cells. However, only in female cells does an active epigenetic landscape emerge at the Xist locus, concomitant with high Xist expression. Several regions within and around the Xic show unsuspected chromatin changes, and we define a series of unusual loci containing highly enriched H3K27me3. Genome-wide ChIP-seq analyses show a female-specific quantitative increase of H3K27me3 across the X chromosome as XCI proceeds in differentiating female ES cells. Using female ES cells with nonrandom XCI and polymorphic X chromosomes, we demonstrate that this increase is specific to the Xi by allele-specific SNP mapping of the ChIP-seq tags. H3K27me3 becomes evenly associated with the Xi in a chromosome-wide fashion. A selective and robust increase of H3K27me3 and concomitant decrease in H3K4me3 is observed over active genes. This indicates that deposition of H3K27me3 during XCI is tightly associated with the act of silencing of individual genes across the Xi.

Epigenetic regulation of learning and memory by Drosophila EHMT/G9a.

Behavioral phenotyping and genome-wide profiling of the histone modifier EHMT in Drosophila reveals a mechanism through which an epigenetic writer may control cognition.

Mll2 is required for H3K4 trimethylation on bivalent promoters in embryonic stem cells, whereas Mll1 is redundant

Trimethylation of histone H3 lysine 4 (H3K4me3) at the promoters of actively transcribed genes is a universal epigenetic mark and a key product of Trithorax group action. Here, we show that Mll2, one of the six Set1/Trithorax-type H3K4 methyltransferases in mammals, is required for trimethylation of bivalent promoters in mouse embryonic stem cells. Mll2 is bound to bivalent promoters but also to most active promoters, which do not require Mll2 for H3K4me3 or mRNA expression. By contrast, the Set1 complex (Set1C) subunit Cxxc1 is primarily bound to active but not bivalent promoters. This indicates that bivalent promoters rely on Mll2 for H3K4me3 whereas active promoters have more than one bound H3K4 methyltransferase, including Set1C. Removal of Mll1, sister to Mll2, had almost no effect on any promoter unless Mll2 was also removed, indicating functional backup between these enzymes. Except for a subset, loss of H3K4me3 on bivalent promoters did not prevent responsiveness to retinoic acid, thereby arguing against a priming model for bivalency. In contrast, we propose that Mll2 is the pioneer trimethyltransferase for promoter definition in the naïve epigenome and that Polycomb group action on bivalent promoters blocks the premature establishment of active, Set1C-bound, promoters.

Genetic variation among isolates of White spot syndrome virus

Summary.White spot syndrome virus (WSSV), member of a new virus family called Nimaviridae, is a major scourge in worldwide shrimp cultivation. Geographical isolates of WSSV identified so far are very similar in morphology and proteome, and show little difference in restriction fragment length polymorphism (RFLP) pattern. We have mapped the genomic differences between three completely sequenced WSSV isolates, originating from Thailand (WSSV-TH), China (WSSV-CN) and Taiwan (WSSV-TW). Alignment of the genomic sequences of these geographical isolates revealed an overall nucleotide identity of 99.32%. The major difference among the three isolates is a deletion of approximately 13 kb (WSSV-TH) and 1 kb (WSSV-CN), present in the same genomic region, relative to WSSV-TW. A second difference involves a genetically variable region of about 750 bp. All other variations >2 bp between the three isolates are located in repeat regions along the genome. Except for the homologous regions (hr1, hr3, hr8 and hr9), these variable repeat regions are almost exclusively located in ORFs, of which the genomic repeat regions in ORF75, ORF94 and ORF125 can be used for PCR based classification of WSSV isolates in epidemiological studies. Furthermore, the comparison identified highly invariable genomic loci, which may be used for reliable monitoring of WSSV infections and for shrimp health certification.

An RNA-Seq Strategy to Detect the Complete Coding and Non-Coding Transcriptome Including Full-Length Imprinted Macro ncRNAs

Imprinted macro non-protein-coding (nc) RNAs are cis-repressor transcripts that silence multiple genes in at least three imprinted gene clusters in the mouse genome. Similar macro or long ncRNAs are abundant in the mammalian genome. Here we present the full coding and non-coding transcriptome of two mouse tissues: differentiated ES cells and fetal head using an optimized RNA-Seq strategy. The data produced is highly reproducible in different sequencing locations and is able to detect the full length of imprinted macro ncRNAs such as Airn and Kcnq1ot1, whose length ranges between 80–118 kb. Transcripts show a more uniform read coverage when RNA is fragmented with RNA hydrolysis compared with cDNA fragmentation by shearing. Irrespective of the fragmentation method, all coding and non-coding transcripts longer than 8 kb show a gradual loss of sequencing tags towards the 3′ end. Comparisons to published RNA-Seq datasets show that the strategy presented here is more efficient in detecting known functional imprinted macro ncRNAs and also indicate that standardization of RNA preparation protocols would increase the comparability of the transcriptome between different RNA-Seq datasets.

Otx2 and Oct4 drive early enhancer activation during embryonic stem cell transition from naive pluripotency.

Summary Embryonic stem cells (ESCs) are unique in that they have the capacity to differentiate into all of the cell types in the body. We know a lot about the complex transcriptional control circuits that maintain the naive pluripotent state under self-renewing conditions but comparatively less about how cells exit from this state in response to differentiation stimuli. Here, we examined the role of Otx2 in this process in mouse ESCs and demonstrate that it plays a leading role in remodeling the gene regulatory networks as cells exit from ground state pluripotency. Otx2 drives enhancer activation through affecting chromatin marks and the activity of associated genes. Mechanistically, Oct4 is required for Otx2 expression, and reciprocally, Otx2 is required for efficient Oct4 recruitment to many enhancer regions. Therefore, the Oct4-Otx2 regulatory axis actively establishes a new regulatory chromatin landscape during the early events that accompany exit from ground state pluripotency.