Nucleosome retention and the stochastic nature of promoter chromatin remodeling for transcription

Exact burst-size distributions for gene-expression models with complex promoter structure

In prokaryotic and eukaryotic cells, most genes are transcribed in a bursty fashion on one hand and complex gene regulations may lead to complex promoter structure on the other hand. This raises an unsolved issue: how does promoter structure shape transcriptional bursting kinetics characterized by burst size and frequency? Here we analyze stochastic models of gene transcription, which consider complex regulatory mechanisms. Notably, we develop an efficient method to derive exact burst-size distributions. The analytical results show that if the promoter of a gene contains only one active state, the burst size indeed follows a geometric distribution, in agreement with the previous result derived under certain limiting conditions. However, if it contains a multitude of active states, the burst size in general obeys a non-geometric distribution, which is a linearly weighted sum of geometric distributions. This superposition principle reveals the essential feature of bursting kinetics in complex cases of transcriptional regulation although it seems that there has been no direct experimental confirmation. The derived burst-size distributions not only highlight the importance of promoter structure in regulating bursting kinetics, but can be also used in the exact inference of this kinetics based on experimental data.

Review and Evaluate the Bioinformatics Analysis Strategies of ATAC-seq and CUT&Tag Data

Efficient and reliable profiling methods are essential to study epigenetics. Tn5, one of the first identified prokaryotic transposases with high DNA-binding and tagmentation efficiency, is widely adopted in different genomic and epigenomic protocols for high-throughputly exploring the genome and epigenome. Based on Tn5, the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) and the Cleavage Under Targets and Tagmentation (CUT&Tag) were developed to measure chromatin accessibility and detect DNA-protein interactions. These methodologies can be applied to large amounts of biological samples with low-input levels, such as rare tissues, embryos, and sorted single cells. However, fast and proper processing of these epigenomic data has become a bottleneck because massive data production continues to increase quickly. Furthermore, inappropriate data analysis can generate biased or misleading conclusions. Therefore, it is essential to evaluate the performance of Tn5-based ATAC-seq and CUT&Tag data processing bioinformatics tools, many of which were developed mostly for analyzing chromatin immunoprecipitation followed by sequencing (ChIP-seq) data. Here, we conducted a comprehensive benchmarking analysis to evaluate the performance of eight popular software for processing ATAC-seq and CUT&Tag data. We compared the sensitivity, specificity, and peak width distribution for both narrow-type and broad-type peak calling. We also tested the influence of the availability of control IgG input in CUT&Tag data analysis. Finally, we evaluated the differential analysis strategies commonly used for analyzing the CUT&Tag data. Our study provided comprehensive guidance for selecting bioinformatics tools and recommended analysis strategies, which were implemented into Docker/Singularity images for streamlined data analysis.

Genotype inference from aggregated chromatin accessibility data reveals genetic regulatory mechanisms

Background

Understanding the genetic causes for variability in chromatin accessibility can shed light on the molecular mechanisms through which genetic variants may affect complex traits. Thousands of ATAC-seq samples have been collected that hold information about chromatin accessibility across diverse cell types and contexts, but most of these are not paired with genetic information and come from diverse distinct projects and laboratories.

Results

We report here joint genotyping, chromatin accessibility peak calling, and discovery of quantitative trait loci which influence chromatin accessibility (caQTLs), demonstrating the capability of performing caQTL analysis on a large scale in a diverse sample set without pre-existing genotype information. Using 10,293 profiling samples representing 1,454 unique donor individuals across 653 studies from public databases, we catalog 23,381 caQTLs in total. After joint discovery analysis, we cluster samples based on accessible chromatin profiles to identify context-specific caQTLs. We find that caQTLs are strongly enriched for annotations of gene regulatory elements across diverse cell types and tissues and are often strongly linked with genetic variation associated with changes in expression (eQTLs), indicating that caQTLs can mediate genetic effects on gene expression. We demonstrate sharing of causal variants for chromatin accessibility and diverse complex human traits, enabling a more complete picture of the genetic mechanisms underlying complex human phenotypes.

Conclusions

Our work provides a proof of principle for caQTL calling from previously ungenotyped samples, and represents one of the largest, most diverse caQTL resources currently available, informing mechanisms of genetic regulation of gene expression and contribution to disease.

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.

Energy-driven genome regulation by ATP-dependent chromatin remodellers

The packaging of DNA into chromatin in eukaryotes regulates gene transcription, DNA replication and DNA repair. ATP-dependent chromatin remodelling enzymes (re)arrange nucleosomes at the first level of chromatin organization. Their Snf2-type motor ATPases alter histone-DNA interactions through a common DNA translocation mechanism. Whether remodeller activities mainly catalyse nucleosome dynamics or accurately co-determine nucleosome organization remained unclear. In this Review, we discuss the emerging mechanisms of chromatin remodelling: dynamic remodeller architectures and their interactions, the inner workings of the ATPase cycle, allosteric regulation and pathological dysregulation. Recent mechanistic insights argue for a decisive role of remodellers in the energy-driven self-organization of chromatin, which enables both stability and plasticity of genome regulation - for example, during development and stress. Different remodellers, such as members of the SWI/SNF, ISWI, CHD and INO80 families, process (epi)genetic information through specific mechanisms into distinct functional outputs. Combinatorial assembly of remodellers and their interplay with histone modifications, histone variants, DNA sequence or DNA-bound transcription factors regulate nucleosome mobilization or eviction or histone exchange. Such input-output relationships determine specific nucleosome positions and compositions with distinct DNA accessibilities and mediate differential genome regulation. Finally, remodeller genes are often mutated in diseases characterized by genome dysregulation, notably in cancer, and we discuss their physiological relevance.

Genome-scale chromatin interaction dynamic measurements for key components of the RNA Pol II general transcription machinery

Background

A great deal of work has revealed in structural detail the components of the machinery responsible for mRNA gene transcription initiation. These include the general transcription factors (GTFs), which assemble at promoters along with RNA Polymerase II (Pol II) to form a preinitiation complex (PIC) aided by the activities of cofactors and site-specific transcription factors (TFs). However, less well understood are the in vivo PIC assembly pathways and their kinetics, an understanding of which is vital for determining on a mechanistic level how rates of in vivo RNA synthesis are established and how cofactors and TFs impact them.

Results

We used competition ChIP to obtain genome-scale estimates of the residence times for five GTFs: TBP, TFIIA, TFIIB, TFIIE and TFIIF in budding yeast. While many GTF-chromatin interactions were short-lived (< 1 min), there were numerous interactions with residence times in the several minutes range. 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.

Conclusions

The datasets and results reported here provide kinetic information for most of the Pol II-driven genes in this organism and therefore offer a rich resource for exploring the mechanistic relationships between PIC assembly, gene regulation, and transcription. The relationships between gene function and GTF dynamics suggest that shared sets of TFs tune PIC assembly kinetics to ensure appropriate levels of expression.

Gene expression model inference from snapshot RNA data using Bayesian non-parametrics

Gene expression models, which are key towards understanding cellular regulatory response, underlie observations of single-cell transcriptional dynamics. Although RNA expression data encode information on gene expression models, existing computational frameworks do not perform simultaneous Bayesian inference of gene expression models and parameters from such data. Rather, gene expression models-composed of gene states, their connectivities and associated parameters-are currently deduced by pre-specifying gene state numbers and connectivity before learning associated rate parameters. Here we propose a method to learn full distributions over gene states, state connectivities and associated rate parameters, simultaneously and self-consistently from single-molecule RNA counts. We propagate noise from fluctuating RNA counts over models by treating models themselves as random variables. We achieve this within a Bayesian non-parametric paradigm. We demonstrate our method on the Escherichia coli lacZ pathway and the Saccharomyces cerevisiae STL1 pathway, and verify its robustness on synthetic data.

UL135 and UL136 Epistasis Controls Reactivation of Human Cytomegalovirus

Human cytomegalovirus (HCMV) is beta herpesvirus that persists indefinitely in the human host through a protracted, latent infection. The polycistronic UL133-UL138 gene locus of HCMV encodes genes regulating latency and reactivation. While UL138 is pro-latency, restricting virus replication in CD34+ hematopoietic progenitor cells (HPCs), UL135 overcomes this restriction for reactivation. By contrast, UL136 is expressed with later kinetics and encodes multiple protein isoforms with differential roles in latency and reactivation. Like UL135, the largest UL136 isoform, UL136p33, is required for reactivation from latency in hematopoietic cells. Furthermore, UL136p33 is unstable, and its instability is important for the establishment of latency and sufficient accumulation of UL136p33 is a checkpoint for reactivation. We hypothesized that stabilizing UL136p33 might overcome the requirement of UL135 for reactivation. To test this, we generated recombinant viruses lacking UL135 that expressed a stabilized variant of UL136p33. Stabilizing UL136p33 did not impact replication of the UL135-mutant virus in fibroblasts. However, in the context of infection in hematopoietic cells, stabilization of UL136p33 strikingly compensated for the loss of UL135, resulting in increased replication in CD34+ HPCs and in humanized NOD- scid IL2Rγ _c ^null (NSG) mice. This finding suggests that while UL135 is essential for reactivation, it functions at steps preceding the accumulation of UL136p33 and that stabilized expression of UL136p33 largely overcomes the requirement for UL135 in reactivation. Taken together, our genetic evidence indicates an epistatic relationship between UL136p33 and UL135 whereby UL135 may initiate events early in reactivation that will result in the accumulation of UL136p33 to a threshold required for productive reactivation.

SIGNIFICANCE

Human cytomegalovirus (HCMV) is one of nine human herpesviruses and a significant human pathogen. While HCMV establishes a life-long latent infection that is typically asymptomatic in healthy individuals, its reactivation from latency can have devastating consequences in the immune compromised. Defining virus-host and virus-virus interactions important for HCMV latency, reactivation and replication is critical to defining the molecular basis of latent and replicative states and in controlling infection and CMV disease. Here we define a genetic relationship between two viral genes in controlling virus reactivation from latency using primary human hematopoietic progenitor cell and humanized mouse models.

Molecular perspectives in hypertrophic heart disease: An epigenetic approach from chromatin modification

Epigenetic changes induced by environmental factors are increasingly relevant in cardiovascular diseases. The most frequent molecular component in cardiac hypertrophy is the reactivation of fetal genes caused by various pathologies, including obesity, arterial hypertension, aortic valve stenosis, and congenital causes. Despite the multiple investigations performed to achieve information about the molecular components of this pathology, its influence on therapeutic strategies is relatively scarce. Recently, new information has been taken about the proteins that modify the expression of fetal genes reactivated in cardiac hypertrophy. These proteins modify the DNA covalently and induce changes in the structure of chromatin. The relationship between histones and DNA has a recognized control in the expression of genes conditioned by the environment and induces epigenetic variations. The epigenetic modifications that regulate pathological cardiac hypertrophy are performed through changes in genomic stability, chromatin architecture, and gene expression. Histone 3 trimethylation at lysine 4, 9, or 27 (H3-K4; -K9; -K27me3) and histone demethylation at lysine 9 and 79 (H3-K9; -K79) are mediators of reprogramming in pathologic hypertrophy. Within the chromatin architecture modifiers, histone demethylases are a group of proteins that have been shown to play an essential role in cardiac cell differentiation and may also be components in the development of cardiac hypertrophy. In the present work, we review the current knowledge about the influence of epigenetic modifications in the expression of genes involved in cardiac hypertrophy and its possible therapeutic approach.

Kinetic Proofreading

Biochemistry and molecular biology rely on the recognition of structural complementarity between molecules. Molecular interactions must be both quickly reversible, i.e., tenuous, and specific. How the cell reconciles these conflicting demands is the subject of this article. The problem and its theoretical solution are discussed within the wider theoretical context of the thermodynamics of stochastic processes (stochastic thermodynamics). The solution-an irreversible reaction cycle that decreases internal error at the expense of entropy export into the environment-is shown to be widely employed by biological processes that transmit genetic and regulatory information. Expected final online publication date for the Annual Review of Biochemistry, Volume 91 is June 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The ATP-dependent SWI/SNF and RSC chromatin remodelers cooperatively induce unfolded protein response genes during endoplasmic reticulum stress

The SWI/SNF subfamily remodelers (SWI/SNF and RSC) generally promote gene expression by displacing or evicting nucleosomes at the promoter regions. Their action creates a nucleosome-depleted region where transcription machinery accesses the DNA. Their function has been shown critical for inducing stress-responsive transcription programs. Although the role of SWI/SNF and RSC complexes in transcription regulation of heat shock responsive genes is well studied, their involvement in other pathways such as unfolded protein response (UPR) and protein quality control (PQC) is less known. This study shows that SWI/SNF occupies the promoters of UPR, HSP and PQC genes in response to unfolded protein stress, and its recruitment at UPR promoters depends on Hac1 transcription factor and other epigenetic factors like Ada2 and Ume6. Disruption of SWI/SNF's activity does not affect the remodeling of these promoters or gene expression. However, inactivation of RSC and SWI/SNF together diminishes induction of most of the UPR, HSP and PQC genes tested. Furthermore, RSC and SWI/SNF colocalize at these promoters, suggesting that these two remodelers functionally cooperate to induce stress-responsive genes under proteotoxic conditions.

Sophisticated Conversations between Chromatin and Chromatin Remodelers, and Dissonances in Cancer

The establishment and maintenance of genome packaging into chromatin contribute to define specific cellular identity and function. Dynamic regulation of chromatin organization and nucleosome positioning are critical to all DNA transactions-in particular, the regulation of gene expression-and involve the cooperative action of sequence-specific DNA-binding factors, histone modifying enzymes, and remodelers. Remodelers are molecular machines that generate various chromatin landscapes, adjust nucleosome positioning, and alter DNA accessibility by using ATP binding and hydrolysis to perform DNA translocation, which is highly regulated through sophisticated structural and functional conversations with nucleosomes. In this review, I first present the functional and structural diversity of remodelers, while emphasizing the basic mechanism of DNA translocation, the common regulatory aspects, and the hand-in-hand progressive increase in complexity of the regulatory conversations between remodelers and nucleosomes that accompanies the increase in challenges of remodeling processes. Next, I examine how, through nucleosome positioning, remodelers guide the regulation of gene expression. Finally, I explore various aspects of how alterations/mutations in remodelers introduce dissonance into the conversations between remodelers and nucleosomes, modify chromatin organization, and contribute to oncogenesis.

TATA and paused promoters active in differentiated tissues have distinct expression characteristics

Core promoter types differ in the extent to which RNA polymerase II (Pol II) pauses after initiation, but how this affects their tissue-specific gene expression characteristics is not well understood. While promoters with Pol II pausing elements are active throughout development, TATA promoters are highly active in differentiated tissues. We therefore used a genomics approach on late-stage Drosophila embryos to analyze the properties of promoter types. Using tissue-specific Pol II ChIP-seq, we found that paused promoters have high levels of paused Pol II throughout the embryo, even in tissues where the gene is not expressed, while TATA promoters only show Pol II occupancy when the gene is active. The promoter types are associated with different chromatin accessibility in ATAC-seq data and have different expression characteristics in single-cell RNA-seq data. The two promoter types may therefore be optimized for different properties: paused promoters show more consistent expression when active, while TATA promoters have lower background expression when inactive. We propose that tissue-specific genes have evolved to use two different strategies for their differential expression across tissues.

Spanning the gap: unraveling RSC dynamics in vivo

Multiple reports over the past 2 years have provided the first complete structural analyses for the essential yeast chromatin remodeler, RSC, providing elaborate molecular details for its engagement with the nucleosome. However, there still remain gaps in resolution, particularly within the many RSC subunits that harbor histone binding domains.

Solving contacts at these interfaces is crucial because they are regulated by posttranslational modifications that control remodeler binding modes and function. Modifications are dynamic in nature often corresponding to transcriptional activation states and cell cycle stage, highlighting not only a need for enriched spatial resolution but also temporal understanding of remodeler engagement with the nucleosome. Our recent work sheds light on some of those gaps by exploring the binding interface between the RSC catalytic motor protein, Sth1, and the nucleosome, in the living nucleus. Using genetically encoded photo-activatable amino acids incorporated into histones of living yeast we are able to monitor the nucleosomal binding of RSC, emphasizing the regulatory roles of histone modifications in a spatiotemporal manner. We observe that RSC prefers to bind H2B SUMOylated nucleosomes in vivo and interacts with neighboring nucleosomes via H3K14ac. Additionally, we establish that RSC is constitutively bound to the nucleosome and is not ejected during mitotic chromatin compaction but alters its binding mode as it progresses through the cell cycle. Our data offer a renewed perspective on RSC mechanics under true physiological conditions.

Gene Regulation in and out of Equilibrium

Determining whether and how a gene is transcribed are two of the central processes of life. The conceptual basis for understanding such gene regulation arose from pioneering biophysical studies in eubacteria. However, eukaryotic genomes exhibit vastly greater complexity, which raises questions not addressed by this bacterial paradigm. First, how is information integrated from many widely separated binding sites to determine how a gene is transcribed? Second, does the presence of multiple energy-expending mechanisms, which are absent from eubacterial genomes, indicate that eukaryotes are capable of improved forms of genetic information processing? An updated biophysical foundation is needed to answer such questions. We describe the linear framework, a graph-based approach to Markov processes, and show that it can accommodate many previous studies in the field. Under the assumption of thermodynamic equilibrium, we introduce a language of higher-order cooperativities and show how it can rigorously quantify gene regulatory properties suggested by experiment. We point out that fundamental limits to information processing arise at thermodynamic equilibrium and can only be bypassed through energy expenditure. Finally, we outline some of the mathematical challenges that must be overcome to construct an improved biophysical understanding of gene regulation.

Identifying chromatin features that regulate gene expression distribution

Gene expression variability, differences in the number of mRNA per cell across a population of cells, is ubiquitous across diverse organisms with broad impacts on cellular phenotypes. The role of chromatin in regulating average gene expression has been extensively studied. However, what aspects of the chromatin contribute to gene expression variability is still underexplored. Here we addressed this problem by leveraging chromatin diversity and using a systematic investigation of randomly integrated expression reporters to identify what aspects of chromatin microenvironment contribute to gene expression variability. Using DNA barcoding and split-pool decoding, we created a large library of isogenic reporter clones and identified reporter integration sites in a massive and parallel manner. By mapping our measurements of reporter expression at different genomic loci with multiple epigenetic profiles including the enrichment of transcription factors and the distance to different chromatin states, we identified new factors that impact the regulation of gene expression distributions.

Microbial metabolic noise

From the time a cell was first placed under the microscope, it became apparent that identifying two clonal cells that "look" identical is extremely challenging. Since then, cell-to-cell differences in shape, size, and protein content have been carefully examined, informing us of the ultimate limits that hinder two cells from occupying an identical phenotypic state. Here, we present recent experimental and computational evidence that similar limits emerge also in cellular metabolism. These limits pertain to stochastic metabolic dynamics and, thus, cell-to-cell metabolic variability, including the resulting adapting benefits. We review these phenomena with a focus on microbial metabolism and conclude with a brief outlook on the potential relationship between metabolic noise and adaptive evolution. This article is categorized under: Metabolic Diseases > Computational Models Metabolic Diseases > Biomedical Engineering.

Control strategies for the timing of intracellular events

While the timing of intracellular events is essential for many cellular processes, gene expression inside a single cell can exhibit substantial cell-to-cell variability, raising the question of how cells ensure precision in event timing despite such stochasticity. We address this question by analyzing a biologically reasonable model of gene expression in the context of first passage time (FPT), focusing on two experimentally measurable statistics: mean FPT (MFPT) and timing variability (TV). We show that (1) transcriptional burst size (BS) and burst frequency (BF) can minimize the TV; (2) translational BS monotonically reduces the MFPT to a nonzero low bound; (3) the timescale of promoter kinetics can minimize both the MFPT and the TV, depending on the ratio of the on-switching rate over the off-switching rate; and (4) positive feedback regulation of any form can all minimize the TV, whereas negative feedback regulation of transcriptional BF or BS always enhances the TV. These control strategies can have broad implications for diverse cellular processes relying on precise temporal triggering of events.

Nucleosomal proofreading of activator-promoter interactions

Specificity in transcriptional regulation is imparted by transcriptional activators that bind to specific DNA sequences from which they stimulate transcription. Specificity may be increased by slowing down the kinetics of regulation: by increasing the energy for dissociation of the activator-DNA complex or decreasing activator concentration. In general, higher dissociation energies imply longer DNA dwell times of the activator; the activator-bound gene may not readily turn off again. Lower activator concentrations entail longer pauses between binding events; the activator-unbound gene is not easily turned on again and activated transcription occurs in stochastic bursts. We show that kinetic proofreading of activator-DNA recognition-insertion of an energy-dissipating delay step into the activation pathway for transcription-reconciles high specificity of transcriptional regulation with fast regulatory kinetics. We show that kinetic proofreading results from the stochastic removal and reformation of promoter nucleosomes, at a distance from equilibrium.

Does genome surveillance explain the global discrepancy between binding and effect of chromatin factors?

Knocking out a chromatin factor often does not alter the transcription of its binding targets. What explains the observed disconnect between binding and effect? We hypothesize that this discrepancy could be associated with the role of chromatin factors in maintaining genetic and epigenetic integrity at promoters, and not necessarily with transcription. Through re-analysis of published datasets, we present several lines of evidence that support our hypothesis and deflate the popular assumptions. We also tested the hypothesis through mutation accumulation assays on yeast knockouts of chromatin factors. Altogether, the proposed hypothesis presents a simple explanation for the global discord between chromatin factor binding and effect. Future work in this direction might fortify the hypothesis and elucidate the underlying mechanisms.

Stochastic transcription in the p53-mediated response to DNA damage is modulated by burst frequency

Discontinuous transcription has been described for different mammalian cell lines and numerous promoters. However, our knowledge of how the activity of individual promoters is adjusted by dynamic signaling inputs from transcription factors is limited. To address this question, we characterized the activity of selected target genes that are regulated by pulsatile accumulation of the tumor suppressor p53 in response to ionizing radiation. We performed time-resolved measurements of gene expression at the single-cell level by smFISH and used the resulting data to inform a mathematical model of promoter activity. We found that p53 target promoters are regulated by frequency modulation of stochastic bursting and can be grouped along three archetypes of gene expression. The occurrence of these archetypes cannot solely be explained by nuclear p53 abundance or promoter binding of total p53. Instead, we provide evidence that the time-varying acetylation state of p53's C-terminal lysine residues is critical for gene-specific regulation of stochastic bursting.

Transcriptional processes: Models and inference

Many biochemical events involve multistep reactions. One of the most important biological processes that involve multistep reaction is the transcriptional process. Models for multistep reaction necessarily need multiple states and it is a challenge to compute model parameters that best agree with experimental data. Therefore, the aim of this work is to design a multistep promoter model which accurately characterizes transcriptional bursting and is consistent with observed data. To address this issue, we develop a model for promoters with several OFF states and a single ON state using Erlang distribution. To explore the combined effects of model and data, we combine Monte Carlo extension of Expectation Maximization (MCEM) and delay Stochastic Simulation Algorithm (DSSA) and call the resultant algorithm as delay Bursty MCEM. We apply this algorithm to time-series data of endogenous mouse glutaminase promoter to validate the model assumptions and infer the kinetic parameters. Our results show that with multiple OFF states, we are able to infer and produce a model which is more consistent with experimental data. Our results also show that delay Bursty MCEM inference is more efficient.

Fold-Change Detection of NF-κB at Target Genes with Different Transcript Outputs

The transcription factor nuclear factor (NF)-κB promotes inflammatory and stress-responsive gene transcription across a range of cell types in response to the cytokine tumor necrosis factor (TNF). Although NF-κB signaling exhibits significant variability across single cells, some target genes supporting high levels of TNF-inducible transcription exhibit fold-change detection of NF-κB, which may buffer against stochastic variation in signaling molecules. It is unknown whether fold-change detection is maintained at NF-κB target genes with low levels of TNF-inducible transcription, for which stochastic promoter events may be more pronounced. Here, we used a microfluidic cell-trapping device to measure how TNF-induced activation of NF-κB controls transcription in single Jurkat T cells at the promoters of integrated HIV and the endogenous cytokine gene IL6, which produce only a few transcripts per cell. We tracked TNF-stimulated NF-κB RelA nuclear translocation by live-cell imaging and then quantified transcript number by RNA FISH in the same cell. We found that TNF-induced transcript abundance at 2 h for low- and high-abundance target genes correlates with similar strength with the fold change in nuclear NF-κB. A computational model of TNF-NF-κB signaling, which implements fold-change detection from competition for binding to κB motifs, could reproduce fold-change detection across the experimentally measured range of transcript outputs. However, multiple model parameters affecting transcription had to be simultaneously varied across promoters to maintain fold-change detection while also matching other trends in the single-cell data for low-abundance transcripts. Our results suggest that cells use multiple biological mechanisms to tune transcriptional output while maintaining robustness of NF-κB fold-change detection.

Modulation of first-passage time for bursty gene expression via random signals

The stochastic nature of cell-specific signal molecules (such as transcription factor, ribosome, etc.) and the intrinsic stochastic nature of gene expression process result in cell-to-cell variations at protein levels. Increasing experimental evidences suggest that cell phenotypic variations often depend on the accumulation of some special proteins. Hence, a natural and fundamental question is: How does input signal affect the timing of protein count up to a given threshold? To this end, we study effects of input signal on the first-passage time (FPT), the time at which the number of proteins crosses a given threshold. Input signal is distinguished into two types: constant input signal and random input signal, regulating only burst frequency (or burst size) of gene expression. Firstly, we derive analytical formulae for FPT moments in each case of constant signal regulation and random signal regulation. Then, we find that random input signal tends to increases the mean and noise of FPT compared with constant input signal. Finally, we observe that different regulation ways of random signal have different effects on FPT, that is, burst size modulation tends to decrease the mean of FPT and increase the noise of FPT compared with burst frequency modulation. Our findings imply a fundamental mechanism that random fluctuating environment may prolong FPT. This can provide theoretical guidance for studies of some cellular key events such as latency of HIV and lysis time of bacteriophage λ. In conclusion, our results reveal impacts of external signal on FPT and aid understanding the regulation mechanism of gene expression.

Fine-tuning of noise in gene expression with nucleosome remodeling

Engineering stochastic fluctuations of gene expression (or “noise”) is integral to precisely bias cellular-fate decisions and statistical phenotypes in both single-cell and multi-cellular systems. Epigenetic regulation has been shown to constitute a large source of noise, and thus, engineering stochasticity is deeply intertwined with epigenetics. Here, utilizing chromatin remodeling, we report that Caffeic acid phenethyl ester (CA) and Pyrimethamine (PYR), two inhibitors of BAF250a, a subunit of the Brahma-associated factor (BAF) nucleosome remodeling complex, enable differential and tunable control of noise in transcription and translation from the human immunodeficiency virus long terminal repeat promoter in a dose and time-dependent manner. CA conserves noise levels while increasing mean abundance, resulting in direct tuning of the transcriptional burst size, while PYR strictly increases transcriptional initiation frequency while conserving a constant transcriptional burst size. Time-dependent treatment with CA reveals non-continuous tuning with noise oscillating at a constant mean abundance at early time points and the burst size increasing for treatments after 5 h. Treatments combining CA and Protein Kinase C agonists result in an even larger increase of abundance while conserving noise levels with a highly non-linear increase in variance of up to 63× untreated controls. Finally, drug combinations provide non-antagonistic combinatorial tuning of gene expression noise and map a noise phase space for future applications with viral and synthetic gene vectors. Active remodeling of nucleosomes and BAF-mediated control of gene expression noise expand a toolbox for the future design and engineering of stochasticity in living systems.

Simulating biological processes: stochastic physics from whole cells to colonies

The last few decades have revealed the living cell to be a crowded spatially heterogeneous space teeming with biomolecules whose concentrations and activities are governed by intrinsically random forces. It is from this randomness, however, that a vast array of precisely timed and intricately coordinated biological functions emerge that give rise to the complex forms and behaviors we see in the biosphere around us. This seemingly paradoxical nature of life has drawn the interest of an increasing number of physicists, and recent years have seen stochastic modeling grow into a major subdiscipline within biological physics. Here we review some of the major advances that have shaped our understanding of stochasticity in biology. We begin with some historical context, outlining a string of important experimental results that motivated the development of stochastic modeling. We then embark upon a fairly rigorous treatment of the simulation methods that are currently available for the treatment of stochastic biological models, with an eye toward comparing and contrasting their realms of applicability, and the care that must be taken when parameterizing them. Following that, we describe how stochasticity impacts several key biological functions, including transcription, translation, ribosome biogenesis, chromosome replication, and metabolism, before considering how the functions may be coupled into a comprehensive model of a 'minimal cell'. Finally, we close with our expectation for the future of the field, focusing on how mesoscopic stochastic methods may be augmented with atomic-scale molecular modeling approaches in order to understand life across a range of length and time scales.

Frequency Modulation of Transcriptional Bursting Enables Sensitive and Rapid Gene Regulation

Gene regulation is a complex non-equilibrium process. Here, we show that quantitating the temporal regulation of key gene states (transcriptionally inactive, active, and refractory) provides a parsimonious framework for analyzing gene regulation. Our theory makes two non-intuitive predictions. First, for transcription factors (TFs) that regulate transcription burst frequency, as opposed to amplitude or duration, weak TF binding is sufficient to elicit strong transcriptional responses. Second, refractoriness of a gene after a transcription burst enables rapid responses to stimuli. We validate both predictions experimentally by exploiting the natural, optogenetic-like responsiveness of the Neurospora GATA-type TF White Collar Complex (WCC) to blue light. Further, we demonstrate that differential regulation of WCC target genes is caused by different gene activation rates, not different TF occupancy, and that these rates are tuned by both the core promoter and the distance between TF-binding site and core promoter. In total, our work demonstrates the relevance of a kinetic, non-equilibrium framework for understanding transcriptional regulation.

An Investigation on the Expression Level of Interleukin 10 (IL-10) in the Healthy and Mastitic Holstein Cows and the Bioinformatics Analysis of Nucleosome Profile

Cytokines are immune regulators that play an essential role in regulating immune response against various infections. The present study focused on the possible association between the expression level of Interleukin 10 (IL-10) in blood and milk samples of 25 healthy and 25 mastitic cows in Fars province, Iran, using a quantitative real-time PCR assay. The experimental groups were categorized according to the number of calvings. The expression level of IL-10 was significantly higher in the blood and milk samples of mastitic cows compared to the healthy ones. Concomitant to increasing the number of calving, a numerical elevation in the expression of IL-10 in blood was observed (P < 0.05). The bioinformatics analysis of IL-10 gene revealed the promoter, exon-intron regions, and nucleosome profile. The nucleosome occupancy site was finally predicted using NUPOP software. Our result indicated that the promoter was not exactly placed in the nucleosome region, which was finally aimed to predict the position and expression of IL-10 gene in the mastitic cows.

Mechanisms of action and regulation of ATP-dependent chromatin-remodelling complexes

Cells utilize diverse ATP-dependent nucleosome-remodelling complexes to carry out histone sliding, ejection or the incorporation of histone variants, suggesting that different mechanisms of action are used by the various chromatin-remodelling complex subfamilies. However, all chromatin-remodelling complex subfamilies contain an ATPase-translocase 'motor' that translocates DNA from a common location within the nucleosome. In this Review, we discuss (and illustrate with animations) an alternative, unifying mechanism of chromatin remodelling, which is based on the regulation of DNA translocation. We propose the 'hourglass' model of remodeller function, in which each remodeller subfamily utilizes diverse specialized proteins and protein domains to assist in nucleosome targeting or to differentially detect nucleosome epitopes. These modules converge to regulate a common DNA translocation mechanism, to inform the conserved ATPase 'motor' on whether and how to apply DNA translocation, which together achieve the various outcomes of chromatin remodelling: nucleosome assembly, chromatin access and nucleosome editing.

Real-time observation of light-controlled transcription in living cells

Gene expression is tightly regulated in space and time. To dissect this process with high temporal resolution, we introduce an optogenetic tool termed BLInCR (Blue Light-Induced Chromatin Recruitment) that combines rapid and reversible light-dependent recruitment of effector proteins with a real-time readout for transcription. We used BLInCR to control the activity of a reporter gene cluster in the human osteosarcoma cell line U2OS by reversibly recruiting the viral transactivator VP16. RNA production was detectable ∼2 minutes after VP16 recruitment and readily decreased when VP16 dissociated from the cluster in the absence of light. Quantitative assessment of the activation process revealed biphasic activation kinetics with a pronounced early phase in cells treated with the histone deacetylase inhibitor SAHA. Comparison with kinetic models for transcription activation suggests that the gene cluster undergoes a maturation process when activated. We anticipate that BLInCR will facilitate the study of transcription dynamics in living cells.

Theoretical estimates of exposure timescales of protein binding sites on DNA regulated by nucleosome kinetics

It is being increasingly realized that nucleosome organization on DNA crucially regulates DNA-protein interactions and the resulting gene expression. While the spatial character of the nucleosome positioning on DNA has been experimentally and theoretically studied extensively, the temporal character is poorly understood. Accounting for ATPase activity and DNA-sequence effects on nucleosome kinetics, we develop a theoretical method to estimate the time of continuous exposure of binding sites of non-histone proteins (e.g. transcription factors and TATA binding proteins) along any genome. Applying the method to Saccharomyces cerevisiae, we show that the exposure timescales are determined by cooperative dynamics of multiple nucleosomes, and their behavior is often different from expectations based on static nucleosome occupancy. Examining exposure times in the promoters of GAL1 and PHO5, we show that our theoretical predictions are consistent with known experiments. We apply our method genome-wide and discover huge gene-to-gene variability of mean exposure times of TATA boxes and patches adjacent to TSS (+1 nucleosome region); the resulting timescale distributions have non-exponential tails.

Role of transcription factor-mediated nucleosome disassembly in PHO5 gene expression

Studying nucleosome dynamics in promoter regions is crucial for understanding gene regulation. Nucleosomes regulate gene expression by sterically occluding transcription factors (TFs) and other non–histone proteins accessing genomic DNA. How the binding competition between nucleosomes and TFs leads to transcriptionally compatible promoter states is an open question. Here, we present a computational study of the nucleosome dynamics and organization in the promoter region of PHO5 gene in Saccharomyces cerevisiae. Introducing a model for nucleosome kinetics that takes into account ATP-dependent remodeling activity, DNA sequence effects, and kinetics of TFs (Pho4p), we compute the probability of obtaining different “promoter states” having different nucleosome configurations. Comparing our results with experimental data, we argue that the presence of local remodeling activity (LRA) as opposed to basal remodeling activity (BRA) is crucial in determining transcriptionally active promoter states. By modulating the LRA and Pho4p binding rate, we obtain different mRNA distributions—Poisson, bimodal, and long-tail. Through this work we explain many features of the PHO5 promoter such as sequence-dependent TF accessibility and the role of correlated dynamics between nucleosomes and TFs in opening/coverage of the TATA box. We also obtain possible ranges for TF binding rates and the magnitude of LRA.

Effects of DNA replication on mRNA noise

There are several sources of fluctuations in gene expression. Here we study the effects of time-dependent DNA replication, itself a tightly controlled process, on noise in mRNA levels. Stochastic simulations of constitutive and regulated gene expression are used to analyze the time-averaged mean and variation in each case. The simulations demonstrate that to capture mRNA distributions correctly, chromosome replication must be realistically modeled. Slow relaxation of mRNA from the low copy number steady state before gene replication to the high steady state after replication is set by the transcript's half-life and contributes significantly to the shape of the mRNA distribution. Consequently both the intrinsic kinetics and the gene location play an important role in accounting for the mRNA average and variance. Exact analytic expressions for moments of the mRNA distributions that depend on the DNA copy number, gene location, cell doubling time, and the rates of transcription and degradation are derived for the case of constitutive expression and subsequently extended to provide approximate corrections for regulated expression and RNA polymerase variability. Comparisons of the simulated models and analytical expressions to experimentally measured mRNA distributions show that they better capture the physics of the system than previous theories.

From Structural Variation of Gene Molecules to Chromatin Dynamics and Transcriptional Bursting

Transcriptional activation of eukaryotic genes is accompanied, in general, by a change in the sensitivity of promoter chromatin to endonucleases. The structural basis of this alteration has remained elusive for decades; but the change has been viewed as a transformation of one structure into another, from "closed" to "open" chromatin. In contradistinction to this static and deterministic view of the problem, a dynamical and probabilistic theory of promoter chromatin has emerged as its solution. This theory, which we review here, explains observed variation in promoter chromatin structure at the level of single gene molecules and provides a molecular basis for random bursting in transcription-the conjecture that promoters stochastically transition between transcriptionally conducive and inconducive states. The mechanism of transcriptional regulation may be understood only in probabilistic terms.

Nucleosomes, transcription, and probability

Speaking of current measurements on single ion channel molecules, David Colquhoun wrote in 2006, "Individual molecules behave randomly, so suddenly we had to learn how to deal with stochastic processes." Here I describe theoretical efforts to understand recent experimental observations on the chromatin structure of single gene molecules, a molecular biologist's path toward probabilistic theories.

Chromatin structure analysis of single gene molecules by psoralen cross-linking and electron microscopy

Nucleosomes occupy a central role in regulating eukaryotic gene expression by blocking access of transcription factors to their target sites on chromosomal DNA. Analysis of chromatin structure and function has mostly been performed by probing DNA accessibility with endonucleases. Such experiments average over large numbers of molecules of the same gene, and more recently, over entire genomes. However, both digestion and averaging erase the structural variation between molecules indicative of dynamic behavior, which must be reconstructed for any theory of regulation. Solution of this problem requires the structural analysis of single gene molecules. In this chapter, we describe a method by which single gene molecules are purified from the yeast Saccharomyces cerevisiae and cross-linked with psoralen, allowing the determination of nucleosome configurations by transmission electron microscopy. We also provide custom analysis software that semi-automates the analysis of micrograph data. This single-gene technique enables detailed examination of chromatin structure at any genomic locus in yeast.

Dynamic regulation of transcription factors by nucleosome remodeling

The chromatin landscape and promoter architecture are dominated by the interplay of nucleosome and transcription factor (TF) binding to crucial DNA sequence elements. However, it remains unclear whether nucleosomes mobilized by chromatin remodelers can influence TFs that are already present on the DNA template. In this study, we investigated the interplay between nucleosome remodeling, by either yeast ISW1a or SWI/SNF, and a bound TF. We found that a TF serves as a major barrier to ISW1a remodeling, and acts as a boundary for nucleosome repositioning. In contrast, SWI/SNF was able to slide a nucleosome past a TF, with concurrent eviction of the TF from the DNA, and the TF did not significantly impact the nucleosome positioning. Our results provide direct evidence for a novel mechanism for both nucleosome positioning regulation by bound TFs and TF regulation via dynamic repositioning of nucleosomes.

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

Cells contain thousands of genes that are encoded by molecules of DNA. In yeast and other eukaryotic organisms, this DNA is wrapped around proteins called histones to make structures called nucleosomes. This compacts the DNA and allows it to fit inside the tiny nucleus within the cell. The positioning of the nucleosomes influences how tightly packed the DNA is, which in turn influences the activity of genes. Less active genes tend to be found within regions of DNA that are tightly packed, while more active genes are found in less tightly packed regions.

To activate a gene, proteins called transcription factors bind to a section of DNA within the gene called the promoter. Enzymes known as ‘chromatin remodelers’ can alter the locations of nucleosomes on DNA to allow the transcription factors access to the promoters of particular genes. In yeast, the SWI/SNF family of chromatin remodelers can disassemble nucleosomes to promote gene activity, while the ISW1 family organises nucleosomes into closely spaced groups to repress gene activity. However, it is not clear if, or how, chromatin remodelers can influence transcription factors that are already bound to DNA.

Here, Li et al. studied the interactions between a transcription factor and the chromatin remodelers in yeast. The experiment used a piece of DNA that contained a bound transcription factor and a single nucleosome. Li et al. used a technique called ‘single molecule DNA unzipping’, which enabled them to precisely locate the position of the nucleosome and transcription factor before and after the nucleosome was remodeled. The experiments found that a chromatin remodeler called ISW1a moved the nucleosome away from the transcription factor, while a SWI/SNF chromatin remodeler moved the nucleosome towards it.

Significantly, Li et al. also found that a transcription factor is a major barrier to ISW1a's remodeling activity, suggesting that ISW1a may use transcription factors as reference points to position nucleosomes. In contrast, SWI/SNF was able to slide a nucleosome past the transcription factor, which led to the transcription factor falling off the DNA. Therefore, SWI/SNF is able to move transcription factors out of the way to deactivate genes.

Li et al. propose a new model for how chromatin remodelers can move nucleosomes and regulate transcription factors to alter gene activity. A future challenge will be to observe these types of activities in living cells.

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

Effects of promoter leakage on dynamics of gene expression

Background

Quantitative analysis of simple molecular networks is an important step forward understanding fundamental intracellular processes. As network motifs occurring recurrently in complex biological networks, gene auto-regulatory circuits have been extensively studied but gene expression dynamics remain to be fully understood, e.g., how promoter leakage affects expression noise is unclear.

Results

In this work, we analyze a gene model with auto regulation, where the promoter is assumed to have one active state with highly efficient transcription and one inactive state with very lowly efficient transcription (termed as promoter leakage). We first derive the analytical distribution of gene product, and then analyze effects of promoter leakage on expression dynamics including bursting kinetics. Interestingly, we find that promoter leakage always reduces expression noise and that increasing the leakage rate tends to simplify phenotypes. In addition, higher leakage results in fewer bursts.

Conclusions

Our results reveal the essential role of promoter leakage in controlling expression dynamics and further phenotype. Specifically, promoter leakage is a universal mechanism of reducing expression noise, controlling phenotypes in different environments and making the gene produce generate fewer bursts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0157-z) contains supplementary material, which is available to authorized users.

Dynamical analysis of mCAT2 gene models with CTN-RNA nuclear retention

As an experimentally well-studied nuclear-retained RNA, CTN-RNA plays a significant role in many aspects of mouse cationic amino acid transporter 2 (mCAT2) gene expression, but relevant dynamical mechanisms have not been completely clarified. Here we first show that CTN-RNA nuclear retention can not only reduce pre-mCAT2 RNA noise but also mediate its coding partner noise. Then, by collecting experimental observations, we conjecture a heterodimer formed by two proteins, p54(nrb) and PSP1, named p54(nrb)-PSP1, by which CTN-RNA can positively regulate the expression of nuclear mCAT2 RNA. Therefore, we construct a sequestration model at the molecular level. By analyzing the dynamics of this model system, we demonstrate why most nuclear-retained CTN-RNAs stabilize at the periphery of paraspeckles, how CTN-RNA regulates its protein-coding partner, and how the mCAT2 gene can maintain a stable expression. In particular, we obtain results that can easily explain the experimental phenomena observed in two cases, namely, when cells are stressed and unstressed. Our entire analysis not only reveals that CTN-RNA nuclear retention may play an essential role in indirectly preventing diseases but also lays the foundation for further study of other members of the nuclear-regulatory RNA family with more complicated molecular mechanisms.

Histone chaperones: assisting histone traffic and nucleosome dynamics

The functional organization of eukaryotic DNA into chromatin uses histones as components of its building block, the nucleosome. Histone chaperones, which are proteins that escort histones throughout their cellular life, are key actors in all facets of histone metabolism; they regulate the supply and dynamics of histones at chromatin for its assembly and disassembly. Histone chaperones can also participate in the distribution of histone variants, thereby defining distinct chromatin landscapes of importance for genome function, stability, and cell identity. Here, we discuss our current knowledge of the known histone chaperones and their histone partners, focusing on histone H3 and its variants. We then place them into an escort network that distributes these histones in various deposition pathways. Through their distinct interfaces, we show how they affect dynamics during DNA replication, DNA damage, and transcription, and how they maintain genome integrity. Finally, we discuss the importance of histone chaperones during development and describe how misregulation of the histone flow can link to disease.

Essential gene disruptions reveal complex relationships between phenotypic robustness, pleiotropy, and fitness

The concept of robustness in biology has gained much attention recently, but a mechanistic understanding of how genetic networks regulate phenotypic variation has remained elusive. One approach to understand the genetic architecture of variability has been to analyze dispensable gene deletions in model organisms; however, the most important genes cannot be deleted. Here, we have utilized two systems in yeast whereby essential genes have been altered to reduce expression. Using high-throughput microscopy and image analysis, we have characterized a large number of morphological phenotypes, and their associated variation, for the majority of essential genes in yeast. Our results indicate that phenotypic robustness is more highly dependent upon the expression of essential genes than on the presence of dispensable genes. Morphological robustness appears to be a general property of a genotype that is closely related to pleiotropy. While the fitness profile across a range of expression levels is idiosyncratic to each gene, the global pattern indicates that there is a window in which phenotypic variation can be released before fitness effects are observable.

Reinitiation enhances reliable transcriptional responses in eukaryotes

Gene transcription is a noisy process carried out by the transcription machinery recruited to the promoter. Noise reduction is a fundamental requirement for reliable transcriptional responses which in turn are crucial for signal transduction. Compared with the relatively simple transcription initiation in prokaryotes, eukaryotic transcription is more complex partially owing to its additional reinitiation mechanism. By theoretical analysis, we showed that reinitiation reduces noise in eukaryotic transcription independent of the transcription level. Besides, a higher reinitiation rate enables a stable scaffold complex an advantage in noise reduction. Finally, we showed that the coupling between scaffold formation and transcription can further reduce transcription noise independent of the transcription level. Furthermore, compared with the reinitiation mechanism, the noise reduction effect of the coupling can be of more significance in the case that the transcription level is low and the intrinsic noise dominates. Our results uncover a mechanistic route which eukaryotes may use to facilitate a more reliable response in the noisy transcription process.

The Yeast Nucleosome Atlas (YNA) database: an integrative gene mining platform for studying chromatin structure and its regulation in yeast

Background

Histone modification and remodeling play crucial roles in regulating gene transcription. These post-translational modifications of histones function in a combinatorial fashion and can be recognized by specific histone-binding proteins, thus regulating gene transcription. Therefore, understanding the combinatorial patterns of the histone code is vital to understanding the associated biological processes. However, most of the datasets regarding histone modification and chromatin regulation are scattered across various studies, and no comprehensive search and query tool has yet been made available to retrieve genes bearing specific histone modification patterns and regulatory proteins.

Description

For this reason, we developed the Yeast Nucleosome Atlas database, or the YNA database, which integrates the available experimental data on nucleosome occupancy, histone modifications, the binding occupancy of regulatory proteins, and gene expression data, and provides the genome-wide gene miner to retrieve genes with a specific combination of these chromatin-related datasets. Moreover, the biological significance analyzer, which analyzes the enrichments of histone modifications, binding occupancy, transcription rate, and functionality of the retrieved genes, was constructed to help researchers to gain insight into the correlation among chromatin regulation and transcription.

Conclusions

Compared to previously established genome browsing databases, YNA provides a powerful gene mining and retrieval interface, and is an investigation tool that can assist users to generate testable hypotheses for studying chromatin regulation during transcription. YNA is available online at http://cosbi3.ee.ncku.edu.tw/yna/.

Nucleosomal promoter variation generates gene expression noise

Gene product molecule numbers fluctuate over time and between cells, confounding deterministic expectations. The molecular origins of this noise of gene expression remain unknown. Recent EM analysis of single PHO5 gene molecules of yeast indicated that promoter molecules stochastically assume alternative nucleosome configurations at steady state, including the fully nucleosomal and nucleosome-free configuration. Given that distinct configurations are unequally conducive to transcription, the nucleosomal variation of promoter molecules may constitute a source of gene expression noise. This notion, however, implies an untested conjecture, namely that the nucleosomal variation arises de novo or intrinsically (i.e., that it cannot be explained as the result of the promoter's deterministic response to variation in its molecular surroundings). Here, we show--by microscopically analyzing the nucleosome configurations of two juxtaposed physically linked PHO5 promoter copies--that the configurational variation, indeed, is intrinsically stochastic and thus, a cause of gene expression noise rather than its effect.

Promoter-mediated transcriptional dynamics

Genes in eukaryotic cells are typically regulated by complex promoters containing multiple binding sites for a variety of transcription factors, but how promoter dynamics affect transcriptional dynamics has remained poorly understood. In this study, we analyze gene models at the transcriptional regulation level, which incorporate the complexity of promoter structure (PS) defined as transcriptional exits (i.e., ON states of the promoter) and the transition pattern (described by a matrix consisting of transition rates among promoter activity states). We show that multiple exits of transcription are the essential origin of generating multimodal distributions of mRNA, but promoters with the same transition pattern can lead to multimodality of different modes, depending on the regulation of transcriptional factors. In turn, for similar mRNA distributions in the models, the mean ON or OFF time distributions may exhibit different characteristics, thus providing the supplemental information on PS. In addition, we demonstrate that the transcriptional noise can be characterized by a nonlinear function of mean ON and OFF times. These results not only reveal essential characteristics of promoter-mediated transcriptional dynamics but also provide signatures useful for inferring PS based on characteristics of transcriptional outputs.

The yeast PHO5 promoter: from single locus to systems biology of a paradigm for gene regulation through chromatin

Chromatin dynamics crucially contributes to gene regulation. Studies of the yeast PHO5 promoter were key to establish this nowadays accepted view and continuously provide mechanistic insight in chromatin remodeling and promoter regulation, both on single locus as well as on systems level. The PHO5 promoter is a context independent chromatin switch module where in the repressed state positioned nucleosomes occlude transcription factor sites such that nucleosome remodeling is a prerequisite for and not consequence of induced gene transcription. This massive chromatin transition from positioned nucleosomes to an extensive hypersensitive site, together with respective transitions at the co-regulated PHO8 and PHO84 promoters, became a prime model for dissecting how remodelers, histone modifiers and chaperones co-operate in nucleosome remodeling upon gene induction. This revealed a surprisingly complex cofactor network at the PHO5 promoter, including five remodeler ATPases (SWI/SNF, RSC, INO80, Isw1, Chd1), and demonstrated for the first time histone eviction in trans as remodeling mode in vivo. Recently, the PHO5 promoter and the whole PHO regulon were harnessed for quantitative analyses and computational modeling of remodeling, transcription factor binding and promoter input-output relations such that this rewarding single-locus model becomes a paradigm also for theoretical and systems approaches to gene regulatory networks.

Threonine-4 of the budding yeast RNAP II CTD couples transcription with Htz1-mediated chromatin remodeling

The C-terminal domain (CTD) of the largest subunit of RNA polymerase II (RNAP II) consists of repeated YSPTSPS heptapeptides and connects transcription with cotranscriptional events. Threonine-4 (Thr4) of the CTD repeats has been shown to function in histone mRNA 3'-end processing in chicken cells and in transcriptional elongation in human cells. Here, we demonstrate that, in budding yeast, Thr4, although dispensable for growth in rich media, is essential in phosphate-depleted or galactose-containing media. Thr4 is required to maintain repression of phosphate-regulated (PHO) genes under normal growth conditions and for full induction of PHO5 and the galactose-induced GAL1 and GAL7 genes. We identify genetic links between Thr4 and the histone variant Htz1 and show that Thr4, as well as the Ino80 chromatin remodeler, is required for activation-associated eviction of Htz1 specifically from promoters of the Thr4-dependent genes. Our study uncovers a connection between transcription and chromatin remodeling linked by Thr4 of the CTD.

Dependence of the sperm/oocyte decision on the nucleosome remodeling factor complex was acquired during recent Caenorhabditis briggsae evolution

The major families of chromatin remodelers have been conserved throughout eukaryotic evolution. Because they play broad, pleiotropic roles in gene regulation, it was not known if their functions could change rapidly. Here, we show that major alterations in the use of chromatin remodelers are possible, because the nucleosome remodeling factor (NURF) complex has acquired a unique role in the sperm/oocyte decision of the nematode Caenorhabditis briggsae. First, lowering the activity of C. briggsae NURF-1 or ISW-1, the core components of the NURF complex, causes germ cells to become oocytes rather than sperm. This observation is based on the analysis of weak alleles and null mutations that were induced with TALENs and on RNA interference. Second, qRT-polymerase chain reaction data show that the C. briggsae NURF complex promotes the expression of Cbr-fog-1 and Cbr-fog-3, two genes that control the sperm/oocyte decision. This regulation occurs in the third larval stage and affects the expression of later spermatogenesis genes. Third, double mutants reveal that the NURF complex and the transcription factor TRA-1 act independently on Cbr-fog-1 and Cbr-fog-3. TRA-1 binds both promoters, and computer analyses predict that these binding sites are buried in nucleosomes, so we suggest that the NURF complex alters chromatin structure to allow TRA-1 access to Cbr-fog-1 and Cbr-fog-3. Finally, lowering NURF activity by mutation or RNA interference does not affect this trait in other nematodes, including the sister species C. nigoni, so it must have evolved recently. We conclude that altered chromatin remodeling could play an important role in evolutionary change.

Single-cell nucleosome mapping reveals the molecular basis of gene expression heterogeneity

Nucleosomes, the basic unit of chromatin, have a critical role in the control of gene expression. Nucleosome positions have generally been determined by examining bulk populations of cells and then correlated with overall gene expression. Here, we describe a technique to determine nucleosome positioning in single cells by virtue of the ability of the nucleosome to protect DNA from GpC methylation. In the acid phosphatase inducible PHO5 gene, we find that there is significant cell-to-cell variation in nucleosome positions and shifts in nucleosome positioning correlate with changes in gene expression. However, nucleosome positioning is not absolute, and even with major shifts in gene expression, some cells fail to change nucleosome configuration. Mutations of the PHO5 promoter that introduce a poly(dA:dT) tract-stimulated gene expression under nonpermissive conditions led to shifts of positioned nucleosomes similar to induction of PHO5. By contrast, mutations that altered AA/TT/AT periodicity reduced gene expression upon PHO5 induction and stabilized nucleosomes in most cells, suggesting that enhanced nucleosome affinity for DNA antagonizes chromatin remodelers. Finally, we determined nucleosome positioning in two regions described as "fuzzy" or nucleosome-free when examined in a bulk assay. These regions consisted of distinct nucleosomes with a larger footprint for potential location and an increase population of cells lacking a nucleosome altogether. These data indicate an underlying complexity of nucleosome positioning that may contribute to the flexibility and heterogeneity of gene expression.