Genome-wide protein-DNA binding dynamics suggest a molecular clutch for transcription factor function

Transcriptional silencing in Saccharomyces cerevisiae: known unknowns

Transcriptional silencing in Saccharomyces cerevisiae is a persistent and highly stable form of gene repression. It involves DNA silencers and repressor proteins that bind nucleosomes. The silenced state is influenced by numerous factors including the concentration of repressors, nature of activators, architecture of regulatory elements, modifying enzymes and the dynamics of chromatin.Silencers function to increase the residence time of repressor Sir proteins at silenced domains while clustering of silenced domains enables increased concentrations of repressors and helps facilitate long-range interactions. The presence of an accessible NDR at the regulatory regions of silenced genes, the cycling of chromatin configurations at regulatory sites, the mobility of Sir proteins, and the non-uniform distribution of the Sir proteins across the silenced domain, all result in silenced chromatin that only stably silences weak promoters and enhancers via changes in transcription burst duration and frequency.These data collectively suggest that silencing is probabilistic and the robustness of silencing is achieved through sub-optimization of many different nodes of action such that a stable expression state is generated and maintained even though individual constituents are in constant flux.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

Dual-role transcription factors stabilize intermediate expression levels

Precise control of gene expression levels is essential for normal cell functions, yet how they are defined and tightly maintained, particularly at intermediate levels, remains elusive. Here, using a series of newly developed sequencing, imaging, and functional assays, we uncover a class of transcription factors with dual roles as activators and repressors, referred to as condensate-forming level-regulating dual-action transcription factors (TFs). They reduce high expression but increase low expression to achieve stable intermediate levels. Dual-action TFs directly exert activating and repressing functions via condensate-forming domains that compartmentalize core transcriptional unit selectively. Clinically relevant mutations in these domains, which are linked to a range of developmental disorders, impair condensate selectivity and dual-action TF activity. These results collectively address a fundamental question in expression regulation and demonstrate the potential of level-regulating dual-action TFs as powerful effectors for engineering controlled expression levels.

Genome-scale chromatin binding dynamics of RNA Polymerase II general transcription machinery components

A great deal of work has revealed, in structural detail, the components of the preinitiation complex (PIC) machinery required for initiation of mRNA gene transcription by RNA polymerase II (Pol II). However, less-well understood are the in vivo PIC assembly pathways and their kinetics, an understanding of which is vital for determining how rates of in vivo RNA synthesis are established. We used competition ChIP in budding yeast to obtain genome-scale estimates of the residence times for five general transcription factors (GTFs): TBP, TFIIA, TFIIB, TFIIE and TFIIF. While many GTF-chromatin interactions were short-lived ( < 1 min), there were numerous interactions with residence times in the range of several minutes. Sets of genes with a shared function also shared similar patterns of GTF kinetic behavior. TFIIE, a GTF that enters the PIC late in the assembly process, had residence times correlated with RNA synthesis rates. The datasets and results reported here provide kinetic information for most of the Pol II-driven genes in this organism, offering a rich resource for exploring the mechanistic relationships between PIC assembly, gene regulation, and transcription.

PML::RARA and GATA2 proteins interact via DNA templates to induce aberrant self-renewal in mouse and human hematopoietic cells

The underlying mechanism(s) by which the PML::RARA fusion protein initiates acute promyelocytic leukemia is not yet clear. We defined the genomic binding sites of PML::RARA in primary mouse and human hematopoietic progenitor cells with V5-tagged PML::RARA, using anti-V5-PML::RARA chromatin immunoprecipitation sequencing and CUT&RUN approaches. Most genomic PML::RARA binding sites were found in regions that were already chromatin-accessible (defined by ATAC-seq) in unmanipulated, wild-type promyelocytes, suggesting that these regions are "open" prior to PML::RARA expression. We found that GATA binding motifs, and the direct binding of the chromatin "pioneering factor" GATA2, were significantly enriched near PML::RARA binding sites. Proximity labeling studies revealed that PML::RARA interacts with ~250 proteins in primary mouse hematopoietic cells; GATA2 and 33 others require PML::RARA binding to DNA for the interaction to occur, suggesting that binding to their cognate DNA target motifs may stabilize their interactions. In the absence of PML::RARA, Gata2 overexpression induces many of the same epigenetic and transcriptional changes as PML::RARA. These findings suggested that PML::RARA may indirectly initiate its transcriptional program by activating Gata2 expression: Indeed, we demonstrated that inactivation of Gata2 prior to PML::RARA expression prevented its ability to induce self-renewal. These data suggested that GATA2 binding creates accessible chromatin regions enriched for both GATA and Retinoic Acid Receptor Element motifs, where GATA2 and PML::RARA can potentially bind and interact with each other. In turn, PML::RARA binding to DNA promotes a feed-forward transcriptional program by positively regulating Gata2 expression. Gata2 may therefore be required for PML::RARA to establish its transcriptional program.

Yeast heterochromatin stably silences only weak regulatory elements by altering burst duration

Transcriptional silencing in Saccharomyces cerevisiae involves the generation of a chromatin state that stably represses transcription. Using multiple reporter assays, a diverse set of upstream activating sequence enhancers and core promoters were investigated for their susceptibility to silencing. We show that heterochromatin stably silences only weak and stress-induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements, and the partial repression of these elements did not result in bistable expression states. Permutation analysis of enhancers and promoters indicates that both elements are targets of repression. Chromatin remodelers help specific regulatory elements to resist repression, most probably by altering nucleosome mobility and changing transcription burst duration. The strong enhancers/promoters can be repressed if silencer-bound Sir1 is increased. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating-type gene regulatory elements but not strong housekeeping gene regulatory sequences.

Transcription factor exchange enables prolonged transcriptional bursts

Single-molecule imaging inside living cells has revealed that transcription factors (TFs) bind to DNA transiently, but a long-standing question is how this transient binding is related to transcription activation. Here, we devised a microscopy method to simultaneously measure transient TF binding at a single locus and the effect of these binding events on transcription. We show that DNA binding of the yeast TF Gal4 activates transcription of a target gene within a few seconds, with at least ∼20% efficiency and with a high initiation rate of ∼1 RNA/s. Gal4 DNA dissociation decreases transcription rapidly. Moreover, at a gene with multiple binding sites, individual Gal4 molecules only rarely stay bound throughout the entire burst but instead frequently exchange during a burst to increase the transcriptional burst duration. Our results suggest a mechanism for enhancer regulation in more complex eukaryotes, where TF cooperativity and exchange enable robust and responsive transcription regulation.

Characterizing Different Modes of Interplay Between Rap1 and H3 Using Inducible H3-depletion Yeast

Histones and transcription factors (TFs) are two important DNA-binding proteins that interact, compete, and together regulate transcriptional processes in response to diverse internal and external stimuli. Condition-specific depletion of histones in Saccharomyces cerevisiae using a galactose-inducible H3 promoter provides a suitable framework for examining transcriptional alteration resulting from reduced nucleosome content. However, the effect on DNA binding activities of TFs is yet to be fully explored. In this work, we combine ChIP-seq of H3 with RNA-seq to elucidate the genome-scale relationships between H3 occupancy patterns and transcriptional dynamics before and after global H3 depletion. ChIP-seq of Rap1 is also conducted in the H3-depletion and control treatments, to investigate the interplay between this master regulator TF and nucleosomal H3, and to explore the impact on diverse transcriptional responses of different groups of target genes and functions. Ultimately, we propose a working model and testable hypotheses regarding the impact of global and local H3 depletion on transcriptional changes. We also demonstrate different potential modes of interaction between Rap1 and H3, which sheds light on the potential multifunctional regulatory capabilities of Rap1 and potentially other pioneer factors.

Active compensation for changes in TDH3 expression mediated by direct regulators of TDH3 in Saccharomyces cerevisiae

Genetic networks are surprisingly robust to perturbations caused by new mutations. This robustness is conferred in part by compensation for loss of a gene's activity by genes with overlapping functions, such as paralogs. Compensation occurs passively when the normal activity of one paralog can compensate for the loss of the other, or actively when a change in one paralog's expression, localization, or activity is required to compensate for loss of the other. The mechanisms of active compensation remain poorly understood in most cases. Here we investigate active compensation for the loss or reduction in expression of the Saccharomyces cerevisiae gene TDH3 by its paralog TDH2. TDH2 is upregulated in a dose-dependent manner in response to reductions in TDH3 by a mechanism requiring the shared transcriptional regulators Gcr1p and Rap1p. TDH1, a second and more distantly related paralog of TDH3, has diverged in its regulation and is upregulated by another mechanism. Other glycolytic genes regulated by Rap1p and Gcr1p show changes in expression similar to TDH2, suggesting that the active compensation by TDH3 paralogs is part of a broader homeostatic response mediated by shared transcriptional regulators.

Yeast Heterochromatin Only Stably Silences Weak Regulatory Elements by Altering Burst Duration

The interplay between nucleosomes and transcription factors leads to programs of gene expression. Transcriptional silencing involves the generation of a chromatin state that represses transcription and is faithfully propagated through DNA replication and cell division. Using multiple reporter assays, including directly visualizing transcription in single cells, we investigated a diverse set of UAS enhancers and core promoters for their susceptibility to heterochromatic gene silencing. These results show that heterochromatin only stably silences weak and stress induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements and the partial repression did not result in bistable expression states. Permutation analysis of different UAS enhancers and core promoters indicate that both elements function together to determine the susceptibility of regulatory sequences to repression. Specific histone modifiers and chromatin remodellers function in an enhancer specific manner to aid these elements to resist repression suggesting that Sir proteins likely function in part by reducing nucleosome mobility. We also show that the strong housekeeping regulatory elements can be repressed if silencer bound Sir1 is increased, suggesting that Sir1 is a limiting component in silencing. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating type gene regulatory elements but not strong housekeeping gene regulatory sequences which could help explain why these genes are often found at the boundaries of silenced domains.

Context-dependent function of the transcriptional regulator Rap1 in gene silencing and activation in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, heterochromatin is formed through interactions between site-specific DNA-binding factors, including the transcriptional activator Repressor Activator Protein (Rap1), and Sir proteins. Despite an understanding of the establishment and maintenance of Sir-silenced chromatin, the mechanism of gene silencing by Sir proteins has remained a mystery. Utilizing high-resolution chromatin immunoprecipitation, we found that Rap1, the native activator of the bidirectional HMLα promoter, bound its recognition sequence in silenced chromatin, and its binding was enhanced by the presence of Sir proteins. In contrast to prior results, various components of transcription machinery were not able to access HMLα in the silenced state. These findings disproved the long-standing model of indiscriminate steric occlusion by Sir proteins and led to investigation of the role of the transcriptional activator Rap1 in Sir-silenced chromatin. Using a highly sensitive assay that monitors loss-of-silencing events, we identified a role for promoter-bound Rap1 in the maintenance of silent chromatin through interactions with the Sir complex. We also found that promoter-bound Rap1 activated HMLα when in an expressed state, and aided in the transition from transcription initiation to elongation. Highlighting the importance of epigenetic context in transcription factor function, these results point toward a model in which the duality of Rap1 function was mediated by local chromatin environment rather than binding-site availability.

Acetylation reprograms MITF target selectivity and residence time

The ability of transcription factors to discriminate between different classes of binding sites associated with specific biological functions underpins effective gene regulation in development and homeostasis. How this is achieved is poorly understood. The microphthalmia-associated transcription factor MITF is a lineage-survival oncogene that plays a crucial role in melanocyte development and melanoma. MITF suppresses invasion, reprograms metabolism and promotes both proliferation and differentiation. How MITF distinguishes between differentiation and proliferation-associated targets is unknown. Here we show that compared to many transcription factors MITF exhibits a very long residence time which is reduced by p300/CBP-mediated MITF acetylation at K206. While K206 acetylation also decreases genome-wide MITF DNA-binding affinity, it preferentially directs DNA binding away from differentiation-associated CATGTG motifs toward CACGTG elements. The results reveal an acetylation-mediated switch that suppresses differentiation and provides a mechanistic explanation of why a human K206Q MITF mutation is associated with Waardenburg syndrome.

A spatio-temporally constrained gene regulatory network directed by PBX1/2 acquires limb patterning specificity via HAND2

A lingering question in developmental biology has centered on how transcription factors with widespread distribution in vertebrate embryos can perform tissue-specific functions. Here, using the murine hindlimb as a model, we investigate the elusive mechanisms whereby PBX TALE homeoproteins, viewed primarily as HOX cofactors, attain context-specific developmental roles despite ubiquitous presence in the embryo. We first demonstrate that mesenchymal-specific loss of PBX1/2 or the transcriptional regulator HAND2 generates similar limb phenotypes. By combining tissue-specific and temporally controlled mutagenesis with multi-omics approaches, we reconstruct a gene regulatory network (GRN) at organismal-level resolution that is collaboratively directed by PBX1/2 and HAND2 interactions in subsets of posterior hindlimb mesenchymal cells. Genome-wide profiling of PBX1 binding across multiple embryonic tissues further reveals that HAND2 interacts with subsets of PBX-bound regions to regulate limb-specific GRNs. Our research elucidates fundamental principles by which promiscuous transcription factors cooperate with cofactors that display domain-restricted localization to instruct tissue-specific developmental programs.

Context dependent function of the transcriptional regulator Rap1 in gene silencing and activation in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, heterochromatin is formed through interactions between site-specific DNA-binding factors, including the transcriptional activator Rap1, and Sir proteins. Despite a vast understanding of the establishment and maintenance of Sir-silenced chromatin, the mechanism of gene silencing by Sir proteins has remained a mystery. Utilizing high resolution chromatin immunoprecipitation, we found that Rap1, the native activator of the bi-directional HML α promoter, bound its recognition sequence in silenced chromatin and its binding was enhanced by the presence of Sir proteins. In contrast to prior results, various components of transcription machinery were not able to access HML α in the silenced state. These findings disproved the long-standing model of indiscriminate steric occlusion by Sir proteins and led to investigation of the transcriptional activator Rap1 in Sir-silenced chromatin. Using a highly sensitive assay that monitors loss-of-silencing events, we identified a novel role for promoter-bound Rap1 in the maintenance of silent chromatin through interactions with the Sir complex. We also found that promoter-bound Rap1 activated HML α when in an expressed state, and aided in the transition from transcription initiation to elongation. Highlighting the importance of epigenetic context in transcription factor function, these results point toward a model in which the duality of Rap1 function was mediated by local chromatin environment rather than binding-site availability.

Significance Statement

The coarse partitioning of the genome into regions of active euchromatin and repressed heterochromatin is an important, and conserved, level gene expression regulation in eukaryotes. Repressor Activator Protein (Rap1) is a transcription factor that promotes the activation of genes when recruited to promoters, and aids in the establishment of heterochromatin through interactions with silencer elements. Here, we investigate the role of Rap1 when bound to a promoter in silent chromatin and dissect the context-specific epigenetic cues that regulate the dual properties of this transcription factor. Together, our data highlight the importance of protein-protein interactions and local chromatin state on transcription factor function.

Transcription factor localization dynamics and DNA binding drive distinct promoter interpretations

Environmental information may be encoded in the temporal dynamics of transcription factor (TF) activation and subsequently decoded by gene promoters to enact stimulus-specific gene expression programs. Previous studies of this behavior focused on the encoding and decoding of information in TF nuclear localization dynamics, yet cells control the activity of TFs in myriad ways, including by regulating their ability to bind DNA. Here, we use light-controlled mutants of the yeast TF Msn2 as a model system to investigate how promoter decoding of TF localization dynamics is affected by changes in the ability of the TF to bind DNA. We find that yeast promoters directly decode the light-controlled localization dynamics of Msn2 and that the effects of changing Msn2 affinity on that decoding behavior are highly promoter dependent, illustrating how cells could regulate TF localization dynamics and DNA binding in concert for improved control of gene expression.

Dynamics of nuclear architecture during early embryonic development and lessons from liveimaging

Nuclear organization has emerged as a potential key regulator of genome function. During development, the deployment of transcriptional programs must be tightly coordinated with cell division and is often accompanied by major changes in the repertoire of expressed genes. These transcriptional and developmental events are paralleled by changes in the chromatin landscape. Numerous studies have revealed the dynamics of nuclear organization underlying them. In addition, advances in live-imaging-based methodologies enable the study of nuclear organization with high spatial and temporal resolution. In this Review, we summarize the current knowledge of the changes in nuclear architecture in the early embryogenesis of various model systems. Furthermore, to highlight the importance of integrating fixed-cell and live approaches, we discuss how different live-imaging techniques can be applied to examine nuclear processes and their contribution to our understanding of transcription and chromatin dynamics in early development. Finally, we provide future avenues for outstanding questions in this field.

Active compensation for changes in TDH3 expression mediated by direct regulators of TDH3 in Saccharomyces cerevisiae

Genetic networks are surprisingly robust to perturbations caused by new mutations. This robustness is conferred in part by compensation for loss of a gene's activity by genes with overlapping functions, such as paralogs. Compensation occurs passively when the normal activity of one paralog can compensate for the loss of the other, or actively when a change in one paralog's expression, localization, or activity is required to compensate for loss of the other. The mechanisms of active compensation remain poorly understood in most cases. Here we investigate active compensation for the loss or reduction in expression of the Saccharomyces cerevisiae gene TDH3 by its paralogs TDH1 and TDH2. TDH1 and TDH2 are upregulated in a dose-dependent manner in response to reductions in TDH3 by a mechanism requiring the shared transcriptional regulators Gcr1p and Rap1p. Other glycolytic genes regulated by Rap1p and Gcr1p show changes in expression similar to TDH2, suggesting that the active compensation by TDH3 paralogs is part of a broader homeostatic response mediated by shared transcriptional regulators.

TF-Prioritizer: a Java pipeline to prioritize condition-specific transcription factors

BACKGROUND

Eukaryotic gene expression is controlled by cis-regulatory elements (CREs), including promoters and enhancers, which are bound by transcription factors (TFs). Differential expression of TFs and their binding affinity at putative CREs determine tissue- and developmental-specific transcriptional activity. Consolidating genomic datasets can offer further insights into the accessibility of CREs, TF activity, and, thus, gene regulation. However, the integration and analysis of multimodal datasets are hampered by considerable technical challenges. While methods for highlighting differential TF activity from combined chromatin state data (e.g., chromatin immunoprecipitation [ChIP], ATAC, or DNase sequencing) and RNA sequencing data exist, they do not offer convenient usability, have limited support for large-scale data processing, and provide only minimal functionality for visually interpreting results.

RESULTS

We developed TF-Prioritizer, an automated pipeline that prioritizes condition-specific TFs from multimodal data and generates an interactive web report. We demonstrated its potential by identifying known TFs along with their target genes, as well as previously unreported TFs active in lactating mouse mammary glands. Additionally, we studied a variety of ENCODE datasets for cell lines K562 and MCF-7, including 12 histone modification ChIP sequencing as well as ATAC and DNase sequencing datasets, where we observe and discuss assay-specific differences.

CONCLUSION

TF-Prioritizer accepts ATAC, DNase, or ChIP sequencing and RNA sequencing data as input and identifies TFs with differential activity, thus offering an understanding of genome-wide gene regulation, potential pathogenesis, and therapeutic targets in biomedical research.

Structural insights into p300 regulation and acetylation-dependent genome organisation

Histone modifications are deposited by chromatin modifying enzymes and read out by proteins that recognize the modified state. BRD4-NUT is an oncogenic fusion protein of the acetyl lysine reader BRD4 that binds to the acetylase p300 and enables formation of long-range intra- and interchromosomal interactions. We here examine how acetylation reading and writing enable formation of such interactions. We show that NUT contains an acidic transcriptional activation domain that binds to the TAZ2 domain of p300. We use NMR to investigate the structure of the complex and found that the TAZ2 domain has an autoinhibitory role for p300. NUT-TAZ2 interaction or mutations found in cancer that interfere with autoinhibition by TAZ2 allosterically activate p300. p300 activation results in a self-organizing, acetylation-dependent feed-forward reaction that enables long-range interactions by bromodomain multivalent acetyl-lysine binding. We discuss the implications for chromatin organisation, gene regulation and dysregulation in disease.

Transcription factor binding process is the primary driver of noise in gene expression

Noise in expression of individual genes gives rise to variations in activity of cellular pathways and generates heterogeneity in cellular phenotypes. Phenotypic heterogeneity has important implications for antibiotic persistence, mutation penetrance, cancer growth and therapy resistance. Specific molecular features such as the presence of the TATA box sequence and the promoter nucleosome occupancy have been associated with noise. However, the relative importance of these features in noise regulation is unclear and how well these features can predict noise has not yet been assessed. Here through an integrated statistical model of gene expression noise in yeast we found that the number of regulating transcription factors (TFs) of a gene was a key predictor of noise, whereas presence of the TATA box and the promoter nucleosome occupancy had poor predictive power. With an increase in the number of regulatory TFs, there was a rise in the number of cooperatively binding TFs. In addition, an increased number of regulatory TFs meant more overlaps in TF binding sites, resulting in competition between TFs for binding to the same region of the promoter. Through modeling of TF binding to promoter and application of stochastic simulations, we demonstrated that competition and cooperation among TFs could increase noise. Thus, our work uncovers a process of noise regulation that arises out of the dynamics of gene regulation and is not dependent on any specific transcription factor or specific promoter sequence.

Systematic analysis of low-affinity transcription factor binding site clusters in vitro and in vivo establishes their functional relevance

Binding to binding site clusters has yet to be characterized in depth, and the functional relevance of low-affinity clusters remains uncertain. We characterized transcription factor binding to low-affinity clusters in vitro and found that transcription factors can bind concurrently to overlapping sites, challenging the notion of binding exclusivity. Furthermore, small clusters with binding sites an order of magnitude lower in affinity give rise to high mean occupancies at physiologically-relevant transcription factor concentrations. To assess whether the observed in vitro occupancies translate to transcriptional activation in vivo, we tested low-affinity binding site clusters in a synthetic and native gene regulatory network in S. cerevisiae. In both systems, clusters of low-affinity binding sites generated transcriptional output comparable to single or even multiple consensus sites. This systematic characterization demonstrates that clusters of low-affinity binding sites achieve substantial occupancies, and that this occupancy can drive expression in eukaryotic promoters.

How to Completely Squeeze a Fungus-Advanced Genome Mining Tools for Novel Bioactive Substances

Fungal species have the capability of producing an overwhelming diversity of bioactive substances that can have beneficial but also detrimental effects on human health. These so-called secondary metabolites naturally serve as antimicrobial "weapon systems", signaling molecules or developmental effectors for fungi and hence are produced only under very specific environmental conditions or stages in their life cycle. However, as these complex conditions are difficult or even impossible to mimic in laboratory settings, only a small fraction of the true chemical diversity of fungi is known so far. This also implies that a large space for potentially new pharmaceuticals remains unexplored. We here present an overview on current developments in advanced methods that can be used to explore this chemical space. We focus on genetic and genomic methods, how to detect genes that harbor the blueprints for the production of these compounds (i.e., biosynthetic gene clusters, BGCs), and ways to activate these silent chromosomal regions. We provide an in-depth view of the chromatin-level regulation of BGCs and of the potential to use the CRISPR/Cas technology as an activation tool.

Concentration-dependent transcriptional switching through a collective action of cis-elements

Gene expression specificity of homeobox transcription factors has remained paradoxical. WUSCHEL activates and represses CLAVATA3 transcription at lower and higher concentrations, respectively. We use computational modeling and experimental analysis to investigate the properties of the cis-regulatory module. We find that intrinsically each cis-element can only activate CLAVATA3 at a higher WUSCHEL concentration. However, together, they repress CLAVATA3 at higher WUSCHEL and activate only at lower WUSCHEL, showing that the concentration-dependent interactions among cis-elements regulate both activation and repression. Biochemical experiments show that two adjacent functional cis-elements bind WUSCHEL with higher affinity and dimerize at relatively lower levels. Moreover, increasing the distance between cis-elements prolongs WUSCHEL monomer binding window, resulting in higher CLAVATA3 activation. Our work showing a constellation of optimally spaced cis-elements of defined affinities determining activation and repression thresholds in regulating CLAVATA3 transcription provides a previously unknown mechanism of cofactor-independent regulation of transcription factor binding in mediating gene expression specificity.

Profiling the quantitative occupancy of myriad transcription factors across conditions by modeling chromatin accessibility data

Over a thousand different transcription factors (TFs) bind with varying occupancy across the human genome. Chromatin immunoprecipitation (ChIP) can assay occupancy genome-wide, but only one TF at a time, limiting our ability to comprehensively observe the TF occupancy landscape, let alone quantify how it changes across conditions. We developed TF occupancy profiler (TOP), a Bayesian hierarchical regression framework, to profile genome-wide quantitative occupancy of numerous TFs using data from a single chromatin accessibility experiment (DNase- or ATAC-seq). TOP is supervised, and its hierarchical structure allows it to predict the occupancy of any sequence-specific TF, even those never assayed with ChIP. We used TOP to profile the quantitative occupancy of hundreds of sequence-specific TFs at sites throughout the genome and examined how their occupancies changed in multiple contexts: in approximately 200 human cell types, through 12 h of exposure to different hormones, and across the genetic backgrounds of 70 individuals. TOP enables cost-effective exploration of quantitative changes in the landscape of TF binding.

Facilitated Dissociation of Nucleoid Associated Proteins from DNA in the Bacterial Confinement

Transcription machinery depends on the temporal formation of protein-DNA complexes. Recent experiments demonstrated that not only the formation but also the lifetime of such complexes can affect the transcriptional machinery. In parallel, in vitro single-molecule studies showed that nucleoid-associated proteins (NAPs) leave the DNA rapidly as the bulk concentration of the protein increases via facilitated dissociation (FD). Nevertheless, whether such a concentration-dependent mechanism is functional in a bacterial cell, in which NAP levels and the 3d chromosomal structure are often coupled, is not clear a priori. Here, by using extensive coarse-grained molecular simulations, we model the unbinding of specific and nonspecific dimeric NAPs from a high-molecular-weight circular DNA molecule in a cylindrical structure mimicking the cellular confinement of a bacterial chromosome. Our simulations confirm that physiologically relevant peak protein levels (tens of micromolar) lead to highly compact chromosomal structures. This compaction results in rapid off rates (shorter DNA residence times) for specifically DNA-binding NAPs, such as the factor for inversion stimulation, which mostly dissociate via a segmental jump mechanism. Contrarily, for nonspecific NAPs, which are more prone to leave their binding sites via 1d sliding, the off rates decrease as the protein levels increase. The simulations with restrained chromosome models reveal that chromosome compaction is in favor of faster dissociation but only for specific proteins, and nonspecific proteins are not affected by the chromosome compaction. Overall, our results suggest that the cellular concentration level of a structural DNA-binding protein can be highly intermingled with its DNA residence time.

DNA-induced spatial entrapment of general transcription machinery can stabilize gene expression in a nondividing cell

How differentiated cells such as muscle or nerve maintain their gene expression for prolonged times is currently elusive. Here, using Xenopus oocyte, we have shown that the stability of gene expression in nondividing cells may arise due to the local entrapment of transcriptional machinery to specific gene transcription start sites. We found that within the same nucleus active versus inactive versions of the same gene are spatially segregated through liquid–liquid phase separation. We further observe that silent genes are associated with RNA-Pol-II phosphorylated on Ser5 but fails to attract RNA-Pol-II elongation factors. We propose that liquid–liquid phase separation mediated entrapment of limiting transcriptional machinery factors maintain stable expression of some genes in nondividing cells.

An important characteristic of cell differentiation is its stability. Only rarely do cells or their stem cell progenitors change their differentiation pathway. If they do, it is often accompanied by a malfunction such as cancer. A mechanistic understanding of the stability of differentiated states would allow better prospects of alleviating the malfunctioning. However, such complete information is yet elusive. Earlier experiments performed in Xenopus oocytes to address this question suggest that a cell may maintain its gene expression by prolonged binding of cell type–specific transcription factors. Here, using DNA competition experiments, we show that the stability of gene expression in a nondividing cell could be caused by the local entrapment of part of the general transcription machinery in transcriptionally active regions. Strikingly, we found that transcriptionally active and silent forms of the same DNA template can stably coexist within the same nucleus. Both DNA templates are associated with the gene-specific transcription factor Ascl1, the core factor TBP2, and the polymerase II (Pol-II) ser5 C-terminal domain (CTD) phosphorylated form, while Pol-II ser2 CTD phosphorylation is restricted to the transcriptionally dominant template. We discover that the active and silent DNA forms are physically separated in the oocyte nucleus through partition into liquid–liquid phase-separated condensates. Altogether, our study proposes a mechanism of transcriptional regulation involving a spatial entrapment of general transcription machinery components to stabilize the active form of a gene in a nondividing cell.

The initiation and maintenance of a differentiated state in development

Proper function of the body is maintained by an intricate interaction and communication among cells. during the animal development how these cells are formed and maintained is an important yet elusive. Understanding of how cells such as muscle and nerve cells maintain their identities would enable us to control the diseases which include malfunctioning in cellular identities such as cancer. In this article, we describe how the concept of formation and maintenance of cell identities has changed over the last 100 years. We will also briefly describe our current experimental work which includes transcriptional dynamics, and protein-protein interaction and how they are bringing new molecular insights. We also describe liquid-liquid phase separation as a potential new mechanism for the stability of gene expression in the non dvididng specialised cells of Xenopus oocytes.

DeNOPA: decoding nucleosome positions sensitively with sparse ATAC-seq data

As the basal bricks, the dynamics and arrangement of nucleosomes orchestrate the higher architecture of chromatin in a fundamental way, thereby affecting almost all nuclear biology processes. Thanks to its rather simple protocol, assay for transposase-accessible chromatin using sequencing (ATAC)-seq has been rapidly adopted as a major tool for chromatin-accessible profiling at both bulk and single-cell levels; however, to picture the arrangement of nucleosomes per se remains a challenge with ATAC-seq. In the present work, we introduce a novel ATAC-seq analysis toolkit, named decoding nucleosome organization profile based on ATAC-seq data (deNOPA), to predict nucleosome positions. Assessments showed that deNOPA outperformed state-of-the-art tools with ultra-sparse ATAC-seq data, e.g. no more than 0.5 fragment per base pair. The remarkable performance of deNOPA was fueled by the short fragment reads, which compose nearly half of sequenced reads in the ATAC-seq libraries and are commonly discarded by state-of-the-art nucleosome positioning tools. However, we found that the short fragment reads enrich information on nucleosome positions and that the linker regions were predicted by reads from both short and long fragments using Gaussian smoothing. Last, using deNOPA, we showed that the dynamics of nucleosome organization may not directly couple with chromatin accessibility in the cis-regulatory regions when human cells respond to heat shock stimulation. Our deNOPA provides a powerful tool with which to analyze the dynamics of chromatin at nucleosome position level with ultra-sparse ATAC-seq data.

TBX2 controls a proproliferative gene expression program in melanoma

Senescence shapes embryonic development, plays a key role in aging, and is a critical barrier to cancer initiation, yet how senescence is regulated remains incompletely understood. TBX2 is an antisenescence T-box family transcription repressor implicated in embryonic development and cancer. However, the repertoire of TBX2 target genes, its cooperating partners, and how TBX2 promotes proliferation and senescence bypass are poorly understood. Here, using melanoma as a model, we show that TBX2 lies downstream from PI3K signaling and that TBX2 binds and is required for expression of E2F1, a key antisenescence cell cycle regulator. Remarkably, TBX2 binding in vivo is associated with CACGTG E-boxes, present in genes down-regulated by TBX2 depletion, more frequently than the consensus T-element DNA binding motif that is restricted to Tbx2 repressed genes. TBX2 is revealed to interact with a wide range of transcription factors and cofactors, including key components of the BCOR/PRC1.1 complex that are recruited by TBX2 to the E2F1 locus. Our results provide key insights into how PI3K signaling modulates TBX2 function in cancer to drive proliferation.

Different SP1 binding dynamics at individual genomic loci in human cells

Using a tamoxifen-inducible time-course ChIP-sequencing (ChIP-seq) approach, we show that the ubiquitous transcription factor SP1 has different binding dynamics at its target sites in the human genome. SP1 very rapidly reaches maximal binding levels at some sites, but binding kinetics at other sites is biphasic, with rapid half-maximal binding followed by a considerably slower increase to maximal binding. While ∼70% of SP1 binding sites are located at promoter regions, loci with slow SP1 binding kinetics are enriched in enhancer and Polycomb-repressed regions. Unexpectedly, SP1 sites with fast binding kinetics tend to have higher quality and more copies of the SP1 sequence motif. Different cobinding factors associate near SP1 binding sites depending on their binding kinetics and on their location at promoters or enhancers. For example, NFY and FOS are preferentially associated near promoter-bound SP1 sites with fast binding kinetics, whereas DNA motifs of ETS and homeodomain proteins are preferentially observed at sites with slow binding kinetics. At promoters but not enhancers, proteins involved in sumoylation and PML bodies associate more strongly with slow SP1 binding sites than with the fast binding sites. The speed of SP1 binding is not associated with nucleosome occupancy, and it is not necessarily coupled to higher transcriptional activity. These results with SP1 are in contrast to those of human TBP, indicating that there is no common mechanism affecting transcription factor binding kinetics. The biphasic kinetics at some SP1 target sites suggest the existence of distinct chromatin states at these loci in different cells within the overall population.

Single molecule characterization of the binding kinetics of a transcription factor and its modulation by DNA sequence and methylation

The interaction of transcription factors with their response elements in DNA is emerging as a highly complex process, whose characterization requires measuring the full distribution of binding and dissociation times in a well-controlled assay. Here, we present a single-molecule assay that exploits the thermal fluctuations of a DNA hairpin to detect the association and dissociation of individual, unlabeled transcription factors. We demonstrate this new approach by following the binding of Egr1 to its consensus motif and the three binding sites found in the promoter of the Lhb gene, and find that both association and dissociation are modulated by the 9 bp core motif and the sequences around it. In addition, CpG methylation modulates the dissociation kinetics in a sequence and position-dependent manner, which can both stabilize or destabilize the complex. Together, our findings show how variations in sequence and methylation patterns synergistically extend the spectrum of a protein's binding properties, and demonstrate how the proposed approach can provide new insights on the function of transcription factors.

Transcription factor binding kinetics and transcriptional bursting: What do we really know?

In eukaryotes, transcription is a discontinuous process with mRNA being generated in bursts, after the binding of transcription factors (TFs) to regulatory elements on the genome. Live-cell single-molecule microscopy has highlighted that transcriptional bursting can be controlled by tuning TF/DNA binding kinetics. Yet the timescales of these two processes seem disconnected with TF/DNA interactions typically lasting orders of magnitude shorter than transcriptional bursts. To test models that could reconcile these discrepancies, reliable measurements of TF binding kinetics are needed, also accounting for the current limitations in performing these single-molecule measurements at specific regulatory elements. Here, we review the recent studies linking TF binding kinetics to transcriptional bursting and outline some current and future challenges that need to be addressed to provide a microscopic description of transcriptional regulation kinetics.

Exact Probability Landscapes of Stochastic Phenotype Switching in Feed-Forward Loops: Phase Diagrams of Multimodality

Feed-forward loops (FFLs) are among the most ubiquitously found motifs of reaction networks in nature. However, little is known about their stochastic behavior and the variety of network phenotypes they can exhibit. In this study, we provide full characterizations of the properties of stochastic multimodality of FFLs, and how switching between different network phenotypes are controlled. We have computed the exact steady-state probability landscapes of all eight types of coherent and incoherent FFLs using the finite-butter Accurate Chemical Master Equation (ACME) algorithm, and quantified the exact topological features of their high-dimensional probability landscapes using persistent homology. Through analysis of the degree of multimodality for each of a set of 10,812 probability landscapes, where each landscape resides over 10⁵–10⁶ microstates, we have constructed comprehensive phase diagrams of all relevant behavior of FFL multimodality over broad ranges of input and regulation intensities, as well as different regimes of promoter binding dynamics. In addition, we have quantified the topological sensitivity of the multimodality of the landscapes to regulation intensities. Our results show that with slow binding and unbinding dynamics of transcription factor to promoter, FFLs exhibit strong stochastic behavior that is very different from what would be inferred from deterministic models. In addition, input intensity play major roles in the phenotypes of FFLs: At weak input intensity, FFL exhibit monomodality, but strong input intensity may result in up to 6 stable phenotypes. Furthermore, we found that gene duplication can enlarge stable regions of specific multimodalities and enrich the phenotypic diversity of FFL networks, providing means for cells toward better adaptation to changing environment. Our results are directly applicable to analysis of behavior of FFLs in biological processes such as stem cell differentiation and for design of synthetic networks when certain phenotypic behavior is desired.

Time-Based Systems Biology Approaches to Capture and Model Dynamic Gene Regulatory Networks

All aspects of transcription and its regulation involve dynamic events. However, capturing these dynamic events in gene regulatory networks (GRNs) offers both a promise and a challenge. The promise is that capturing and modeling the dynamic changes in GRNs will allow us to understand how organisms adapt to a changing environment. The ability to mount a rapid transcriptional response to environmental changes is especially important in nonmotile organisms such as plants. The challenge is to capture these dynamic, genome-wide events and model them in GRNs. In this review, we cover recent progress in capturing dynamic interactions of transcription factors with their targets-at both the local and genome-wide levels-and how they are used to learn how GRNs operate as a function of time. We also discuss recent advances that employ time-based machine learning approaches to forecast gene expression at future time points, a key goal of systems biology.

Runx proteins and transcriptional mechanisms that govern memory CD8 T cell development

Adaptive immunity to intracellular pathogens and tumors is mediated by antigen-experienced CD8 T cells. Individual naive CD8 T cells have the potential to differentiate into a diverse array of antigen-experienced subsets that exhibit distinct effector functions, life spans, anatomic positioning, and potential for regenerating an entirely new immune response during iterative pathogenic exposures. The developmental process by which activated naive cells undergo diversification involves regulation of chromatin structure and transcription but is not entirely understood. This review examines how alterations in chromatin structure, transcription factor binding, extracellular signals, and single-cell gene expression explain the differential development of distinct effector (T_EFF ) and memory (T_MEM ) CD8 T cell subsets. Special emphasis is placed on how Runx proteins function with additional transcription factors to pioneer changes in chromatin accessibility and drive transcriptional programs that establish the core attributes of cytotoxic T lymphocytes, subdivide circulating and non-circulating T_MEM cell subsets, and govern terminal differentiation. The discussion integrates the roles of specific cytokine signals, transcriptional circuits and how regulation of individual nucleosomes and RNA polymerase II activity can contribute to the process of differentiation. A model that integrates many of these features is discussed to conceptualize how activated CD8 T cells arrive at their fates.

Chromatin information content landscapes inform transcription factor and DNA interactions

Interactions between transcription factors and chromatin are fundamental to genome organization and regulation and, ultimately, cell state. Here, we use information theory to measure signatures of organized chromatin resulting from transcription factor-chromatin interactions encoded in the patterns of the accessible genome, which we term chromatin information enrichment (CIE). We calculate CIE for hundreds of transcription factor motifs across human samples and identify two classes: low and high CIE. The 10-20% of common and tissue-specific high CIE transcription factor motifs, associate with higher protein-DNA residence time, including different binding site subclasses of the same transcription factor, increased nucleosome phasing, specific protein domains, and the genetic control of both chromatin accessibility and gene expression. These results show that variations in the information encoded in chromatin architecture reflect functional biological variation, with implications for cell state dynamics and memory.

SUMO and Transcriptional Regulation: The Lessons of Large-Scale Proteomic, Modifomic and Genomic Studies

One major role of the eukaryotic peptidic post-translational modifier SUMO in the cell is transcriptional control. This occurs via modification of virtually all classes of transcriptional actors, which include transcription factors, transcriptional coregulators, diverse chromatin components, as well as Pol I-, Pol II- and Pol III transcriptional machineries and their regulators. For many years, the role of SUMOylation has essentially been studied on individual proteins, or small groups of proteins, principally dealing with Pol II-mediated transcription. This provided only a fragmentary view of how SUMOylation controls transcription. The recent advent of large-scale proteomic, modifomic and genomic studies has however considerably refined our perception of the part played by SUMO in gene expression control. We review here these developments and the new concepts they are at the origin of, together with the limitations of our knowledge. How they illuminate the SUMO-dependent transcriptional mechanisms that have been characterized thus far and how they impact our view of SUMO-dependent chromatin organization are also considered.

Learning causal networks using inducible transcription factors and transcriptome-wide time series

We present IDEA (the Induction Dynamics gene Expression Atlas), a dataset constructed by independently inducing hundreds of transcription factors (TFs) and measuring timecourses of the resulting gene expression responses in budding yeast. Each experiment captures a regulatory cascade connecting a single induced regulator to the genes it causally regulates. We discuss the regulatory cascade of a single TF, Aft1, in detail; however, IDEA contains > 200 TF induction experiments with 20 million individual observations and 100,000 signal‐containing dynamic responses. As an application of IDEA, we integrate all timecourses into a whole‐cell transcriptional model, which is used to predict and validate multiple new and underappreciated transcriptional regulators. We also find that the magnitudes of coefficients in this model are predictive of genetic interaction profile similarities. In addition to being a resource for exploring regulatory connectivity between TFs and their target genes, our modeling approach shows that combining rapid perturbations of individual genes with genome‐scale time‐series measurements is an effective strategy for elucidating gene regulatory networks.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Genome-wide off-rates reveal how DNA binding dynamics shape transcription factor function

Protein-DNA interactions are dynamic, and these dynamics are an important aspect of chromatin-associated processes such as transcription or replication. Due to a lack of methods to study on- and off-rates across entire genomes, protein-DNA interaction dynamics have not been studied extensively. Here, we determine in vivo off-rates for the Saccharomyces cerevisiae chromatin organizing factor Abf1, at 191 sites simultaneously across the yeast genome. Average Abf1 residence times span a wide range, varying between 4.2 and 33 min. Sites with different off-rates are associated with different functional characteristics. This includes their transcriptional dependency on Abf1, nucleosome positioning and the size of the nucleosome-free region, as well as the ability to roadblock RNA polymerase II for termination. The results show how off-rates contribute to transcription factor function and that DIVORSEQ (Determining In Vivo Off-Rates by SEQuencing) is a meaningful way of investigating protein-DNA binding dynamics genome-wide.

Single molecule studies reveal that p53 tetramers dynamically bind response elements containing one or two half sites

The tumor suppressor protein p53 is critical for cell fate decisions, including apoptosis, senescence, and cell cycle arrest. p53 is a tetrameric transcription factor that binds DNA response elements to regulate transcription of target genes. p53 response elements consist of two decameric half-sites, and data suggest one p53 dimer in the tetramer binds to each half-site. Despite a broad literature describing p53 binding DNA, unanswered questions remain, due partly to the need for more quantitative and structural studies with full length protein. Here we describe a single molecule fluorescence system to visualize full length p53 tetramers binding DNA in real time. The data revealed a dynamic interaction in which tetrameric p53/DNA complexes assembled and disassembled without a dimer/DNA intermediate. On a wild type DNA containing two half sites, p53/DNA complexes existed in two kinetically distinct populations. p53 tetramers bound response elements containing only one half site to form a single population of complexes with reduced kinetic stability. Altering the spacing and helical phasing between two half sites affected both the population distribution of p53/DNA complexes and their kinetic stability. Our real time single molecule measurements of full length p53 tetramers binding DNA reveal the parameters that define the stability of p53/DNA complexes, and provide insight into the pathways by which those complexes assemble.

Optogenetic Control Reveals Differential Promoter Interpretation of Transcription Factor Nuclear Translocation Dynamics

Gene expression is thought to be affected not only by the concentration of transcription factors (TFs) but also the dynamics of their nuclear translocation. Testing this hypothesis requires direct control of TF dynamics. Here, we engineer CLASP, an optogenetic tool for rapid and tunable translocation of a TF of interest. Using CLASP fused to Crz1, we observe that, for the same integrated concentration of nuclear TF over time, changing input dynamics changes target gene expression: pulsatile inputs yield higher expression than continuous inputs, or vice versa, depending on the target gene. Computational modeling reveals that a dose-response saturating at low TF input can yield higher gene expression for pulsatile versus continuous input, and that multi-state promoter activation can yield the opposite behavior. Our integrated tool development and modeling approach characterize promoter responses to Crz1 nuclear translocation dynamics, extracting quantitative features that may help explain the differential expression of target genes.

CLASP is a modular optogenetic strategy to control the nuclear localization of transcription factors (TFs) and elicit gene expression from their cognate promoters. CLASP control of Crz1 nuclear localization, coupled with computational modeling, revealed how promoters can differentially decode dynamic transcription factor signals. The integrated strategy of CLASP development and modeling presents a generalized approach to causally investigate the transcriptional consequences of dynamic TF nuclear shuttling.

Long-term association of a transcription factor with its chromatin binding site can stabilize gene expression and cell fate commitment

Some kinds of transcription factor proteins are very important in initiating and guiding cell fate differentiation. Overexpression of these factors can force many other kinds of cells to become muscle or nerve. Examples are MyoD for muscle and Ascl1 for nerve. It is not known how long such a factor must remain bound to its binding site for it to have its function; this could be seconds, minutes, hours, or days. We have developed a procedure to determine the required residence time for the Ascl1 nerve factor to have its function. This factor remains closely associated with its chromatin binding site for hours or days. This may be a general characteristic of such factors in nondividing (adult) cells.

Some lineage-determining transcription factors are overwhelmingly important in directing embryonic cells to a particular differentiation pathway, such as Ascl1 for nerve. They also have an exceptionally strong ability to force cells to change from an unrelated pathway to one preferred by their action. Transcription factors are believed to have a very short residence time of only a few seconds on their specific DNA or chromatin-binding sites. We have developed a procedure in which DNA containing one copy of the binding site for the neural-inducing factor Ascl1 is injected directly into a Xenopus oocyte nucleus which has been preloaded with a limiting amount of the Ascl1 transcription factor protein. This is followed by a further injection of DNA as a competitor, either in a plasmid or in chromosomal DNA, containing the same binding site but with a different reporter. Importantly, expression of the reporter provides a measure of the function of the transcription factor in addition to its residence time. The same long residence time and resistance to competition are seen with the estrogen receptor and its DNA response elements. We find that in this nondividing oocyte, the nerve-inducing factor Ascl1 can remain bound to a specific chromatin site for hours or days and thereby help to stabilize gene expression. This stability of transcription factor binding to chromatin is a necessary part of its action because removal of this factor causes discontinuation of its effect on gene expression. Stable transcription factor binding may be a characteristic of nondividing cells.

Tuning Transcription Factor Availability through Acetylation-Mediated Genomic Redistribution

It is widely assumed that decreasing transcription factor DNA-binding affinity reduces transcription initiation by diminishing occupancy of sequence-specific regulatory elements. However, in vivo transcription factors find their binding sites while confronted with a large excess of low-affinity degenerate motifs. Here, using the melanoma lineage survival oncogene MITF as a model, we show that low-affinity binding sites act as a competitive reservoir in vivo from which transcription factors are released by mitogen-activated protein kinase (MAPK)-stimulated acetylation to promote increased occupancy of their regulatory elements. Consequently, a low-DNA-binding-affinity acetylation-mimetic MITF mutation supports melanocyte development and drives tumorigenesis, whereas a high-affinity non-acetylatable mutant does not. The results reveal a paradoxical acetylation-mediated molecular clutch that tunes transcription factor availability via genome-wide redistribution and couples BRAF to tumorigenesis. Our results further suggest that p300/CREB-binding protein-mediated transcription factor acetylation may represent a common mechanism to control transcription factor availability.

•

Reducing transcription factor DNA-binding affinity increases activity in vivo

•

Acetylation is triggered by MAPK signaling

•

Acetylation leads to genome-wide transcription factor redistribution

•

Acetylation of MITF drives tumorigenesis and melanocyte development

Short Residence Times of DNA-Bound Transcription Factors Can Reduce Gene Expression Noise and Increase the Transmission of Information in a Gene Regulation System

Gene expression noise is not just ubiquitous but also variable, and we still do not understand some of the most elementary factors that affect it. Among them is the residence time of a transcription factor (TF) on DNA, the mean time that a DNA-bound TF remains bound. Here, we use a stochastic model of transcriptional regulation to study how residence time affects the gene expression noise that arises when a TF induces gene expression. We find that the effect of residence time on gene expression noise depends on the TF’s concentration and its affinity to DNA, which determine the level of induction of the gene. At high levels of induction, residence time has no effect on gene expression noise. However, as the level of induction decreases, short residence times reduce gene expression noise. The reason is that fast on-off TF binding dynamics prevent long periods where proteins are predominantly synthesized or degraded, which can cause excessive fluctuations in gene expression. As a consequence, short residence times can help a gene regulation system acquire information about the cellular environment it operates in. Our predictions are consistent with the observation that experimentally measured residence times are usually modest and lie between seconds to minutes.

Tuning up Transcription Factors for Therapy

The recent developments in the delivery and design of transcription factors put their therapeutic applications within reach, exemplified by cell replacement, cancer differentiation and T-cell based cancer therapies. The success of such applications depends on the efficacy and precision in the action of transcription factors. The biophysical and genetic characterization of the paradigmatic prokaryotic repressors, LacI and TetR and the designer transcription factors, transcription activator-like effector (TALE) and CRISPR-dCas9 revealed common principles behind their efficacy, which can aid the optimization of transcriptional activators and repressors. Further studies will be required to analyze the linkage between dissociation constants and enzymatic activity, the role of phase separation and squelching in activation and repression and the long-range interaction of transcription factors with epigenetic regulators in the context of the chromosomes. Understanding these mechanisms will help to tailor natural and synthetic transcription factors to the needs of specific applications.

From reads to insight: a hitchhiker's guide to ATAC-seq data analysis

Assay of Transposase Accessible Chromatin sequencing (ATAC-seq) is widely used in studying chromatin biology, but a comprehensive review of the analysis tools has not been completed yet. Here, we discuss the major steps in ATAC-seq data analysis, including pre-analysis (quality check and alignment), core analysis (peak calling), and advanced analysis (peak differential analysis and annotation, motif enrichment, footprinting, and nucleosome position analysis). We also review the reconstruction of transcriptional regulatory networks with multiomics data and highlight the current challenges of each step. Finally, we describe the potential of single-cell ATAC-seq and highlight the necessity of developing ATAC-seq specific analysis tools to obtain biologically meaningful insights.

Rap1-mediated nucleosome displacement can regulate gene expression in senescent cells without impacting the pace of senescence

Cell senescence is accompanied, and in part mediated, by changes in chromatin, including histone losses, but underlying mechanisms are not well understood. We reported previously that during yeast cell senescence driven by telomere shortening, the telomeric protein Rap1 plays a major role in reprogramming gene expression by relocalizing hundreds of new target genes (called NRTS, for new Rap1 targets at senescence) to the promoters. This leads to two types of histone loss: Rap1 lowers histone level globally by repressing histone gene expression, and it also causes local nucleosome displacement at the promoters of upregulated NRTS. Here, we present evidence of direct binding between Rap1 and histone H3/H4 heterotetramers, and map amino acids involved in the interaction within the Rap1 SANT domain to amino acids 392-394 (SHY). Introduction of a point mutation within the native RAP1 locus that converts these residues to alanines (RAP1^SHY ), and thus disrupts Rap1-H3/H4 interaction, does not interfere with Rap1 relocalization to NRTS at senescence, but prevents full nucleosome displacement and gene upregulation, indicating direct Rap1-H3/H4 contacts are involved in nucleosome displacement. Consistent with this, the histone H3/H4 chaperone Asf1 is similarly unnecessary for Rap1 localization to NRTS but is required for full Rap1-mediated nucleosome displacement and gene activation. Remarkably, RAP1^SHY does not affect the pace of senescence-related cell cycle arrest, indicating that some changes in gene expression at senescence are not coupled to this arrest.

Lessons from eRNAs: understanding transcriptional regulation through the lens of nascent RNAs

Nascent transcription assays, such as global run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq), have uncovered a myriad of unstable RNAs being actively produced from numerous sites genome-wide. These transcripts provide a more complete and immediate picture of the impact of regulatory events. Transcription factors recruit RNA polymerase II, effectively initiating the process of transcription; repressors inhibit polymerase recruitment. Efficiency of recruitment is dictated by sequence elements in and around the RNA polymerase loading zone. A combination of sequence elements and RNA binding proteins subsequently influence the ultimate stability of the resulting transcript. Some of these transcripts are capable of providing feedback on the process, influencing subsequent transcription. By monitoring RNA polymerase activity, nascent assays provide insights into every step of the regulated process of transcription.