In search of noise-induced bimodality
© Kim and Sauro; licensee BioMed Central Ltd. 2012
Received: 20 September 2012
Accepted: 6 November 2012
Published: 7 November 2012
Many biological studies are carried out on large populations of cells, often in order to obtain enough material to make measurements. However, we now know that noise is endemic in biological systems and this results in cell-to-cell variability in what appears to be a population of identical cells. Although often neglected, this noise can have a dramatic effect on system responses to environmental cues with significant and often counter-intuitive biological outcomes. A recent study in BMC Systems Biology provides an example of this, documenting a bimodal distribution of activated extracellular signal-regulated kinase in a population of cells exposed to epidermal growth factor and demonstrating that the observed bimodality of the response is induced purely by noise.
See research article: http://0-www.biomedcentral.com.brum.beds.ac.uk/1752-0509/6/109
Noise in biological systems is endemic and can contribute to biological phenotypes. For example, noise can affect fate determination in virus-infected cells by randomly switching between latency and reactivation, and it can also cause Escherichia coli to switch between competency and non-competency for DNA uptake. The origin of biological noise is attributed to randomness in biological reaction events, and leads to cell-to-cell variability in genetically identical cell populations. To understand noise-induced phenotypes, it is important to have the capability to observe behaviors at the single-cell level.
The paper by Birtwistle et al.  highlights an interesting issue with respect to noise in the mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase (ERK) protein signaling pathway. This pathway receives environmental cues, such as changes in epidermal growth factor (EGF) that leads to an amplified signal in the form of activated ERK. The activation of ERK in turn results in diverse cellular responses, such as proliferation, differentiation, and apoptosis. In their study, flow cytometry was used to collect data in the form of fluorescence signals emitted from individual single cells. The signals indicated the activated ERK levels and showed a bimodal histogram in the cell population, which might otherwise have been obscured without adopting single-cell-level measurement techniques.
Bimodality can be the result of positive feedback mechanisms resulting in bistablity, a state very much like a light toggle switch. Such systems have been found in a number of natural systems so it was an obvious choice for Birtwistle et al. to investigate this possibility. However, although the ERK activation pathway has been found to exhibit bistability in some types of cell, it does not in others, and the experimental observations made by Birtwistle el al. in their cell line were not compatible with the mechanisms that used positive feedback. Closer investigation reveals something completely different and novel.
The observed bimodal distributions  were shown by computer modeling to be induced purely by noise without being related to positive feedback. There have been many theoretical studies to understand the mechanisms for such noise-induced bimodality. Many of these are related to the close interplay between nonlinearities in the system and noise in biological signals. The signal noise can be processed via a nonlinear input-output response, causing sufficient signal distortion to transform unimodal signal distributions into bimodal ones . In the case of the signaling pathway studied by Birtwistle et al., the EGF signal is transmitted across the membrane to generate activated Ras, and the Ras signal, considered unimodal in the cell line used, activates the MAPK/ERK cascade after passing a threshold. The existence of this threshold means that the system is not purely linear, and signal distortion due to the system nonlinearity can explain the bimodality that emerges. This provides a first in vivo example in protein signaling networks that shows noise-induced bimodality due to its inherent nonlinear signal processing without positive feedback.
Careful attention is, however, needed when computing distributions, more precisely histograms, depending on the choice of the scale of the x-axis. As in the case of flow cytometry, a signal is often visualized in the log scale to show the broad ranges of the signal values. This leads researchers to use the histogram of log-transformed signal values. One pitfall of this procedure is that it is possible to generate transformation artifacts. One example is shown in Figure 1, where a unimodal (or very weakly bimodal) distribution in a linear scale (Figure 1h) becomes a bimodal distribution in a log scale (Figure 1g). This is because the log-scale representation causes larger values of × to be compressed more visually, so that a greater number of samples will be taken at larger × values and fewer at smaller × values (Figure 1c,g). Therefore, flow cytometry data may need to be treated carefully. Birtwistle et al. represented their data in a log scale using Kaplan-Meier empirical cumulative distribution functions, which may reduce the problem. Although this method can visually amplify the bimodal distribution as the conventional log scale representation does, it was confirmed via personal communication that the bimodality robustly appears in the wide range of input doses (0.5 nM and 1 nM EGF in Figure 1B in Birtwistle et al. ) and the transformation artifact does not eliminate the bimodality except for the case of 0.1 nM EGF in the same figure.
From these simulations, we can conclude that the variability in the input signal and the activation threshold both individually enhanced bimodality in the output signal. This is because the value of y is determined by the distance between the value of × and the activation threshold. Thus, the variability in the threshold has the same effect as the input variability by changing the distance between the value of × and the threshold. What happens if variability appears in both × and the threshold simultaneously? If the input signal and the threshold can be assumed to fluctuate independently due to the fact that they can be processed through sufficiently different biological systems, then the presence of the variability in both will enhance the bimodality further.
Finally, we can consider the variability in the saturation levels of the response curves. This variability smoothed the sharp peak appearing at the saturation level (y = 1) for the case that the saturation level was fixed (Figure 2f). Thus, the variability in the saturation level did not have any significant effect on bimodality.
The lesson from this work, and one that we see more and more often, is that the interaction of noise and the underlying deterministic dynamics can result in non-intuitive behavior. We are only beginning to understand how noise is exploited by nature  and furthermore by system designers like synthetic biologists , but the influence of noise is likely to be subtle and counter-intuitive to our normal deterministic view of the world.
This work was supported by a National Science Foundation (NSF) Molecular Cell Biology Grant 1158573 and NIGMS grant GM081070.
- Birtwistle MR, Rauch J, Kiyatkin A, Aksamitiene E, Dobrzyński M, Hoek JB, Kolch W, Ogunnaike BA, Kholodenko BN: Emergence of bimodal cell population responses from the interplay between analog single-cell signaling and protein expression noise. BMC Syst Biol. 2012, 6: 109-10.1186/1752-0509-6-109.PubMedPubMed CentralView ArticleGoogle Scholar
- Ochab-Marcinek A, Tabaka M: Bimodal gene expression in noncooperative regulatory systems. Proc Natl Acad Sci USA. 2010, 107: 22096-22101. 10.1073/pnas.1008965107.PubMedPubMed CentralView ArticleGoogle Scholar
- Taniguchi Y, Choi PJ, Li GW, Chen H, Babu M, Hearn J, Emili A, Xie XS: Quantifying E. coli proteome and transcriptome with single-molecule sensitivity in single cells. Science. 2010, 329: 533-538. 10.1126/science.1188308.PubMedPubMed CentralView ArticleGoogle Scholar
- Eldar A, Elowitz MB: Functional roles for noise in genetic circuits. Nature. 2010, 467: 167-173. 10.1038/nature09326.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim KH, Sauro HM: Adjusting phenotypes by noise control. PLoS Comp Biol. 2012, 8: e1002344-10.1371/journal.pcbi.1002344.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.