Phenotypic heterogeneity promotes adaptive evolution

Evolutionary regain of lost gene circuit function

Natural or synthetic genetic modules can lose their function over long-term evolution if the function is costly. How populations can evolve to restore such broken function is poorly understood. To test the reversibility of evolutionary breakdown, we use yeast cell populations with a chromosomally integrated synthetic gene circuit. In previous evolution experiments the gene circuit lost its costly function through various mutations. By exposing such mutant populations to conditions where regaining gene circuit function would be beneficial, we find adaptation scenarios with or without repairing lost gene circuit function. These results are important for drug resistance or future synthetic biology applications where evolutionary loss and regain of function play a significant role.

Evolutionary reversibility—the ability to regain a lost function—is an important problem both in evolutionary and synthetic biology, where repairing natural or synthetic systems broken by evolutionary processes may be valuable. Here, we use a synthetic positive-feedback (PF) gene circuit integrated into haploid Saccharomyces cerevisiae cells to test if the population can restore lost PF function. In previous evolution experiments, mutations in a gene eliminated the fitness costs of PF activation. Since PF activation also provides drug resistance, exposing such compromised or broken mutants to both drug and inducer should create selection pressure to regain drug resistance and possibly PF function. Indeed, evolving 7 PF mutant strains in the presence of drug revealed 3 adaptation scenarios through genomic, PF-external mutations that elevate PF basal expression, possibly by affecting transcription, translation, degradation, and other fundamental cellular processes. Nonfunctional mutants gained drug resistance without ever developing high expression, while quasifunctional and dysfunctional PF mutants developed high expression nongenetically, which then diminished, although more slowly for dysfunctional mutants where revertant clones arose. These results highlight how intracellular context, such as the growth rate, can affect regulatory network dynamics and evolutionary dynamics, which has important consequences for understanding the evolution of drug resistance and developing future synthetic biology applications.

Integrated evolutionary analysis reveals antimicrobial peptides with limited resistance

Antimicrobial peptides (AMPs) are promising antimicrobials, however, the potential of bacterial resistance is a major concern. Here we systematically study the evolution of resistance to 14 chemically diverse AMPs and 12 antibiotics in Escherichia coli. Our work indicates that evolution of resistance against certain AMPs, such as tachyplesin II and cecropin P1, is limited. Resistance level provided by point mutations and gene amplification is very low and antibiotic-resistant bacteria display no cross-resistance to these AMPs. Moreover, genomic fragments derived from a wide range of soil bacteria confer no detectable resistance against these AMPs when introduced into native host bacteria on plasmids. We have found that simple physicochemical features dictate bacterial propensity to evolve resistance against AMPs. Our work could serve as a promising source for the development of new AMP-based therapeutics less prone to resistance, a feature necessary to avoid any possible interference with our innate immune system.

Optogenetic Repressors of Gene Expression in Yeasts Using Light-Controlled Nuclear Localization

INTRODUCTION

Controlling gene expression is a fundamental goal of basic and synthetic biology because it allows insight into cellular function and control of cellular activity. We explored the possibility of generating an optogenetic repressor of gene expression in the model organism Saccharomyces cerevisiae by using light to control the nuclear localization of nuclease-dead Cas9, dCas9.

METHODS

The dCas9 protein acts as a repressor for a gene of interest when localized to the nucleus in the presence of an appropriate guide RNA (sgRNA). We engineered dCas9, the mammalian transcriptional repressor Mxi1, and an optogenetic tool to control nuclear localization (LINuS) as parts in an existing yeast optogenetic toolkit. This allowed expression cassettes containing novel dCas9 repressor configurations and guide RNAs to be rapidly constructed and integrated into yeast.

RESULTS

Our library of repressors displays a range of basal repression without the need for inducers or promoter modification. Populations of cells containing these repressors can be combined to generate a heterogeneous population of yeast with a 100-fold expression range. We find that repression can be dialed modestly in a light dose- and intensity-dependent manner. We used this library to repress expression of the lanosterol 14-alpha-demethylase Erg11, generating yeast with a range of sensitivity to the important antifungal drug fluconazole.

CONCLUSIONS

This toolkit will be useful for spatiotemporal perturbation of gene expression in Saccharomyces cerevisiae. Additionally, we believe that the simplicity of our scheme will allow these repressors to be easily modified to control gene expression in medically relevant fungi, such as pathogenic yeasts.

Intraspecies polymorphisms of Kluyveromyces marxianus strains from Yaghnob valley

In this study, 29 strains of Kluyveromyces marxianus with peculiar genetic and phenotypic traits previously isolated from a fermented goat milk of Yaghnob valley were investigated for chromosome length polymorphism (CLP) by PFGE, adhesion properties and carbon usage by Biolog analysis. Obtained data showed that strains differed in terms of number and size of chromosome bands. The number of bands ranged from 5 to 7, suggesting a probable genome size from 1.4 to 2.6 Mb. Strains showed a certain level of cell surface hydrophobicity ranging from 32% to 77.7%. Strains were also tested for their ability to form a biofilm on polystyrene plates: planktonic cells ranged from 6.3 cfu/mL to 7.95 cfu/mL, while sessile from 7.11 cfu/mL to 8.6 cfu/mL. The strains able to adhere to polystyrene plates were also able to form a mature MAT. Biolog analysis revealed that almost all strains were able to use putrescine, malic acid, α-D lactose, phenylethylamine, β-methyl D-gucoside and xylose; 5 strains were able to grow on cellobiose and 3 were able to catabolise α-ketobutyric. The obtained data highlighted a number of interesting features underlying the peculiar capacities of these strains for industrial applications.

Rapid Evolution of Reduced Susceptibility against a Balanced Dual-Targeting Antibiotic through Stepping-Stone Mutations

Multitargeting antibiotics, i.e., single compounds capable of inhibiting two or more bacterial targets, are generally considered to be a promising therapeutic strategy against resistance evolution. The rationale for this theory is that multitargeting antibiotics demand the simultaneous acquisition of multiple mutations at their respective target genes to achieve significant resistance. The theory presumes that individual mutations provide little or no benefit to the bacterial host. Here, we propose that such individual stepping-stone mutations can be prevalent in clinical bacterial isolates, as they provide significant resistance to other antimicrobial agents. To test this possibility, we focused on gepotidacin, an antibiotic candidate that selectively inhibits both bacterial DNA gyrase and topoisomerase IV. In a susceptible organism, Klebsiella pneumoniae, a combination of two specific mutations in these target proteins provide an >2,000-fold reduction in susceptibility, while individually, none of these mutations affect resistance significantly. Alarmingly, strains with decreased susceptibility against gepotidacin are found to be as virulent as the wild-type Klebsiella pneumoniae strain in a murine model. Moreover, numerous pathogenic isolates carry mutations which could promote the evolution of clinically significant reduction of susceptibility against gepotidacin in the future. As might be expected, prolonged exposure to ciprofloxacin, a clinically widely employed gyrase inhibitor, coselected for reduced susceptibility against gepotidacin. We conclude that extensive antibiotic usage could select for mutations that serve as stepping-stones toward resistance against antimicrobial compounds still under development. Our research indicates that even balanced multitargeting antibiotics are prone to resistance evolution.

A Quantitative Model Explains Single-Cell Dynamics of the Adaptive Response in Escherichia coli

DNA damage caused by alkylating chemicals induces an adaptive response in Escherichia coli that increases the tolerance of cells to further damage. Signaling of the response occurs through irreversible methylation of the Ada protein, which acts as a DNA repair protein and damage sensor. Methylated Ada induces its own gene expression through a positive feedback loop. However, random fluctuations in the abundance of Ada jeopardize the reliability of the induction signal. I developed a quantitative model to test how gene expression noise and feedback amplification affect the fidelity of the adaptive response. A remarkably simple model accurately reproduced experimental observations from single-cell measurements of gene expression dynamics in a microfluidic device. Stochastic simulations showed that delays in the adaptive response are a direct consequence of the very low number of Ada molecules present to signal DNA damage. For cells that have zero copies of Ada, response activation becomes a memoryless process that is dictated by an exponential waiting time distribution between basal Ada expression events. Experiments also confirmed the model prediction that the strength of the adaptive response drops with an increasing growth rate of cells.

Every which way? On predicting tumor evolution using cancer progression models

Successful prediction of the likely paths of tumor progression is valuable for diagnostic, prognostic, and treatment purposes. Cancer progression models (CPMs) use cross-sectional samples to identify restrictions in the order of accumulation of driver mutations and thus CPMs encode the paths of tumor progression. Here we analyze the performance of four CPMs to examine whether they can be used to predict the true distribution of paths of tumor progression and to estimate evolutionary unpredictability. Employing simulations we show that if fitness landscapes are single peaked (have a single fitness maximum) there is good agreement between true and predicted distributions of paths of tumor progression when sample sizes are large, but performance is poor with the currently common much smaller sample sizes. Under multi-peaked fitness landscapes (i.e., those with multiple fitness maxima), performance is poor and improves only slightly with sample size. In all cases, detection regime (when tumors are sampled) is a key determinant of performance. Estimates of evolutionary unpredictability from the best performing CPM, among the four examined, tend to overestimate the true unpredictability and the bias is affected by detection regime; CPMs could be useful for estimating upper bounds to the true evolutionary unpredictability. Analysis of twenty-two cancer data sets shows low evolutionary unpredictability for several of the data sets. But most of the predictions of paths of tumor progression are very unreliable, and unreliability increases with the number of features analyzed. Our results indicate that CPMs could be valuable tools for predicting cancer progression but that, currently, obtaining useful predictions of paths of tumor progression from CPMs is dubious, and emphasize the need for methodological work that can account for the probably multi-peaked fitness landscapes in cancer.

Developmental noise and ecological opportunity across space can release constraints on the evolution of plasticity

Phenotypic plasticity is a potentially definitive solution to environment heterogeneity, driving biologists to understand why it is not ubiquitous in nature. While costs and constraints may limit the success of plasticity, we are still far from a complete theory of when these limitations actually proscribe adaptive plasticity. Here I use a simple model of plasticity incorporating developmental noise to explore the competitive and evolutionary relationships of specialist and generalist genotypes spreading across a heterogeneous landscape. Results show that plasticity can arise in the context of specialism, preadapting genotypes to later evolve toward plastic generalism. Developmental noise helps a mutant with imperfect plasticity successfully compete against its ancestor, providing an evolutionary path by which subsequent mutations can refine plasticity toward its optimum. These results address how the complex selection pressures across a heterogeneous environment can help evolution find paths around constraints arising from developmental mechanisms.

Role of network-mediated stochasticity in mammalian drug resistance

A major challenge in biology is that genetically identical cells in the same environment can display gene expression stochasticity (noise), which contributes to bet-hedging, drug tolerance, and cell-fate switching. The magnitude and timescales of stochastic fluctuations can depend on the gene regulatory network. Currently, it is unclear how gene expression noise of specific networks impacts the evolution of drug resistance in mammalian cells. Answering this question requires adjusting network noise independently from mean expression. Here, we develop positive and negative feedback-based synthetic gene circuits to decouple noise from the mean for Puromycin resistance gene expression in Chinese Hamster Ovary cells. In low Puromycin concentrations, the high-noise, positive-feedback network delays long-term adaptation, whereas it facilitates adaptation under high Puromycin concentration. Accordingly, the low-noise, negative-feedback circuit can maintain resistance by acquiring mutations while the positive-feedback circuit remains mutation-free and regains drug sensitivity. These findings may have profound implications for chemotherapeutic inefficiency and cancer relapse.

The role of gene expression noise in the evolution of drug resistance in mammalian cells is unclear. Here, by uncoupling noise from mean expression of a drug resistance gene in CHO cells the authors show that noisy expression aids adaptation to high drug levels, but delays it at low drug levels.

Gene Expression Noise Produces Cell-to-Cell Heterogeneity in Eukaryotic Homologous Recombination Rate

Variation in gene expression among genetically identical individual cells (called gene expression noise) directly contributes to phenotypic diversity. Whether such variation can impact genome stability and lead to variation in genotype remains poorly explored. We addressed this question by investigating whether noise in the expression of genes affecting homologous recombination (HR) activity either directly (RAD52) or indirectly (RAD27) confers cell-to-cell heterogeneity in HR rate in Saccharomyces cerevisiae. Using cell sorting to isolate subpopulations with various expression levels, we show that spontaneous HR rate is highly heterogeneous from cell-to-cell in clonal populations depending on the cellular amount of proteins affecting HR activity. Phleomycin-induced HR is even more heterogeneous, showing that RAD27 expression variation strongly affects the rate of recombination from cell-to-cell. Strong variations in HR rate between subpopulations are not correlated to strong changes in cell cycle stage. Moreover, this heterogeneity occurs even when simultaneously sorting cells at equal expression level of another gene involved in DNA damage response (BMH1) that is upregulated by DNA damage, showing that the initiating DNA damage is not responsible for the observed heterogeneity in HR rate. Thus gene expression noise seems mainly responsible for this phenomenon. Finally, HR rate non-linearly scales with Rad27 levels showing that total amount of HR cannot be explained solely by the time- or population-averaged Rad27 expression. Altogether, our data reveal interplay between heterogeneity at the gene expression and genetic levels in the production of phenotypic diversity with evolutionary consequences from microbial to cancer cell populations.

Gene expression variability: the other dimension in transcriptome analysis

Transcriptome sequencing is a powerful technique to study molecular changes that underlie the differences in physiological conditions and disease progression. A typical question that is posed in such studies is finding genes with significant changes between sample groups. In this respect expression variability is regarded as a nuisance factor that is primarily of technical origin and complicates the data analysis. However, it is becoming apparent that the biological variation in gene expression might be an important molecular phenotype that can affect physiological parameters. In this review we explore the recent literature on technical and biological variability in gene expression, sources of expression variability, (epi-)genetic hallmarks, and evolutionary constraints in genes with robust and variable gene expression. We provide an overview of recent findings on effects of external cues, such as diet and aging, on expression variability and on other biological phenomena that can be linked to it. We discuss metrics and tools that were developed for quantification of expression variability and highlight the importance of future studies in this direction. To assist the adoption of expression variability analysis, we also provide a detailed description and computer code, which can easily be utilized by other researchers. We also provide a reanalysis of recently published data to highlight the value of the analysis method.

Intratumoral Heterogeneity: More Than Just Mutations

Most human tumors are composed of genetically and phenotypically heterogeneous cancer cell populations, which poses a major challenge for the clinical management of cancer patients. Advances of single-cell technologies have allowed the profiling of tumors at unprecedented depth, which, in combination with newly developed computational tools, enable the dissection of tumor evolution with increasing precision. However, our understanding of mechanisms that regulate intratumoral heterogeneity and our ability to modulate it has been lagging behind. Recent data demonstrate that epigenetic regulators, including histone demethylases, may control the cell-to-cell variability of transcriptomes and chromatin profiles and they may modulate therapeutic responses via this function. Thus, the therapeutic targeting of epigenetic enzymes may be used to decrease intratumoral cellular heterogeneity and treatment resistance, when used in combination with other types of agents.

Modeling cell population dynamics

 Quantitative modeling is quickly becoming an integral part of biology, due to the ability of mathematical models and computer simulations to generate insights and predict the behavior of living systems. Single-cell models can be incapable or misleading for inferring population dynamics, as they do not consider the interactions between cells via metabolites or physical contact, nor do they consider competition for limited resources such as nutrients or space. Here we examine methods that are commonly used to model and simulate cell populations. First, we cover simple models where analytic solutions are available, and then move on to more complex scenarios where computational methods are required. Overall, we present a summary of mathematical models used to describe cell population dynamics, which may aid future model development and highlights the importance of population modeling in biology.

Modeling cell population dynamics

 Quantitative modeling is quickly becoming an integral part of biology, due to the ability of mathematical models and computer simulations to generate insights and predict the behavior of living systems. Single-cell models can be incapable or misleading for inferring population dynamics, as they do not consider the interactions between cells via metabolites or physical contact, nor do they consider competition for limited resources such as nutrients or space. Here we examine methods that are commonly used to model and simulate cell populations. First, we cover simple models where analytic solutions are available, and then move on to more complex scenarios where computational methods are required. Overall, we present a summary of mathematical models used to describe cell population dynamics, which may aid future model development and highlights the importance of population modeling in biology.

Control of nongenetic heterogeneity in growth rate and stress tolerance of Saccharomyces cerevisiae by cyclic AMP-regulated transcription factors

Genetically identical cells exhibit extensive phenotypic variation even under constant and benign conditions. This so-called nongenetic heterogeneity has important clinical implications: within tumors and microbial infections, cells show nongenetic heterogeneity in growth rate and in susceptibility to drugs or stress. The budding yeast, Saccharomyces cerevisiae, shows a similar form of nongenetic heterogeneity in which growth rate correlates positively with susceptibility to acute heat stress at the single-cell level. Using genetic and chemical perturbations, combined with high-throughput single-cell assays of yeast growth and gene expression, we show here that heterogeneity in intracellular cyclic AMP (cAMP) levels acting through the conserved Ras/cAMP/protein kinase A (PKA) pathway and its target transcription factors, Msn2 and Msn4, underlies this nongenetic heterogeneity. Lower levels of cAMP correspond to slower growth, as shown by direct comparison of cAMP concentration in subpopulations enriched for slower vs. faster growing cells. Concordantly, an endogenous reporter of this pathway’s activity correlates with growth in individual cells. The paralogs Msn2 and Msn4 differ in their roles in nongenetic heterogeneity in a way that demonstrates slow growth and stress tolerance are not inevitably linked. Heterogeneity in growth rate requires each, whereas only Msn2 is required for heterogeneity in expression of Tsl1, a subunit of trehalose synthase that contributes to acute-stress tolerance. Perturbing nongenetic heterogeneity by mutating genes in this pathway, or by culturing wild-type cells with the cell-permeable cAMP analog 8-bromo-cAMP or the PKA inhibitor H89, significantly impacts survival of acute heat stress. Perturbations that increase intracellular cAMP levels reduce the slower-growing subpopulation and increase susceptibility to acute heat stress, whereas PKA inhibition slows growth and decreases susceptibility to acute heat stress. Loss of Msn2 reduces, but does not completely eliminate, the correlation in individual cells between growth rate and acute-stress survival, suggesting a major role for the Msn2 pathway in nongenetic heterogeneity but also a residual benefit of slow growth. Our results shed light on the genetic control of nongenetic heterogeneity and suggest a possible means of defeating bet-hedging pathogens or tumor cells by making them more uniformly susceptible to treatment.

Nongenetic heterogeneity exists when a trait differs among individuals that have identical genotypes and environments. A clonal population can maximize its long-term success in an uncertain environment by diversifying its phenotypes via nongenetic heterogeneity: the currently unfavored ones may become the favored ones when conditions change. Nongenetic heterogeneity has clinical relevance. For example, populations of tumor cells or infectious microbes show cell-to-cell differences in growth and in drug or stress tolerance. This heterogeneity hampers efficient treatment and can potentiate harmful evolution of a tumor or pathogen. We show that in budding yeast, heterogeneity in intracellular cyclic AMP levels acting through the conserved Ras/cAMP/protein kinase A (PKA) pathway and its target transcription factors, Msn2 and Msn4, underlies the nongenetic heterogeneity of both single-cell growth rate and acute heat-stress tolerance. Perturbations of this pathway significantly affect population survival upon acute heat stress. These results illuminate a mechanism of nongenetic heterogeneity and suggest the potential value of antitumor or antifungal treatment strategies that target nongenetic heterogeneity to render the tumor or pathogen population more uniformly susceptible to a second drug that aims to kill.

NanoSIMS single cell analyses reveal the contrasting nitrogen sources for small phytoplankton

Nitrogen (N) is a limiting nutrient in vast regions of the world's oceans, yet the sources of N available to various phytoplankton groups remain poorly understood. In this study, we investigated inorganic carbon (C) fixation rates and nitrate (NO₃^-), ammonium (NH₄⁺) and urea uptake rates at the single cell level in photosynthetic pico-eukaryotes (PPE) and the cyanobacteria Prochlorococcus and Synechococcus. To that end, we used dual ¹⁵N and ¹³C-labeled incubation assays coupled to flow cytometry cell sorting and nanoSIMS analysis on samples collected in the North Pacific Subtropical Gyre (NPSG) and in the California Current System (CCS). Based on these analyses, we found that photosynthetic growth rates (based on C fixation) of PPE were higher in the CCS than in the NSPG, while the opposite was observed for Prochlorococcus. Reduced forms of N (NH₄⁺ and urea) accounted for the majority of N acquisition for all the groups studied. NO₃^- represented a reduced fraction of total N uptake in all groups but was higher in PPE (17.4 ± 11.2% on average) than in Prochlorococcus and Synechococcus (4.5 ± 6.5 and 2.9 ± 2.1% on average, respectively). This may in part explain the contrasting biogeography of these picoplankton groups. Moreover, single cell analyses reveal that cell-to-cell heterogeneity within picoplankton groups was significantly greater for NO₃^- uptake than for C fixation and NH₄⁺ uptake. We hypothesize that cellular heterogeneity in NO₃^- uptake within groups facilitates adaptation to the fluctuating availability of NO₃^- in the environment.

Synthetic circuits reveal how mechanisms of gene regulatory networks constrain evolution

Phenotypic variation is the raw material of adaptive Darwinian evolution. The phenotypic variation found in organismal development is biased towards certain phenotypes, but the molecular mechanisms behind such biases are still poorly understood. Gene regulatory networks have been proposed as one cause of constrained phenotypic variation. However, most pertinent evidence is theoretical rather than experimental. Here, we study evolutionary biases in two synthetic gene regulatory circuits expressed in Escherichia coli that produce a gene expression stripe-a pivotal pattern in embryonic development. The two parental circuits produce the same phenotype, but create it through different regulatory mechanisms. We show that mutations cause distinct novel phenotypes in the two networks and use a combination of experimental measurements, mathematical modelling and DNA sequencing to understand why mutations bring forth only some but not other novel gene expression phenotypes. Our results reveal that the regulatory mechanisms of networks restrict the possible phenotypic variation upon mutation. Consequently, seemingly equivalent networks can indeed be distinct in how they constrain the outcome of further evolution.

Duplications and losses of genes encoding known elements of the stress defence system of the Aspergilli contribute to the evolution of these filamentous fungi but do not directly influence their environmental stress tolerance

The contribution of stress protein duplication and deletion events to the evolution of the Aspergilli was studied. We performed a large-scale homology analysis of stress proteins and generated and analysed three stress defence system models based on Saccharomyces cerevisiae, Schizosaccharomyces pombe and Aspergillus nidulans. Although both yeast-based and A. nidulans-based models were suitable to trace evolutionary changes, the A. nidulans-based model performed better in mapping stress protein radiations. The strong Mantel correlation found between the positions of species in the phylogenetic tree on the one hand and either in the A. nidulans-based or S. cerevisiae-based models on the other hand demonstrated that stress protein expansions and reductions contributed significantly to the evolution of the Aspergilli. Interestingly, stress tolerance attributes correlated well with the number of orthologs only for a few stress proteins. Notable examples are Ftr1 iron permease and Fet3 ferro-O2-oxidoreductase, elements of the reductive iron assimilation pathway, in the S. cerevisiae-based model, as well as MpkC, a HogA-like mitogen activated protein kinase in the A. nidulans-based model. In the case of the iron assimilation proteins, the number of orthologs showed a positive correlation with H2O2-induced stress tolerance while the number of MpkC orthologs correlated positively with Congo Red induced cell wall stress, sorbitol induced osmotic stress and H2O2 induced oxidative stress tolerances. For most stress proteins, changes in the number of orthologs did not correlate well with any stress tolerance attributes. As a consequence, stress tolerance patterns of the studied Aspergilli did not correlate with either the sets of stress response proteins in general or with the phylogeny of the species studied. These observations suggest that stress protein duplication and deletion events significantly contributed to the evolution of stress tolerance attributes of Aspergilli. In contrast, there are other processes, which may counterbalance the effects of stress gene duplications or deletions including (i) alterations in the structures of stress proteins leading to changes in their biological activities, (ii) varying biosynthesis of stress proteins, (iii) rewiring stress response regulatory networks or even (iv) acquiring new stress response genes by horizontal gene transfer. All these multilevel changes are indispensable for the successful adaptation of filamentous fungi to altering environmental conditions, especially when these organisms are entering new ecological niches.

Experimental Design, Population Dynamics, and Diversity in Microbial Experimental Evolution

In experimental evolution, laboratory-controlled conditions select for the adaptation of species, which can be monitored in real time. Despite the current popularity of such experiments, nature's most pervasive biological force was long believed to be observable only on time scales that transcend a researcher's life-span, and studying evolution by natural selection was therefore carried out solely by comparative means. Eventually, microorganisms' propensity for fast evolutionary changes proved us wrong, displaying strong evolutionary adaptations over a limited time, nowadays massively exploited in laboratory evolution experiments. Here, we formulate a guide to experimental evolution with microorganisms, explaining experimental design and discussing evolutionary dynamics and outcomes and how it is used to assess ecoevolutionary theories, improve industrially important traits, and untangle complex phenotypes. Specifically, we give a comprehensive overview of the setups used in experimental evolution. Additionally, we address population dynamics and genetic or phenotypic diversity during evolution experiments and expand upon contributing factors, such as epistasis and the consequences of (a)sexual reproduction. Dynamics and outcomes of evolution are most profoundly affected by the spatiotemporal nature of the selective environment, where changing environments might lead to generalists and structured environments could foster diversity, aided by, for example, clonal interference and negative frequency-dependent selection. We conclude with future perspectives, with an emphasis on possibilities offered by fast-paced technological progress. This work is meant to serve as an introduction to those new to the field of experimental evolution, as a guide to the budding experimentalist, and as a reference work to the seasoned expert.

An Artificial Yeast Genetic Circuit Enables Deep Mutational Scanning of an Antimicrobial Resistance Protein

Understanding the molecular mechanisms underlying antibiotic resistance requires concerted efforts in enzymology and medicinal chemistry. Here we describe a new synthetic biology approach to antibiotic development, where the presence of tetracycline antibiotics is linked to a life-death selection in Saccharomyces cerevisiae. This artificial genetic circuit allowed the deep mutational scanning of the tetracycline inactivating enzyme TetX, revealing key functional residues. We used both positive and negative selections to confirm the importance of different residues for TetX activity, and profiled activity hotspots for different tetracyclines to reveal substrate-specific activity determinants. We found that precise positioning of FAD and hydrophobic shielding of the tetracycline are critical for enzymatic inactivation of doxycycline. However, positioning of FAD is suboptimal in the case of anhydrotetracycline, potentially explaining its comparatively poor degradation and potential as an inhibitor for this family of enzymes. By combining artificial genetic circuits whose function can be modulated by antimicrobial resistance determinants, we establish a framework to select for the next generation of antibiotics.

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.

Adaptive response of yeast cells to triggered toxicity of phosphoribulokinase

Adjustment of plasmid copy number resulting from the balance between positive and negative impacts of borne synthetic genes, plays a critical role in the global efficiency of multistep metabolic engineering. Differential expression of co-expressed engineered genes is frequently observed depending on growth phases, metabolic status and triggered adjustments of plasmid copy numbers, constituting a dynamic process contributing to minimize global engineering burden. A yeast model involving plasmid based expression of phosphoribulokinase (PRKp), a key enzyme for the reconstruction of synthetic Calvin cycle, was designed to gain further insights into such a mechanism. A conditional PRK expression cassette was cloned either onto a low (ARS-CEN based) or a high (2-micron origin based) copy number plasmid using complementation of a trp1 genomic mutation as constant positive selection. Evolution of plasmid copy numbers, PRKp expressions, and cell growth rates were dynamically monitored following gene de-repression through external doxycycline concentration shifts. In the absence of RubisCO encoding gene permitting metabolic recycling, PRKp expression that led to depletion of ribulose phosphate, a critical metabolite for aromatic amino-acids biosynthesis, and accumulation of the dead-end diphosphate product contribute to toxicity. Triggered copy number adjustment was found to be a dynamic process depending both on plasmid types and levels of PRK induction. With the ARS-CEN plasmid, cell growth was abruptly affected only when level PRKp expression exceeded a threshold value. In contrast, a proportional relationship was observed with the 2-micron plasmid consistent with large copy number adjustments. Micro-compartment partitioning of bulk cultures by embedding individual cells into inverse culture medium/oil droplets, revealed the presence of slow and fast growing subpopulations that differ in relative proportions for low and high copy number plasmids.

CpG island composition differences are a source of gene expression noise indicative of promoter responsiveness

BACKGROUND

Population phenotypic variation can arise from genetic differences between individuals, or from cellular heterogeneity in an isogenic group of cells or organisms. The emergence of gene expression differences between genetically identical cells is referred to as gene expression noise, the sources of which are not well understood.

RESULTS

In this work, by studying gene expression noise between multiple cell lineages and mammalian species, we find consistent evidence of a role for CpG islands as sources of gene expression noise. Variation in noise among CpG island promoters can be partially attributed to differences in island size, in which short islands have noisier gene expression. Building on these findings, we investigate the potential for short CpG islands to act as fast response elements to environmental stimuli. Specifically, we find that these islands are enriched amongst primary response genes in SWI/SNF-independent stimuli, suggesting that expression noise is an indicator of promoter responsiveness.

CONCLUSIONS

Thus, through the integration of single-cell RNA expression profiling, chromatin landscape and temporal gene expression dynamics, we have uncovered a role for short CpG island promoters as fast response elements.

Shift from stochastic to spatially-ordered expression of serine-glycine synthesis enzymes in 3D microtumors

Cell-to-cell differences in protein expression in normal tissues and tumors are a common phenomenon, but the underlying principles that govern this heterogeneity are largely unknown. Here, we show that in monolayer cancer cell-line cultures, the expression of the five metabolic enzymes of serine-glycine synthesis (SGS), including its rate-limiting enzyme, phosphoglycerate dehydrogenase (PHGDH), displays stochastic cell-to-cell variation. By contrast, in cancer cell line-derived three-dimensional (3D) microtumors PHGDH expression is restricted to the outermost part of the microtumors' outer proliferative cell layer, while the four other SGS enzymes display near uniform expression throughout the microtumor. A mathematical model suggests that metabolic stress in the microtumor core activates factors that restrict PHGDH expression. Thus, intracellular enzyme expression in growing cell ecosystems can shift to spatially ordered patterns in 3D structured environments due to emergent cell-cell communication, with potential implications for the design of effective anti-metabolic cancer therapies.

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.

Antibiotic-resistant bacteria show widespread collateral sensitivity to antimicrobial peptides

Antimicrobial peptides are promising alternative antimicrobial agents. However, little is known about whether resistance to small-molecule antibiotics leads to cross-resistance (decreased sensitivity) or collateral sensitivity (increased sensitivity) to antimicrobial peptides. We systematically addressed this question by studying the susceptibilities of a comprehensive set of 60 antibiotic-resistant Escherichia coli strains towards 24 antimicrobial peptides. Strikingly, antibiotic-resistant bacteria show a high frequency of collateral sensitivity to antimicrobial peptides, whereas cross-resistance is relatively rare. We identify clinically relevant multidrug-resistance mutations that increase bacterial sensitivity to antimicrobial peptides. Collateral sensitivity in multidrug-resistant bacteria arises partly through regulatory changes shaping the lipopolysaccharide composition of the bacterial outer membrane. These advances allow the identification of antimicrobial peptide-antibiotic combinations that enhance antibiotic activity against multidrug-resistant bacteria and slow down de novo evolution of resistance. In particular, when co-administered as an adjuvant, the antimicrobial peptide glycine-leucine-amide caused up to 30-fold decrease in the antibiotic resistance level of resistant bacteria. Our work provides guidelines for the development of efficient peptide-based therapies of antibiotic-resistant infections.

Behavioral Variability and Phenotypic Diversity in Bacterial Chemotaxis

Living cells detect and process external signals using signaling pathways that are affected by random fluctuations. These variations cause the behavior of individual cells to fluctuate over time (behavioral variability) and generate phenotypic differences between genetically identical individuals (phenotypic diversity). These two noise sources reduce our ability to predict biological behavior because they diversify cellular responses to identical signals. Here, we review recent experimental and theoretical advances in understanding the mechanistic origin and functional consequences of such variation in Escherichia coli chemotaxis-a well-understood model of signal transduction and behavior. After briefly summarizing the architecture and logic of the chemotaxis system, we discuss determinants of behavior and chemotactic performance of individual cells. Then, we review how cell-to-cell differences in protein abundance map onto differences in individual chemotactic abilities and how phenotypic variability affects the performance of the population. We conclude with open questions to be addressed by future research.

Ergodicity, hidden bias and the growth rate gain

Many single-cell observables are highly heterogeneous. A part of this heterogeneity stems from age-related phenomena: the fact that there is a nonuniform distribution of cells with different ages. This has led to a renewed interest in analytic methodologies including use of the 'von Foerster equation' for predicting population growth and cell age distributions. Here we discuss how some of the most popular implementations of this machinery assume a strong condition on the ergodicity of the cell cycle duration ensemble. We show that one common definition for the term ergodicity, 'a single individual observed over many generations recapitulates the behavior of the entire ensemble' is implied by the other, 'the probability of observing any state is conserved across time and over all individuals' in an ensemble with a fixed number of individuals but that this is not true when the ensemble is growing. We further explore the impact of generational correlations between cell cycle durations on the population growth rate. Finally, we explore the 'growth rate gain'-the phenomenon that variations in the cell cycle duration leads to an improved population-level growth rate-in this context. We highlight that, fundamentally, this effect is due to asymmetric division.

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.

Epigenetic and Transcriptional Variability Shape Phenotypic Plasticity

Epigenetic and transcriptional variability contribute to the vast diversity of cellular and organismal phenotypes and are key in human health and disease. In this review, we describe different types, sources, and determinants of epigenetic and transcriptional variability, enabling cells and organisms to adapt and evolve to a changing environment. We highlight the latest research and hypotheses on how chromatin structure and the epigenome influence gene expression variability. Further, we provide an overview of challenges in the analysis of biological variability. An improved understanding of the molecular mechanisms underlying epigenetic and transcriptional variability, at both the intra- and inter-individual level, provides great opportunity for disease prevention, better therapeutic approaches, and personalized medicine.

Impact of Morphological Sectors on Antifungal Susceptibility Testing and Virulence Studies

Morphological heterogeneity of Aspergillus terreus cultures was observed during continued cultivation of amphotericin B (AMB)-resistant isolates on drug-free medium. Outgrowth leads to the emergence of multiple sectors that might result from increased growth rates at drug-free conditions. We evaluated the differences in AMB susceptibility and virulence between sector subcultures (ATSec), AMB-resistant (ATR) strains, and AMB-susceptible (ATS) strains. By comparing A. terreus AMB-resistant (ATR) strains and A. terreus sector (ATSec) cultures we observed a highly significant reduction of AMB MICs in ATSec (ATR MIC, 2 to 32 μg/ml; ATSec MIC, 0.12 to 5 μg/ml). Furthermore, Galleria mellonella survival studies revealed an enhanced virulence of ATSec, which was comparable with that of AMB-sensitive Aspergillus terreus strains (median survival rates for ATS isolates, 72 h; for ATSec isolate ATSec^G1, 84 h; for ATR isolates, 144 h). Our findings clearly demonstrate that spontaneous culture degeneration occurs in A. terreus and, most importantly, crucially impacts drug efficacy and virulence.

Transcription termination factor Rho and microbial phenotypic heterogeneity

Populations of genetically identical microorganisms exhibit high degree of cell-to-cell phenotypic diversity even when grown in uniform environmental conditions. Heterogeneity is a genetically determined trait, which ensures bacterial adaptation and survival in the ever changing environmental conditions. Fluctuations in gene expression (noise) at the level of transcription initiation largely contribute to cell-to-cell variability within population. Not surprisingly, the analyses of the mechanisms driving phenotypic heterogeneity are mainly focused on the activity of promoters and transcriptional factors. Less attention is currently given to a role of intrinsic and factor-dependent transcription terminators. Here, we discuss recent advances in understanding the regulatory role of the multi-functional transcription termination factor Rho, the major inhibitor of pervasive transcription in bacteria and the emerging global regulator of gene expression. We propose that termination activity of Rho might be among the mechanisms by which cells manage the intensity of transcriptional noise, thus affecting population heterogeneity.

From transcriptional complexity to cellular phenotypes: Lessons from yeast

Pervasive transcription has been reported in many eukaryotic organisms, revealing a highly interleaved transcriptome organization that involves thousands of coding and non-coding RNAs. However, to date, the biological impact of transcriptome complexity is still poorly understood. Here I will review how subtle variations of the transcriptome can lead to divergent cellular phenotypes by fine-tuning both its coding potential and regulation. I will discuss strategies that can be used to link molecular variations with divergent biological outcomes. Finally, I will explore the implication of transcriptional complexity for our understanding of gene expression in the context of cell-to-cell phenotypic variability. Copyright © 2017 John Wiley & Sons, Ltd.