Search

L1TD1 is a cytoplasmic RNA-binding protein specifically expressed in pluripotent stem cells and, unlike its mouse ortholog, is essential for the maintenance of stemness in human cells. Although L1TD1 is the only known protein-coding gene domesticated from a LINE-1 (L1) retroelement, the functional legacy of its ancestral protein, ORF1p of L1, and how it is manifested in L1TD1 are still unknown. Here, we determined RNAs associated with L1TD1 and found that, like ORF1p, L1TD1 binds L1 RNAs and localizes to high-density ribonucleoprotein (RNP) condensates. Unexpectedly, L1TD1 enhanced the translation of a subset of mRNAs enriched in the condensates. L1TD1 depletion promoted the formation of stress granules in embryonic stem cells. In HeLa cells, ectopically expressed L1TD1 facilitated the dissolution of stress granules and granules formed by pathological mutations of TDP-43 and FUS. The glutamate-rich domain and the ORF1-homology domain of L1TD1 facilitated dispersal of the RNPs and induced autophagy, respectively. These results provide insights into how L1TD1 regulates gene expression in pluripotent stem cells. We propose that the ability of L1TD1 to dissolve stress granules may provide novel opportunities for treatment of neurodegenerative diseases caused by disturbed stress granule dynamics., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Bromodomain-containing protein 9 (BRD9) is a specific subunit of the non-canonical SWI/SNF (ncBAF) chromatin-remodeling complex, whose function in human embryonic stem cells (hESCs) remains unclear. Here, we demonstrate that impaired BRD9 function reduces the self-renewal capacity of hESCs and alters their differentiation potential. Specifically, BRD9 depletion inhibits meso-endoderm differentiation while promoting neural ectoderm differentiation. Notably, supplementation of NODAL, TGF-β, Activin A or WNT3A rescues the differentiation defects caused by BRD9 loss. Mechanistically, BRD9 forms a complex with BRD4, SMAD2/3, β-CATENIN and P300, which regulates the expression of pluripotency genes and the activity of TGF-β/Nodal/Activin and Wnt signaling pathways. This is achieved by regulating the deposition of H3K27ac on associated genes, thus maintaining and directing hESC differentiation. BRD9-mediated regulation of the TGF-β/Activin/Nodal pathway is also demonstrated in the development of pancreatic and breast cancer cells. In summary, our study highlights the crucial role of BRD9 in the regulation of hESC self-renewal and differentiation, as well as its participation in the progression of pancreatic and breast cancers., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Dot1l is a histone methyltransferase without a SET domain and is responsible for H3K79 methylation, which marks active transcription. In contradiction, Dot1l also participates in silencing gene expression. The target regions and mechanism of Dot1l in repressing transcription remain enigmatic. Here, we show that Dot1l represses endogenous retroviruses in embryonic stem cells (ESCs). Specifically, the absence of Dot1l led to the activation of MERVL, which is a marker of 2-cell-like cells. In addition, Dot1l deletion activated the 2-cell-like state and predisposed ESCs to differentiate into trophectoderm lineage. Transcriptome analysis revealed activation of 2-cell genes and meiotic genes by Dot1l deletion. Mechanistically, Dot1l interacted with and co-localized with Npm1 on MERVL, and depletion of Npm1 similarly augmented MERVL expression. The catalytic activity and AT-hook domain of Dot1l are important to suppress MERVL. Notably, Dot1l-Npm1 restricts MERVL by regulating protein level and deposition of histone H1. Furthermore, Dot1l is critical for Npm1 to efficiently interact with histone H1 and inhibit ubiquitination of H1 whereas Npm1 is essential for Dot1l to interact with MERVL. Altogether, we discover that Dot1l represses MERVL through chaperoning H1 by collaborating with Npm1. Importantly, our findings shed light on the non-canonical transcriptional repressive role of Dot1l in ESCs., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

An engineered SOX17 variant with point mutations within its DNA binding domain termed SOX17FNV is a more potent pluripotency inducer than SOX2, yet the underlying mechanism remains unclear. Although wild-type SOX17 was incapable of inducing pluripotency, SOX17FNV outperformed SOX2 in mouse and human pluripotency reprogramming. In embryonic stem cells, SOX17FNV could replace SOX2 to maintain pluripotency despite considerable sequence differences and upregulated genes expressed in cleavage-stage embryos. Mechanistically, SOX17FNV co-bound OCT4 more cooperatively than SOX2 in the context of the canonical SoxOct DNA element. SOX2, SOX17, and SOX17FNV were all able to bind nucleosome core particles in vitro, which is a prerequisite for pioneer transcription factors. Experiments using purified proteins and in cellular contexts showed that SOX17 variants phase-separated more efficiently than SOX2, suggesting an enhanced ability to self-organise. Systematic deletion analyses showed that the N-terminus of SOX17FNV was dispensable for its reprogramming activity. However, the C-terminus encodes essential domains indicating multivalent interactions that drive transactivation and reprogramming. We defined a minimal SOX17FNV (miniSOX) that can support reprogramming with high activity, reducing the payload of reprogramming cassettes. This study uncovers the mechanisms behind SOX17FNV-induced pluripotency and establishes engineered SOX factors as powerful cell engineering tools., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The establishment of cellular identity is driven by transcriptional and epigenetic regulators of the chromatin proteome - the chromatome. Comprehensive analyses of the chromatome composition and dynamics can therefore greatly improve our understanding of gene regulatory mechanisms. Here, we developed an accurate mass spectrometry (MS)-based proteomic method called Chromatin Aggregation Capture (ChAC) followed by Data-Independent Acquisition (DIA) and analyzed chromatome reorganizations during major phases of pluripotency. This enabled us to generate a comprehensive atlas of proteomes, chromatomes, and chromatin affinities for the ground, formative and primed pluripotency states, and to pinpoint the specific binding and rearrangement of regulatory components. These comprehensive datasets combined with extensive analyses identified phase-specific factors like QSER1 and JADE1/2/3 and provide a detailed foundation for an in-depth understanding of mechanisms that govern the phased progression of pluripotency. The technical advances reported here can be readily applied to other models in development and disease., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The histone H3 variant, H3.3, is localized at specific regions in the genome, especially promoters and active enhancers, and has been shown to play important roles in development. A lysine to methionine substitution in position 27 (H3.3K27M) is a main cause of Diffuse Intrinsic Pontine Glioma (specifically Diffuse Midline Glioma, K27M-mutant), a lethal type of pediatric cancer. H3.3K27M has a dominant-negative effect by inhibiting the Polycomb Repressor Complex 2 (PRC2) activity. Here, we studied the immediate, genome-wide, consequences of the H3.3K27M mutation independent of PRC2 activity. We developed Doxycycline (Dox)-inducible mouse embryonic stem cells (ESCs) carrying a single extra copy of WT-H3.3, H3.3K27M and H3.3K27L, all fused to HA. We performed RNA-Seq and ChIP-Seq at different times following Dox induction in undifferentiated and differentiated ESCs. We find increased binding of H3.3 around transcription start sites in cells expressing both H3.3K27M and H3.3K27L compared with WT, but not in cells treated with PRC2 inhibitors. Differentiated cells carrying either H3.3K27M or H3.3K27L retain expression of ESC-active genes, in expense of expression of genes related to neuronal differentiation. Taken together, our data suggest that a modifiable H3.3K27 is required for proper histone incorporation and cellular maturation, independent of PRC2 activity., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

RNA helicases are involved in multiple steps of RNA metabolism to direct their roles in gene expression, yet their functions in pluripotency control remain largely unexplored. Starting from an RNA interference (RNAi) screen of RNA helicases, we identified that eIF4A3, a DEAD-box (Ddx) helicase component of the exon junction complex (EJC), is essential for the maintenance of embryonic stem cells (ESCs). Mechanistically, we show that eIF4A3 post-transcriptionally controls the pluripotency-related cell cycle regulators and that its depletion causes the loss of pluripotency via cell cycle dysregulation. Specifically, eIF4A3 is required for the efficient nuclear export of Ccnb1 mRNA, which encodes Cyclin B1, a key component of the pluripotency-promoting pathway during the cell cycle progression of ESCs. Our results reveal a previously unappreciated role for eIF4A3 and its associated EJC in maintaining stem cell pluripotency through post-transcriptional control of the cell cycle., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Neural differentiation of embryonic stem cells (ESCs) requires precisely orchestrated gene regulation, a process governed in part by changes in 3D chromatin structure. How these changes regulate gene expression in this context remains unclear. In this study, we observed enrichment of the transcription factor KLF4 at some poised or closed enhancers at TSS-linked regions of genes associated with neural differentiation. Combination analysis of ChIP, HiChIP and RNA-seq data indicated that KLF4 loss in ESCs induced changes in 3D chromatin structure, including increased chromatin interaction loops between neural differentiation-associated genes and active enhancers, leading to upregulated expression of neural differentiation-associated genes and therefore early neural differentiation. This study suggests KLF4 inhibits early neural differentiation by regulation of 3D chromatin structure, which is a new mechanism of early neural differentiation., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Although extended pluripotent stem cells (EPSCs) have the potential to form both embryonic and extraembryonic lineages, how their transcriptional regulatory mechanism differs from that of embryonic stem cells (ESCs) remains unclear. Here, we discovered that YY1 binds to specific open chromatin regions in EPSCs. Yy1 depletion in EPSCs leads to a gene expression pattern more similar to that of ESCs than control EPSCs. Moreover, Yy1 depletion triggers a series of epigenetic crosstalk activities, including changes in DNA methylation, histone modifications and high-order chromatin structures. Yy1 depletion in EPSCs disrupts the enhancer-promoter (EP) interactions of EPSC-specific genes, including Dnmt3l. Yy1 loss results in DNA hypomethylation and dramatically reduces the enrichment of H3K4me3 and H3K27ac on the promoters of EPSC-specific genes by upregulating the expression of Kdm5c and Hdac6 through facilitating the formation of CCCTC-binding factor (CTCF)-mediated EP interactions surrounding their loci. Furthermore, single-cell RNA sequencing (scRNA-seq) experiments revealed that YY1 is required for the derivation of extraembryonic endoderm (XEN)-like cells from EPSCs in vitro. Together, this study reveals that YY1 functions as a key regulator of multidimensional epigenetic crosstalk associated with extended pluripotency., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Mammalian early epiblasts at different phases are characterized by naïve, formative, and primed pluripotency states, involving extensive transcriptome changes. Here, we report that deadenylase Cnot8 of Ccr4-Not complex plays essential roles during the transition from naïve to formative state. Knock out (KO) Cnot8 resulted in early embryonic lethality in mice, but Cnot8 KO embryonic stem cells (ESCs) could be established. Compared with the cells differentiated from normal ESCs, Cnot8 KO cells highly expressed a great many genes during their differentiation into the formative state, including several hundred naïve-like genes enriched in lipid metabolic process and gene expression regulation that may form the naïve regulation networks. Knockdown expression of the selected genes of naïve regulation networks partially rescued the differentiation defects of Cnot8 KO ESCs. Cnot8 depletion led to the deadenylation defects of its targets, increasing their poly(A) tail lengths and half-life, eventually elevating their expression levels. We further found that Cnot8 was involved in the clearance of targets through its deadenylase activity and the binding of Ccr4-Not complex, as well as the interacting with Tob1 and Pabpc1. Our results suggest that Cnot8 eliminates naïve regulation networks through mRNA clearance, and is essential for naïve-to-formative pluripotency transition., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Tet enzymes (Tet1/2/3) oxidize 5-methylcytosine to promote DNA demethylation and partner with chromatin modifiers to regulate gene expression. Tet1 is highly expressed in embryonic stem cells (ESCs), but its enzymatic and non-enzymatic roles in gene regulation are not dissected. We have generated Tet1 catalytically inactive (Tet1m/m) and knockout (Tet1-/-) ESCs and mice to study these functions. Loss of Tet1, but not loss of its catalytic activity, caused aberrant upregulation of bivalent (H3K4me3+; H3K27me3+) developmental genes, leading to defects in differentiation. Wild-type and catalytic-mutant Tet1 occupied similar genomic loci which overlapped with H3K27 tri-methyltransferase PRC2 and the deacetylase complex Sin3a at promoters of bivalent genes and with the helicase Chd4 at active genes. Loss of Tet1, but not loss of its catalytic activity, impaired enrichment of PRC2 and Sin3a at bivalent promoters leading to reduced H3K27 trimethylation and deacetylation, respectively, in absence of any changes in DNA methylation. Tet1-/-, but not Tet1m/m, embryos expressed higher levels of Gata6 and were developmentally delayed. Thus, the critical functions of Tet1 in ESCs and early development are mediated through its non-catalytic roles in regulating H3K27 modifications to silence developmental genes, and are more important than its catalytic functions in DNA demethylation., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

CMTR1 (cap methyltransferase 1) catalyses methylation of the first transcribed nucleotide of RNAPII transcripts (N1 2'-O-Me), creating part of the mammalian RNA cap structure. In addition to marking RNA as self, N1 2'-O-Me has ill-defined roles in RNA expression and translation. Here, we investigated the gene specificity of CMTR1 and its impact on RNA expression in embryonic stem cells. Using chromatin immunoprecipitation, CMTR1 was found to bind to transcription start sites (TSS) correlating with RNAPII levels, predominantly binding at histone genes and ribosomal protein (RP) genes. Repression of CMTR1 expression resulted in repression of RNAPII binding at the TSS and repression of RNA expression, particularly of histone and RP genes. In correlation with regulation of histones and RP genes, CMTR1 repression resulted in repression of translation and induction of DNA replication stress and damage. Indicating a direct role for CMTR1 in transcription, addition of recombinant CMTR1 to purified nuclei increased transcription of the histone and RP genes. CMTR1 was found to be upregulated during neural differentiation and there was an enhanced requirement for CMTR1 for gene expression and proliferation during this process. We highlight the distinct roles of the cap methyltransferases RNMT and CMTR1 in target gene expression and differentiation., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Transposable elements are an abundant source of transcription factor binding sites, and favorable genomic integration may lead to their recruitment by the host genome for gene regulatory functions. However, it is unclear how frequent co-option of transposable elements as regulatory elements is, to which regulatory programs they contribute and how they compare to regulatory elements devoid of transposable elements. Here, we report a transcription initiation-centric, in-depth characterization of the transposon-derived regulatory landscape of mouse embryonic stem cells. We demonstrate that a substantial number of transposable element insertions, in particular endogenous retroviral elements, are associated with open chromatin regions that are divergently transcribed into unstable RNAs in a cell-type specific manner, and that these elements contribute to a sizable proportion of active enhancers and gene promoters. We further show that transposon subfamilies contribute differently and distinctly to the pluripotency regulatory program through their repertoires of transcription factor binding site sequences, shedding light on the formation of regulatory programs and the origins of regulatory elements., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Precise genomic modification using prime editing (PE) holds enormous potential for research and clinical applications. In this study, we generated all-in-one prime editing (PEA1) constructs that carry all the components required for PE, along with a selection marker. We tested these constructs (with selection) in HEK293T, K562, HeLa and mouse embryonic stem (ES) cells. We discovered that PE efficiency in HEK293T cells was much higher than previously observed, reaching up to 95% (mean 67%). The efficiency in K562 and HeLa cells, however, remained low. To improve PE efficiency in K562 and HeLa, we generated a nuclease prime editor and tested this system in these cell lines as well as mouse ES cells. PE-nuclease greatly increased prime editing initiation, however, installation of the intended edits was often accompanied by extra insertions derived from the repair template. Finally, we show that zygotic injection of the nuclease prime editor can generate correct modifications in mouse fetuses with up to 100% efficiency., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Whole genome base-resolution methylome sequencing allows for the most comprehensive analysis of DNA methylation, however, the considerable sequencing cost often limits its applications. While reduced representation sequencing can be an affordable alternative, over 80% of CpGs in the genome are not covered. Building on our recently developed TET-assisted pyridine borane sequencing (TAPS) method, we here described endonuclease enrichment TAPS (eeTAPS), which utilizes dihydrouracil (DHU)-cleaving endonuclease digestion of TAPS-converted DNA to enrich methylated CpG sites (mCpGs). eeTAPS can accurately detect 87% of mCpGs in the mouse genome with a sequencing depth equivalent to 4× whole genome sequencing. In comparison, reduced representation TAPS (rrTAPS) detected less than 4% of mCpGs with 2.5× sequencing depth. Our results demonstrate eeTAPS to be a new strategy for cost-effective genome-wide methylation analysis at single-CpG resolution that can fill the gap between whole-genome and reduced representation sequencing., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Heterochromatin binding protein HP1β plays an important role in chromatin organization and cell differentiation, however the underlying mechanisms remain unclear. Here, we generated HP1β-/- embryonic stem cells and observed reduced heterochromatin clustering and impaired differentiation. We found that during stem cell differentiation, HP1β is phosphorylated at serine 89 by CK2, which creates a binding site for the pluripotency regulator KAP1. This phosphorylation dependent sequestration of KAP1 in heterochromatin compartments causes a downregulation of pluripotency factors and triggers pluripotency exit. Accordingly, HP1β-/- and phospho-mutant cells exhibited impaired differentiation, while ubiquitination-deficient KAP1-/- cells had the opposite phenotype with enhanced differentiation. These results suggest that KAP1 regulates pluripotency via its ubiquitination activity. We propose that the formation of subnuclear membraneless heterochromatin compartments may serve as a dynamic reservoir to trap or release cellular factors. The sequestration of essential regulators defines a novel and active role of heterochromatin in gene regulation and represents a dynamic mode of remote control to regulate cellular processes like cell fate decisions., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Using human embryonic, adult and cancer stem cells/stem cell-like cells (SCs), we demonstrate that DNA replication speed differs in SCs and their differentiated counterparts. While SCs decelerate DNA replication, differentiated cells synthesize DNA faster and accumulate DNA damage. Notably, both replication phenotypes depend on p53 and polymerase iota (POLι). By exploring protein interactions and newly synthesized DNA, we show that SCs promote complex formation of p53 and POLι at replication sites. Intriguingly, in SCs the translocase ZRANB3 is recruited to POLι and required for slow-down of DNA replication. The known role of ZRANB3 in fork reversal suggests that the p53-POLι complex mediates slow but safe bypass of replication barriers in SCs. In differentiated cells, POLι localizes more transiently to sites of DNA synthesis and no longer interacts with p53 facilitating fast POLι-dependent DNA replication. In this alternative scenario, POLι associates with the p53 target p21, which antagonizes PCNA poly-ubiquitination and, thereby potentially disfavors the recruitment of translocases. Altogether, we provide evidence for diametrically opposed DNA replication phenotypes in SCs and their differentiated counterparts putting DNA replication-based strategies in the spotlight for the creation of therapeutic opportunities targeting SCs., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Nuclear architecture influences gene regulation and cell identity by controlling the three-dimensional organization of genes and their distal regulatory sequences, which may be far apart in linear space. The genome is functionally and spatially segregated in the eukaryotic nucleus with transcriptionally active regions in the nuclear interior separated from repressive regions, including those at the nuclear periphery. Here, we describe the identification of a novel type of nuclear peripheral chromatin domain that is enriched for tissue-specific transcriptional enhancers. Like other chromatin at the nuclear periphery, these regions are marked by H3K9me2. But unlike the nuclear peripheral Lamina-Associated Domains (LADs), these novel, enhancer-rich domains have limited Lamin B interaction. We therefore refer to them as H3K9me2-Only Domains (KODs). In mouse embryonic stem cells, KODs are found in Hi-C-defined A compartments and feature relatively accessible chromatin. KODs are characterized by low gene expression and enhancers located in these domains bear the histone marks of an inactive or poised state. These results indicate that KODs organize a subset of inactive, tissue-specific enhancers at the nuclear periphery. We hypothesize that KODs may play a role in facilitating and perhaps constraining the enhancer-promoter interactions underlying spatiotemporal regulation of gene expression programs in differentiation and development., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Epigenetics, especially histone marks, functions beyond the DNA sequences to regulate gene expression. Depletion of NSD1, which catalyzes H3K36me2, leads to both up- and down-regulation of gene expression, indicating NSD1 is associated with both active and repressed gene expression. It's known that NSD1 regulates the deposition and expansion of H3K27me3, a repressive mark for gene expression, to keep active gene transcription. However, how NSD1 functions to repress gene expression is largely unknown. Here, we find that, when NSD1 is knocked out in mouse embryonic stem cells (mESCs), H3K27ac increases correlatively with the decrease of H3K36me2 at active enhancers, which is associated with mesoderm differentiation genes, leading to elevated gene expression. Mechanistically, NSD1 recruits HDAC1, the deacetylase of H3K27ac, to chromatin. Moreover, HDAC1 knockout (KO) recapitulates the increase of H3K27ac at active enhancers as the NSD1 depletion. Together, we propose that NSD1 deposits H3K36me2 and recruits HDAC1 at active enhancers to serve as a 'safeguard', preventing further activation of active enhancer-associated genes., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The genetic code of mammalian cells can be expanded to allow the incorporation of non-canonical amino acids (ncAAs) by suppressing in-frame amber stop codons (UAG) with an orthogonal pyrrolysyl-tRNA synthetase (PylRS)/tRNAPylCUA (PylT) pair. However, the feasibility of this approach is substantially hampered by unpredictable variations in incorporation efficiencies at different stop codon positions within target proteins. Here, we apply a proteomics-based approach to quantify ncAA incorporation rates at hundreds of endogenous amber stop codons in mammalian cells. With these data, we compute iPASS (Identification of Permissive Amber Sites for Suppression; available at www.bultmannlab.eu/tools/iPASS), a linear regression model to predict relative ncAA incorporation efficiencies depending on the surrounding sequence context. To verify iPASS, we develop a dual-fluorescence reporter for high-throughput flow-cytometry analysis that reproducibly yields context-specific ncAA incorporation efficiencies. We show that nucleotides up- and downstream of UAG synergistically influence ncAA incorporation efficiency independent of cell line and ncAA identity. Additionally, we demonstrate iPASS-guided optimization of ncAA incorporation rates by synonymous exchange of codons flanking the amber stop codon. This combination of in silico analysis followed by validation in living mammalian cells substantially simplifies identification as well as adaptation of sites within a target protein to confer high ncAA incorporation rates., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Maintenance of stem-cell identity requires proper regulation of enhancer activity. Both transcription factors OCT4/SOX2/NANOG and histone methyltransferase complexes MLL/SET1 were shown to regulate enhancer activity, but how they are regulated in embryonic stem cells (ESCs) remains further studies. Here, we report a transcription factor BACH1, which directly interacts with OCT4/SOX2/NANOG (OSN) and MLL/SET1 methyltransferase complexes and maintains pluripotency in mouse ESCs (mESCs). BTB domain and bZIP domain of BACH1 are required for these interactions and pluripotency maintenance. Loss of BACH1 reduced the interaction between NANOG and MLL1/SET1 complexes, and decreased their occupancy on chromatin, and further decreased H3 lysine 4 trimethylation (H3K4me3) level on gene promoters and (super-) enhancers, leading to decreased enhancer activity and transcription activity, especially on stemness-related genes. Moreover, BACH1 recruited NANOG through chromatin looping and regulated remote NANOG binding, fine-tuning enhancer-promoter activity and gene expression. Collectively, these observations suggest that BACH1 maintains pluripotency in ESCs by recruiting NANOG and MLL/SET1 complexes to chromatin and maintaining the trimethylated state of H3K4 and enhancer-promoter activity, especially on stemness-related genes., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Long noncoding RNAs (lncRNAs) play a wide range of roles in the epigenetic regulation of crucial biological processes, but the functions of lncRNAs in cortical development are poorly understood. Using human embryonic stem cell (hESC)-based 2D neural differentiation approach and 3D cerebral organoid system, we identified that the lncRNA PAUPAR, which is adjacent to PAX6, plays essential roles in cortical differentiation by interacting with PAX6 to regulate the expression of a large number of neural genes. Mechanistic studies showed that PAUPAR confers PAX6 proper binding sites on the target neural genes by directly binding the genomic regions of these genes. Moreover, PAX6 recruits the histone methyltransferase NSD1 through its C-terminal PST enrichment domain, then regulate H3K36 methylation and the expression of target genes. Collectively, our data reveal that the PAUPAR/PAX6/NSD1 complex plays a critical role in the epigenetic regulation of hESC cortical differentiation and highlight the importance of PAUPAR as an intrinsic regulator of cortical differentiation., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Enhancers play important roles in controlling gene expression in a choreographed spatial and temporal manner during development. However, it is unclear how these regulatory regions are established during differentiation. Here we investigated the genome-wide binding profile of the forkhead transcription factor FOXK2 in human embryonic stem cells (ESCs) and downstream cell types. This transcription factor is bound to thousands of regulatory regions in human ESCs, and binding at many sites is maintained as cells differentiate to mesendodermal and neural precursor cell (NPC) types, alongside the emergence of new binding regions. FOXK2 binding is generally associated with active histone marks in any given cell type. Furthermore newly acquired, or retained FOXK2 binding regions show elevated levels of activating histone marks following differentiation to NPCs. In keeping with this association with activating marks, we demonstrate a role for FOXK transcription factors in gene activation during NPC differentiation. FOXK2 occupancy in ESCs is therefore an early mark for delineating the regulatory regions, which become activated in later lineages., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Gene-editing experiments commonly elicit the error-prone non-homologous end joining for DNA double-strand break (DSB) repair. Microhomology-mediated end joining (MMEJ) can generate more predictable outcomes for functional genomic and somatic therapeutic applications. We compared three DSB repair prediction algorithms - MENTHU, inDelphi, and Lindel - in identifying MMEJ-repaired, homogeneous genotypes (PreMAs) in an independent dataset of 5,885 distinct Cas9-mediated mouse embryonic stem cell DSB repair events. MENTHU correctly identified 46% of all PreMAs available, a ∼2- and ∼60-fold sensitivity increase compared to inDelphi and Lindel, respectively. In contrast, only Lindel correctly predicted predominant single-base insertions. We report the new algorithm MENdel, a combination of MENTHU and Lindel, that achieves the most predictive coverage of homogeneous out-of-frame mutations in this large dataset. We then estimated the frequency of Cas9-targetable homogeneous frameshift-inducing DSBs in vertebrate coding regions for gene discovery using MENdel. 47 out of 54 genes (87%) contained at least one early frameshift-inducing DSB and 49 out of 54 (91%) did so when also considering Cas12a-mediated deletions. We suggest that the use of MENdel helps researchers use MMEJ at scale for reverse genetics screenings and with sufficient intra-gene density rates to be viable for nearly all loss-of-function based gene editing therapeutic applications., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Mammalian de novo DNA methyltransferases (DNMT) are responsible for the establishment of cell-type-specific DNA methylation in healthy and diseased tissues. Through genome-wide analysis of de novo methylation activity in murine stem cells we uncover that DNMT3A prefers to methylate CpGs followed by cytosines or thymines, while DNMT3B predominantly methylates CpGs followed by guanines or adenines. These signatures are further observed at non-CpG sites, resembling methylation context observed in specialised cell types, including neurons and oocytes. We further show that these preferences result from structural differences in the catalytic domains of the two de novo DNMTs and are not a consequence of differential recruitment to the genome. Molecular dynamics simulations suggest that, in case of human DNMT3A, the preference is due to favourable polar interactions between the flexible Arg836 side chain and the guanine that base-pairs with the cytosine following the CpG. By exchanging arginine to a lysine, the corresponding side chain in DNMT3B, the sequence preference is reversed, confirming the requirement for arginine at this position. This context-dependent enzymatic activity provides additional insights into the complex regulation of DNA methylation patterns., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The heterochromatin protein 1 (HP1) family members are canonical effectors and propagators of gene repression mediated by histone H3 lysine 9 (H3K9) methylation. HP1γ exhibits an increased interaction with active transcription elongation-associated factors in embryonic stem cells (ESCs) compared to somatic cells. However, whether this association has a functional consequence remains elusive. Here we find that genic HP1γ colocalizes and enhances enrichment of transcription elongation-associated H3K36me3 rather than H3K9me3. Unexpectedly, sustained H3K36me3 deposition is dependent on HP1γ. HP1γ-deleted ESCs display reduced H3K36me3 enrichment, concomitant with decreased expression at shared genes which function to maintain cellular homeostasis. Both the H3K9me3-binding chromodomain and histone binding ability of HP1γ are dispensable for maintaining H3K36me3 levels. Instead, the chromoshadow together with the hinge domain of HP1γ that confer protein and nucleic acid-binding ability are sufficient because they retain the ability to interact with NSD1, an H3K36 methyltransferase. HP1γ-deleted ESCs have a slower self-renewal rate and an impaired ability to differentiate towards cardiac mesoderm. Our findings reveal a requirement for HP1γ in faithful establishment of transcription elongation in ESCs, which regulates pluripotency., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

In vertebrates, DNA methylation predominantly occurs at CG dinucleotides however, widespread non-CG methylation (mCH) has been reported in mammalian embryonic stem cells and in the brain. In mammals, mCH is found at CAC trinucleotides in the nervous system, where it is associated with transcriptional repression, and at CAG trinucleotides in embryonic stem cells, where it positively correlates with transcription. Moreover, CAC methylation appears to be a conserved feature of adult vertebrate brains. Unlike any of those methylation signatures, here we describe a novel form of mCH that occurs in the TGCT context within zebrafish mosaic satellite repeats. TGCT methylation is inherited from both male and female gametes, remodelled during mid-blastula transition, and re-established during gastrulation in all embryonic layers. Moreover, we identify DNA methyltransferase 3ba (Dnmt3ba) as the primary enzyme responsible for the deposition of this mCH mark. Finally, we observe that TGCT-methylated repeats are specifically associated with H3K9me3-marked heterochromatin suggestive of a functional interplay between these two gene-regulatory marks. Altogether, this work provides insight into a novel form of vertebrate mCH and highlights the substrate diversity of vertebrate DNA methyltransferases., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is used to identify genome-wide DNA regions bound by proteins. Given one ChIP-seq experiment with replicates, binding sites not observed in all the replicates will usually be interpreted as noise and discarded. However, the recent discovery of high-occupancy target (HOT) regions suggests that there are regions where binding of multiple transcription factors can be identified. To investigate ChIP-seq variability, we developed a reproducibility score and a method that identifies cell-specific variable regions in ChIP-seq data by integrating replicated ChIP-seq experiments for multiple protein targets on a particular cell type. Using our method, we found variable regions in human cell lines K562, GM12878, HepG2, MCF-7 and in mouse embryonic stem cells (mESCs). These variable-occupancy target regions (VOTs) are CG dinucleotide rich, and show enrichment at promoters and R-loops. They overlap significantly with HOT regions, but are not blacklisted regions producing non-specific binding ChIP-seq peaks. Furthermore, in mESCs, VOTs are conserved among placental species suggesting that they could have a function important for this taxon. Our method can be useful to point to such regions along the genome in a given cell type of interest, to improve the downstream interpretative analysis before follow-up experiments., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The miRNA pathway has three segments-biogenesis, targeting and downstream regulatory effectors. We aimed to better understand their cellular control by exploring the miRNA-mRNA-targeting relationships. We first used human evolutionarily conserved sites. Strikingly, AGOs 1-3 are all among the top 14 mRNAs with the highest miRNA site counts, along with ANKRD52, the phosphatase regulatory subunit of the recently identified AGO phosphorylation cycle; and the AGO phosphorylation cycle mRNAs share much more than expected miRNA sites. The mRNAs for TNRC6, which acts with AGOs to channel miRNA-mediated regulatory actions onto specific mRNAs, are also heavily miRNA-targeted. In contrast, upstream miRNA biogenesis mRNAs are not, and neither are downstream regulatory effectors. In short, binding site enrichment in miRNA targeting machinery mRNAs, but neither upstream biogenesis nor downstream effector mRNAs, was observed, endowing a cellular capacity for intensive and specific feedback control of the targeting activity. The pattern was confirmed with experimentally determined miRNA-mRNA target relationships. Moreover, genetic experiments demonstrated cellular utilization of this capacity. Thus, we uncovered a capacity for intensive, and specific, feedback-regulation of miRNA targeting activity directly by miRNAs themselves, i.e. segment-specific feedback auto-regulation of miRNA pathway, complementing miRNAs pairing with transcription factors to form hybrid feedback-loop., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

We have previously shown that the highly prevalent acute myeloid leukemia (AML) mutation, Arg882His, in DNMT3A disrupts its cooperative mechanism and leads to reduced enzymatic activity, thus explaining the genomic hypomethylation in AML cells. However, the underlying cause of the oncogenic effect of Arg882His in DNMT3A is not fully understood. Here, we discovered that DNMT3A WT enzyme under conditions that favor non-cooperative kinetic mechanism as well as DNMT3A Arg882His variant acquire CpG flanking sequence preference akin to that of DNMT3B, which is non-cooperative. We tested if DNMT3A Arg882His could preferably methylate DNMT3B-specific target sites in vivo. Rescue experiments in Dnmt3a/3b double knockout mouse embryonic stem cells show that the corresponding Arg878His mutation in mouse DNMT3A severely impairs its ability to methylate major satellite DNA, a DNMT3A-preferred target, but has no overt effect on the ability to methylate minor satellite DNA, a DNMT3B-preferred target. We also observed a previously unappreciated CpG flanking sequence bias in major and minor satellite repeats that is consistent with DNMT3A and DNMT3B specificity suggesting that DNA methylation patterns are guided by the sequence preference of these enzymes. We speculate that aberrant methylation of DNMT3B target sites could contribute to the oncogenic potential of DNMT3A AML variant., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Retinoic acid (RA) induces rapid differentiation of embryonic stem cells (ESCs), partly by activating expression of the transcription factor Hoxa1, which regulates downstream target genes that promote ESCs differentiation. However, mechanisms of RA-induced Hoxa1 expression and ESCs early differentiation remain largely unknown. Here, we identify a distal enhancer interacting with the Hoxa1 locus through a long-range chromatin loop. Enhancer deletion significantly inhibited expression of RA-induced Hoxa1 and endoderm master control genes such as Gata4 and Gata6. Transcriptome analysis revealed that RA-induced early ESCs differentiation was blocked in Hoxa1 enhancer knockout cells, suggesting a requirement for the enhancer. Restoration of Hoxa1 expression partly rescued expression levels of ∼40% of genes whose expression changed following enhancer deletion, and ∼18% of promoters of those rescued genes were directly bound by Hoxa1. Our data show that a distal enhancer maintains Hoxa1 expression through long-range chromatin loop and that Hoxa1 directly regulates downstream target genes expression and then orchestrates RA-induced early differentiation of ESCs. This discovery reveals mechanisms of a novel enhancer regulating RA-induced Hoxa genes expression and early ESCs differentiation., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Chromatin immunoprecipitation (ChIP) is the most widely used approach for identification of genome-associated proteins and their modifications. We have previously introduced a microplate-based ChIP platform, Matrix ChIP, where the entire ChIP procedure is done on the same plate without sample transfers. Compared to conventional ChIP protocols, the Matrix ChIP assay is faster and has increased throughput. However, even with microplate ChIP assays, sample preparation and chromatin fragmentation (which is required to map genomic locations) remains a major bottleneck. We have developed a novel technology (termed 'PIXUL') utilizing an array of ultrasound transducers for simultaneous shearing of samples in standard 96-well microplates. We integrated PIXUL with Matrix ChIP ('PIXUL-ChIP'), that allows for fast, reproducible, low-cost and high-throughput sample preparation and ChIP analysis of 96 samples (cell culture or tissues) in one day. Further, we demonstrated that chromatin prepared using PIXUL can be used in an existing ChIP-seq workflow. Thus, the high-throughput capacity of PIXUL-ChIP provides the means to carry out ChIP-qPCR or ChIP-seq experiments involving dozens of samples. Given the complexity of epigenetic processes, the use of PIXUL-ChIP will advance our understanding of these processes in health and disease, as well as facilitate screening of epigenetic drugs., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Cancer-associated mutations of the core splicing factor 3 B1 (SF3B1) result in selection of novel 3' splice sites (3'SS), but precise molecular mechanisms of oncogenesis remain unclear. SF3B1 stabilizes the interaction between U2 snRNP and branch point (BP) on the pre-mRNA. It has hence been speculated that a change in BP selection is the basis for novel 3'SS selection. Direct quantitative determination of BP utilization is however technically challenging. To define BP utilization by SF3B1-mutant spliceosomes, we used an overexpression approach in human cells as well as a complementary strategy using isogenic murine embryonic stem cells with monoallelic K700E mutations constructed via CRISPR/Cas9-based genome editing and a dual vector homology-directed repair methodology. A synthetic minigene library with degenerate regions in 3' intronic regions (3.4 million individual minigenes) was used to compare BP usage of SF3B1K700E and SF3B1WT. Using this model, we show that SF3B1K700E spliceosomes utilize non-canonical sequence variants (at position -1 relative to BP adenosine) more frequently than wild-type spliceosomes. These predictions were confirmed using minigene splicing assays. Our results suggest a model of BP utilization by mutant SF3B1 wherein it is able to utilize non-consensus alternative BP sequences by stabilizing weaker U2-BP interactions.

The controlled and stepwise oxidation of 5mC to 5hmC, 5fC and 5caC by Tet enzymes is influencing the chemical and biological properties of cytosine. Besides direct effects on gene regulation, oxidised forms influence the dynamics of demethylation and re-methylation processes. So far, no combined methods exist which allow to precisely determine the strand specific localisation of cytosine modifications along with their CpG symmetric distribution. Here we describe a comprehensive protocol combining conventional hairpin bisulfite with oxidative bisulfite sequencing (HPoxBS) to determine the strand specific distribution of 5mC and 5hmC at base resolution. We apply this method to analyse the contribution of local oxidative effects on DNA demethylation in mouse ES cells. Our method includes the HPoxBS workflow and subsequent data analysis using our developed software tools. Besides a precise estimation and display of strand specific 5mC and 5hmC levels at base resolution we apply the data to predict region specific activities of Dnmt and Tet enzymes. Our experimental and computational workflow provides a precise double strand display of 5mC and 5hmC modifications at single base resolution. Based on our data we predict region specific Tet and Dnmt enzyme efficiencies shaping the distinct locus levels and patterns of 5hmC and 5mC.

It is becoming increasingly clear that chromosome organization plays an important role in gene regulation. High-resolution methods such as 4C, Capture-C and promoter capture Hi-C (PCHiC) enable the study of chromatin loops such as those formed between promoters and enhancers or CTCF/cohesin binding sites. An important aspect of 4C/Capture-C/PCHiC analyses is the reliable identification of chromatin loops, preferably not based on visual inspection of a DNA contact profile, but on reproducible statistical analysis that robustly scores interaction peaks in the non-uniform contact background. Here, we present peakC, an R package for the analysis of 4C/Capture-C/PCHiC data. We generated 4C data for 13 viewpoints in two tissues in at least triplicate to test our methods. We developed a non-parametric peak caller based on rank-products. Sampling analysis shows that not read depth but template quality is the most important determinant of success in 4C experiments. By performing peak calling on single experiments we show that the peak calling results are similar to the replicate experiments, but that false positive rates are significantly reduced by performing replicates. Our software is user-friendly and enables robust peak calling for one-vs-all chromosome capture experiments. peakC is available at: https://github.com/deWitLab/peakC.

Transcription factors and chromatin remodeling proteins control the transcriptional variability for ESC lineage commitment. During ESC differentiation, chromatin modifiers are recruited to the regulatory regions by transcription factors, thereby activating the lineage-specific genes or silencing the transcription of active ESC genes. However, the underlying mechanisms that link transcription factors to exit from pluripotency are yet to be identified. In this study, we show that the Ctbp2-interacting zinc finger proteins, Zfp217 and Zfp516, function as linkers for the chromatin regulators during ESC differentiation. CRISPR-Cas9-mediated knock-outs of both Zfp217 and Zfp516 in ESCs prevent the exit from pluripotency. Both zinc finger proteins regulate the Ctbp2-mediated recruitment of the NuRD complex and polycomb repressive complex 2 (PRC2) to active ESC genes, subsequently switching the H3K27ac to H3K27me3 during ESC differentiation for active gene silencing. We therefore suggest that some zinc finger proteins orchestrate to control the concise epigenetic states on active ESC genes during differentiation, resulting in natural lineage commitment.

Cyclin-dependent kinase 1 (Cdk1) is indispensable for embryonic stem cell (ESC) maintenance and embryo development. Even though some reports have described a connection between Cdk1 and Oct4, there is no evidence that Cdk1 activity is directly linked to the ESC pluripotency transcription program. We recently reported that Aurkb/PP1-mediated Oct4 resetting is important to cell cycle maintenance and pluripotency in mouse ESCs (mESCs). In this study, we show that Cdk1 is an upstream regulator of the Oct4 phosphorylation state during cell cycle progression, and it coordinates the chromatin associated state of Oct4 for pluripotency-related gene expression within the cell cycle. Upon entry into mitosis, Aurkb in the chromosome passenger complex becomes fully activated and PP1 activity is inhibited downstream of Cdk1 activation, leading to sustaining Oct4(S229) phosphorylation and dissociation of Oct4 from chromatin during the mitotic phase. Cdk1 inhibition at the mitotic phase abnormally results in Oct4 dephosphorylation, chromosome decondensation and chromatin association of Oct4, even in replicated chromosome. Our study results suggest a molecular mechanism by which Cdk1 directly links the cell cycle to the pluripotency transcription program in mESCs.

Sin3a is a core component of histone-deacetylation-activity-associated transcriptional repressor complex, playing important roles in early embryo development. Here, we reported that down-regulation of Sin3a led to the loss of embryonic stem cell (ESC) self-renewal and skewed differentiation into mesendoderm lineage. We found that Sin3a functioned as a transcriptional coactivator of the critical Nodal antagonist Lefty1 through interacting with Tet1 to de-methylate the Lefty1 promoter. Further studies showed that two amino acid residues (Phe147, Phe182) in the PAH1 domain of Sin3a are essential for Sin3a-Tet1 interaction and its activity in regulating pluripotency. Furthermore, genome-wide analyses of Sin3a, Tet1 and Pol II ChIP-seq and of 5mC MeDIP-seq revealed that Sin3a acted with Tet1 to facilitate the transcription of a set of their co-target genes. These results link Sin3a to epigenetic DNA modifications in transcriptional activation and have implications for understanding mechanisms underlying versatile functions of Sin3a in mouse ESCs.

The new web resource EviNet provides an easily run interface to network enrichment analysis for exploration of novel, experimentally defined gene sets. The major advantages of this analysis are (i) applicability to any genes found in the global network rather than only to those with pathway/ontology term annotations, (ii) ability to connect genes via different molecular mechanisms rather than within one high-throughput platform, and (iii) statistical power sufficient to detect enrichment of very small sets, down to individual genes. The users' gene sets are either defined prior to upload or derived interactively from an uploaded file by differential expression criteria. The pathways and networks used in the analysis can be chosen from the collection menu. The calculation is typically done within seconds or minutes and the stable URL is provided immediately. The results are presented in both visual (network graphs) and tabular formats using jQuery libraries. Uploaded data and analysis results are kept in separated project directories not accessible by other users. EviNet is available at https://www.evinet.org/.

Embryonic stem cells (ESCs) and meiosis are featured by relatively higher frequent homologous recombination associated with DNA double strand breaks (DSB) repair. Here, we show that Pold3 plays important roles in DSB repair, telomere maintenance and genomic stability of both ESCs and spermatocytes in mice. By attempting to generate Pold3 deficient mice using CRISPR/Cas9 or transcription activator-like effector nucleases, we show that complete loss of Pold3 (Pold3-/-) resulted in early embryonic lethality at E6.5. Rapid DNA damage response and massive apoptosis occurred in both outgrowths of Pold3-null (Pold3-/-) blastocysts and Pold3 inducible knockout (iKO) ESCs. While Pold3-/- ESCs were not achievable, Pold3 iKO led to increased DNA damage response, telomere loss and chromosome breaks accompanied by extended S phase. Meanwhile, loss of Pold3 resulted in replicative stress, micronucleation and aneuploidy. Also, DNA repair was impaired in Pold3+/- or Pold3 knockdown ESCs. Moreover, Pold3 mediates DNA replication and repair by regulating 53BP1, RIF1, ATR and ATM pathways. Furthermore, spermatocytes of Pold3 haploinsufficient (Pold3+/-) mice with increasing age displayed impaired DSB repair, telomere shortening and loss, and chromosome breaks, like Pold3 iKO ESCs. These data suggest that Pold3 maintains telomere integrity and genomic stability of both ESCs and meiosis by suppressing replicative stress.

Transposable elements, including endogenous retroviruses (ERVs), constitute a large fraction of the mammalian genome. They are transcriptionally silenced during early development to protect genome integrity and aberrant transcription. However, the mechanisms that control their repression are not fully understood. To systematically study ERV repression, we carried out an RNAi screen in mouse embryonic stem cells (ESCs) and identified a list of novel regulators. Among them, Rif1 displays the strongest effect. Rif1 depletion by RNAi or gene deletion led to increased transcription and increased chromatin accessibility at ERV regions and their neighboring genes. This transcriptional de-repression becomes more severe when DNA methylation is lost. On the mechanistic level, Rif1 directly occupies ERVs and is required for repressive histone mark H3K9me3 and H3K27me3 assembly and DNA methylation. It interacts with histone methyltransferases and facilitates their recruitment to ERV regions. Importantly, Rif1 represses ERVs in human ESCs as well, and the evolutionally-conserved HEAT-like domain is essential for its function. Finally, Rif1 acts as a barrier during somatic cell reprogramming, and its depletion significantly enhances reprogramming efficiency. Together, our study uncovered many previously uncharacterized repressors of ERVs, and defined an essential role of Rif1 in the epigenetic defense against ERV activation., (Published by Oxford University Press on behalf of Nucleic Acids Research 2017.)

Alternative splicing (AS) is a genetically and epigenetically regulated pre-mRNA processing to increase transcriptome and proteome diversity. Comprehensively decoding these regulatory mechanisms holds promise in getting deeper insights into a variety of biological contexts involving in AS, such as development and diseases. We assembled splicing (epi)genetic code, DeepCode, for human embryonic stem cell (hESC) differentiation by integrating heterogeneous features of genomic sequences, 16 histone modifications with a multi-label deep neural network. With the advantages of epigenetic features, DeepCode significantly improves the performance in predicting the splicing patterns and their changes during hESC differentiation. Meanwhile, DeepCode reveals the superiority of epigenomic features and their dominant roles in decoding AS patterns, highlighting the necessity of including the epigenetic properties when assembling a more comprehensive splicing code. Moreover, DeepCode allows the robust predictions across cell lineages and datasets. Especially, we identified a putative H3K36me3-regulated AS event leading to a nonsense-mediated mRNA decay of BARD1. Reduced BARD1 expression results in the attenuation of ATM/ATR signalling activities and further the hESC differentiation. These results suggest a novel candidate mechanism linking histone modifications to hESC fate decision. In addition, when trained in different contexts, DeepCode can be expanded to a variety of biological and biomedical fields., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Histone variants and their chaperones are key regulators of eukaryotic transcription, and are critical for normal development. The histone variant H3.3 has been shown to play important roles in pluripotency and differentiation, and although its genome-wide patterns have been investigated, little is known about the role of its dynamic turnover in transcriptional regulation. To elucidate the role of H3.3 dynamics in embryonic stem cell (ESC) biology, we generated mouse ESC lines carrying a single copy of a doxycycline (Dox)-inducible HA-tagged version of H3.3 and monitored the rate of H3.3 incorporation by ChIP-seq at varying time points following Dox induction, before and after RA-induced differentiation. Comparing H3.3 turnover profiles in ESCs and RA-treated cells, we identified a hyperdynamic H3.3-containing nucleosome at the -1 position in promoters of genes expressed in ESCs. This dynamic nucleosome is restricted and shifted downstream into the +1 position following differentiation. We suggest that histone turnover dynamics provides an additional mechanism involved in expression regulation, and that a hyperdynamic -1 nucleosome marks promoters in ESCs. Our data provide evidence for regional regulation of H3.3 turnover in ESC promoters, and calls for testing, in high resolution, the dynamic behavior of additional histone variants and other structural chromatin proteins., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

We report here a robust, tunable, and reversible transcription control system for endogenous genes. The REMOTE-control system (Reversible Manipulation of Transcription at Endogenous loci) employs enhanced lac repression and tet activation systems. With this approach, we show in mouse embryonic stem cells that endogenous Dnmt1 gene transcription could be up- or downregulated in a tunable, inducible, and reversible manner across nearly two orders of magnitude. Transcriptional repression of Dnmt1 by REMOTE-control was potent enough to cause embryonic lethality in mice, reminiscent of a genetic knockout of Dnmt1 and could substantially suppress intestinal polyp formation when applied to an ApcMin model. Binding by the enhanced lac repressor was sufficiently tight to allow strong attenuation of transcriptional elongation, even at operators located many kilobases downstream of the transcription start site and to produce invariably tight repression of all of the strong viral/mammalian promoters tested. Our approach of targeting tet transcriptional activators to the endogenous Dnmt1 promoter resulted in robust upregulation of this highly expressed housekeeping gene. Our system provides exquisite control of the level, timing, and cell-type specificity of endogenous gene expression, and the potency and versatility of the system will enable high resolution in vivo functional analyses., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Transposable elements (TEs) compose about 40% of the murine genome. Retrotransposition of active TEs such as LINE-1 (L1) tremendously impacts genetic diversification and genome stability. Therefore, transcription and transposition activities of retrotransposons are tightly controlled. Here, we show that the Krüppel-like zinc finger protein Zfp281 directly binds and suppresses a subset of retrotransposons, including the active young L1 repeat elements, in mouse embryonic stem (ES) cells. In addition, we find that Zfp281-regulated L1s are highly enriched for 5-hydroxymethylcytosine (5hmC) and H3K4me3. The COMPASS-like H3K4 methyltransferase Mll2 is the major H3K4me3 methylase at the Zfp281-regulated L1s and required for their proper expression. Our studies also reveal that Zfp281 functions partially through recruiting the L1 regulators DNA hydroxymethylase Tet1 and Sin3A, and restricting Mll2 at these active L1s, leading to their balanced expression. In summary, our data indicate an instrumental role of Zfp281 in suppressing the young active L1s in mouse ES cells., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Mouse embryonic stem (ES) cells are a popular model system to study biological processes, though uncovering recessive phenotypes requires inactivating both alleles. Building upon resources from the International Knockout Mouse Consortium (IKMC), we developed a targeting vector for second allele inactivation in conditional-ready IKMC 'knockout-first' ES cell lines. We applied our technology to several epigenetic regulators, recovering bi-allelic targeted clones with a high efficiency of 60% and used Flp recombinase to restore expression in two null cell lines to demonstrate how our system confirms causality through mutant phenotype reversion. We designed our strategy to select against re-targeting the 'knockout-first' allele and identify essential genes in ES cells, including the histone methyltransferase Setdb1. For confirmation, we exploited the flexibility of our system, enabling tamoxifen inducible conditional gene ablation while controlling for genetic background and tamoxifen effects. Setdb1 ablated ES cells exhibit severe growth inhibition, which is not rescued by exogenous Nanog expression or culturing in naive pluripotency '2i' media, suggesting that the self-renewal defect is mediated through pluripotency network independent pathways. Our strategy to generate null mutant mouse ES cells is applicable to thousands of genes and repurposes existing IKMC Intermediate Vectors., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Telomerase is highly expressed in cancer and embryonic stem cells (ESCs) and implicated in controlling genome integrity, cancer formation and stemness. Previous studies identified that Krüppel-like transcription factor 4 (KLF4) activates telomerase reverse transcriptase (TERT) expression and contributes to the maintenance of self-renewal in ESCs. However, little is known about how KLF4 regulates TERT expression. Here, we discover poly(ADP-ribose) polymerase 1 (PARP1) as a novel KLF4-interacting partner. Knockdown of PARP1 reduces TERT expression and telomerase activity not only in cancer cells, but also in human and mouse ESCs. Recruitment of KLF4 to TERT promoter is reduced in PARP1-suppressed cells. The poly(ADP-ribose) polymerase activity is dispensable, while the oligo(ADP-ribose) polymerase activity is required for the PARP1- and KLF4-mediated TERT activation. Repression of Parp1 in mouse ESCs decreases expression of pluripotent markers and induces differentiation. These results suggest that PARP1 recruits KLF4 to activate telomerase expression and stem cell pluripotency, indicating a positive regulatory role of the PARP1-KLF4 complex in telomerase expression in cancer and stem cells., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Faithful expression of transgenes in cell cultures and mice is often challenged by locus dependent epigenetic silencing. We investigated silencing of Tet-controlled expression cassettes within the mouse ROSA26 locus. We observed pronounced DNA methylation of the Tet promoter concomitant with loss of expression in mES cells as well as in differentiated cells and transgenic animals. Strikingly, the ROSA26 promoter remains active and methylation free indicating that this silencing mechanism specifically affects the transgene, but does not spread to the host's chromosomal neighborhood. To reactivate Tet cassettes a synthetic fusion protein was constructed and expressed in silenced cells. This protein includes the enzymatic domains of ten eleven translocation methylcytosine dioxygenase 1 (TET-1) as well as the Tet repressor DNA binding domain. Expression of the synthetic fusion protein and Doxycycline treatment allowed targeted demethylation of the Tet promoter in the ROSA26 locus and in another genomic site, rescuing transgene expression in cells and transgenic mice. Thus, inducible, reversible and site-specific epigenetic modulation is a promising strategy for reactivation of silenced transgene expression, independent of the integration site., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

TET1 oxidizes methylated cytosine into 5-hydroxymethylcytosine (5hmC), resulting in regulation of DNA methylation and gene expression. Full length TET1 (TET1FL) has a CXXC domain that binds to unmethylated CpG islands (CGIs). This CXXC domain allows TET1 to protect CGIs from aberrant methylation, but it also limits its ability to regulate genes outside of CGIs. Here, we report a novel isoform of TET1 (TET1ALT) that has a unique transcription start site from an alternate promoter in intron 2, yielding a protein with a unique translation start site. Importantly, TET1ALT lacks the CXXC domain but retains the catalytic domain. TET1ALT is repressed in embryonic stem cells (ESCs) but becomes activated in embryonic and adult tissues while TET1FL is expressed in ESCs, but repressed in adult tissues. Overexpression of TET1ALT shows production of 5hmC with distinct (and weaker) effects on DNA methylation or gene expression when compared to TET1FL. TET1ALT is aberrantly activated in multiple cancer types including breast, uterine and glioblastoma, and TET1 activation is associated with a worse overall survival in breast, uterine and ovarian cancers. Our data suggest that the predominantly activated isoform of TET1 in cancer cells does not protect from CGI methylation and likely mediates dynamic site-specific demethylation outside of CGIs., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Designer nucleases like CRISPR/Cas9 enable fluent site-directed damage or small mutations in many genomes. Strategies for their use to achieve more complex tasks like regional exchanges for gene humanization or the establishment of conditional alleles are still emerging. To optimize Cas9-assisted targeting, we measured the relationship between targeting frequency and homology length in targeting constructs using a hypoxanthine-guanine phosphoribosyl-transferase assay in mouse embryonic stem cells. Targeting frequency with supercoiled plasmids improved steeply up to 2 kb total homology and continued to increase with even longer homology arms, thereby implying that Cas9-assisted targeting efficiencies can be improved using homology arms of 1 kb or greater. To humanize the Kmt2d gene, we built a hybrid mouse/human targeting construct in a bacterial artificial chromosome by recombineering. To simplify the possible outcomes, we employed a single Cas9 cleavage strategy and best achieved the intended 42 kb regional exchange with a targeting construct including a very long homology arm to recombine ∼42 kb away from the cleavage site. We recommend the use of long homology arm targeting constructs for accurate and efficient complex genome engineering, particularly when combined with the simplifying advantages of using just one Cas9 cleavage at the genome target site., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)