Trade-offs, trade-ups, and high mutational parallelism underlie microbial adaptation during extreme cycles of feast and famine.

Microbes are evolutionarily robust organisms capable of rapid adaptation to complex stress, which enables them to colonize harsh environments. In nature, microbes are regularly challenged by starvation, which is a particularly complex stress because resource limitation often co-occurs with changes in pH, osmolarity, and toxin accumulation created by metabolic waste. Often overlooked are the additional complications introduced by eventual resource replenishment, as successful microbes must withstand rapid environmental shifts before swiftly capitalizing on replenished resources to avoid invasion by competing species. To understand how microbes navigate trade-offs between growth and survival, ultimately adapting to thrive in environments with extreme fluctuations, we experimentally evolved 16 Escherichia coli populations for 900 days in repeated feast/famine conditions with cycles of 100-day starvation before resource replenishment. Using longitudinal population-genomic analysis, we found that evolution in response to extreme feast/famine is characterized by narrow adaptive trajectories with high mutational parallelism and notable mutational order. Genetic reconstructions reveal that early mutations result in trade-offs for biofilm and motility but trade-ups for growth and survival, as these mutations conferred positively correlated advantages during both short-term and long-term culture. Our results demonstrate how microbes can navigate the adaptive landscapes of regularly fluctuating conditions and ultimately follow mutational trajectories that confer benefits across diverse environments.

Bacterial adenine cross-feeding stems from a purine salvage bottleneck.

Diverse ecosystems host microbial relationships that are stabilized by nutrient cross-feeding. Cross-feeding can involve metabolites that should hold value for the producer. Externalization of such communally valuable metabolites is often unexpected and difficult to predict. Previously, we discovered purine externalization by Rhodopseudomonas palustris by its ability to rescue an Escherichia coli purine auxotroph. Here we found that an E. coli purine auxotroph can stably coexist with R. palustris due to purine cross-feeding. We identified the cross-fed purine as adenine. Adenine was externalized by R. palustris under diverse growth conditions. Computational modeling suggested that adenine externalization occurs via diffusion across the cytoplasmic membrane. RNAseq analysis led us to hypothesize that adenine accumulation and externalization stem from a salvage pathway bottleneck at the enzyme encoded by apt. Ectopic expression of apt eliminated adenine externalization, supporting our hypothesis. A comparison of 49 R. palustris strains suggested that purine externalization is relatively common, with 16 strains exhibiting the trait. Purine externalization was correlated with the genomic orientation of apt, but apt orientation alone could not always explain purine externalization. Our results provide a mechanistic understanding of how a communally valuable metabolite can participate in cross-feeding. Our findings also highlight the challenge in identifying genetic signatures for metabolite externalization.

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.

Differential adenine methylation analysis reveals increased variability in 6mA in the absence of methyl-directed mismatch repair.

IMPORTANCE: Methylation greatly influences the bacterial genome by guiding DNA repair and regulating pathogenic and stress-response phenotypes. But, the rate of epigenetic changes and their consequences on molecular phenotypes are underexplored. Through a detailed characterization of genome-wide adenine methylation in a commonly used laboratory strain of Escherichia coli, we reveal that mismatch repair deficient populations experience an increase in epimutations resulting in a genome-wide reduction of 6mA methylation in a manner consistent with genetic drift. Our findings highlight how methylation patterns evolve and the constraints on epigenetic evolution due to post-replicative DNA repair, contributing to a deeper understanding of bacterial genome evolution and how epimutations may introduce semi-permanent variation that can influence adaptation.

Case Studies in the Assessment of Microbial Fitness: Seemingly Subtle Changes Can Have Major Effects on Phenotypic Outcomes.

Following the completion of an adaptive evolution experiment, fitness evaluations are routinely conducted to assess the magnitude of adaptation. In doing so, proper consideration should be given when determining the appropriate methods as trade-offs may exist between accuracy and throughput. Here, we present three instances in which small changes in the framework or execution of fitness evaluations significantly impacted the outcomes. The first case illustrates that discrepancies in fitness conclusions can arise depending on the approach to evaluating fitness, the culture vessel used, and the sampling method. The second case reveals that variations in environmental conditions can occur associated with culture vessel material. Specifically, these subtle changes can greatly affect microbial physiology leading to changes in the culture pH and distorting fitness measurements. Finally, the last case reports that heterogeneity in CFU formation time can result in inaccurate fitness conclusions. Based on each case, considerations and recommendations are presented for future adaptive evolution experiments.

Rapid evolution of mutation rate and spectrum in response to environmental and population-genetic challenges.

Ecological and demographic factors can significantly shape the evolution of microbial populations both directly and indirectly, as when changes in the effective population size affect the efficiency of natural selection on the mutation rate. However, it remains unclear how rapidly the mutation-rate responds evolutionarily to the entanglement of ecological and population-genetic factors over time. Here, we directly assess the mutation rate and spectrum of Escherichia coli clones isolated from populations evolving in response to 1000 days of different transfer volumes and resource-replenishment intervals. The evolution of mutation rates proceeded rapidly in response to demographic and/or environmental changes, with substantial bidirectional shifts observed as early as 59 generations. These results highlight the remarkable rapidity by which mutation rates are shaped in asexual lineages in response to environmental and population-genetic forces, and are broadly consistent with the drift-barrier hypothesis for the evolution of mutation rates, while also highlighting situations in which mutator genotypes may be promoted by positive selection.

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.

Evolutionary Dynamics of Asexual Hypermutators Adapting to a Novel Environment.

How microbes adapt to a novel environment is a central question in evolutionary biology. Although adaptive evolution must be fueled by beneficial mutations, whether higher mutation rates facilitate the rate of adaptive evolution remains unclear. To address this question, we cultured Escherichia coli hypermutating populations, in which a defective methyl-directed mismatch repair pathway causes a 140-fold increase in single-nucleotide mutation rates. In parallel with wild-type E. coli, populations were cultured in tubes containing Luria-Bertani broth, a complex medium known to promote the evolution of subpopulation structure. After 900 days of evolution, in three transfer schemes with different population-size bottlenecks, hypermutators always exhibited similar levels of improved fitness as controls. Fluctuation tests revealed that the mutation rates of hypermutator lines converged evolutionarily on those of wild-type populations, which may have contributed to the absence of fitness differences. Further genome-sequence analysis revealed that, although hypermutator populations have higher rates of genomic evolution, this largely reflects strong genetic linkage. Despite these linkage effects, the evolved population exhibits parallelism in fixed mutations, including those potentially related to biofilm formation, transcription regulation, and mutation-rate evolution. Together, these results are generally inconsistent with a hypothesized positive relationship between the mutation rate and the adaptive speed of evolution, and provide insight into how clonal adaptation occurs in novel environments.

Multi-omic Characterization of Intraspecies Variation in Laboratory and Natural Environments.

Investigation of microbial communities has led to many advances in our understanding of ecosystem function, whether that ecosystem is a subglacial lake or the human gut. Within these communities, much emphasis has been placed on interspecific variation and between-species relationships. However, with current advances in sequencing technology resulting in both the reduction in sequencing costs and the rise of shotgun metagenomic sequencing, the importance of intraspecific variation and within-species relationships is becoming realized. Our group conducts multi-omic analyses to understand how spatial structure and resource availability influence diversification within a species and the potential for long-term coexistence of multiple ecotypes within a microbial community. Here, we present examples of ecotypic variation observed in the lab and in the wild, current challenges faced when investigating intraspecies diversity, and future developments that we expect to define the field over the next 5 years.

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.

Enhanced nutrient uptake is sufficient to drive emergent cross-feeding between bacteria in a synthetic community.

Interactive microbial communities are ubiquitous, influencing biogeochemical cycles and host health. One widespread interaction is nutrient exchange, or cross-feeding, wherein metabolites are transferred between microbes. Some cross-fed metabolites, such as vitamins, amino acids, and ammonium (NH4+), are communally valuable and impose a cost on the producer. The mechanisms that enforce cross-feeding of communally valuable metabolites are not fully understood. Previously we engineered a cross-feeding coculture between N2-fixing Rhodopseudomonas palustris and fermentative Escherichia coli. Engineered R. palustris excretes essential nitrogen as NH4+ to E. coli, while E. coli excretes essential carbon as fermentation products to R. palustris. Here, we sought to determine whether a reciprocal cross-feeding relationship would evolve spontaneously in cocultures with wild-type R. palustris, which is not known to excrete NH4+. Indeed, we observed the emergence of NH4+ cross-feeding, but driven by adaptation of E. coli alone. A missense mutation in E. coli NtrC, a regulator of nitrogen scavenging, resulted in constitutive activation of an NH4+ transporter. This activity likely allowed E. coli to subsist on the small amount of leaked NH4+ and better reciprocate through elevated excretion of fermentation products from a larger E. coli population. Our results indicate that enhanced nutrient uptake by recipients, rather than increased excretion by producers, is an underappreciated yet possibly prevalent mechanism by which cross-feeding can emerge.

Evaluation of a fluorescence resonance energy transfer quantitative polymerase chain reaction assay for identification of gyrA mutations conferring enrofloxacin resistance in canine urinary Escherichia coli isolates and canine urine specimens.

OBJECTIVE To evaluate a fluorescence resonance energy transfer quantitative PCR (FRET-qPCR) assay for detection of gyrA mutations conferring fluoroquinolone resistance in canine urinary Escherichia coli isolates and canine urine specimens. SAMPLE 264 canine urinary E coli isolates and 283 clinical canine urine specimens. PROCEDURES The E coli isolates were used to validate the FRET-qPCR assay. Urine specimens were evaluated by bacterial culture and identification, isolate enrofloxacin susceptibility testing, and FRET-qPCR assay. Sensitivity and specificity of the FRET-qPCR assay for detection of gyrA mutations in urine specimens and in E coli isolated from urine specimens were computed, with results of enrofloxacin susceptibility testing used as the reference standard. RESULTS The validated FRET-qPCR assay discriminated between enrofloxacin-resistant and enrofloxacin-susceptible E coli isolates with an area under the receiver operating characteristic curve of 0.92. The assay accurately identified 25 of 40 urine specimens as containing enrofloxacin-resistant isolates (sensitivity, 62.5%) and 226 of 243 urine specimens as containing enrofloxacin-susceptible isolates (specificity, 93.0%). When the same assay was performed on E coli isolates recovered from these specimens, sensitivity (77.8%) and specificity (94.8%) increased. Moderate agreement was achieved between results of the FRET-qPCR assay and enrofloxacin susceptibility testing for E coli isolates recovered from urine specimens. CONCLUSIONS AND CLINICAL RELEVANCE The FRET-qPCR assay was able to rapidly distinguish between enrofloxacin-resistant and enrofloxacin-susceptible E coli in canine clinical urine specimens through detection of gyrA mutations. Therefore, the assay may be useful in clinical settings to screen such specimens for enrofloxacin-resistant E coli to avoid inappropriate use of enrofloxacin and contributing to antimicrobial resistance.

An Escherichia coli Nitrogen Starvation Response Is Important for Mutualistic Coexistence with Rhodopseudomonas palustris.

Microbial mutualistic cross-feeding interactions are ubiquitous and can drive important community functions. Engaging in cross-feeding undoubtedly affects the physiology and metabolism of individual species involved. However, the nature in which an individual species' physiology is influenced by cross-feeding and the importance of those physiological changes for the mutualism have received little attention. We previously developed a genetically tractable coculture to study bacterial mutualisms. The coculture consists of fermentative Escherichia coli and phototrophic Rhodopseudomonas palustris In this coculture, E. coli anaerobically ferments sugars into excreted organic acids as a carbon source for R. palustris In return, a genetically engineered R. palustris strain constitutively converts N2 into NH4+, providing E. coli with essential nitrogen. Using transcriptome sequencing (RNA-seq) and proteomics, we identified transcript and protein levels that differ in each partner when grown in coculture versus monoculture. When in coculture with R. palustris, E. coli gene expression changes resembled a nitrogen starvation response under the control of the transcriptional regulator NtrC. By genetically disrupting E. coli NtrC, we determined that a nitrogen starvation response is important for a stable coexistence, especially at low R. palustris NH4+ excretion levels. Destabilization of the nitrogen starvation regulatory network resulted in variable growth trends and, in some cases, extinction. Our results highlight that alternative physiological states can be important for survival within cooperative cross-feeding relationships.IMPORTANCE Mutualistic cross-feeding between microbes within multispecies communities is widespread. Studying how mutualistic interactions influence the physiology of each species involved is important for understanding how mutualisms function and persist in both natural and applied settings. Using a bacterial mutualism consisting of Rhodopseudomonas palustris and Escherichia coli growing cooperatively through bidirectional nutrient exchange, we determined that an E. coli nitrogen starvation response is important for maintaining a stable coexistence. The lack of an E. coli nitrogen starvation response ultimately destabilized the mutualism and, in some cases, led to community collapse after serial transfers. Our findings thus inform on the potential necessity of an alternative physiological state for mutualistic coexistence with another species compared to the physiology of species grown in isolation.

Escherichia coli cultures maintain stable subpopulation structure during long-term evolution.

How genetic variation is generated and maintained remains a central question in evolutionary biology. When presented with a complex environment, microbes can take advantage of genetic variation to exploit new niches. Here we present a massively parallel experiment where WT and repair-deficient (∆mutL) Escherichia coli populations have evolved over 3 y in a spatially heterogeneous and nutritionally complex environment. Metagenomic sequencing revealed that these initially isogenic populations evolved and maintained stable subpopulation structure in just 10 mL of medium for up to 10,000 generations, consisting of up to five major haplotypes with many minor haplotypes. We characterized the genomic, transcriptomic, exometabolomic, and phenotypic differences between clonal isolates, revealing subpopulation structure driven primarily by spatial segregation followed by differential utilization of nutrients. In addition to genes regulating the import and catabolism of nutrients, major polymorphisms of note included insertion elements transposing into fimE (regulator of the type I fimbriae) and upstream of hns (global regulator of environmental-change and stress-response genes), both known to regulate biofilm formation. Interestingly, these genes have also been identified as critical to colonization in uropathogenic E. coli infections. Our findings illustrate the complexity that can arise and persist even in small cultures, raising the possibility that infections may often be promoted by an evolving and complex pathogen population.

Selection on Position of Nonsense Codons in Introns.

Introns occasionally remain in mature messenger RNAs (mRNAs) due to splicing errors and the translated, aberrant proteins that result represent a metabolic cost and may have other deleterious consequences. The nonsense-mediated decay (NMD) pathway degrades aberrant mRNAs, which it recognizes by the presence of an in-frame premature termination codon (PTC). We investigated whether selection has shaped the location of PTCs in introns to reduce waste and facilitate NMD. We found across seven model organisms, that in both first and last introns, PTCs occur earlier in introns than expected by chance, suggesting that selection favors earlier position. This pattern is more pronounced in species with larger effective population sizes. The pattern does not hold for last introns in the two mammal species, however, perhaps because in these species NMD is not initiated from 3'-terminal introns. We conclude that there is compelling evidence that the location of PTCs is shaped by selection for reduced waste and efficient degradation of aberrant mRNAs.

The repeatability of genome-wide mutation rate and spectrum estimates.

Over the last decade, mutation studies have grown in popularity due to the affordability and accessibility of whole genome sequencing. As the number of species in which spontaneous mutation has been directly estimated approaches 20 across two domains of life, questions arise over the repeatability of results in such experiments. Five species were identified in which duplicate mutation studies have been performed. Across these studies the difference in estimated spontaneous mutation rate is at most, weakly significant (p < 0.01). However, a highly significant (p < 10(-5)), threefold difference in the rate of insertions/deletions (indels) exists between two recent studies in Schizosaccharomyces pombe. Upon investigation of the ancestral genome sequence for both studies, a possible anti-mutator allele was identified. The observed variation in indel rate may imply that the use of indel markers, such as microsatellites, for the investigation of genetic diversity within and among populations may be inappropriate because of the assumption of uniform mutation rate within a species.

Similar Mutation Rates but Highly Diverse Mutation Spectra in Ascomycete and Basidiomycete Yeasts.

Yeast species are extremely diverse and not monophyletic. Because the majority of yeast research focuses on ascomycetes, the mutational determinants of genetic diversity across yeast species are not well understood. By combining mutation-accumulation techniques with whole-genome sequencing, we resolved the genomic mutation rate and spectrum of the oleaginous (oil-producing) 'red yeast' Rhodotorula toruloides, the first such study in the fungal phylum Basidiomycota. We find that the mutation spectrum is quite different from what has been observed in all other studied unicellular eukaryotes, but similar to that in most bacteria­a predominance of transitions relative to transversions. Rhodotorula toruloides has a significantly higher A:T→G:C transition rate­possibly elevated by the abundant flanking G/C nucleotides in the GC-rich genome, as well as a much lower G:C→T:A transversion rate. In spite of these striking differences, there are substantial consistencies between R. toruloides and the ascomycete model yeasts: a spontaneous base-substitution mutation rate of 1.90 × 10 −10 per site per cell division as well as an elevated mutation rate at non-methylated 5'CpG3' sites. These results imply the evolution of variable mutation spectra in the face of similar mutation rates in yeasts.

Genome-Wide Estimates of Mutation Rates and Spectrum in Schizosaccharomyces pombe Indicate CpG Sites are Highly Mutagenic Despite the Absence of DNA Methylation.

We accumulated mutations for 1952 generations in 79 initially identical, haploid lines of the fission yeast Schizosaccharomyces pombe, and then performed whole-genome sequencing to determine the mutation rates and spectrum. We captured 696 spontaneous mutations across the 79 mutation accumulation (MA) lines. We compared the mutation spectrum and rate to a recently published equivalent experiment on the same species, and to another model ascomycetous yeast, the budding yeast Saccharomyces cerevisiae. While the two species are approximately 600 million years diverged from each other, they share similar life histories, genome size and genomic G/C content. We found that Sc. pombe and S. cerevisiae have similar mutation rates, but Sc. pombe exhibits a stronger insertion bias. Intriguingly, we observed an increased mutation rate at cytosine nucleotides, specifically CpG nucleotides, which is also seen in S. cerevisiae. However, the absence of methylation in Sc. pombe and the pattern of mutation at these sites, primarily C → A as opposed to C → T, strongly suggest that the increased mutation rate is not caused by deamination of methylated cytosines. This result implies that the high mutability of CpG dinucleotides in other species may be caused in part by a methylation-independent mechanism. Many of our findings mirror those seen in the recent study, despite the use of different passaging conditions, indicating that MA is a reliable method for estimating mutation rates and spectra.