Dynamics of Adaptation During Three Years of Evolution Under Long-Term Stationary Phase

Different RNA recognition by ProQ and FinO depends on the sequence surrounding intrinsic terminator hairpins

Escherichia coli ProQ and FinO proteins both have RNA-binding FinO domains, which bind to intrinsic transcription terminators, but each protein recognizes distinct RNAs. To explore how ProQ and FinO discriminate between RNAs, we transplanted sequences surrounding terminator hairpins between RNAs specific for each protein, and compared their binding to ProQ, the isolated FinO domain of ProQ (ProQNTD), and FinO. The results showed that the binding specificity of chimeric RNAs toward ProQ, ProQNTD, or FinO was determined by the origin of the transplanted sequence. Further analysis showed that the sequence surrounding the terminator hairpin, including a purine-purine mismatch, in natural RNA ligands of FinO and in chimeric RNAs, weakened their binding by ProQNTD Overall, our studies suggest that RNA sequence elements surrounding the intrinsic terminator hairpin contribute to the discrimination between RNAs by ProQ and FinO.

Does strain-level persistence of lactobacilli in long-term back-slopped sourdoughs inform on domestication of food-fermenting lactic acid bacteria?

Sourdoughs are maintained by back-slopping over long time periods. To determine strain-level persistence of bacteria, we characterized four sourdoughs from bakeries over a period of 3.3, 11.0, 18.0, and 19.0 years. One sourdough included isolates of Levilactobacillus spp. and Fructilactobacillus spp. that differed by fewer than 10 single-nucleotide polymorphisms (SNPs) from the isolates obtained 3.3 years earlier and thus likely represent the same strain. Isolates of Levilactobacillus parabrevis differed by 200-300 SNPs; their genomes were under positive selection, indicating transmission from an external source. In two other sourdoughs, isolates of Fructilactobacillus sanfranciscensis that were obtained 11 and 18 years apart differed by 19 and 29 SNPs, respectively, again indicating repeated isolation of the same strain. The isolate of Fl. sanfranciscensis from the fourth sourdough differed by 45 SNPs from the isolate obtained 19 years previously. We thus identified strain-level persistence in three out of four long-term back-slopped sourdoughs, making it possible that strains persisted over periods that are long enough to allow bacterial speciation and domestication.IMPORTANCEThe assembly of microbial communities in sourdough is shaped by dispersal and selection. Speciation and domestication of fermentation microbes in back-slopped food fermentations have been documented for food-fermenting fungi including sourdough yeasts but not for bacteria, which evolve at a slower rate. Bacterial speciation in food fermentations requires strain-level persistence of fermentation microbes over hundreds or thousands of years. By documenting strain-level persistence in three out of four sourdoughs over a period of up to 18 years, we demonstrate that persistence over hundreds or thousands of years is possible, if not likely. We thus not only open a new perspective on fermentation control in bakeries but also support the possibility that all humans, despite their cultural diversity, share the same fermentation microbes.

Evolution of pH-sensitive transcription termination in Escherichia coli during adaptation to repeated long-term starvation

Fluctuating environments that consist of regular cycles of co-occurring stress are a common challenge faced by cellular populations. For a population to thrive in constantly changing conditions, an ability to coordinate a rapid cellular response is essential. Here, we identify a mutation conferring an arginine-to-histidine (Arg to His) substitution in the transcription terminator Rho. The rho R109H mutation frequently arose in Escherichia coli populations experimentally evolved under repeated long-term starvation conditions, during which the accumulation of metabolic waste followed by transfer into fresh media results in drastic environmental pH fluctuations associated with feast and famine. Metagenomic sequencing revealed that populations containing the rho mutation also possess putative loss-of-function mutations in ydcI, which encodes a recently characterized transcription factor associated with pH homeostasis. Genetic reconstructions of these mutations show that the rho allele confers plasticity via an alkaline-induced reduction of Rho function that, when found in tandem with a ΔydcI allele, leads to intracellular alkalization and genetic assimilation of Rho mutant function. We further identify Arg to His substitutions at analogous sites in rho alleles from species that regularly experience neutral to alkaline pH fluctuations in their environments. Our results suggest that Arg to His substitutions in Rho may serve to rapidly coordinate complex physiological responses through pH sensing and shed light on how cellular populations use environmental cues to coordinate rapid responses to complex, fluctuating environments.

Pseudomonas putida Dynamics of Adaptation under Prolonged Resource Exhaustion

Many nonsporulating bacterial species survive prolonged resource exhaustion, by entering a state termed long-term stationary phase. Here, we performed long-term stationary phase evolutionary experiments on the bacterium Pseudomonas putida, followed by whole-genome sequencing of evolved clones. We show that P. putida is able to persist and adapt genetically under long-term stationary phase. We observed an accumulation of mutations within the evolving P. putida populations. Within each population, independently evolving lineages are established early on and persist throughout the 4-month-long experiment. Mutations accumulate in a highly convergent manner, with similar loci being mutated across independently evolving populations. Across populations, mutators emerge, that due to mutations within mismatch repair genes developed a much higher rate of mutation than other clones with which they coexisted within their respective populations. While these general dynamics of the adaptive process are quite similar to those we previously observed in the model bacterium Escherichia coli, the specific loci that are involved in adaptation only partially overlap between P. putida and E. coli.

Escherichia coli adaptation under prolonged resource exhaustion is characterized by extreme parallelism and frequent historical contingency

Like many other non-sporulating bacterial species, Escherichia coli is able to survive prolonged periods of resource exhaustion, by entering a state of growth called long-term stationary phase (LTSP). In July 2015, we initiated a set of evolutionary experiments aimed at characterizing the dynamics of E. coli adaptation under LTSP. In these experiments populations of E. coli were allowed to initially grow on fresh rich media, but were not provided with any new external growth resources since their establishment. Utilizing whole genome sequencing data obtained for hundreds of clones sampled at 12 time points spanning the first six years of these experiments, we reveal several novel aspects of the dynamics of adaptation. First, we show that E. coli continuously adapts genetically, up to six years under resource exhaustion, through the highly convergent accumulation of mutations. We further show that upon entry into LTSP, long-lasting lineages are established. This lineage structure is in itself convergent, with similar lineages arising across independently evolving populations. The high parallelism with which adaptations occur under LTSP, combined with the LTSP populations' lineage structure, enable us to screen for pairs of loci displaying a significant association in the occurrence of mutations, suggestive of a historical contingency. We find that such associations are highly frequent and that a third of convergently mutated loci are involved in at least one such association. Combined our results demonstrate that LTSP adaptation is characterized by remarkably high parallelism and frequent historical contingency.

Evolution of pH-sensitive transcription termination during adaptation to repeated long-term starvation

Fluctuating environments that consist of regular cycles of co-occurring stress are a common challenge faced by cellular populations. For a population to thrive in constantly changing conditions, an ability to coordinate a rapid cellular response is essential. Here, we identify a mutation conferring an arginine-to-histidine (Arg to His) substitution in the transcription terminator Rho. The rho R109H mutation frequently arose in E. coli populations experimentally evolved under repeated long-term starvation conditions, during which feast and famine result in drastic environmental pH fluctuations. Metagenomic sequencing revealed that populations containing the rho mutation also possess putative loss-of-function mutations in ydcI, which encodes a recently characterized transcription factor associated with pH homeostasis. Genetic reconstructions of these mutations show that the rho allele confers a plastic alkaline-induced reduction of Rho function that, when found in tandem with a ΔydcI allele, leads to intracellular alkalinization and genetic assimilation of Rho mutant function. We further identify Arg to His substitutions at analogous sites in rho alleles from species originating from fluctuating alkaline environments. Our results suggest that Arg to His substitutions in global regulators of gene expression can serve to rapidly coordinate complex responses through pH sensing and shed light on how cellular populations across the tree of life use environmental cues to coordinate rapid responses to complex, fluctuating environments.

Linking genotypic and phenotypic changes in the E. coli long-term evolution experiment using metabolomics

Changes in an organism's environment, genome, or gene expression patterns can lead to changes in its metabolism. The metabolic phenotype can be under selection and contributes to adaptation. However, the networked and convoluted nature of an organism's metabolism makes relating mutations, metabolic changes, and effects on fitness challenging. To overcome this challenge, we use the long-term evolution experiment (LTEE) with E. coli as a model to understand how mutations can eventually affect metabolism and perhaps fitness. We used mass spectrometry to broadly survey the metabolomes of the ancestral strains and all 12 evolved lines. We combined this metabolic data with mutation and expression data to suggest how mutations that alter specific reaction pathways, such as the biosynthesis of nicotinamide adenine dinucleotide, might increase fitness in the system. Our work provides a better understanding of how mutations might affect fitness through the metabolic changes in the LTEE and thus provides a major step in developing a complete genotype-phenotype map for this experimental system.

Biochemical and genetic dissection of the RNA-binding surface of the FinO domain of Escherichia coli ProQ

RNA-binding proteins play important roles in bacterial gene regulation through interactions with both coding and noncoding RNAs. ProQ is a FinO-domain protein that binds a large set of RNAs in Escherichia coli, though the details of how ProQ binds these RNAs remain unclear. In this study, we used a combination of in vivo and in vitro binding assays to confirm key structural features of E. coli ProQ's FinO domain and explore its mechanism of RNA interactions. Using a bacterial three-hybrid assay, we performed forward genetic screens to confirm the importance of the concave face of ProQ in RNA binding. Using gel shift assays, we directly probed the contributions of ten amino acids on ProQ binding to seven RNA targets. Certain residues (R58, Y70, and R80) were found to be essential for binding of all seven RNAs, while substitutions of other residues (K54 and R62) caused more moderate binding defects. Interestingly, substitutions of two amino acids (K35, R69), which are evolutionarily variable but adjacent to conserved residues, showed varied effects on the binding of different RNAs; these may arise from the differing sequence context around each RNA's terminator hairpin. Together, this work confirms many of the essential RNA-binding residues in ProQ initially identified in vivo and supports a model in which residues on the conserved concave face of the FinO domain such as R58, Y70, and R80 form the main RNA-binding site of E. coli ProQ, while additional contacts contribute to the binding of certain RNAs.

Mutation bias and the predictability of evolution

Predicting evolutionary outcomes is an important research goal in a diversity of contexts. The focus of evolutionary forecasting is usually on adaptive processes, and efforts to improve prediction typically focus on selection. However, adaptive processes often rely on new mutations, which can be strongly influenced by predictable biases in mutation. Here, we provide an overview of existing theory and evidence for such mutation-biased adaptation and consider the implications of these results for the problem of prediction, in regard to topics such as the evolution of infectious diseases, resistance to biochemical agents, as well as cancer and other kinds of somatic evolution. We argue that empirical knowledge of mutational biases is likely to improve in the near future, and that this knowledge is readily applicable to the challenges of short-term prediction. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Biochemical and genetic dissection of the RNA-binding surface of the FinO domain of Escherichia coli ProQ

RNA-binding proteins play important roles in bacterial gene regulation through interactions with both coding and non-coding RNAs. ProQ is a FinO-domain protein that binds a large set of RNAs in Escherichia coli , though the details of how ProQ binds these RNAs remain unclear. In this study, we used a combination of in vivo and in vitro binding assays to confirm key structural features of E. coli ProQ's FinO domain and explore its mechanism of RNA interactions. Using a bacterial three-hybrid assay, we performed forward genetic screens to confirm the importance of the concave face of ProQ in RNA binding. Using gel shift assays, we directly probed the contributions of ten amino acids on ProQ binding to seven RNA targets. Certain residues (R58, Y70, and R80) were found to be essential for binding of all seven RNAs, while substitutions of other residues (K54 and R62) caused more moderate binding defects. Interestingly, substitutions of two amino acids (K35, R69), which are evolutionarily variable but adjacent to conserved residues, showed varied effects on the binding of different RNAs; these may arise from the differing sequence context around each RNA's terminator hairpin. Together, this work confirms many of the essential RNA-binding residues in ProQ initially identified in vivo and supports a model in which residues on the conserved concave face of the FinO domain such as R58, Y70 and R80 form the main RNA-binding site of E. coli ProQ, while additional contacts contribute to the binding of certain RNAs.

The RNA-Binding Protein ProQ Promotes Antibiotic Persistence in Salmonella

Bacterial populations can survive exposure to antibiotics through transient phenotypic and gene expression changes. These changes can be attributed to a small subpopulation of bacteria, giving rise to antibiotic persistence. Although this phenomenon has been known for decades, much remains to be learned about the mechanisms that drive persister formation. The RNA-binding protein ProQ has recently emerged as a global regulator of gene expression. Here, we show that ProQ impacts persister formation in Salmonella. In vitro, ProQ contributes to growth arrest in a subset of cells that are able to survive treatment at high concentrations of different antibiotics. The underlying mechanism for ProQ-dependent persister formation involves the activation of metabolically costly processes, including the flagellar pathway and the type III protein secretion system encoded on Salmonella pathogenicity island 2. Importantly, we show that the ProQ-dependent phenotype is relevant during macrophage infection and allows Salmonella to survive the combined action of host immune defenses and antibiotics. Together, our data highlight the importance of ProQ in Salmonella persistence and pathogenesis. IMPORTANCE Bacteria can avoid eradication by antibiotics through a phenomenon known as persistence. Persister cells arise through phenotypic heterogeneity and constitute a small fraction of dormant cells within a population of actively growing bacteria, which is susceptible to antibiotic killing. In this study, we show that ProQ, an RNA-binding protein and global regulator of gene expression, promotes persisters in the human pathogen Salmonella enterica serovar Typhimurium. Bacteria lacking the proQ gene outcompete wild-type bacteria under laboratory conditions, are less prone to enter growth dormancy, and form fewer persister cells. The basis for these phenotypes lies in ProQ's ability to activate energy-consuming cellular processes, including flagellar motility and protein secretion. Importantly, we show that ProQ contributes to the persister phenotype during Salmonella infection of macrophages, indicating an important role of this global regulator in Salmonella pathogenesis.

Rapid Adaptation Often Occurs through Mutations to the Most Highly Conserved Positions of the RNA Polymerase Core Enzyme

Mutations to the genes encoding the RNA polymerase core enzyme (RNAPC) and additional housekeeping regulatory genes were found to be involved in adaptation, in the context of numerous evolutionary experiments, in which bacteria were exposed to diverse selective pressures. This provides a conundrum, as the housekeeping genes that were so often mutated in response to these diverse selective pressures tend to be among the genes that are most conserved in their sequences across the bacterial phylogeny. In order to further examine this apparent discrepancy, we characterized the precise positions of the RNAPC involved in adaptation to a large variety of selective pressures. We found that RNAPC lab adaptations tended to occur at positions displaying traits associated with higher selective constraint. Specifically, compared to other RNAPC positions, positions involved in adaptation tended to be more conserved in their sequences within bacteria, were more often located within defined protein domains, and were located closer to the complex's active site. Higher sequence conservation was also found for resource exhaustion adaptations occurring within additional housekeeping genes. Combined, our results demonstrate that the positions that change most readily in response to well-defined selective pressures exerted in lab environments are often also those that evolve most slowly in nature.

Adaptation to novel spatially-structured environments is driven by the capsule and alters virulence-associated traits

The extracellular capsule is a major virulence factor, but its ubiquity in free-living bacteria with large environmental breadths suggests that it shapes adaptation to novel niches. Yet, how it does so, remains unexplored. Here, we evolve three Klebsiella strains and their capsule mutants in parallel. Their comparison reveals different phenotypic and genotypic evolutionary changes that alter virulence-associated traits. Non-capsulated populations accumulate mutations that reduce exopolysaccharide production and increase biofilm formation and yield, whereas most capsulated populations become hypermucoviscous, a signature of hypervirulence. Hence, adaptation to novel environments primarily occurs by fine-tuning expression of the capsular locus. The same evolutionary conditions selecting for mutations in the capsular gene wzc leading to hypermucoviscosity also result in increased susceptibility to antibiotics by mutations in the ramA regulon. This implies that general adaptive processes outside the host can affect capsule evolution and its role in virulence and infection outcomes may be a by-product of such adaptation.

Phenotypic and genotypic evolution in worrisome Klebsiella spp. is influenced by the capsule. Here the authors show that adaptation outside the host can impact virulence-associated traits, including de novo emergence of hypermucoviscosity.

Complex Ecotype Dynamics Evolve in Response to Fluctuating Resources

Ecotypic diversification and its associated cooperative behaviors are frequently observed in natural microbial populations whose access to resources is often sporadic. However, the extent to which fluctuations in resource availability influence the emergence of cooperative ecotypes is not fully understood. To determine how exposure to repeated resource limitation affects the establishment and long-term maintenance of ecotypes in a structured environment, we followed 32 populations of Escherichia coli evolving to either 1-day or 10-day feast/famine cycles for 900 days. Population-level analysis revealed that compared to populations evolving to 1-day cycles, 10-day populations evolved increased biofilm density, higher parallelism in mutational targets, and increased mutation rates. As previous investigations of evolution in structured environments have identified biofilm formation as the earliest observable phenotype associated with diversification of ecotypes, we revived cultures midway through the evolutionary process and conducted additional genomic, transcriptional, and phenotypic analyses of clones isolated from these evolving populations. We found not only that 10-day feast/famine cycles support multiple ecotypes but also that these ecotypes exhibit cooperative behavior. Consistent with the black queen hypothesis, or evolution of cooperation by gene loss, transcriptomic evidence suggests the evolution of bidirectional cross-feeding behaviors based on essential resources. These results provide insight into how analogous cooperative relationships may emerge in natural microbial communities. IMPORTANCE Despite regular feast and famine conditions representing an environmental pressure that is commonly encountered by microbial communities, the evolutionary outcomes of repeated cycles of feast and famine have been less studied. By experimentally evolving initially isogenic Escherichia coli populations to 10-day feast/famine cycles, we observed rapid diversification into ecotypes with evidence of bidirectional cross-feeding on costly resources and frequency-dependent fitness. Although unidirectional cross-feeding has been repeatedly observed to evolve in laboratory culture, most investigations of bidirectional cooperative behaviors in microbial populations have been conducted in engineered communities. This work demonstrates the de novo evolution of black queen relationships in a microbial population originating from a single ancestor, providing a model for investigation of the eco-evolutionary processes leading to mutualistic cooperation.

The role of sigma factor competition in bacterial adaptation under prolonged starvation

The study of adaptive microbial evolution in the laboratory can illuminate the genetic mechanisms of gaining fitness under a pre-defined set of selection factors. Laboratory evolution of bacteria under long-term starvation has gained importance in recent years because of its ability to uncover adaptive strategies that overcome prolonged nutrient limitation, a condition often encountered by natural microbes. In this evolutionary paradigm, bacteria are maintained in an energy-restricted environment in a growth phase called long-term stationary phase (LTSP). This phase is characterized by a stable, viable population size and highly dynamic genetic changes. Multiple independent iterations of LTSP evolution experiments have given rise to mutants that are slow-growing compared to the ancestor. Although the antagonistic regulation between rapid growth and the stress response is well-known in bacteria (especially Escherichia coli), the growth deficit of many LTSP-adapted mutants has not been explored in detail. In this review, I pinpoint the trade-off between growth and stress response as a dominant driver of evolutionary strategies under prolonged starvation. Focusing on mainly E. coli-based research, I discuss the various affectors and regulators of the competition between sigma factors to occupy their targets on the genome, and assess its effect on growth advantage in stationary phase (GASP). Finally, I comment on some crucial issues that hinder the progress of the field, including identification of novel metabolites in nutrient-depleted media, and the importance of using multidisciplinary research to resolve them.

Mutation bias shapes the spectrum of adaptive substitutions

How do mutational biases influence the process of adaptation? A common assumption is that selection alone determines the course of adaptation from abundant preexisting variation. Yet, theoretical work shows broad conditions under which the mutation rate to a given type of variant strongly influences its probability of contributing to adaptation. Here we introduce a statistical approach to analyzing how mutation shapes protein sequence adaptation. Using large datasets from three different species, we show that the mutation spectrum has a proportional influence on the types of changes fixed in adaptation. We also show via computer simulations that a variety of factors can influence how closely the spectrum of adaptive substitutions reflects the spectrum of variants introduced by mutation.

Evolutionary adaptation often occurs by the fixation of beneficial mutations. This mode of adaptation can be characterized quantitatively by a spectrum of adaptive substitutions, i.e., a distribution for types of changes fixed in adaptation. Recent work establishes that the changes involved in adaptation reflect common types of mutations, raising the question of how strongly the mutation spectrum shapes the spectrum of adaptive substitutions. We address this question with a codon-based model for the spectrum of adaptive amino acid substitutions, applied to three large datasets covering thousands of amino acid changes identified in natural and experimental adaptation in Saccharomyces cerevisiae, Escherichia coli, and Mycobacterium tuberculosis. Using species-specific mutation spectra based on prior knowledge, we find that the mutation spectrum has a proportional influence on the spectrum of adaptive substitutions in all three species. Indeed, we find that by inferring the mutation rates that best explain the spectrum of adaptive substitutions, we can accurately recover the species-specific mutation spectra. However, we also find that the predictive power of the model differs substantially between the three species. To better understand these differences, we use population simulations to explore the factors that influence how closely the spectrum of adaptive substitutions mirrors the mutation spectrum. The results show that the influence of the mutation spectrum decreases with increasing mutational supply (Nμ) and that predictive power is strongly affected by the number and diversity of beneficial mutations.

Primordial mimicry induces morphological change in Escherichia coli

The morphology of primitive cells has been the subject of extensive research. A spherical form was commonly presumed in prebiotic studies but lacked experimental evidence in living cells. Whether and how the shape of living cells changed are unclear. Here we exposed the rod-shaped bacterium Escherichia coli to a resource utilization regime mimicking a primordial environment. Oleate was given as an easy-to-use model prebiotic nutrient, as fatty acid vesicles were likely present on the prebiotic Earth and might have been used as an energy resource. Six evolutionary lineages were generated under glucose-free but oleic acid vesicle (OAV)-rich conditions. Intriguingly, fitness increase was commonly associated with the morphological change from rod to sphere and the decreases in both the size and the area-to-volume ratio of the cell. The changed cell shape was conserved in either OAVs or glucose, regardless of the trade-offs in carbon utilization and protein abundance. Highly differentiated mutations present in the genome revealed two distinct strategies of adaption to OAV-rich conditions, i.e., either directly targeting the cell wall or not. The change in cell morphology of Escherichia coli for adapting to fatty acid availability supports the assumption of the primitive spherical form.

Lu et al. investigate the evolution of the shape of living cells by generating six experimental lineages of the rod-shaped E. coli under glucose-free conditions in the presence of oleic acid mimicking a primordial environment. The authors show that the morphological changes from rod to sphere accompanied fitness increases and adaptation amongst fatty acid availability supports the assumption of a primitive spherical form.

Microbial population dynamics and evolutionary outcomes under extreme energy limitation

Microorganisms commonly inhabit energy-limited ecosystems where cellular maintenance and reproduction is highly constrained. To gain insight into how individuals persist under such conditions, we derived demographic parameters from a collection of 21 heterotrophic bacterial taxa by censusing 100 populations in an effectively closed system for 1,000 d. All but one taxon survived prolonged resource scarcity, yielding estimated times to extinction ranging over four orders of magnitude from 100 to 105 y. Our findings corroborate reports of long-lived bacteria recovered from ancient environmental samples, while providing insight into mechanisms of persistence. As death rates declined over time, lifespan was extended through the scavenging of dead cells. Although reproduction was suppressed in the absence of exogenous resources, populations continued to evolve. Hundreds of mutations were acquired, contributing to genome-wide signatures of purifying selection as well as molecular signals of adaptation. Consistent ecological and evolutionary dynamics indicate that distantly related bacteria respond to energy limitation in a similar and predictable manner, which likely contributes to the stability and robustness of microbial life.

Endurance of extremely prolonged nutrient prevention across kingdoms of life

Numerous observations demonstrate that microorganisms can survive very long periods of nutrient deprivation and starvation. Moreover, it is evident that prolonged periods of starvation are a feature of many habitats, and many cells in all kingdoms of life are found in prolonged starvation conditions. Bacteria exhibit a range of responses to long-term starvation. These include genetic adaptations such as the long-term stationary phase and the growth advantage in stationary phase phenotypes characterized by mutations in stress-signaling genes and elevated mutation rates. Here, we suggest using the term "endurance of prolonged nutrient prevention" (EPNP phase), to describe this phase, which was also recently described in eukaryotes. Here, we review this literature and describe the current knowledge about the adaptations to very long-term starvation conditions in bacteria and eukaryotes, its conceptual and structural conservation across all kingdoms of life, and point out possible directions that merit further research.