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.

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.

Transposable elements (TEs) comprise almost half of primate genomes and their aberrant regulation can result in deleterious effects. In pluripotent stem cells, rapidly evolving KRAB-ZNF genes target TEs for silencing by H3K9me3. To investigate the evolution of TE silencing, we performed H3K9me3 ChIP-seq experiments in induced pluripotent stem cells from 10 human and 7 chimpanzee individuals. We identified four million orthologous TEs and found the SVA and ERV families to be marked most frequently by H3K9me3. We found little evidence of inter-species differences in TE silencing, with as many as 82% of putatively silenced TEs marked at similar levels in humans and chimpanzees. TEs that are preferentially silenced in one species are a similar age to those silenced in both species and are not more likely to be associated with expression divergence of nearby orthologous genes. Our data suggest limited species-specificity of TE silencing across 6 million years of primate evolution.

RUNX1 is a transcription factor that plays key roles in hematopoietic development and in hematopoiesis and lymphopoiesis. In this article, we report that RUNX1 regulates a gene expression program in naive mouse B cells that affects the dynamics of cell cycle entry in response to stimulation of the BCR. Conditional knockout of Runx1 in mouse resting B cells resulted in accelerated entry into S-phase after BCR engagement. Our results indicate that Runx1 regulates the cyclin D2 (Ccnd2) gene, the immediate early genes Fosl2, Atf3, and Egr2, and the Notch pathway gene Rbpj in mouse B cells, reducing the rate at which transcription of these genes increases after BCR stimulation. RUNX1 interacts with the chromatin remodeler SNF-2-related CREB-binding protein activator protein (SRCAP), recruiting it to promoter and enhancer regions of the Ccnd2 gene. BCR-mediated activation triggers switching between binding of RUNX1 and its paralog RUNX3 and between SRCAP and the switch/SNF remodeling complex member BRG1. Binding of BRG1 is increased at the Ccnd2 and Rbpj promoters in the Runx1 knockout cells after BCR stimulation. We also find that RUNX1 exerts positive or negative effects on a number of genes that affect the activation response of mouse resting B cells. These include Cd22 and Bank1, which act as negative regulators of the BCR, and the IFN receptor subunit gene Ifnar1 The hyperresponsiveness of the Runx1 knockout B cells to BCR stimulation and its role in regulating genes that are associated with immune regulation suggest that RUNX1 could be involved in regulating B cell tolerance.

Cohesin is a regulator of genome architecture with roles in sister chromatid cohesion and chromosome compaction. The recruitment and mobility of cohesin complexes on DNA is restricted by nucleosomes. Here, we show that the role of cohesin in chromosome organization requires the histone chaperone FACT ('facilitates chromatin transcription') in Saccharomyces cerevisiae. We find that FACT interacts directly with cohesin, and is dynamically required for its localization on chromatin. Depletion of FACT in metaphase cells prevents cohesin accumulation at pericentric regions and causes reduced binding on chromosome arms. Using the Hi-C technique, we show that cohesin-dependent TAD (topological associated domain)-like structures in G1 and metaphase chromosomes are reduced in the absence of FACT. Sister chromatid cohesion is intact in FACT-depleted cells, although chromosome segregation failure is observed. Our data show that FACT contributes to the formation of cohesin-dependent TADs, thus uncovering a new role for this complex in nuclear organization during interphase and mitotic chromosome folding.

Innate immune responses rely on inducible gene expression programmes which, in contrast to steady-state transcription, are highly dependent on cohesin. Here we address transcriptional parameters underlying this cohesin-dependence by single-molecule RNA-FISH and single-cell RNA-sequencing. We show that inducible innate immune genes are regulated predominantly by an increase in the probability of active transcription, and that probabilities of enhancer and promoter transcription are coordinated. Cohesin has no major impact on the fraction of transcribed inducible enhancers, or the number of mature mRNAs produced per transcribing cell. Cohesin is, however, required for coupling the probabilities of enhancer and promoter transcription. Enhancer-promoter coupling may not be explained by spatial proximity alone, and at the model locus Il12b can be disrupted by selective inhibition of the cohesinopathy-associated BET bromodomain BD2. Our data identify discrete steps in enhancer-mediated inducible gene expression that differ in cohesin-dependence, and suggest that cohesin and BD2 may act on shared pathways.

DNA methylation is a heritable epigenetic modification involved in gene silencing, imprinting, and the suppression of retrotransposons. Global DNA demethylation occurs in the early embryo and the germ line, and may be mediated by Tet (ten eleven translocation) enzymes, which convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC). Tet enzymes have been studied extensively in mouse embryonic stem (ES) cells, which are generally cultured in the absence of vitamin C, a potential cofactor for Fe(II) 2-oxoglutarate dioxygenase enzymes such as Tet enzymes. Here we report that addition of vitamin C to mouse ES cells promotes Tet activity, leading to a rapid and global increase in 5hmC. This is followed by DNA demethylation of many gene promoters and upregulation of demethylated germline genes. Tet1 binding is enriched near the transcription start site of genes affected by vitamin C treatment. Importantly, vitamin C, but not other antioxidants, enhances the activity of recombinant Tet1 in a biochemical assay, and the vitamin-C-induced changes in 5hmC and 5mC are entirely suppressed in Tet1 and Tet2 double knockout ES cells. Vitamin C has a stronger effect on regions that gain methylation in cultured ES cells compared to blastocysts, and in vivo are methylated only after implantation. In contrast, imprinted regions and intracisternal A particle retroelements, which are resistant to demethylation in the early embryo, are resistant to vitamin-C-induced DNA demethylation. Collectively, the results of this study establish vitamin C as a direct regulator of Tet activity and DNA methylation fidelity in ES cells.

Several tools have been developed to enable biologists to perform initial browsing and exploration of sequencing data. However the computational tool set for further analyses often requires significant computational expertise to use and many of the biologists with the knowledge needed to interpret these data must rely on programming experts.

The development of cell therapy for repairing damaged or diseased skeletal muscle has been hindered by the inability to significantly expand immature, transplantable myogenic stem cells (MuSCs) in culture. To overcome this limitation, a deeper understanding of the mechanisms regulating the transition between activated, proliferating MuSCs and differentiation-primed, poorly engrafting progenitors is needed. Here, we show that methyltransferase Setd7 facilitates such transition by regulating the nuclear accumulation of β-catenin in proliferating MuSCs. Genetic or pharmacological inhibition of Setd7 promotes in vitro expansion of MuSCs and increases the yield of primary myogenic cell cultures. Upon transplantation, both mouse and human MuSCs expanded with a Setd7 small-molecule inhibitor are better able to repopulate the satellite cell niche, and treated mouse MuSCs show enhanced therapeutic potential in preclinical models of muscular dystrophy. Thus, Setd7 inhibition may help bypass a key obstacle in the translation of cell therapy for muscle disease.

The "arms race" relationship between transposable elements (TEs) and their host has promoted a series of epigenetic silencing mechanisms directed against TEs. Retrotransposons, a class of TEs, are often located in repressed regions and are thought to induce heterochromatin formation and spreading. However, direct evidence for TE-induced local heterochromatin in mammals is surprisingly scarce. To examine this phenomenon, we chose two mouse embryonic stem (ES) cell lines that possess insertionally polymorphic retrotransposons (IAP, ETn/MusD, and LINE elements) at specific loci in one cell line but not the other. Employing ChIP-seq data for these cell lines, we show that IAP elements robustly induce H3K9me3 and H4K20me3 marks in flanking genomic DNA. In contrast, such heterochromatin is not induced by LINE copies and only by a minority of polymorphic ETn/MusD copies. DNA methylation is independent of the presence of IAP copies, since it is present in flanking regions of both full and empty sites. Finally, such spreading into genes appears to be rare, since the transcriptional start sites of very few genes are less than one Kb from an IAP. However, the B3galtl gene is subject to transcriptional silencing via IAP-induced heterochromatin. Hence, although rare, IAP-induced local heterochromatin spreading into nearby genes may influence expression and, in turn, host fitness.

Retrotransposition of endogenous retroviruses (ERVs) poses a substantial threat to genome stability. Transcriptional silencing of a subset of these parasitic elements in early mouse embryonic and germ cell development is dependent upon the lysine methyltransferase SETDB1, which deposits H3K9 trimethylation (H3K9me3) and the co-repressor KAP1, which binds SETDB1 when SUMOylated. Here we identified the transcription co-factor hnRNP K as a novel binding partner of the SETDB1/KAP1 complex in mouse embryonic stem cells (mESCs) and show that hnRNP K is required for ERV silencing. RNAi-mediated knockdown of hnRNP K led to depletion of H3K9me3 at ERVs, concomitant with de-repression of proviral reporter constructs and specific ERV subfamilies, as well as a cohort of germline-specific genes directly targeted by SETDB1. While hnRNP K recruitment to ERVs is dependent upon KAP1, SETDB1 binding at these elements requires hnRNP K. Furthermore, an intact SUMO conjugation pathway is necessary for SETDB1 recruitment to proviral chromatin and depletion of hnRNP K resulted in reduced SUMOylation at ERVs. Taken together, these findings reveal a novel regulatory hierarchy governing SETDB1 recruitment and in turn, transcriptional silencing in mESCs.

Chromosome movements and programmed DNA double-strand breaks (DSBs) promote homologue pairing and initiate recombination at meiosis onset. Meiotic progression involves checkpoint-controlled termination of these events when all homologue pairs achieve synapsis and form crossover precursors. Exploiting the temporo-spatial organisation of the C. elegans germline and time-resolved methods of protein removal, we show that surveillance of the synaptonemal complex (SC) controls meiotic progression. In nuclei with fully synapsed homologues and crossover precursors, removing different meiosis-specific cohesin complexes, which are individually required for SC stability, or a SC central region component causes functional redeployment of the chromosome movement and DSB machinery, triggering whole-nucleus reorganisation. This apparent reversal of the meiotic programme requires CHK-2 kinase reactivation via signalling from chromosome axes containing HORMA proteins, but occurs in the absence of transcriptional changes. Our results uncover an unexpected plasticity of the meiotic programme and show how chromosome signalling orchestrates nuclear organisation and meiotic progression.

CD4 and CD8 mark helper and cytotoxic T cell lineages, respectively, and serve as coreceptors for MHC-restricted TCR recognition. How coreceptor expression is matched with TCR specificity is central to understanding CD4/CD8 lineage choice, but visualising coreceptor gene activity in individual selection intermediates has been technically challenging. It therefore remains unclear whether the sequence of coreceptor gene expression in selection intermediates follows a stereotypic pattern, or is responsive to signaling. Here we use single cell RNA sequencing (scRNA-seq) to classify mouse thymocyte selection intermediates by coreceptor gene expression. In the unperturbed thymus, Cd4+Cd8a- selection intermediates appear before Cd4-Cd8a+ selection intermediates, but the timing of these subsets is flexible according to the strength of TCR signals. Our data show that selection intermediates discriminate MHC class prior to the loss of coreceptor expression and suggest a model where signal strength informs the timing of coreceptor gene activity and ultimately CD4/CD8 lineage choice.

Remnants of ancient transposable elements (TEs) are abundant in mammalian genomes. These sequences contain multiple regulatory motifs and hence are capable of influencing expression of host genes. TEs are known to be released from epigenetic repression and can become transcriptionally active in cancer. Such activation could also lead to lineage-inappropriate activation of oncogenes, as previously described in lymphomas. However, there are few reports of this mechanism occurring in non-blood cancers. Here, we re-analyzed whole transcriptome data from a large cohort of patients with colon cancer, compared to matched normal colon control samples, to detect genes or transcripts ectopically expressed through activation of TE promoters. Among many such transcripts, we identified six where the affected gene has described role in cancer and where the TE-driven gene mRNA is expressed in primary colon cancer, but not normal matched tissue, and confirmed expression in colon cancer-derived cell lines. We further characterized a TE-gene chimeric transcript involving the Interleukin 33 (IL-33) gene (termed LTR-IL-33), that is ectopically expressed in a subset of colon cancer samples through the use of an endogenous retroviral long terminal repeat (LTR) promoter of the MSTD family. The LTR-IL-33 chimeric transcript encodes a novel shorter isoform of the protein, which is missing the initial N-terminus (including many conserved residues) of Native IL-33. In vitro studies showed that LTR-IL-33 expression is required for optimal CRC cell line growth as 3D colonospheres. Taken together, these data demonstrate the significance of TEs as regulators of aberrant gene expression in colon cancer.

Phosphorylation of histone H3 at serine 10 (H3S10ph) by Aurora kinases plays an important role in mitosis; however, H3S10ph also marks regulatory regions of inducible genes in interphase mammalian cells, implicating mitosis-independent functions. Using the fluorescent ubiquitin-mediated cell cycle indicator (FUCCI), we found that 30% of the genome in interphase mouse embryonic stem cells (ESCs) is marked with H3S10ph. H3S10ph broadly demarcates gene-rich regions in G1 and is positively correlated with domains of early DNA replication timing (RT) but negatively correlated with H3K9me2 and lamin-associated domains (LADs). Consistent with mitosis-independent kinase activity, this pattern was preserved in ESCs treated with Hesperadin, a potent inhibitor of Aurora B/C kinases. Disruption of H3S10ph by expression of nonphosphorylatable H3.3S10A results in ectopic spreading of H3K9me2 into adjacent euchromatic regions, mimicking the phenotype observed in Drosophila JIL-1 kinase mutants. Conversely, interphase H3S10ph domains expand in Ehmt1 (also known as Glp) null ESCs, revealing that H3S10ph deposition is restricted by H3K9me2. Strikingly, spreading of H3S10ph at RT transition regions (TTRs) is accompanied by aberrant transcription initiation of genes co-oriented with the replication fork in Ehmt1-/- and Ehmt2-/- ESCs, indicating that establishment of repressive chromatin on the leading strand following DNA synthesis may depend upon these lysine methyltransferases. H3S10ph is also anti-correlated with H3K9me2 in interphase murine embryonic fibroblasts (MEFs) and is restricted to intragenic regions of actively transcribing genes by EHMT2. Taken together, these observations reveal that H3S10ph may play a general role in restricting the spreading of repressive chromatin in interphase mammalian cells.

DNA methylation and histone H3 lysine 9 trimethylation (H3K9me3) play important roles in silencing of genes and retroelements. However, a comprehensive comparison of genes and repetitive elements repressed by these pathways has not been reported. Here we show that in mouse embryonic stem cells (mESCs), the genes upregulated after deletion of the H3K9 methyltransferase Setdb1 are distinct from those derepressed in mESC deficient in the DNA methyltransferases Dnmt1, Dnmt3a, and Dnmt3b, with the exception of a small number of primarily germline-specific genes. Numerous endogenous retroviruses (ERVs) lose H3K9me3 and are concomitantly derepressed exclusively in SETDB1 knockout mESCs. Strikingly, ~15% of upregulated genes are induced in association with derepression of promoter-proximal ERVs, half in the context of "chimeric" transcripts that initiate within these retroelements and splice to genic exons. Thus, SETDB1 plays a previously unappreciated yet critical role in inhibiting aberrant gene transcription by suppressing the expression of proximal ERVs.

Cell competition is emerging as a quality-control mechanism that eliminates unfit cells in a wide range of settings from development to the adult. However, the nature of the cells normally eliminated by cell competition and what triggers their elimination remains poorly understood. In mice, 35% of epiblast cells are eliminated before gastrulation. Here we show that cells with mitochondrial defects are eliminated by cell competition during early mouse development. Using single-cell transcriptional profiling of eliminated mouse epiblast cells, we identify hallmarks of cell competition and mitochondrial defects. We demonstrate that mitochondrial defects are common to a range of different loser cell types and that manipulating mitochondrial function triggers cell competition. Moreover, we show that in the mouse embryo, cell competition eliminates cells with sequence changes in mt-Rnr1 and mt-Rnr2, and that even non-pathological changes in mitochondrial DNA sequences can induce cell competition. Our results suggest that cell competition is a purifying selection that optimizes mitochondrial performance before gastrulation.

In mouse embryonic stem cells (mESCs), transcriptional silencing of numerous class I and II endogenous retroviruses (ERVs), including IAP, ETn and MMERVK10C, is dependent upon the H3K9 methyltransferase (KMTase) SETDB1/ESET and its binding partner KAP1/TRIM28. In contrast, the H3K9 KMTases G9a and GLP and HP1 proteins are dispensable for this process. Intriguingly, MERVL retroelements are actively transcribed exclusively in the two-cell (2C) embryo, but the molecular basis of silencing of these class III ERVs at later developmental stages has not been systematically addressed.

Genetic variation in the enzymes that catalyze posttranslational modification of proteins is a potentially important source of phenotypic variation during evolution. Ubiquitination is one such modification that affects turnover of virtually all of the proteins in the cell in addition to roles in signaling and epigenetic regulation. UBE2D3 is a promiscuous E2 enzyme, which acts as an ubiquitin donor for E3 ligases that catalyze ubiquitination of developmentally important proteins. We have used protein sequence comparison of UBE2D3 orthologs to identify a position in the C-terminal α-helical region of UBE2D3 that is occupied by a conserved serine in amniotes and by alanine in anamniote vertebrate and invertebrate lineages. Acquisition of the serine (S138) in the common ancestor to modern amniotes created a phosphorylation site for Aurora B. Phosphorylation of S138 disrupts the structure of UBE2D3 and reduces the level of the protein in mouse embryonic stem cells (ESCs). Substitution of S138 with the anamniote alanine (S138A) increases the level of UBE2D3 in ESCs as well as being a gain of function early embryonic lethal mutation in mice. When mutant S138A ESCs were differentiated into extraembryonic primitive endoderm, levels of the PDGFRα and FGFR1 receptor tyrosine kinases were reduced and primitive endoderm differentiation was compromised. Proximity ligation analysis showed increased interaction between UBE2D3 and the E3 ligase CBL and between CBL and the receptor tyrosine kinases. Our results identify a sequence change that altered the ubiquitination landscape at the base of the amniote lineage with potential effects on amniote biology and evolution.