Interplay between pleiotropy and secondary selection determines rise and fall of mutators in stress response

Dramatic rise of mutators has been found to accompany adaptation of bacteria in response to many kinds of stress. Two views on the evolutionary origin of this phenomenon emerged: the pleiotropic hypothesis positing that it is a byproduct of environmental stress or other specific stress response mechanisms and the second order selection which states that mutators hitchhike to fixation with unrelated beneficial alleles. Conventional population genetics models could not fully resolve this controversy because they are based on certain assumptions about fitness landscape. Here we address this problem using a microscopic multiscale model, which couples physically realistic molecular descriptions of proteins and their interactions with population genetics of carrier organisms without assuming any a priori fitness landscape. We found that both pleiotropy and second order selection play a crucial role at different stages of adaptation: the supply of mutators is provided through destabilization of error correction complexes or fluctuations of production levels of prototypic mismatch repair proteins (pleiotropic effects), while rise and fixation of mutators occur when there is a sufficient supply of beneficial mutations in replication-controlling genes. This general mechanism assures a robust and reliable adaptation of organisms to unforeseen challenges. This study highlights physical principles underlying physical biological mechanisms of stress response and adaptation

A Genome-Wide Analysis of Promoter-Mediated Phenotypic Noise in Escherichia coli

Gene expression is subject to random perturbations that lead to fluctuations in the rate of protein production. As a consequence, for any given protein, genetically identical organisms living in a constant environment will contain different amounts of that particular protein, resulting in different phenotypes. This phenomenon is known as “phenotypic noise.” In bacterial systems, previous studies have shown that, for specific genes, both transcriptional and translational processes affect phenotypic noise. Here, we focus on how the promoter regions of genes affect noise and ask whether levels of promoter-mediated noise are correlated with genes' functional attributes, using data for over 60% of all promoters in Escherichia coli. We find that essential genes and genes with a high degree of evolutionary conservation have promoters that confer low levels of noise. We also find that the level of noise cannot be attributed to the evolutionary time that different genes have spent in the genome of E. coli. In contrast to previous results in eukaryotes, we find no association between promoter-mediated noise and gene expression plasticity. These results are consistent with the hypothesis that, in bacteria, natural selection can act to reduce gene expression noise and that some of this noise is controlled through the sequence of the promoter region alon

The role of input noise in transcriptional regulation

Even under constant external conditions, the expression levels of genes fluctuate. Much emphasis has been placed on the components of this noise that are due to randomness in transcription and translation; here we analyze the role of noise associated with the inputs to transcriptional regulation, the random arrival and binding of transcription factors to their target sites along the genome. This noise sets a fundamental physical limit to the reliability of genetic control, and has clear signatures, but we show that these are easily obscured by experimental limitations and even by conventional methods for plotting the variance vs. mean expression level. We argue that simple, global models of noise dominated by transcription and translation are inconsistent with the embedding of gene expression in a network of regulatory interactions. Analysis of recent experiments on transcriptional control in the early Drosophila embryo shows that these results are quantitatively consistent with the predicted signatures of input noise, and we discuss the experiments needed to test the importance of input noise more generally.Comment: 11 pages, 5 figures minor correction

Effect of promoter architecture on the cell-to-cell variability in gene expression

According to recent experimental evidence, the architecture of a promoter, defined as the number, strength and regulatory role of the operators that control the promoter, plays a major role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effort that addresses the question of how changes in promoter architecture affect noise in gene expression in a systematic rather than case-by-case fashion. In this article, we make such a systematic investigation, based on a simple microscopic model of gene regulation that incorporates stochastic effects. In particular, we show how operator strength and operator multiplicity affect this variability. We examine different modes of transcription factor binding to complex promoters (cooperative, independent, simultaneous) and how each of these affects the level of variability in transcription product from cell-to-cell. We propose that direct comparison between in vivo single-cell experiments and theoretical predictions for the moments of the probability distribution of mRNA number per cell can discriminate between different kinetic models of gene regulation.Comment: 35 pages, 6 figures, Submitte

Accelerated expansion from ghost-free bigravity: a statistical analysis with improved generality

We study the background cosmology of the ghost-free, bimetric theory of gravity. We perform an extensive statistical analysis of the model using both frequentist and Bayesian frameworks and employ the constraints on the expansion history of the Universe from the observations of supernovae, the cosmic microwave background and the large scale structure to estimate the model's parameters and test the goodness of the fits. We explore the parameter space of the model with nested sampling to find the best-fit chi-square, obtain the Bayesian evidence, and compute the marginalized posteriors and mean likelihoods. We mainly focus on a class of sub-models with no explicit cosmological constant (or vacuum energy) term to assess the ability of the theory to dynamically cause a late-time accelerated expansion. The model behaves as standard gravity without a cosmological constant at early times, with an emergent extra contribution to the energy density that converges to a cosmological constant in the far future. The model can in most cases yield very good fits and is in perfect agreement with the data. This is because many points in the parameter space of the model exist that give rise to time-evolution equations that are effectively very similar to those of the $\Lambda$CDM. This similarity makes the model compatible with observations as in the $\Lambda$CDM case, at least at the background level. Even though our results indicate a slightly better fit for the $\Lambda$CDM concordance model in terms of the $p$-value and evidence, none of the models is statistically preferred to the other. However, the parameters of the bigravity model are in general degenerate. A similar but perturbative analysis of the model as well as more data will be required to break the degeneracies and constrain the parameters, in case the model will still be viable compared to the $\Lambda$CDM.Comment: 42 pages, 9 figures; typos corrected in equations (2.12), (2.13), (3.7), (3.8) and (3.9); more discussions added (footnotes 5, 8, 10 and 13) and abstract, sections 4.2, 4.3 and 5 (conclusions) modified in response to referee's comments; references added; acknowledgements modified; all results completely unchanged; matches version accepted for publication in JHE

A synthetic gene network for tuning protein degradation in Saccharomyces cerevisiae

Protein decay rates are regulated by degradation machinery that clears unnecessary housekeeping proteins and maintains appropriate dynamic resolution for transcriptional regulators. Turnover rates are also crucial for fluorescence reporters that must strike a balance between sufficient fluorescence for signal detection and temporal resolution for tracking dynamic responses. Here, we use components of the Escherichia coli degradation machinery to construct a Saccharomyces cerevisiae strain that allows for tunable degradation of a tagged protein. Using a microfluidic platform tailored for single-cell fluorescence measurements, we monitor protein decay rates after repression using an ssrA-tagged fluorescent reporter. We observe a half-life ranging from 91 to 22 min, depending on the level of activation of the degradation genes. Computational modeling of the underlying set of enzymatic reactions leads to GFP decay curves that are in excellent agreement with the observations, implying that degradation is governed by Michaelis–Menten-type interactions. In addition to providing a reporter with tunable dynamic resolution, our findings set the stage for explorations of the effect of protein degradation on gene regulatory and signalling pathways