Antibiotic collateral sensitivity is contingent on the repeatability of evolution

Distinguishing mutants that resist drugs via different mechanisms by examining fitness tradeoffs

There is growing interest in designing multidrug therapies that leverage tradeoffs to combat resistance. Tradeoffs are common in evolution and occur when, for example, resistance to one drug results in sensitivity to another. Major questions remain about the extent to which tradeoffs are reliable, specifically, whether the mutants that provide resistance to a given drug all suffer similar tradeoffs. This question is difficult because the drug-resistant mutants observed in the clinic, and even those evolved in controlled laboratory settings, are often biased towards those that provide large fitness benefits. Thus, the mutations (and mechanisms) that provide drug resistance may be more diverse than current data suggests. Here, we perform evolution experiments utilizing lineage-tracking to capture a fuller spectrum of mutations that give yeast cells a fitness advantage in fluconazole, a common antifungal drug. We then quantify fitness tradeoffs for each of 774 evolved mutants across 12 environments, finding these mutants group into 6 classes with characteristically different tradeoffs. Their unique tradeoffs may imply that each group of mutants affects fitness through different underlying mechanisms. Some of the groupings we find are surprising. For example, we find some mutants that resist single drugs do not resist their combination, while others do. And some mutants to the same gene have different tradeoffs than others. These findings, on one hand, demonstrate the difficulty in relying on consistent or intuitive tradeoffs when designing multidrug treatments. On the other hand, by demonstrating that hundreds of adaptive mutations can be reduced to a few groups with characteristic tradeoffs, our findings may yet empower multidrug strategies that leverage tradeoffs to combat resistance. More generally speaking, by grouping mutants that likely affect fitness through similar underlying mechanisms, our work guides efforts to map the phenotypic effects of mutation.

Evolutionary druggability for low-dimensional fitness landscapes toward new metrics for antimicrobial applications

The term 'druggability' describes the molecular properties of drugs or targets in pharmacological interventions and is commonly used in work involving drug development for clinical applications. There are no current analogues for this notion that quantify the drug-target interaction with respect to a given target variant's sensitivity across a breadth of drugs in a panel, or a given drug's range of effectiveness across alleles of a target protein. Using data from low-dimensional empirical fitness landscapes composed of 16 β-lactamase alleles and 7 β-lactam drugs, we introduce two metrics that capture (i) the average susceptibility of an allelic variant of a drug target to any available drug in a given panel ('variant vulnerability'), and (ii) the average applicability of a drug (or mixture) across allelic variants of a drug target ('drug applicability'). Finally, we (iii) disentangle the quality and magnitude of interactions between loci in the drug target and the seven drug environments in terms of their mutation by mutation by environment (G x G x E) interactions, offering mechanistic insight into the variant variability and drug applicability metrics. Summarizing, we propose that our framework can be applied to other datasets and pathogen-drug systems to understand which pathogen variants in a clinical setting are the most concerning (low variant vulnerability), and which drugs in a panel are most likely to be effective in an infection defined by standing genetic variation in the pathogen drug target (high drug applicability).

Threshold-awareness in adaptive cancer therapy

Although adaptive cancer therapy shows promise in integrating evolutionary dynamics into treatment scheduling, the stochastic nature of cancer evolution has seldom been taken into account. Various sources of random perturbations can impact the evolution of heterogeneous tumors, making performance metrics of any treatment policy random as well. In this paper, we propose an efficient method for selecting optimal adaptive treatment policies under randomly evolving tumor dynamics. The goal is to improve the cumulative "cost" of treatment, a combination of the total amount of drugs used and the total treatment time. As this cost also becomes random in any stochastic setting, we maximize the probability of reaching the treatment goals (tumor stabilization or eradication) without exceeding a pre-specified cost threshold (or a "budget"). We use a novel Stochastic Optimal Control formulation and Dynamic Programming to find such "threshold-aware" optimal treatment policies. Our approach enables an efficient algorithm to compute these policies for a range of threshold values simultaneously. Compared to treatment plans shown to be optimal in a deterministic setting, the new "threshold-aware" policies significantly improve the chances of the therapy succeeding under the budget, which is correlated with a lower general drug usage. We illustrate this method using two specific examples, but our approach is far more general and provides a new tool for optimizing adaptive therapies based on a broad range of stochastic cancer models.

A stochastic programming approach to the antibiotics time machine problem

Antibiotics Time Machine is an important problem to understand antibiotic resistance and how it can be reversed. Mathematically, it can be modeled as follows: Consider a set of genotypes, each of which contain a set of mutated and unmutated genes. Suppose that a set of growth rate measurements of each genotype under a set of antibiotics is given. The transition probabilities of a 'realization' of a Markov chain associated with each arc under each antibiotic are computable via a predefined function given the growth rate realizations. The aim is to maximize the expected probability of reaching to the genotype with all unmutated genes given the initial genotype in a predetermined number of transitions, considering the following two sources of uncertainties: (i) the randomness in growth rates, (ii) the randomness in transition probabilities, which are functions of growth rates. We develop stochastic mixed-integer linear programming and dynamic programming approaches to solve static and dynamic versions of the Antibiotics Time Machine Problem under the aforementioned uncertainties. We adapt a Sample Average Approximation approach that exploits the special structure of the problem and provide accurate solutions that perform very well in an out-of-sample analysis.

Reinforcement learning informs optimal treatment strategies to limit antibiotic resistance

Antimicrobial resistance was estimated to be associated with 4.95 million deaths worldwide in 2019. It is possible to frame the antimicrobial resistance problem as a feedback-control problem. If we could optimize this feedback-control problem and translate our findings to the clinic, we could slow, prevent, or reverse the development of high-level drug resistance. Prior work on this topic has relied on systems where the exact dynamics and parameters were known a priori. In this study, we extend this work using a reinforcement learning (RL) approach capable of learning effective drug cycling policies in a system defined by empirically measured fitness landscapes. Crucially, we show that it is possible to learn effective drug cycling policies despite the problems of noisy, limited, or delayed measurement. Given access to a panel of 15 [Formula: see text]-lactam antibiotics with which to treat the simulated Escherichia coli population, we demonstrate that RL agents outperform two naive treatment paradigms at minimizing the population fitness over time. We also show that RL agents approach the performance of the optimal drug cycling policy. Even when stochastic noise is introduced to the measurements of population fitness, we show that RL agents are capable of maintaining evolving populations at lower growth rates compared to controls. We further tested our approach in arbitrary fitness landscapes of up to 1,024 genotypes. We show that minimization of population fitness using drug cycles is not limited by increasing genome size. Our work represents a proof-of-concept for using AI to control complex evolutionary processes.

In Vitro Microevolution and Co-Selection Assessment of Amoxicillin and Cefotaxime Impact on Escherichia coli Resistance Development

The global spread of antimicrobial resistance has become a prominent issue in both veterinary and public health in the 21st century. The extensive use of amoxicillin, a beta-lactam antibiotic, and consequent resistance development are particularly alarming in food-producing animals, with a focus on the swine and poultry sectors. Another beta-lactam, cefotaxime, is widely utilized in human medicine, where the escalating resistance to third- and fourth-generation cephalosporins is a major concern. The aim of this study was to simulate the development of phenotypic and genotypic resistance to beta-lactam antibiotics, focusing on amoxicillin and cefotaxime. The investigation of the minimal inhibitory concentrations (MIC) of antibiotics was performed at 1×, 10×, 100×, and 1000× concentrations using the modified microbial evolution and growth arena (MEGA-plate) method. Our results indicate that amoxicillin significantly increased the MIC values of several tested antibiotics, except for oxytetracycline and florfenicol. In the case of cefotaxime, this increase was observed in all classes. A total of 44 antimicrobial resistance genes were identified in all samples. Chromosomal point mutations, particularly concerning cefotaxime, revealed numerous complex mutations, deletions, insertions, and single nucleotide polymorphisms (SNPs) that were not experienced in the case of amoxicillin. The findings suggest that, regarding amoxicillin, the point mutation of the acrB gene could explain the observed MIC value increases due to the heightened activity of the acrAB-tolC efflux pump system. However, under the influence of cefotaxime, more intricate processes occurred, including complex amino acid substitutions in the ampC gene promoter region, increased enzyme production induced by amino acid substitutions and SNPs, as well as mutations in the acrR and robA repressor genes that heightened the activity of the acrAB-tolC efflux pump system. These changes may contribute to the significant MIC increases observed for all tested antibiotics. The results underscore the importance of understanding cross-resistance development between individual drugs when choosing clinical alternative drugs. The point mutations in the mdtB and emrR genes may also contribute to the increased activity of the mdtABC-tolC and emrAB-tolC pump systems against all tested antibiotics. The exceptionally high mutation rate induced by cephalosporins justifies further investigations to clarify the exact mechanism behind.

The evolution of diverse antimicrobial responses in vancomycin-intermediate Staphylococcus aureus and its therapeutic implications

Staphylococcus aureus causes endocarditis, osteomyelitis, and bacteremia. Clinicians often prescribe vancomycin as an empiric therapy to account for methicillin-resistant S. aureus (MRSA) and narrow treatment based on culture susceptibility results. However, these results reflect a single time point before empiric treatment and represent a limited subset of the total bacterial population within the patient. Thus, while they may indicate that the infection is susceptible to a particular drug, this recommendation may no longer be accurate during therapy. Here, we addressed how antibiotic susceptibility changes over time by accounting for evolution. We evolved 18 methicillin-susceptible S. aureus (MSSA) populations under increasing vancomycin concentrations until they reached intermediate resistance levels. Sequencing revealed parallel mutations that affect cell membrane stress response and cell-wall biosynthesis. The populations exhibited repeated cross-resistance to daptomycin and varied responses to meropenem, gentamicin, and nafcillin. We accounted for this variability by deriving likelihood estimates that express a population's probability of exhibiting a drug response following vancomycin treatment. Our results suggest antistaphylococcal penicillins are preferable first-line treatments for MSSA infections but also highlight the inherent uncertainty that evolution poses to effective therapies. Infections may take varied evolutionary paths; therefore, considering evolution as a probabilistic process should inform our therapeutic choices.

Sub-inhibitory antibiotic treatment selects for enhanced metabolic efficiency

Bacterial growth and metabolic rates are often closely related. However, under antibiotic selection, a paradox in this relationship arises: antibiotic efficacy decreases when bacteria are metabolically dormant, yet antibiotics select for resistant cells that grow fastest during treatment. That is, antibiotic selection counterintuitively favors bacteria with fast growth but slow metabolism. Despite this apparent contradiction, antibiotic resistant cells have historically been characterized primarily in the context of growth, whereas the extent of analogous changes in metabolism is comparatively unknown. Here, we observed that previously evolved antibiotic-resistant strains exhibited a unique relationship between growth and metabolism whereby nutrient utilization became more efficient, regardless of the growth rate. To better understand this unexpected phenomenon, we used a simplified model to simulate bacterial populations adapting to sub-inhibitory antibiotic selection through successive bottlenecking events. Simulations predicted that sub-inhibitory bactericidal antibiotic concentrations could select for enhanced metabolic efficiency, defined based on nutrient utilization: drug-adapted cells are able to achieve the same biomass while utilizing less substrate, even in the absence of treatment. Moreover, simulations predicted that restoring metabolic efficiency would re-sensitize resistant bacteria exhibiting metabolic-dependent resistance; we confirmed this result using adaptive laboratory evolutions of Escherichia coli under carbenicillin treatment. Overall, these results indicate that metabolic efficiency is under direct selective pressure during antibiotic treatment and that differences in evolutionary context may determine both the efficacy of different antibiotics and corresponding re-sensitization approaches.IMPORTANCEThe sustained emergence of antibiotic-resistant pathogens combined with the stalled drug discovery pipelines highlights the critical need to better understand the underlying evolution mechanisms of antibiotic resistance. To this end, bacterial growth and metabolic rates are often closely related, and resistant cells have historically been characterized exclusively in the context of growth. However, under antibiotic selection, antibiotics counterintuitively favor cells with fast growth, and slow metabolism. Through an integrated approach of mathematical modeling and experiments, this study thereby addresses the significant knowledge gap of whether antibiotic selection drives changes in metabolism that complement, and/or act independently, of antibiotic resistance phenotypes.

Diverse mutant selection windows shape spatial heterogeneity in evolving populations

Mutant selection windows (MSWs), the range of drug concentrations that select for drug-resistant mutants, have long been used as a model for predicting drug resistance and designing optimal dosing strategies in infectious disease. The canonical MSW model offers comparisons between two subtypes at a time: drug-sensitive and drug-resistant. In contrast, the fitness landscape model with N alleles, which maps genotype to fitness, allows comparisons between N genotypes simultaneously, but does not encode continuous drug response data. In clinical settings, there may be a wide range of drug concentrations selecting for a variety of genotypes in both cancer and infectious diseases. Therefore, there is a need for a more robust model of the pathogen response to therapy to predict resistance and design new therapeutic approaches. Fitness seascapes, which model genotype-by-environment interactions, permit multiple MSW comparisons simultaneously by encoding genotype-specific dose-response data. By comparing dose-response curves, one can visualize the range of drug concentrations where one genotype is selected over another. In this work, we show how N-allele fitness seascapes allow for N * 2N-1 unique MSW comparisons. In spatial drug diffusion models, we demonstrate how fitness seascapes reveal spatially heterogeneous MSWs, extending the MSW model to more fully reflect the selection of drug resistant genotypes. Furthermore, using synthetic data and empirical dose-response data in cancer, we find that the spatial structure of MSWs shapes the evolution of drug resistance in an agent-based model. By simulating a tumor treated with cyclic drug therapy, we find that mutant selection windows introduced by drug diffusion promote the proliferation of drug resistant cells. Our work highlights the importance and utility of considering dose-dependent fitness seascapes in evolutionary medicine.

Global surveillance of antimicrobial resistance in food animals using priority drugs maps

Antimicrobial resistance (AMR) in food animals is a growing threat to animal health and potentially to human health. In resource-limited settings, allocating resources to address AMR can be guided with maps. Here, we mapped AMR prevalence in 7 antimicrobials in Escherichia coli and nontyphoidal Salmonella species across low- and middle-income countries (LIMCs), using 1088 point-prevalence surveys in combination with a geospatial model. Hotspots of AMR were predicted in China, India, Brazil, Chile, and part of central Asia and southeastern Africa. The highest resistance prevalence was for tetracycline (59% for E. coli and 54% for nontyphoidal Salmonella, average across LMICs) and lowest for cefotaxime (33% and 19%). We also identified the antimicrobial with the highest probability of resistance exceeding critical levels (50%) in the future (1.7-12.4 years) for each 10 × 10 km pixel on the map. In Africa and South America, 78% locations were associated with penicillins or tetracyclines crossing 50% resistance in the future. In contrast, in Asia, 77% locations were associated with penicillins or sulphonamides. Our maps highlight diverging geographic trends of AMR prevalence across antimicrobial classes, and can be used to target AMR surveillance in AMR hotspots for priority antimicrobial classes.

Laboratory Evolution of Antimicrobial Resistance in Bacteria to Develop Rational Treatment Strategies

Laboratory evolution studies, particularly with Escherichia coli, have yielded invaluable insights into the mechanisms of antimicrobial resistance (AMR). Recent investigations have illuminated that, with repetitive antibiotic exposures, bacterial populations will adapt and eventually become tolerant and resistant to the drugs. Through intensive analyses, these inquiries have unveiled instances of convergent evolution across diverse antibiotics, the pleiotropic effects of resistance mutations, and the role played by loss-of-function mutations in the evolutionary landscape. Moreover, a quantitative analysis of multidrug combinations has shed light on collateral sensitivity, revealing specific drug combinations capable of suppressing the acquisition of resistance. This review article introduces the methodologies employed in the laboratory evolution of AMR in bacteria and presents recent discoveries concerning AMR mechanisms derived from laboratory evolution. Additionally, the review outlines the application of laboratory evolution in endeavors to formulate rational treatment strategies.

The Pseudomonas aeruginosa Resistome: Permanent and Transient Antibiotic Resistance, an Overview

One of the most concerning characteristics of Pseudomonas aeruginosa is its low susceptibility to several antibiotics of common use in clinics, as well as its facility to acquire increased resistance levels. Consequently, the study of the antibiotic resistance mechanisms of this bacterium is of relevance for human health. For such a study, different types of resistance should be distinguished. The intrinsic resistome is composed of a set of genes, present in the core genome of P. aeruginosa, which contributes to its characteristic, species-specific, phenotype of susceptibility to antibiotics. Acquired resistance refers to those genetic events, such as the acquisition of mutations or antibiotic resistance genes that reduce antibiotic susceptibility. Finally, antibiotic resistance can be transiently acquired in the presence of specific compounds or under some growing conditions. The current article provides information on methods currently used to analyze intrinsic, mutation-driven, and transient antibiotic resistance in P. aeruginosa.

The evolution of resistance to synergistic multi-drug combinations is more complex than evolving resistance to each individual drug component

Multidrug antibiotic resistance is an urgent public health concern. Multiple strategies have been suggested to alleviate this problem, including the use of antibiotic combinations and cyclic therapies. We examine how adaptation to (1) combinations of drugs affects resistance to individual drugs, and to (2) individual drugs alters responses to drug combinations. To evaluate this, we evolved multiple strains of drug resistant Staphylococcus epidermidis in the lab. We show that evolving resistance to four highly synergistic combinations does not result in cross-resistance to all of their components. Likewise, prior resistance to one antibiotic in a combination does not guarantee survival when exposed to the combination. We also identify four 3-step and four 2-step treatments that inhibit bacterial growth and confer collateral sensitivity with each step, impeding the development of multidrug resistance. This study highlights the importance of considering higher-order drug combinations in sequential therapies and how antibiotic interactions can influence the evolutionary trajectory of bacterial populations.

Evolutionary Epidemiology Consequences of Trait-Dependent Control of Heterogeneous Parasites

AbstractDisease control can induce both demographic and evolutionary responses in host-parasite systems. Foreseeing the outcome of control therefore requires knowledge of the eco-evolutionary feedback between control and system. Previous work has assumed that control strategies have a homogeneous effect on the parasite population. However, this is not true when control targets those traits that confer to the parasite heterogeneous levels of resistance, which can additionally be related to other key parasite traits through evolutionary trade-offs. In this work, we develop a minimal model coupling epidemiological and evolutionary dynamics to explore possible trait-dependent effects of control strategies. In particular, we consider a parasite expressing continuous levels of a trait-determining resource exploitation and a control treatment that can be either positively or negatively correlated with that trait. We demonstrate the potential of trait-dependent control by considering that the decision maker may want to minimize both the damage caused by the disease and the use of treatment, due to possible environmental or economic costs. We identify efficient strategies showing that the optimal type of treatment depends on the amount applied. Our results pave the way for the study of control strategies based on evolutionary constraints, such as collateral sensitivity and resistance costs, which are receiving increasing attention for both public health and agricultural purposes.

Fitness Costs of Antibiotic Resistance Impede the Evolution of Resistance to Other Antibiotics

Antibiotic resistance is a major threat to global health, claiming the lives of millions every year. With a nearly dry antibiotic development pipeline, novel strategies are urgently needed to combat resistant pathogens. One emerging strategy is the use of sequential antibiotic therapy, postulated to reduce the rate at which antibiotic resistance evolves. Here, we use the soft agar gradient evolution (SAGE) system to carry out high-throughput in vitro bacterial evolution against antibiotic pressure. We find that evolution of resistance to the antibiotic chloramphenicol (CHL) severely affects bacterial fitness, slowing the rate at which resistance to the antibiotics nitrofurantoin and streptomycin emerges. In vitro acquisition of compensatory mutations in the CHL-resistant cells markedly improves fitness and nitrofurantoin adaptation rates but fails to restore rates to wild-type levels against streptomycin. Genome sequencing reveals distinct evolutionary paths to resistance in fitness-impaired populations, suggesting resistance trade-offs in favor of mitigation of fitness costs. We show that the speed of bacterial fronts in SAGE plates is a reliable indicator of adaptation rates and evolutionary trajectories to resistance. Identification of antibiotics whose mutational resistance mechanisms confer stable impairments may help clinicians prescribe sequential antibiotic therapies that are less prone to resistance evolution.

Translating eco-evolutionary biology into therapy to tackle antibiotic resistance

Antibiotic resistance is currently one of the most important public health problems. The golden age of antibiotic discovery ended decades ago, and new approaches are urgently needed. Therefore, preserving the efficacy of the antibiotics currently in use and developing compounds and strategies that specifically target antibiotic-resistant pathogens is critical. The identification of robust trends of antibiotic resistance evolution and of its associated trade-offs, such as collateral sensitivity or fitness costs, is invaluable for the design of rational evolution-based, ecology-based treatment approaches. In this Review, we discuss these evolutionary trade-offs and how such knowledge can aid in informing combination or alternating antibiotic therapies against bacterial infections. In addition, we discuss how targeting bacterial metabolism can enhance drug activity and impair antibiotic resistance evolution. Finally, we explore how an improved understanding of the original physiological function of antibiotic resistance determinants, which have evolved to reach clinical resistance after a process of historical contingency, may help to tackle antibiotic resistance.

Treatment of evolving cancers will require dynamic decision support

Cancer research has traditionally focused on developing new agents, but an underexplored question is that of the dose and frequency of existing drugs. Based on the modus operandi established in the early days of chemotherapies, most drugs are administered according to predetermined schedules that seek to deliver the maximum tolerated dose and are only adjusted for toxicity. However, we believe that the complex, evolving nature of cancer requires a more dynamic and personalized approach. Chronicling the milestones of the field, we show that the impact of schedule choice crucially depends on processes driving treatment response and failure. As such, cancer heterogeneity and evolution dictate that a one-size-fits-all solution is unlikely-instead, each patient should be mapped to the strategy that best matches their current disease characteristics and treatment objectives (i.e. their 'tumorscape'). To achieve this level of personalization, we need mathematical modeling. In this perspective, we propose a five-step 'Adaptive Dosing Adjusted for Personalized Tumorscapes (ADAPT)' paradigm to integrate data and understanding across scales and derive dynamic and personalized schedules. We conclude with promising examples of model-guided schedule personalization and a call to action to address key outstanding challenges surrounding data collection, model development, and integration.

Nothing makes sense in drug resistance except in the light of evolution

Our ability to fight infectious diseases is being increasingly compromised due to the emergence and spread of pathogens that become resistant to one or several drugs. This phenomenon is ubiquitous among pathogens and has parallels in cancer treatment. Given the urgency of the problem, there is a need for a paradigm shift in drug therapy toward one in which the objective to prevent the evolution of drug resistance is considered alongside the main objective of eliminating the infection or tumor. Here, I stress the importance of considering an evolutionary perspective to achieve this goal, and review recent advances in this direction, including therapies that exploit the fitness trade-offs of resistance.

Diverse evolutionary pathways challenge the use of collateral sensitivity as a strategy to suppress resistance

Drug resistance remains a major obstacle to malaria control and eradication efforts, necessitating the development of novel therapeutic strategies to treat this disease. Drug combinations based on collateral sensitivity, wherein resistance to one drug causes increased sensitivity to the partner drug, have been proposed as an evolutionary strategy to suppress the emergence of resistance in pathogen populations. In this study, we explore collateral sensitivity between compounds targeting the Plasmodium dihydroorotate dehydrogenase (DHODH). We profiled the cross-resistance and collateral sensitivity phenotypes of several DHODH mutant lines to a diverse panel of DHODH inhibitors. We focus on one compound, TCMDC-125334, which was active against all mutant lines tested, including the DHODH C276Y line, which arose in selections with the clinical candidate DSM265. In six selections with TCMDC-125334, the most common mechanism of resistance to this compound was copy number variation of the dhodh locus, although we did identify one mutation, DHODH I263S, which conferred resistance to TCMDC-125334 but not DSM265. We found that selection of the DHODH C276Y mutant with TCMDC-125334 yielded additional genetic changes in the dhodh locus. These double mutant parasites exhibited decreased sensitivity to TCMDC-125334 and were highly resistant to DSM265. Finally, we tested whether collateral sensitivity could be exploited to suppress the emergence of resistance in the context of combination treatment by exposing wildtype parasites to both DSM265 and TCMDC-125334 simultaneously. This selected for parasites with a DHODH V532A mutation which were cross-resistant to both compounds and were as fit as the wildtype parent in vitro. The emergence of these cross-resistant, evolutionarily fit parasites highlights the mutational flexibility of the DHODH enzyme.

Diverse mutant selection windows shape spatial heterogeneity in evolving populations

Mutant selection windows (MSWs), the range of drug concentrations that select for drug-resistant mutants, have long been used as a model for predicting drug resistance and designing optimal dosing strategies in infectious disease. The canonical MSW model offers comparisons between two subtypes at a time: drug-sensitive and drug-resistant. In contrast, the fitness landscape model with N alleles, which maps genotype to fitness, allows comparisons between N genotypes simultaneously, but does not encode continuous drug response data. In clinical settings, there may be a wide range of drug concentrations selecting for a variety of genotypes. Therefore, there is a need for a more robust model of the pathogen response to therapy to predict resistance and design new therapeutic approaches. Fitness seascapes, which model genotype-by-environment interactions, permit multiple MSW comparisons simultaneously by encoding genotype-specific dose-response data. By comparing dose-response curves, one can visualize the range of drug concentrations where one genotype is selected over another. In this work, we show how N-allele fitness seascapes allow for N*2N-1 unique MSW comparisons. In spatial drug diffusion models, we demonstrate how fitness seascapes reveal spatially heterogeneous MSWs, extending the MSW model to more accurately reflect the selection fo drug resistant genotypes. Furthermore, we find that the spatial structure of MSWs shapes the evolution of drug resistance in an agent-based model. Our work highlights the importance and utility of considering dose-dependent fitness seascapes in evolutionary medicine.

Evolutionary druggability: leveraging low-dimensional fitness landscapes towards new metrics for antimicrobial applications

The term "druggability" describes the molecular properties of drugs or targets in pharmacological interventions and is commonly used in work involving drug development for clinical applications. There are no current analogues for this notion that quantify the drug-target interaction with respect to a given target variant's sensitivity across a breadth of drugs in a panel, or a given drug's range of effectiveness across alleles of a target protein. Using data from low-dimensional empirical fitness landscapes composed of 16 β-lactamase alleles and seven β-lactam drugs, we introduce two metrics that capture (i) the average susceptibility of an allelic variant of a drug target to any available drug in a given panel ("variant vulnerability"), and (ii) the average applicability of a drug (or mixture) across allelic variants of a drug target ("drug applicability"). Finally, we (iii) disentangle the quality and magnitude of interactions between loci in the drug target and the seven drug environments in terms of their mutation by mutation by environment (G × G × E) interactions, offering mechanistic insight into the variant variability and drug applicability metrics. Summarizing, we propose that our framework can be applied to other datasets and pathogen-drug systems to understand which pathogen variants in a clinical setting are the most concerning (low variant vulnerability), and which drugs in a panel are most likely to be effective in an infection defined by standing genetic variation in the pathogen drug target (high drug applicability).

Evolution of antimicrobial resistance in E. coli biofilm treated with high doses of ciprofloxacin

The evolution of antimicrobial resistance (AMR) has mainly been studied in planktonic bacteria exposed to sub-inhibitory antimicrobial (AM) concentrations. However, in a number of infections that are treated with AMs the bacteria are located in biofilms where they tolerate high doses of AM. In the present study, we continuously exposed biofilm residing E. coli at body temperature to high ciprofloxacin (CIP) concentrations increasing from 4 to 130 times the minimal inhibitory concentration (MIC), i.e., from 0.06 to 2.0 mg/L. After 1 week, the biofilms were full of CIP resistant bacteria. The evolutionary trajectory observed was the same as described in the literature for planktonic bacteria, i.e., starting with a single mutation in the target gene gyrA followed by mutations in parC, gyrB, and parE, as well as in genes for regulation of multidrug efflux pump systems and outer membrane porins. Strains with higher numbers of these mutations also displayed higher MIC values. Furthermore, the evolution of CIP resistance was more rapid, and resulted in strains with higher MIC values, when the bacteria were biofilm residing than when they were in a planktonic suspension. These results may indicate that extensive clinical AM treatment of biofilm-residing bacteria may not only fail to eradicate the infection but also pose an increased risk of AMR development.

Pervasive genotype-by-environment interactions shape the fitness effects of antibiotic resistance mutations

The fitness effects of antibiotic resistance mutations are a major driver of resistance evolution. While the nutrient environment affects bacterial fitness, experimental studies of resistance typically measure fitness of mutants in a single environment only. We explored how the nutrient environment affected the fitness effects of rifampicin-resistant rpoB mutations in Escherichia coli under several conditions critical for the emergence and spread of resistance-the presence of primary or secondary antibiotic, or the absence of any antibiotic. Pervasive genotype-by-environment (GxE) interactions determined fitness in all experimental conditions, with rank order of fitness in the presence and absence of antibiotics being strongly dependent on the nutrient environment. GxE interactions also affected the magnitude and direction of collateral effects of secondary antibiotics, in some cases so drastically that a mutant that was highly sensitive in one nutrient environment exhibited cross-resistance to the same antibiotic in another. It is likely that the mutant-specific impact of rpoB mutations on the global transcriptome underpins the observed GxE interactions. The pervasive, mutant-specific GxE interactions highlight the importance of doing what is rarely done when studying the evolution and spread of resistance in experimental and clinical work: assessing fitness of antibiotic-resistant mutants across a range of relevant environments.

Graph 'texture' features as novel metrics that can summarize complex biological graphs

Objective.Image texture features, such as those derived by Haralicket al, are a powerful metric for image classification and are used across fields including cancer research. Our aim is to demonstrate how analogous texture features can be derived for graphs and networks. We also aim to illustrate how these new metrics summarize graphs, may aid comparative graph studies, may help classify biological graphs, and might assist in detecting dysregulation in cancer.Approach.We generate the first analogies of image texture for graphs and networks. Co-occurrence matrices for graphs are generated by summing over all pairs of neighboring nodes in the graph. We generate metrics for fitness landscapes, gene co-expression and regulatory networks, and protein interaction networks. To assess metric sensitivity we varied discretization parameters and noise. To examine these metrics in the cancer context we compare metrics for both simulated and publicly available experimental gene expression and build random forest classifiers for cancer cell lineage.Main results.Our novel graph 'texture' features are shown to be informative of graph structure and node label distributions. The metrics are sensitive to discretization parameters and noise in node labels. We demonstrate that graph texture features vary across different biological graph topologies and node labelings. We show how our texture metrics can be used to classify cell line expression by lineage, demonstrating classifiers with 82% and 89% accuracy.Significance.New metrics provide opportunities for better comparative analyzes and new models for classification. Our texture features are novel second-order graph features for networks or graphs with ordered node labels. In the complex cancer informatics setting, evolutionary analyses and drug response prediction are two examples where new network science approaches like this may prove fruitful.

Collateral sensitivity: An evolutionary trade-off between antibiotic resistance mechanisms, attractive for dealing with drug-resistance crisis

Background

The discovery and development of antimicrobial drugs were one of the most significant advances in medicine, but the evolution of microbial resistance limited the efficiency of these drugs.

Aim

This paper reviews the collateral sensitivity in bacteria and its potential and limitation as a new target for treating infections.

Results and Discussion

Knowledge mechanisms of resistance to antimicrobial agents are useful to trace a practical approach to treat and control of resistant pathogens. The effect of a resistance mechanism to certain antibiotics on the susceptibility or resistance to other drugs is a key point that may be helpful for applying a strategy to control resistance challenges. In an evolutionary trade-off known as collateral sensitivity, the resistance mechanism to a certain drug may be mediated by the hypersensitivity to other drugs. Collateral sensitivity has been described for different drugs in various bacteria, but the molecular mechanisms affecting susceptibility are not well demonstrated. Collateral sensitivity could be studied to detect its potential in the battle against resistance crisis as well as in the treatment of pathogens adapting to antibiotics. Collateral sensitivity-based antimicrobial therapy may have the potential to limit the emergence of antibiotic resistance.

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.

Exploiting convergent phenotypes to derive a pan-cancer cisplatin response gene expression signature

Precision medicine offers remarkable potential for the treatment of cancer, but is largely focused on tumors that harbor actionable mutations. Gene expression signatures can expand the scope of precision medicine by predicting response to traditional (cytotoxic) chemotherapy agents without relying on changes in mutational status. We present a new signature extraction method, inspired by the principle of convergent phenotypes, which states that tumors with disparate genetic backgrounds may evolve similar phenotypes independently. This evolutionary-informed method can be utilized to produce consensus signatures predictive of response to over 200 chemotherapeutic drugs found in the Genomics of Drug Sensitivity in Cancer (GDSC) Database. Here, we demonstrate its use by extracting the Cisplatin Response Signature (CisSig). We show that this signature can predict cisplatin response within carcinoma-based cell lines from the GDSC database, and expression of the signatures aligns with clinical trends seen in independent datasets of tumor samples from The Cancer Genome Atlas (TCGA) and Total Cancer Care (TCC) database. Finally, we demonstrate preliminary validation of CisSig for use in muscle-invasive bladder cancer, predicting overall survival in a small cohort of patients who undergo cisplatin-containing chemotherapy. This methodology can be used to produce robust signatures that, with further clinical validation, may be used for the prediction of traditional chemotherapeutic response, dramatically increasing the reach of personalized medicine in cancer.

Collateral sensitivity profiling in drug-resistant Escherichia coli identifies natural products suppressing cephalosporin resistance

The rapid emergence of antimicrobial resistance presents serious health challenges to the management of infectious diseases, a problem that is further exacerbated by slowing rates of antimicrobial drug discovery in recent years. The phenomenon of collateral sensitivity (CS), whereby resistance to one drug is accompanied by increased sensitivity to another, provides new opportunities to address both these challenges. Here, we present a high-throughput screening platform termed Collateral Sensitivity Profiling (CSP) to map the difference in bioactivity of large chemical libraries across 29 drug-resistant strains of E. coli. CSP screening of 80 commercial antimicrobials demonstrated multiple CS interactions. Further screening of a 6195-member natural product library revealed extensive CS relationships in nature. In particular, we report the isolation of known and new analogues of borrelidin A with potent CS activities against cephalosporin-resistant strains. Co-dosing ceftazidime with borrelidin A slows broader cephalosporin resistance with no recognizable resistance to borrelidin A itself.

Tackling antibiotic resistance by inducing transient and robust collateral sensitivity

Collateral sensitivity (CS) is an evolutionary trade-off traditionally linked to the mutational acquisition of antibiotic resistance (AR). However, AR can be temporally induced, and the possibility that this causes transient, non-inherited CS, has not been addressed. Mutational acquisition of ciprofloxacin resistance leads to robust CS to tobramycin in pre-existing antibiotic-resistant mutants of Pseudomonas aeruginosa. Further, the strength of this phenotype is higher when nfxB mutants, over-producing the efflux pump MexCD-OprJ, are selected. Here, we induce transient nfxB-mediated ciprofloxacin resistance by using the antiseptic dequalinium chloride. Notably, non-inherited induction of AR renders transient tobramycin CS in the analyzed antibiotic-resistant mutants and clinical isolates, including tobramycin-resistant isolates. Further, by combining tobramycin with dequalinium chloride we drive these strains to extinction. Our results support that transient CS could allow the design of new evolutionary strategies to tackle antibiotic-resistant infections, avoiding the acquisition of AR mutations on which inherited CS depends.

The evolution of antibiotic resistance is associated with collateral drug phenotypes in Mycobacterium tuberculosis

The increasing incidence of drug resistance in Mycobacterium tuberculosis has diminished the efficacy of almost all available antibiotics, complicating efforts to combat the spread of this global health burden. Alongside the development of new drugs, optimised drug combinations are needed to improve treatment success and prevent the further spread of antibiotic resistance. Typically, antibiotic resistance leads to reduced sensitivity, yet in some cases the evolution of drug resistance can lead to enhanced sensitivity to unrelated drugs. This phenomenon of collateral sensitivity is largely unexplored in M. tuberculosis but has the potential to identify alternative therapeutic strategies to combat drug-resistant strains that are unresponsive to current treatments. Here, by using drug susceptibility profiling, genomics and evolutionary studies we provide evidence for the existence of collateral drug sensitivities in an isogenic collection M. tuberculosis drug-resistant strains. Furthermore, in proof-of-concept studies, we demonstrate how collateral drug phenotypes can be exploited to select against and prevent the emergence of drug-resistant strains. This study highlights that the evolution of drug resistance in M. tuberculosis leads to collateral drug responses that can be exploited to design improved drug regimens.

Rapid Phenotypic Convergence towards Collateral Sensitivity in Clinical Isolates of Pseudomonas aeruginosa Presenting Different Genomic Backgrounds

Collateral sensitivity (CS) is an evolutionary trade-off by which acquisition of resistance to an antibiotic leads to increased susceptibility to another. This Achilles' heel of antibiotic resistance could be exploited to design evolution-based strategies for treating bacterial infections. To date, most studies in the field have focused on the identification of CS patterns in model strains. However, one of the main requirements for the clinical application of this trade-off is that it must be robust and has to emerge in different genomic backgrounds, including preexisting drug-resistant isolates, since infections are frequently caused by pathogens already resistant to antibiotics. Here, we report the first analysis of CS robustness in clinical strains of Pseudomonas aeruginosa presenting different ab initio mutational resistomes. We identified a robust CS pattern associated with short-term evolution in the presence of ciprofloxacin of clinical P. aeruginosa isolates, including representatives of high-risk epidemic clones belonging to sequence type (ST) 111, ST175, and ST244. We observed the acquisition of different ciprofloxacin resistance mutations in strains presenting varied STs and different preexisting mutational resistomes. Importantly, despite these genetic differences, the use of ciprofloxacin led to a robust CS to aztreonam and tobramycin. In addition, we describe the possible application of this evolutionary trade-off to drive P. aeruginosa infections to extinction by using the combination of ciprofloxacin-tobramycin or ciprofloxacin-aztreonam. Our results support the notion that the identification of robust patterns of CS may establish the basis for developing evolution-informed treatment strategies to tackle bacterial infections, including those due to antibiotic-resistant pathogens. IMPORTANCE Collateral sensitivity (CS) is a trade-off of antibiotic resistance evolution that could be exploited to design strategies for treating bacterial infections. Clinical application of CS requires it to robustly emerge in different genomic backgrounds. In this study, we performed an analysis to identify robust patterns of CS associated with the use of ciprofloxacin in clinical isolates of P. aeruginosa presenting different mutational resistomes and including high-risk epidemic clones (ST111, ST175, and ST244). We demonstrate the robustness of CS to tobramycin and aztreonam and the potential application of this evolutionary observation to drive P. aeruginosa infections to extinction. Our results support the notion that the identification of robust CS patterns may establish the basis for developing evolutionary strategies to tackle bacterial infections, including those due to antibiotic-resistant pathogens.

Resistance-resistant antibacterial treatment strategies

Antibiotic resistance is a major danger to public health that threatens to claim the lives of millions of people per year within the next few decades. Years of necessary administration and excessive application of antibiotics have selected for strains that are resistant to many of our currently available treatments. Due to the high costs and difficulty of developing new antibiotics, the emergence of resistant bacteria is outpacing the introduction of new drugs to fight them. To overcome this problem, many researchers are focusing on developing antibacterial therapeutic strategies that are "resistance-resistant"-regimens that slow or stall resistance development in the targeted pathogens. In this mini review, we outline major examples of novel resistance-resistant therapeutic strategies. We discuss the use of compounds that reduce mutagenesis and thereby decrease the likelihood of resistance emergence. Then, we examine the effectiveness of antibiotic cycling and evolutionary steering, in which a bacterial population is forced by one antibiotic toward susceptibility to another antibiotic. We also consider combination therapies that aim to sabotage defensive mechanisms and eliminate potentially resistant pathogens by combining two antibiotics or combining an antibiotic with other therapeutics, such as antibodies or phages. Finally, we highlight promising future directions in this field, including the potential of applying machine learning and personalized medicine to fight antibiotic resistance emergence and out-maneuver adaptive pathogens.

Antibiotic Cycling Affects Resistance Evolution Independently of Collateral Sensitivity

Antibiotic cycling has been proposed as a promising approach to slow down resistance evolution against currently employed antibiotics. It remains unclear, however, to which extent the decreased resistance evolution is the result of collateral sensitivity, an evolutionary trade-off where resistance to one antibiotic enhances the sensitivity to the second, or due to additional effects of the evolved genetic background, in which mutations accumulated during treatment with a first antibiotic alter the emergence and spread of resistance against a second antibiotic via other mechanisms. Also, the influence of antibiotic exposure patterns on the outcome of drug cycling is unknown. Here, we systematically assessed the effects of the evolved genetic background by focusing on the first switch between two antibiotics against Salmonella Typhimurium, with cefotaxime fixed as the first and a broad variety of other drugs as the second antibiotic. By normalizing the antibiotic concentrations to eliminate the effects of collateral sensitivity, we demonstrated a clear contribution of the evolved genetic background beyond collateral sensitivity, which either enhanced or reduced the adaptive potential depending on the specific drug combination. We further demonstrated that the gradient strength with which cefotaxime was applied affected both cefotaxime resistance evolution and adaptation to second antibiotics, an effect that was associated with higher levels of clonal interference and reduced cost of resistance in populations evolved under weaker cefotaxime gradients. Overall, our work highlights that drug cycling can affect resistance evolution independently of collateral sensitivity, in a manner that is contingent on the antibiotic exposure pattern.

Antimicrobial Peptide Combination Can Hinder Resistance Evolution

Antibiotic-resistant microbial pathogens are becoming a major threat to human health. Therefore, there is an urgent need to develop new alternatives to conventional antibiotics. One such promising alternative is antimicrobial peptides (AMPs), which are produced by virtually all organisms and typically inhibit bacteria via membrane disruption. However, previous studies demonstrated that bacteria can rapidly develop AMP resistance. Here, we study whether combination therapy, known to be able to inhibit the evolution of resistance to conventional antibiotics, can also hinder the evolution of AMP resistance. To do so, we evolved the opportunistic pathogen Staphylococcus aureus in the presence of individual AMPs, AMP pairs, and a combinatorial antimicrobial peptide library. Treatment with some AMP pairs indeed hindered the evolution of resistance compared with individual AMPs. In particular, resistance to pairs was delayed when resistance to the individual AMPs came at a cost of impaired bacterial growth and did not confer cross-resistance to other tested AMPs. The lowest level of resistance evolved during treatment with the combinatorial antimicrobial peptide library termed random antimicrobial peptide mixture, which contains more than a million different peptides. A better understanding of how AMP combinations affect the evolution of resistance is a crucial step in order to design "resistant proof" AMP cocktails that will offer a sustainable treatment option for antibiotic-resistant pathogens. IMPORTANCE The main insights gleaned from this study are the following. (i) AMP combination treatment can delay the evolution of resistance in S. aureus. Treatment with some AMP pairs resulted in significantly lower resistance then treatment with either of the individual AMPs. Treatment with a random AMP library resulted in no detectable resistance. (ii) The rate at which resistance to combination arises correlates with the cost of resistance to individual AMPs and their cross-resistance. In particular, combinations to which the least resistance arose involved AMPs with high fitness cost of resistance and low cross-resistance. (iii) No broad-range AMP resistance evolved. Strains that evolved resistance to some AMPs typically remained sensitive to other AMPs, alleviating concerns regarding the evolution of resistance to immune system AMPs in response to AMP treatment.

Challenges and potential solutions for studying the genetic and phenotypic architecture of adaptation in microbes

All organisms are defined by the makeup of their DNA. Over billions of years, the structure and information contained in that DNA, often referred to as genetic architecture, have been honed by a multitude of evolutionary processes. Mutations that cause genetic elements to change in a way that results in beneficial phenotypic change are more likely to survive and propagate through the population in a process known as adaptation. Recent work reveals that the genetic targets of adaptation are varied and can change with genetic background. Further, seemingly similar adaptive mutations, even within the same gene, can have diverse and unpredictable effects on phenotype. These challenges represent major obstacles in predicting adaptation and evolution. In this review, we cover these concepts in detail and identify three emerging synergistic solutions: higher-throughput evolution experiments combined with updated genotype-phenotype mapping strategies and physiological models. Our review largely focuses on recent literature in yeast, and the field seems to be on the cusp of a new era with regard to studying the predictability of evolution.

Adaptive laboratory evolution and independent component analysis disentangle complex vancomycin adaptation trajectories

Infections with methicillin-resistant Staphylococcus aureus (MRSA) are associated with significant morbidity and mortality. Vancomycin is a last-line antibiotic used to treat MRSA infections; however, strains with decreased susceptibility to vancomycin (vancomycin-intermediate S. aureus [VISA]) have been spreading, and VISA infections are associated with prolonged therapeutic treatment and treatment failure. To map out the evolutionary trajectory behind VISA development, we characterized the mutational, transcriptional, and phenotypic landscape of 10 lineages of S. aureus USA300 strain JE2 that evolved in parallel to vancomycin. We demonstrate that MRSA strains adapt to vancomycin by divergent pathways leading to high or low oxacillin susceptibility characterized by mutational or transcriptional profiles. Our results point to diagnostic possibilities that may support personalized antibiotic treatment regimes.

Human infections with methicillin-resistant Staphylococcus aureus (MRSA) are commonly treated with vancomycin, and strains with decreased susceptibility, designated as vancomycin-intermediate S. aureus (VISA), are associated with treatment failure. Here, we profiled the phenotypic, mutational, and transcriptional landscape of 10 VISA strains adapted by laboratory evolution from one common MRSA ancestor, the USA300 strain JE2. Using functional and independent component analysis, we found that: 1) despite the common genetic background and environmental conditions, the mutational landscape diverged between evolved strains and included mutations previously associated with vancomycin resistance (in vraT, graS, vraFG, walKR, and rpoBCD) as well as novel adaptive mutations (SAUSA300_RS04225, ssaA, pitAR, and sagB); 2) the first wave of mutations affected transcriptional regulators and the second affected genes involved in membrane biosynthesis; 3) expression profiles were predominantly strain-specific except for sceD and lukG, which were the only two genes significantly differentially expressed in all clones; 4) three independent virulence systems (φSa3, SaeR, and T7SS) featured as the most transcriptionally perturbed gene sets across clones; 5) there was a striking variation in oxacillin susceptibility across the evolved lineages (from a 10-fold increase to a 63-fold decrease) that also arose in clinical MRSA isolates exposed to vancomycin and correlated with susceptibility to teichoic acid inhibitors; and 6) constitutive expression of the VraR regulon explained cross-susceptibility, while mutations in walK were associated with cross-resistance. Our results show that adaptation to vancomycin involves a surprising breadth of mutational and transcriptional pathways that affect antibiotic susceptibility and possibly the clinical outcome of infections.

Evolution of Habitat-Dependent Antibiotic Resistance in Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic human pathogen that usually causes difficult-to-treat infections due to its low intrinsic antibiotic susceptibility and outstanding capacity for becoming resistant to antibiotics. In addition, it has a remarkable metabolic versatility, being able to grow in different habitats, from natural niches to different and changing inpatient environments. Study of the environmental conditions that shape genetic and phenotypic changes of P. aeruginosa toward antibiotic resistance supposes a novelty, since experimental evolution assays are usually performed with well-defined antibiotics in regular laboratory growth media. Therefore, in this work we address the extent to which the nutrients’ availability may constrain the evolution of antibiotic resistance. We determined that P. aeruginosa genetic trajectories toward resistance to tobramycin, ceftazidime, and ceftazidime-avibactam are different when evolving in laboratory rich medium, urine, or synthetic sputum. Furthermore, our study, linking genotype with phenotype, showed a clear impact of each analyzed environment on both the fitness and resistance level associated with particular resistance mutations. This indicates that the phenotype associated with specific resistance mutations is variable and dependent on the bacterial metabolic state in each particular habitat. Our results support that the design of evolution-based strategies to tackle P. aeruginosa infections should be based on robust patterns of evolution identified within each particular infection and body location.

IMPORTANCE Predicting evolution toward antibiotic resistance (AR) and its associated trade-offs, such as collateral sensitivity, is important to design evolution-based strategies to tackle AR. However, the effect of nutrients' availability on such evolution, particularly those that can be found under in vivo infection conditions, has been barely addressed. We analyzed the evolutionary patterns of P. aeruginosa in the presence of antibiotics in different media, including urine and synthetic sputum, whose compositions are similar to the ones in infections, finding that AR evolution differs, depending on growth conditions. Furthermore, the representative mutants isolated under each condition tested render different AR levels and fitness costs, depending on nutrients’ availability, supporting the idea that environmental constraints shape the phenotypes associated with specific AR mutations. Consequently, the selection of AR mutations that render similar phenotypes is environment dependent. The analysis of evolution patterns toward AR requires studying growth conditions mimicking those that bacteria face during in vivo evolution.

Prediction of Antibiotic Resistance Evolution by Growth Measurement of All Proximal Mutants of Beta-Lactamase

The antibiotic resistance crisis continues to threaten human health. Better predictions of the evolution of antibiotic resistance genes could contribute to the design of more sustainable treatment strategies. However, comprehensive prediction of antibiotic resistance gene evolution via laboratory approaches remains challenging. By combining site-specific integration and high-throughput sequencing, we quantified relative growth under the respective selection of cefotaxime or ceftazidime selection in ∼23,000 Escherichia coli MG1655 strains that each carried a unique, single-copy variant of the extended-spectrum β-lactamase gene blaCTX-M-14 at the chromosomal att HK022 site. Significant synergistic pleiotropy was observed within four subgenic regions, suggesting key regions for the evolution of resistance to both antibiotics. Moreover, we propose PEARP and PEARR, two deep-learning models with strong clinical correlations, for the prospective and retrospective prediction of blaCTX-M-14 evolution, respectively. Single to quintuple mutations of blaCTX-M-14 predicted to confer resistance by PEARP were significantly enriched among the clinical isolates harboring blaCTX-M-14 variants, and the PEARR scores matched the minimal inhibitory concentrations obtained for the 31 intermediates in all hypothetical trajectories. Altogether, we conclude that the measurement of local fitness landscape enables prediction of the evolutionary trajectories of antibiotic resistance genes, which could be useful for a broad range of clinical applications, from resistance prediction to designing novel treatment strategies.

Cancer's second genome: Microbial cancer diagnostics and redefining clonal evolution as a multispecies process: Humans and their tumors are not aseptic, and the multispecies nature of cancer modulates clinical care and clonal evolution: Humans and their tumors are not aseptic, and the multispecies nature of cancer modulates clinical care and clonal evolution

The presence and role of microbes in human cancers has come full circle in the last century. Tumors are no longer considered aseptic, but implications for cancer biology and oncology remain underappreciated. Opportunities to identify and build translational diagnostics, prognostics, and therapeutics that exploit cancer's second genome-the metagenome-are manifold, but require careful consideration of microbial experimental idiosyncrasies that are distinct from host-centric methods. Furthermore, the discoveries of intracellular and intra-metastatic cancer bacteria necessitate fundamental changes in describing clonal evolution and selection, reflecting bidirectional interactions with non-human residents. Reconsidering cancer clonality as a multispecies process similarly holds key implications for understanding metastasis and prognosing therapeutic resistance while providing rational guidance for the next generation of bacterial cancer therapies. Guided by these new findings and challenges, this Review describes opportunities to exploit cancer's metagenome in oncology and proposes an evolutionary framework as a first step towards modeling multispecies cancer clonality. Also see the video abstract here: https://youtu.be/-WDtIRJYZSs.

Collateral sensitivity to pleuromutilins in vancomycin-resistant Enterococcus faecium

The acquisition of resistance to one antibiotic sometimes leads to collateral sensitivity to a second antibiotic. Here, we show that vancomycin resistance in Enterococcus faecium is associated with a remarkable increase in susceptibility to pleuromutilin antibiotics (such as lefamulin), which target the bacterial ribosome. The trade-off between vancomycin and pleuromutilins is mediated by epistasis between the van gene cluster and msrC, encoding an ABC-F protein that protects bacterial ribosomes from antibiotic targeting. In mouse models of vancomycin-resistant E. faecium colonization and septicemia, pleuromutilin treatment reduces colonization and improves survival more effectively than standard therapy (linezolid). Our findings suggest that pleuromutilins may be useful for the treatment of vancomycin-resistant E. faecium infections.

Genomic heterogeneity underlies multidrug resistance in Pseudomonas aeruginosa: A population-level analysis beyond susceptibility testing

BACKGROUND

Pseudomonas aeruginosa is a persistent and difficult-to-treat pathogen in many patients, especially those with Cystic Fibrosis (CF). Herein, we describe a longitudinal analysis of a series of multidrug resistant (MDR) P. aeruginosa isolates recovered in a 17-month period, from a young female CF patient who underwent double lung transplantation. Our goal was to understand the genetic basis of the observed resistance phenotypes, establish the genomic population diversity, and define the nature of sequence evolution over time.

METHODS

Twenty-two sequential P. aeruginosa isolates were obtained within a 17-month period, before and after a double-lung transplant. At the end of the study period, antimicrobial susceptibility testing, whole genome sequencing (WGS), phylogenetic analyses and RNAseq were performed in order to understand the genetic basis of the observed resistance phenotypes, establish the genomic population diversity, and define the nature of sequence changes over time.

RESULTS

The majority of isolates were resistant to almost all tested antibiotics. A phylogenetic reconstruction revealed 3 major clades representing a genotypically and phenotypically heterogeneous population. The pattern of mutation accumulation and variation of gene expression suggested that a group of closely related strains was present in the patient prior to transplantation and continued to change throughout the course of treatment. A trend toward accumulation of mutations over time was observed. Different mutations in the DNA mismatch repair gene mutL consistent with a hypermutator phenotype were observed in two clades. RNAseq performed on 12 representative isolates revealed substantial differences in the expression of genes associated with antibiotic resistance and virulence traits.

CONCLUSIONS

The overwhelming current practice in the clinical laboratories setting relies on obtaining a pure culture and reporting the antibiogram from a few isolated colonies to inform therapy decisions. Our analyses revealed significant underlying genomic heterogeneity and unpredictable evolutionary patterns that were independent of prior antibiotic treatment, highlighting the need for comprehensive sampling and population-level analysis when gathering microbiological data in the context of CF P. aeruginosa chronic infection. Our findings challenge the applicability of antimicrobial stewardship programs based on single-isolate resistance profiles for the selection of antibiotic regimens in chronic infections such as CF.

Collateral responses to classical cytotoxic chemotherapies are heterogeneous and sensitivities are sparse

Chemotherapy resistance is a major obstacle to curing cancer patients. Combination drug regimens have shown promise as a method to overcome resistance; however, to date only some cancers have been cured with this method. Collateral sensitivity-the phenomenon whereby resistance to one drug is co-occurrent with sensitivity to a second drug-has been gaining traction as a promising new concept to guide rational design of combination regimens. Here we evolved over 100 subclones of the Eµ-Myc; p19^ARF-/- cell line to be resistant to one of four classical chemotherapy agents: doxorubicin, vincristine, paclitaxel, and cisplatin. We then surveyed collateral responses to acquisition of resistance to these agents. Although numerous collateral sensitivities have been documented for antibiotics and targeted cancer therapies, we observed only one collateral sensitivity: half of cell lines that acquired resistance to paclitaxel also acquired a collateral sensitivity to verapamil. However, we found that the mechanism of this collateral sensitivity was unrelated to the mechanism of paclitaxel resistance. Interestingly, we observed heterogeneity in the phenotypic response to acquisition of resistance to most of the drugs we tested, most notably for paclitaxel, suggesting the existence of multiple different states of resistance. Surprisingly, this phenotypic heterogeneity in paclitaxel resistant cell lines was unrelated to transcriptomic heterogeneity among those cell lines. These features of phenotypic and transcriptomic heterogeneity must be taken into account in future studies of treated tumor subclones and in design of chemotherapy combinations.

Rapid decline of ceftazidime resistance in antibiotic-free and sub-lethal environments is contingent on genetic background

Trade-offs of antibiotic resistance evolution, such as fitness cost and collateral sensitivity (CS), could be exploited to drive evolution toward antibiotic susceptibility. Decline of resistance may occur when resistance to other drug leads to CS to the first one and when compensatory mutations, or genetic reversion of the original ones, reduce fitness cost. Here we describe the impact of antibiotic-free and sublethal environments on declining ceftazidime resistance in different Pseudomonas aeruginosa resistant mutants. We determined that decline of ceftazidime resistance occurs within 450 generations, which is caused by newly acquired mutations and not by reversion of the original ones, and that the original CS of these mutants is preserved. In addition, we observed that the frequency and degree of this decline is contingent on genetic background. Our results are relevant to implement evolution-based therapeutic approaches, as well as to redefine global policies of antibiotic use, such as drug cycling.

Quantification of long-term doxorubicin response dynamics in breast cancer cell lines to direct treatment schedules

While acquired chemoresistance is recognized as a key challenge to treating many types of cancer, the dynamics with which drug sensitivity changes after exposure are poorly characterized. Most chemotherapeutic regimens call for repeated dosing at regular intervals, and if drug sensitivity changes on a similar time scale then the treatment interval could be optimized to improve treatment performance. Theoretical work suggests that such optimal schedules exist, but experimental confirmation has been obstructed by the difficulty of deconvolving the simultaneous processes of death, adaptation, and regrowth taking place in cancer cell populations. Here we present a method of optimizing drug schedules in vitro through iterative application of experimentally calibrated models, and demonstrate its ability to characterize dynamic changes in sensitivity to the chemotherapeutic doxorubicin in three breast cancer cell lines subjected to treatment schedules varying in concentration, interval between pulse treatments, and number of sequential pulse treatments. Cell populations are monitored longitudinally through automated imaging for 600–800 hours, and this data is used to calibrate a family of cancer growth models, each consisting of a system of ordinary differential equations, derived from the bi-exponential model which characterizes resistant and sensitive subpopulations. We identify a model incorporating both a period of growth arrest in surviving cells and a delay in the death of chemosensitive cells which outperforms the original bi-exponential growth model in Akaike Information Criterion based model selection, and use the calibrated model to quantify the performance of each drug schedule. We find that the inter-treatment interval is a key variable in determining the performance of sequential dosing schedules and identify an optimal retreatment time for each cell line which extends regrowth time by 40%-239%, demonstrating that the time scale of changes in chemosensitivity following doxorubicin exposure allows optimization of drug scheduling by varying this inter-treatment interval.

Acquired chemoresistance is a common cause of treatment failure in cancer. The scheduling of a multi-dose course of chemotherapeutic treatment may influence the dynamics of acquired chemoresistance, and drug schedule optimization may increase the duration of effectiveness of a particular chemotherapeutic agent for a particular patient. Here we present a method for experimentally optimizing an in vitro drug schedule through iterative rounds of experimentation and computational analysis, and demonstrate the method’s ability to improve the performance of doxorubicin treatment in three breast carcinoma cell lines. Specifically, we find that the interval between drug exposures can be optimized while holding drug concentration and number of treatments constant, suggesting that this may be a key variable to explore in future drug schedule optimization efforts. We further use this method’s model calibration and selection process to extract information about the underlying biology of the doxorubicin response, and find that the incorporation of delays on both cell death and regrowth are necessary for accurate parameterization of cell growth data.

Stability of scRNA-Seq Analysis Workflows is Susceptible to Preprocessing and is Mitigated by Regularized or Supervised Approaches

Background:

Statistical methods developed to address various questions in single-cell datasets show increased variability to different parameter regimes. In order to delineate further the robustness of commonly utilized methods for single-cell RNA-Seq, we aimed to comprehensively review scRNA-Seq analysis workflows in the setting of dimension reduction, clustering, and trajectory inference.

Methods:

We utilized datasets with temporal single-cell transcriptomics profiles from public repositories. Combining multiple methods at each level of the workflow, we have performed over 6 k analysis and evaluated the results of clustering and pseudotime estimation using adjusted rand index and rank correlation metrics. We have further integrated neural network methods to assess whether models with increased complexity can show increased bias/variance trade-off.

Results:

Combinatorial workflows showed that utilizing non-linear dimension reduction techniques such as t-SNE and UMAP are sensitive to initial preprocessing steps hence clustering results on dimension reduced space of single-cell datasets should be utilized carefully. Similarly, pseudotime estimation methods that depend on previous non-linear dimension reduction steps can result in highly variable trajectories. In contrast, methods that avoid non-linearity such as WOT can result in repeatable inferences of temporal gene expression dynamics. Furthermore, imputation methods do not improve clustering or trajectory inference results substantially in terms of repeatability. In contrast, the selection of the normalization method shows an increased effect on downstream analysis where ScTransform reduces variability overall.

Revisiting Antibiotic Resistance: Mechanistic Foundations to Evolutionary Outlook

Antibiotics are the pivotal pillar of contemporary healthcare and have contributed towards its advancement over the decades. Antibiotic resistance emerged as a critical warning to public wellbeing because of unsuccessful management efforts. Resistance is a natural adaptive tool that offers selection pressure to bacteria, and hence cannot be stopped entirely but rather be slowed down. Antibiotic resistance mutations mostly diminish bacterial reproductive fitness in an environment without antibiotics; however, a fraction of resistant populations 'accidentally' emerge as the fittest and thrive in a specific environmental condition, thus favouring the origin of a successful resistant clone. Therefore, despite the time-to-time amendment of treatment regimens, antibiotic resistance has evolved relentlessly. According to the World Health Organization (WHO), we are rapidly approaching a 'post-antibiotic' era. The knowledge gap about antibiotic resistance and room for progress is evident and unified combating strategies to mitigate the inadvertent trends of resistance seem to be lacking. Hence, a comprehensive understanding of the genetic and evolutionary foundations of antibiotic resistance will be efficacious to implement policies to force-stop the emergence of resistant bacteria and treat already emerged ones. Prediction of possible evolutionary lineages of resistant bacteria could offer an unswerving impact in precision medicine. In this review, we will discuss the key molecular mechanisms of resistance development in clinical settings and their spontaneous evolution.

Collateral Sensitivity Interactions between Antibiotics Depend on Local Abiotic Conditions

Mutations conferring resistance to one antibiotic can increase (cross-resistance) or decrease (collateral sensitivity) resistance to others. Antibiotic combinations displaying collateral sensitivity could be used in treatments that slow resistance evolution. However, lab-to-clinic translation requires understanding whether collateral effects are robust across different environmental conditions. Here, we isolated and characterized resistant mutants of Escherichia coli using five antibiotics, before measuring collateral effects on resistance to other paired antibiotics. During both isolation and phenotyping, we varied conditions in ways relevant in nature (pH, temperature, and bile). This revealed that local abiotic conditions modified expression of resistance against both the antibiotic used during isolation and other antibiotics. Consequently, local conditions influenced collateral sensitivity in two ways: by favoring different sets of mutants (with different collateral sensitivities) and by modifying expression of collateral effects for individual mutants. These results place collateral sensitivity in the context of environmental variation, with important implications for translation to real-world applications. IMPORTANCE When bacteria become resistant to an antibiotic, the genetic changes involved sometimes increase (cross-resistance) or decrease (collateral sensitivity) their resistance to other antibiotics. Antibiotic combinations showing repeatable collateral sensitivity could be used in treatment to slow resistance evolution. However, collateral sensitivity interactions may depend on the local environmental conditions that bacteria experience, potentially reducing repeatability and clinical application. Here, we show that variation in local conditions (pH, temperature, and bile salts) can influence collateral sensitivity in two ways: by favoring different sets of mutants during bacterial resistance evolution (with different collateral sensitivities to other antibiotics) and by modifying expression of collateral effects for individual mutants. This suggests that translation from the lab to the clinic of new approaches exploiting collateral sensitivity will be influenced by local abiotic conditions.

The population genetics of collateral resistance and sensitivity

Resistance mutations against one drug can elicit collateral sensitivity against other drugs. Multi-drug treatments exploiting such trade-offs can help slow down the evolution of resistance. However, if mutations with diverse collateral effects are available, a treated population may evolve either collateral sensitivity or collateral resistance. How to design treatments robust to such uncertainty is unclear. We show that many resistance mutations in Escherichia coli against various antibiotics indeed have diverse collateral effects. We propose to characterize such diversity with a joint distribution of fitness effects (JDFE) and develop a theory for describing and predicting collateral evolution based on simple statistics of the JDFE. We show how to robustly rank drug pairs to minimize the risk of collateral resistance and how to estimate JDFEs. In addition to practical applications, these results have implications for our understanding of evolution in variable environments.

Drugs known as antibiotics are the main treatment for most serious infections caused by bacteria. However, many bacteria are acquiring genetic mutations that make them resistant to the effects of one or more types of antibiotics, making them harder to eliminate.

One way to tackle drug-resistant bacteria is to develop new types of antibiotics; however, in recent years, the rate at which new antibiotics have become available has been dwindling. Using two or more existing drugs, one after another, can also be an effective way to eliminate resistant bacteria. The success of any such ‘multi-drug’ treatment lies in being able to predict whether mutations that make the bacteria resistant to one drug simultaneously make it sensitive to another, a phenomenon known as collateral sensitivity.

Different resistance mutations may have different collateral effects: some may increase the bacteria’s sensitivity to the second drug, while others might make the bacteria more resistant. However, it is currently unclear how to design robust multi-drug treatments that take this diversity of collateral effects into account.

Here, Ardell and Kryazhimskiy used a concept called JDFE (short for the joint distribution of fitness effects) to describe the diversity of collateral effects in a population of bacteria exposed to a single drug. This information was then used to mathematically model how collateral effects evolved in the population over time.

Ardell and Kryazhimskiy showed that this approach can predict how likely a population is to become collaterally sensitive or collaterally resistant to a second antibiotic. Drug pairs can then be ranked according to the risk of collateral resistance emerging, so long as information on the variety of resistance mutations available to the bacteria are included in the model.

Each year, more than 700,000 people die from infections caused by bacteria that are resistant to one or more antibiotics. The findings of Ardell and Kryazhimskiy may eventually help clinicians design multi-drug treatments that effectively eliminate bacterial infections and help to prevent more bacteria from evolving resistance to antibiotics. However, to achieve this goal, more research is needed to fully understand the range collateral effects caused by resistance mutations.

Hidden role of mutations in the evolutionary process

Mutations not only alter allele frequencies in a genetic pool but may also determine the fate of an evolutionary process. Here we study which allele fixes in a one-step, one-way model including the wild type and two adaptive mutations. We study the effect of the four basic evolutionary mechanisms-genetic drift, natural selection, mutation, and gene flow-on mutant fixation and its kinetics. Determining which allele is more likely to fix is not simply a question of comparing fitnesses and mutation rates. For instance, if the allele of interest is less fit than the other, then not only must it have a greater mutation rate, but also its mutation rate must exceed a specific threshold for it to prevail. We find exact expressions for such conditions. Our conclusions are based on the mathematical description of two extreme but important regimes, as well as on simulations.

Allogenous Selection of Mutational Collateral Resistance: Old Drugs Select for New Resistance Within Antibiotic Families

Allogeneous selection occurs when an antibiotic selects for resistance to more advanced members of the same family. The mechanisms of allogenous selection are (a) collateral expansion, when the antibiotic expands the gene and gene-containing bacterial populations favoring the emergence of other mutations, inactivating the more advanced antibiotics; (b) collateral selection, when the old antibiotic selects its own resistance but also resistance to more modern drugs; (c) collateral hyper-resistance, when resistance to the old antibiotic selects in higher degree for populations resistant to other antibiotics of the family than to itself; and (d) collateral evolution, when the simultaneous or sequential use of antibiotics of the same family selects for new mutational combinations with novel phenotypes in this family, generally with higher activity (higher inactivation of the antibiotic substrates) or broader spectrum (more antibiotics of the family are inactivated). Note that in some cases, collateral selection derives from collateral evolution. In this article, examples of allogenous selection are provided for the major families of antibiotics. Improvements in minimal inhibitory concentrations with the newest drugs do not necessarily exclude “old” antibiotics of the same family of retaining some selective power for resistance to the newest agents. If this were true, the use of older members of the same drug family would facilitate the emergence of mutational resistance to the younger drugs of the family, which is frequently based on previously established resistance traits. The extensive use of old drugs (particularly in low-income countries and in farming) might be significant for the emergence and selection of resistance to the novel members of the family, becoming a growing source of variation and selection of resistance to the whole family. In terms of future research, it could be advisable to focus antimicrobial drug discovery more on the identification of new targets and new (unique) classes of antimicrobial agents, than on the perpetual chemical exploitation of classic existing ones.