Bacterial resistance response and resource availability mediate viral coexistence.

Viruses that infect bacteria, known as bacteriophages or phages, are the most prevalent entities on Earth. Their genetic diversity in nature is well documented, and members of divergent lineages can be found sharing the same ecological niche. This viral diversity can be influenced by a number of factors, including productivity, spatial structuring of the environment, and host-range trade-offs. Rapid evolution is also known to promote diversity by buffering ecological systems from extinction. There is, however, little known about the impact of coevolution on the maintenance of viral diversity within a microbial community. To address this, we developed a 4 species experimental system where two bacterial hosts, a generalist and a specialist phage, coevolved in a spatially homogenous environment over time. We observed the persistence of both viruses if the resource availability was sufficiently high. This coexistence occurred in the absence of any detectable host-range trade-offs that are costly for generalists and thus known to promote viral diversity. However, the coexistence was lost if two bacteria were not permitted to evolve alongside the phages or if two phages coevolved with a single bacterial host. Our findings indicate that a host's resistance response in mixed-species communities plays a significant role in maintaining viral diversity in the environment.

Ribosome-binding antibiotics increase bacterial longevity and growth efficiency.

Antibiotics, by definition, reduce bacterial growth rates in optimal culture conditions; however, the real-world environments bacteria inhabit see rapid growth punctuated by periods of low nutrient availability. How antibiotics mediate population decline during these periods is poorly understood. Bacteria cannot optimize for all environmental conditions because a growth-longevity tradeoff predicts faster growth results in faster population decline, and since bacteriostatic antibiotics slow growth, they should also mediate longevity. We quantify how antibiotics, their targets, and resistance mechanisms influence longevity using populations of Escherichia coli and, as the tradeoff predicts, populations are maintained for longer if they encounter ribosome-binding antibiotics doxycycline and erythromycin, a finding that is not observed using antibiotics with alternative cellular targets. This tradeoff also predicts resistance mechanisms that increase growth rates during antibiotic treatment could be detrimental during nutrient stresses, and indeed, we find resistance by ribosomal protection removes benefits to longevity provided by doxycycline. We therefore liken ribosomal protection to a "Trojan horse" because it provides protection from an antibiotic but, during nutrient stresses, it promotes the demise of the bacteria. Seeking mechanisms to support these observations, we show doxycycline promotes efficient metabolism and reduces the concentration of reactive oxygen species. Seeking generality, we sought another mechanism that affects longevity and we found the number of doxycycline targets, namely, the ribosomal RNA operons, mediates growth and longevity even without antibiotics. We conclude that slow growth, as observed during antibiotic treatment, can help bacteria overcome later periods of nutrient stress.

Canonical Host-Pathogen Trade-Offs Subverted by Mutations with Dual Benefits.

AbstractHost-parasite coevolution is expected to drive the evolution of genetic diversity because the traits used in arms races-namely, host range and parasite resistance-are hypothesized to trade off with traits used in resource competition. We therefore tested data for several trade-offs among 93 isolates of bacteriophage λ and 51 Escherichia coli genotypes that coevolved during a laboratory experiment. Surprisingly, we found multiple trade-ups (positive trait correlations) but little evidence of several canonical trade-offs. For example, some bacterial genotypes evaded a trade-off between phage resistance and absolute fitness, instead evolving simultaneous improvements in both traits. This was surprising because our experimental design was predicted to expose resistance-fitness trade-offs by culturing E. coli in a medium where the phage receptor, LamB, is also used for nutrient acquisition. On reflection, LamB mediates not one but many trade-offs, allowing for more complex trait interactions than just pairwise trade-offs. Here, we report that mathematical reasoning and laboratory data highlight how trade-ups should exist whenever an evolutionary system exhibits multiple interacting trade-offs. Does this mean that coevolution should not promote genetic diversity? No, quite the contrary. We deduce that whenever positive trait correlations are observed in multidimensional traits, other traits may trade off and so provide the right circumstances for diversity maintenance. Overall, this study reveals that there are predictive limits when data account only for pairwise trait correlations, and it argues that a wider range of circumstances than previously anticipated can promote genetic and species diversity.

Metabolic efficiency reshapes the seminal relationship between pathogen growth rate and virulence.

A cornerstone of classical virulence evolution theories is the assumption that pathogen growth rate is positively correlated with virulence, the amount of damage pathogens inflict on their hosts. Such theories are key for incorporating evolutionary principles into sustainable disease management strategies. Yet, empirical evidence raises doubts over this central assumption underpinning classical theories, thus undermining their generality and predictive power. In this paper, we identify a key component missing from current theories which redefines the growth-virulence relationship in a way that is consistent with data. By modifying the activity of a single metabolic gene, we engineered strains of Magnaporthe oryzae with different nutrient acquisition and growth rates. We conducted in planta infection studies and uncovered an unexpected non-monotonic relationship between growth rate and virulence that is jointly shaped by how growth rate and metabolic efficiency interact. This novel mechanistic framework paves the way for a much-needed new suite of virulence evolution theories.

Analysis of Pneumocystis Transcription Factor Evolution and Implications for Biology and Lifestyle.

Pneumocystis jirovecii kills hundreds of thousands of immunocompromised patients each year. Yet many aspects of the biology of this obligate pathogen remain obscure because it is not possible to culture the fungus in vitro independently of its host. Consequently, our understanding of Pneumocystis pathobiology is heavily reliant upon bioinformatic inferences. We have exploited a powerful combination of genomic and phylogenetic approaches to examine the evolution of transcription factors in Pneumocystis species. We selected protein families (Pfam families) that correspond to transcriptional regulators and used bioinformatic approaches to compare these families in the seven Pneumocystis species that have been sequenced to date with those from other yeasts, other human and plant pathogens, and other obligate parasites. Some Pfam families of transcription factors have undergone significant reduction during their evolution in the Pneumocystis genus, and other Pfam families have been lost or appear to be in the process of being lost. Meanwhile, other transcription factor families have been retained in Pneumocystis species, and some even appear to have undergone expansion. On this basis, Pneumocystis species seem to have retained transcriptional regulators that control chromosome maintenance, ribosomal gene regulation, RNA processing and modification, and respiration. Meanwhile, regulators that promote the assimilation of alternative carbon sources, amino acid, lipid, and sterol biosynthesis, and iron sensing and homeostasis appear to have been lost. Our analyses of transcription factor retention, loss, and gain provide important insights into the biology and lifestyle of Pneumocystis. IMPORTANCE Pneumocystis jirovecii is a major fungal pathogen of humans that infects healthy individuals, colonizing the lungs of infants. In immunocompromised and transplant patients, this fungus causes life-threatening pneumonia, and these Pneumocystis infections remain among the most common and serious infections in HIV/AIDS patients. Yet we remain remarkably ignorant about the biology and epidemiology of Pneumocystis due to the inability to culture this fungus in vitro. Our analyses of transcription factor retentions, losses, and gains in sequenced Pneumocystis species provide valuable new views of their specialized biology, suggesting the retention of many metabolic and stress regulators and the loss of others that are essential in free-living fungi. Given the lack of in vitro culture methods for Pneumocystis, this powerful bioinformatic approach has advanced our understanding of the lifestyle of P. jirovecii and the nature of its dependence on the host for survival.

Seeking patterns of antibiotic resistance in ATLAS, an open, raw MIC database with patient metadata.

Antibiotic resistance represents a growing medical concern where raw, clinical datasets are under-exploited as a means to track the scale of the problem. We therefore sought patterns of antibiotic resistance in the Antimicrobial Testing Leadership and Surveillance (ATLAS) database. ATLAS holds 6.5M minimal inhibitory concentrations (MICs) for 3,919 pathogen-antibiotic pairs isolated from 633k patients in 70 countries between 2004 and 2017. We show most pairs form coherent, although not stationary, timeseries whose frequencies of resistance are higher than other databases, although we identified no systematic bias towards including more resistant strains in ATLAS. We sought data anomalies whereby MICs could shift for methodological and not clinical or microbiological reasons and found artefacts in over 100 pathogen-antibiotic pairs. Using an information-optimal clustering methodology to classify pathogens into low and high antibiotic susceptibilities, we used ATLAS to predict changes in resistance. Dynamics of the latter exhibit complex patterns with MIC increases, and some decreases, whereby subpopulations' MICs can diverge. We also identify pathogens at risk of developing clinical resistance in the near future.

Would that it were so simple: Interactions between multiple traits undermine classical single-trait-based predictions of microbial community function and evolution.

Understanding how microbial traits affect the evolution and functioning of microbial communities is fundamental for improving the management of harmful microorganisms, while promoting those that are beneficial. Decades of evolutionary ecology research has focused on examining microbial cooperation, diversity, productivity and virulence but with one crucial limitation. The traits under consideration, such as public good production and resistance to antibiotics or predation, are often assumed to act in isolation. Yet, in reality, multiple traits frequently interact, which can lead to unexpected and undesired outcomes for the health of macroorganisms and ecosystem functioning. This is because many predictions generated in a single-trait context aimed at promoting diversity, reducing virulence or controlling antibiotic resistance can fail for systems where multiple traits interact. Here, we provide a much needed discussion and synthesis of the most recent research to reveal the widespread and diverse nature of multi-trait interactions and their consequences for predicting and controlling microbial community dynamics. Importantly, we argue that synthetic microbial communities and multi-trait mathematical models are powerful tools for managing the beneficial and detrimental impacts of microbial communities, such that past mistakes, like those made regarding the stewardship of antimicrobials, are not repeated.

The Antibiotic Dosage of Fastest Resistance Evolution: Gene Amplifications Underpinning the Inverted-U.

To determine the dosage at which antibiotic resistance evolution is most rapid, we treated Escherichia coli in vitro, deploying the antibiotic erythromycin at dosages ranging from zero to high. Adaptation was fastest just below erythromycin's minimal inhibitory concentration (MIC) and genotype-phenotype correlations determined from whole genome sequencing revealed the molecular basis: simultaneous selection for copy number variation in three resistance mechanisms which exhibited an "inverted-U" pattern of dose-dependence, as did several insertion sequences and an integron. Many genes did not conform to this pattern, however, reflecting changes in selection as dose increased: putative media adaptation polymorphisms at zero antibiotic dosage gave way to drug target (ribosomal RNA operon) amplification at mid dosages whereas prophage-mediated drug efflux amplifications dominated at the highest dosages. All treatments exhibited E. coli increases in the copy number of efflux operons acrAB and emrE at rates that correlated with increases in population density. For strains where the inverted-U was no longer observed following the genetic manipulation of acrAB, it could be recovered by prolonging the antibiotic treatment at subMIC dosages.

Predicting microbial growth dynamics in response to nutrient availability.

Developing mathematical models to accurately predict microbial growth dynamics remains a key challenge in ecology, evolution, biotechnology, and public health. To reproduce and grow, microbes need to take up essential nutrients from the environment, and mathematical models classically assume that the nutrient uptake rate is a saturating function of the nutrient concentration. In nature, microbes experience different levels of nutrient availability at all environmental scales, yet parameters shaping the nutrient uptake function are commonly estimated for a single initial nutrient concentration. This hampers the models from accurately capturing microbial dynamics when the environmental conditions change. To address this problem, we conduct growth experiments for a range of micro-organisms, including human fungal pathogens, baker's yeast, and common coliform bacteria, and uncover the following patterns. We observed that the maximal nutrient uptake rate and biomass yield were both decreasing functions of initial nutrient concentration. While a functional form for the relationship between biomass yield and initial nutrient concentration has been previously derived from first metabolic principles, here we also derive the form of the relationship between maximal nutrient uptake rate and initial nutrient concentration. Incorporating these two functions into a model of microbial growth allows for variable growth parameters and enables us to substantially improve predictions for microbial dynamics in a range of initial nutrient concentrations, compared to keeping growth parameters fixed.

Evolution of drug-resistant and virulent small colonies in phenotypically diverse populations of the human fungal pathogen Candida glabrata.

Antimicrobial resistance frequently carries a fitness cost to a pathogen, measured as a reduction in growth rate compared to the sensitive wild-type, in the absence of antibiotics. Existing empirical evidence points to the following relationship between cost of resistance and virulence. If a resistant pathogen suffers a fitness cost in terms of reduced growth rate it commonly has lower virulence compared to the sensitive wild-type. If this cost is absent so is the reduction in virulence. Here we show, using experimental evolution of drug resistance in the fungal human pathogen Candida glabrata, that reduced growth rate of resistant strains need not result in reduced virulence. Phenotypically heterogeneous populations were evolved in parallel containing highly resistant sub-population small colony variants (SCVs) alongside sensitive sub-populations. Despite their low growth rate in the absence of an antifungal drug, the SCVs did not suffer a marked alteration in virulence compared with the wild-type ancestral strain, or their co-isolated sensitive strains. This contrasts with classical theory that assumes growth rate to positively correlate with virulence. Our work thus highlights the complexity of the relationship between resistance, basic life-history traits and virulence.

Predicting community dynamics of antibiotic-sensitive and -resistant species in fluctuating environments.

Microbes occupy almost every niche within and on their human hosts. Whether colonizing the gut, mouth or bloodstream, microorganisms face temporal fluctuations in resources and stressors within their niche but we still know little of how environmental fluctuations mediate certain microbial phenotypes, notably antimicrobial-resistant ones. For instance, do rapid or slow fluctuations in nutrient and antimicrobial concentrations select for, or against, resistance? We tackle this question using an ecological approach by studying the dynamics of a synthetic and pathogenic microbial community containing two species, one sensitive and the other resistant to an antibiotic drug where the community is exposed to different rates of environmental fluctuation. We provide mathematical models, supported by experimental data, to demonstrate that simple community outcomes, such as competitive exclusion, can shift to coexistence and ecosystem bistability as fluctuation rates vary. Theory gives mechanistic insight into how these dynamical regimes are related. Importantly, our approach highlights a fundamental difference between resistance in single-species populations, the context in which it is usually assayed, and that in communities. While fast environmental changes are known to select against resistance in single-species populations, here we show that they can promote the resistant species in mixed-species communities. Our theoretical observations are verified empirically using a two-species Candida community.

Privatization of public goods can cause population decline.

Microbes commonly deploy a risky strategy to acquire nutrients from their environment, involving the production of costly public goods that can be exploited by neighbouring individuals. Why engage in such a strategy when an exploitation-free alternative is readily available whereby public goods are kept private? We address this by examining metabolism of Saccharomyces cerevisiae in its native form and by creating a new three-strain synthetic community deploying different strategies of sucrose metabolism. Public-metabolizers digest resources externally, private-metabolizers internalize resources before digestion, and cheats avoid the metabolic costs of digestion but exploit external products generated by competitors. A combination of mathematical modelling and ecological experiments reveal that private-metabolizers invade and take over an otherwise stable community of public-metabolizers and cheats. However, owing to the reduced growth rate of private-metabolizers and population bottlenecks that are frequently associated with microbial communities, privatizing public goods can become unsustainable, leading to population decline.

Author Correction: Drug-mediated metabolic tipping between antibiotic resistant states in a mixed-species community.

In the version of this Article originally published, the following sentence was missing from the Acknowledgements: "R.E.B. is an EPSRC Healthcare Technologies Impact Fellow EP/N033671/1; I.G. is funded by ERC Consolidator grant 647292 MathModExp; A.J.P.B., N.A.R.G. and A.T. were funded by BBSRC grant BB/F00513X/1; K.H., I.G., S.N. and E.C. were funded by BBSRC grant BB/F005210/2." This text has now been added.

Drug-mediated metabolic tipping between antibiotic resistant states in a mixed-species community.

Microbes rarely exist in isolation, rather, they form intricate multi-species communities that colonize our bodies and inserted medical devices. However, the efficacy of antimicrobials is measured in clinical laboratories exclusively using microbial monocultures. Here, to determine how multi-species interactions mediate selection for resistance during antibiotic treatment, particularly following drug withdrawal, we study a laboratory community consisting of two microbial pathogens. Single-species dose responses are a poor predictor of community dynamics during treatment so, to better understand those dynamics, we introduce the concept of a dose-response mosaic, a multi-dimensional map that indicates how species' abundance is affected by changes in abiotic conditions. We study the dose-response mosaic of a two-species community with a 'Gene × Gene × Environment × Environment' ecological interaction whereby Candida glabrata, which is resistant to the antifungal drug fluconazole, competes for survival with Candida albicans, which is susceptible to fluconazole. The mosaic comprises several zones that delineate abiotic conditions where each species dominates. Zones are separated by loci of bifurcations and tipping points that identify what environmental changes can trigger the loss of either species. Observations of the laboratory communities corroborated theory, showing that changes in both antibiotic concentration and nutrient availability can push populations beyond tipping points, thus creating irreversible shifts in community composition from drug-sensitive to drug-resistant species. This has an important consequence: resistant species can increase in frequency even if an antibiotic is withdrawn because, unwittingly, a tipping point was passed during treatment.

When increasing population density can promote the evolution of metabolic cooperation.

Microbial cooperation drives ecological and epidemiological processes and is affected by the ecology and demography of populations. Population density influences the selection for cooperation, with spatial structure and the type of social dilemma, namely public-goods production or self-restraint, shaping the outcome. While existing theories predict that in spatially structured environments increasing population density can select either for or against cooperation, experimental studies with both public-goods production and self-restraint systems have only ever shown that increasing population density favours cheats. We suggest that the disparity between theory and empirical studies results from experimental procedures not capturing environmental conditions predicted by existing theories to influence the outcome. Our study resolves this issue and provides the first experimental evidence that high population density can favour cooperation in spatially structured environments for both self-restraint and public-goods production systems. Moreover, using a multi-trait mathematical model supported by laboratory experiments we extend this result to systems where the self-restraint and public-goods social dilemmas interact. We thus provide a systematic understanding of how the strength of interaction between the two social dilemmas and the degree of spatial structure within an environment affect selection for cooperation. These findings help to close the current gap between theory and experiments.

The unconstrained evolution of fast and efficient antibiotic-resistant bacterial genomes.

Evolutionary trajectories are constrained by trade-offs when mutations that benefit one life history trait incur fitness costs in other traits. As resistance to tetracycline antibiotics by increased efflux can be associated with an increase in length of the Escherichia coli chromosome of 10% or more, we sought costs of resistance associated with doxycycline. However, it was difficult to identify any because the growth rate (r), carrying capacity (K) and drug efflux rate of E. coli increased during evolutionary experiments where the species was exposed to doxycycline. Moreover, these improvements remained following drug withdrawal. We sought mechanisms for this seemingly unconstrained adaptation, particularly as these traits ought to trade-off according to rK selection theory. Using prokaryote and eukaryote microorganisms, including clinical pathogens, we show that r and K can trade-off, but need not, because of 'rK trade-ups'. r and K trade-off only in sufficiently carbon-rich environments where growth is inefficient. We then used E. coli ribosomal RNA (rRNA) knockouts to determine specific mutations, namely changes in rRNA operon (rrn) copy number, than can simultaneously maximize r and K. The optimal genome has fewer operons, and therefore fewer functional ribosomes, than the ancestral strain. It is, therefore, unsurprising for r-adaptation in the presence of a ribosome-inhibiting antibiotic, doxycycline, to also increase population size. We found two costs for this improvement: an elongated lag phase and the loss of stress protection genes.

Density-Dependent Recycling Promotes the Long-Term Survival of Bacterial Populations during Periods of Starvation.

The amount of natural resources in the Earth's environment is in flux, which can trigger catastrophic collapses of ecosystems. How populations survive under nutrient-poor conditions is a central question in ecology. Curiously, some bacteria persist for a long time in nutrient-poor environments. Although this survival may be accomplished through cell death and the recycling of dead cells, the importance of these processes and the mechanisms underlying the survival of the populations have not been quantitated. Here, we use microbial laboratory experiments and mathematical models to demonstrate that death and recycling are essential activities for the maintenance of cell survival. We also show that the behavior of the survivors is governed by population density feedback, wherein growth is limited not only by the available resources but also by the population density. The numerical simulations suggest that population density-dependent recycling could be an advantageous behavior under starvation conditions. IMPORTANCE: How organisms survive after exhaustion of resources is a central question in ecology. Starving Escherichia coli constitute a model system to understand survival mechanisms during long-term starvation. Although death and the recycling of dead cells might play a key role in the maintenance of long-term survival, their mechanisms and importance have not been quantitated. Here, we verified the significance of social recycling of dead cells for long-term survival. We also show that the survivors restrained their recycling and did not use all available nutrients released from dead cells, which may be advantageous under starvation conditions. These results indicate that not only the utilization of dead cells but also restrained recycling coordinate the effective utilization of limited resources for long-term survival under starvation.

Harbouring public good mutants within a pathogen population can increase both fitness and virulence.

Existing theory, empirical, clinical and field research all predict that reducing the virulence of individuals within a pathogen population will reduce the overall virulence, rendering disease less severe. Here, we show that this seemingly successful disease management strategy can fail with devastating consequences for infected hosts. We deploy cooperation theory and a novel synthetic system involving the rice blast fungus Magnaporthe oryzae. In vivo infections of rice demonstrate that M. oryzae virulence is enhanced, quite paradoxically, when a public good mutant is present in a population of high-virulence pathogens. We reason that during infection, the fungus engages in multiple cooperative acts to exploit host resources. We establish a multi-trait cooperation model which suggests that the observed failure of the virulence reduction strategy is caused by the interference between different social traits. Multi-trait cooperative interactions are widespread, so we caution against the indiscriminant application of anti-virulence therapy as a disease-management strategy.

Stability of Cross-Feeding Polymorphisms in Microbial Communities.

Cross-feeding, a relationship wherein one organism consumes metabolites excreted by another, is a ubiquitous feature of natural and clinically-relevant microbial communities and could be a key factor promoting diversity in extreme and/or nutrient-poor environments. However, it remains unclear how readily cross-feeding interactions form, and therefore our ability to predict their emergence is limited. In this paper we developed a mathematical model parameterized using data from the biochemistry and ecology of an E. coli cross-feeding laboratory system. The model accurately captures short-term dynamics of the two competitors that have been observed empirically and we use it to systematically explore the stability of cross-feeding interactions for a range of environmental conditions. We find that our simple system can display complex dynamics including multi-stable behavior separated by a critical point. Therefore whether cross-feeding interactions form depends on the complex interplay between density and frequency of the competitors as well as on the concentration of resources in the environment. Moreover, we find that subtly different environmental conditions can lead to dramatically different results regarding the establishment of cross-feeding, which could explain the apparently unpredictable between-population differences in experimental outcomes. We argue that mathematical models are essential tools for disentangling the complexities of cross-feeding interactions.

Kinase Inhibition Leads to Hormesis in a Dual Phosphorylation-Dephosphorylation Cycle.

Many antimicrobial and anti-tumour drugs elicit hormetic responses characterised by low-dose stimulation and high-dose inhibition. While this can have profound consequences for human health, with low drug concentrations actually stimulating pathogen or tumour growth, the mechanistic understanding behind such responses is still lacking. We propose a novel, simple but general mechanism that could give rise to hormesis in systems where an inhibitor acts on an enzyme. At its core is one of the basic building blocks in intracellular signalling, the dual phosphorylation-dephosphorylation motif, found in diverse regulatory processes including control of cell proliferation and programmed cell death. Our analytically-derived conditions for observing hormesis provide clues as to why this mechanism has not been previously identified. Current mathematical models regularly make simplifying assumptions that lack empirical support but inadvertently preclude the observation of hormesis. In addition, due to the inherent population heterogeneities, the presence of hormesis is likely to be masked in empirical population-level studies. Therefore, examining hormetic responses at single-cell level coupled with improved mathematical models could substantially enhance detection and mechanistic understanding of hormesis.