Mechanisms by which transcription factors gain access to target sequence elements in chromatin

Chance promoter activities illuminate the origins of eukaryotic intergenic transcriptions

It is debated whether the pervasive intergenic transcription from eukaryotic genomes has functional significance or simply reflects the promiscuity of RNA polymerases. We approach this question by comparing chance promoter activities with the expression levels of intergenic regions in the model eukaryote Saccharomyces cerevisiae. We build a library of over 105 strains, each carrying a 120-nucleotide, chromosomally integrated, completely random sequence driving the potential transcription of a barcode. Quantifying the RNA concentration of each barcode in two environments reveals that 41-63% of random sequences have significant, albeit usually low, promoter activities. Therefore, even in eukaryotes, where the presence of chromatin is thought to repress transcription, chance transcription is prevalent. We find that only 1-5% of yeast intergenic transcriptions are unattributable to chance promoter activities or neighboring gene expressions, and these transcriptions exhibit higher-than-expected environment-specificity. These findings suggest that only a minute fraction of intergenic transcription is functional in yeast.

Nuclear export restricts Gdown1 to a mitotic function

Approximately half of purified mammalian RNA polymerase II (Pol II) is associated with a tightly interacting sub-stoichiometric subunit, Gdown1. Previous studies have established that Gdown1 inhibits transcription initiation through competitive interactions with general transcription factors and blocks the Pol II termination activity of transcription termination factor 2 (TTF2). However, the biological functions of Gdown1 remain poorly understood. Here, we utilized genetic, microscopic, and multi-omics approaches to functionally characterize Gdown1 in three human cell lines. Acute depletion of Gdown1 caused minimal direct effects on transcription. We show that Gdown1 resides predominantly in the cytoplasm of interphase cells, shuttles between the cytoplasm and nucleus, and is regulated by nuclear export. Gdown1 enters the nucleus at the onset of mitosis. Consistently, genetic ablation of Gdown1 is associated with partial de-repression of mitotic transcription, and Gdown1 KO cells present with evidence of aberrant mitoses coupled to p53 pathway activation. Evidence is presented demonstrating that Gdown1 modulates the combined functions of purified productive elongation factors PAF1C, RTF1, SPT6, DSIF and P-TEFb in vitro. Collectively, our findings support a model wherein the Pol II-regulatory function of Gdown1 occurs during mitosis and is required for genome integrity.

MNase Digestion Protection Patterns of the Linker DNA in Chromatosomes

The compact nucleosomal structure limits DNA accessibility and regulates DNA-dependent cellular activities. Linker histones bind to nucleosomes and compact nucleosomal arrays into a higher-order chromatin structure. Recent developments in high throughput technologies and structural computational studies provide nucleosome positioning at a high resolution and contribute to the information of linker histone location within a chromatosome. However, the precise linker histone location within the chromatin fibre remains unclear. Using monomer extension, we mapped core particle and chromatosomal positions over a core histone-reconstituted, 1.5 kb stretch of DNA from the chicken adult β-globin gene, after titration with linker histones and linker histone globular domains. Our results show that, although linker histone globular domains and linker histones display a wide variation in their binding affinity for different positioned nucleosomes, they do not alter nucleosome positions or generate new nucleosome positions. Furthermore, the extra ~20 bp of DNA protected in a chromatosome is usually symmetrically distributed at each end of the core particle, suggesting linker histones or linker histone globular domains are located close to the nucleosomal dyad axis.

Histone H3K4me1 and H3K27ac play roles in nucleosome eviction and eRNA transcription, respectively, at enhancers

Histone H3K4me1 and H3K27ac are enhancer-specific modifications and are required for enhancers to activate transcription of target genes. However, the reciprocal effects of these histone modifications on each other and their roles in enhancers are not clear. Here to comparatively analyze the role of these modifications, we inhibited H3K4me1 and H3K27ac by deleting the SET domains of histone methyltransferases MLL3 and MLL4 and the HAT domain of histone acetyltransferase p300, respectively, in erythroid K562 cells. The loss of H3K4me1 reduced H3K27ac at the β-globin enhancer LCR HSs, but H3K27ac reduction did not affect H3K4me1. This unequal relationship between two modifications was revealed in putative enhancers by genome-wide analysis using ChIP-seq. Histone H3 eviction at putative enhancers was weakened by the loss of H3K4me1 but not by the loss of H3K27ac. Chromatin remodeling complexes were recruited into the β-globin LCR HSs in a H3K4me1-dependent manner. In contrast, H3K27ac was required for enhancer RNA (eRNA) transcription, and H3K4me1 was not enough for it. Forced H3K27ac-induced eRNA transcription without affecting H3K4me1 at the β-globin LCR HSs. These results indicate that H3K4me1 and H3K27ac affect each other in different ways and play more direct roles in nucleosome eviction and eRNA transcription, respectively, at enhancers.

Nuclear Organization during Hepatogenesis in Zebrafish Requires Uhrf1

Acquisition of cellular fate during development is initiated and maintained by well-coordinated patterns of gene expression that are dictated by the epigenetic landscape and genome organization in the nucleus. While the epigenetic marks that mediate developmental gene expression patterns during organogenesis have been well studied, less is known about how epigenetic marks influence nuclear organization during development. This study examines the relationship between nuclear structure, chromatin accessibility, DNA methylation, and gene expression during hepatic outgrowth in zebrafish larvae. We investigate the relationship between these features using mutants that lack DNA methylation. Hepatocyte nuclear morphology was established coincident with hepatocyte differentiation at 80 h post-fertilization (hpf), and nuclear shape and size continued to change until the conclusion of outgrowth and morphogenesis at 120 hpf. Integrating ATAC-Seq analysis with DNA methylation profiling of zebrafish livers at 120 hpf showed that closed and highly methylated chromatin occupies most transposable elements and that open chromatin correlated with gene expression. DNA hypomethylation, due to mutation of genes encoding ubiquitin-like, containing PHD and RING Finger Domains 1 (uhrf1) and DNA methyltransferase (dnmt1), did not block hepatocyte differentiation, but had dramatic effects on nuclear organization. Hepatocytes in uhrf1 mutants have large, deformed nuclei with multiple nucleoli, downregulation of nucleolar genes, and a complete lack of the nuclear lamina. Loss of lamin B2 staining was phenocopied by dnmt1 mutation. Together, these data show that hepatocyte nuclear morphogenesis coincides with organ morphogenesis and outgrowth, and that DNA methylation directs chromatin organization, and, in turn, hepatocyte nuclear shape and size during liver development.

Chromatin Profiling Techniques: Exploring the Chromatin Environment and Its Contributions to Complex Traits

The genetic architecture of complex traits is multifactorial. Genome-wide association studies (GWASs) have identified risk loci for complex traits and diseases that are disproportionately located at the non-coding regions of the genome. On the other hand, we have just begun to understand the regulatory roles of the non-coding genome, making it challenging to precisely interpret the functions of non-coding variants associated with complex diseases. Additionally, the epigenome plays an active role in mediating cellular responses to fluctuations of sensory or environmental stimuli. However, it remains unclear how exactly non-coding elements associate with epigenetic modifications to regulate gene expression changes and mediate phenotypic outcomes. Therefore, finer interrogations of the human epigenomic landscape in associating with non-coding variants are warranted. Recently, chromatin-profiling techniques have vastly improved our understanding of the numerous functions mediated by the epigenome and DNA structure. Here, we review various chromatin-profiling techniques, such as assays of chromatin accessibility, nucleosome distribution, histone modifications, and chromatin topology, and discuss their applications in unraveling the brain epigenome and etiology of complex traits at tissue homogenate and single-cell resolution. These techniques have elucidated compositional and structural organizing principles of the chromatin environment. Taken together, we believe that high-resolution epigenomic and DNA structure profiling will be one of the best ways to elucidate how non-coding genetic variations impact complex diseases, ultimately allowing us to pinpoint cell-type targets with therapeutic potential.

Local states of chromatin compaction at transcription start sites control transcription levels

The ‘open’ and ‘compact’ regions of chromatin are considered to be regions of active and silent transcription, respectively. However, individual genes produce transcripts at different levels, suggesting that transcription output does not depend on the simple open-compact conversion of chromatin, but on structural variations in chromatin itself, which so far have remained elusive. In this study, weakly crosslinked chromatin was subjected to sedimentation velocity centrifugation, which fractionated the chromatin according to its degree of compaction. Open chromatin remained in upper fractions, while compact chromatin sedimented to lower fractions depending on the level of nucleosome assembly. Although nucleosomes were evenly detected in all fractions, histone H1 was more highly enriched in the lower fractions. H1 was found to self-associate and crosslinked to histone H3, suggesting that H1 bound to H3 interacts with another H1 in an adjacent nucleosome to form compact chromatin. Genome-wide analyses revealed that nearly the entire genome consists of compact chromatin without differences in compaction between repeat and non-repeat sequences; however, active transcription start sites (TSSs) were rarely found in compact chromatin. Considering the inverse correlation between chromatin compaction and RNA polymerase binding at TSSs, it appears that local states of chromatin compaction determine transcription levels.

Mobility connects: transposable elements wire new transcriptional networks by transferring transcription factor binding motifs

Transposable elements (TEs) constitute major fractions of plant genomes. Their potential to be mobile provides them with the capacity to cause major genome rearrangements. Those effects are potentially deleterious and enforced the evolution of epigenetic suppressive mechanisms controlling TE activity. However, beyond their deleterious effects, TE insertions can be neutral or even advantageous for the host, leading to long-term retention of TEs in the host genome. Indeed, TEs are increasingly recognized as major drivers of evolutionary novelties by regulating the expression of nearby genes. TEs frequently contain binding motifs for transcription factors and capture binding motifs during transposition, which they spread through the genome by transposition. Thus, TEs drive the evolution and diversification of gene regulatory networks by recruiting lineage-specific targets under the regulatory control of specific transcription factors. This process can explain the rapid and repeated evolution of developmental novelties, such as C₄ photosynthesis and a wide spectrum of stress responses in plants. It also underpins the convergent evolution of embryo nourishing tissues, the placenta in mammals and the endosperm in flowering plants. Furthermore, the gene regulatory network underlying flower development has also been largely reshaped by TE-mediated recruitment of regulatory elements; some of them being preserved across long evolutionary timescales. In this review, we highlight the potential role of TEs as evolutionary toolkits in plants by showcasing examples of TE-mediated evolutionary novelties.

RNA aptamer capture of macromolecular complexes for mass spectrometry analysis

Specific genomic functions are dictated by macromolecular complexes (MCs) containing multiple proteins. Affinity purification of these complexes, often using antibodies, followed by mass spectrometry (MS) has revolutionized our ability to identify the composition of MCs. However, conventional immunoprecipitations suffer from contaminating antibody/serum-derived peptides that limit the sensitivity of detection for low-abundant interacting partners using MS. Here, we present AptA-MS (aptamer affinity-mass spectrometry), a robust strategy primarily using a specific, high-affinity RNA aptamer against Green Fluorescent Protein (GFP) to identify interactors of a GFP-tagged protein of interest by high-resolution MS. Utilizing this approach, we have identified the known molecular chaperones that interact with human Heat Shock Factor 1 (HSF1), and observed an increased association with several proteins upon heat shock, including translation elongation factors and histones. HSF1 is known to be regulated by multiple post-translational modifications (PTMs), and we observe both known and new sites of modifications on HSF1. We show that AptA-MS provides a dramatic target enrichment and detection sensitivity in evolutionarily diverse organisms and allows identification of PTMs without the need for modification-specific enrichments. In combination with the expanding libraries of GFP-tagged cell lines, this strategy offers a general, inexpensive, and high-resolution alternative to conventional approaches for studying MCs.

Developmental Transcriptional Enhancers: A Subtle Interplay between Accessibility and Activity: Considering Quantitative Accessibility Changes between Different Regulatory States of an Enhancer Deconvolutes the Complex Relationship between Accessibility and Activity

Measurements of open chromatin in specific cell types are widely used to infer the spatiotemporal activity of transcriptional enhancers. How reliable are these predictions? In this review, it is argued that the relationship between the accessibility and activity of an enhancer is insufficiently described by simply considering open versus closed chromatin, or active versus inactive enhancers. Instead, recent studies focusing on the quantitative nature of accessibility signal reveal subtle differences between active enhancers and their different inactive counterparts: the closed silenced state and the accessible primed and repressed states. While the open structure as such is not a specific indicator of enhancer activity, active enhancers display a higher degree of accessibility than the primed and repressed states. Molecular mechanisms that may account for these quantitative differences are discussed. A model that relates molecular events at an enhancer to changes in its activity and accessibility in a developing tissue is also proposed.

Zscan4 binds nucleosomal microsatellite DNA and protects mouse two-cell embryos from DNA damage

Zinc finger protein Zscan4 is selectively expressed in mouse two-cell (2C) embryos undergoing zygotic genome activation (ZGA) and in a rare subpopulation of embryonic stem cells with 2C-like features. Here, we show that Zscan4 specifically recognizes a subset of (CA)n microsatellites, repeat sequences prone to genomic instability. Zscan4-associated microsatellite regions are characterized by low nuclease sensitivity and high histone occupancy. In vitro, Zscan4 binds nucleosomes and protects them from disassembly upon torsional strain. Furthermore, Zscan4 depletion leads to elevated DNA damage in 2C mouse embryos in a transcription-dependent manner. Together, our results identify Zscan4 as a DNA sequence-dependent microsatellite binding factor and suggest a developmentally regulated mechanism, which protects fragile genomic regions from DNA damage at a time of embryogenesis associated with high transcriptional burden and genomic stress.

Common Regulatory Targets of NFIA, NFIX and NFIB during Postnatal Cerebellar Development

Transcriptional regulation plays a central role in controlling neural stem and progenitor cell proliferation and differentiation during neurogenesis. For instance, transcription factors from the nuclear factor I (NFI) family have been shown to co-ordinate neural stem and progenitor cell differentiation within multiple regions of the embryonic nervous system, including the neocortex, hippocampus, spinal cord and cerebellum. Knockout of individual Nfi genes culminates in similar phenotypes, suggestive of common target genes for these transcription factors. However, whether or not the NFI family regulates common suites of genes remains poorly defined. Here, we use granule neuron precursors (GNPs) of the postnatal murine cerebellum as a model system to analyse regulatory targets of three members of the NFI family: NFIA, NFIB and NFIX. By integrating transcriptomic profiling (RNA-seq) of Nfia- and Nfix-deficient GNPs with epigenomic profiling (ChIP-seq against NFIA, NFIB and NFIX, and DNase I hypersensitivity assays), we reveal that these transcription factors share a large set of potential transcriptional targets, suggestive of complementary roles for these NFI family members in promoting neural development.

The Role of Chromatin Accessibility in cis-Regulatory Evolution

Transcription factor (TF) binding is determined by sequence as well as chromatin accessibility. Although the role of accessibility in shaping TF-binding landscapes is well recorded, its role in evolutionary divergence of TF binding, which in turn can alter cis-regulatory activities, is not well understood. In this work, we studied the evolution of genome-wide binding landscapes of five major TFs in the core network of mesoderm specification, between Drosophila melanogaster and Drosophila virilis, and examined its relationship to accessibility and sequence-level changes. We generated chromatin accessibility data from three important stages of embryogenesis in both Drosophila melanogaster and Drosophila virilis and recorded conservation and divergence patterns. We then used multivariable models to correlate accessibility and sequence changes to TF-binding divergence. We found that accessibility changes can in some cases, for example, for the master regulator Twist and for earlier developmental stages, more accurately predict binding change than is possible using TF-binding motif changes between orthologous enhancers. Accessibility changes also explain a significant portion of the codivergence of TF pairs. We noted that accessibility and motif changes offer complementary views of the evolution of TF binding and developed a combined model that captures the evolutionary data much more accurately than either view alone. Finally, we trained machine learning models to predict enhancer activity from TF binding and used these functional models to argue that motif and accessibility-based predictors of TF-binding change can substitute for experimentally measured binding change, for the purpose of predicting evolutionary changes in enhancer activity.

Clock-controlled rhythmic transcription: is the clock enough and how does it work?

Circadian clocks regulate the rhythmic expression of thousands of genes underlying the daily oscillations of biological functions. Here, we discuss recent findings showing that circadian clock rhythmic transcriptional outputs rely on additional mechanisms than just clock gene DNA binding, which may ultimately contribute to the plasticity of circadian transcriptional programs.

PAX8 activates metabolic genes via enhancer elements in Renal Cell Carcinoma

Transcription factor networks shape the gene expression programs responsible for normal cell identity and pathogenic state. Using Core Regulatory Circuitry analysis (CRC), we identify PAX8 as a candidate oncogene in Renal Cell Carcinoma (RCC) cells. Validation of large-scale functional genomic screens confirms that PAX8 silencing leads to decreased proliferation of RCC cell lines. Epigenomic analyses of PAX8-dependent cistrome demonstrate that PAX8 largely occupies active enhancer elements controlling genes involved in various metabolic pathways. We selected the ferroxidase Ceruloplasmin (CP) as an exemplary gene to dissect PAX8 molecular functions. PAX8 recruits histone acetylation activity at bound enhancers looping onto the CP promoter. Importantly, CP expression correlates with sensitivity to PAX8 silencing and identifies a subset of RCC cases with poor survival. Our data identifies PAX8 as a candidate oncogene in RCC and provides a potential biomarker to monitor its activity.

Transcription factors are critical regulators of cell identity. Here, the authors use computational and functional genomic approaches to show an oncogenic role of PAX8 in renal cancer. Mechanistic dissection of PAX8 functions reveal its role in activating genes associated with metabolic pathways.

LeNup: learning nucleosome positioning from DNA sequences with improved convolutional neural networks

Motivation

Nucleosome positioning plays significant roles in proper genome packing and its accessibility to execute transcription regulation. Despite a multitude of nucleosome positioning resources available on line including experimental datasets of genome-wide nucleosome occupancy profiles and computational tools to the analysis on these data, the complex language of eukaryotic Nucleosome positioning remains incompletely understood.

Results

Here, we address this challenge using an approach based on a state-of-the-art machine learning method. We present a novel convolutional neural network (CNN) to understand nucleosome positioning. We combined Inception-like networks with a gating mechanism for the response of multiple patterns and long term association in DNA sequences. We developed the open-source package LeNup based on the CNN to predict nucleosome positioning in Homo sapiens, Caenorhabditis elegans, Drosophila melanogaster as well as Saccharomyces cerevisiae genomes. We trained LeNup on four benchmark datasets. LeNup achieved greater predictive accuracy than previously published methods.

Availability and implementation

LeNup is freely available as Python and Lua script source code under a BSD style license from https://github.com/biomedBit/LeNup.

Supplementary information

Supplementary data are available at Bioinformatics online.

Low-Affinity Binding Sites and the Transcription Factor Specificity Paradox in Eukaryotes

Eukaryotic transcription factors (TFs) from the same structural family tend to bind similar DNA sequences, despite the ability of these TFs to execute distinct functions in vivo. The cell partly resolves this specificity paradox through combinatorial strategies and the use of low-affinity binding sites, which are better able to distinguish between similar TFs. However, because these sites have low affinity, it is challenging to understand how TFs recognize them in vivo. Here, we summarize recent findings and technological advancements that allow for the quantification and mechanistic interpretation of TF recognition across a wide range of affinities. We propose a model that integrates insights from the fields of genetics and cell biology to provide further conceptual understanding of TF binding specificity. We argue that in eukaryotes, target specificity is driven by an inhomogeneous 3D nuclear distribution of TFs and by variation in DNA binding affinity such that locally elevated TF concentration allows low-affinity binding sites to be functional.

ATAC-seq reveals regional differences in enhancer accessibility during the establishment of spatial coordinates in the Drosophila blastoderm

Establishment of spatial coordinates during Drosophila embryogenesis relies on differential regulatory activity of axis patterning enhancers. Concentration gradients of activator and repressor transcription factors (TFs) provide positional information to each enhancer, which in turn promotes transcription of a target gene in a specific spatial pattern. However, the interplay between an enhancer regulatory activity and its accessibility as determined by local chromatin organization is not well understood. We profiled chromatin accessibility with ATAC-seq in narrow, genetically tagged domains along the antero-posterior axis in the Drosophila blastoderm. We demonstrate that one-quarter of the accessible genome displays significant regional variation in its ATAC-seq signal immediately after zygotic genome activation. Axis patterning enhancers are enriched among the most variable intervals, and their accessibility changes correlate with their regulatory activity. In an embryonic domain where an enhancer receives a net activating TF input and promotes transcription, it displays elevated accessibility in comparison to a domain where it receives a net repressive input. We propose that differential accessibility is a signature of patterning cis-regulatory elements in the Drosophila blastoderm and discuss potential mechanisms by which accessibility of enhancers may be modulated by activator and repressor TFs.

Sumoylation of DNA-bound transcription factor Sko1 prevents its association with nontarget promoters

Sequence-specific transcription factors (TFs) represent one of the largest groups of proteins that is targeted for SUMO post-translational modification, in both yeast and humans. SUMO modification can have diverse effects, but recent studies showed that sumoylation reduces the interaction of multiple TFs with DNA in living cells. Whether this relates to a general role for sumoylation in TF binding site selection, however, has not been fully explored because few genome-wide studies aimed at studying such a role have been reported. To address this, we used genome-wide analysis to examine how sumoylation regulates Sko1, a yeast bZIP TF with hundreds of known binding sites. We find that Sko1 is sumoylated at Lys 567 and, although many of its targets are osmoresponse genes, the level of Sko1 sumoylation is not stress-regulated and the modification does not depend or impinge on its phosphorylation by the osmostress kinase Hog1. We show that Sko1 mutants that cannot bind DNA are not sumoylated, but attaching a heterologous DNA binding domain restores the modification, implicating DNA binding as a major determinant for Sko1 sumoylation. Genome-wide chromatin immunoprecipitation (ChIP-seq) analysis shows that a sumoylation-deficient Sko1 mutant displays increased occupancy levels at its numerous binding sites, which inhibits the recruitment of the Hog1 kinase to some induced osmostress genes. This strongly supports a general role for sumoylation in reducing the association of TFs with chromatin. Extending this result, remarkably, sumoylation-deficient Sko1 binds numerous additional promoters that are not normally regulated by Sko1 but contain sequences that resemble the Sko1 binding motif. Our study points to an important role for sumoylation in modulating the interaction of a DNA-bound TF with chromatin to increase the specificity of TF-DNA interactions.

Granule neuron precursor cell proliferation is regulated by NFIX and intersectin 1 during postnatal cerebellar development

Cerebellar granule neurons are the most numerous neuronal subtype in the central nervous system. Within the developing cerebellum, these neurons are derived from a population of progenitor cells found within the external granule layer of the cerebellar anlage, namely the cerebellar granule neuron precursors (GNPs). The timely proliferation and differentiation of these precursor cells, which, in rodents occurs predominantly in the postnatal period, is tightly controlled to ensure the normal morphogenesis of the cerebellum. Despite this, our understanding of the factors mediating how GNP differentiation is controlled remains limited. Here, we reveal that the transcription factor nuclear factor I X (NFIX) plays an important role in this process. Mice lacking Nfix exhibit reduced numbers of GNPs during early postnatal development, but elevated numbers of these cells at postnatal day 15. Moreover, Nfix^-/- GNPs exhibit increased proliferation when cultured in vitro, suggestive of a role for NFIX in promoting GNP differentiation. At a mechanistic level, profiling analyses using both ChIP-seq and RNA-seq identified the actin-associated factor intersectin 1 as a downstream target of NFIX during cerebellar development. In support of this, mice lacking intersectin 1 also displayed delayed GNP differentiation. Collectively, these findings highlight a key role for NFIX and intersectin 1 in the regulation of cerebellar development.

TIP55, a splice isoform of the KAT5 acetyltransferase, is essential for developmental gene regulation and organogenesis

Regulation of chromatin structure is critical for cell type-specific gene expression. Many chromatin regulatory complexes exist in several different forms, due to alternative splicing and differential incorporation of accessory subunits. However, in vivo studies often utilize mutations that eliminate multiple forms of complexes, preventing assessment of the specific roles of each. Here we examined the developmental roles of the TIP55 isoform of the KAT5 histone acetyltransferase. In contrast to the pre-implantation lethal phenotype of mice lacking all four Kat5 transcripts, mice specifically deficient for Tip55 die around embryonic day 11.5 (E11.5). Prior to developmental arrest, defects in heart and neural tube were evident in Tip55 mutant embryos. Specification of cardiac and neural cell fates appeared normal in Tip55 mutants. However, cell division and survival were impaired in heart and neural tube, respectively, revealing a role for TIP55 in cellular proliferation. Consistent with these findings, transcriptome profiling revealed perturbations in genes that function in multiple cell types and developmental pathways. These findings show that Tip55 is dispensable for the pre- and early post-implantation roles of Kat5, but is essential during organogenesis. Our results raise the possibility that isoform-specific functions of other chromatin regulatory proteins may play important roles in development.

Polycomb repressive complex 1 shapes the nucleosome landscape but not accessibility at target genes

Polycomb group (PcG) proteins are transcriptional repressors that play important roles in regulating gene expression during animal development. In vitro experiments have shown that PcG protein complexes can compact chromatin to limit the activity of chromatin remodeling enzymes and access of the transcriptional machinery to DNA. In fitting with these ideas, gene promoters associated with PcG proteins have been reported to be less accessible than other gene promoters. However, it remains largely untested in vivo whether PcG proteins define chromatin accessibility or other chromatin features. To address this important question, we examine the chromatin accessibility and nucleosome landscape at PcG protein-bound promoters in mouse embryonic stem cells using the assay for transposase accessible chromatin (ATAC)-seq. Combined with genetic ablation strategies, we unexpectedly discover that although PcG protein-occupied gene promoters exhibit reduced accessibility, this does not rely on PcG proteins. Instead, the Polycomb repressive complex 1 (PRC1) appears to play a unique role in driving elevated nucleosome occupancy and decreased nucleosomal spacing in Polycomb chromatin domains. Our new genome-scale observations argue, in contrast to the prevailing view, that PcG proteins do not significantly affect chromatin accessibility and highlight an underappreciated complexity in the relationship between chromatin accessibility, the nucleosome landscape, and PcG-mediated transcriptional repression.

Zbtb7a is a transducer for the control of promoter accessibility by NF-kappa B and multiple other transcription factors

Gene expression in eukaryotes is controlled by DNA sequences at promoter and enhancer regions, whose accessibility for binding by regulatory proteins dictates their specific patterns of activity. Here, we identify the protein Zbtb7a as a factor required for inducible changes in accessibility driven by transcription factors (TFs). We show that Zbtb7a binds to a significant fraction of genomic promoters and enhancers, encompassing many target genes of nuclear factor kappa B (NFκB) p65 and a variety of other TFs. While Zbtb7a binding is not alone sufficient to directly activate promoters, it is required to enable TF-dependent control of accessibility and normal gene expression. Using p65 as a model TF, we show that Zbtb7a associates with promoters independently of client TF binding. Moreover, the presence of prebound Zbtb7a can specify promoters that are amenable to TF-induced changes in accessibility. Therefore, Zbtb7a represents a widely used promoter factor that transduces signals from other TFs to enable control of accessibility and regulation of gene expression.

Gene activation is driven by the binding of regulatory proteins to the specific DNA sequences that control each gene. However, these sequences are not always accessible for binding in every type of cell, and so differences in their accessibility can underlie the range of cell types in which particular genes can be activated. Although several cellular processes can alter the accessibilities of these sequences, it is still often unclear how these processes are directed to act at specific genes. We have discovered that the protein Zbtb7a binds near numerous gene-regulatory sequences throughout the genome and that it enables other DNA-binding proteins to trigger changes in their accessibility and to activate nearby genes. However, unlike many other factors that control gene activation, the binding of Zbtb7a alone does not seem to be sufficient to switch on gene expression; instead, its function is required for activation of genes that are independently bound by a specific set of transcription factors (TFs), and it could therefore be considered to “transduce” their gene-regulatory activities. The implication of this is that the presence or absence of Zbtb7a at any gene in a particular cell type may represent one of the aspects that can determine whether that gene is able to be activated or not.

Genome-wide determinants of sequence-specific DNA binding of general regulatory factors

General regulatory factors (GRFs), such as Reb1, Abf1, Rap1, Mcm1, and Cbf1, positionally organize yeast chromatin through interactions with a core consensus DNA sequence. It is assumed that sequence recognition via direct base readout suffices for specificity and that spurious nonfunctional sites are rendered inaccessible by chromatin. We tested these assumptions through genome-wide mapping of GRFs in vivo and in purified biochemical systems at near–base pair (bp) resolution using several ChIP-exo–based assays. We find that computationally predicted DNA shape features (e.g., minor groove width, helix twist, base roll, and propeller twist) that are not defined by a unique consensus sequence are embedded in the nonunique portions of GRF motifs and contribute critically to sequence-specific binding. This dual source specificity occurs at GRF sites in promoter regions where chromatin organization starts. Outside of promoter regions, strong consensus sites lack the shape component and consequently lack an intrinsic ability to bind cognate GRFs, without regard to influences from chromatin. However, sites having a weak consensus and low intrinsic affinity do exist in these regions but are rendered inaccessible in a chromatin environment. Thus, GRF site-specificity is achieved through integration of favorable DNA sequence and shape readouts in promoter regions and by chromatin-based exclusion from fortuitous weak sites within gene bodies. This study further revealed a severe G/C nucleotide cross-linking selectivity inherent in all formaldehyde-based ChIP assays, which includes ChIP-seq. However, for most tested proteins, G/C selectivity did not appreciably affect binding site detection, although it does place limits on the quantitativeness of occupancy levels.

Resolving systematic errors in widely used enhancer activity assays in human cells

The identification of transcriptional enhancers in the human genome is a prime goal in biology. Enhancers are typically predicted via chromatin marks, yet their function is primarily assessed with plasmid-based reporter assays. Here, we show that such assays are rendered unreliable by two previously reported phenomena relating to plasmid transfection into human cells: (i) the bacterial plasmid origin of replication (ORI) functions as a conflicting core promoter and (ii) a type I interferon (IFN-I) response is activated. These cause confounding false positives and negatives in luciferase assays and STARR-seq screens. We overcome both problems by employing the ORI as core promoter and by inhibiting two IFN-I-inducing kinases, enabling genome-wide STARR-seq screens in human cells. In HeLa-S3 cells, we uncover strong enhancers, IFN-I-induced enhancers, and enhancers endogenously silenced at the chromatin level. Our findings apply to all episomal enhancer activity assays in mammalian cells and are key to the characterization of human enhancers.

Genomic copy number variation analysis in multiple system atrophy

Genomic variation includes single-nucleotide variants, small insertions or deletions (indels), and copy number variants (CNVs). CNVs affect gene expression by altering the genome structure and transposable elements within a region. CNVs are greater than 1 kb in size; hence, CNVs can produce more variation than can individual single-nucleotide variations that are detected by next-generation sequencing. Multiple system atrophy (MSA) is an α-synucleinopathy adult-onset disorder. Pathologically, it is characterized by insoluble aggregation of filamentous α-synuclein in brain oligodendrocytes. Generally, MSA is sporadic, although there are rare cases of familial MSA. In addition, the frequencies of the clinical phenotypes differ considerably among countries. Reports indicate that genetic factors play roles in the mechanisms involved in the pathology and onset of MSA. To evaluate the genetic background of this disorder, we attempted to determine whether there are differences in CNVs between patients with MSA and normal control subjects. We found that the number of CNVs on chromosomes 5, 22, and 4 was increased in MSA; 3 CNVs in non-coding regions were considered risk factors for MSA. Our results show that CNVs in non-coding regions influence the expression of genes through transcription-related mechanisms and potentially increase subsequent structural alterations of chromosomes. Therefore, these CNVs likely play roles in the molecular mechanisms underlying MSA.

Structural and biochemical analyses of monoubiquitinated human histones H2B and H4

Monoubiquitination is a major histone post-translational modification. In humans, the histone H2B K120 and histone H4 K31 residues are monoubiquitinated and may form transcriptionally active chromatin. In this study, we reconstituted nucleosomes containing H2B monoubiquitinated at position 120 (H2Bub₁₂₀) and/or H4 monoubiquitinated at position 31 (H4ub₃₁). We found that the H2Bub₁₂₀ and H4ub₃₁ monoubiquitinations differently affect nucleosome stability: the H2Bub₁₂₀ monoubiquitination enhances the H2A–H2B association with the nucleosome, while the H4ub₃₁ monoubiquitination decreases the H3–H4 stability in the nucleosome, when compared with the unmodified nucleosome. The H2Bub₁₂₀ and H4ub₃₁ monoubiquitinations both antagonize the Mg²⁺-dependent compaction of a poly-nucleosome, suggesting that these monoubiquitinations maintain more relaxed conformations of chromatin. In the crystal structure, the H2Bub₁₂₀ and H4ub₃₁ monoubiquitinations do not change the structure of the nucleosome core particle and the ubiquitin molecules were flexibly disordered in the H2Bub₁₂₀/H4ub₃₁ nucleosome structure. These results revealed the differences and similarities of the H2Bub₁₂₀ and H4ub₃₁ monoubiquitinations at the mono- and poly-nucleosome levels and provide novel information to clarify the roles of monoubiquitination in chromatin.

Identification of a core TP53 transcriptional program with highly distributed tumor suppressive activity

The tumor suppressor TP53 is the most frequently mutated gene product in human cancer. Close to half of all solid tumors carry inactivating mutations in the TP53 gene, while in the remaining cases, TP53 activity is abrogated by other oncogenic events, such as hyperactivation of its endogenous repressors MDM2 or MDM4. Despite identification of hundreds of genes regulated by this transcription factor, it remains unclear which direct target genes and downstream pathways are essential for the tumor suppressive function of TP53. We set out to address this problem by generating multiple genomic data sets for three different cancer cell lines, allowing the identification of distinct sets of TP53-regulated genes, from early transcriptional targets through to late targets controlled at the translational level. We found that although TP53 elicits vastly divergent signaling cascades across cell lines, it directly activates a core transcriptional program of ∼100 genes with diverse biological functions, regardless of cell type or cellular response to TP53 activation. This core program is associated with high-occupancy TP53 enhancers, high levels of paused RNA polymerases, and accessible chromatin. Interestingly, two different shRNA screens failed to identify a single TP53 target gene required for the anti-proliferative effects of TP53 during pharmacological activation in vitro. Furthermore, bioinformatics analysis of thousands of cancer genomes revealed that none of these core target genes are frequently inactivated in tumors expressing wild-type TP53. These results support the hypothesis that TP53 activates a genetically robust transcriptional program with highly distributed tumor suppressive functions acting in diverse cellular contexts.

Epigenetic Landscape during Coronavirus Infection

Coronaviruses (CoV) comprise a large group of emerging human and animal pathogens, including the highly pathogenic severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) strains. The molecular mechanisms regulating emerging coronavirus pathogenesis are complex and include virus–host interactions associated with entry, replication, egress and innate immune control. Epigenetics research investigates the genetic and non-genetic factors that regulate phenotypic variation, usually caused by external and environmental factors that alter host expression patterns and performance without any change in the underlying genotype. Epigenetic modifications, such as histone modifications, DNA methylation, chromatin remodeling, and non-coding RNAs, function as important regulators that remodel host chromatin, altering host expression patterns and networks in a highly flexible manner. For most of the past two and a half decades, research has focused on the molecular mechanisms by which RNA viruses antagonize the signaling and sensing components that regulate induction of the host innate immune and antiviral defense programs upon infection. More recently, a growing body of evidence supports the hypothesis that viruses, even lytic RNA viruses that replicate in the cytoplasm, have developed intricate, highly evolved, and well-coordinated processes that are designed to regulate the host epigenome, and control host innate immune antiviral defense processes, thereby promoting robust virus replication and pathogenesis. In this article, we discuss the strategies that are used to evaluate the mechanisms by which viruses regulate the host epigenome, especially focusing on highly pathogenic respiratory RNA virus infections as a model. By combining measures of epigenome reorganization with RNA and proteomic datasets, we articulate a spatial-temporal data integration approach to identify regulatory genomic clusters and regions that play a crucial role in the host’s innate immune response, thereby defining a new viral antagonism mechanism following emerging coronavirus infection.

NucTools: analysis of chromatin feature occupancy profiles from high-throughput sequencing data

Background

Biomedical applications of high-throughput sequencing methods generate a vast amount of data in which numerous chromatin features are mapped along the genome. The results are frequently analysed by creating binary data sets that link the presence/absence of a given feature to specific genomic loci. However, the nucleosome occupancy or chromatin accessibility landscape is essentially continuous. It is currently a challenge in the field to cope with continuous distributions of deep sequencing chromatin readouts and to integrate the different types of discrete chromatin features to reveal linkages between them.

Results

Here we introduce the NucTools suite of Perl scripts as well as MATLAB- and R-based visualization programs for a nucleosome-centred downstream analysis of deep sequencing data. NucTools accounts for the continuous distribution of nucleosome occupancy. It allows calculations of nucleosome occupancy profiles averaged over several replicates, comparisons of nucleosome occupancy landscapes between different experimental conditions, and the estimation of the changes of integral chromatin properties such as the nucleosome repeat length. Furthermore, NucTools facilitates the annotation of nucleosome occupancy with other chromatin features like binding of transcription factors or architectural proteins, and epigenetic marks like histone modifications or DNA methylation. The applications of NucTools are demonstrated for the comparison of several datasets for nucleosome occupancy in mouse embryonic stem cells (ESCs) and mouse embryonic fibroblasts (MEFs).

Conclusions

The typical workflows of data processing and integrative analysis with NucTools reveal information on the interplay of nucleosome positioning with other features such as for example binding of a transcription factor CTCF, regions with stable and unstable nucleosomes, and domains of large organized chromatin K9me2 modifications (LOCKs). As potential limitations and problems we discuss how inter-replicate variability of MNase-seq experiments can be addressed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3580-2) contains supplementary material, which is available to authorized users.

Systematic dissection of genomic features determining transcription factor binding and enhancer function

Enhancers regulate gene expression through the binding of sequence-specific transcription factors (TFs) to cognate motifs. Various features influence TF binding and enhancer function-including the chromatin state of the genomic locus, the affinities of the binding site, the activity of the bound TFs, and interactions among TFs. However, the precise nature and relative contributions of these features remain unclear. Here, we used massively parallel reporter assays (MPRAs) involving 32,115 natural and synthetic enhancers, together with high-throughput in vivo binding assays, to systematically dissect the contribution of each of these features to the binding and activity of genomic regulatory elements that contain motifs for PPARγ, a TF that serves as a key regulator of adipogenesis. We show that distinct sets of features govern PPARγ binding vs. enhancer activity. PPARγ binding is largely governed by the affinity of the specific motif site and higher-order features of the larger genomic locus, such as chromatin accessibility. In contrast, the enhancer activity of PPARγ binding sites depends on varying contributions from dozens of TFs in the immediate vicinity, including interactions between combinations of these TFs. Different pairs of motifs follow different interaction rules, including subadditive, additive, and superadditive interactions among specific classes of TFs, with both spatially constrained and flexible grammars. Our results provide a paradigm for the systematic characterization of the genomic features underlying regulatory elements, applicable to the design of synthetic regulatory elements or the interpretation of human genetic variation.

Noncommutative Biology: Sequential Regulation of Complex Networks

Single-cell variability in gene expression is important for generating distinct cell types, but it is unclear how cells use the same set of regulatory molecules to specifically control similarly regulated genes. While combinatorial binding of transcription factors at promoters has been proposed as a solution for cell-type specific gene expression, we found that such models resulted in substantial information bottlenecks. We sought to understand the consequences of adopting sequential logic wherein the time-ordering of factors informs the final outcome. We showed that with noncommutative control, it is possible to independently control targets that would otherwise be activated simultaneously using combinatorial logic. Consequently, sequential logic overcomes the information bottleneck inherent in complex networks. We derived scaling laws for two noncommutative models of regulation, motivated by phosphorylation/neural networks and chromosome folding, respectively, and showed that they scale super-exponentially in the number of regulators. We also showed that specificity in control is robust to the loss of a regulator. Lastly, we connected these theoretical results to real biological networks that demonstrate specificity in the context of promiscuity. These results show that achieving a desired outcome often necessitates roundabout steps.

DNA is the blueprint of life. Yet the order in which a cell follows these instructions makes it capable of generating thousands of different fates. How this information is extracted from underlying gene regulatory networks is unclear, especially given that biological networks are highly interconnected, and that the number of signaling pathways is relatively small (approximately 5–10). The conventional approach for increasing the information capacity of a limited set of regulators is to use them in combination. Surprisingly, combinatorial logic does not increase the diversity of target configurations or cell fates, but instead causes information bottlenecks. A different approach, called sequential logic, uses noncommutative sequences of a small set of regulators to drive networks to a large number of novel configurations. If certain targets are first protected, then even promiscuous regulators can activate specific subsets of lineage-specific targets. In this paper we show how sequential logic outperforms combinatorial logic, and argue that noncommutative sequences underlie a number of cases of biological regulation, e.g. how a small number of signaling pathways generates a large diversity of cell types in development. In addition to explaining biological networks, sequential logic may be a general experimental design strategy in synthetic and single-cell biology.

Expansion of GA Dinucleotide Repeats Increases the Density of CLAMP Binding Sites on the X-Chromosome to Promote Drosophila Dosage Compensation

Dosage compensation is an essential process that equalizes transcript levels of X-linked genes between sexes by forming a domain of coordinated gene expression. Throughout the evolution of Diptera, many different X-chromosomes acquired the ability to be dosage compensated. Once each newly evolved X-chromosome is targeted for dosage compensation in XY males, its active genes are upregulated two-fold to equalize gene expression with XX females. In Drosophila melanogaster, the CLAMP zinc finger protein links the dosage compensation complex to the X-chromosome. However, the mechanism for X-chromosome identification has remained unknown. Here, we combine biochemical, genomic and evolutionary approaches to reveal that expansion of GA-dinucleotide repeats likely accumulated on the X-chromosome over evolutionary time to increase the density of CLAMP binding sites, thereby driving the evolution of dosage compensation. Overall, we present new insight into how subtle changes in genomic architecture, such as expansions of a simple sequence repeat, promote the evolution of coordinated gene expression.

Regulation of chaperone binding and nucleosome dynamics by key residues within the globular domain of histone H3

BACKGROUND

Nucleosomes have an important role in modulating access of DNA by regulatory factors. The role specific histone residues have in this process has been shown to be an important mechanism of transcription regulation. Previously, we identified eight amino acids in histones H3 and H4 that are required for nucleosome occupancy over highly transcribed regions of the genome.

RESULTS

We investigate the mechanism through which three of these previously identified histone H3 amino acids regulate nucleosome architecture. We find that histone H3 K122, Q120, and R49 are required for Spt2, Spt6, and Spt16 occupancies at genomic locations where transcription rates are high, but not over regions of low transcription rates. Furthermore, substitution at one residue, K122, located on the dyad axis of the nucleosome, results in improper reassembly and disassembly of nucleosomes, likely accounting for the transcription rate-dependent regulation by these mutant histones.

CONCLUSIONS

These data show that when specific amino acids of histone proteins are substituted, Spt2, Spt6, and Spt16 occupancies are reduced and nucleosome dynamics are altered. Therefore, these data support a mechanism for histone chaperone binding where these factors interact with histone proteins to promote their activities during transcription.

A SAM oligomerization domain shapes the genomic binding landscape of the LEAFY transcription factor

Deciphering the mechanisms directing transcription factors (TFs) to specific genome regions is essential to understand and predict transcriptional regulation. TFs recognize short DNA motifs primarily through their DNA-binding domain. Some TFs also possess an oligomerization domain suspected to potentiate DNA binding but for which the genome-wide influence remains poorly understood. Here we focus on the LEAFY transcription factor, a master regulator of flower development in angiosperms. We have determined the crystal structure of its conserved amino-terminal domain, revealing an unanticipated Sterile Alpha Motif oligomerization domain. We show that this domain is essential to LEAFY floral function. Moreover, combined biochemical and genome-wide assays suggest that oligomerization is required for LEAFY to access regions with low-affinity binding sites or closed chromatin. This finding shows that domains that do not directly contact DNA can nevertheless have a profound impact on the DNA binding landscape of a TF.

The LEAFY transcription factor is a master regulator of flower development in plants. Here the authors describe the structure of a LEAFY oligomerization domain and show that mutations that disrupt oligomerization alter its capacity to bind low affinity and poorly accessible target sites.

The Effects of Sequence Variation on Genome-wide NRF2 Binding--New Target Genes and Regulatory SNPs

Transcription factor binding specificity is crucial for proper target gene regulation. Motif discovery algorithms identify the main features of the binding patterns, but the accuracy on the lower affinity sites is often poor. Nuclear factor E2-related factor 2 (NRF2) is a ubiquitous redox-activated transcription factor having a key protective role against endogenous and exogenous oxidant and electrophile stress. Herein, we decipher the effects of sequence variation on the DNA binding sequence of NRF2, in order to identify both genome-wide binding sites for NRF2 and disease-associated regulatory SNPs (rSNPs) with drastic effects on NRF2 binding. Interactions between NRF2 and DNA were studied using molecular modelling, and NRF2 chromatin immunoprecipitation-sequence datasets together with protein binding microarray measurements were utilized to study binding sequence variation in detail. The binding model thus generated was used to identify genome-wide binding sites for NRF2, and genomic binding sites with rSNPs that have strong effects on NRF2 binding and reside on active regulatory elements in human cells. As a proof of concept, miR-126–3p and -5p were identified as NRF2 target microRNAs, and a rSNP (rs113067944) residing on NRF2 target gene (Ferritin, light polypeptide, FTL) promoter was experimentally verified to decrease NRF2 binding and result in decreased transcriptional activity.

Nucleosome positioning: resources and tools online

Nucleosome positioning is an important process required for proper genome packing and its accessibility to execute the genetic program in a cell-specific, timely manner. In the recent years hundreds of papers have been devoted to the bioinformatics, physics and biology of nucleosome positioning. The purpose of this review is to cover a practical aspect of this field, namely, to provide a guide to the multitude of nucleosome positioning resources available online. These include almost 300 experimental datasets of genome-wide nucleosome occupancy profiles determined in different cell types and more than 40 computational tools for the analysis of experimental nucleosome positioning data and prediction of intrinsic nucleosome formation probabilities from the DNA sequence. A manually curated, up to date list of these resources will be maintained at http://generegulation.info.

The Ku subunit of telomerase binds Sir4 to recruit telomerase to lengthen telomeres in S. cerevisiae

In Saccharomyces cerevisiae and in humans, the telomerase RNA subunit is bound by Ku, a ring-shaped protein heterodimer best known for its function in DNA repair. Ku binding to yeast telomerase RNA promotes telomere lengthening and telomerase recruitment to telomeres, but how this is achieved remains unknown. Using telomere-length analysis and chromatin immunoprecipitation, we show that Sir4 – a previously identified Ku-binding protein that is a component of telomeric silent chromatin – is required for Ku-mediated telomere lengthening and telomerase recruitment. We also find that specifically tethering Sir4 directly to Ku-binding-defective telomerase RNA restores otherwise-shortened telomeres to wild-type length. These findings suggest that Sir4 is the telomere-bound target of Ku-mediated telomerase recruitment and provide one mechanism for how the Sir4-competing Rif1 and Rif2 proteins negatively regulate telomere length in yeast.

DOI:http://dx.doi.org/10.7554/eLife.07750.001

Inside a cell's nucleus, DNA is packaged into structures called chromosomes. The ends of every chromosome are capped by repeating sequences of DNA known as telomeres, which protect the chromosomes from damage. Every time a cell divides, the telomeres shorten. If telomere length falls below a critical level, the cell can die or enter a state in which it can no longer divide.

During cell division, an enzyme called telomerase normally restores telomeres to their original length. Telomerase is made up of several proteins and an RNA molecule. In yeast and humans, a protein called Ku is one part of the telomerase enzyme. Ku binds to the RNA subunit of telomerase and helps the enzyme find and interact with the telomeres. Previous research has shown that Ku is unable to work alone to recruit telomerase to the chromosome.

A protein called Sir4 binds to telomeres and cells lacking it have short telomeres, but the reason behind this was not known. Hass and Zappulla confirmed previous reports that Ku binds to Sir4 using a biochemical approach. Additional experiments provided genetic evidence that this binding interaction is important for telomerase to lengthen telomeres appropriately.

Cells in which the RNA subunit of telomerase is unable to bind effectively to Ku have short telomeres. Hass and Zappulla directly tethered Sir4 to this defective RNA and found this restored the shortened telomeres to a normal length, indicating that Sir4 normally binds Ku to recruit telomerase. Discovering this mode of recruitment also helps to explain how two other telomeric proteins (Rif1 and 2) limit telomere lengthening; they compete with Ku-Sir4 recruitment to form a length-regulating system.

Taken together, Hass and Zappulla's results provide strong evidence that Sir4 cooperates with Ku to control the lengthening of chromosome ends. Future research will hopefully reveal the precise space and time requirements for this telomerase-controlling system in yeast. Additionally, because Ku has been reported to be a subunit of human telomerase, future studies could also explore whether human cells use a similar strategy.

DOI:http://dx.doi.org/10.7554/eLife.07750.002

Unraveling the 3D genome: genomics tools for multiscale exploration

A decade of rapid method development has begun to yield exciting insights into the 3D architecture of the metazoan genome and the roles it may play in regulating transcription. Here we review core methods and new tools in the modern genomicist's toolbox at three length scales, ranging from single base pairs to megabase-scale chromosomal domains, and discuss the emerging picture of the 3D genome that these tools have revealed. Blind spots remain, especially at intermediate length scales spanning a few nucleosomes, but thanks in part to new technologies that permit targeted alteration of chromatin states and time-resolved studies, the next decade holds great promise for hypothesis-driven research into the mechanisms that drive genome architecture and transcriptional regulation.

Activity of heat shock genes' promoters in thermally contrasting animal species

Heat shock gene promoters represent a highly conserved and universal system for the rapid induction of transcription after various stressful stimuli. We chose pairs of mammalian and insect species that significantly differ in their thermoresistance and constitutive levels of Hsp70 to compare hsp promoter strength under normal conditions and after heat shock (HS). The first pair includes the HSPA1 gene promoter of camel (Camelus dromedarius) and humans. It was demonstrated that the camel HSPA1A and HSPA1L promoters function normally in vitro in human cell cultures and exceed the strength of orthologous human promoters under basal conditions. We used the same in vitro assay for Drosophila melanogaster Schneider-2 (S2) cells to compare the activity of the hsp70 and hsp83 promoters of the second species pair represented by Diptera, i.e., Stratiomys singularior and D. melanogaster, which dramatically differ in thermoresistance and the pattern of Hsp70 accumulation. Promoter strength was also monitored in vivo in D. melanogaster strains transformed with constructs containing the S. singularior hsp70 ORF driven either by its own promoter or an orthologous promoter from the D. melanogaster hsp70Aa gene. Analysis revealed low S. singularior hsp70 promoter activity in vitro and in vivo under basal conditions and after HS in comparison with the endogenous promoter in D. melanogaster cells, which correlates with the absence of canonical GAGA elements in the promoters of the former species. Indeed, the insertion of GAGA elements into the S. singularior hsp70 regulatory region resulted in a dramatic increase in promoter activity in vitro but only modestly enhanced the promoter strength in the larvae of the transformed strains. In contrast with hsp70 promoters, hsp83 promoters from both of the studied Diptera species demonstrated high conservation and universality.

Interactions of chromatin context, binding site sequence content, and sequence evolution in stress-induced p53 occupancy and transactivation

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10⁻⁷) and correlated with stronger p53RE sequences (p<10⁻¹¹⁰) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.

It is well established that p53 binds DNA elements near p53 target genes to regulate the response to cellular stress. To assess factors influencing binding to response elements and subsequent gene expression, we have analyzed 2932 p53-occupied response elements (p53REs) in the context of genome-wide chromatin state, DNA accessibility and dynamics, and considered roles for binding-sequence specificity and evolutionary conservation. While p53 occupancy level shows little apparent direct relationship to gene expression change, after grouping expressed genes by their chromatin status at baseline, a relationship between occupancy of p53REs and gene expression change emerged. Analysis of p53RE sequences demonstrated that p53 occupancy was strongly correlated with sequence similarity to p53RE consensus, but surprisingly, was inversely correlated with nucleosome accessibility (DHS) and evolutionary conservation. These data revealed a systematic interaction between p53RE content and chromatin context that affects both quantitative p53 occupancy and the induced transactivation response to exposure. Moreover, this interaction appears to have been tuned via evolutionary events involving transposable elements, which strongly bind p53, but in only a few instances affect gene expression levels. Models of p53-regulated gene expression response that consider both chromatin state and sequence context may prove useful in guiding strategies for cancer prevention or therapy.

p53 binding to human genome: crowd control navigation in chromatin context

p53 is the most studied human protein because of its role in maintaining genomic stability. Binding to genomic targets is essential for transcription-dependent p53 tumor suppression, but how p53 selects targets remains unclear. Here, the impact of chromatin context on p53 genome-wide binding and targets selection is discussed. It is proposed that p53 genomic binding serves not only to regulate transcription, but to sense epigenomic changes threatening the genomic integrity. The problem of p53 navigating the human genome is discussed with respect to the degenerate p53 binding motif. This discussion relates to the fundamental problem of DNA binding factors navigating large genomes in search for cognate binding sites.

Single molecule fluorescence methodologies for investigating transcription factor binding kinetics to nucleosomes and DNA

Site specific DNA binding complexes must bind their DNA target sites and then reside there for a sufficient amount of time for proper regulation of DNA processing including transcription, replication and DNA repair. In eukaryotes, the occupancy of DNA binding complexes at their target sites is regulated by chromatin structure and dynamics. Methodologies that probe both the binding and dissociation kinetics of DNA binding proteins with naked and nucleosomal DNA are essential for understanding the mechanisms by which these complexes function. Here, we describe single-molecule fluorescence methodologies for quantifying the binding and dissociation kinetics of transcription factors at a target site within DNA, nucleosomes and nucleosome arrays. This approach allowed for the unexpected observation that nucleosomes impact not only binding but also dissociation kinetics of transcription factors and is well-suited for the investigation of numerous DNA processing complexes that directly interact with DNA organized into chromatin.

DNase footprint signatures are dictated by factor dynamics and DNA sequence

Genomic footprinting has emerged as an unbiased discovery method for transcription factor (TF) occupancy at cognate DNA in vivo. A basic premise of footprinting is that sequence-specific TF-DNA interactions are associated with localized resistance to nucleases, leaving observable signatures of cleavage within accessible chromatin. This phenomenon is interpreted to imply protection of the critical nucleotides by the stably bound protein factor. However, this model conflicts with previous reports of many TFs exchanging with specific binding sites in living cells on a timescale of seconds. We show that TFs with short DNA residence times have no footprints at bound motif elements. Moreover, the nuclease cleavage profile within a footprint originates from the DNA sequence in the factor-binding site, rather than from the protein occupying specific nucleotides. These findings suggest a revised understanding of TF footprinting and reveal limitations in comprehensive reconstruction of the TF regulatory network using this approach.

Targeted H3R26 deimination specifically facilitates estrogen receptor binding by modifying nucleosome structure

Transcription factor binding to DNA in vivo causes the recruitment of chromatin modifiers that can cause changes in chromatin structure, including the modification of histone tails. We previously showed that estrogen receptor (ER) target gene activation is facilitated by peptidylarginine deiminase 2 (PAD2)-catalyzed histone H3R26 deimination (H3R26Cit). Here we report that the genomic distributions of ER and H3R26Cit in breast cancer cells are strikingly coincident, linearly correlated, and observed as early as 2 minutes following estradiol treatment. The H3R26Cit profile is unlike that of previously described histone modifications and is characterized by sharp, narrow peaks. Paired-end MNase ChIP-seq indicates that the charge-neutral H3R26Cit modification facilitates ER binding to DNA by altering the fine structure of the nucleosome. Clinically, we find that PAD2 and H3R26Cit levels correlate with ER expression in breast tumors and that high PAD2 expression is associated with increased survival in ER+ breast cancer patients. These findings provide insight into how transcription factors gain access to nucleosomal DNA and implicate PAD2 as a novel therapeutic target for ER+ breast cancer.

Transcription factors bind to DNA to activate and repress gene transcription. Many transcription factors, particularly nuclear receptors, associate with their cognate DNA element in a highly dynamic manner in vivo. Highly acetylated histone tails and DNase sensitive chromatin are amenable to the initial binding of transcription factors. Upon binding to DNA, transcription factor binding recruits remodelers and coactivators that can cause a concomitant increase in accessibility and acetylation. Herein, we show that estrogen receptor recruitment of a histone deiminase causes the positively charged H3R26 residue to be neutralized. This modification changes the fine structure of the nucleosome particle and facilitates estrogen receptor binding. Lastly, we find that high deiminase expression is associated with increased survival in estrogen receptor-positive breast cancer patients.

Transient estrogen receptor binding and p300 redistribution support a squelching mechanism for estradiol-repressed genes

Proper gene regulation is essential for proper organismal development and appropriate responses to external stimuli. Specialized factors, termed master regulators, are often responsible for orchestrating the molecular events that result from signaling cascades. Master regulators coordinate the activation and repression of specific gene classes. Estrogen receptor α (ER) precipitates the signaling cascade that results from endogenous or exogenous estrogen hormones. ER is a classic transcriptional activator and the mechanisms by which ER coordinates gene activation are well characterized. However, it remains unclear how ER coordinates the immediate repression of genes. We integrated genomic transcription, chromosome looping, transcription factor binding, and chromatin structure data to analyze the molecular cascade that results from estradiol (E2)-induced signaling in human MCF-7 breast cancer cells and addressed the context-specific nature of gene regulation. We defined a class of genes that are immediately repressed upon estrogen stimulation, and we compared and contrasted the molecular characteristics of these repressed genes vs activated and unregulated genes. The most striking and unique feature of the repressed gene class is transient binding of ER at early time points after estrogen stimulation. We also found that p300, a coactivator and acetyltransferase, quantitatively redistributes from non-ER enhancers to ER enhancers after E2 treatment. These data support an extension of the classic physiological squelching model, whereby ER hijacks coactivators from repressed genes and redistributes the coactivators to ER enhancers that activate transcription.