Population diversification in a yeast metabolic program promotes anticipation of environmental shifts

The type of carbon source not the growth rate it supports can determine diauxie in Saccharomyces cerevisiae

How cells choose between carbon sources is a classic example of cellular decision-making. Microbes often prioritise glucose, but there has been little investigation of whether other sugars are also preferred. Here we study budding yeast growing on mixtures of sugars with palatinose, a sucrose isomer that cells catabolise with the MAL regulon. We find that the decision-making involves more than carbon flux-sensing: yeast prioritise galactose over palatinose, but sucrose and fructose weakly if at all despite each allowing faster growth than palatinose. With genetic perturbations and transcriptomics, we show that the regulation is active with repression of the MAL genes via Gal4, the GAL regulon's master regulator. We argue, using mathematical modelling, that cells enforce their preference for galactose through weakening the MAL regulon's positive feedback. They do so through decreasing intracellular palatinose by repressing MAL11, the palatinose transporter, and expressing the isomaltases IMA1 and IMA5. Supporting these predictions, we show that deleting IMA1 abolishes diauxie. Our results demonstrate that budding yeast actively prioritises carbon sources other than glucose and that such priorities need not reflect differences in growth rates. They imply that carbon-sensing strategies even in model organisms are more complex than previously thought.

A biological circuit to anticipate trend

Organisms gain by anticipating future changes in the environment. Those environmental changes often follow stochastic trends. The steeper the slope of the trend, the more likely the trend's momentum carries the future trend in the same direction. This article presents a simple biological circuit that measures the momentum, providing a prediction about future trend. The circuit calculates the momentum by the difference between a short-term and a long-term exponential moving average. The time lengths of the two moving averages can be adjusted by changing the decay rates of state variables. Different time lengths for those averages trade off between errors caused by noise and errors caused by lags in predicting a change in the direction of the trend. Prior studies have emphasized circuits that make similar calculations about trends. However, those prior studies embedded their analyses in the details of particular applications, obscuring the simple generality and wide applicability of the approach. The model here contributes to the topic by clarifying the great simplicity and generality of anticipation for stochastic trends. This article also notes that, in financial analysis, the difference between moving averages is widely used to predict future trends in asset prices. The financial measure is called the moving average convergence-divergence indicator. Connecting the biological problem to financial analysis opens the way for future studies in biology to exploit the variety of highly developed trend models in finance.

Shaping of microbial phenotypes by trade-offs

Growth rate maximization is an important fitness strategy for microbes. However, the wide distribution of slow-growing oligotrophic microbes in ecosystems suggests that rapid growth is often not favored across ecological environments. In many circumstances, there exist trade-offs between growth and other important traits (e.g., adaptability and survival) due to physiological and proteome constraints. Investments on alternative traits could compromise growth rate and microbes need to adopt bet-hedging strategies to improve fitness in fluctuating environments. Here we review the mechanistic role of trade-offs in controlling bacterial growth and further highlight its ecological implications in driving the emergences of many important ecological phenomena such as co-existence, population heterogeneity and oligotrophic/copiotrophic lifestyles.

Biodiversity is enhanced by sequential resource utilization and environmental fluctuations via emergent temporal niches

How natural communities maintain their remarkable biodiversity and which species survive in complex communities are central questions in ecology. Resource competition models successfully explain many phenomena but typically predict only as many species as resources can coexist. Here, we demonstrate that sequential resource utilization, or diauxie, with periodic growth cycles can support many more species than resources. We explore how communities modify their own environments by sequentially depleting resources to form sequences of temporal niches, or intermediately depleted environments. Biodiversity is enhanced when community-driven or environmental fluctuations modulate the resource depletion order and produce different temporal niches on each growth cycle. Community-driven fluctuations under constant environmental conditions are rare, but exploring them illuminates the temporal niche structure that emerges from sequential resource utilization. With environmental fluctuations, we find most communities have more stably coexisting species than resources with survivors accurately predicted by the same temporal niche structure and each following a distinct optimal strategy. Our results thus present a new niche-based approach to understanding highly diverse fluctuating communities.

Trade-offs between the instantaneous growth rate and long-term fitness: Consequences for microbial physiology and predictive computational models

Microbial systems biology has made enormous advances in relating microbial physiology to the underlying biochemistry and molecular biology. By meticulously studying model microorganisms, in particular Escherichia coli and Saccharomyces cerevisiae, increasingly comprehensive computational models predict metabolic fluxes, protein expression, and growth. The modeling rationale is that cells are constrained by a limited pool of resources that they allocate optimally to maximize fitness. As a consequence, the expression of particular proteins is at the expense of others, causing trade-offs between cellular objectives such as instantaneous growth, stress tolerance, and capacity to adapt to new environments. While current computational models are remarkably predictive for E. coli and S. cerevisiae when grown in laboratory environments, this may not hold for other growth conditions and other microorganisms. In this contribution, we therefore discuss the relationship between the instantaneous growth rate, limited resources, and long-term fitness. We discuss uses and limitations of current computational models, in particular for rapidly changing and adverse environments, and propose to classify microbial growth strategies based on Grimes's CSR framework.

Advances in automated and reactive flow cytometry for synthetic biotechnology

Automated flow cytometry (FC) has been initially considered for bioprocess monitoring and optimization. More recently, new physical and software interfaces have been made available, facilitating the access to this technology for labs and industries. It also comes with new capabilities, such as being able to act on the cultivation conditions based on population data. This approach, known as reactive FC, extended the range of applications of automated FC to bioprocess control and the stabilization of cocultures, but also to the broad field of synthetic and systems biology for the characterization of gene circuits. However, several issues must be addressed before automated and reactive FC can be considered standard and modular technologies.

Single-cell analysis of multiple invertible promoters reveals differential inversion rates as a strong determinant of bacterial population heterogeneity

Population heterogeneity can promote bacterial fitness in response to unpredictable environmental conditions. A major mechanism of phenotypic variability in the human gut symbiont Bacteroides spp. involves the inversion of promoters that drive the expression of capsular polysaccharides, which determine the architecture of the cell surface. High-throughput single-cell sequencing reveals substantial population heterogeneity generated through combinatorial promoter inversion regulated by a broadly conserved serine recombinase. Exploiting control over population diversification, we show that populations with different initial compositions converge to a similar composition over time. Combining our data with stochastic computational modeling, we demonstrate that the differential rates of promoter inversion are a major mechanism shaping population dynamics. More broadly, our approach could be used to interrogate single-cell combinatorial phase variable states of diverse microbes including bacterial pathogens.

Divergence of alternative sugar preferences through modulation of the expression and activity of the Gal3 sensor in yeast

Optimized nutrient utilization is crucial for the progression of microorganisms in competing communities. Here we investigate how different budding yeast species and ecological isolates have established divergent preferences for two alternative sugar substrates: Glucose, which is fermented preferentially by yeast, and galactose, which is alternatively used upon induction of the relevant GAL metabolic genes. We quantified the dose-dependent induction of the GAL1 gene encoding the central galactokinase enzyme and found that a very large diversification exists between different yeast ecotypes and species. The sensitivity of GAL1 induction correlates with the growth performance of the respective yeasts with the alternative sugar. We further define some of the mechanisms, which have established different glucose/galactose consumption strategies in representative yeast strains by modulating the activity of the Gal3 inducer. (1) Optimal galactose consumers, such as Saccharomyces uvarum, contain a hyperactive GAL3 promoter, sustaining highly sensitive GAL1 expression, which is not further improved upon repetitive galactose encounters. (2) Desensitized galactose consumers, such as S. cerevisiae Y12, contain a less sensitive Gal3 sensor, causing a shift of the galactose response towards higher sugar concentrations even in galactose experienced cells. (3) Galactose insensitive sugar consumers, such as S. cerevisiae DBVPG6044, contain an interrupted GAL3 gene, causing extremely reluctant galactose consumption, which is, however, improved upon repeated galactose availability. In summary, different yeast strains and natural isolates have evolved galactose utilization strategies, which cover the whole range of possible sensitivities by modulating the expression and/or activity of the inducible galactose sensor Gal3.

Meta-analysis of three-stressor combinations on population-level fitness reveal substantial higher-order interactions

Although natural populations are typically subjected to multiple stressors, most past research has focused on single-stressor and two-stressor interactions, with little attention paid to higher-order interactions among three or more stressors. However, higher-order interactions increasingly appear to be widespread. Consequently, we used a recently introduced and improved framework to re-analyze higher-order ecological interactions. We conducted a literature review of the last 100 years (1920-2020) and reanalyzed 142 ecological three-stressor interactions on species' populations from 38 published papers; the vast majority of these studies were from the past 10 years. We found that 95.8 % (n = 136) of the three-stressor combinations had either not been categorized before or resulted in different interactions than previously reported. We also found substantial levels of emergent properties-interactions that are not due to strong pairwise interactions within the combination but rather uniquely due to all three stressors being combined. Calculating net interactions-the overall accounting for all possible interactions within a combination including the emergent and all pairwise interactions-we found that the most prevalent interaction type is antagonism, corresponding to a smaller than expected effect based on single stressor effects. In contrast, for emergent interactions, the most prevalent interaction type is synergistic, resulting in a larger than expected effect based on single stressor effects. Additionally, we found that hidden suppressive interactions-where a pairwise interaction is suppressed by a third stressor-are found in the majority of combinations (74 %). Collectively, understanding multiple stressor interactions through applying an appropriate framework is crucial for answering fundamental questions in ecology and has implications for conservation biology and population management. Crucially, identifying emergent properties can reveal hidden suppressive interactions that could be particularly important for the ecological management of at-risk populations.

Recognition of galactose by a scaffold protein recruits a transcriptional activator for the GAL regulon induction in Candida albicans

The GAL pathway of yeasts has long served as a model system for understanding of how regulatory mode of eukaryotic metabolic pathways evolves. While Gal4 mode has been well-characterized in Saccharomycetaceae clade, little is known about the regulation of the GAL pathway in other yeasts. Here, we find that Rep1, a Ndt80-like family transcription factor, serves as a galactose sensor in the commensal-pathogenic fungus Candida albicans. It is presented at the GAL gene promoters independent of the presence of galactose. Rep1 recognizes galactose via a direct physical interaction. The net result of this interaction is the recruitment of a transcriptional activator Cga1 (Candida galactose gene activator, orf19.4959) and transcription of the GAL genes proceeds. Rep1 and Cga1 are conserved across the CTG species. Rep1 itself does not possess transcriptional activity. Instead, it provides a scaffold to recruit different factors for transcriptional regulation. Rep1-Cga1 mode of regulation represents a new example of network rewiring in fungi, which provides insight into how C. albicans evolves transcriptional programs to colonize diverse host niches.

Spatiotemporal Patterning Enabled by Gene Regulatory Networks

Spatiotemporal pattern formation plays a key role in various biological phenomena including embryogenesis and neural network formation. Though the reaction-diffusion systems enabling pattern formation have been studied phenomenologically, the biomolecular mechanisms behind these processes have not been modeled in detail. Here, we study the emergence of spatiotemporal patterns due to simple, synthetic and commonly observed two- and three-node gene regulatory network motifs coupled with their molecular diffusion in one- and two-dimensional space. We investigate the patterns formed due to the coupling of inherent multistable and oscillatory behavior of the toggle switch, toggle switch with double self-activation, toggle triad, and repressilator with the effect of spatial diffusion of these molecules. We probe multiple parameter regimes corresponding to different regions of stability (monostable, multistable, oscillatory) and assess the impact of varying diffusion coefficients. This analysis offers valuable insights into the design principles of pattern formation facilitated by these network motifs, and it suggests the mechanistic underpinnings of biological pattern formation.

Stringent response ensures the timely adaptation of bacterial growth to nutrient downshift

Timely adaptation to nutrient downshift is crucial for bacteria to maintain fitness during feast and famine cycle in the natural niche. However, the molecular mechanism that ensures the timely adaption of bacterial growth to nutrient downshift remains poorly understood. Here, we quantitatively investigated the adaptation of Escherichia coli to various kinds of nutrient downshift. We found that relA deficient strain, which is devoid of stringent response, exhibits a significantly longer growth lag than wild type strain during adapting to both amino acid downshift and carbon downshift. Quantitative proteomics show that increased (p)ppGpp level promotes the growth adaption of bacteria to amino acid downshift via triggering the proteome resource re-allocation from ribosome synthesis to amino acid biosynthesis. Such type of proteome re-allocation is significantly delayed in the relA-deficient strain, which underlies its longer lag than wild type strain during amino acid downshift. During carbon downshift, a lack of stringent response in relA deficient strain leads to disruption of the transcription-translation coordination, thus compromising the transcription processivity and further the timely expression of related catabolic operons for utilizing secondary carbon sources. Our studies shed light on the fundamental strategy of bacteria to maintain fitness under nutrient-fluctuating environments.

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.

Escherichia coli robustly expresses ATP synthase at growth rate‐maximizing concentrations

Fitness‐enhancing adaptations of protein expression and its regulation are an important aspect of bacterial evolution. A key question is whether evolution has led to optimal protein expression that maximizes immediate growth rate (short‐term fitness) in a robust manner (consistently across diverse conditions). Alternatively, they could display suboptimal short‐term fitness, because they cannot do better or because they instead strive for long‐term fitness maximization by, for instance, preparing for future conditions. To address this question, we focus on the ATP‐producing enzyme F₁F₀ H⁺‐ATPase, which is an abundant enzyme and ubiquitously expressed across conditions. Its expression is highly regulated and dependent on growth rate and nutrient conditions. For instance, during growth on sugars, when metabolism is overflowing acetate, glycolysis supplies most ATP, while H⁺‐ATPase is the main source of ATP synthesis during growth on acetate. We tested the optimality of H⁺‐ATPase expression in Escherichia coli across different nutrient conditions. In all tested conditions, wild‐type E. coli expresses its H⁺‐ATPase remarkably close (within a few per cent) to optimal concentrations that maximize immediate growth rate. This work indicates that bacteria can indeed achieve robust optimal protein expression for immediate growth‐rate maximization.

Diauxic lags explain unexpected coexistence in multi-resource environments

How the coexistence of species is affected by the presence of multiple resources is a major question in microbial ecology. We experimentally demonstrate that differences in diauxic lags, which occur as species deplete their own environments and adapt their metabolisms, allow slow-growing microbes to stably coexist with faster-growing species in multi-resource environments despite being excluded in single-resource environments. In our focal example, an Acinetobacter species (Aci2) competitively excludes Pseudomonas aurantiaca (Pa) on alanine and on glutamate. However, they coexist on the combination of both resources. Experiments reveal that Aci2 grows faster but Pa has shorter diauxic lags. We establish a tradeoff between Aci2's fast growth and Pa's short lags as their mechanism for coexistence. We model this tradeoff to accurately predict how environmental changes affect community composition. We extend our work by surveying a large set of competitions and observe coexistence nearly four times as frequently when the slow-grower is the fast-switcher. Our work illustrates a simple mechanism, based entirely on supplied-resource growth dynamics, for the emergence of multi-resource coexistence.

Identification of genes related to hydrolysis and assimilation of Agave fructans in Candida apicola NRRL Y-50540 and Torulaspora delbrueckii NRRL Y-50541 by de novo transcriptome analysis

Fructans are the main sugar in agave pine used by yeasts during mezcal fermentation processes, from which Candida apicola NRRL Y-50540 and Torulaspora delbrueckii NRRL Y-50541 were isolated. De novo transcriptome analysis was carried out to identify genes involved in the hydrolysis and assimilation of Agave fructans (AF). We identified a transcript annotated as SUC2, which is related to β-fructofuranosidase activity, and several differential expressed genes involved in the transcriptional regulation of SUC2 such as: MIG1, MTH1, SNF1, SNF5, REG1, SSN6, SIP1, SIP2, SIP5, GPR1, RAS2, and PKA. Some of these genes were specifically expressed in some of the yeasts according to their fructans assimilation metabolism. Different hexose transporters that could be related to the assimilation of fructose and glucose were found in both the transcriptomes. Our findings provide a better understanding of AF assimilation in these yeasts and provide resources for further metabolic engineering and biotechnology applications.

Polysaccharide utilization loci in Bacteroides determine population fitness and community-level interactions

Polysaccharide utilization loci (PULs) are co-regulated bacterial genes that sense nutrients and enable glycan digestion. Human gut microbiome members, notably Bacteroides, contain numerous PULs that enable glycan utilization and shape ecological dynamics. To investigate the role of PULs on fitness and inter-species interactions, we develop a CRISPR-based genome editing tool to study 23 PULs in Bacteroides uniformis (BU). BU PULs show distinct glycan-degrading functions and transcriptional coordination that enables the population to adapt upon loss of other PULs. Exploiting a BU mutant barcoding strategy, we demonstrate that in vitro fitness and BU colonization in the murine gut are enhanced by deletion of specific PULs and modulated by glycan availability. PULs mediate glycan-dependent interactions with butyrate producers that depend on the degradation mechanism and glycan utilization ability of the butyrate producer. Thus, PULs determine community dynamics and butyrate production and provide a selective advantage or disadvantage depending on the nutritional landscape.

Robustness: linking strain design to viable bioprocesses

Microbial cell factories are becoming increasingly popular for the sustainable production of various chemicals. Metabolic engineering has led to the design of advanced cell factories; however, their long-term yield, titer, and productivity falter when scaled up and subjected to industrial conditions. This limitation arises from a lack of robustness - the ability to maintain a constant phenotype despite the perturbations of such processes. This review describes predictable and stochastic industrial perturbations as well as state-of-the-art technologies to counter process variability. Moreover, we distinguish robustness from tolerance and discuss the potential of single-cell studies for improving system robustness. Finally, we highlight ways of achieving consistent and comparable quantification of robustness that can guide the selection of strains for industrial bioprocesses.

Experimental Evolution of Anticipatory Regulation in Escherichia coli

Environmental cues in an ecological niche are often temporal in nature. For instance, in temperate climates, temperature is higher in daytime compared to during night. In response to these temporal cues, bacteria have been known to exhibit anticipatory regulation, whereby triggering response to a yet to appear cue. Such an anticipatory response in known to enhance Darwinian fitness, and hence, is likely an important feature of regulatory networks in microorganisms. However, the conditions under which an anticipatory response evolves as an adaptive response are not known. In this work, we develop a quantitative model to study response of a population to two temporal environmental cues, and predict variables which are likely important for evolution of anticipatory regulatory response. We follow this with experimental evolution of Escherichia coli in alternating environments of rhamnose and paraquat for ∼850 generations. We demonstrate that growth in this cyclical environment leads to evolution of anticipatory regulation. As a result, pre-exposure to rhamnose leads to a greater fitness in paraquat environment. Genome sequencing reveals that this anticipatory regulation is encoded via mutations in global regulators. Overall, our study contributes to understanding of how environment shapes the topology of regulatory networks in an organism.

Suboptimal resource allocation in changing environments constrains response and growth in bacteria

To respond to fluctuating conditions, microbes typically need to synthesize novel proteins. As this synthesis relies on sufficient biosynthetic precursors, microbes must devise effective response strategies to manage depleting precursors. To better understand these strategies, we investigate the active response of Escherichia coli to changes in nutrient conditions, connecting transient gene expression to growth phenotypes. By synthetically modifying gene expression during changing conditions, we show how the competition by genes for the limited protein synthesis capacity constrains cellular response. Despite this constraint cells substantially express genes that are not required, trapping them in states where precursor levels are low and the genes needed to replenish the precursors are outcompeted. Contrary to common modeling assumptions, our findings highlight that cells do not optimize growth under changing environments but rather exhibit hardwired response strategies that may have evolved to promote fitness in their native environment. The constraint and the suboptimality of the cellular response uncovered provide a conceptual framework relevant for many research applications, from the prediction of evolution to the improvement of gene circuits in biotechnology.

Computational analysis of GAL pathway pinpoints mechanisms underlying natural variation

Quantitative traits are measurable phenotypes that show continuous variation over a wide phenotypic range. Enormous effort has recently been put into determining the genetic influences on a variety of quantitative traits with mixed success. We identified a quantitative trait in a tractable model system, the GAL pathway in yeast, which controls the uptake and metabolism of the sugar galactose. GAL pathway activation depends both on galactose concentration and on the concentrations of competing, preferred sugars such as glucose. Natural yeast isolates show substantial variation in the behavior of the pathway. All studied yeast strains exhibit bimodal responses relative to external galactose concentration, i.e. a set of galactose concentrations existed at which both GAL-induced and GAL-repressed subpopulations were observed. However, these concentrations differed in different strains. We built a mechanistic model of the GAL pathway and identified parameters that are plausible candidates for capturing the phenotypic features of a set of strains including standard lab strains, natural variants, and mutants. In silico perturbation of these parameters identified variation in the intracellular galactose sensor, Gal3p, the negative feedback node within the GAL regulatory network, Gal80p, and the hexose transporters, HXT, as the main sources of the bimodal range variation. We were able to switch the phenotype of individual yeast strains in silico by tuning parameters related to these three elements. Determining the basis for these behavioral differences may give insight into how the GAL pathway processes information, and into the evolution of nutrient metabolism preferences in different strains. More generally, our method of identifying the key parameters that explain phenotypic variation in this system should be generally applicable to other quantitative traits.

Microbes adopt elaborate strategies for the preferred uptake and use of nutrients to cope with complex and fluctuating environments. As a result, yeast strains originating from different ecological niches show significant variation in the way they induce genes in the galactose metabolism (GAL) pathway in response to nutrient signals. To identify the mechanistic sources of this variation, we built a mathematical model to simulate the dynamics of the galactose metabolic regulation network, and studied how parameters with different biological implications contributed to the natural variation. We found that variations in the behavior of the galactose sensor Gal3p, the negative feedback node Gal80p, and the hexose transporters HXT were critical elements in the GAL pathway response. Tuning single parameters in silico was sufficient to achieve phenotype switching between different yeast strains. Our computational approach should be generally useful to help pinpoint the genetic and molecular bases of natural variation in other systems.

Cis-regulatory variants affect gene expression dynamics in yeast

Evolution of cis-regulatory sequences depends on how they affect gene expression and motivates both the identification and prediction of cis-regulatory variants responsible for expression differences within and between species. While much progress has been made in relating cis-regulatory variants to expression levels, the timing of gene activation and repression may also be important to the evolution of cis-regulatory sequences. We investigated allele-specific expression (ASE) dynamics within and between Saccharomyces species during the diauxic shift and found appreciable cis-acting variation in gene expression dynamics. Within-species ASE is associated with intergenic variants, and ASE dynamics are more strongly associated with insertions and deletions than ASE levels. To refine these associations, we used a high-throughput reporter assay to test promoter regions and individual variants. Within the subset of regions that recapitulated endogenous expression, we identified and characterized cis-regulatory variants that affect expression dynamics. Between species, chimeric promoter regions generate novel patterns and indicate constraints on the evolution of gene expression dynamics. We conclude that changes in cis-regulatory sequences can tune gene expression dynamics and that the interplay between expression dynamics and other aspects of expression is relevant to the evolution of cis-regulatory sequences.

Variation in the modality of a yeast signaling pathway is mediated by a single regulator

Bimodal gene expression by genetically identical cells is a pervasive feature of signaling networks and has been suggested to allow organisms to hedge their ‘bets’ in uncertain conditions. In the galactose-utilization (GAL) pathway of Saccharomyces cerevisiae, gene induction is unimodal or bimodal depending on natural genetic variation and pre-induction conditions. Here, we find that this variation in modality arises from regulation of two features of the pathway response: the fraction of cells that show induction and their level of expression. GAL3, the galactose sensor, controls the fraction of induced cells, and titrating its expression is sufficient to control modality; moreover, all the observed differences in modality between different pre-induction conditions and among natural isolates can be explained by changes in GAL3’s regulation and activity. The ability to switch modality by tuning the activity of a single protein may allow rapid adaptation of bet hedging to maximize fitness in complex environments.

Gene Regulation and Cellular Metabolism: An Essential Partnership

It is recognized that cell metabolism is tightly connected to other cellular processes such as regulation of gene expression. Metabolic pathways not only provide the precursor molecules necessary for gene expression, but they also provide ATP, the primary fuel driving gene expression. However, metabolic conditions are highly variable since nutrient uptake is not a uniform process. Thus, cells must continually calibrate gene expression to their changing metabolite and energy budgets. This review discusses recent advances in understanding the molecular and biophysical mechanisms that connect metabolism and gene regulation as cells navigate their growth, proliferation, and differentiation. Particular focus is given to these mechanisms in the context of organismal development.

Variation, Variegation and Heritable Gene Repression in S. cerevisiae

Phenotypic heterogeneity provides growth advantages for a population upon changes of the environment. In S. cerevisiae, such heterogeneity has been observed as "on/off" states in the expression of individual genes in individual cells. These variations can persist for a limited or extended number of mitotic divisions. Such traits are known to be mediated by heritable chromatin structures, by the mitotic transmission of transcription factors involved in gene regulatory circuits or by the cytoplasmic partition of prions or other unstructured proteins. The significance of such epigenetic diversity is obvious, however, we have limited insight into the mechanisms that generate it. In this review, we summarize the current knowledge of epigenetically maintained heterogeneity of gene expression and point out similarities and converging points between different mechanisms. We discuss how the sharing of limiting repression or activation factors can contribute to cell-to-cell variations in gene expression and to the coordination between short- and long- term epigenetic strategies. Finally, we discuss the implications of such variations and strategies in adaptation and aging.

Insights on the Control of Yeast Single-Cell Growth Variability by Members of the Trehalose Phosphate Synthase (TPS) Complex

Single-cell variability of growth is a biological phenomenon that has attracted growing interest in recent years. Important progress has been made in the knowledge of the origin of cell-to-cell heterogeneity of growth, especially in microbial cells. To better understand the origins of such heterogeneity at the single-cell level, we developed a new methodological pipeline that coupled cytometry-based cell sorting with automatized microscopy and image analysis to score the growth rate of thousands of single cells. This allowed investigating the influence of the initial amount of proteins of interest on the subsequent growth of the microcolony. As a preliminary step to validate this experimental setup, we referred to previous findings in yeast where the expression level of Tsl1, a member of the Trehalose Phosphate Synthase (TPS) complex, negatively correlated with cell division rate. We unfortunately could not find any influence of the initial TSL1 expression level on the growth rate of the microcolonies. We also analyzed the effect of the natural variations of trehalose-6-phosphate synthase (TPS1) expression on cell-to-cell growth heterogeneity, but we did not find any correlation. However, due to the already known altered growth of the tps1Δ mutants, we tested this strain at the single-cell level on a permissive carbon source. This mutant showed an outstanding lack of reproducibility of growth rate distributions as compared to the wild-type strain, with variable proportions of non-growing cells between cultivations and more heterogeneous microcolonies in terms of individual growth rates. Interestingly, this variable behavior at the single-cell level was reminiscent to the high variability that is also stochastically suffered at the population level when cultivating this tps1Δ strain, even when using controlled bioreactors.

Slow Adaptive Response of Budding Yeast Cells to Stable Conditions of Continuous Culture Can Occur without Genome Modifications

Continuous cultures assure the invariability of environmental conditions and the metabolic state of cultured microorganisms, whereas batch-cultured cells undergo constant changes in nutrients availability. For that reason, continuous culture is sometimes employed in the whole transcriptome, whole proteome, or whole metabolome studies. However, the typical method for establishing uniform growth of a cell population, i.e., by limited chemostat, results in the enrichment of the cell population gene pool with mutations adaptive for starvation conditions. These adaptive changes can skew the results of large-scale studies. It is commonly assumed that these adaptations reflect changes in the genome, and this assumption has been confirmed experimentally in rare cases. Here we show that in a population of budding yeast cells grown for over 200 generations in continuous culture in non-limiting minimal medium and therefore not subject to selection pressure, remodeling of transcriptome occurs, but not as a result of the accumulation of adaptive mutations. The observed changes indicate a shift in the metabolic balance towards catabolism, a decrease in ribosome biogenesis, a decrease in general stress alertness, reorganization of the cell wall, and transactions occurring at the cell periphery. These adaptive changes signify the acquisition of a new lifestyle in a stable nonstressful environment. The absence of underlying adaptive mutations suggests these changes may be regulated by another mechanism.

Microbial metabolic noise

From the time a cell was first placed under the microscope, it became apparent that identifying two clonal cells that "look" identical is extremely challenging. Since then, cell-to-cell differences in shape, size, and protein content have been carefully examined, informing us of the ultimate limits that hinder two cells from occupying an identical phenotypic state. Here, we present recent experimental and computational evidence that similar limits emerge also in cellular metabolism. These limits pertain to stochastic metabolic dynamics and, thus, cell-to-cell metabolic variability, including the resulting adapting benefits. We review these phenomena with a focus on microbial metabolism and conclude with a brief outlook on the potential relationship between metabolic noise and adaptive evolution. This article is categorized under: Metabolic Diseases > Computational Models Metabolic Diseases > Biomedical Engineering.

A Putative Bet-Hedging Strategy Buffers Budding Yeast against Environmental Instability

To grow and divide, cells must extract resources from dynamic and unpredictable environments. Many organisms use different metabolic strategies for distinct contexts. Budding yeast can produce ATP from carbon sources by mechanisms that prioritize either speed (fermentation) or yield (respiration). Withdrawing glucose from exponentially growing cells reveals variability in their ability to switch from fermentation to respiration. We observe two subpopulations of glucose-starved cells: recoverers, which rapidly adapt and resume growth, and arresters, which enter a shock state characterized by deformation of many cellular structures, including mitochondria. These states are heritable, and on high glucose, arresters grow and divide faster than recoverers. Recoverers have a fitness advantage during a carbon source shift but are less fit in a constant, high-glucose environment, and we observe natural variation in the frequency of the two states across wild yeast strains. These experiments suggest that bet hedging has evolved in budding yeast.

A Similar Speciation Process Relying on Cellular Stochasticity in Microbial and Cancer Cell Populations

Similarities between microbial and cancer cells were noticed in recent years and serve as a basis for an atavism theory of cancer. Cancer cells would rely on the reactivation of an ancestral "genetic program" that would have been repressed in metazoan cells. Here we argue that cancer cells resemble unicellular organisms mainly in their similar way to exploit cellular stochasticity to produce cell specialization and maximize proliferation. Indeed, the relationship between low stochasticity, specialization, and quiescence found in normal differentiated metazoan cells is lost in cancer. On the contrary, low stochasticity and specialization are associated with high proliferation among cancer cells, as it is observed for the "specialist" cells in microbial populations that fully exploit nutritional resources to maximize proliferation. Thus, we propose a model where the appearance of cancer phenotypes can be solely due to an adaptation and a speciation process based on initial increase in cellular stochasticity.

Wide lag time distributions break a trade-off between reproduction and survival in bacteria

Many microorganisms face a fundamental trade-off between reproduction and survival: Rapid growth boosts population size but makes microorganisms sensitive to external stressors. Here, we show that starved bacteria encountering new resources can break this trade-off by evolving phenotypic heterogeneity in lag time. We quantify the distribution of single-cell lag times of populations of starved Escherichia coli and show that population growth after starvation is primarily determined by the cells with shortest lag due to the exponential nature of bacterial population dynamics. As a consequence, cells with long lag times have no substantial effect on population growth resumption. However, we observe that these cells provide tolerance to stressors such as antibiotics. This allows an isogenic population to break the trade-off between reproduction and survival. We support this argument with an evolutionary model which shows that bacteria evolve wide lag time distributions when both rapid growth resumption and survival under stressful conditions are under selection. Our results can explain the prevalence of antibiotic tolerance by lag and demonstrate that the benefits of phenotypic heterogeneity in fluctuating environments are particularly high when minorities with extreme phenotypes dominate population dynamics.

Direct Observation of the Dynamics of Single-Cell Metabolic Activity during Microbial Diauxic Growth

Population-level analyses are rapidly becoming inadequate to answer many of biomedical science and microbial ecology's most pressing questions. The role of microbial populations within ecosystems and the evolutionary selective pressure on individuals depend fundamentally on the metabolic activity of single cells. Yet, many existing single-cell technologies provide only indirect evidence of metabolic specialization because they rely on correlations between transcription and phenotype established at the level of the population to infer activity. In this study, we take a top-down approach using isotope labels and secondary ion mass spectrometry to track the uptake of carbon and nitrogen atoms from different sources into biomass and directly observe dynamic changes in anabolic specialization at the level of single cells. We investigate the classic microbiological phenomenon of diauxic growth at the single-cell level in the model methylotroph Methylobacterium extorquens In nature, this organism inhabits the phyllosphere, where it experiences diurnal changes in the available carbon substrates, necessitating an overhaul of central carbon metabolism. We show that the population exhibits a unimodal response to the changing availability of viable substrates, a conclusion that supports the canonical model but has thus far been supported by only indirect evidence. We anticipate that the ability to monitor the dynamics of anabolism in individual cells directly will have important applications across the fields of ecology, medicine, and biogeochemistry, especially where regulation downstream of transcription has the potential to manifest as heterogeneity that would be undetectable with other existing single-cell approaches.IMPORTANCE Understanding how genetic information is realized as the behavior of individual cells is a long-term goal of biology but represents a significant technological challenge. In clonal microbial populations, variation in gene regulation is often interpreted as metabolic heterogeneity. This follows the central dogma of biology, in which information flows from DNA to RNA to protein and ultimately manifests as activity. At present, DNA and RNA can be characterized in single cells, but the abundance and activity of proteins cannot. Inferences about metabolic activity usually therefore rely on the assumption that transcription reflects activity. By tracking the atoms from which they build their biomass, we make direct observations of growth rate and substrate specialization in individual cells throughout a period of growth in a changing environment. This approach allows the flow of information from DNA to be constrained from the distal end of the regulatory cascade and will become an essential tool in the rapidly advancing field of single-cell metabolism.

Microbial phenotypic heterogeneity in response to a metabolic toxin: Continuous, dynamically shifting distribution of formaldehyde tolerance in Methylobacterium extorquens populations

While microbiologists often make the simplifying assumption that genotype determines phenotype in a given environment, it is becoming increasingly apparent that phenotypic heterogeneity (in which one genotype generates multiple phenotypes simultaneously even in a uniform environment) is common in many microbial populations. The importance of phenotypic heterogeneity has been demonstrated in a number of model systems involving binary phenotypic states (e.g., growth/non-growth); however, less is known about systems involving phenotype distributions that are continuous across an environmental gradient, and how those distributions change when the environment changes. Here, we describe a novel instance of phenotypic diversity in tolerance to a metabolic toxin within wild-type populations of Methylobacterium extorquens, a ubiquitous phyllosphere methylotroph capable of growing on the methanol periodically released from plant leaves. The first intermediate in methanol metabolism is formaldehyde, a potent cellular toxin that is lethal in high concentrations. We have found that at moderate concentrations, formaldehyde tolerance in M. extorquens is heterogeneous, with a cell's minimum tolerance level ranging between 0 mM and 8 mM. Tolerant cells have a distinct gene expression profile from non-tolerant cells. This form of heterogeneity is continuous in terms of threshold (the formaldehyde concentration where growth ceases), yet binary in outcome (at a given formaldehyde concentration, cells either grow normally or die, with no intermediate phenotype), and it is not associated with any detectable genetic mutations. Moreover, tolerance distributions within the population are dynamic, changing over time in response to growth conditions. We characterized this phenomenon using bulk liquid culture experiments, colony growth tracking, flow cytometry, single-cell time-lapse microscopy, transcriptomics, and genome resequencing. Finally, we used mathematical modeling to better understand the processes by which cells change phenotype, and found evidence for both stochastic, bidirectional phenotypic diversification and responsive, directed phenotypic shifts, depending on the growth substrate and the presence of toxin.

Scientists tend to appreciate microbes for their simplicity and predictability: a population of genetically identical cells inhabiting a uniform environment is expected to behave in a uniform way. However, counter-examples to this assumption are frequently being discovered, forcing a re-examination of the relationship between genotype and phenotype. In most such examples, bacterial cells are found to split into two discrete populations, for instance growing and non-growing. Here, we report the discovery of a novel example of microbial phenotypic heterogeneity in which cells are distributed along a gradient of phenotypes, ranging from low to high tolerance of a toxic chemical. Furthermore, we demonstrate that the distribution of phenotypes changes in different growth conditions, and we use mathematical modeling to show that cells may change their phenotype either randomly or in a particular direction in response to the environment. Our work expands our understanding of how a bacterial cell's genome, family history, and environment all contribute to its behavior, with implications for the diverse situations in which we care to understand the growth of any single-celled populations.

The ecology and evolution of microbial metabolic strategies

Free-living microbes are generally capable of growing on multiple different nutrients. Some of those nutrients are used simultaneously, while others are used sequentially. The pattern of nutrient preferences and co-utilization defines the metabolic strategy of a microorganism. Metabolic strategies can substantially affect ecological interactions between species, but their evolution and distribution across the tree of life remain poorly characterized. We discuss how the confluence of better computational models of genotype-phenotype maps and high-throughput experimental tools can help us fill gaps in our knowledge and incorporate metabolic strategies into quantitative predictive models of microbial consortia.

MMHelper: An automated framework for the analysis of microscopy images acquired with the mother machine

Live-cell imaging in microfluidic devices now allows the investigation of cellular heterogeneity within microbial populations. In particular, the mother machine technology developed by Wang et al. has been widely employed to investigate single-cell physiological parameters including gene expression, growth rate, mutagenesis, and response to antibiotics. One of the advantages of the mother machine technology is the ability to generate vast amounts of images; however, the time consuming analysis of these images constitutes a severe bottleneck. Here we overcome this limitation by introducing MMHelper ( https://doi.org/10.5281/zenodo.3254394 ), a publicly available custom software implemented in Python which allows the automated analysis of brightfield or phase contrast, and any associated fluorescence, images of bacteria confined in the mother machine. We show that cell data extracted via MMHelper from tens of thousands of individual cells imaged in brightfield are consistent with results obtained via semi-automated image analysis based on ImageJ. Furthermore, we benchmark our software capability in processing phase contrast images from other laboratories against other publicly available software. We demonstrate that MMHelper has over 90% detection efficiency for brightfield and phase contrast images and provides a new open-source platform for the extraction of single-bacterium data, including cell length, area, and fluorescence intensity.

Evolution: The Genetics of Milk Champagne Yeast

A newly characterized group of budding yeast found in fermented milk drinks in West and Central Asia far outshines its relatives in the ability to use galactose, a milk sugar byproduct.

Scalable nonlinear programming framework for parameter estimation in dynamic biological system models

We present a nonlinear programming (NLP) framework for the scalable solution of parameter estimation problems that arise in dynamic modeling of biological systems. Such problems are computationally challenging because they often involve highly nonlinear and stiff differential equations as well as many experimental data sets and parameters. The proposed framework uses cutting-edge modeling and solution tools which are computationally efficient, robust, and easy-to-use. Specifically, our framework uses a time discretization approach that: i) avoids repetitive simulations of the dynamic model, ii) enables fully algebraic model implementations and computation of derivatives, and iii) enables the use of computationally efficient nonlinear interior point solvers that exploit sparse and structured linear algebra techniques. We demonstrate these capabilities by solving estimation problems for synthetic human gut microbiome community models. We show that an instance with 156 parameters, 144 differential equations, and 1,704 experimental data points can be solved in less than 3 minutes using our proposed framework (while an off-the-shelf simulation-based solution framework requires over 7 hours). We also create large instances to show that the proposed framework is scalable and can solve problems with up to 2,352 parameters, 2,304 differential equations, and 20,352 data points in less than 15 minutes. The proposed framework is flexible and easy-to-use, can be broadly applied to dynamic models of biological systems, and enables the implementation of sophisticated estimation techniques to quantify parameter uncertainty, to diagnose observability/uniqueness issues, to perform model selection, and to handle outliers.

Constructing and validating dynamic models of biological systems spanning biomolecular networks to ecological systems is a challenging problem. Here we present a scalable computational framework to rapidly infer parameters in complex dynamic models of biological systems from large-scale experimental data. The framework was applied to infer parameters of a synthetic microbial community model from large-scale time series data. We also demonstrate that this framework can be used to analyze parameter uncertainty, to diagnose whether the experimental data are sufficient to uniquely determine the parameters, to determine the model that best describes the data, and to infer parameters in the face of data outliers.

Growth strategy of microbes on mixed carbon sources

A classic problem in microbiology is that bacteria display two types of growth behavior when cultured on a mixture of two carbon sources: the two sources are sequentially consumed one after another (diauxie) or they are simultaneously consumed (co-utilization). The search for the molecular mechanism of diauxie led to the discovery of the lac operon. However, questions remain as why microbes would bother to have different strategies of taking up nutrients. Here we show that diauxie versus co-utilization can be understood from the topological features of the metabolic network. A model of optimal allocation of protein resources quantitatively explains why and how the cell makes the choice. In case of co-utilization, the model predicts the percentage of each carbon source in supplying the amino acid pools, which is quantitatively verified by experiments. Our work solves a long-standing puzzle and provides a quantitative framework for the carbon source utilization of microbes.

Bacteria grown on two carbon sources either consume both sources simultaneously or consume them sequentially. Here the authors use a metabolic network model of E. coli to show that optimal protein resource allocation and topological features of the network can explain the choice of carbon acquisition.

Sensitive high-throughput single-cell RNA-seq reveals within-clonal transcript correlations in yeast populations

Single-cell RNA-seq has revealed extensive cellular heterogeneity within many organisms, but few methods have been developed for microbial clonal populations. The yeast genome displays unusually dense transcript spacing, with interleaved and overlapping transcription from both strands, resulting in a minuscule but complex pool of RNA protected by a resilient cell wall. Here, we have developed a sensitive, scalable, and inexpensive yeast single-cell RNA-seq (yscRNA-seq) method that digitally counts transcript start sites in a strand- and isoform-specific manner. YscRNA-seq detects the expression of low-abundant, non-coding RNAs, and at least half of the protein-coding genome in each cell. Within clonal cells, we observed a negative correlation for the expression of sense/antisense pairs, while paralogs and divergent transcripts co-express. Combining yscRNA-seq with index sorting, we uncovered a linear relationship between cell size and RNA content. Although we detected an average of ~3.5 molecules/gene, the number of expressed isoforms are restricted at the single-cell level. Remarkably, the expression of metabolic genes is highly variable, while their stochastic expression primes cells for increased fitness towards the corresponding environmental challenge. These findings suggest that functional transcript diversity acts as a mechanism for providing a selective advantage to individual cells within otherwise transcriptionally heterogeneous populations.

The Crabtree Effect Shapes the Saccharomyces cerevisiae Lag Phase during the Switch between Different Carbon Sources

When faced with environmental changes, microbes often enter a temporary growth arrest during which they reprogram the expression of specific genes to adapt to the new conditions. A prime example of such a lag phase occurs when microbes need to switch from glucose to other, less-preferred carbon sources. Despite its industrial relevance, the genetic network that determines the duration of the lag phase has not been studied in much detail. Here, we performed a genome-wide Bar-Seq screen to identify genetic determinants of the Saccharomyces cerevisiae glucose-to-galactose lag phase. The results show that genes involved in respiration, and specifically those encoding complexes III and IV of the electron transport chain, are needed for efficient growth resumption after the lag phase. Anaerobic growth experiments confirmed the importance of respiratory energy conversion in determining the lag phase duration. Moreover, overexpression of the central regulator of respiration, HAP4, leads to significantly shorter lag phases. Together, these results suggest that the glucose-induced repression of respiration, known as the Crabtree effect, is a major determinant of microbial fitness in fluctuating carbon environments.IMPORTANCE The lag phase is arguably one of the prime characteristics of microbial growth. Longer lag phases result in lower competitive fitness in variable environments, and the duration of the lag phase is also important in many industrial processes where long lag phases lead to sluggish, less efficient fermentations. Despite the immense importance of the lag phase, surprisingly little is known about the exact molecular processes that determine its duration. Our study uses the molecular toolbox of S. cerevisiae combined with detailed growth experiments to reveal how the transition from fermentative to respirative metabolism is a key bottleneck for cells to overcome the lag phase. Together, our findings not only yield insight into the key molecular processes and genes that influence lag duration but also open routes to increase the efficiency of industrial fermentations and offer an experimental framework to study other types of lag behavior.

An ensemble of mathematical models showing diauxic growth behaviour

Background

Carbon catabolite repression (CCR) controls the order in which different carbon sources are metabolised. Although this system is one of the paradigms of regulation in bacteria, the underlying mechanisms remain controversial. CCR involves the coordination of different subsystems of the cell - responsible for the uptake of carbon sources, their breakdown for the production of energy and precursors, and the conversion of the latter to biomass. The complexity of this integrated system, with regulatory mechanisms cutting across metabolism, gene expression, and signalling, has motivated important modelling efforts over the past four decades, especially in the enterobacterium Escherichia coli.

Results

Starting from a simple core model with only four intracellular metabolites, we develop an ensemble of model variants, all showing diauxic growth behaviour during a batch process. The model variants fall into one of the four categories: flux balance models, kinetic models with growth dilution, kinetic models with regulation, and resource allocation models. The model variants differ from one another in only a single aspect, each breaking the symmetry between the two substrate assimilation pathways in a different manner, and can be quantitatively compared using a so-called diauxic growth index. For each of the model variants, we predict the behaviour in two new experimental conditions, namely a glucose pulse for a culture growing in minimal medium with lactose and a batch culture with different initial concentrations of the components of the transport systems. When qualitatively comparing these predictions with experimental data for these two conditions, a number of models can be excluded while other model variants are still not discriminable. The best-performing model variants are based on inducer inclusion and activation of enzymatic genes by a global transcription factor, but the other proposed factors may complement these well-known regulatory mechanisms.

Conclusions

The model ensemble presented here offers a better understanding of the variety of mechanisms that have been proposed to play a role in CCR. In addition, it provides an educational resource for systems biology, as it gives an introduction to a broad range of modeling approaches in the context of a simple but biologically relevant example.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0604-8) contains supplementary material, which is available to authorized users.

A GRX1 Promoter Variant Confers Constitutive Noisy Bimodal Expression That Increases Oxidative Stress Resistance in Yeast

Higher noise in the expression of stress-related genes was previously shown to confer better resistance in selective conditions. Thus, evolving the promoter of such genes toward higher transcriptional noise appears to be an attractive strategy to engineer microbial strains with enhanced stress resistance. Here we generated hundreds of promoter variants of the GRX1 gene involved in oxidative stress resistance in Saccharomyces cerevisiae and created a yeast library by replacing the native GRX1 promoter by these variants at the native locus. An outlier clone with very strong increase in noise (6-times) at the same mean expression level as the native strain was identified whereas the other noisiest clones were only 3-times increased. This variant provides constitutive bimodal expression and consists in 3 repeated but differently mutated copies of the GRX1 promoter. In spite of the multi-factorial oxidative stress-response in yeast, replacement of the native promoter by this variant is sufficient alone to confer strongly enhanced resistance to H₂O₂ and cumene hydroperoxide. New replacement of this variant by the native promoter in the resistant strain suppresses the resistance. This work shows that increasing noise of target genes in a relevant strategy to engineer microbial strains toward better stress resistance. Multiple promoter replacement could synergize the effect observed here with the sole GRX1 promoter replacement. Finally this work suggests that combining several mutated copies of the target promoter could allow enhancing transcriptional-mediated noise at higher levels than mutating a single copy by providing constitutive bimodal and highly heterogeneous expression distribution.

The Culture Environment Influences Both Gene Regulation and Phenotypic Heterogeneity in Escherichia coli

Microorganisms shape the composition of the medium they are growing in, which in turn has profound consequences on the reprogramming of the population gene-expression profile. In this paper, we investigate the progressive changes in pH and sugar availability in the medium of a growing Escherichia coli (E. coli) culture. We show how these changes have an effect on both the cellular heterogeneity within the microbial community and the gene-expression profile of the microbial population. We measure the changes in gene-expression as E. coli moves from lag, to exponential, and finally into stationary phase. We found that pathways linked to the changes in the medium composition such as ribosomal, tricarboxylic acid cycle (TCA), transport, and metabolism pathways are strongly regulated during the different growth phases. In order to quantify the corresponding temporal changes in the population heterogeneity, we measure the fraction of E. coli persisters surviving different antibiotic treatments during the various phases of growth. We show that the composition of the medium in which β-lactams or quinolones, but not aminoglycosides, are dissolved strongly affects the measured phenotypic heterogeneity within the culture. Our findings contribute to a better understanding on how the composition of the culture medium influences both the reprogramming in the population gene-expression and the emergence of phenotypic variants.

Assigning function to natural allelic variation via dynamic modeling of gene network induction

More and more natural DNA variants are being linked to physiological traits. Yet, understanding what differences they make on molecular regulations remains challenging. Important properties of gene regulatory networks can be captured by computational models. If model parameters can be “personalized” according to the genotype, their variation may then reveal how DNA variants operate in the network. Here, we combined experiments and computations to visualize natural alleles of the yeast GAL3 gene in a space of model parameters describing the galactose response network. Alleles altering the activation of Gal3p by galactose were discriminated from those affecting its activity (production/degradation or efficiency of the activated protein). The approach allowed us to correctly predict that a non‐synonymous SNP would change the binding affinity of Gal3p with the Gal80p transcriptional repressor. Our results illustrate how personalizing gene regulatory models can be used for the mechanistic interpretation of genetic variants.

The emergence of metabolic heterogeneity and diverse growth responses in isogenic bacterial cells

Microorganisms adapt to frequent environmental changes through population diversification. Previous studies demonstrated phenotypic diversity in a clonal population and its important effects on microbial ecology. However, the dynamic changes of phenotypic composition have rarely been characterized. Also, cellular variations and environmental factors responsible for phenotypic diversity remain poorly understood. Here, we studied phenotypic diversity driven by metabolic heterogeneity. We characterized metabolic activities and growth kinetics of starved Escherichia coli cells subject to nutrient upshift at single-cell resolution. We observed three subpopulations with distinct metabolic activities and growth phenotypes. One subpopulation was metabolically active and immediately grew upon nutrient upshift. One subpopulation was metabolically inactive and non-viable. The other subpopulation was metabolically partially active, and did not grow upon nutrient upshift. The ratio of these subpopulations changed dynamically during starvation. A long-term observation of cells with partial metabolic activities indicated that their metabolism was later spontaneously restored, leading to growth recovery. Further investigations showed that oxidative stress can induce the emergence of a subpopulation with partial metabolic activities. Our findings reveal the emergence of metabolic heterogeneity and associated dynamic changes in phenotypic composition. In addition, the results shed new light on microbial dormancy, which has important implications in microbial ecology and biomedicine.

Molecular and cellular bases of adaptation to a changing environment in microorganisms

Environmental heterogeneity constitutes an evolutionary challenge for organisms. While evolutionary dynamics under variable conditions has been explored for decades, we still know relatively little about the cellular and molecular mechanisms involved. It is of paramount importance to examine these molecular bases because they may play an important role in shaping the course of evolution. In this review, we examine the diversity of adaptive mechanisms in the face of environmental changes. We exploit the recent literature on microbial systems because those have benefited the most from the recent emergence of genetic engineering and experimental evolution followed by genome sequencing. We identify four emerging trends: (i) an adaptive molecular change in a pathway often results in fitness trade-off in alternative environments but the effects are dependent on a mutation's genetic background; (ii) adaptive changes often modify transcriptional and signalling pathways; (iii) several adaptive changes may occur within the same molecular pathway but be associated with pleiotropy of different signs across environments; (iv) because of their large associated costs, macromolecular changes such as gene amplification and aneuploidy may be a rapid mechanism of adaptation in the short-term only. The course of adaptation in a variable environment, therefore, depends on the complexity of the environment but also on the molecular relationships among the genes involved and between the genes and the phenotypes under selection.

Genetic and Epigenetic Strategies Potentiate Gal4 Activation to Enhance Fitness in Recently Diverged Yeast Species

Certain genes show more rapid reactivation for several generations following repression, a conserved phenomenon called epigenetic transcriptional memory. Following previous growth in galactose, GAL gene transcriptional memory confers a strong fitness benefit in Saccharomyces cerevisiae adapting to growth in galactose for up to 8 generations. A genetic screen for mutants defective for GAL gene memory revealed new insights into the molecular mechanism, adaptive consequences, and evolutionary history of memory. A point mutation in the Gal1 co-activator that disrupts the interaction with the Gal80 inhibitor specifically and completely disrupted memory. This mutation confirms that cytoplasmically inherited Gal1 produced during previous growth in galactose directly interferes with Gal80 repression to promote faster induction of GAL genes. This mitotically heritable mode of regulation is recently evolved; in a diverged Saccharomyces species, GAL genes show constitutively faster activation due to genetically encoded basal expression of Gal1. Thus, recently diverged species utilize either epigenetic or genetic strategies to regulate the same molecular mechanism. The screen also revealed that the central domain of the Gal4 transcription factor both regulates the stochasticity of GAL gene expression and potentiates stronger GAL gene activation in the presence of Gal1. The central domain is critical for GAL gene transcriptional memory; Gal4 lacking the central domain fails to potentiate GAL gene expression and is unresponsive to previous Gal1 expression.