Individual cells in genetically homogeneous populations have been found to express different numbers of molecules of specific proteins. We investigated the origins of these variations in mammalian cells by counting individual molecules of mRNA produced from a reporter gene that was stably integrated into the cell’s genome. We found that there are massive variations in the number of mRNA molecules present in each cell. These variations occur because mRNAs are synthesized in short but intense bursts of transcription beginning when the gene transitions from an inactive to an active state and ending when they transition back to the inactive state. We show that these transitions are intrinsically random and not due to global, extrinsic factors such as the levels of transcriptional activators. Moreover, the gene activation causes burst-like expression of all genes within a wider genomic locus. We further found that bursts are also exhibited in the synthesis of natural genes. The bursts of mRNA expression can be buffered at the protein level by slow protein degradation rates. A stochastic model of gene activation and inactivation was developed to explain the statistical properties of the bursts. The model showed that increasing the level of transcription factors increases the average size of the bursts rather than their frequency. These results demonstrate that gene expression in mammalian cells is subject to large, intrinsically random fluctuations and raise questions about how cells are able to function in the face of such noise.

Summary. Motivated by the intriguing complexity of biochemical circuitry within individual cells we study Stochastic Chemical Reaction Networks (SCRNs), a formal model that considers a set of chemical reactions acting on a finite number of molecules in a well-stirred solution according to standard chemical kinetics equations. SCRNs have been widely used for describing naturally occurring (bio)chemical systems, and with the advent of synthetic biology they become a promising language for the design of artificial biochemical circuits. Our interest here is the computational power of SCRNs and how they relate to more conventional models of computation. We survey known connections and give new connections between SCRNs and

The partitioning and subsequent inheritance of cellular factors like proteins and RNAs is a ubiquitous feature of cell division. However, direct quantitative measures of how such nongenetic inheritance affects subsequent changes in gene expression have been lacking. We tracked families of the yeast Saccharomyces cerevisiae as they switch between two semi-stable epigenetic states. We found that long after two cells have divided, they continued to switch in a synchronized manner, whereas individual cells have exponentially distributed switching times. By comparing these results to a Poisson process, we show that the time evolution of an epigenetic state depends initially on inherited factors, with stochastic processes requiring several generations to decorrelate closely related cells. Finally, a simple stochastic model demonstrates that a single fluctuating regulatory protein that is synthesized in large bursts can explain the bulk of our results.

nvolves lower levels of noise than does the production of most other genes. Our results support the hypothesis that noise in gene expression is a biologically important variable, is generally detrimental to organismal fitness, and is subject to natural selection. Introduction Stochasticity is a ubiquitous characteristic of life. Such apparent randomness, or "noise," can be observed in a wide range of organisms, resulting in phenomena ranging from progressive loss of cell-cycle synchronization in an initially synchronized population of microbes to the pattern of hair coloration in female calico cats. An important source of stochasticity in biological systems is the random noise of transcription and translation, which can result in very different rates of synthesis of a specific protein in genetically identical cells in essentially identical environments (Elowitz et al. 2002; Ozbudak et al. 2002; Blake et al. 2003). Understanding how stochasticity contributes to cellular phenotypes is

Here, we describe a bicistronic reporter system for the analysis of promoter activity in a variety of Gram-negative bacteria at both the population and single-cell levels. This synthetic genetic tool utilizes an artificial operon comprising the gfp and lacZ genes that are assembled in a suicide vector, which is integrated at specific sites within the chromosome of the target bacterium, thereby creating a monocopy reporter system. This tool was instrumental for the complete in vivo characterization of two promoters, Pb and Pc, that drive the expression of the benzoate and catechol degradation pathways, respectively, of the soil bacterium Pseudomonas putida KT2440. The parameterization of these promoters in a population (using b-galactosidase assays) and in single cells (using flow cytometry) was necessary to examine the basic numerical features of these systems, such as the basal and maximal levels and the induction kinetics in response to an inducer (benzoate). Remarkably, GFP afforded a view of the process at a much higher resolution compared with standard lacZ tests; changes in fluorescence faithfully reflected variations in the transcriptional regimes of individual bacteria. The broad host range of the vector/reporter platform is an asset for the characterization of promoters in different bacteria, thereby

Biological oscillators naturally exhibit stochastic fluctuations in period and amplitude due to the random nature of molecular reactions. Accurately measuring the precision of noisy oscillators and the heterogeneity in period and strength of rhythmicity across a population of cells requires single-cell recordings of sufficient length to fully represent the variability of oscillations. We found persistent, independent circadian oscillations of clock gene expression in 6-week-long bioluminescence recordings of 80 primary fibroblast cells dissociated from PER2::LUC mice and kept in vitro for 6 months. Due to the stochastic nature of rhythmicity, the proportion of cells appearing rhythmic increases with the length of interval examined, with 100 % of cells found to be rhythmic when using 3-week windows. Mean period and amplitude are remarkably stable throughout the 6-week recordings, with precision improving over time. For individual cells, precision of period and amplitude are correlated with cell size and rhythm amplitude, but not with period, and period exhibits much less cycle-to-cycle variability (CV 7.3%) than does amplitude (CV 37%). The time series are long enough to distinguish stochastic fluctuations within each cell from differences among cells, and we conclude that the cells do exhibit significant heterogeneity in period and strength of rhythmicity, which we measure using a novel statistical metric. Furthermore, stochastic modeling suggests that these single-cell clocks operate near a Hopf bifurcation, such that intrinsic noise

Over the last few years, experimental data on the fluctuations in gene activity between individual cells and within the same cell over time have confirmed that gene expression is a ‘‘noisy’ ’ process. This variation is in part due to the small number of molecules taking part in some of the key reactions that are involved in gene expression. One of the consequences of this is that protein production often occurs in bursts, each due to a single promoter or transcription factor binding event. Recently, the distribution of the number of proteins produced in such bursts has been experimentally measured, offering a unique opportunity to study the relative importance of different sources of noise in gene expression. Here, we provide a derivation of the theoretical probability distribution of these bursts for a wide variety of different models of gene expression. We show that there is a good fit between our theoretical distribution and that obtained from two different published experimental datasets. We then prove that, irrespective of the details of the model, the burst size distribution is always geometric and hence determined by a single parameter. Many different combinations of the biochemical rates for the constituent reactions of both transcription and translation will therefore lead to the same experimentally observed burst size distribution. It is thus impossible to identify different sources of fluctuations purely from protein burst size data or to use such data to estimate all of the model parameters. We explore methods of inferring these values when additional types of experimental

Experiments in recent years have vividly demonstrated that gene expression can be highly stochastic. How protein concentration fluctuations affect the growth rate of a population of cells is, however, a wide-open question. We present a mathematical model that makes it possible to quantify the effect of protein concentration fluctuations on the growth rate of a population of genetically identical cells. The model predicts that the population’s growth rate depends on how the growth rate of a single cell varies with protein concentration, the variance of the protein concentration fluctuations, and the correlation time of these fluctuations. The model also predicts that when the average concentration of a protein is close to the value that maximizes the growth rate, fluctuations in its concentration always reduce the growth rate. However, when the average protein concentration deviates sufficiently from the optimal level, fluctuations can enhance the growth rate of the population, even when the growth rate of a cell depends linearly on the protein concentration. The model also shows that the ensemble or population average of a quantity, such as the average protein expression level or its variance, is in general not equal to its time average as obtained from tracing a single cell and its descendants. We apply our model to perform a cost-benefit analysis of gene regulatory control. Our analysis predicts that the optimal expression level of a gene regulatory protein is determined by the trade-off between the cost of synthesizing the regulatory protein and the benefit of minimizing the fluctuations in the expression of its target gene. We discuss possible experiments that could test our

A model-driven discovery process, Computing Life, is used to identify an ensemble of genetic networks that describe the biological clock. A clock mechanism involving the genes white-collar-1 and white-collar-2 (wc-1 and wc-2) that encode a transcriptional activator (as well as a blue-light receptor) and an oscillator frequency (frq) that encodes a cyclin that deactivates the activator is used to guide this discovery process through three cycles of microarray experiments. Central to this discovery process is a new methodology for the rational design of a Maximally Informative Next Experiment (MINE), based on the genetic network ensemble. In each experimentation cycle, the MINE approach is used to select the most informative new experiment in order to mine for clock-controlled genes, the outputs of the clock. As much as 25 % of the N. crassa transcriptome appears to be under clock-control. Clock outputs include genes with products in DNA metabolism, ribosome biogenesis in RNA metabolism, cell cycle, protein metabolism, transport, carbon metabolism, isoprenoid (including carotenoid) biosynthesis, development, and varied signaling processes. Genes under the transcription factor complex WCC ( = WC-1/WC-2) control were resolved into four classes, circadian only (612 genes), light-responsive only (396), both circadian and light-responsive (328), and neither circadian nor light-responsive (987). In each of three cycles of microarray experiments data support that wc-1 and wc-2 are auto-regulated by WCC. Among 11,000 N. crassa genes a total of 295 genes, including a large fraction of phosphatases/kinases, appear to be under the immediate control of the FRQ

Many cellular systems rely on the ability to interpret spatial heterogeneities in chemoattractant concentration to direct cell migration. The accuracy of this process is limited by stochastic fluctuations in the concentration of the external signal and in the internal signaling components. Here we use information theory to determine the optimal scheme to detect the location of an external chemoattractant source in the presence of noise. We compute the minimum amount of mutual information needed between the chemoattractant gradient and the internal signal to achieve a prespecified chemotactic accuracy. We show that more accurate chemotaxis requires greater mutual information. We also demonstrate that a priori information can improve chemotaxis efficiency. We compare the optimal signaling schemes with existing experimental measurements and models of eukaryotic gradient sensing. Remarkably, there is good quantitative agreement between the optimal response when no a priori assumption is made about the location of the existing source, and the observed experimental response of unpolarized Dictyostelium discoideum cells. In contrast, the measured response of polarized D. discoideum cells matches closely the optimal scheme, assuming prior knowledge of the external gradient—for example, through prolonged chemotaxis in a given direction. Our results demonstrate that different observed classes of responses in cells (polarized and unpolarized) are optimal under varying information assumptions. Citation: Andrews BW, Iglesias PA (2007) An information-theoretic characterization of the optimal gradient sensing response of cells. PLoS Comput Biol 3(8): e153. doi:10.