Genome-wide nucleosome specificity and function of chromatin remodellers in ES cells

Structural insights into the human NuA4/TIP60 acetyltransferase and chromatin remodeling complex

The human nucleosome acetyltransferase of histone H4 (NuA4)/Tat-interactive protein, 60 kilodalton (TIP60) coactivator complex, a fusion of the yeast switch/sucrose nonfermentable related 1 (SWR1) and NuA4 complexes, both incorporates the histone variant H2A.Z into nucleosomes and acetylates histones H4, H2A, and H2A.Z to regulate gene expression and maintain genome stability. Our cryo-electron microscopy studies show that, within the NuA4/TIP60 complex, the E1A binding protein P400 (EP400) subunit serves as a scaffold holding the different functional modules in specific positions, creating a distinct arrangement of the actin-related protein (ARP) module. EP400 interacts with the transformation/transcription domain-associated protein (TRRAP) subunit by using a footprint that overlaps with that of the Spt-Ada-Gcn5 acetyltransferase (SAGA) complex, preventing the formation of a hybrid complex. Loss of the TRRAP subunit leads to mislocalization of NuA4/TIP60, resulting in the redistribution of H2A.Z and its acetylation across the genome, emphasizing the dual functionality of NuA4/TIP60 as a single macromolecular assembly.

Multifunctional histone variants in genome function

Histones are integral components of eukaryotic chromatin that have a pivotal role in the organization and function of the genome. The dynamic regulation of chromatin involves the incorporation of histone variants, which can dramatically alter its structural and functional properties. Contrary to an earlier view that limited individual histone variants to specific genomic functions, new insights have revealed that histone variants exert multifaceted roles involving all aspects of genome function, from governing patterns of gene expression at precise genomic loci to participating in genome replication, repair and maintenance. This conceptual change has led to a new understanding of the intricate interplay between chromatin and DNA-dependent processes and how this connection translates into normal and abnormal cellular functions.

guidedNOMe-seq quantifies chromatin states at single allele resolution for hundreds of custom regions in parallel

Since the introduction of next generation sequencing technologies, the field of epigenomics has evolved rapidly. However, most commonly used assays are enrichment-based methods and thus only semi-quantitative. Nucleosome occupancy and methylome sequencing (NOMe-seq) allows for quantitative inference of chromatin states with single locus resolution, but this requires high sequencing depth and is therefore prohibitively expensive to routinely apply to organisms with large genomes. To overcome this limitation, we introduce guidedNOMe-seq, where we combine NOMe profiling with large scale sgRNA synthesis and Cas9-mediated region-of-interest (ROI) liberation. To facilitate quantitative comparisons between multiple samples, we additionally develop an R package to standardize differential analysis of any type of NOMe-seq data. We extensively benchmark guidedNOMe-seq in a proof-of-concept study, dissecting the interplay of ChAHP and CTCF on chromatin. In summary we present a cost-effective, scalable, and customizable target enrichment extension to the existing NOMe-seq protocol allowing genome-scale quantification of nucleosome occupancy and transcription factor binding at single allele resolution.

ChAHP2 and ChAHP control diverse retrotransposons by complementary activities

Retrotransposon control in mammals is an intricate process that is effectuated by a broad network of chromatin regulatory pathways. We previously discovered ChAHP, a protein complex with repressive activity against short interspersed element (SINE) retrotransposons that is composed of the transcription factor ADNP, chromatin remodeler CHD4, and HP1 proteins. Here we identify ChAHP2, a protein complex homologous to ChAHP, in which ADNP is replaced by ADNP2. ChAHP2 is predominantly targeted to endogenous retroviruses (ERVs) and long interspersed elements (LINEs) via HP1β-mediated binding of H3K9 trimethylated histones. We further demonstrate that ChAHP also binds these elements in a manner mechanistically equivalent to that of ChAHP2 and distinct from DNA sequence-specific recruitment at SINEs. Genetic ablation of ADNP2 alleviates ERV and LINE1 repression, which is synthetically exacerbated by additional depletion of ADNP. Together, our results reveal that the ChAHP and ChAHP2 complexes function to control both nonautonomous and autonomous retrotransposons by complementary activities, further adding to the complexity of mammalian transposon control.

CARM1 hypermethylates the NuRD chromatin remodeling complex to promote cell cycle gene expression and breast cancer development

Protein arginine methyltransferase CARM1 has been shown to methylate a large number of non-histone proteins, and play important roles in gene transcriptional activation, cell cycle progress, and tumorigenesis. However, the critical substrates through which CARM1 exerts its functions remain to be fully characterized. Here, we reported that CARM1 directly interacts with the GATAD2A/2B subunit in the nucleosome remodeling and deacetylase (NuRD) complex, expanding the activities of NuRD to include protein arginine methylation. CARM1 and NuRD bind and activate a large cohort of genes with implications in cell cycle control to facilitate the G1 to S phase transition. This gene activation process requires CARM1 to hypermethylate GATAD2A/2B at a cluster of arginines, which is critical for the recruitment of the NuRD complex. The clinical significance of this gene activation mechanism is underscored by the high expression of CARM1 and NuRD in breast cancers, and the fact that knockdown CARM1 and NuRD inhibits cancer cell growth in vitro and tumorigenesis in vivo. Targeting CARM1-mediated GATAD2A/2B methylation with CARM1 specific inhibitors potently inhibit breast cancer cell growth in vitro and tumorigenesis in vivo. These findings reveal a gene activation program that requires arginine methylation established by CARM1 on a key chromatin remodeler, and targeting such methylation might represent a promising therapeutic avenue in the clinic.

cBAF generates subnucleosomes that expand OCT4 binding and function beyond DNA motifs at enhancers

The canonical BRG/BRM-associated factor (cBAF) complex is essential for chromatin opening at enhancers in mammalian cells. However, the nature of the open chromatin remains unclear. Here, we show that, in addition to producing histone-free DNA, cBAF generates stable hemisome-like subnucleosomal particles containing the four core histones associated with 50-80 bp of DNA. Our genome-wide analysis indicates that cBAF makes these particles by targeting and splitting fragile nucleosomes. In mouse embryonic stem cells, these subnucleosomes become an in vivo binding substrate for the master transcription factor OCT4 independently of the presence of OCT4 DNA motifs. At enhancers, the OCT4-subnucleosome interaction increases OCT4 occupancy and amplifies the genomic interval bound by OCT4 by up to one order of magnitude compared to the region occupied on histone-free DNA. We propose that cBAF-dependent subnucleosomes orchestrate a molecular mechanism that projects OCT4 function in chromatin opening beyond its DNA motifs.

Identification of a molecular network regulated by multiple ASD high risk genes

Genetic sequencing has identified high-confidence ASD risk genes with loss-of-function mutations. How the haploinsufficiency of distinct ASD risk genes causes ASD remains to be elucidated. In this study, we examined the role of four top-ranking ASD risk genes, ADNP, KDM6B, CHD2, and MED13, in gene expression regulation. ChIP-seq analysis reveals that gene targets with the binding of these ASD risk genes at promoters are enriched in RNA processing and DNA repair. Many of these targets are found in ASD gene database (SFARI), and are involved in transcription regulation and chromatin remodeling. Common gene targets of these ASD risk genes include a network of high confidence ASD genes associated with gene expression regulation, such as CTNNB1 and SMARCA4. We further directly examined the transcriptional impact of the deficiency of these ASD risk genes. Our mRNA profiling with qPCR assays in cells with the knockdown of Adnp, Kdm6b, Chd2 or Med13 has revealed an intricate pattern of their cross-regulation, as well as their influence on the expression of other ASD genes. In addition, some synaptic genes, such as Snap25 and Nrxn1, are strongly regulated by deficiency of the four ASD risk genes, which could be through the direct binding at promoters or indirectly through the targets like Ctnnb1 or Smarca4. The identification of convergent and divergent gene targets that are regulated by multiple ASD risk genes will help to understand the molecular mechanisms underlying common and unique phenotypes associated with haploinsufficiency of ASD-associated genes.

Targeting dependency on a paralog pair of CBP/p300 against de-repression of KREMEN2 in SMARCB1-deficient cancers

SMARCB1, a subunit of the SWI/SNF chromatin remodeling complex, is the causative gene of rhabdoid tumors and epithelioid sarcomas. Here, we identify a paralog pair of CBP and p300 as a synthetic lethal target in SMARCB1-deficient cancers by using a dual siRNA screening method based on the "simultaneous inhibition of a paralog pair" concept. Treatment with CBP/p300 dual inhibitors suppresses growth of cell lines and tumor xenografts derived from SMARCB1-deficient cells but not from SMARCB1-proficient cells. SMARCB1-containing SWI/SNF complexes localize with H3K27me3 and its methyltransferase EZH2 at the promotor region of the KREMEN2 locus, resulting in transcriptional downregulation of KREMEN2. By contrast, SMARCB1 deficiency leads to localization of H3K27ac, and recruitment of its acetyltransferases CBP and p300, at the KREMEN2 locus, resulting in transcriptional upregulation of KREMEN2, which cooperates with the SMARCA1 chromatin remodeling complex. Simultaneous inhibition of CBP/p300 leads to transcriptional downregulation of KREMEN2, followed by apoptosis induction via monomerization of KREMEN1 due to a failure to interact with KREMEN2, which suppresses anti-apoptotic signaling pathways. Taken together, our findings indicate that simultaneous inhibitors of CBP/p300 could be promising therapeutic agents for SMARCB1-deficient cancers.

Acute depletion of BRG1 reveals its primary function as an activator of transcription

The mammalian SWI/SNF-like BAF complexes play critical roles during animal development and pathological conditions. Previous gene deletion studies and characterization of human gene mutations implicate that the complexes both repress and activate a large number of genes. However, the direct function of the complexes in cells remains largely unclear due to the relatively long-term nature of gene deletion or natural mutation. Here we generate a mouse line by knocking in the auxin-inducible degron tag (AID) to the Smarca4 gene, which encodes BRG1, the essential ATPase subunit of the BAF complexes. We show that the tagged BRG1 can be efficiently depleted by osTIR1 expression and auxin treatment for 6 to 10 h in CD4 + T cells, hepatocytes, and fibroblasts isolated from the knock-in mice. The acute depletion of BRG1 leads to decreases in nascent RNAs and RNA polymerase II binding at a large number of genes, which are positively correlated with the loss of BRG1. Further, these changes are correlated with diminished accessibility at DNase I Hypersensitive Sites (DHSs) and p300 binding. The acute BRG1 depletion results in three major patterns of nucleosome shifts leading to narrower nucleosome spacing surrounding transcription factor motifs and at enhancers and transcription start sites (TSSs), which are correlated with loss of BRG1, decreased chromatin accessibility and decreased nascent RNAs. Acute depletion of BRG1 severely compromises the Trichostatin A (TSA) -induced histone acetylation, suggesting a substantial interplay between the chromatin remodeling activity of BRG1 and histone acetylation. Our data suggest BRG1 mainly plays a direct positive role in chromatin accessibility, RNAPII binding, and nascent RNA production by regulating nucleosome positioning and facilitating transcription factor binding to their target sites.

Epigenetic reprogramming in the transition from pluripotency to totipotency

Mammalian development commences with the zygote, which can differentiate into both embryonic and extraembryonic tissues, a capability known as totipotency. Only the zygote and embryos around zygotic genome activation (ZGA) (two-cell embryo stage in mice and eight-cell embryo in humans) are totipotent cells. Epigenetic modifications undergo extremely extensive changes during the acquisition of totipotency and subsequent development of differentiation. However, the underlying molecular mechanisms remain elusive. Recently, the discovery of mouse two-cell embryo-like cells, human eight-cell embryo-like cells, extended pluripotent stem cells and totipotent-like stem cells with extra-embryonic developmental potential has greatly expanded our understanding of totipotency. Experiments with these in vitro models have led to insights into epigenetic changes in the reprogramming of pluri-to-totipotency, which have informed the exploration of preimplantation development. In this review, we highlight the recent findings in understanding the mechanisms of epigenetic remodeling during totipotency capture, including RNA splicing, DNA methylation, chromatin configuration, histone modifications, and nuclear organization.

Chromatin Modifier EP400 Regulates Oocyte Quality and Zygotic Genome Activation in Mice

Epigenetic modifiers that accumulate in oocytes, play a crucial role in steering the developmental program of cleavage embryos and initiating life. However, the identification of key maternal epigenetic regulators remains elusive. In the findings, the essential role of maternal Ep400, a chaperone for H3.3, in oocyte quality and early embryo development in mice is highlighted. Depletion of Ep400 in oocytes resulted in a decline in oocyte quality and abnormalities in fertilization. Preimplantation embryos lacking maternal Ep400 exhibited reduced major zygotic genome activation (ZGA) and experienced developmental arrest at the 2-to-4-cell stage. The study shows that EP400 forms protein complex with NFYA, occupies promoters of major ZGA genes, modulates H3.3 distribution between euchromatin and heterochromatin, promotes transcription elongation, activates the expression of genes regulating mitochondrial functions, and facilitates the expression of rate-limiting enzymes of the TCA cycle. This intricate process driven by Ep400 ensures the proper execution of the developmental program, emphasizing its critical role in maternal-to-embryonic transition.

Widespread impact of nucleosome remodelers on transcription at cis-regulatory elements

Nucleosome remodeling complexes and other regulatory factors work in concert to build a chromatin environment that directs the expression of a distinct set of genes in each cell using cis-regulatory elements (CREs), such as promoters and enhancers, that drive transcription of both mRNAs and CRE-associated non-coding RNAs (ncRNAs). Two classes of CRE-associated ncRNAs include upstream antisense RNAs (uaRNAs), which are transcribed divergently from a shared mRNA promoter, and enhancer RNAs (eRNAs), which are transcribed bidirectionally from active enhancers. The complicated network of CRE regulation by nucleosome remodelers remains only partially explored, with a focus on a select, limited number of remodelers. We endeavored to elucidate a remodeler-based regulatory network governing CRE-associated transcription (mRNA, eRNA, and uaRNA) in murine embryonic stem (ES) cells to test the hypothesis that many SNF2-family nucleosome remodelers collaborate to regulate the coding and non-coding transcriptome via alteration of underlying nucleosome architecture. Using depletion followed by transient transcriptome sequencing (TT-seq), we identified thousands of misregulated mRNAs and CRE-associated ncRNAs across the remodelers examined, identifying novel contributions by understudied remodelers in the regulation of coding and noncoding transcription. Our findings suggest that mRNA and eRNA transcription are coordinately co-regulated, while mRNA and uaRNAs sharing a common promoter are independently regulated. Subsequent mechanistic studies suggest that while remodelers SRCAP and CHD8 modulate transcription through classical mechanisms such as transcription factors and histone variants, a broad set of remodelers including SMARCAL1 indirectly contribute to transcriptional regulation through maintenance of genomic stability and proper Integrator complex localization. This study systematically examines the contribution of SNF2-remodelers to the CRE-associated transcriptome, identifying at least two classes for remodeler action.

Genome-wide ATAC-see screening identifies TFDP1 as a modulator of global chromatin accessibility

Chromatin accessibility is a hallmark of active regulatory regions and is functionally linked to transcriptional networks and cell identity. However, the molecular mechanisms and networks that govern chromatin accessibility have not been thoroughly studied. Here we conducted a genome-wide CRISPR screening combined with an optimized ATAC-see protocol to identify genes that modulate global chromatin accessibility. In addition to known chromatin regulators like CREBBP and EP400, we discovered a number of previously unrecognized proteins that modulate chromatin accessibility, including TFDP1, HNRNPU, EIF3D and THAP11 belonging to diverse biological pathways. ATAC-seq analysis upon their knockouts revealed their distinct and specific effects on chromatin accessibility. Remarkably, we found that TFDP1, a transcription factor, modulates global chromatin accessibility through transcriptional regulation of canonical histones. In addition, our findings highlight the manipulation of chromatin accessibility as an approach to enhance various cell engineering applications, including genome editing and induced pluripotent stem cell reprogramming.

Dancr-BRG1 regulates Nfatc1 transcription and Pgc1β-dependent metabolic shifts in osteoclastogenesis

Long non-coding RNA (lncRNA) serves as a vital regulator of bone metabolism, but its role in pathologically overactive osteoclast differentiation remains elusive. Here, we identify lncRNA Dancr (Differentiation Antagonizing Non-protein Coding RNA) as a critical suppressor of osteoclastogenesis and bone resorption, which is down-regulated in response to estrogen deficiency. Global or osteoclast-specific Dancr Knockout mice display significant trabecular bone deterioration and enhanced osteoclast activity, but minimal alteration of bone formation. Moreover, the bone-targeted delivery of Dancr by Adeno-associated viral remarkably attenuates ovariectomy-induced osteopenia in mice. Mechanistically, Dancr establishes a direct interaction with Brahma-related gene 1 to prevent its binding and preserve H3K27me3 enrichment at the nuclear factor of activated T cells 1 and proliferator-activated receptor gamma coactivator 1-beta promoters, thereby maintaining appropriate expression of osteoclastic genes and metabolic programs during osteoclastogenesis. These results demonstrate that Dancr is a key molecule maintaining proper osteoclast differentiation and bone homeostasis under physiological conditions, and Dancr overexpression constitutes a potential strategy for treating osteoporosis.

The ncBAF Complex Regulates Transcription in AML Through H3K27ac Sensing by BRD9

The non-canonical BAF complex (ncBAF) subunit BRD9 is essential for acute myeloid leukemia (AML) cell viability but has an unclear role in leukemogenesis. Because BRD9 is required for ncBAF complex assembly through its DUF3512 domain, precise bromodomain inhibition is necessary to parse the role of BRD9 as a transcriptional regulator from that of a scaffolding protein. To understand the role of BRD9 bromodomain function in regulating AML, we selected a panel of five AML cell lines with distinct driver mutations, disease classifications, and genomic aberrations and subjected these cells to short-term BRD9 bromodomain inhibition. We examined the bromodomain-dependent growth of these cell lines, identifying a dependency in AML cell lines but not HEK293T cells. To define a mechanism through which BRD9 maintains AML cell survival, we examined nascent transcription, chromatin accessibility, and ncBAF complex binding genome-wide after bromodomain inhibition. We identified extensive regulation of transcription by BRD9 bromodomain activity, including repression of myeloid maturation factors and tumor suppressor genes, while standard AML chemotherapy targets were repressed by inhibition of the BRD9 bromodomain. BRD9 bromodomain activity maintained accessible chromatin at both gene promoters and gene-distal putative enhancer regions, in a manner that qualitatively correlated with enrichment of BRD9 binding. Furthermore, we identified reduced chromatin accessibility at GATA, ETS, and AP-1 motifs and increased chromatin accessibility at SNAIL-, HIC-, and TP53-recognized motifs after BRD9 inhibition. These data suggest a role for BRD9 in regulating AML cell differentiation through modulation of accessibility at hematopoietic transcription factor binding sites.

SIGNIFICANCE

The bromodomain-containing protein BRD9 is essential for AML cell viability, but it is unclear whether this requirement is due to the protein's role as an epigenetic reader. We inhibited this activity and identified altered gene-distal chromatin regulation and transcription consistent with a more mature myeloid cell state.

The Nucleosome Remodelling and Deacetylation complex coordinates the transcriptional response to lineage commitment in pluripotent cells

As cells exit the pluripotent state and begin to commit to a specific lineage they must activate genes appropriate for that lineage while silencing genes associated with pluripotency and preventing activation of lineage-inappropriate genes. The Nucleosome Remodelling and Deacetylation (NuRD) complex is essential for pluripotent cells to successfully undergo lineage commitment. NuRD controls nucleosome density at regulatory sequences to facilitate transcriptional responses, and also has been shown to prevent unscheduled transcription (transcriptional noise) in undifferentiated pluripotent cells. How these activities combine to ensure cells engage a gene expression program suitable for successful lineage commitment has not been determined. Here, we show that NuRD is not required to silence all genes. Rather, it restricts expression of genes primed for activation upon exit from the pluripotent state, but maintains them in a transcriptionally permissive state in self-renewing conditions, which facilitates their subsequent activation upon exit from naïve pluripotency. We further show that NuRD coordinates gene expression changes, which acts to maintain a barrier between different stable states. Thus NuRD-mediated chromatin remodelling serves multiple functions, including reducing transcriptional noise, priming genes for activation and coordinating the transcriptional response to facilitate lineage commitment.

Chromatin gatekeeper and modifier CHD proteins in development, and in autism and other neurological disorders

Chromatin, a protein-DNA complex, is a dynamic structure that stores genetic information within the nucleus and responds to molecular/cellular changes in its structure, providing conditional access to the genetic machinery. ATP-dependent chromatin modifiers regulate access of transcription factors and RNA polymerases to DNA by either "opening" or "closing" the structure of chromatin, and its aberrant regulation leads to a variety of neurodevelopmental disorders. The chromodomain helicase DNA-binding (CHD) proteins are ATP-dependent chromatin modifiers involved in the organization of chromatin structure, act as gatekeepers of genomic access, and deposit histone variants required for gene regulation. In this review, we first discuss the structural and functional domains of the CHD proteins, and their binding sites, and phosphorylation, acetylation, and methylation sites. The conservation of important amino acids in SWItch/sucrose non-fermenting (SWI/SNF) domains, and their protein and mRNA tissue expression profiles are discussed. Next, we convey the important binding partners of CHD proteins, their protein complexes and activities, and their involvements in epigenetic regulation. We also show the ChIP-seq binding dynamics for CHD1, CHD2, CHD4, and CHD7 proteins at promoter regions of histone genes, as well as several genes that are critical for neurodevelopment. The role of CHD proteins in development is also discussed. Finally, this review provides information about CHD protein mutations reported in autism and neurodevelopmental disorders, and their pathogenicity. Overall, this review provides information on the progress of research into CHD proteins, their structural and functional domains, epigenetics, and their role in stem cell, development, and neurological disorders.

Global identification of SWI/SNF targets reveals compensation by EP400

Mammalian SWI/SNF chromatin remodeling complexes move and evict nucleosomes at gene promoters and enhancers to modulate DNA access. Although SWI/SNF subunits are commonly mutated in disease, therapeutic options are limited by our inability to predict SWI/SNF gene targets and conflicting studies on functional significance. Here, we leverage a fast-acting inhibitor of SWI/SNF remodeling to elucidate direct targets and effects of SWI/SNF. Blocking SWI/SNF activity causes a rapid and global loss of chromatin accessibility and transcription. Whereas repression persists at most enhancers, we uncover a compensatory role for the EP400/TIP60 remodeler, which reestablishes accessibility at most promoters during prolonged loss of SWI/SNF. Indeed, we observe synthetic lethality between EP400 and SWI/SNF in cancer cell lines and human cancer patient data. Our data define a set of molecular genomic features that accurately predict gene sensitivity to SWI/SNF inhibition in diverse cancer cell lines, thereby improving the therapeutic potential of SWI/SNF inhibitors.

Nucleosome density shapes kilobase-scale regulation by a mammalian chromatin remodeler

Nearly all essential nuclear processes act on DNA packaged into arrays of nucleosomes. However, our understanding of how these processes (for example, DNA replication, RNA transcription, chromatin extrusion and nucleosome remodeling) occur on individual chromatin arrays remains unresolved. Here, to address this deficit, we present SAMOSA-ChAAT: a massively multiplex single-molecule footprinting approach to map the primary structure of individual, reconstituted chromatin templates subject to virtually any chromatin-associated reaction. We apply this method to distinguish between competing models for chromatin remodeling by the essential imitation switch (ISWI) ATPase SNF2h: nucleosome-density-dependent spacing versus fixed-linker-length nucleosome clamping. First, we perform in vivo single-molecule nucleosome footprinting in murine embryonic stem cells, to discover that ISWI-catalyzed nucleosome spacing correlates with the underlying nucleosome density of specific epigenomic domains. To establish causality, we apply SAMOSA-ChAAT to quantify the activities of ISWI ATPase SNF2h and its parent complex ACF on reconstituted nucleosomal arrays of varying nucleosome density, at single-molecule resolution. We demonstrate that ISWI remodelers operate as density-dependent, length-sensing nucleosome sliders, whose ability to program DNA accessibility is dictated by single-molecule nucleosome density. We propose that the long-observed, context-specific regulatory effects of ISWI complexes can be explained in part by the sensing of nucleosome density within epigenomic domains. More generally, our approach promises molecule-precise views of the essential processes that shape nuclear physiology.

FACT regulates pluripotency through proximal and distal regulation of gene expression in murine embryonic stem cells

BACKGROUND

The FACT complex is a conserved histone chaperone with critical roles in transcription and histone deposition. FACT is essential in pluripotent and cancer cells, but otherwise dispensable for most mammalian cell types. FACT deletion or inhibition can block induction of pluripotent stem cells, yet the mechanism through which FACT regulates cell fate decisions remains unclear.

RESULTS

To explore the mechanism for FACT function, we generated AID-tagged murine embryonic cell lines for FACT subunit SPT16 and paired depletion with nascent transcription and chromatin accessibility analyses. We also analyzed SPT16 occupancy using CUT&RUN and found that SPT16 localizes to both promoter and enhancer elements, with a strong overlap in binding with OCT4, SOX2, and NANOG. Over a timecourse of SPT16 depletion, nucleosomes invade new loci, including promoters, regions bound by SPT16, OCT4, SOX2, and NANOG, and TSS-distal DNaseI hypersensitive sites. Simultaneously, transcription of Pou5f1 (encoding OCT4), Sox2, Nanog, and enhancer RNAs produced from these genes' associated enhancers are downregulated.

CONCLUSIONS

We propose that FACT maintains cellular pluripotency through a precise nucleosome-based regulatory mechanism for appropriate expression of both coding and non-coding transcripts associated with pluripotency.

BAZ2A-RNA mediated association with TOP2A and KDM1A represses genes implicated in prostate cancer

BAZ2A represses rRNA genes (rDNA) that are transcribed by RNA polymerase I. In prostate cancer (PCa), BAZ2A function goes beyond this role because it represses genes frequently silenced in metastatic disease. However, the mechanisms of this BAZ2A-mediated repression remain elusive. Here, we show that BAZ2A represses genes through its RNA-binding TAM domain using mechanisms differing from rDNA silencing. Although the TAM domain mediates BAZ2A recruitment to rDNA, in PCa, this is not required for BAZ2A association with target genes. Instead, the BAZ2A-TAM domain in association with RNA mediates the interaction with topoisomerase 2A (TOP2A) and histone demethylase KDM1A, whose expression positively correlates with BAZ2A levels in localized and metastatic PCa. TOP2A and KDM1A pharmacological inhibition up-regulate BAZ2A-repressed genes that are regulated by inactive enhancers bound by BAZ2A, whereas rRNA genes are not affected. Our findings showed a novel RNA-based mechanism of gene regulation in PCa. Furthermore, we determined that RNA-mediated interactions between BAZ2A and TOP2A and KDM1A repress genes critical to PCa and may prove to be useful to stratify prostate cancer risk and treatment in patients.

Mitotic bookmarking by SWI/SNF subunits

For cells to initiate and sustain a differentiated state, it is necessary that a 'memory' of this state is transmitted through mitosis to the daughter cells1-3. Mammalian switch/sucrose non-fermentable (SWI/SNF) complexes (also known as Brg1/Brg-associated factors, or BAF) control cell identity by modulating chromatin architecture to regulate gene expression4-7, but whether they participate in cell fate memory is unclear. Here we provide evidence that subunits of SWI/SNF act as mitotic bookmarks to safeguard cell identity during cell division. The SWI/SNF core subunits SMARCE1 and SMARCB1 are displaced from enhancers but are bound to promoters during mitosis, and we show that this binding is required for appropriate reactivation of bound genes after mitotic exit. Ablation of SMARCE1 during a single mitosis in mouse embryonic stem cells is sufficient to disrupt gene expression, impair the occupancy of several established bookmarks at a subset of their targets and cause aberrant neural differentiation. Thus, SWI/SNF subunit SMARCE1 has a mitotic bookmarking role and is essential for heritable epigenetic fidelity during transcriptional reprogramming.

The esBAF and ISWI nucleosome remodeling complexes influence occupancy of overlapping dinucleosomes and fragile nucleosomes in murine embryonic stem cells

BACKGROUND

Nucleosome remodeling factors regulate the occupancy and positioning of nucleosomes genome-wide through ATP-driven DNA translocation. While many nucleosomes are consistently well-positioned, some nucleosomes and alternative nucleosome structures are more sensitive to nuclease digestion or are transitory. Fragile nucleosomes are nucleosome structures that are sensitive to nuclease digestion and may be composed of either six or eight histone proteins, making these either hexasomes or octasomes. Overlapping dinucleosomes are composed of two merged nucleosomes, lacking one H2A:H2B dimer, creating a 14-mer wrapped by ~ 250 bp of DNA. In vitro studies of nucleosome remodeling suggest that the collision of adjacent nucleosomes by sliding stimulates formation of overlapping dinucleosomes.

RESULTS

To better understand how nucleosome remodeling factors regulate alternative nucleosome structures, we depleted murine embryonic stem cells of the transcripts encoding remodeler ATPases BRG1 or SNF2H, then performed MNase-seq. We used high- and low-MNase digestion to assess the effects of nucleosome remodeling factors on nuclease-sensitive or "fragile" nucleosome occupancy. In parallel we gel-extracted MNase-digested fragments to enrich for overlapping dinucleosomes. We recapitulate prior identification of fragile nucleosomes and overlapping dinucleosomes near transcription start sites, and identify enrichment of these features around gene-distal DNaseI hypersensitive sites, CTCF binding sites, and pluripotency factor binding sites. We find that BRG1 stimulates occupancy of fragile nucleosomes but restricts occupancy of overlapping dinucleosomes.

CONCLUSIONS

Overlapping dinucleosomes and fragile nucleosomes are prevalent within the ES cell genome, occurring at hotspots of gene regulation beyond their characterized existence at promoters. Although neither structure is fully dependent on either nucleosome remodeling factor, both fragile nucleosomes and overlapping dinucleosomes are affected by knockdown of BRG1, suggesting a role for the complex in creating or removing these structures.

Single-cell multi-scale footprinting reveals the modular organization of DNA regulatory elements

Cis-regulatory elements control gene expression and are dynamic in their structure, reflecting changes to the composition of diverse effector proteins over time1-3. Here we sought to connect the structural changes at cis-regulatory elements to alterations in cellular fate and function. To do this we developed PRINT, a computational method that uses deep learning to correct sequence bias in chromatin accessibility data and identifies multi-scale footprints of DNA-protein interactions. We find that multi-scale footprints enable more accurate inference of TF and nucleosome binding. Using PRINT with single-cell multi-omics, we discover wide-spread changes to the structure and function of candidate cis-regulatory elements (cCREs) across hematopoiesis, wherein nucleosomes slide, expose DNA for TF binding, and promote gene expression. Activity segmentation using the co-variance across cell states identifies "sub-cCREs" as modular cCRE subunits of regulatory DNA. We apply this single-cell and PRINT approach to characterize the age-associated alterations to cCREs within hematopoietic stem cells (HSCs). Remarkably, we find a spectrum of aging alterations among HSCs corresponding to a global gain of sub-cCRE activity while preserving cCRE accessibility. Collectively, we reveal the functional importance of cCRE structure across cell states, highlighting changes to gene regulation at single-cell and single-base-pair resolution.

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

The autism risk factor CHD8 is a chromatin activator in human neurons and functionally dependent on the ERK-MAPK pathway effector ELK1

The chromodomain helicase DNA-binding protein CHD8 is the most frequently mutated gene in autism spectrum disorder. Despite its prominent disease involvement, little is known about its molecular function in the human brain. CHD8 is a chromatin regulator which binds to the promoters of actively transcribed genes through genomic targeting mechanisms which have yet to be fully defined. By generating a conditional loss-of-function and an endogenously tagged allele in human pluripotent stem cells, we investigated the molecular function and the interaction of CHD8 with chromatin in human neurons. Chromatin accessibility analysis and transcriptional profiling revealed that CHD8 functions as a transcriptional activator at its target genes in human neurons. Furthermore, we found that CHD8 chromatin targeting is cell context-dependent. In human neurons, CHD8 preferentially binds at ETS motif-enriched promoters. This enrichment is particularly prominent on the promoters of genes whose expression significantly changes upon the loss of CHD8. Indeed, among the ETS transcription factors, we identified ELK1 as being most highly correlated with CHD8 expression in primary human fetal and adult cortical neurons and most highly expressed in our stem cell-derived neurons. Remarkably, ELK1 was necessary to recruit CHD8 specifically to ETS motif-containing sites. These findings imply that ELK1 and CHD8 functionally cooperate to regulate gene expression and chromatin states at MAPK/ERK target genes in human neurons. Our results suggest that the MAPK/ERK/ELK1 axis potentially contributes to the pathogenesis caused by CHD8 mutations in human neurodevelopmental disorders.

The chromatin remodeler CHD6 promotes colorectal cancer development by regulating TMEM65-mediated mitochondrial dynamics via EGF and Wnt signaling

Chromodomain helicase DNA binding protein (CHD) family plays critical roles in regulating gene transcription. The family is linked to cancer disease, but the family member's role in tumorigenesis remains largely unknown. Here, we report that CHD6 is highly expressed in colorectal cancer (CRC). CHD6 knockdown inhibited cancer cell proliferation, migration, invasion, and tumorigenesis. Consistently, Villin-specific Chd6 knockout in mice attenuates cancer formation in AOM/DSS model. We found that aberrant EGF signals promoted the stability of CHD6 by diminishing ubiquitin-mediated degradation. EGF signal inhibits GSK3β activity, which in turn prevents phosphodegron formation of CHD6, thereby hindering E3 ligase FBXW7-mediated CHD6 ubiquitination and degradation. CHD6's chromatin remodeler activity engages in binding Wnt signaling transcription factor TCF4 to facilitate the transcriptional expression of TMEM65, a mitochondrial inner membrane protein involved in ATP production and mitochondrial dynamics. In addition, Wnt signaling is also an upstream regulator of CHD6. CHD6 promoter contains TCF4 and β-catenin binding site, and CHD6 can be transcriptionally activated by Wnt ligand to facilitate TMEM65 transcription. Thus CHD6-TMEM65 axis can be regulated by both EGF and Wnt signaling pathways through two different mechanisms. We further illustrate that CHD6-TMEM65 axis is deregulated in cancer and that co-administration of Wnt inhibitor LGK974 and the anti-EGFR monoclonal antibody cetuximab largely restricted the growth of patient-derived xenografts of CRC. Targeting CHD6-TMEM65 axis may be effective for cancer intervention.

Cofactors: a new layer of specificity to enhancer regulation

Cofactors are essential effectors of the transcription control machinery. How this functionally diverse group of factors is used in the genome remains elusive. A recent study by Neumayr, Haberle et al. sheds light on this question, showing that enhancers depend on defined combinations of cofactors for their activation.

CHD6 promotes broad nucleosome eviction for transcriptional activation in prostate cancer cells

Despite being a member of the chromodomain helicase DNA-binding protein family, little is known about the exact role of CHD6 in chromatin remodeling or cancer disease. Here we show that CHD6 binds to chromatin to promote broad nucleosome eviction for transcriptional activation of many cancer pathways. By integrating multiple patient cohorts for bioinformatics analysis of over a thousand prostate cancer datasets, we found CHD6 expression elevated in prostate cancer and associated with poor prognosis. Further comprehensive experiments demonstrated that CHD6 regulates oncogenicity of prostate cancer cells and tumor development in a murine xenograft model. ChIP-Seq for CHD6, along with MNase-Seq and RNA-Seq, revealed that CHD6 binds on chromatin to evict nucleosomes from promoters and gene bodies for transcriptional activation of oncogenic pathways. These results demonstrated a key function of CHD6 in evicting nucleosomes from chromatin for transcriptional activation of prostate cancer pathways.

Decoding histone ubiquitylation

Histone ubiquitylation is a critical part of both active and repressed transcriptional states, and lies at the heart of DNA damage repair signaling. The histone residues targeted for ubiquitylation are often highly conserved through evolution, and extensive functional studies of the enzymes that catalyze the ubiquitylation and de-ubiquitylation of histones have revealed key roles linked to cell growth and division, development, and disease in model systems ranging from yeast to human cells. Nonetheless, the downstream consequences of these modifications have only recently begun to be appreciated on a molecular level. Here we review the structure and function of proteins that act as effectors or “readers” of histone ubiquitylation. We highlight lessons learned about how ubiquitin recognition lends specificity and function to intermolecular interactions in the context of transcription and DNA repair, as well as what this might mean for how we think about histone modifications more broadly.

Regulation of sensory perception and motor abilities by brain-specific action of chromatin remodeling factor CHD1

The ATP-dependent chromatin remodeling factor CHD1 (chromodomain-helicase-DNA binding protein 1) is involved in both the de novo assembly and the remodeling of chromatin. Recently, we discovered a crucial role of CHD1 in the incorporation of the histone variant H3.3 in the fly brain illustrated by widespread transcriptional upregulation and shortened lifespan in Chd1-mutant animals. Because many genes linked to sensory perception were dysregulated in Chd1-mutant heads, we studied the role of CHD1 in these processes. Here we show that Chd1-mutant flies have severe defects in their response behavior to olfactory and gustatory but not visual stimuli. Further analyses suggested that poor performance in gustatory response assays was caused by reduced motivation for foraging and feeding rather than defects in taste perception. Moreover, we show that shortened lifespan of Chd1-mutant flies is accompanied by indications of premature functional aging as suggested by defects in negative geotaxis and exploratory walking assays. The latter phenotype was rescued by neuronal re-expression of Chd1, while the olfactory defects were not. Interestingly, we found evidence for indirect regulation of the non-neuronal expression of odorant binding proteins (Obp) by neuronal expression of Chd1. Together, these results emphasize the crucial role of CHD1 activity controlling diverse neuronal processes thereby affecting healthy lifespan.

The Role of the Histone Variant H2A.Z in Metazoan Development

During the emergence and radiation of complex multicellular eukaryotes from unicellular ancestors, transcriptional systems evolved by becoming more complex to provide the basis for this morphological diversity. The way eukaryotic genomes are packaged into a highly complex structure, known as chromatin, underpins this evolution of transcriptional regulation. Chromatin structure is controlled by a variety of different epigenetic mechanisms, including the major mechanism for altering the biochemical makeup of the nucleosome by replacing core histones with their variant forms. The histone H2A variant H2A.Z is particularly important in early metazoan development because, without it, embryos cease to develop and die. However, H2A.Z is also required for many differentiation steps beyond the stage that H2A.Z-knockout embryos die. H2A.Z can facilitate the activation and repression of genes that are important for pluripotency and differentiation, and acts through a variety of different molecular mechanisms that depend upon its modification status, its interaction with histone and nonhistone partners, and where it is deposited within the genome. In this review, we discuss the current knowledge about the different mechanisms by which H2A.Z regulates chromatin function at various developmental stages and the chromatin remodeling complexes that determine when and where H2A.Z is deposited.

Ranking reprogramming factors for cell differentiation

Transcription factor over-expression is a proven method for reprogramming cells to a desired cell type for regenerative medicine and therapeutic discovery. However, a general method for the identification of reprogramming factors to create an arbitrary cell type is an open problem. Here we examine the success rate of methods and data for differentiation by testing the ability of nine computational methods (CellNet, GarNet, EBseq, AME, DREME, HOMER, KMAC, diffTF and DeepAccess) to discover and rank candidate factors for eight target cell types with known reprogramming solutions. We compare methods that use gene expression, biological networks and chromatin accessibility data, and comprehensively test parameter and preprocessing of input data to optimize performance. We find the best factor identification methods can identify an average of 50-60% of reprogramming factors within the top ten candidates, and methods that use chromatin accessibility perform the best. Among the chromatin accessibility methods, complex methods DeepAccess and diffTF have higher correlation with the ranked significance of transcription factor candidates within reprogramming protocols for differentiation. We provide evidence that AME and diffTF are optimal methods for transcription factor recovery that will allow for systematic prioritization of transcription factor candidates to aid in the design of new reprogramming protocols.

Chromatin accessibility shapes meiotic recombination in mouse primordial germ cells through assisting double-strand breaks and loop formation

Meiotic recombination is a driver of evolution, and aberrant recombination is a major contributor to aneuploidy in mammals. Mechanism of recombination remains elusive yet. Here, we present a computational analysis to explore recombination-related dynamics of chromatin accessibility in mouse primordial germ cells (PGCs). Our data reveals that: (1) recombination hotspots which get accessible at meiosis-specific DNase I-hypersensitive sites (DHSs) only when PGCs enter meiosis are located preferentially in intronic and distal intergenic regions; (2) stable DHSs maintained stably across PGC differentiation are enriched by CTCF motifs and CTCF binding and mediate chromatin loop formation; (3) compared with the specific DHSs aroused at meiotic stage, stable DHSs are largely encoded in DNA sequence and also enriched by epigenetic marks; (4) PRDM9 is likely to target nucleosome-occupied hotspot regions and remodels local chromatin structure to make them accessible for recombination machinery; and (5) cells undergoing meiotic recombination are deficient in TAD structure and chromatin loop arrays are organized regularly along the axis formed between homologous chromosomes. Taken together, by analyzing DHS-related DNA features, epigenetic marks and 3D genome structure, we revealed some specific roles of chromatin accessibility in recombination, which would expand our understanding of recombination mechanism.

Autism-associated protein POGZ controls ESCs and ESC neural induction by association with esBAF

Background

The POGZ gene has been found frequently mutated in neurodevelopmental disorders (NDDs), particularly autism spectrum disorder (ASD) and intellectual disability (ID). However, little is known about its roles in embryonic stem cells (ESCs), neural development and diseases.

Methods

We generated Pogz−/− ESCs and directed ESC differentiation toward a neural fate. We performed biochemistry, ChIP-seq, ATAC-seq, and bioinformatics analyses to understand the role of POGZ.

Results

We show that POGZ is required for the maintenance of ESC identity and the up-regulation of neural genes during ESC differentiation toward a neural fate. Genome-wide binding analysis shows that POGZ is primarily localized to gene promoter and enhancer regions. POGZ functions as both a transcriptional activator and repressor, and its loss leads to deregulation of differentiation genes, including neural genes. POGZ physically associates with the SWI-SNF (esBAF) chromatin remodeler complex, and together they modulate enhancer activities via epigenetic modifications such as chromatin remodeling and histone modification. During ESC neural induction, POGZ-mediated recruitment of esBAF/BRG1 and H3K27ac are important for proper expression of neural progenitor genes.

Limitations

The genotype and allele relevant to human neurodevelopmental disorders is heterozygous loss of function. This work is designed to study the effects of loss of POGZ function on ESCs and during ESC neural induction. Also, this work lacks of in vivo validation using animal models.

Conclusions

The data suggest that POGZ is both a transcription factor and a genome regulator, and its loss leads to defects in neural induction and neurogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13229-022-00502-9.

Inositol polyphosphate multikinase physically binds to the SWI/SNF complex and modulates BRG1 occupancy in mouse embryonic stem cells

Inositol polyphosphate multikinase (IPMK), a key enzyme in inositol polyphosphate (IP) metabolism, is a pleiotropic signaling factor involved in major biological events, including transcriptional control. In the yeast, IPMK and its IP products promote the activity of the chromatin remodeling complex SWI/SNF, which plays a critical role in gene expression by regulating chromatin accessibility. However, the direct link between IPMK and chromatin remodelers remains unclear, raising the question of how IPMK contributes to transcriptional regulation in mammals. By employing unbiased screening approaches and in vivo/in vitro immunoprecipitation, here we demonstrate that mammalian IPMK physically interacts with the SWI/SNF complex by directly binding to SMARCB1, BRG1, and SMARCC1. Furthermore, we identified the specific domains required for IPMK-SMARCB1 binding. Notably, using CUT&RUN and ATAC-seq assays, we discovered that IPMK co-localizes with BRG1 and regulates BRG1 localization as well as BRG1-mediated chromatin accessibility in a genome-wide manner in mouse embryonic stem cells. Together, these findings show that IPMK regulates the promoter targeting of the SWI/SNF complex, thereby contributing to SWI/SNF-meditated chromatin accessibility, transcription, and differentiation in mouse embryonic stem cells.

CHD chromatin remodeling protein diversification yields novel clades and domains absent in classic model organisms

Chromatin remodelers play a fundamental role in the assembly of chromatin, regulation of transcription, and DNA repair. Biochemical and functional characterizations of the CHD family of chromatin remodelers from a variety of model organisms have shown that these remodelers participate in a wide range of activities. However, because the evolutionary history of CHD homologs is unclear, it is difficult to predict which of these activities are broadly conserved and which have evolved more recently in individual eukaryotic lineages. Here, we performed a comprehensive phylogenetic analysis of 8,042 CHD homologs from 1,894 species to create a model for the evolution of this family across eukaryotes with a particular focus on the timing of duplications that gave rise to the diverse copies observed in plants, animals, and fungi. Our analysis confirms that the three major subfamilies of CHD remodelers originated in the eukaryotic last common ancestor, and subsequent losses occurred independently in different lineages. Improved taxon sampling identified several subfamilies of CHD remodelers in plants that were absent or highly divergent in the model plant Arabidopsis thaliana. Whereas the timing of CHD subfamily expansions in vertebrates corresponds to whole genome duplication events, the mechanisms underlying CHD diversification in land plants appear more complicated. Analysis of protein domains reveals that CHD remodeler diversification has been accompanied by distinct transitions in domain architecture, contributing to the functional differences observed between these remodelers. This study demonstrates the importance of proper taxon sampling when studying ancient evolutionary events to prevent misinterpretation of subsequent lineage-specific changes and provides an evolutionary framework for functional and comparative analysis of this critical chromatin remodeler family across eukaryotes.

Distinct transcription kinetics of pluripotent cell states

Mouse embryonic stem cells (mESCs) can adopt naïve, ground, and paused pluripotent states that give rise to unique transcriptomes. Here, we use transient transcriptome sequencing (TT-seq) to define both coding and non-coding transcription units (TUs) in these three pluripotent states and combine TT-seq with RNA polymerase II occupancy profiling to unravel the kinetics of RNA metabolism genome-wide. Compared to the naïve state (serum), RNA synthesis and turnover rates are globally reduced in the ground state (2i) and the paused state (mTORi). The global reduction in RNA synthesis goes along with a genome-wide decrease of polymerase elongation velocity, which is related to epigenomic features and alterations in the Pol II termination window. Our data suggest that transcription activity is the main determinant of steady state mRNA levels in the naïve state and that genome-wide changes in transcription kinetics invoke ground and paused pluripotent states.

CHD4 Conceals Aberrant CTCF-Binding Sites at TAD Interiors by Regulating Chromatin Accessibility in Mouse Embryonic Stem Cells

CCCTC-binding factor (CTCF) critically contributes to 3D chromatin organization by determining topologically associated domain (TAD) borders. Although CTCF primarily binds at TAD borders, there also exist putative CTCF-binding sites within TADs, which are spread throughout the genome by retrotransposition. However, the detailed mechanism responsible for masking the putative CTCF-binding sites remains largely elusive. Here, we show that the ATP-dependent chromatin remodeler, chromodomain helicase DNA-binding 4 (CHD4), regulates chromatin accessibility to conceal aberrant CTCF-binding sites embedded in H3K9me3-enriched heterochromatic B2 short interspersed nuclear elements (SINEs) in mouse embryonic stem cells (mESCs). Upon CHD4 depletion, these aberrant CTCF-binding sites become accessible and aberrant CTCF recruitment occurs within TADs, resulting in disorganization of local TADs. RNA-binding intrinsically disordered domains (IDRs) of CHD4 are required to prevent this aberrant CTCF binding, and CHD4 is critical for the repression of B2 SINE transcripts. These results collectively reveal that a CHD4-mediated mechanism ensures appropriate CTCF binding and associated TAD organization in mESCs.

In vivo targeted DamID identifies CHD8 genomic targets in fetal mouse brain

Genetic studies of autism have revealed causal roles for chromatin remodeling gene mutations. Chromodomain helicase DNA binding protein 8 (CHD8) encodes a chromatin remodeler with significant de novo mutation rates in sporadic autism. However, relationships between CHD8 genomic function and autism-relevant biology remain poorly elucidated. Published studies utilizing ChIP-seq to map CHD8 protein-DNA interactions have high variability, consistent with technical challenges and limitations associated with this method. Thus, complementary approaches are needed to establish CHD8 genomic targets and regulatory functions in developing brain. We used in utero CHD8 Targeted DamID followed by sequencing (TaDa-seq) to characterize CHD8 binding in embryonic mouse cortex. CHD8 TaDa-seq reproduced interaction patterns observed from ChIP-seq and further highlighted CHD8 distal interactions associated with neuronal loci. This study establishes TaDa-seq as a useful alternative for mapping protein-DNA interactions in vivo and provides insights into the regulatory targets of CHD8 and autism-relevant pathophysiology associated with CHD8 mutations.

Comprehensive understanding of Tn5 insertion preference improves transcription regulatory element identification

Tn5 transposase, which can efficiently tagment the genome, has been widely adopted as a molecular tool in next-generation sequencing, from short-read sequencing to more complex methods such as assay for transposase-accessible chromatin using sequencing (ATAC-seq). Here, we systematically map Tn5 insertion characteristics across several model organisms, finding critical parameters that affect its insertion. On naked genomic DNA, we found that Tn5 insertion is not uniformly distributed or random. To uncover drivers of these biases, we used a machine learning framework, which revealed that DNA shape cooperatively works with DNA motif to affect Tn5 insertion preference. These intrinsic insertion preferences can be modeled using nucleotide dependence information from DNA sequences, and we developed a computational pipeline to correct for these biases in ATAC-seq data. Using our pipeline, we show that bias correction improves the overall performance of ATAC-seq peak detection, recovering many potential false-negative peaks. Furthermore, we found that these peaks are bound by transcription factors, underscoring the biological relevance of capturing this additional information. These findings highlight the benefits of an improved understanding and precise correction of Tn5 insertion preference.

Acute depletion of the ARID1A subunit of SWI/SNF complexes reveals distinct pathways for activation and repression of transcription

The ARID1A subunit of SWI/SNF chromatin remodeling complexes is a potent tumor suppressor. Here, a degron is applied to detect rapid loss of chromatin accessibility at thousands of loci where ARID1A acts to generate accessible minidomains of nucleosomes. Loss of ARID1A also results in the redistribution of the coactivator EP300. Co-incident EP300 dissociation and lost chromatin accessibility at enhancer elements are highly enriched adjacent to rapidly downregulated genes. In contrast, sites of gained EP300 occupancy are linked to genes that are transcriptionally upregulated. These chromatin changes are associated with a small number of genes that are differentially expressed in the first hours following loss of ARID1A. Indirect or adaptive changes dominate the transcriptome following growth for days after loss of ARID1A and result in strong engagement with cancer pathways. The identification of this hierarchy suggests sites for intervention in ARID1A-driven diseases.

•

Degradation of ARID1A disrupts nucleosomes flanking pluripotency transcription factors

•

EP300 is rapidly redistributed with increased occupancy adjacent to upregulated genes

•

These changes are associated with misregulation of a few hundred genes within 2 h

•

During subsequent days, widespread indirect changes mimic a premalignant state

Sentinels of chromatin: chromodomain helicase DNA-binding proteins in development and disease

Chromatin is highly dynamic, undergoing continuous global changes in its structure and type of histone and DNA modifications governed by processes such as transcription, repair, replication, and recombination. Members of the chromodomain helicase DNA-binding (CHD) family of enzymes are ATP-dependent chromatin remodelers that are intimately involved in the regulation of chromatin dynamics, altering nucleosomal structure and DNA accessibility. Genetic studies in yeast, fruit flies, zebrafish, and mice underscore essential roles of CHD enzymes in regulating cellular fate and identity, as well as proper embryonic development. With the advent of next-generation sequencing, evidence is emerging that these enzymes are subjected to frequent DNA copy number alterations or mutations and show aberrant expression in malignancies and other human diseases. As such, they might prove to be valuable biomarkers or targets for therapeutic intervention.

Advances in Chromodomain Helicase DNA-Binding (CHD) Proteins Regulating Stem Cell Differentiation and Human Diseases

Background: Regulation of gene expression is critical for stem cell differentiation, tissue development, and human health maintenance. Recently, epigenetic modifications of histone and chromatin remodeling have been verified as key controllers of gene expression and human diseases.

Objective: In this study, we review the role of chromodomain helicase DNA-binding (CHD) proteins in stem cell differentiation, cell fate decision, and several known human developmental disorders and cancers.

Conclusion: CHD proteins play a crucial role in stem cell differentiation and human diseases.

CHD1 controls H3.3 incorporation in adult brain chromatin to maintain metabolic homeostasis and normal lifespan

The ATP-dependent chromatin remodeling factor CHD1 is essential for the assembly of variant histone H3.3 into paternal chromatin during sperm chromatin remodeling in fertilized eggs. It remains unclear, however, if CHD1 has a similar role in normal diploid cells. Using a specifically tailored quantitative mass spectrometry approach, we show that Chd1 disruption results in reduced H3.3 levels in heads of Chd1 mutant flies. Chd1 deletion perturbs brain chromatin structure in a similar way as H3.3 deletion and leads to global de-repression of transcription. The physiological consequences are reduced food intake, metabolic alterations, and shortened lifespan. Notably, brain-specific CHD1 expression rescues these phenotypes. We further demonstrate a strong genetic interaction between Chd1 and H3.3 chaperone Hira. Thus, our findings establish CHD1 as a factor required for the assembly of H3.3-containing chromatin in adult cells and suggest a crucial role for CHD1 in the brain as a regulator of organismal health and longevity.

BAZ2A-mediated repression via H3K14ac-marked enhancers promotes prostate cancer stem cells

Prostate cancer (PCa) is one of the most prevalent cancers in men. Cancer stem cells are thought to be associated with PCa relapse. Here, we show that BAZ2A is required for PCa cells with a cancer stem‐like state. BAZ2A genomic occupancy in PCa cells coincides with H3K14ac‐enriched chromatin regions. This association is mediated by BAZ2A‐bromodomain (BAZ2A‐BRD) that specifically binds H3K14ac. BAZ2A associates with inactive enhancers marked by H3K14ac and repressing transcription of genes frequently silenced in aggressive and poorly differentiated PCa. BAZ2A‐mediated repression is also linked to EP300 that acetylates H3K14ac. BAZ2A‐BRD mutations or treatment with inhibitors abrogating BAZ2A‐BRD/H3K14ac interaction impair PCa stem cells. Furthermore, pharmacological inactivation of BAZ2A‐BRD impairs Pten‐loss oncogenic transformation of prostate organoids. Our findings indicate a role of BAZ2A‐BRD in PCa stem cell features and suggest potential epigenetic‐reader therapeutic strategies to target BAZ2A in aggressive PCa.

INO80 promotes H2A.Z occupancy to regulate cell fate transition in pluripotent stem cells

The INO80 chromatin remodeler is involved in many chromatin-dependent cellular functions. However, its role in pluripotency and cell fate transition is not fully defined. We examined the impact of Ino80 deletion in the naïve and primed pluripotent stem cells. We found that Ino80 deletion had minimal effect on self-renewal and gene expression in the naïve state, but led to cellular differentiation and de-repression of developmental genes in the transition toward and maintenance of the primed state. In the naïve state, INO80 pre-marked gene promoters that would adopt bivalent histone modifications by H3K4me3 and H3K27me3 upon transition into the primed state. In the primed state, in contrast to its known role in H2A.Z exchange, INO80 promoted H2A.Z occupancy at these bivalent promoters and facilitated H3K27me3 installation and maintenance as well as downstream gene repression. Together, our results identified an unexpected function of INO80 in H2A.Z deposition and gene regulation. We showed that INO80-dependent H2A.Z occupancy is a critical licensing step for the bivalent domains, and thereby uncovered an epigenetic mechanism by which chromatin remodeling, histone variant deposition and histone modification coordinately control cell fate.

AIAP: A Quality Control and Integrative Analysis Package to Improve ATAC-seq Data Analysis

Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) is a technique widely used to investigate genome-wide chromatin accessibility. The recently published Omni-ATAC-seq protocol substantially improves the signal/noise ratio and reduces the input cell number. High-quality data are critical to ensure accurate analysis. Several tools have been developed for assessing sequencing quality and insertion size distribution for ATAC-seq data; however, key quality control (QC) metrics have not yet been established to accurately determine the quality of ATAC-seq data. Here, we optimized the analysis strategy for ATAC-seq and defined a series of QC metrics for ATAC-seq data, including reads under peak ratio (RUPr), background (BG), promoter enrichment (ProEn), subsampling enrichment (SubEn), and other measurements. We incorporated these QC tests into our recently developed ATAC-seq Integrative Analysis Package (AIAP) to provide a complete ATAC-seq analysis system, including quality assurance, improved peak calling, and downstream differential analysis. We demonstrated a significant improvement of sensitivity (20%–60%) in both peak calling and differential analysis by processing paired-end ATAC-seq datasets using AIAP. AIAP is compiled into Docker/Singularity, and it can be executed by one command line to generate a comprehensive QC report. We used ENCODE ATAC-seq data to benchmark and generate QC recommendations, and developed qATACViewer for the user-friendly interaction with the QC report. The software, source code, and documentation of AIAP are freely available at https://github.com/Zhang-lab/ATAC-seq_QC_analysis.

mSWI/SNF promotes Polycomb repression both directly and through genome-wide redistribution

The mammalian SWI/SNF complex, or BAF complex, has a conserved and direct role in antagonizing Polycomb-mediated repression. Yet, BAF also promotes repression by Polycomb in stem cells and cancer. How BAF both antagonizes and promotes Polycomb-mediated repression remains unknown. Here, we utilize targeted protein degradation to dissect the BAF-Polycomb axis in mouse embryonic stem cells on short timescales. We report that rapid BAF depletion redistributes Polycomb repressive complexes PRC1 and PRC2 from highly occupied domains, like Hox clusters, to weakly occupied sites normally opposed by BAF. Polycomb redistribution from highly repressed domains results in their decompaction, gain of active epigenomic features and transcriptional derepression. Surprisingly, through dose-dependent degradation of PRC1 and PRC2, we identify a conventional role for BAF in Polycomb-mediated repression, in addition to global Polycomb redistribution. These findings provide new mechanistic insight into the highly dynamic state of the Polycomb-Trithorax axis.

Overarching control of autophagy and DNA damage response by CHD6 revealed by modeling a rare human pathology

Members of the chromodomain-helicase-DNA binding (CHD) protein family are chromatin remodelers implicated in human pathologies, with CHD6 being one of its least studied members. We discovered a de novo CHD6 missense mutation in a patient clinically presenting the rare Hallermann-Streiff syndrome (HSS). We used genome editing to generate isogenic iPSC lines and model HSS in relevant cell types. By combining genomics with functional in vivo and in vitro assays, we show that CHD6 binds a cohort of autophagy and stress response genes across cell types. The HSS mutation affects CHD6 protein folding and impairs its ability to recruit co-remodelers in response to DNA damage or autophagy stimulation. This leads to accumulation of DNA damage burden and senescence-like phenotypes. We therefore uncovered a molecular mechanism explaining HSS onset via chromatin control of autophagic flux and genotoxic stress surveillance.