Bet-hedging in stochastically switching environments

TOR regulates variability of protein synthesis rates

Cellular processes are subject to inherent variability, but the extent to which cells can regulate this variability has received little investigation. Here, we explore the characteristics of the rate of cellular protein synthesis in single cells of the eukaryote fission yeast. Strikingly, this rate is highly variable despite protein synthesis being dependent on hundreds of reactions which might be expected to average out at the overall cellular level. The rate is variable over short time scales, and exhibits homoeostatic behaviour at the population level. Cells can regulate the level of variability through processes involving the TOR pathway, suggesting there is an optimal level of variability conferring a selective advantage. While this could be an example of bet-hedging, but we propose an alternative explanation: regulated 'loose' control of complex processes of overall cellular metabolism such as protein synthesis, may lead to this variability. This could ensure cells are fluid in control and agile in response to changing conditions, and may constitute a novel organisational principle of complex metabolic cellular systems.

Hierarchic regulation of a metabolic pathway: H-NS, CRP, and SsrB control myo-inositol utilization by Salmonella enterica

IMPORTANCE

The capacity to utilize myo-inositol (MI) as sole carbon and energy source is widespread among bacteria, among them the intestinal pathogen S. Typhimurium. This study elucidates the complex and hierarchical regulation that underlies the utilization of MI by S. Typhimurium under substrate limitation. A total of seven regulatory factors have been identified so far, allowing the pathogen an environment-dependent, efficient, and fine-tuned regulation of a metabolic property that provides growth advantages in different environments.

The Migration and the Fate of Dental Pulp Stem Cells

Human dental pulp stem cells (hDPSCs) are adult mesenchymal stem cells (MSCs) obtained from dental pulp and derived from the neural crest. They can differentiate into odontoblasts, osteoblasts, chondrocytes, adipocytes and nerve cells, and they play a role in tissue repair and regeneration. In fact, DPSCs, depending on the microenvironmental signals, can differentiate into odontoblasts and regenerate dentin or, when transplanted, replace/repair damaged neurons. Cell homing depends on recruitment and migration, and it is more effective and safer than cell transplantation. However, the main limitations of cell homing are the poor cell migration of MSCs and the limited information we have on the regulatory mechanism of the direct differentiation of MSCs. Different isolation methods used to recover DPSCs can yield different cell types. To date, most studies on DPSCs use the enzymatic isolation method, which prevents direct observation of cell migration. Instead, the explant method allows for the observation of single cells that can migrate at two different times and, therefore, could have different fates, for example, differentiation and self-renewal. DPSCs use mesenchymal and amoeboid migration modes with the formation of lamellipodia, filopodia and blebs, depending on the biochemical and biophysical signals of the microenvironment. Here, we present current knowledge on the possible intriguing role of cell migration, with particular attention to microenvironmental cues and mechanosensing properties, in the fate of DPSCs.

Do microbes have a memory? History-dependent behavior in the adaptation to variable environments

Microbes are constantly confronted with changes and challenges in their environment. A proper response to these environmental cues is needed for optimal cellular functioning and fitness. Interestingly, past exposure to environmental cues can accelerate or boost the response when this condition returns, even in daughter cells that have not directly encountered the initial cue. Moreover, this behavior is mostly epigenetic and often goes hand in hand with strong heterogeneity in the strength and speed of the response between isogenic cells of the same population, which might function as a bet-hedging strategy. In this review, we discuss examples of history-dependent behavior (HDB) or “memory,” with a specific focus on HDB in fluctuating environments. In most examples discussed, the lag time before the response to an environmental change is used as an experimentally measurable proxy for HDB. We highlight different mechanisms already implicated in HDB, and by using HDB in fluctuating carbon conditions as a case study, we showcase how the metabolic state of a cell can be a key determining factor for HDB. Finally, we consider possible evolutionary causes and consequences of such HDB.

Life-History traits and the replicator equation

Due to the relevance for conservation biology, there is an increasing interest to extend evolutionary genomics models to plant, animal or microbial species. However, this requires to understand the effect of life-history traits absent in humans on genomic evolution. In this context, it is fundamentally of interest to generalize the replicator equation, which is at the heart of most population genomics models. However, as the inclusion of life-history traits generates models with a large state space, the analysis becomes involving. We focus, here, on quiescence and seed banks, two features common to many plant, invertebrate and microbial species. We develop a method to obtain a low-dimensional replicator equation in the context of evolutionary game theory, based on two assumptions: (1) the life-history traits are per se neutral, and (2) frequency-dependent selection is weak. We use the results to investigate the evolution and maintenance of cooperation based on the Prisoner's dilemma and the snowdrift game. We first consider the generalized replicator equation, and then refine the investigation using adaptive dynamics. It turns out that, depending on the structure and timing of the quiescence/dormancy life-history trait, cooperation in a homogeneous population can be stabilized. We finally discuss and highlight the relevance of these results for plant, invertebrate and microbial communities.

The role of memory in non-genetic inheritance and its impact on cancer treatment resistance

Intra-tumour heterogeneity is a leading cause of treatment failure and disease progression in cancer. While genetic mutations have long been accepted as a primary mechanism of generating this heterogeneity, the role of phenotypic plasticity is becoming increasingly apparent as a driver of intra-tumour heterogeneity. Consequently, understanding the role of this plasticity in treatment resistance and failure is a key component of improving cancer therapy. We develop a mathematical model of stochastic phenotype switching that tracks the evolution of drug-sensitive and drug-tolerant subpopulations to clarify the role of phenotype switching on population growth rates and tumour persistence. By including cytotoxic therapy in the model, we show that, depending on the strategy of the drug-tolerant subpopulation, stochastic phenotype switching can lead to either transient or permanent drug resistance. We study the role of phenotypic heterogeneity in a drug-resistant, genetically homogeneous population of non-small cell lung cancer cells to derive a rational treatment schedule that drives population extinction and avoids competitive release of the drug-tolerant sub-population. This model-informed therapeutic schedule results in increased treatment efficacy when compared against periodic therapy, and, most importantly, sustained tumour decay without the development of resistance.

Persistence as an Optimal Hedging Strategy

Bacteria invest in a slow-growing subpopulation, called persisters, to ensure survival in the face of uncertainty. This hedging strategy is remarkably similar to financial hedging, where diversifying an investment portfolio protects against economic uncertainty. We provide a new, to our knowledge, theoretical foundation for understanding cellular hedging by unifying the study of biological population dynamics and the mathematics of financial risk management through optimal control theory. Motivated by the widely accepted role of volatility in the emergence of persistence, we consider several models of environmental volatility described by continuous-time stochastic processes. This allows us to study an emergent cellular hedging strategy that maximizes the expected per capita growth rate of the population. Analytical and simulation results probe the optimal persister strategy, revealing results that are consistent with experimental observations and suggest new opportunities for experimental investigation and design. Overall, we provide a new, to our knowledge, way of conceptualizing and modeling cellular decision making in volatile environments by explicitly unifying theory from mathematical biology and finance.

Comparative study between discrete and continuum models for the evolution of competing phenotype-structured cell populations in dynamical environments

Deterministic continuum models formulated as nonlocal partial differential equations for the evolutionary dynamics of populations structured by phenotypic traits have been used recently to address open questions concerning the adaptation of asexual species to periodically fluctuating environmental conditions. These models are usually defined on the basis of population-scale phenomenological assumptions and cannot capture adaptive phenomena that are driven by stochastic variability in the evolutionary paths of single individuals. In light of these considerations, in this paper we develop a stochastic individual-based model for the coevolution of two competing phenotype-structured cell populations that are exposed to time-varying nutrient levels and undergo spontaneous, heritable phenotypic changes with different probabilities. Here, the evolution of every cell is described by a set of rules that result in a discrete-time branching random walk on the space of phenotypic states, and nutrient levels are governed by a difference equation in which a sink term models nutrient consumption by the cells. We formally show that the deterministic continuum counterpart of this model comprises a system of nonlocal partial differential equations for the cell population density functions coupled with an ordinary differential equation for the nutrient concentration. We compare the individual-based model and its continuum analog, focusing on scenarios whereby the predictions of the two models differ. The results obtained clarify the conditions under which significant differences between the two models can emerge due to bottleneck effects that bring about both lower regularity of the density functions of the two populations and more pronounced demographic stochasticity. In particular, bottleneck effects emerge in the presence of lower probabilities of phenotypic variation and are more apparent when the two populations are characterized by lower fitness initial mean phenotypes and smaller initial levels of phenotypic heterogeneity. The emergence of these effects, and thus the agreement between the two modeling approaches, is also dependent on the initial proportions of the two populations. As an illustrative example, we demonstrate the implications of these results in the context of the mathematical modeling of the early stage of metastatic colonization of distant organs.

Evolutionary dynamics of competing phenotype-structured populations in periodically fluctuating environments

Living species, ranging from bacteria to animals, exist in environmental conditions that exhibit spatial and temporal heterogeneity which requires them to adapt. Risk-spreading through spontaneous phenotypic variations is a known concept in ecology, which is used to explain how species may survive when faced with the evolutionary risks associated with temporally varying environments. In order to support a deeper understanding of the adaptive role of spontaneous phenotypic variations in fluctuating environments, we consider a system of non-local partial differential equations modelling the evolutionary dynamics of two competing phenotype-structured populations in the presence of periodically oscillating nutrient levels. The two populations undergo heritable, spontaneous phenotypic variations at different rates. The phenotypic state of each individual is represented by a continuous variable, and the phenotypic landscape of the populations evolves in time due to variations in the nutrient level. Exploiting the analytical tractability of our model, we study the long-time behaviour of the solutions to obtain a detailed mathematical depiction of the evolutionary dynamics. The results suggest that when nutrient levels undergo small and slow oscillations, it is evolutionarily more convenient to rarely undergo spontaneous phenotypic variations. Conversely, under relatively large and fast periodic oscillations in the nutrient levels, which bring about alternating cycles of starvation and nutrient abundance, higher rates of spontaneous phenotypic variations confer a competitive advantage. We discuss the implications of our results in the context of cancer metabolism.

Understanding the role of phenotypic switching in cancer drug resistance

The emergence of acquired drug resistance in cancer represents a major barrier to treatment success. While research has traditionally focused on genetic sources of resistance, recent findings suggest that cancer cells can acquire transient resistant phenotypes via epigenetic modifications and other non-genetic mechanisms. Although these resistant phenotypes are eventually relinquished by individual cells, they can temporarily 'save' the tumor from extinction and enable the emergence of more permanent resistance mechanisms. These observations have generated interest in the potential of epigenetic therapies for long-term tumor control or eradication. In this work, we develop a mathematical model to study how phenotypic switching at the single-cell level affects resistance evolution in cancer. We highlight unique features of non-genetic resistance, probe the evolutionary consequences of epigenetic drugs and explore potential therapeutic strategies. We find that even short-term epigenetic modifications and stochastic fluctuations in gene expression can drive long-term drug resistance in the absence of any bona fide resistance mechanisms. We also find that an epigenetic drug that slightly perturbs the average retention of the resistant phenotype can turn guaranteed treatment failure into guaranteed success. Lastly, we find that combining an epigenetic drug with an anti-cancer agent can significantly outperform monotherapy, and that treatment outcome is heavily affected by drug sequencing.

Model genotype-phenotype mappings and the algorithmic structure of evolution

Cancers are complex dynamic systems that undergo evolution and selection. Personalized medicine approaches in the clinic increasingly rely on predictions of tumour response to one or more therapies; these predictions are complicated by the inevitable evolution of the tumour. Despite enormous amounts of data on the mutational status of cancers and numerous therapies developed in recent decades to target these mutations, many of these treatments fail after a time due to the development of resistance in the tumour. The emergence of these resistant phenotypes is not easily predicted from genomic data, since the relationship between genotypes and phenotypes, termed the genotype-phenotype (GP) mapping, is neither injective nor functional. We present a review of models of this mapping within a generalized evolutionary framework that takes into account the relation between genotype, phenotype, environment and fitness. Different modelling approaches are described and compared, and many evolutionary results are shown to be conserved across studies despite using different underlying model systems. In addition, several areas for future work that remain understudied are identified, including plasticity and bet-hedging. The GP-mapping provides a pathway for understanding the potential routes of evolution taken by cancers, which will be necessary knowledge for improving personalized therapies.

Are the better cooperators dormant or quiescent?

Despite the wealth of empirical and theoretical studies, the origin and maintenance of cooperation is still an evolutionary riddle. In this context, ecological life-history traits which affect the efficiency of selection may play a role despite being often ignored. We consider here species such as bacteria, fungi, invertebrates and plants which exhibit resting stages in the form of a quiescent state or a seed bank. When quiescent, individuals are inactive and reproduce upon activation, while under seed bank parents produce offspring remaining dormant for different amount of time. We assume weak frequency-dependent selection modeled using game-theory and the prisoner's dilemma (cooperation/defect) as payoff matrix. The cooperators and defectors are allowed to evolve different quiescence or dormancy times. By means of singular perturbation theory we reduce the model to a one-dimensional equation resembling the well known replicator equation, in which the gain functions are scaled with lumped parameters reflecting the time scale of the resting state of the cooperators and defectors. If both time scales are identical cooperation cannot persist in a homogeneous population. If, however, the time scale of the cooperator is distinctively different from that of the defector, cooperation may become a locally asymptotically stable strategy. Interestingly enough, in the seed bank case the cooperator needs to become active faster than the defector, while in the quiescent case the cooperator has to be slower. We use adaptive dynamics to identify situations where cooperation may evolve and form a convergent stable ESS. We conclude by highlighting the relevance of these results for many non-model species and the maintenance of cooperation in microbial, invertebrate or plant populations.

Mechanisms behind persistence of a fire-sensitive alternative stable state system in the Gibson Desert, Western Australia

Alternative vegetation types that switch from one to another under contrasting fire regimes are termed fire-mediated alternative stable states (FMASS). Typically, pyrophylic communities (i.e., vegetation assemblages favored by burning) dominate under high frequencies or intensities of fire. Conversely, fire-sensitive (pyrophobic) vegetation types persist under long fire-free conditions. As the persistence traits of plants of FMASS systems are generally poorly researched, threshold levels of pyric disturbance that trigger 'state-switching' are often unknown. Dense thickets of the obligate-seeder shrub waputi (Aluta maisonneuvei ssp. maisonneuvei [Myrtaceae]) form fire-retarding woody islands within highly flammable spinifex (Triodia spp.) grasslands in arid Australia. To examine the tolerance of Aluta thickets to burning, we investigated: (1) the influence of post-fire rainfall and fire severity on recruitment (a field study); (2) soil seedbank densities (a field study); and (3) fire-related dormancy cues in seeds (a germination trial). We found a positive relationship between recruitment and post-fire rainfall volume, and much higher mean recruitment at sites with high- (5.9 seedlings/m2) than low-severity-burnt (2.2 seedlings/m2) and unburnt shrubs (0.03 seedlings/m2). Post-fire regeneration was mediated by dense soil-borne seedbanks, and the germination trial indicated that smoke promoted germination. Although Aluta shrubs are invariably fire-killed, high-severity fires are unlikely to lead to state shifts from shrubland to grassland because of the ability of mature stands to regenerate from dense, fire-cued seedbanks. Nevertheless, given that Aluta seedlings are exceptionally slow-growing, post-fire droughts combined with fire-return intervals less than the Aluta primary juvenile period of c. 5 years could drive conversion from Aluta- to Triodia-dominated vegetation.

Evolutionary Stability of Salmonella Competition with the Gut Microbiota: How the Environment Fosters Heterogeneity in Exploitative and Interference Competition

Following ingestion, gastrointestinal pathogens compete against the gastrointestinal microbiota and overcome host immune defenses in order to cause infections. Besides employing direct killing mechanisms, the commensal microbiota occupies metabolic niches to outcompete invading pathogens. Salmonella enterica serovar Typhimurium (S. Typhimurium) uses several strategies to successfully colonize the gut and establish infection, of which an increasing number is based on phenotypic heterogeneity within the S. Typhimurium population. The utilization of myo-inositol (MI) and the production of colicin confer a selective advantage over the microbiota in terms of exploitative and interference competition, respectively. In this review, we summarize the genetic basis underlying bistability of MI catabolism and colicin production. As demonstrated by single-cell analyses, a stochastic switch in the expression of the genes responsible for colicin production and MI degradation constitutes the heterogeneity of the two phenotypes. Both genetic systems are tightly regulated to avoid their expression under non-appropriate conditions and possible detrimental effects on bacterial fitness. Moreover, evolutionary mechanisms underlying formation and stability of these phenotypes in S. Typhimurium are discussed. We propose that both MI catabolism and colicin production create a bet-hedging strategy, which provides an adaptive benefit for S. Typhimurium in the fluctuating environment of the mammalian gut.

Evolutionary model for the unequal segregation of high copy plasmids

Plasmids are extrachromosomal DNA elements of microorganisms encoding beneficial genetic information. They were thought to be equally distributed to daughter cells during cell division. Here we use mathematical modeling to investigate the evolutionary stability of plasmid segregation for high-copy plasmids—plasmids that are present in up to several hundred copies per cell—carrying antibiotic resistance genes. Evolutionary stable strategies (ESS) are determined by numerical analysis of a plasmid-load structured population model. The theory predicts that the evolutionary stable segregation strategy of a cell depends on the plasmid copy number: For low and medium plasmid load, both daughters receive in average an equal share of plasmids, while in case of high plasmid load, one daughter obtains distinctively and systematically more plasmids. These findings are in good agreement with recent experimental results. We discuss the interpretation and practical consequences.

In the last years, it becomes more and more clear that heterogeneity in isogenic bacterial populations is rather the rule than the exception. This observation is interesting as it reveals the complex social life of bacteria, and also because of tremendous practical implications in medicine, biotechnology, and ecology. The central questions in this field are the identification of the underlying proximate causes (molecular mechanisms) on the one hand and on the other hand the identification of ultimate causes (evolutionary forces) that shape the social life of bacteria. We focus on plasmid dynamics, in particular on plasmid segregation. Recent experiments showed that plasmid segregation depends on the plasmid load. We identify possible evolutionary factors that shaped this process. It turns out that the ambivalence in the effect of plasmids—advantageous if present in low copy numbers, a metabolic burden if present in high copy numbers—is able to explain the experimental observations. The experimental findings can be interpreted as a variant of the principle of division of labor, as it is well known from e.g. persister cells or sporulation. Our model extends the theory of unequal segregation of damage to the ambivalent role of plasmids. Similarly, it is known that certain gene regulatory proteins are acting in a dose-dependent manner. Due to differences in their cellular concentrations and in their affinities to various target promoters, differential gene expression patterns are achieved. Consequently, tight concentration control is observed [1]. Another example is the strict copy-dependent utilization of autolysins during cell division. These enzymes carefully open the bacterial cell wall to allow for its extension [2]. Overproduction of these enzymes leads to cell lysis [3]. These molecules are in principle also candidates for a segregation characteristic similar to that of high copy plasmids described here.

Shortsighted Evolution Constrains the Efficacy of Long-Term Bet Hedging

To survive unpredictable environmental change, many organisms adopt bet-hedging strategies that are initially costly but provide a long-term fitness benefit. The temporal extent of these deferred fitness benefits determines whether bet-hedging organisms can survive long enough to realize them. In this article, we examine a model of microbial bet hedging in which there are two paths to extinction: unpredictable environmental change and demographic stochasticity. In temporally correlated environments, these drivers of extinction select for different switching strategies. Rapid phenotype switching ensures survival in the face of unpredictable environmental change, while slower-switching organisms become extinct. However, when both switching strategies are present in the same population, then demographic stochasticity-enforced by a limited population size-leads to extinction of the faster-switching organism. As a result, we find a novel form of evolutionary suicide whereby selection in a fluctuating environment can favor bet-hedging strategies that ultimately increase the risk of extinction. Population structures with multiple subpopulations and dispersal can reduce the risk of extinction from unpredictable environmental change and shift the balance so as to facilitate the evolution of slower-switching organisms.

Exploration-exploitation tradeoffs dictate the optimal distributions of phenotypes for populations subject to fitness fluctuations

We study a minimal model for the growth of a phenotypically heterogeneous population of cells subject to a fluctuating environment in which they can replicate (by exploiting available resources) and modify their phenotype within a given landscape (thereby exploring novel configurations). The model displays an exploration-exploitation trade-off whose specifics depend on the statistics of the environment. Most notably, the phenotypic distribution corresponding to maximum population fitness (i.e., growth rate) requires a nonzero exploration rate when the magnitude of environmental fluctuations changes randomly over time, while a purely exploitative strategy turns out to be optimal in two-state environments, independently of the statistics of switching times. We obtain analytical insight into the limiting cases of very fast and very slow exploration rates by directly linking population growth to the features of the environment.

Effect of a locally adapted genome on environmentally induced epigenetic variation

Both genetic variation and environmentally induced epigenetic changes allow organisms to persist through the heterogeneity of their habitats. Selection on genetic variation can promote local adaptation of populations. However, in absence of genetic variation, clonal organisms mostly rely on epigenetics to respond to environmental heterogeneity. We used the potential of unisexual organisms in incorporating their host genome, to empirically assess whether the presence of a locally adapted genome affects environmentally induced epigenetic changes in clonal organisms. We addressed this problematic by using unisexual lineages of the kleptogen vertebrate Ambystoma laterale-jeffersonianum complex that can optionally incorporate genetic material from locally adapted sexual hosts through genomic exchanges. More specifically, we compared environmentally induced epigenetic changes between lineages strictly reproducing clonally vs. those incorporating a locally adapted genome. The results revealed that both lineage and sample site components, as well as their interaction, affected epigenetic variation. When lineages were analysed separately, differences among sample sites were only detected in lineages impervious to genomic exchanges. Sample sites had no significant effect on the epigenetic variation of lineages that performed genomic exchanges. These results suggest that environmentally induced epigenetic variation among sites depends more on the lack of locally adapted alleles than on the level of genetic variation.

Eco-evolutionary causes and consequences of temporal changes in intratumoural blood flow

Temporal changes in blood flow are commonly observed in malignant tumours, but the evolutionary causes and consequences are rarely considered. We propose that stochastic temporal variations in blood flow and microenvironmental conditions arise from the eco-evolutionary dynamics of tumour angiogenesis in which cancer cells, as individual units of selection, can influence and respond only to local environmental conditions. This leads to new vessels arising from the closest available vascular structure regardless of the size or capacity of this parental vessel. These dynamics produce unstable vascular networks with unpredictable spatial and temporal variations in blood flow and microenvironmental conditions. Adaptations of evolving populations to temporally varying environments in nature include increased diversity, greater motility and invasiveness, and highly plastic phenotypes, allowing for broad metabolic adaptability and rapid shifts to high rates of proliferation and profound quiescence. These adaptive strategies, when adopted in cancer cells, promote many commonly observed phenotypic properties including those found in the stem phenotype and in epithelial-to-mesenchymal transition. Temporal variations in intratumoural blood flow, which occur through the promotion of cancer cell phenotypes that facilitate both metastatic spread and resistance to therapy, may have substantial clinical consequences.

Links between evolutionary processes and phenotypic robustness in microbes

The costs and benefits of random phenotypic heterogeneity in microbes have been vigorously debated and experimental tested for decades; yet, this conversation is largely independent from discussion of phenotypic robustness in other disciplines. In this review I connect microbial examples of stochasticity with studies on the ecological and population-genetic consequences of phenotypic variability. These topics illustrate the complexity of selection pressures on phenotypic robustness and provide inspiration that this complexity can be parsed with theoretical advances and the experimental power of microbial systems.

Simulation-Based Exploration of Quorum Sensing Triggered Resistance of Biofilms to Antibiotics

We present a mathematical model of biofilm response to antibiotics, controlled by a quorum sensing system that confers increased resistance. The model is a highly nonlinear system of partial differential equations that we investigate in computer simulations. Our results suggest that an adaptive, quorum sensing-controlled, mechanism to switch between modes of fast growth with little protection and protective modes of slow growth may confer benefits to biofilm populations. It enhances the formation of micro-niches in the inner regions of the biofilm in which bacteria are not easily reached by antibiotics. Whereas quorum sensing inhibitors can delay the onset of increased resistance, their advantage is lost after up-regulation. This emphasizes the importance of timing for treatment of biofilms with antibiotics.

Population heterogeneity in microbial bioprocesses: origin, analysis, mechanisms, and future perspectives

Population heterogeneity is omnipresent in all bioprocesses even in homogenous environments. Its origin, however, is only so well understood that potential strategies like bet-hedging, noise in gene expression and division of labour that lead to population heterogeneity can be derived from experimental studies simulating the dynamics in industrial scale bioprocesses. This review aims at summarizing the current state of the different parts of single cell studies in bioprocesses. This includes setups to visualize different phenotypes of single cells, computational approaches connecting single cell physiology with environmental influence and special cultivation setups like scale-down reactors that have been proven to be useful to simulate large-scale conditions. A step in between investigation of populations and single cells is studying subpopulations with distinct properties that differ from the rest of the population with sub-omics methods which are also presented here. Moreover, the current knowledge about population heterogeneity in bioprocesses is summarized for relevant industrial production hosts and mixed cultures, as they provide the unique opportunity to distribute metabolic burden and optimize production processes in a way that is impossible in traditional monocultures. In the end, approaches to explain the underlying mechanism of population heterogeneity and the evidences found to support each hypothesis are presented. For instance, population heterogeneity serving as a bet-hedging strategy that is used as coordinated action against bioprocess-related stresses while at the same time spreading the risk between individual cells as it ensures the survival of least a part of the population in any environment the cells encounter.

Phenotypic Diversity Using Bimodal and Unimodal Expression of Stress Response Proteins

Populations of cells need to express proteins to survive the sudden appearance of stressors. However, these mechanisms may be taxing. Populations can introduce diversity, allowing individual cells to stochastically switch between fast-growing and stress-tolerant states. One way to achieve this is to use genetic networks coupled with noise to generate bimodal distributions with two distinct subpopulations, each adapted to a stress condition. Another survival strategy is to rely on random fluctuations in gene expression to produce continuous, unimodal distributions of the stress response protein. To quantify the environmental conditions where bimodal versus unimodal expression is beneficial, we used a differential evolution algorithm to evolve optimal distributions of stress response proteins given environments with sudden fluctuations between low and high stress. We found that bimodality evolved for a large range of environmental conditions. However, we asked whether these findings were an artifact of considering two well-defined stress environments (low and high stress). As noise in the environment increases, or when there is an intermediate environment (medium stress), the benefits of bimodality decrease. Our results indicate that under realistic conditions, a continuum of resistance phenotypes generated through a unimodal distribution is sufficient to ensure survival without a high cost to the population.

Delayed bet-hedging resilience strategies under environmental fluctuations

Many biological populations, such as bacterial colonies, have developed through evolution a protection mechanism, called bet hedging, to increase their probability of survival under stressful environmental fluctuation. In this context, the concept of preadaptation refers to a common type of bet-hedging protection strategy in which a relatively small number of individuals in a population stochastically switch their phenotypes to a dormant metabolic state in which they increase their probability of survival against potential environmental shocks. Hence, if an environmental shock took place at some point in time, preadapted organisms would be better adapted to survive and proliferate once the shock is over. In many biological populations, the mechanisms of preadaptation and proliferation present delays whose influence in the fitness of the population are not well understood. In this paper, we propose a rigorous mathematical framework to analyze the role of delays in both preadaptation and proliferation mechanisms in the survival of biological populations, with an emphasis on bacterial colonies. Our theoretical framework allows us to analytically quantify the average growth rate of a bet-hedging bacterial colony with stochastically delayed reactions with arbitrary precision. We verify the accuracy of the proposed method by numerical simulations and conclude that the growth rate of a bet-hedging population shows a nontrivial dependency on their preadaptation and proliferation delays. Contrary to the current belief, our results show that faster reactions do not, in general, increase the overall fitness of a biological population.

Stochasticity in the Genotype-Phenotype Map: Implications for the Robustness and Persistence of Bet-Hedging

Nongenetic variation in phenotypes, or bet-hedging, has been observed as a driver of drug resistance in both bacterial infections and cancers. Here, we study how bet-hedging emerges in genotype-phenotype (GP) mapping through a simple interaction model: a molecular switch. We use simple chemical reaction networks to implement stochastic switches that map gene products to phenotypes, and investigate the impact of structurally distinct mappings on the evolution of phenotypic heterogeneity. Bet-hedging naturally emerges within this model, and is robust to evolutionary loss through mutations to both the expression of individual genes, and to the network itself. This robustness explains an apparent paradox of bet-hedging-why does it persist in environments where natural selection necessarily acts to remove it? The structure of the underlying molecular mechanism, itself subject to selection, can slow the evolutionary loss of bet-hedging to ensure a survival mechanism against environmental catastrophes even when they are rare. Critically, these properties, taken together, have profound implications for the use of treatment-holidays to combat bet-hedging-driven resistant disease, as the efficacy of breaks from treatment will ultimately be determined by the structure of the GP mapping.

The lag-phase during diauxic growth is a trade-off between fast adaptation and high growth rate

Bi-phasic or diauxic growth is often observed when microbes are grown in a chemically defined medium containing two sugars (for example glucose and lactose). Typically, the two growth stages are separated by an often lengthy phase of arrested growth, the so-called lag-phase. Diauxic growth is usually interpreted as an adaptation to maximise population growth in multi-nutrient environments. However, the lag-phase implies a substantial loss of growth during the switch-over. It therefore remains unexplained why the lag-phase is adaptive. Here we show by means of a stochastic simulation model based on the bacterial PTS system that it is not possible to shorten the lag-phase without incurring a permanent growth-penalty. Mechanistically, this is due to the inherent and well established limitations of biological sensors to operate efficiently at a given resource cost. Hence, there is a trade-off between lost growth during the diauxic switch and the long-term growth potential of the cell. Using simulated evolution we predict that the lag-phase will evolve depending on the distribution of conditions experienced during adaptation. In environments where switching is less frequently required, the lag-phase will evolve to be longer whereas, in frequently changing environments, the lag-phase will evolve to be shorter.

The dynamic balance of import and export of zinc in Escherichia coli suggests a heterogeneous population response to stress

Zinc is essential for life, but toxic in excess. Thus all cells must control their internal zinc concentration. We used a systems approach, alternating rounds of experiments and models, to further elucidate the zinc control systems in Escherichia coli. We measured the response to zinc of the main specific zinc import and export systems in the wild-type, and a series of deletion mutant strains. We interpreted these data with a detailed mathematical model and Bayesian model fitting routines. There are three key findings: first, that alternate, non-inducible importers and exporters are important. Second, that an internal zinc reservoir is essential for maintaining the internal zinc concentration. Third, our data fitting led us to propose that the cells mount a heterogeneous response to zinc: some respond effectively, while others die or stop growing. In a further round of experiments, we demonstrated lower viable cell counts in the mutant strain tested exposed to excess zinc, consistent with this hypothesis. A stochastic model simulation demonstrated considerable fluctuations in the cellular levels of the ZntA exporter protein, reinforcing this proposal. We hypothesize that maintaining population heterogeneity could be a bet-hedging response allowing a population of cells to survive in varied and fluctuating environments.

Variability in thermal and phototactic preferences in Drosophila may reflect an adaptive bet-hedging strategy

Organisms use various strategies to cope with fluctuating environmental conditions. In diversified bet‐hedging, a single genotype exhibits phenotypic heterogeneity with the expectation that some individuals will survive transient selective pressures. To date, empirical evidence for bet‐hedging is scarce. Here, we observe that individual Drosophila melanogaster flies exhibit striking variation in light‐ and temperature‐preference behaviors. With a modeling approach that combines real world weather and climate data to simulate temperature preference‐dependent survival and reproduction, we find that a bet‐hedging strategy may underlie the observed interindividual behavioral diversity. Specifically, bet‐hedging outcompetes strategies in which individual thermal preferences are heritable. Animals employing bet‐hedging refrain from adapting to the coolness of spring with increased warm‐seeking that inevitably becomes counterproductive in the hot summer. This strategy is particularly valuable when mean seasonal temperatures are typical, or when there is considerable fluctuation in temperature within the season. The model predicts, and we experimentally verify, that the behaviors of individual flies are not heritable. Finally, we model the effects of historical weather data, climate change, and geographic seasonal variation on the optimal strategies underlying behavioral variation between individuals, characterizing the regimes in which bet‐hedging is advantageous.

Slow and steady wins the race: a bacterial exploitative competition strategy in fluctuating environments

One promising frontier for synthetic biology is the development of synthetic ecologies, whereby interacting species form an additional layer of connectivity for engineered gene circuits. Toward this goal, an important step is to understand different types of bacterial interactions in natural settings, among which competition is the most prevalent. By constructing a two-species population dynamics model, here, we mimicked bacterial growth in nature with resource-limited fluctuating environments and searched for optimal strategies for bacterial exploitative competition. In a simple game with two strategy options (constant or susceptible growth), we found that the species playing the constant growth strategy always outplays or is evenly matched with its competitor, suggesting that constant growth is a "no-loss" good bet. We also showed that adoption of sophisticated strategies enables a species to maximize its fitness when its competitor grows susceptibly. The pursuit of fitness maximization is, however, associated with potential loss if both species are capable of strategy adjustment, indicating an intrinsic risk-return trade-off. These findings offer new insights into bacterial competition and may also facilitate the engineering of microbial consortia for synthetic biology applications.

Proteome variance differences within populations of European whitefish (Coregonus lavaretus) originating from contrasting salinity environments

Variation in gene expression is an important component of the phenotypic differences observed in nature. Gene expression variance across biological groups and environmental conditions has been studied extensively and has revealed specific genes and molecular mechanisms of interest. However, little is known regarding the importance of within-population gene expression variation to environmental adaptation. To address this issue, we quantified the proteomes of individuals of European whitefish (Coregonus lavaretus) from populations that have previously been shown to have adapted during early development to freshwater and brackishwater salinity environments. Using MS-based label-free proteomics, we studied 955 proteins in eight hatch-stage fish embryos from each population that had been reared in either freshwater or brackishwater salinity conditions. By comparing the levels of within-population protein expression variance over individuals and per protein between populations, we found that fish embryos from the population less affected by salinity level had also markedly higher levels of expression variance. Gene Ontologies and molecular pathways associated with osmoregulation showed the most significant difference of within-population proteome variance between populations. Several new candidate genes for salinity adaptation were identified, emphasising the added value of combining assessments of within-population gene expression variation with standard gene expression analysis practices for better understanding the mechanisms of environmental adaptation.

BIOLOGICAL SIGNIFICANCE

We demonstrate the benefits of studying within-population gene expression variance together with more typical methods of gene expression profiling. Proteome variance differences within European whitefish populations originating from different salinity environments allowed us to identify several new candidate genes for salinity adaptation in teleost fish and generate many further hypotheses to be tested. This article is part of a Special Issue entitled: Proteomics of non-model organisms.