The impact of the embryonic DNA methylation program on CTCF-mediated genome regulation.

During mammalian embryogenesis, both the 5-cytosine DNA methylation (5meC) landscape and three dimensional (3D) chromatin architecture are profoundly remodeled during a process known as 'epigenetic reprogramming.' An understudied aspect of epigenetic reprogramming is how the 5meC flux, per se, affects the 3D genome. This is pertinent given the 5meC-sensitivity of DNA binding for a key regulator of chromosome folding: CTCF. We profiled the CTCF binding landscape using a mouse embryonic stem cell (ESC) differentiation protocol that models embryonic 5meC dynamics. Mouse ESCs lacking DNA methylation machinery are able to exit naive pluripotency, thus allowing for dissection of subtle effects of CTCF on gene expression. We performed CTCF HiChIP in both wild-type and mutant conditions to assess gained CTCF-CTCF contacts in the absence of 5meC. We performed H3K27ac HiChIP to determine the impact that ectopic CTCF binding has on cis-regulatory contacts. Using 5meC epigenome editing, we demonstrated that the methyl-mark is able to impair CTCF binding at select loci. Finally, a detailed dissection of the imprinted Zdbf2 locus showed how 5meC-antagonism of CTCF allows for proper gene regulation during differentiation. This work provides a comprehensive overview of how 5meC impacts the 3D genome in a relevant model for early embryonic events.

A genetic screen identifies BEND3 as a regulator of bivalent gene expression and global DNA methylation.

Epigenetic mechanisms are essential to establish and safeguard cellular identities in mammals. They dynamically regulate the expression of genes, transposable elements and higher-order chromatin structures. Consequently, these chromatin marks are indispensable for mammalian development and alterations often lead to disease, such as cancer. Bivalent promoters are especially important during differentiation and development. Here we used a genetic screen to identify new regulators of a bivalent repressed gene. We identify BEND3 as a regulator of hundreds of bivalent promoters, some of which it represses, and some of which it activates. We show that BEND3 is recruited to a CpG-containg consensus site that is present in multiple copies in many bivalent promoters. Besides having direct effect on the promoters it binds, the loss of BEND3 leads to genome-wide gains of DNA methylation, which are especially marked at regions normally protected by the TET enzymes. DNA hydroxymethylation is reduced in Bend3 mutant cells, possibly as consequence of altered gene expression leading to diminished alpha-ketoglutarate production, thus lowering TET activity. Our results clarify the direct and indirect roles of an important chromatin regulator, BEND3, and, more broadly, they shed light on the regulation of bivalent promoters.

Divergent transcriptional and transforming properties of PAX3-FOXO1 and PAX7-FOXO1 paralogs.

The hallmarks of the alveolar subclass of rhabdomyosarcoma are chromosomal translocations that generate chimeric PAX3-FOXO1 or PAX7-FOXO1 transcription factors. Overexpression of either PAX-FOXO1s results in related cell transformation in animal models. Yet, in patients the two structural genetic aberrations they derived from are associated with distinct pathological manifestations. To assess the mechanisms underlying these differences, we generated isogenic fibroblast lines expressing either PAX-FOXO1 paralog. Mapping of their genomic recruitment using CUT&Tag revealed that the two chimeric proteins have distinct DNA binding preferences. In addition, PAX7-FOXO1 binding results in greater recruitment of the H3K27ac activation mark than PAX3-FOXO1 binding and is accompanied by greater transcriptional activation of neighbouring genes. These effects are associated with a PAX-FOXO1-specific alteration in the expression of genes regulating cell shape and the cell cycle. Consistently, PAX3-FOXO1 accentuates fibroblast cellular traits associated with contractility and surface adhesion and limits entry into S phase. In contrast, PAX7-FOXO1 drives cells to adopt an amoeboid shape, reduces entry into M phase, and causes increased DNA damage. Altogether, our results argue that the diversity of rhabdomyosarcoma manifestation arises, in part, from the divergence between the genomic occupancy and transcriptional activity of PAX3-FOXO1 and PAX7-FOXO1.

The p38/HOG stress-activated protein kinase network couples growth to division in Candida albicans.

Cell size is a complex trait that responds to developmental and environmental cues. Quantitative size analysis of mutant strain collections disrupted for protein kinases and transcriptional regulators in the pathogenic yeast Candida albicans uncovered 66 genes that altered cell size, few of which overlapped with known size genes in the budding yeast Saccharomyces cerevisiae. A potent size regulator specific to C. albicans was the conserved p38/HOG MAPK module that mediates the osmostress response. Basal HOG activity inhibited the SBF G1/S transcription factor complex in a stress-independent fashion to delay the G1/S transition. The HOG network also governed ribosome biogenesis through the master transcriptional regulator Sfp1. Hog1 bound to the promoters and cognate transcription factors for ribosome biogenesis regulons and interacted genetically with the SBF G1/S machinery, and thereby directly linked cell growth and division. These results illuminate the evolutionary plasticity of size control and identify the HOG module as a nexus of cell cycle and growth regulation.

Development and application of an integrated allele-specific pipeline for methylomic and epigenomic analysis (MEA).

BACKGROUND: Allele-specific transcriptional regulation, including of imprinted genes, is essential for normal mammalian development. While the regulatory regions controlling imprinted genes are associated with DNA methylation (DNAme) and specific histone modifications, the interplay between transcription and these epigenetic marks at allelic resolution is typically not investigated genome-wide due to a lack of bioinformatic packages that can process and integrate multiple epigenomic datasets with allelic resolution. In addition, existing ad-hoc software only consider SNVs for allele-specific read discovery. This limitation omits potentially informative INDELs, which constitute about one fifth of the number of SNVs in mice, and introduces a systematic reference bias in allele-specific analyses. RESULTS: Here, we describe MEA, an INDEL-aware Methylomic and Epigenomic Allele-specific analysis pipeline which enables user-friendly data exploration, visualization and interpretation of allelic imbalance. Applying MEA to mouse embryonic datasets yields robust allele-specific DNAme maps and low reference bias. We validate allele-specific DNAme at known differentially methylated regions and show that automated integration of such methylation data with RNA- and ChIP-seq datasets yields an intuitive, multidimensional view of allelic gene regulation. MEA uncovers numerous novel dynamically methylated loci, highlighting the sensitivity of our pipeline. Furthermore, processing and visualization of epigenomic datasets from human brain reveals the expected allele-specific enrichment of H3K27ac and DNAme at imprinted as well as novel monoallelically expressed genes, highlighting MEA's utility for integrating human datasets of distinct provenance for genome-wide analysis of allelic phenomena. CONCLUSIONS: Our novel pipeline for standardized allele-specific processing and visualization of disparate epigenomic and methylomic datasets enables rapid analysis and navigation with allelic resolution. MEA is freely available as a Docker container at https://github.com/julienrichardalbert/MEA .

A functional portrait of Med7 and the mediator complex in Candida albicans.

Mediator is a multi-subunit protein complex that regulates gene expression in eukaryotes by integrating physiological and developmental signals and transmitting them to the general RNA polymerase II machinery. We examined, in the fungal pathogen Candida albicans, a set of conditional alleles of genes encoding Mediator subunits of the head, middle, and tail modules that were found to be essential in the related ascomycete Saccharomyces cerevisiae. Intriguingly, while the Med4, 8, 10, 11, 14, 17, 21 and 22 subunits were essential in both fungi, the structurally highly conserved Med7 subunit was apparently non-essential in C. albicans. While loss of CaMed7 did not lead to loss of viability under normal growth conditions, it dramatically influenced the pathogen's ability to grow in different carbon sources, to form hyphae and biofilms, and to colonize the gastrointestinal tracts of mice. We used epitope tagging and location profiling of the Med7 subunit to examine the distribution of the DNA sites bound by Mediator during growth in either the yeast or the hyphal form, two distinct morphologies characterized by different transcription profiles. We observed a core set of 200 genes bound by Med7 under both conditions; this core set is expanded moderately during yeast growth, but is expanded considerably during hyphal growth, supporting the idea that Mediator binding correlates with changes in transcriptional activity and that this binding is condition specific. Med7 bound not only in the promoter regions of active genes but also within coding regions and at the 3' ends of genes. By combining genome-wide location profiling, expression analyses and phenotyping, we have identified different Med7p-influenced regulons including genes related to glycolysis and the Filamentous Growth Regulator family. In the absence of Med7, the ribosomal regulon is de-repressed, suggesting Med7 is involved in central aspects of growth control.

DNA methylation shapes the Polycomb landscape during the exit from naive pluripotency.

In mammals, 5-methylcytosine (5mC) and Polycomb repressive complex 2 (PRC2)-deposited histone 3 lysine 27 trimethylation (H3K27me3) are generally mutually exclusive at CpG-rich regions. As mouse embryonic stem cells exit the naive pluripotent state, there is massive gain of 5mC concomitantly with restriction of broad H3K27me3 to 5mC-free, CpG-rich regions. To formally assess how 5mC shapes the H3K27me3 landscape, we profiled the epigenome of naive and differentiated cells in the presence and absence of the DNA methylation machinery. Surprisingly, we found that 5mC accumulation is not required to restrict most H3K27me3 domains. Instead, this 5mC-independent H3K27me3 restriction is mediated by aberrant expression of the PRC2 antagonist Ezhip (encoding EZH inhibitory protein). At the subset of regions where 5mC appears to genuinely supplant H3K27me3, we identified 163 candidate genes that appeared to require 5mC deposition and/or H3K27me3 depletion for their activation in differentiated cells. Using site-directed epigenome editing to directly modulate 5mC levels, we demonstrated that 5mC deposition is sufficient to antagonize H3K27me3 deposition and confer gene activation at individual candidates. Altogether, we systematically measured the antagonistic interplay between 5mC and H3K27me3 in a system that recapitulates early embryonic dynamics. Our results suggest that H3K27me3 restraint depends on 5mC, both directly and indirectly. Our study also implies a noncanonical role of 5mC in gene activation, which may be important not only for normal development but also for cancer progression, as oncogenic cells frequently exhibit dynamic replacement of 5mC for H3K27me3 and vice versa.

Conservation and divergence of canonical and non-canonical imprinting in murids.

BACKGROUND: Genomic imprinting affects gene expression in a parent-of-origin manner and has a profound impact on complex traits including growth and behavior. While the rat is widely used to model human pathophysiology, few imprinted genes have been identified in this murid. To systematically identify imprinted genes and genomic imprints in the rat, we use low input methods for genome-wide analyses of gene expression and DNA methylation to profile embryonic and extraembryonic tissues at allele-specific resolution. RESULTS: We identify 14 and 26 imprinted genes in these tissues, respectively, with 10 of these genes imprinted in both tissues. Comparative analyses with mouse reveal that orthologous imprinted gene expression and associated canonical DNA methylation imprints are conserved in the embryo proper of the Muridae family. However, only 3 paternally expressed imprinted genes are conserved in the extraembryonic tissue of murids, all of which are associated with non-canonical H3K27me3 imprints. The discovery of 8 novel non-canonical imprinted genes unique to the rat is consistent with more rapid evolution of extraembryonic imprinting. Meta-analysis of novel imprinted genes reveals multiple mechanisms by which species-specific imprinted expression may be established, including H3K27me3 deposition in the oocyte, the appearance of ZFP57 binding motifs, and the insertion of endogenous retroviral promoters. CONCLUSIONS: In summary, we provide an expanded list of imprinted loci in the rat, reveal the extent of conservation of imprinted gene expression, and identify potential mechanisms responsible for the evolution of species-specific imprinting.

Maternal DNMT3A-dependent de novo methylation of the paternal genome inhibits gene expression in the early embryo.

De novo DNA methylation (DNAme) during mammalian spermatogenesis yields a densely methylated genome, with the exception of CpG islands (CGIs), which are hypomethylated in sperm. While the paternal genome undergoes widespread DNAme loss before the first S-phase following fertilization, recent mass spectrometry analysis revealed that the zygotic paternal genome is paradoxically also subject to a low level of de novo DNAme. However, the loci involved, and impact on transcription were not addressed. Here, we employ allele-specific analysis of whole-genome bisulphite sequencing data and show that a number of genomic regions, including several dozen CGI promoters, are de novo methylated on the paternal genome by the 2-cell stage. A subset of these promoters maintains DNAme through development to the blastocyst stage. Consistent with paternal DNAme acquisition, many of these loci are hypermethylated in androgenetic blastocysts but hypomethylated in parthenogenetic blastocysts. Paternal DNAme acquisition is lost following maternal deletion of Dnmt3a, with a subset of promoters, which are normally transcribed from the paternal allele in blastocysts, being prematurely transcribed at the 4-cell stage in maternal Dnmt3a knockout embryos. These observations uncover a role for maternal DNMT3A activity in post-fertilization epigenetic reprogramming and transcriptional silencing of the paternal genome.

LTR retrotransposons transcribed in oocytes drive species-specific and heritable changes in DNA methylation.

De novo DNA methylation (DNAme) during mouse oogenesis occurs within transcribed regions enriched for H3K36me3. As many oocyte transcripts originate in long terminal repeats (LTRs), which are heterogeneous even between closely related mammals, we examined whether species-specific LTR-initiated transcription units (LITs) shape the oocyte methylome. Here we identify thousands of syntenic regions in mouse, rat, and human that show divergent DNAme associated with private LITs, many of which initiate in lineage-specific LTR retrotransposons. Furthermore, CpG island (CGI) promoters methylated in mouse and/or rat, but not human oocytes, are embedded within rodent-specific LITs and vice versa. Notably, at a subset of such CGI promoters, DNAme persists on the maternal genome in fertilized and parthenogenetic mouse blastocysts or in human placenta, indicative of species-specific epigenetic inheritance. Polymorphic LITs are also responsible for disparate DNAme at promoter CGIs in distantly related mouse strains, revealing that LITs also promote intra-species divergence in CGI DNAme.