CpG islands and GC content dictate nucleosome depletion in a transcription-independent manner at mammalian promoters

The GC-content at the 5' ends of human protein-coding genes is undergoing mutational decay

BACKGROUND

In vertebrates, most protein-coding genes have a peak of GC-content near their 5' transcriptional start site (TSS). This feature promotes both the efficient nuclear export and translation of mRNAs. Despite the importance of GC-content for RNA metabolism, its general features, origin, and maintenance remain mysterious. We investigate the evolutionary forces shaping GC-content at the transcriptional start site (TSS) of genes through both comparative genomic analysis of nucleotide substitution rates between different species and by examining human de novo mutations.

RESULTS

Our data suggests that GC-peaks at TSSs were present in the last common ancestor of amniotes, and likely that of vertebrates. We observe that in apes and rodents, where recombination is directed away from TSSs by PRDM9, GC-content at the 5' end of protein-coding gene is currently undergoing mutational decay. In canids, which lack PRDM9 and perform recombination at TSSs, GC-content at the 5' end of protein-coding is increasing. We show that these patterns extend into the 5' end of the open reading frame, thus impacting synonymous codon position choices.

CONCLUSIONS

Our results indicate that the dynamics of this GC-peak in amniotes is largely shaped by historic patterns of recombination. Since decay of GC-content towards the mutation rate equilibrium is the default state for non-functional DNA, the observed decrease in GC-content at TSSs in apes and rodents indicates that the GC-peak is not being maintained by selection on most protein-coding genes in those species.

HMGA1 orchestrates chromatin compartmentalization and sequesters genes into 3D networks coordinating senescence heterogeneity

HMGA1 is an abundant non-histone chromatin protein that has been implicated in embryonic development, cancer, and cellular senescence, but its specific role remains elusive. Here, we combine functional genomics approaches with graph theory to investigate how HMGA1 genomic deposition controls high-order chromatin networks in an oncogene-induced senescence model. While the direct role of HMGA1 in gene activation has been described previously, we find little evidence to support this. Instead, we show that the heterogeneous linear distribution of HMGA1 drives a specific 3D chromatin organization. HMGA1-dense loci form highly interactive networks, similar to, but independent of, constitutive heterochromatic loci. This, coupled with the exclusion of HMGA1-poor chromatin regions, leads to coordinated gene regulation through the repositioning of genes. In the absence of HMGA1, the whole process is largely reversed, but many regulatory interactions also emerge, amplifying the inflammatory senescence-associated secretory phenotype. Such HMGA1-mediated fine-tuning of gene expression contributes to the heterogeneous nature of senescence at the single-cell level. A similar 'buffer' effect of HMGA1 on inflammatory signalling is also detected in lung cancer cells. Our study reveals a mechanism through which HMGA1 modulates chromatin compartmentalization and gene regulation in senescence and beyond.

Genome-wide Methylation Dynamics and Context-dependent Gene Expression Variability in Differentiating Preadipocytes

Obesity, characterized by the accumulation of excess fat, is a complex condition resulting from the combination of genetic and epigenetic factors. Recent studies have found correspondence between DNA methylation and cell differentiation, suggesting a role of the former in cell fate determination. There is a lack of comprehensive understanding concerning the underpinnings of preadipocyte differentiation, specifically when cells are undergoing terminal differentiation (TD). To gain insight into dynamic genome-wide methylation, 3T3 L1 preadipocyte cells were differentiated by a hormone cocktail. The genomic DNA was isolated from undifferentiated cells and 4 hours, 2 days postdifferentiated cells, and 15 days TD cells. We employed whole-genome bisulfite sequencing (WGBS) to ascertain global genomic DNA methylation alterations at single base resolution as preadipocyte cells differentiate. The genome-wide distribution of DNA methylation showed similar overall patterns in pre-, post-, and terminally differentiated adipocytes, according to WGBS analysis. DNA methylation decreases at 4 hours after differentiation initiation, followed by methylation gain as cells approach TD. Studies revealed novel differentially methylated regions (DMRs) associated with adipogenesis. DMR analysis suggested that though DNA methylation is global, noticeable changes are observed at specific sites known as "hotspots." Hotspots are genomic regions rich in transcription factor (TF) binding sites and exhibit methylation-dependent TF binding. Subsequent analysis indicated hotspots as part of DMRs. The gene expression profile of key adipogenic genes in differentiating adipocytes is context-dependent, as we found a direct and inverse relationship between promoter DNA methylation and gene expression.

Pioneer Transcription Factors: The First Domino in Zygotic Genome Activation

Zygotic genome activation (ZGA) is a pivotal event in mammalian embryogenesis, marking the transition from maternal to zygotic control of development. During the ZGA process that is characterized by the intricate cascade of gene expression, who tipped the first domino in a meticulously arranged sequence is a subject of paramount interest. Recently, Dux, Obox and Nr5a2 were identified as pioneer transcription factors that reside at the top of transcriptional hierarchy. Through co-option of retrotransposon elements as hubs for transcriptional activation, these pioneer transcription factors rewire the gene regulatory network, thus initiating ZGA. In this review, we provide a snapshot of the mechanisms underlying the functions of these pioneer transcription factors. We propose that ZGA is the starting point where the embryo's own genome begins to influence development trajectory, therefore in-depth dissecting the functions of pioneer transcription factors during ZGA will form a cornerstone of our understanding for early embryonic development, which will pave the way for advancing our grasp of mammalian developmental biology and optimizing in vitro production (IVP) techniques.

A transcription factor (TF) inference method that broadly measures TF activity and identifies mechanistically distinct TF networks

TF profiler is a method of inferring transcription factor regulatory activity, i.e. when a TF is present and actively regulating transcription, directly directly from nascent sequencing assays such as PRO-seq and GRO-seq. Transcription factors orchestrate transcription and play a critical role in cellular maintenance, identity and response to external stimuli. While ChIP assays have measured DNA localization, they fall short of identifying when and where transcription factors are actively regulating transcription. Our method, on the other hand, uses RNA polymerase activity to infer TF activity across hundreds of data sets and transcription factors. Based on these classifications we identify three distinct classes of transcription factors: ubiquitous factors that play roles in cellular homeostasis, driving basal gene programs across tissues and cell types, tissue specific factors that act almost exclusively at enhancers and are themselves regulated at transcription, and stimulus responsive TFs which are regulated post-transcriptionally but act predominantly at enhancers. TF profiler is broadly applicable, providing regulatory insights on any PRO-seq sample for any transcription factor with a known binding motif.

Transcriptional regulation of FACT involves Coordination of chromatin accessibility and CTCF binding

Histone chaperone FACT (facilitates chromatin transcription) is well known to promote chromatin recovery during transcription. However, the mechanism how FACT regulates genome-wide chromatin accessibility and transcription factor binding has not been fully elucidated. Through loss-of-function studies, we show here that FACT component Ssrp1 is required for DNA replication and DNA damage repair and is also essential for progression of cell phase transition and cell proliferation in mouse embryonic fibroblast cells. On the molecular level, absence of the Ssrp1 leads to increased chromatin accessibility, enhanced CTCF binding, and a remarkable change in dynamic range of gene expression. Our study thus unequivocally uncovers a unique mechanism by which FACT complex regulates transcription by coordinating genome-wide chromatin accessibility and CTCF binding.

Single-cell methylation analysis of brain tissue prioritizes mutations that alter transcription

Relating genetic variants to behavior remains a fundamental challenge. To assess the utility of DNA methylation marks in discovering causative variants, we examined their relationship to genetic variation by generating single-nucleus methylomes from the hippocampus of eight inbred mouse strains. At CpG sequence densities under 40 CpG/Kb, cells compensate for loss of methylated sites by methylating additional sites to maintain methylation levels. At higher CpG sequence densities, the exact location of a methylated site becomes more important, suggesting that variants affecting methylation will have a greater effect when occurring in higher CpG densities than in lower. We found this to be true for a variant's effect on transcript abundance, indicating that candidate variants can be prioritized based on CpG sequence density. Our findings imply that DNA methylation influences the likelihood that mutations occur at specific sites in the genome, supporting the view that the distribution of mutations is not random.

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.

Tuning Methylation-Dependent Silencing Dynamics by Synthetic Modulation of CpG Density

Methylation of cytosines in CG dinucleotides (CpGs) within promoters has been shown to lead to gene silencing in mammals in natural contexts. Recently, engineered recruitment of methyltransferases (DNMTs) at specific loci was shown to be sufficient to silence synthetic and endogenous gene expression through this mechanism. A critical parameter for DNA methylation-based silencing is the distribution of CpGs within the target promoter. However, how the number or density of CpGs in the target promoter affects the dynamics of silencing by DNMT recruitment has remained unclear. Here, we constructed a library of promoters with systematically varying CpG content, and analyzed the rate of silencing in response to recruitment of DNMT. We observed a tight correlation between silencing rate and CpG content. Further, methylation-specific analysis revealed a constant accumulation rate of methylation at the promoter after DNMT recruitment. We identified a single CpG site between TATA box and transcription start site (TSS) that accounted for a substantial part of the difference in silencing rates between promoters with differing CpG content, indicating that certain residues play disproportionate roles in controlling silencing. Together, these results provide a library of promoters for synthetic epigenetic and gene regulation applications, as well as insights into the regulatory link between CpG content and silencing rate.

RNA polymerase II CTD is dispensable for transcription and required for termination in human cells

The largest subunit of RNA polymerase (Pol) II harbors an evolutionarily conserved C-terminal domain (CTD), composed of heptapeptide repeats, central to the transcriptional process. Here, we analyze the transcriptional phenotypes of a CTD-Δ5 mutant that carries a large CTD truncation in human cells. Our data show that this mutant can transcribe genes in living cells but displays a pervasive phenotype with impaired termination, similar to but more severe than previously characterized mutations of CTD tyrosine residues. The CTD-Δ5 mutant does not interact with the Mediator and Integrator complexes involved in the activation of transcription and processing of RNAs. Examination of long-distance interactions and CTCF-binding patterns in CTD-Δ5 mutant cells reveals no changes in TAD domains or borders. Our data demonstrate that the CTD is largely dispensable for the act of transcription in living cells. We propose a model in which CTD-depleted Pol II has a lower entry rate onto DNA but becomes pervasive once engaged in transcription, resulting in a defect in termination.

High-sensitive nascent transcript sequencing reveals BRD4-specific control of widespread enhancer and target gene transcription

Gene transcription by RNA polymerase II (Pol II) is under control of promoters and distal regulatory elements known as enhancers. Enhancers are themselves transcribed by Pol II correlating with their activity. How enhancer transcription is regulated and coordinated with transcription at target genes has remained unclear. Here, we developed a high-sensitive native elongating transcript sequencing approach, called HiS-NET-seq, to provide an extended high-resolution view on transcription, especially at lowly transcribed regions such as enhancers. HiS-NET-seq uncovers new transcribed enhancers in human cells. A multi-omics analysis shows that genome-wide enhancer transcription depends on the BET family protein BRD4. Specifically, BRD4 co-localizes to enhancer and promoter-proximal gene regions, and is required for elongation activation at enhancers and their genes. BRD4 keeps a set of enhancers and genes in proximity through long-range contacts. From these studies BRD4 emerges as a general regulator of enhancer transcription that may link transcription at enhancers and genes.

OBOX regulates mouse zygotic genome activation and early development

Zygotic genome activation (ZGA) activates the quiescent genome to enable the maternal-to-zygotic transition1,2. However, the identity of transcription factors that underlie mammalian ZGA in vivo remains elusive. Here we show that OBOX, a PRD-like homeobox domain transcription factor family (OBOX1-OBOX8)3-5, are key regulators of mouse ZGA. Mice deficient for maternally transcribed Obox1/2/5/7 and zygotically expressed Obox3/4 had a two-cell to four-cell arrest, accompanied by impaired ZGA. The Obox knockout defects could be rescued by restoring either maternal and zygotic OBOX, which suggests that maternal and zygotic OBOX redundantly support embryonic development. Chromatin-binding analysis showed that Obox knockout preferentially affected OBOX-binding targets. Mechanistically, OBOX facilitated the 'preconfiguration' of RNA polymerase II, as the polymerase relocated from the initial one-cell binding targets to ZGA gene promoters and distal enhancers. Impaired polymerase II preconfiguration in Obox mutants was accompanied by defective ZGA and chromatin accessibility transition, as well as aberrant activation of one-cell polymerase II targets. Finally, ectopic expression of OBOX activated ZGA genes and MERVL repeats in mouse embryonic stem cells. These data thus demonstrate that OBOX regulates mouse ZGA and early embryogenesis.

G4access identifies G-quadruplexes and their associations with open chromatin and imprinting control regions

Metazoan promoters are enriched in secondary DNA structure-forming motifs, such as G-quadruplexes (G4s). Here we describe 'G4access', an approach to isolate and sequence G4s associated with open chromatin via nuclease digestion. G4access is antibody- and crosslinking-independent and enriches for computationally predicted G4s (pG4s), most of which are confirmed in vitro. Using G4access in human and mouse cells, we identify cell-type-specific G4 enrichment correlated with nucleosome exclusion and promoter transcription. G4access allows measurement of variations in G4 repertoire usage following G4 ligand treatment, HDAC and G4 helicases inhibitors. Applying G4access to cells from reciprocal hybrid mouse crosses suggests a role for G4s in the control of active imprinting regions. Consistently, we also observed that G4access peaks are unmethylated, while methylation at pG4s correlates with nucleosome repositioning on DNA. Overall, our study provides a new tool for studying G4s in cellular dynamics and highlights their association with open chromatin, transcription and their antagonism to DNA methylation.

Transcriptional and epigenomic profiling identifies YAP signaling as a key regulator of intestinal epithelium maturation

During intestinal organogenesis, equipotent epithelial progenitors mature into phenotypically distinct stem cells that are responsible for lifelong maintenance of the tissue. While the morphological changes associated with the transition are well characterized, the molecular mechanisms underpinning the maturation process are not fully understood. Here, we leverage intestinal organoid cultures to profile transcriptional, chromatin accessibility, DNA methylation, and three-dimensional (3D) chromatin conformation landscapes in fetal and adult epithelial cells. We observed prominent differences in gene expression and enhancer activity, which are accompanied by local changes in 3D organization, DNA accessibility, and methylation between the two cellular states. Using integrative analyses, we identified sustained Yes-Associated Protein (YAP) transcriptional activity as a major gatekeeper of the immature fetal state. We found the YAP-associated transcriptional network to be regulated at various levels of chromatin organization and likely to be coordinated by changes in extracellular matrix composition. Together, our work highlights the value of unbiased profiling of regulatory landscapes for the identification of key mechanisms underlying tissue maturation.

The Interplay between Dysregulated Metabolism and Epigenetics in Cancer

Cellular metabolism (or energetics) and epigenetics are tightly coupled cellular processes. It is arguable that of all the described cancer hallmarks, dysregulated cellular energetics and epigenetics are the most tightly coregulated. Cellular metabolic states regulate and drive epigenetic changes while also being capable of influencing, if not driving, epigenetic reprogramming. Conversely, epigenetic changes can drive altered and compensatory metabolic states. Cancer cells meticulously modify and control each of these two linked cellular processes in order to maintain their tumorigenic potential and capacity. This review aims to explore the interplay between these two processes and discuss how each affects the other, driving and enhancing tumorigenic states in certain contexts.

Tuning methylation-dependent silencing dynamics by synthetic modulation of CpG density

Methylation of cytosines in CG dinucleotides (CpGs) within promoters has been shown to lead to gene silencing in mammals in natural contexts. Recently, engineered recruitment of methyltransferases (DNMTs) at specific loci was shown to be sufficient to silence synthetic and endogenous gene expression through this mechanism. A critical parameter for DNA methylation-based silencing is the distribution of CpGs within the target promoter. However, how the number or density of CpGs in the target promoter affects the dynamics of silencing by DNMT recruitment has remained unclear. Here we constructed a library of promoters with systematically varying CpG content, and analyzed the rate of silencing in response to recruitment of DNMT. We observed a tight correlation between silencing rate and CpG content. Further, methylation-specific analysis revealed a constant accumulation rate of methylation at the promoter after DNMT recruitment. We identified a single CpG site between TATA box and transcription start site (TSS) that accounted for a substantial part of the difference in silencing rates between promoters with differing CpG content, indicating that certain residues play disproportionate roles in controlling silencing. Together, these results provide a library of promoters for synthetic epigenetic and gene regulation applications, as well as insights into the regulatory link between CpG content and silencing rate.

Alternative promoters in CpG depleted regions are prevalently associated with epigenetic misregulation of liver cancer transcriptomes

Transcriptional regulation is commonly governed by alternative promoters. However, the regulatory architecture in alternative and reference promoters, and how they differ, remains elusive. In 100 CAGE-seq libraries from hepatocellular carcinoma patients, here we annotate 4083 alternative promoters in 2926 multi-promoter genes, which are largely undetected in normal livers. These genes are enriched in oncogenic processes and predominantly show association with overall survival. Alternative promoters are narrow nucleosome depleted regions, CpG island depleted, and enriched for tissue-specific transcription factors. Globally tumors lose DNA methylation. We show hierarchical retention of intragenic DNA methylation with CG-poor regions rapidly losing methylation, while CG-rich regions retain it, a process mediated by differential SETD2, H3K36me3, DNMT3B, and TET1 binding. This mechanism is validated in SETD2 knockdown cells and SETD2-mutated patients. Selective DNA methylation loss in CG-poor regions makes the chromatin accessible for alternative transcription. We show alternative promoters can control tumor transcriptomes and their regulatory architecture.

Systematic evaluation of chromatin immunoprecipitation sequencing to study histone occupancy in dormancy transitions of grapevine buds

The regulation of DNA accessibility by histone modification has emerged as a paradigm of developmental and environmental programming. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a versatile tool to investigate in vivo protein-DNA interaction and has enabled advances in mechanistic understanding of physiologies. The technique has been successfully demonstrated in several plant species and tissues; however, it has remained challenging in woody tissues, in particular complex structures such as perennating buds. Here we developed a ChIP method specifically for mature dormant buds of grapevine (Vitis vinifera cv. Cabernet Sauvignon). Each step of the protocol was systematically optimized, including crosslinking, chromatin extraction, sonication and antibody validation. Analysis of histone H3-enriched DNA was performed to evaluate the success of the protocol and identify occupancy of histone H3 along grapevine bud chromatin. To our best knowledge, this is the first ChIP experiment protocol optimized for the grapevine bud system.

Normalization benchmark of ATAC-seq datasets shows the importance of accounting for GC-content effects

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) allows the study of epigenetic regulation of gene expression by assessing chromatin configuration for an entire genome. Despite its popularity, there have been limited studies investigating the analytical challenges related to ATAC-seq data, with most studies leveraging tools developed for bulk transcriptome sequencing. Here, we show that GC-content effects are omnipresent in ATAC-seq datasets. Since the GC-content effects are sample specific, they can bias downstream analyses such as clustering and differential accessibility analysis. We introduce a normalization method based on smooth-quantile normalization within GC-content bins and evaluate it together with 11 different normalization procedures on 8 public ATAC-seq datasets. Accounting for GC-content effects in the normalization is crucial for common downstream ATAC-seq data analyses, improving accuracy and interpretability. Through case studies, we show that exploratory data analysis is essential to guide the choice of an appropriate normalization method for a given dataset.

Uncovering the mechanisms of transcription elongation by eukaryotic RNA polymerases I, II, and III

Eukaryotes express three nuclear RNA polymerases (Pols I, II, and III) that are essential for cell survival. Despite extensive investigation of the three Pols, significant knowledge gaps regarding their biochemical properties remain because each Pol has been evaluated independently under disparate experimental conditions and methodologies. To advance our understanding of the Pols, we employed identical in vitro transcription assays for direct comparison of their elongation rates, elongation complex (EC) stabilities, and fidelities. Pol I is the fastest, most likely to misincorporate, forms the least stable EC, and is most sensitive to alterations in reaction buffers. Pol II is the slowest of the Pols, forms the most stable EC, and negligibly misincorporated an incorrect nucleotide. The enzymatic properties of Pol III were intermediate between Pols I and II in all assays examined. These results reveal unique enzymatic characteristics of the Pols that provide new insights into their evolutionary divergence.

Contribution of CRISPRable DNA to human complex traits

CRISPR-Cas is a powerful genome editing tool for various species and human cell lines, widely used in many research areas including studying the mechanisms, targets, and gene therapies of human diseases. Recent developments have even allowed high-throughput genetic screening using the CRISPR system. However, due to the practical and ethical limitations in human gene editing research, little is known about whether CRISPR-editable DNA segments could influence human complex traits or diseases. Here, we investigated the human genomic regions condensed with different CRISPR Cas enzymes’ protospacer-adjacent motifs (PAMs). We found that Cas enzymes with GC-rich PAMs could interfere more with the genomic regions that harbor enriched heritability for human complex traits and diseases. The results linked GC content across the genome to the functional genomic elements in the heritability enrichment of human complex traits. We provide a genetic overview of the effects of high-throughput genome editing on human complex traits.

An analysis of different CRISPR protospacer-adjacent motifs (PAMs) from various Cas enzymes shows that GC-rich PAMs are more abundant in genomic regions that harbour enriched heritability for human complex traits.

False-positive IRESes from Hoxa9 and other genes resulting from errors in mammalian 5' UTR annotations

Hyperconserved genomic sequences have great promise for understanding core biological processes. It has been recently proposed that scores of hyperconserved 5' untranslated regions (UTRs), also known as transcript leaders (hTLs), encode internal ribosome entry sites (IRESes) that drive cap-independent translation, in part, via interactions with ribosome expansion segments. However, the direct functional significance of such interactions has not yet been definitively demonstrated. We provide evidence that the putative IRESes previously reported in Hox gene hTLs are rarely included in transcript leaders. Instead, these regions function independently as transcriptional promoters. In addition, we find the proposed RNA structure of the putative Hoxa9 IRES is not conserved. Instead, sequences previously shown to be essential for putative IRES activity encode a hyperconserved transcription factor binding site (E-box) that contributes to its promoter activity and is bound by several transcription factors, including USF1 and USF2. Similar E-box sequences enhance the promoter activities of other putative Hoxa gene IRESes. Moreover, we provide evidence that the vast majority of hTLs with putative IRES activity overlap transcriptional promoters, enhancers, and 3' splice sites that are most likely responsible for their reported IRES activities. These results argue strongly against recently reported widespread IRES-like activities from hTLs and contradict proposed interactions between ribosomal expansion segment ES9S and putative IRESes. Furthermore, our work underscores the importance of accurate transcript annotations, controls in bicistronic reporter assays, and the power of synthesizing publicly available data from multiple sources.

An enhancer located in a Pde6c intron drives transient expression in the cone photoreceptors of developing mouse and human retinas

How cone photoreceptors are formed during retinal development is only partially known. This is in part because we do not fully understand the gene regulatory network responsible for cone genesis. We reasoned that cis-regulatory elements (enhancers) active in nascent cones would be regulated by the same upstream network that controls cone formation. To dissect this network, we searched for enhancers active in developing cones. By electroporating enhancer-driven fluorescent reporter plasmids, we observed that a sequence within an intron of the cone-specific Pde6c gene acted as an enhancer in developing mouse cones. Similar fluorescent reporter plasmids were used to generate stable transgenic human induced pluripotent stem cells that were then grown into three-dimensional human retinal organoids. These organoids contained fluorescently labeled cones, demonstrating that the Pde6c enhancer was also active in human cones. We observed that enhancer activity was transient and labeled a minor population of developing rod photoreceptors in both mouse and human systems. This cone-enriched pattern argues that the Pde6c enhancer is activated in cells poised between rod and cone fates. Additionally, it suggests that the Pde6c enhancer is activated by the same regulatory network that selects or stabilizes cone fate choice. To further understand this regulatory network, we identified essential enhancer sequence regions through a series of mutagenesis experiments. This suggested that the Pde6c enhancer was regulated by transcription factor binding at five or more locations. Binding site predictions implicated transcription factor families known to control photoreceptor formation and families not previously associated with cone development. These results provide a framework for deciphering the gene regulatory network that controls cone genesis in both human and mouse systems. Our new transgenic human stem cell lines provide a tool for determining which cone developmental mechanisms are shared and distinct between mice and humans.

The accessible promoter-mediated supplementary effect of host factors provides new insight into the tropism of SARS-CoV-2

In the past year, the rapid spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resulted in the worldwide coronavirus disease 2019 (COVID-19) pandemic. Yet our understanding of the SARS-CoV-2 tropism mechanism is still insufficient. In this study, we examined the chromatin accessibility at the promoters of host factor genes (ACE2, TMPRSS2, NRP1, BSG, CTSL, and FURIN) in 14 tissue types, 23 tumor types, and 189 cell lines. We showed that the promoters of ACE2 and TMPRSS2 were accessible in a tissue- and cell-specific pattern, which is accordant with previous clinical research on SARS-CoV-2 tropism. We were able to further verify that type I interferon (IFN) could induce angiotensin-converting enzyme 2 (ACE2) expression in Caco-2 cells by enhancing the binding of HNF1A, the transcription factor of ACE2, to ACE2 promoter without changing chromatin accessibility. We then performed transcription factor (TF)-gene interactions network and pathway analyses and discovered that the TFs regulating host factor genes are enriched in pathways associated with viral infection. Finally, we established a novel model that suggests that open chromatin at the promoter mediates the host factors’ supplementary effect and ensures SARS-CoV-2 entry. Our work uncovers the relationship between epigenetic regulation and SARS-CoV-2 tropism and provides clues for further investigation of COVID-19 pathogenesis.

CpG island reconfiguration for the establishment and synchronization of polycomb functions upon exit from naive pluripotency

Polycomb group (PcG) proteins are essential for post-implantation development by depositing repressive histone modifications at promoters, mainly CpG islands (CGIs), of developmental regulator genes. However, promoter PcG marks are erased after fertilization and de novo established in peri-implantation embryos, coinciding with the transition from naive to primed pluripotency. Nevertheless, the molecular basis for this establishment remains unknown. In this study, we show that the expression of the long KDM2B isoform (KDM2BLF), which contains the demethylase domain, is specifically induced at peri-implantation and that its H3K36me2 demethylase activity is required for PcG enrichment at CGIs. Moreover, KDM2BLF interacts with BRG1/BRM-associated factor (BAF) and stabilizes BAF occupancy at CGIs for subsequent gain of accessibility, which precedes PcG enrichment. Consistently, KDM2BLF inactivation results in significantly delayed post-implantation development. In summary, our data unveil dynamic chromatin configuration of CGIs during exit from naive pluripotency and provide a conceptual framework for the spatiotemporal establishment of PcG functions.

Global approaches for profiling transcription initiation

Transcription start site (TSS) selection influences transcript stability and translation as well as protein sequence. Alternative TSS usage is pervasive in organismal development, is a major contributor to transcript isoform diversity in humans, and is frequently observed in human diseases including cancer. In this review, we discuss the breadth of techniques that have been used to globally profile TSSs and the resulting insights into gene regulation, as well as future prospects in this area of inquiry.

Analyses of Promoter , Enhancer, and Nucleosome Organization in Mammalian Cells by MNase-Seq

MNase-Seq is a genome-wide procedure that allows mapping of DNA associated to nucleosomes following micrococcal nuclease digestion. It is a rapid and robust technology useful for the analysis of chromatin properties genome-wide at the resolution of mono-nucleosomes. Here, we describe how to produce high-resolution nucleosome maps of cells grown in suspension or adherent mammalian cells. After only three steps: nuclei or cell preparation, native MNase digestion and DNA purification, libraries for high-throughput sequencing can be prepared. Genome-wide nucleosome maps allow analyzing chromatin opening at promoters or enhancers, nucleosome displacement, or labile nucleosome occupancy depending on the digestion condition used. As presented, MNase-Seq is a versatile tool for investigating chromatin dynamics, regulation, and to define open chromatin regions of regulatory elements in mammalian genomes.

Maize decrease in DNA methylation 1 targets RNA-directed DNA methylation on active chromatin

DNA methylation plays vital roles in repressing transposable element activity and regulating gene expression. The chromatin-remodeling factor Decrease in DNA methylation 1 (DDM1) is crucial for maintaining DNA methylation across diverse plant species, and is required for RNA-directed DNA methylation (RdDM) to maintain mCHH islands in maize (Zea mays). However, the mechanisms by which DDM1 is involved in RdDM are not well understood. In this work, we used chromatin immunoprecipitation coupled with high-throughput sequencing to ascertain the genome-wide occupancy of ZmDDM1 in the maize genome. The results revealed that ZmDDM1 recognized an 8-bp-long GC-rich degenerate DNA sequence motif, which is enriched in transcription start sites and other euchromatic regions. Meanwhile, 24-nucleotide siRNAs and CHH methylation were delineated at the edge of ZmDDM1-occupied sites. ZmDDM1 co-purified with Argonaute 4 (ZmAGO4) proteins, providing further evidence that ZmDDM1 is a component of RdDM complexes in planta. Consistent with this, the vast majority of ZmDDM1-targeted regions co-localized with ZmAGO4-bound genomic sites. Overall, our results suggest a model that ZmDDM1 may be recruited to euchromatic regions via recognition of a GC-rich motif, thereby remodeling chromatin to provide access for RdDM activities in maize.

Nuclear FGFR1 Regulates Gene Transcription and Promotes Antiestrogen Resistance in ER+ Breast Cancer

PURPOSE

FGFR1 overexpression has been associated with endocrine resistance in ER⁺ breast cancer. We found FGFR1 localized in the nucleus of breast cancer cells in primary tumors resistant to estrogen suppression. We investigated a role of nuclear FGFR1 on gene transcription and antiestrogen resistance.

EXPERIMENTAL DESIGN

Tumors from patients treated with letrozole were subjected to Ki67 and FGFR1 IHC. MCF7 cells were transduced with FGFR1(SP-)(NLS) to promote nuclear FGFR1 overexpression. FGFR1 genomic activity in ER⁺/FGFR1-amplified breast cancer cells ± FOXA1 siRNA or ± the FGFR tyrosine kinase inhibitor (TKI) erdafitinib was examined by chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq). The nuclear and chromatin-bound FGFR1 interactome was investigated by mass spectrometry (MS).

RESULTS

High nuclear FGFR1 expression in ER⁺ primary tumors positively correlated with post-letrozole Ki67 values. Nuclear FGFR1 overexpression influenced gene transcription and promoted resistance to estrogen suppression and to fulvestrant in vivo. A gene expression signature induced by nuclear FGFR1 correlated with shorter survival in the METABRIC cohort of patients treated with antiestrogens. ChIP-Seq revealed FGFR1 occupancy at transcription start sites, overlapping with active transcription histone marks. MS analysis of the nuclear FGFR1 interactome identified phosphorylated RNA-Polymerase II and FOXA1, with FOXA1 RNAi impairing FGFR1 recruitment to chromatin. Treatment with erdafitinib did not impair nuclear FGFR1 translocation and genomic activity.

CONCLUSIONS

These data suggest nuclear FGFR1 contributes to endocrine resistance by modulating gene transcription in ER⁺ breast cancer. Nuclear FGFR1 activity was unaffected by FGFR TKIs, thus supporting the development of treatment strategies to inhibit nuclear FGFR1 in ER+/FGFR1 overexpressing breast cancer.

BANP opens chromatin and activates CpG-island-regulated genes

The majority of gene transcripts generated by RNA polymerase II in mammalian genomes initiate at CpG island (CGI) promoters^1,2, yet our understanding of their regulation remains limited. This is in part due to the incomplete information that we have on transcription factors, their DNA-binding motifs and which genomic binding sites are functional in any given cell type^3-5. In addition, there are orphan motifs without known binders, such as the CGCG element, which is associated with highly expressed genes across human tissues and enriched near the transcription start site of a subset of CGI promoters^6-8. Here we combine single-molecule footprinting with interaction proteomics to identify BTG3-associated nuclear protein (BANP) as the transcription factor that binds this element in the mouse and human genome. We show that BANP is a strong CGI activator that controls essential metabolic genes in pluripotent stem and terminally differentiated neuronal cells. BANP binding is repelled by DNA methylation of its motif in vitro and in vivo, which epigenetically restricts most binding to CGIs and accounts for differential binding at aberrantly methylated CGI promoters in cancer cells. Upon binding to an unmethylated motif, BANP opens chromatin and phases nucleosomes. These findings establish BANP as a critical activator of a set of essential genes and suggest a model in which the activity of CGI promoters relies on methylation-sensitive transcription factors that are capable of chromatin opening.

K-mer Content Changes with Node Degree in Promoter-Enhancer Network of Mouse ES Cells

Maps of Hi-C contacts between promoters and enhancers can be analyzed as networks, with cis-regulatory regions as nodes and their interactions as edges. We checked if in the published promoter-enhancer network of mouse embryonic stem (ES) cells the differences in the node type (promoter or enhancer) and the node degree (number of regions interacting with a given promoter or enhancer) are reflected by sequence composition or sequence similarity of the interacting nodes. We used counts of all k-mers (k = 4) to analyze the sequence composition and the Euclidean distance between the k-mer count vectors (k-mer distance) as the measure of sequence (dis)similarity. The results we obtained with 4-mers are interpretable in terms of dinucleotides. Promoters are GC-rich as compared to enhancers, which is known. Enhancers are enriched in scaffold/matrix attachment regions (S/MARs) patterns and depleted of CpGs. Furthermore, we show that promoters are more similar to their interacting enhancers than vice-versa. Most notably, in both promoters and enhancers, the GC content and the CpG count increase with the node degree. As a consequence, enhancers of higher node degree become more similar to promoters, whereas higher degree promoters become less similar to enhancers. We confirmed the key results also for human keratinocytes.

The Long Road to Understanding RNAPII Transcription Initiation and Related Syndromes

In eukaryotes, transcription of protein-coding genes requires the assembly at core promoters of a large preinitiation machinery containing RNA polymerase II (RNAPII) and general transcription factors (GTFs). Transcription is potentiated by regulatory elements called enhancers, which are recognized by specific DNA-binding transcription factors that recruit cofactors and convey, following chromatin remodeling, the activating cues to the preinitiation complex. This review summarizes nearly five decades of work on transcription initiation by describing the sequential recruitment of diverse molecular players including the GTFs, the Mediator complex, and DNA repair factors that support RNAPII to enable RNA synthesis. The elucidation of the transcription initiation mechanism has greatly benefited from the study of altered transcription components associated with human diseases that could be considered transcription syndromes.

Characterizing Promoter and Enhancer Sequences by a Deep Learning Method

Promoters and enhancers are well-known regulatory elements modulating gene expression. As confirmed by high-throughput sequencing technologies, these regulatory elements are bidirectionally transcribed. That is, promoters produce stable mRNA in the sense direction and unstable RNA in the antisense direction, while enhancers transcribe unstable RNA in both directions. Although it is thought that enhancers and promoters share a similar architecture of transcription start sites (TSSs), how the transcriptional machinery distinctly uses these genomic regions as promoters or enhancers remains unclear. To address this issue, we developed a deep learning (DL) method by utilizing a convolutional neural network (CNN) and the saliency algorithm. In comparison with other classifiers, our CNN presented higher predictive performance, suggesting the overarching importance of the high-order sequence features, captured by the CNN. Moreover, our method revealed that there are substantial sequence differences between the enhancers and promoters. Remarkably, the 20–120 bp downstream regions from the center of bidirectional TSSs seemed to contribute to the RNA stability. These regions in promoters tend to have a larger number of guanines and cytosines compared to those in enhancers, and this feature contributed to the classification of the regulatory elements. Our CNN-based method can capture the complex TSS architectures. We found that the genomic regions around TSSs for promoters and enhancers contribute to RNA stability and show GC-biased characteristics as a critical determinant for promoter TSSs.

Not just a writer: PRC2 as a chromatin reader

PRC2 deposits the H3K27me3 repressive mark, which facilitates transcription repression of developmental genes. The decision of whether a particular gene is silenced at a given point during development is heavily dependent on the chromatin context. More than just a simple epigenetic writer, PRC2 employs several distinct chromatin reading capabilities to sense the local chromatin environment and modulate the H3K27me3 writer activity in a context-dependent manner. Here we discuss the complex interplay of PRC2 with the hallmarks of active and repressive chromatin, how it affects H3K27me3 deposition and how it guides transcriptional activity.

Prediction of disease-associated functional variants in noncoding regions through a comprehensive analysis by integrating datasets and features

One of the greatest challenges in human genetics is deciphering the link between functional variants in noncoding sequences and the pathophysiology of complex diseases. To address this issue, many methods have been developed to sort functional single-nucleotide variants (SNVs) for neutral SNVs in noncoding regions. In this study, we integrated well-established features and commonly used datasets and merged them into large-scale datasets based on a random forest model, which yielded promising performance and outperformed some cutting-edge approaches. Our analyses of feature importance and data coverage also provide certain clues for future research in enhancing the prediction of functional noncoding SNVs.

Androgen and glucocorticoid receptor direct distinct transcriptional programs by receptor-specific and shared DNA binding sites

The glucocorticoid (GR) and androgen (AR) receptors execute unique functions in vivo, yet have nearly identical DNA binding specificities. To identify mechanisms that facilitate functional diversification among these transcription factor paralogs, we studied them in an equivalent cellular context. Analysis of chromatin and sequence suggest that divergent binding, and corresponding gene regulation, are driven by different abilities of AR and GR to interact with relatively inaccessible chromatin. Divergent genomic binding patterns can also be the result of subtle differences in DNA binding preference between AR and GR. Furthermore, the sequence composition of large regions (>10 kb) surrounding selectively occupied binding sites differs significantly, indicating a role for the sequence environment in guiding AR and GR to distinct binding sites. The comparison of binding sites that are shared shows that the specificity paradox can also be resolved by differences in the events that occur downstream of receptor binding. Specifically, shared binding sites display receptor-specific enhancer activity, cofactor recruitment and changes in histone modifications. Genomic deletion of shared binding sites demonstrates their contribution to directing receptor-specific gene regulation. Together, these data suggest that differences in genomic occupancy as well as divergence in the events that occur downstream of receptor binding direct functional diversification among transcription factor paralogs.

The "Genomic Code": DNA Pervasively Moulds Chromatin Structures Leaving no Room for "Junk"

The chromatin of the human genome was analyzed at three DNA size levels. At the first, compartment level, two "gene spaces" were found many years ago: A GC-rich, gene-rich "genome core" and a GC-poor, gene-poor "genome desert", the former corresponding to open chromatin centrally located in the interphase nucleus, the latter to closed chromatin located peripherally. This bimodality was later confirmed and extended by the discoveries (1) of LADs, the Lamina-Associated Domains, and InterLADs; (2) of two "spatial compartments", A and B, identified on the basis of chromatin interactions; and (3) of "forests and prairies" characterized by high and low CpG islands densities. Chromatin compartments were shown to be associated with the compositionally different, flat and single- or multi-peak DNA structures of the two, GC-poor and GC-rich, "super-families" of isochores. At the second, sub-compartment, level, chromatin corresponds to flat isochores and to isochore loops (due to compositional DNA gradients) that are susceptible to extrusion. Finally, at the short-sequence level, two sets of sequences, GC-poor and GC-rich, define two different nucleosome spacings, a short one and a long one. In conclusion, chromatin structures are moulded according to a "genomic code" by DNA sequences that pervade the genome and leave no room for "junk".

ATRX limits the accessibility of histone H3-occupied HSV genomes during lytic infection

Histones are rapidly loaded on the HSV genome upon entry into the nucleus of human fibroblasts, but the effects of histone loading on viral replication have not been fully defined. We showed recently that ATRX is dispensable for de novo deposition of H3 to HSV genomes after nuclear entry but restricted infection through maintenance of viral heterochromatin. To further investigate the roles that ATRX and other histone H3 chaperones play in restriction of HSV, we infected human fibroblasts that were systematically depleted of nuclear H3 chaperones. We found that the ATRX/DAXX complex is unique among nuclear H3 chaperones in its capacity to restrict ICP0-null HSV infection. Only depletion of ATRX significantly alleviated restriction of viral replication. Interestingly, no individual nuclear H3 chaperone was required for deposition of H3 onto input viral genomes, suggesting that during lytic infection, H3 deposition may occur through multiple pathways. ChIP-seq for total histone H3 in control and ATRX-KO cells infected with ICP0-null HSV showed that HSV DNA is loaded with high levels of histones across the entire viral genome. Despite high levels of H3, ATAC-seq analysis revealed that HSV DNA is highly accessible, especially in regions of high GC content, and is not organized largely into ordered nucleosomes during lytic infection. ATRX reduced accessibility of viral DNA to the activity of a TN5 transposase and enhanced accumulation of viral DNA fragment sizes associated with nucleosome-like structures. Together, these findings support a model in which ATRX restricts viral infection by altering the structure of histone H3-loaded viral chromatin that reduces viral DNA accessibility for transcription. High GC rich regions of the HSV genome, especially the S component inverted repeats of the HSV-1 genome, show increased accessibility, which may lead to increased ability to transcribe the IE genes encoded in these regions during initiation of infection.

Species-Specific Relationships between DNA and Chromatin Properties of CpG Islands in Embryonic Stem Cells and Differentiated Cells

CpG islands often exhibit low DNA methylation, high histone H3 lysine 4 trimethylation, low nucleosome density, and high DNase I hypersensitivity, yet the rules by which CpG islands are sensed remain poorly understood. In this study, we first evaluated the relationships between the DNA and the chromatin properties of CpG islands in embryonic stem cells using modified bacterial artificial chromosomes. Then, using a bioinformatic approach, we identified strict CpG-island density and length thresholds in mouse embryonic stem and differentiated cells that consistently specify low DNA methylation levels. Surprisingly, the human genome exhibited a dramatically different relationship between DNA properties and DNA methylation levels of CpG islands. Further analysis allowed speculation that this difference is accommodated in part by evolutionary changes in the nucleotide composition of orthologous promoters. Thus, a change in the rules by which CpG-island properties are sensed may have co-evolved with compensatory genome adaptation events during mammalian evolution.

•

Strict DNA property rules can predict low DNA methylation at murine CpG islands

•

Mice and humans differ greatly in how DNA properties influence DNA methylation

•

Evolutionary adaptation of CpG islands compensates for the species differences

Signaling-to-chromatin pathways in the immune system

Complex organisms are able to respond to diverse environmental cues by rapidly inducing specific transcriptional programs comprising a few dozen genes among thousands. The highly complex environment within the nucleus-a crowded milieu containing large genomes tightly condensed with histone proteins in the form of chromatin-makes inducible transcription a challenge for the cell, akin to the proverbial needle in a haystack. The different signaling pathways and transcription factors involved in the transmission of information from the cell surface to the nucleus have been readily explored, but not so much the specific mechanisms employed by the cell to ultimately instruct the chromatin changes necessary for a fast and robust transcription activation. Signaling pathways rely on cascades of protein kinases that, in addition to activating transcription factors can also activate the chromatin template by phosphorylating histone proteins, what we refer to as "signaling-to-chromatin." These pathways appear to be selectively employed and especially critical for driving inducible transcription in macrophages and likely in diverse other immune cell populations. Here, we discuss signaling-to-chromatin pathways with potential relevance in diverse immune cell populations together with chromatin related mechanisms that help to "solve" the needle in a haystack challenge of robust chromatin activation and inducible transcription.

Structural dynamics of DNA depending on methylation pattern

DNA methylation is associated with a number of biological phenomena, and plays crucial roles in epigenetic regulation of eukaryotic gene expression. It is also suggested that DNA methylation alters the mechanical properties of DNA molecules, which is likely to affect epigenetic regulation. However, it has not been systematically investigated how methylation changes the structural and dynamic features of DNA. In this research, to elucidate the effects of methylation on DNA mechanics, a fully atomic molecular dynamics simulation of double-stranded DNA with several methylation patterns was performed. Through the analysis of the relative positioning of the nucleotides (base-step variables), characteristic changes in terms of local flexibility were observed, which further affected the overall DNA geometry and stiffness. These findings may serve as a basis for a discussion on methylation-dependent DNA dynamics in physiological conditions.

The impact of different negative training data on regulatory sequence predictions

Regulatory regions, like promoters and enhancers, cover an estimated 5–15% of the human genome. Changes to these sequences are thought to underlie much of human phenotypic variation and a substantial proportion of genetic causes of disease. However, our understanding of their functional encoding in DNA is still very limited. Applying machine or deep learning methods can shed light on this encoding and gapped k-mer support vector machines (gkm-SVMs) or convolutional neural networks (CNNs) are commonly trained on putative regulatory sequences. Here, we investigate the impact of negative sequence selection on model performance. By training gkm-SVM and CNN models on open chromatin data and corresponding negative training dataset, both learners and two approaches for negative training data are compared. Negative sets use either genomic background sequences or sequence shuffles of the positive sequences. Model performance was evaluated on three different tasks: predicting elements active in a cell-type, predicting cell-type specific elements, and predicting elements' relative activity as measured from independent experimental data. Our results indicate strong effects of the negative training data, with genomic backgrounds showing overall best results. Specifically, models trained on highly shuffled sequences perform worse on the complex tasks of tissue-specific activity and quantitative activity prediction, and seem to learn features of artificial sequences rather than regulatory activity. Further, we observe that insufficient matching of genomic background sequences results in model biases. While CNNs achieved and exceeded the performance of gkm-SVMs for larger training datasets, gkm-SVMs gave robust and best results for typical training dataset sizes without the need of hyperparameter optimization.

CpG Islands Shape the Epigenome Landscape

Epigenetic modifications and nucleosome positioning play an important role in modulating gene expression. However, how the patterns of epigenetic modifications and nucleosome positioning are established around promoters is not well understood. Here, we have addressed these questions in a series of genome-wide experiments coupled to a novel bioinformatic analysis approach. Our data reveal a clear correlation between CpG density, promoter activity and accumulation of active or repressive histone marks. CGI boundaries define the chromatin promoter regions that will be epigenetically modified. CpG-rich promoters are targeted by histone modifications and histone variants, while CpG-poor promoters are regulated by DNA methylation. CGIs boundaries, but not transcriptional activity, are essential determinants of H2A.Z positioning in vicinity of the promoters, suggesting that the presence of H2A.Z is not related to transcriptional control. Accordingly, H2A.Z depletion has no impact on gene expression of arrested mouse embryonic fibroblasts. Therefore, the underlying DNA sequence, the promoter CpG density and, to a lesser extent, transcriptional activity, are key factors implicated in promoter chromatin architecture.

On the irrelevancy of hydroxyl radical to DNA damage from oxidative stress and implications for epigenetics

Contrary to frequent reports in the literature, hydroxyl radical is not a key species participating in endogenous oxidative DNA damage. Instead, carbonate radical anion is formed from the Fenton reaction under cellular conditions and from decomposition of nitrosoperoxycarbonate generated during inflammation. Carbonate radical anion is a potent one-electron oxidant capable of generating base radical cations that can migrate over long distances in duplex DNA, ultimately generating 8-oxo-7,8-dihydroguanine at a redox-sensitive sequence such as GGG. Such a mechanism enables G-quadruplex-forming sequences to act as long-range sensors of oxidative stress, impacting gene expression via the DNA repair mechanism that reads and ultimately erases the oxidized base. With a writing, reading and erasing mechanism in place, oxidative 'damage' to DNA might be relabeled as 'epigenetic' modifications.

The landscape of accessible chromatin in bovine oocytes and early embryos

Chromatin reorganization governs the regulation of gene expression during preimplantation development. However, the landscape of chromatin dynamics in this period has not been explored in bovine. In this study, we constructed a genome-wide map of accessible chromatin in bovine oocytes and early embryos using an improved assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) which revealed unique features of the accessible chromatin during bovine early embryo development. We found that chromatin accessibility is low in oocytes and 2-/4-cell embryos, followed by a significant increase in embryos during major embryonic genome activation (EGA), and peaked in elongating day 14 embryos. Genome-wide characteristics of open chromatin showed that ATAC-seq signals in both transcription start sites (TSS) and transcription end sites (TES) were strong. Additionally, the distal ATAC-seq peaks were enriched in repeat elements in a type-specific and stage-specific manner. We further unveiled a series of transcription factor (TF) motifs with distinct variation of enrichment from distal ATAC-seq peaks. By integrated analysis of chromatin accessibility with transcriptomes and DNA methylomes in bovine early embryos, we showed that promoter accessibility was positively correlated with gene expression, especially during major EGA, and was strongly correlated to DNA methylation and CpG density. Finally, we identified the critical chromatin signatures and TFs that differ between in vivo and in vitro derived blastocysts, which provides insights to the potential mechanisms leading to low quality of embryos produced in vitro. Together, this comprehensive analysis revealed critical features of chromatin landscape and epigenetic reprogramming during bovine preimplantation embryo development.

Simple and efficient profiling of transcription initiation and transcript levels with STRIPE-seq

Accurate mapping of transcription start sites (TSSs) is key for understanding transcriptional regulation. However, current protocols for genome-wide TSS profiling are laborious and/or expensive. We present Survey of TRanscription Initiation at Promoter Elements with high-throughput sequencing (STRIPE-seq), a simple, rapid, and cost-effective protocol for sequencing capped RNA 5' ends from as little as 50 ng total RNA. Including depletion of uncapped RNA and reaction cleanups, a STRIPE-seq library can be constructed in about 5 h. We show application of STRIPE-seq to TSS profiling in yeast and human cells and show that it can also be effectively used for quantification of transcript levels and analysis of differential gene expression. In conjunction with our ready-to-use computational workflows, STRIPE-seq is a straightforward, efficient means by which to probe the landscape of transcriptional initiation.

Role of cell-type specific nucleosome positioning in inducible activation of mammalian promoters

The organization of nucleosomes across functional genomic elements represents a critical layer of control. Here, we present a strategy for high-resolution nucleosome profiling at selected genomic features, and use this to analyse dynamic nucleosome positioning at inducible and cell-type-specific mammalian promoters. We find that nucleosome patterning at inducible promoters frequently resembles that at active promoters, even before stimulus-driven activation. Accordingly, the nucleosome profile at many inactive inducible promoters is sufficient to predict cell-type-specific responsiveness. Induction of gene expression is generally not associated with major changes to nucleosome patterning, and a subset of inducible promoters can be activated without stable nucleosome depletion from their transcription start sites. These promoters are generally dependent on remodelling enzymes for their inducible activation, and exhibit transient nucleosome depletion only at alleles undergoing transcription initiation. Together, these data reveal how the responsiveness of inducible promoters to activating stimuli is linked to cell-type-specific nucleosome patterning.

Nucleosome organisation plays important roles in regulating functional genomic elements. Here, the authors use high-resolution profiling to analyse dynamic nucleosome positioning at inducible and cell-type-specific promoters, providing a global view of chromatin architecture at inducible promoters.

Sensitive Automated Measurement of Histone-DNA Affinities in Nucleosomes

The DNA of eukaryotes is wrapped around histone octamers to form nucleosomes. Although it is well established that the DNA sequence significantly influences nucleosome formation, its precise contribution has remained controversial, partially owing to the lack of quantitative affinity data. Here, we present a method to measure DNA-histone binding free energies at medium throughput and with high sensitivity. Competitive nucleosome formation is achieved through automation, and a modified epifluorescence microscope is used to rapidly and accurately measure the fractions of bound/unbound DNA based on fluorescence anisotropy. The procedure allows us to obtain full titration curves with high reproducibility. We applied this technique to measure the histone-DNA affinities for 47 DNA sequences and analyzed how the affinities correlate with relevant DNA sequence features. We found that the GC content has a significant impact on nucleosome-forming preferences, but 10 bp dinucleotide periodicities and the presence of poly(dA:dT) stretches do not.

•

Robotics permits full titration series to measure histone-DNA binding affinities

•

Fluorescence anisotropy used as a fast, sensitive readout of bound/unbound DNA

•

Free energies span three orders of magnitude, less for naturally occurring sequences

•

GC content is a major determinant of measured binding free energies

Chromatin accessibility established by Pou5f3, Sox19b and Nanog primes genes for activity during zebrafish genome activation

In many organisms, early embryonic development is driven by maternally provided factors until the controlled onset of transcription during zygotic genome activation. The regulation of chromatin accessibility and its relationship to gene activity during this transition remain poorly understood. Here, we generated chromatin accessibility maps with ATAC-seq from genome activation until the onset of lineage specification. During this period, chromatin accessibility increases at regulatory elements. This increase is independent of RNA polymerase II-mediated transcription, with the exception of the hypertranscribed miR-430 locus. Instead, accessibility often precedes the transcription of associated genes. Loss of the maternal transcription factors Pou5f3, Sox19b, and Nanog, which are known to be required for zebrafish genome activation, results in decreased accessibility at regulatory elements. Importantly, the accessibility of regulatory regions, especially when established by Pou5f3, Sox19b and Nanog, is predictive for future transcription. Our results show that the maternally provided transcription factors Pou5f3, Sox19b, and Nanog open up chromatin and prime genes for activity during zygotic genome activation in zebrafish.