Selection to minimise noise in living systems and its implications for the evolution of gene expression

Gene expression, like many biological processes, is subject to noise. This noise has been measured on a global scale, but its general importance to the fitness of an organism is unclear. Here, I show that noise in gene expression in yeast has evolved to prevent harmful stochastic variation in the levels of genes that reduce fitness when their expression levels change. Therefore, there has probably been widespread selection to minimise noise in gene expression. Selection to minimise noise, because it results in gene expression that is stable to stochastic variation in cellular components, may also constrain the ability of gene expression to respond to non‐stochastic variation. I present evidence that this has indeed been the case in yeast. I therefore conclude that gene expression noise is an important biological trait, and one that probably limits the evolvability of complex living systems.     

External noise and feedback regulation: steady-state statistics of auto-regulatory genetic network

The steady-state statistics of a single gene auto-regulatory genetic network with the additive external Gaussian white noises is investigated. The main result shows that the negative feedback will result in that the mRNA noise has a positive contribution to the protein noise, but the positive feedback will result in that the mRNA noise has a negative contribution to the protein noise. If there is no feed back, then the contribution of mRNA noise to protein noise is always positive. On the other hand, the analysis and numerical simulations of linear and nonlinear feedback show that it is possible that the negative feedback increases, but the positive feedback decreases, the protein noise.     

Cost-benefit theory and optimal design of gene regulation functions

Cells respond to the environment by regulating the expression of genes according to environmental signals. The relation between the input signal level and the expression of the gene is called the gene regulation function. It is of interest to understand the shape of a gene regulation function in terms of the environment in which it has evolved and the basic constraints of biological systems. Here we address this by presenting a cost-benefit theory for gene regulation functions that takes into account temporally varying inputs in the environment and stochastic noise in the biological components. We apply this theory to the well-studied lac operon of E. coli. The present theory explains the shape of this regulation function in terms of temporal variation of the input signals, and of minimizing the deleterious effect of cell-cell variability in regulatory protein levels. We also apply the theory to understand the evolutionary tradeoffs in setting the number of regulatory proteins and for selection of feed-forward loops in genetic circuits. The present cost-benefit theory can be used to understand the shape of other gene regulatory functions in terms of environment and noise constraints.     

The Consistent Phylogenetic Signal in Genome Trees Revealed by Reducing the Impact of Noise

Phylogenetic trees based on gene repertoires are remarkably similar to the current consensus of life history. Yet it has been argued that shared gene content is unreliable for phylogenetic reconstruction because of convergence in gene content due to horizontal gene transfer and parallel gene loss. Here we test this argument, by filtering out as noise those orthologous groups that have an inconsistent phylogenetic distribution, using two independent methods. The resulting phylogenies do indeed contain small but significant improvements. More importantly, we find that the majority of orthologous groups contain some phylogenetic signal and that the resulting phylogeny is the only detectable signal present in the gene distribution across genomes. Horizontal gene transfer or parallel gene loss does not cause systematic biases in the gene content tree.     

Noise and Information Transmission in Promoters with Multiple Internal States

Based on the measurements of noise in gene expression performed during the past decade, it has become customary to think of gene regulation in terms of a two-state model, where the promoter of a gene can stochastically switch between an ON and an OFF state. As experiments are becoming increasingly precise and the deviations from the two-state model start to be observable, we ask about the experimental signatures of complex multistate promoters, as well as the functional consequences of this additional complexity. In detail, we i), extend the calculations for noise in gene expression to promoters described by state transition diagrams with multiple states, ii), systematically compute the experimentally accessible noise characteristics for these complex promoters, and iii), use information theory to evaluate the channel capacities of complex promoter architectures and compare them with the baseline provided by the two-state model. We find that adding internal states to the promoter generically decreases channel capacity, except in certain cases, three of which (cooperativity, dual-role regulation, promoter cycling) we analyze in detail.     

An integer optimization algorithm for robust identification of non-linear gene regulatory networks

ABSTRACT: BACKGROUND: Reverse engineering gene networks and identifying regulatory interactions are integral to understanding cellular decision making processes. Advancement in high throughput experimental techniques has initiated innovative data driven analysis of gene regulatory networks. However, inherent noise associated with biological systems requires numerous experimental replicates for reliable conclusions. Furthermore, evidence of robust algorithms directly exploiting basic biological traits are few. Such algorithms are expected to be efficient in their performance and robust in their prediction. RESULTS: We have developed a network identification algorithm to accurately infer both the topology and strength of regulatory interactions from time series gene expression data in the presence of significant experimental noise and non-linear behavior. In this novel formulism, we have addressed data variability in biological systems by integrating network identification with the bootstrap resampling technique, hence predicting robust interactions from limited experimental replicates subjected to noise. Furthermore, we have incorporated non-linearity in gene dynamics using the S-system formulation. The basic network identification formulation exploits the trait of sparsity of biological interactions. Towards that, the identification algorithm is formulated as an integer-programming problem by introducing binary variables for each network component. The objective function is targeted to minimize the network connections subjected to the constraint of maximal agreement between the experimental and predicted gene dynamics. The developed algorithm is validated using both in-silico and experimental data-sets. These studies show that the algorithm can accurately predict the topology and connection strength of the in silico networks, as quantified by high precision and recall, and small discrepancy between the actual and predicted kinetic parameters. Furthermore, in both the in silico and experimental case studies, the predicted gene expression profiles are in very close agreement with the dynamics of the input data. CONCLUSIONS: Our integer programming algorithm effectively utilizes bootstrapping to identify robust gene regulatory networks from noisy, non-linear time-series gene expression data. With significant noise and non-linearities being inherent to biological systems, the present formulism, with the incorporation of network sparsity, is extremely relevant to gene regulatory networks, and while the formulation has been validated against in silico and E. Coli data, it can be applied to any biological system.     

Noise transmission during the dynamic pattern formation in fly embryos

Background: Developmental patterning is highly reproducible and accurate at the single-cell level during fly embryogenesis despite the gene expression noise and external perturbations such as the variation of the embryo length, temperature and genes. To reveal the underlying mechanism, it is very important to characterize the noise transmission during the dynamic pattern formation. Two hypotheses have been proposed. The “channel” scenario requires a highly reproducible input and an accurate interpretation by downstream genes. In contrast, the “filter” scenario proposes a noisy input and a noise filter via the cross-regulation of the downstream network. It has been under great debates which scenario the fly embryogenesis follows. Results: The first 3-h developmental patterning of fly embryos is orchestrated by a hierarchical segmentation gene network, which rewires upon the maternal to zygotic transition. Starting from the highly reproducible maternal gradients, the positional information is refined to the single-cell precision through the highly dynamical evolved zygotic gene expression profiles. Thus the fly embryo development might strictly fit into neither the originally proposed “filter” nor “channel” scenario. The controversy that which scenario the fly embryogenesis follows could be further clarified by combining quantitative measurements and modeling. Conclusions: Fly embryos have become one of the perfect model systems for quantitative systems biology studies. The underlying mechanism discovered from fly embryogenesis will deepen our understanding of the noise control of the gene network, facilitate searching for more efficient and safer methods for cell programming and reprogramming, and have the great potential for tissue engineering and regenerative medicine.