SUMMARY Neuronal activity-inducible gene transcription correlates with rapid and transient increases in histone acetylation at promoters and enhancers of activity-regulated genes. Exactly how histone acetylation modulates transcription of these genes has remained unknown. We used single-cell in situ transcriptional analysis to show that Fos and Npas4 are transcribed in stochastic bursts in mouse neurons and that membrane depolarization increases mRNA expression by increasing burst frequency. We then expressed dCas9-p300 or dCas9-HDAC8 fusion proteins to mimic or block activity-induced histone acetylation locally at enhancers. Adding histone acetylation increased Fos transcription by prolonging burst duration and resulted in higher Fos protein levels and an elevation of resting membrane potential. Inhibiting histone acetylation reduced Fos transcription by reducing burst frequency and impaired experience-dependent Fos protein induction in the hippocampus in vivo. Thus, activity-inducible histone acetylation tunes the transcriptional dynamics of experience-regulated genes to affect selective changes in neuronal gene expression and cellular function.

No abstract

Highlights d Integral control allows perfect adaptation but can be slow or even unstable d Molecular proportional and derivative control modules improve adaptation performance d Control modules can be composed in a modular fashion to achieve design specifications

Single-molecule RNA fluorescence in situ hybridization (smFISH) provides unparalleled resolution in the measurement of the abundance and localization of nascent and mature RNA transcripts in fixed, single cells. We developed a computational pipeline (BayFish) to infer the kinetic parameters of gene expression from smFISH data at multiple time points after gene induction. Given an underlying model of gene expression, BayFish uses a Monte Carlo method to estimate the Bayesian posterior probability of the model parameters and quantify the parameter uncertainty given the observed smFISH data. We tested BayFish on synthetic data and smFISH measurements of the neuronal activity-inducible gene Npas4 in primary neurons.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-017-1297-9) contains supplementary material, which is available to authorized users.

The ability of cells to regulate their function through feedback control is a fundamental underpinning of life. The capability to engineer de novo feedback control with biological molecules is ushering in an era of robust functionality for many applications in biotechnology and medicine. To fulfill their potential, feedback control strategies implemented with biological molecules need to be generalizable, modular and operationally predictable. Proportional-Integral-Derivative (PID) control fulfills this role for technological systems and is a commonly used strategy in engineering. Integral feedback control allows a system to return to an invariant steady-state value after step disturbances, hence enabling its robust operation. Proportional and derivative feedback control used with integral control help sculpt the dynamics of the return to steady-state following perturbation. Recently, a biomolecular implementation of integral control was proposed based on an antithetic motif in which two molecules interact stoichiometrically to annihilate each other's function. In this work, we report how proportional and derivative implementations can be layered on top of this integral architecture to achieve a biochemical PID control design. We illustrate through computational and analytical treatments that the addition of proportional and derivative control improves performance, and discuss practical biomolecular implementations of these control strategies.

Epigenetic switches are bistable, molecular systems built from self-reinforcing feedback loops that can spontaneously switch between heritable phenotypes in the absence of DNA mutation. It has been hypothesized that epigenetic switches first evolved as a mechanism of bet-hedging and adaptation, but the evolutionary trajectories and conditions by which an epigenetic switch can outcompete adaptation through genetic mutation remain unknown. Here, we used computer simulations to evolve a mechanistic, biophysical model of a self-activating genetic circuit, which can both adapt genetically through mutation and exhibit epigenetic switching. We evolved these genetic circuits under a fluctuating environment that alternatively selected for low and high protein expression levels. In all tested conditions, the population first evolved by genetic mutation towards a region of genotypes where genetic adaptation can occur faster after each environmental transition. Once in this region, the self-activating genetic circuit can exhibit epigenetic switching, which starts competing with genetic adaptation. We show a trade-off between either minimizing the adaptation time or increasing the robustness of the phenotype to biochemical noise. Epigenetic switching was superior in a fast fluctuating environment because it adapted faster than genetic mutation after an environmental transition, while still attenuating the effect of biochemical noise on the phenotype. Conversely, genetic adaptation was favored in a slowly fluctuating environment because it maximized the phenotypic robustness to biochemical noise during the constant environment between transitions, even if this resulted in slower adaptation. This simple trade-off predicts the conditions and trajectories under which an epigenetic switch evolved to outcompete genetic adaptation, shedding light on possible mechanisms by which bet-hedging strategies might emerge and persist in natural populations.

Mathematical models can aid the design of genetic circuits, but may yield inaccurate results if individual parts are not modeled at the appropriate resolution. To illustrate the importance of this concept, we study transcriptional cascades consisting of two inducible synthetic transcription factors connected in series. Despite the simplicity of this design, we find that accurate prediction of circuit behavior requires mapping the dose responses of each circuit component along the dimensions of both its expression level and its inducer concentration. Using this multidimensional characterization, we were able to computationally explore the behavior of 16 different circuit designs. We experimentally verified a subset of these predictions and found substantial agreement. This method of biological part characterization enables the use of models to identify (un)­desired circuit behaviors prior to experimental implementation, thus shortening the design–build–test cycle for more complex circuits.

Single-molecule RNA fluorescence in situ hybridization (smFISH) provides unparalleled resolution on the abundance and localization of nascent and mature transcripts in single cells. Gene expression dynamics are typically inferred by measuring mRNA abundance in small numbers of fixed cells sampled from a population at multiple time-points after induction. The sparse data that arise from the small number of cells obtained using smFISH present a challenge for inferring transcription dynamics. Here, we developed a computational pipeline (BayFish) to infer kinetic parameters of gene expression from smFISH data at multiple time points after induction. Given an underlying model of gene expression, BayFish uses a Monte Carlo method to estimate the Bayesian posterior probability of the model parameters and quantify the parameter uncertainty given the observed smFISH data. We tested BayFish on smFISH measurements of the neuronal activity inducible gene Npas4 in primary neurons. We showed that a 2-state promoter model can recapitulate Npas4 dynamics after induction and we inferred that the transition rate from the promoter OFF state to the ON state is increased by the stimulus. Author SummaryGene expression can exhibit cell-to-cell variability due to the stochastic nature of biochemical reactions. Single cell assays (e.g. smFISH) directly quantify stochastic gene expression by measuring the number of active promoters and transcripts per cell in a population of cells. The data are distributions and their shape and time-evolution contain critical information on the underlying process of gene expression. Recent work has combined models of stochastic gene expression with maximum likelihood methods to infer kinetic parameters from smFISH distributions. However, these approaches do not provide a probability distribution or likelihood of model parameters inferred from the smFISH data. This information is useful because it indicates which parameters are loosely constrained by the data and suggests follow up experiments. We developed a suite of MATLAB programs (BayFish) that estimate the Bayesian posterior probability of model parameters from smFISH data. The user specifies an underlying model of PLOS