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

Site-selective modifications by lipid A phosphoethanolamine transferases linked to colistin resistance and bacterial fitness

Genes encoding lipid A modifying phosphoethanolamine transferases (PETs) are genetically diverse and can confer resistance to colistin and antimicrobial peptides. To better understand the functional diversity of PETs, we characterized three canonical mobile colistin resistance (mcr) alleles (mcr-1, -3, -9), one intrinsic pet (eptA), and two mcr-like genes (petB, petC) in Escherichia coli. Using an isogenic expression system, we show that mcr-1 and mcr-3 confer similar phenotypes of decreased colistin susceptibility with low fitness costs. mcr-9, which is phylogenetically closely related to mcr-3, and eptA only provide fitness advantages in the presence of sub-inhibitory concentrations of colistin and significantly reduce fitness in media without colistin. PET-B and PET-C were phenotypically distinct from bonafide PETs; neither impacted colistin susceptibility nor caused considerable fitness cost. Strikingly, we found for the first time that different PETs selectively modify different phosphates of lipid A; MCR-1, MCR-3, and PET-C selectively modify the 4'-phosphate, whereas MCR-9 and EptA modify the 1-phosphate. However, 4'-phosphate modifications facilitated by MCR-1 and -3 are associated with lowered colistin susceptibility and low toxicity. Our results suggest that PETs have a wide phenotypic diversity and that increased colistin resistance is associated with specific lipid A modification patterns that have been largely unexplored thus far.

IMPORTANCE

Rising levels of resistance to increasing numbers of antimicrobials have led to the revival of last resort antibiotic colistin. Unfortunately, resistance to colistin is also spreading in the form of mcr genes, making it essential to (i) improve the identification of resistant bacteria to allow clinicians to prescribe effective drug regimens and (ii) develop new combination therapies effective at targeting resistant bacteria. Our results demonstrate that PETs, including MCR variants, are site-selective in Escherichia coli and that site-selectivity correlates with the level of susceptibility and fitness costs conferred by certain PETs. Site selectivity associated with a given PET may not only help predict colistin resistance phenotypes but may also provide an avenue to (i) improve drug regimens and (ii) develop new combination therapies to better combat colistin-resistant bacteria.

Nutrient-rich spatial refuges buffer against extinction and promote evolutionary rescue in evolving microbial populations

Microbial populations are often exposed to long-term abiotic disturbances, which can reduce population viability and cause local extinction. Eco-evolutionary theory suggests that spatial refuges can facilitate persistence and evolutionary rescue. However, one drawback of spatial refuges is reduced exposure to nutrients such as carbon and oxygen, suggesting the protective effect of refuges depends on the interplay between environmental conditions and the degree of stress. Here, we test this general idea using mathematical modelling, and experimental evolution of the model bacterium Pseudomonas fluorescens SBW25 under salinity stress. As our model predicted, we find that the ability of spatial refuges to rescue evolving populations from extinction crucially depends on nutrient availability. Populations evolving under salinity stress where nutrient-rich spatial refuges were available, harboured clones that displayed enhanced salt resistance, indicating that nutrient-rich spatial refuges can facilitate evolutionary rescue. Furthermore, while control-salinity-evolved populations adapted to spatial structure by evolving enhanced motility (likely through parallel mutations in PFLU_4551, a predicted aerotaxis response regulator), this phenotype was constrained under high salinity, because increased motility negates the benefits of a spatial refuge. Our results reveal a general interplay between spatial refuges and nutrient availability that could be leveraged to reduce extinction risk in natural populations.

Single-cell heterogeneity in ribosome content and the consequences for the growth laws

Across species and environments, the ribosome content of cell populations correlates with population growth rate. The robustness and universality of this correlation have led to its classification as a "growth law." This law has fueled theories about how evolution selects for microbial organisms that maximize their growth rate based on nutrient availability, and it has informed models about how individual cells regulate their growth rates and ribosomal content. However, due to methodological limitations, this growth law has rarely been studied at the level of individual cells. While populations of fast-growing cells tend to have more ribosomes than populations of slow-growing cells, it is unclear whether individual cells tightly regulate their ribosome content to match their environment. Here, we employ recent groundbreaking single-cell RNA sequencing techniques to study this growth law at the single-cell level in two different microbes, S. cerevisiae (a single-celled yeast and eukaryote) and B. subtilis (a bacterium and prokaryote). In both species, we observe significant variation in the ribosomal content of single cells that is not predictive of growth rate. Fast-growing populations include cells exhibiting transcriptional signatures of slow growth and stress, as do cells with the highest ribosome content we survey. Broadening our focus to non-ribosomal transcripts reveals subpopulations of cells in unique transcriptional states suggestive that they have evolved to do things other than maximize their rate of growth. Overall, these results indicate that single-cell ribosome levels are not finely tuned to match population growth rates or nutrient availability and cannot be predicted by a Gaussian process model that assumes measurements are sampled from a normal distribution centered on the population average. This work encourages the expansion of growth law and other models that predict how growth rates are regulated or how they evolve to consider single-cell heterogeneity. To this end, we provide extensive data and analysis of ribosomal and transcriptomic variation across thousands of single cells from multiple conditions, replicates, and species.

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.

In Vitro Resistance-Predicting Studies and In Vitro Resistance-Related Parameters-A Hit-to-Lead Perspective

Despite the urgent need for new antibiotics, very few innovative antibiotics have recently entered clinics or clinical trials. To provide a constant supply of new drug candidates optimized in terms of their potential to select for resistance in natural settings, in vitro resistance-predicting studies need to be improved and scaled up. In this review, the following in vitro parameters are presented: frequency of spontaneous mutant selection (FSMS), mutant prevention concentration (MPC), dominant mutant prevention concentration (MPC-D), inferior-mutant prevention concentration (MPC-F), and minimal selective concentration (MSC). The utility of various adaptive laboratory evolution (ALE) approaches (serial transfer, continuous culture, and evolution in spatiotemporal microenvironments) for comparing hits in terms of the level and time required for multistep resistance to emerge is discussed. We also consider how the hit-to-lead stage can benefit from high-throughput ALE setups based on robotic workstations, do-it-yourself (DIY) continuous cultivation systems, microbial evolution and growth arena (MEGA) plates, soft agar gradient evolution (SAGE) plates, microfluidic chips, or microdroplet technology. Finally, approaches for evaluating the fitness of in vitro-generated resistant mutants are presented. This review aims to draw attention to newly emerged ideas on how to improve the in vitro forecasting of the potential of compounds to select for resistance in natural settings.

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.

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.

Extreme Sensitivity of Fitness to Environmental Conditions: Lessons from #1BigBatch

The phrase "survival of the fittest" has become an iconic descriptor of how natural selection works. And yet, precisely measuring fitness, even for single-celled microbial populations growing in controlled laboratory conditions, remains a challenge. While numerous methods exist to perform these measurements, including recently developed methods utilizing DNA barcodes, all methods are limited in their precision to differentiate strains with small fitness differences. In this study, we rule out some major sources of imprecision, but still find that fitness measurements vary substantially from replicate to replicate. Our data suggest that very subtle and difficult to avoid environmental differences between replicates create systematic variation across fitness measurements. We conclude by discussing how fitness measurements should be interpreted given their extreme environment dependence. This work was inspired by the scientific community who followed us and gave us tips as we live tweeted a high-replicate fitness measurement experiment at #1BigBatch.

Identifying Targets of Selection in Laboratory Evolution Experiments

Adaptive evolution navigates a balance between chance and determinism. The stochastic processes of mutation and drift generate phenotypic variation; however, once mutations reach an appreciable frequency in the population, their fate is governed by the deterministic action of selection, enriching for favorable genotypes and purging the less-favorable ones. The net result is that replicate populations will traverse similar-but not identical-pathways to higher fitness. This parallelism in evolutionary outcomes can be leveraged to identify the genes and pathways under selection. However, distinguishing between beneficial and neutral mutations is challenging because many beneficial mutations will be lost due to drift and clonal interference, and many neutral (and even deleterious) mutations will fix by hitchhiking. Here, we review the best practices that our laboratory uses to identify genetic targets of selection from next-generation sequencing data of evolved yeast populations. The general principles for identifying the mutations driving adaptation will apply more broadly.