Epistasis and the Evolution of Antimicrobial Resistance

Unpredictability of the Fitness Effects of Antimicrobial Resistance Mutations Across Environments in Escherichia coli

The evolution of antimicrobial resistance (AMR) in bacteria is a major public health concern, and antibiotic restriction is often implemented to reduce the spread of resistance. These measures rely on the existence of deleterious fitness effects (i.e. costs) imposed by AMR mutations during growth in the absence of antibiotics. According to this assumption, resistant strains will be outcompeted by susceptible strains that do not pay the cost during the period of restriction. The fitness effects of AMR mutations are generally studied in laboratory reference strains grown in standard growth environments; however, the genetic and environmental context can influence the magnitude and direction of a mutation's fitness effects. In this study, we measure how three sources of variation impact the fitness effects of Escherichia coli AMR mutations: the type of resistance mutation, the genetic background of the host, and the growth environment. We demonstrate that while AMR mutations are generally costly in antibiotic-free environments, their fitness effects vary widely and depend on complex interactions between the mutation, genetic background, and environment. We test the ability of the Rough Mount Fuji fitness landscape model to reproduce the empirical data in simulation. We identify model parameters that reasonably capture the variation in fitness effects due to genetic variation. However, the model fails to accommodate the observed variation when considering multiple growth environments. Overall, this study reveals a wealth of variation in the fitness effects of resistance mutations owing to genetic background and environmental conditions, which will ultimately impact their persistence in natural populations.

Genetic factors associated with acquired phenotypic drug resistance and its compensatory evolution during tuberculosis treatment

OBJECTIVES

We elucidated the factors, evolution, and compensation of antimicrobial resistance (AMR) in Mycobacterium tuberculosis (MTB) isolates under dual pressure from the intra-host environment and anti-tuberculosis (anti-TB) drugs.

METHODS

This retrospective case-control study included 337 patients with pulmonary tuberculosis from 15 clinics in Tianjin, China, with phenotypic drug susceptibility testing results available for at least two time points between January 1, 2009 and December 31, 2016. Patients in the case group exhibited acquired AMR to isoniazid (INH) or rifampicin (RIF), while those in the control group lacked acquired AMR. The whole-genome sequencing (WGS) was conducted on 149 serial longitudinal MTB isolates from 46 patients who acquired or reversed phenotypic INH/RIF-resistance during treatment. The genetic basis, associated factors, and intra-host evolution of acquired phenotypic INH/RIF-resistance were elucidated using a combined analysis.

RESULTS

Anti-TB interruption duration of ≥30 days showed association with acquired phenotypic INH/RIF resistance (aOR = 2·2, 95% CI, 1·0-5·1) and new rpoB mutations (p = 0·024). The MTB evolution was 1·2 (95% CI, 1·02-1·38) single nucleotide polymorphisms per genome per year under dual pressure from the intra-host environment and anti-TB drugs. AMR-associated mutations occurred before phenotypic AMR appearance in cases with acquired phenotypic INH (10 of 16) and RIF (9 of 22) resistances.

DISCUSSION

Compensatory evolution may promote the fixation of INH/RIF-resistance mutations and affect phenotypic AMR. The TB treatment should be adjusted based on gene sequencing results, especially in persistent culture positivity during treatment, which highlights the clinical importance of WGS in identifying reinfection and AMR acquisition before phenotypic drug susceptibility testing.

Global epistasis in plasmid-mediated antimicrobial resistance

Antimicrobial resistance (AMR) in bacteria is a major public health threat and conjugative plasmids play a key role in the dissemination of AMR genes among bacterial pathogens. Interestingly, the association between AMR plasmids and pathogens is not random and certain associations spread successfully at a global scale. The burst of genome sequencing has increased the resolution of epidemiological programs, broadening our understanding of plasmid distribution in bacterial populations. Despite the immense value of these studies, our ability to predict future plasmid-bacteria associations remains limited. Numerous empirical studies have recently reported systematic patterns in genetic interactions that enable predictability, in a phenomenon known as global epistasis. In this perspective, we argue that global epistasis patterns hold the potential to predict interactions between plasmids and bacterial genomes, thereby facilitating the prediction of future successful associations. To assess the validity of this idea, we use previously published data to identify global epistasis patterns in clinically relevant plasmid-bacteria associations. Furthermore, using simple mechanistic models of antibiotic resistance, we illustrate how global epistasis patterns may allow us to generate new hypotheses on the mechanisms associated with successful plasmid-bacteria associations. Collectively, we aim at illustrating the relevance of exploring global epistasis in the context of plasmid biology.

Environmental modulation of global epistasis in a drug resistance fitness landscape

Interactions between mutations (epistasis) can add substantial complexity to genotype-phenotype maps, hampering our ability to predict evolution. Yet, recent studies have shown that the fitness effect of a mutation can often be predicted from the fitness of its genetic background using simple, linear relationships. This phenomenon, termed global epistasis, has been leveraged to reconstruct fitness landscapes and infer adaptive trajectories in a wide variety of contexts. However, little attention has been paid to how patterns of global epistasis may be affected by environmental variation, despite this variation frequently being a major driver of evolution. This is particularly relevant for the evolution of drug resistance, where antimicrobial drugs may change the environment faced by pathogens and shape their adaptive trajectories in ways that can be difficult to predict. By analyzing a fitness landscape of four mutations in a gene encoding an essential enzyme of P. falciparum (a parasite cause of malaria), here we show that patterns of global epistasis can be strongly modulated by the concentration of a drug in the environment. Expanding on previous theoretical results, we demonstrate that this modulation can be quantitatively explained by how specific gene-by-gene interactions are modified by drug dose. Importantly, our results highlight the need to incorporate potential environmental variation into the global epistasis framework in order to predict adaptation in dynamic environments.

Epigenetic effects of ceftriaxone-resistant Neisseria gonorrhoeae FC428 mosaic-like sequences found in PenA sequences unique to Neisseria subflava and related species

OBJECTIVES

The aim of this study was to explore the origin of the PenA mosaic amino acid sequence in the ceftriaxone-resistant Neisseria gonorrhoeae FC428 clone.

METHODS

The penA sequences of 27 Neisseria subflava pharyngeal isolates were determined by the Sanger method and penA sequences of 52 isolates from nine Neisseria species were obtained from the NCBI database. Comparative analysis of each PenA sequence was performed by multiple sequence alignment using ClustalW. In vitro resistance acquisition experiments were conducted to investigate the possibility of selection pressure by cefixime-induced amino acid substitution mutations in PenA.

RESULTS

All N. subflava strains, including two with low susceptibility to expanded-spectrum cephalosporins (ESCs), possessed the majority of the PenA FC428 sequence. Furthermore, a number of strains, but not all, of closely related species of N. subflava showed similar results. PenA FC428 sequences were also found in some strains of distantly related species. No new mutations in the penA sequence were observed in colonies with increased MIC in in vitro resistance acquisition experiments.

CONCLUSIONS

This study provides strong evidence that the FC428 PenA mosaic sequence originated from N. subflava and related species among oral commensal Neisseria species. The results of in vitro resistance acquisition experiments also suggested that one of the PenA FC428-like sequence gene polymorphisms resulted in the expression of ESC resistance. Furthermore, many of the PenA FC428 mosaic sequences were thought to be involved in the so-called epistasis effect that regulates the expression of resistance, without directly contributing to the resistance level itself.

Chance Favors the Prepared Genomes: Horizontal Transfer Shapes the Emergence of Antibiotic Resistance Mutations in Core Genes

Bacterial lineages acquire novel traits at diverse rates in part because the genetic background impacts the successful acquisition of novel genes by horizontal transfer. Yet, how horizontal transfer affects the subsequent evolution of core genes remains poorly understood. Here, we studied the evolution of resistance to quinolones in Escherichia coli accounting for population structure. We found 60 groups of genes whose gain or loss induced an increase in the probability of subsequently becoming resistant to quinolones by point mutations in the gyrase and topoisomerase genes. These groups include functions known to be associated with direct mitigation of the effect of quinolones, with metal uptake, cell growth inhibition, biofilm formation, and sugar metabolism. Many of them are encoded in phages or plasmids. Although some of the chronologies may reflect epidemiological trends, many of these groups encoded functions providing latent phenotypes of antibiotic low-level resistance, tolerance, or persistence under quinolone treatment. The mutations providing resistance were frequent and accumulated very quickly. Their emergence was found to increase the rate of acquisition of other antibiotic resistances setting the path for multidrug resistance. Hence, our findings show that horizontal gene transfer shapes the subsequent emergence of adaptive mutations in core genes. In turn, these mutations further affect the subsequent evolution of resistance by horizontal gene transfer. Given the substantial gene flow within bacterial genomes, interactions between horizontal transfer and point mutations in core genes may be a key to the success of adaptation processes.

An Amino Acid Substitution in Elongation Factor EF-G1A Alters the Antibiotic Susceptibility of Pseudomonas aeruginosa LasR-Null Mutants

The opportunistic bacterium Pseudomonas aeruginosa uses the LasR-I quorum-sensing system to increase resistance to the aminoglycoside antibiotic tobramycin. Paradoxically, lasR-null mutants are commonly isolated from chronic human infections treated with tobramycin, suggesting there may be a mechanism that permits the emergence of lasR-null mutants under tobramycin selection. We hypothesized that some other genetic mutations that emerge in these isolates might modulate the effects of lasR-null mutations on antibiotic resistance. To test this hypothesis, we inactivated lasR in several highly tobramycin-resistant isolates from long-term evolution experiments. In some of these isolates, inactivating lasR further increased resistance, compared with decreasing resistance of the wild-type ancestor. These strain-dependent effects were due to a G61A nucleotide polymorphism in the fusA1 gene encoding amino acid substitution A21T in the translation elongation factor EF-G1A. The EF-G1A mutational effects required the MexXY efflux pump and the MexXY regulator ArmZ. The fusA1 mutation also modulated ΔlasR mutant resistance to two other antibiotics, ciprofloxacin and ceftazidime. Our results identify a gene mutation that can reverse the direction of the antibiotic selection of lasR mutants, a phenomenon known as sign epistasis, and provide a possible explanation for the emergence of lasR-null mutants in clinical isolates. IMPORTANCE One of the most common mutations in Pseudomonas aeruginosa clinical isolates is in the quorum sensing lasR gene. In laboratory strains, lasR disruption decreases resistance to the clinical antibiotic tobramycin. To understand how lasR mutations emerge in tobramycin-treated patients, we mutated lasR in highly tobramycin-resistant laboratory strains and determined the effects on resistance. Disrupting lasR enhanced the resistance of some strains. These strains had a single amino acid substitution in the translation factor EF-G1A. The EF-G1A mutation reversed the selective effects of tobramycin on lasR mutants. These results illustrate how adaptive mutations can lead to the emergence of new traits in a population and are relevant to understanding how genetic diversity contributes to the progression of disease during chronic infections.

Trends and challenges of multi-drug resistance in childhood tuberculosis

Drug-resistant tuberculosis (DR-TB) in children is a growing global health concern, This review provides an overview of the current epidemiology of childhood TB and DR-TB, including prevalence, incidence, and mortality. We discuss the challenges in diagnosing TB and DR-TB in children and the limitations of current diagnostic tools. We summarize the challenges associated with treating multi-drug resistance TB in childhood, including limitations of current treatment options, drug adverse effects, prolonged regimens, and managing and monitoring during treatment. We highlight the urgent need for improved diagnosis and treatment of DR-TB in children. The treatment of children with multidrug-resistant tuberculosis will be expanded to include the evaluation of new drugs or new combinations of drugs. Basic research is needed to support the technological development of biomarkers to assess the phase of therapy, as well as the urgent need for improved diagnostic and treatment options.

Tobramycin adaptation alters the antibiotic susceptibility of Pseudomonas aeruginosa quorum sensing-null mutants

The opportunistic bacterium Pseudomonas aeruginosa uses the LasR-I quorum sensing system to increase resistance to the aminioglycoside antibiotic tobramycin. Paradoxically, lasR-null mutants are commonly isolated from chronic human infections treated with tobramycin, suggesting there may be a mechanism allowing the lasR-null mutants to persist under tobramycin selection. We hypothesized that the effects of inactivating lasR on tobramycin resistance might be dependent on the presence or absence of other gene mutations in that strain, a phenomenon known as epistasis. To test this hypothesis, we inactivated lasR in several highly tobramycin-resistant isolates from long-term evolution experiments. We show that the effects of ΔlasR on tobramycin resistance are strain dependent. The effects can be attributed to a point mutation in the gene encoding the translation elongation factor fusA1 (G61A nucleotide substitution), which confers a strong selective advantage to lasR-null PA14 under tobramycin selection. This fusA1 G61A mutation results in increased activity of the MexXY efflux pump and expression of the mexXY regulator ArmZ. The fusA1 mutation can also modulate ΔlasR mutant resistance to two other antibiotics, ciprofloxacin and ceftazidime. Our results demonstrate the importance of epistatic gene interactions on antibiotic susceptibility of lasR-null mutants. These results support of the idea that gene interactions might play a significant role in the evolution of quorum sensing in P. aeruginosa.

Environmental complexity is more important than mutation in driving the evolution of latent novel traits in E. coli

Recent experiments show that adaptive Darwinian evolution in one environment can lead to the emergence of multiple new traits that provide no immediate benefit in this environment. Such latent non-adaptive traits, however, can become adaptive in future environments. We do not know whether mutation or environment-driven selection is more important for the emergence of such traits. To find out, we evolve multiple wild-type and mutator E. coli populations under two mutation rates in simple (single antibiotic) environments and in complex (multi-antibiotic) environments. We then assay the viability of evolved populations in dozens of new environments and show that all populations become viable in multiple new environments different from those they had evolved in. The number of these new environments increases with environmental complexity but not with the mutation rate. Genome sequencing demonstrates the reason: Different environments affect pleiotropic mutations differently. Our experiments show that the selection pressure provided by an environment can be more important for the evolution of novel traits than the mutational supply experienced by a wild-type and a mutator strain of E. coli.

Novel traits without immediate fitness benefit evolve frequently but we don’t know whether mutation or environment-driven selection drives this evolution. Here, using experimental evolution of E. coli populations, the authors demonstrate the importance of selection in the evolution of latent novel traits.

Intragenic compensation through the lens of deep mutational scanning

A significant fraction of mutations in proteins are deleterious and result in adverse consequences for protein function, stability, or interaction with other molecules. Intragenic compensation is a specific case of positive epistasis when a neutral missense mutation cancels effect of a deleterious mutation in the same protein. Permissive compensatory mutations facilitate protein evolution, since without them all sequences would be extremely conserved. Understanding compensatory mechanisms is an important scientific challenge at the intersection of protein biophysics and evolution. In human genetics, intragenic compensatory interactions are important since they may result in variable penetrance of pathogenic mutations or fixation of pathogenic human alleles in orthologous proteins from related species. The latter phenomenon complicates computational and clinical inference of an allele's pathogenicity. Deep mutational scanning is a relatively new technique that enables experimental studies of functional effects of thousands of mutations in proteins. We review the important aspects of the field and discuss existing limitations of current datasets. We reviewed ten published DMS datasets with quantified functional effects of single and double mutations and described rates and patterns of intragenic compensation in eight of them.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01005-w.

Association between Antimicrobial Consumption and the Prevalence of Nosocomial Carbapenem-Resistant Escherichia coli and Klebsiella pneumoniae in a Tertiary Hospital in Northern Taiwan

Carbapenem-resistant Enterobacteriales has become a threat in Taiwan. This is the first local study focusing on the association between carbapenem-resistant Enterobacteriales and antimicrobial consumption. From January 2012 to December 2020, data were collected in a tertiary care hospital in Taipei, Taiwan. Antimicrobial consumption was estimated by the defined daily dose/1,000 patient-days. During the same period, the prevalence of carbapenem-resistant Escherichia coli (CREC) and carbapenem-resistant Klebsiella pneumoniae (CRKP) were collected through routine surveillance data. The following retrospective analyses were conducted: 1) analysis of antimicrobial consumption over time, (2) analysis and forecast of CREC and CRKP prevalence over time, and 3) analysis of correlation between antimicrobial consumption and the prevalence of CREC and CRKP. The consumption of piperacillin/tazobactam (β = 0.615), fluoroquinolones (β = 0.856), meropenem (β = 0.819), and doripenem (β = 0.891) increased during the observation period (P < 0.001), and the consumption of aminoglycosides (β = -0.852) and imipenem/cilastatin (β = -0.851) decreased (P < 0.001). The prevalence of CRKP rose over time (β = 0.522, P = 0.001) and correlated positively with the consumption of fluoroquinolones, levofloxacin, penicillin/β-lactamase inhibitor, piperacillin/tazobactam, meropenem, and doripenem (P < 0.05). The prevalence of CRKP and CREC both correlated negatively with consumption of aminoglycosides (P < 0.01). The prevalence of CRKP in our hospital increased as the forecast predicted based on an autoregressive integrated moving average model. This study provides alarming messages for members participating in antimicrobial stewardship programs, including the increasing prevalence of CRKP, the increasing consumption of broad-spectrum antibiotics, and the positive correlation between them.

Machine Learning for Antimicrobial Resistance Prediction: Current Practice, Limitations, and Clinical Perspective

Antimicrobial resistance (AMR) is a global health crisis that poses a great threat to modern medicine. Effective prevention strategies are urgently required to slow the emergence and further dissemination of AMR. Given the availability of data sets encompassing hundreds or thousands of pathogen genomes, machine learning (ML) is increasingly being used to predict resistance to different antibiotics in pathogens based on gene content and genome composition. A key objective of this work is to advocate for the incorporation of ML into front-line settings but also highlight the further refinements that are necessary to safely and confidently incorporate these methods. The question of what to predict is not trivial given the existence of different quantitative and qualitative laboratory measures of AMR. ML models typically treat genes as independent predictors, with no consideration of structural and functional linkages; they also may not be accurate when new mutational variants of known AMR genes emerge. Finally, to have the technology trusted by end users in public health settings, ML models need to be transparent and explainable to ensure that the basis for prediction is clear. We strongly advocate that the next set of AMR-ML studies should focus on the refinement of these limitations to be able to bridge the gap to diagnostic implementation.

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.

The Neglected Contribution of Streptomycin to the Tuberculosis Drug Resistance Problem

The airborne pathogen Mycobacterium tuberculosis is responsible for a present major public health problem worsened by the emergence of drug resistance. M. tuberculosis has acquired and developed streptomycin (STR) resistance mechanisms that have been maintained and transmitted in the population over the last decades. Indeed, STR resistant mutations are frequently identified across the main M. tuberculosis lineages that cause tuberculosis outbreaks worldwide. The spread of STR resistance is likely related to the low impact of the most frequent underlying mutations on the fitness of the bacteria. The withdrawal of STR from the first-line treatment of tuberculosis potentially lowered the importance of studying STR resistance. However, the prevalence of STR resistance remains very high, could be underestimated by current genotypic methods, and was found in outbreaks of multi-drug (MDR) and extensively drug (XDR) strains in different geographic regions. Therefore, the contribution of STR resistance to the problem of tuberculosis drug resistance should not be neglected. Here, we review the impact of STR resistance and detail well-known and novel candidate STR resistance mechanisms, genes, and mutations. In addition, we aim to provide insights into the possible role of STR resistance in the development of multi-drug resistant tuberculosis.

Resistance Trend Estimation Using Regression Analysis to Enhance Antimicrobial Surveillance: A Multi-Centre Study in London 2009-2016

In the last years, there has been an increase of antimicrobial resistance rates around the world with the misuse and overuse of antimicrobials as one of the main leading drivers. In response to this threat, a variety of initiatives have arisen to promote the efficient use of antimicrobials. These initiatives rely on antimicrobial surveillance systems to promote appropriate prescription practices and are provided by national or global health care institutions with limited consideration of the variations within hospitals. As a consequence, physicians' adherence to these generic guidelines is still limited. To fill this gap, this work presents an automated approach to performing local antimicrobial surveillance from microbiology data. Moreover, in addition to the commonly reported resistance rates, this work estimates secular resistance trends through regression analysis to provide a single value that effectively communicates the resistance trend to a wider audience. The methods considered for trend estimation were ordinary least squares regression, weighted least squares regression with weights inversely proportional to the number of microbiology records available and autoregressive integrated moving average. Among these, weighted least squares regression was found to be the most robust against changes in the granularity of the time series and presented the best performance. To validate the results, three case studies have been thoroughly compared with the existing literature: (i) Escherichia coli in urine cultures; (ii) Escherichia coli in blood cultures; and (iii) Staphylococcus aureus in wound cultures. The benefits of providing local rather than general antimicrobial surveillance data of a higher quality is two fold. Firstly, it has the potential to stimulate engagement among physicians to strengthen their knowledge and awareness on antimicrobial resistance which might encourage prescribers to change their prescription habits more willingly. Moreover, it provides fundamental knowledge to the wide range of stakeholders to revise and potentially tailor existing guidelines to the specific needs of each hospital.

Structural and Genomic Insights Into Pyrazinamide Resistance in Mycobacterium tuberculosis Underlie Differences Between Ancient and Modern Lineages

Resistance to drugs used to treat tuberculosis disease (TB) continues to remain a public health burden, with missense point mutations in the underlying Mycobacterium tuberculosis bacteria described for nearly all anti-TB drugs. The post-genomics era along with advances in computational and structural biology provide opportunities to understand the interrelationships between the genetic basis and the structural consequences of M. tuberculosis mutations linked to drug resistance. Pyrazinamide (PZA) is a crucial first line antibiotic currently used in TB treatment regimens. The mutational promiscuity exhibited by the pncA gene (target for PZA) necessitates computational approaches to investigate the genetic and structural basis for PZA resistance development. We analysed 424 missense point mutations linked to PZA resistance derived from ∼35K M. tuberculosis clinical isolates sourced globally, which comprised the four main M. tuberculosis lineages (Lineage 1-4). Mutations were annotated to reflect their association with PZA resistance. Genomic measures (minor allele frequency and odds ratio), structural features (surface area, residue depth and hydrophobicity) and biophysical effects (change in stability and ligand affinity) of point mutations on pncA protein stability and ligand affinity were assessed. Missense point mutations within pncA were distributed throughout the gene, with the majority (>80%) of mutations with a destabilising effect on protomer stability and on ligand affinity. Active site residues involved in PZA binding were associated with multiple point mutations highlighting mutational diversity due to selection pressures at these functionally important sites. There were weak associations between genomic measures and biophysical effect of mutations. However, mutations associated with PZA resistance showed statistically significant differences between structural features (surface area and residue depth), but not hydrophobicity score for mutational sites. Most interestingly M. tuberculosis lineage 1 (ancient lineage) exhibited a distinct protein stability profile for mutations associated with PZA resistance, compared to modern lineages.

Compatibility of Evolutionary Responses to Constituent Antibiotics Drive Resistance Evolution to Drug Pairs

Antibiotic combinations are considered a relevant strategy to tackle the global antibiotic resistance crisis since they are believed to increase treatment efficacy and reduce resistance evolution (WHO treatment guidelines for drug-resistant tuberculosis: 2016 update.). However, studies of the evolution of bacterial resistance to combination therapy have focused on a limited number of drugs and have provided contradictory results (Lipsitch, Levin BR. 1997; Hegreness et al. 2008; Munck et al. 2014). To address this gap in our understanding, we performed a large-scale laboratory evolution experiment, adapting eight replicate lineages of Escherichia coli to a diverse set of 22 different antibiotics and 33 antibiotic pairs. We found that combination therapy significantly limits the evolution of de novode novo resistance in E. coli, yet different drug combinations vary substantially in their propensity to select for resistance. In contrast to current theories, the phenotypic features of drug pairs are weak predictors of resistance evolution. Instead, the resistance evolution is driven by the relationship between the evolutionary trajectories that lead to resistance to a drug combination and those that lead to resistance to the component drugs. Drug combinations requiring a novel genetic response from target bacteria compared with the individual component drugs significantly reduce resistance evolution. These data support combination therapy as a treatment option to decelerate resistance evolution and provide a novel framework for selecting optimized drug combinations based on bacterial evolutionary responses.

Predicting plasmid persistence in microbial communities by coarse-grained modeling

Plasmids are a major type of mobile genetic elements (MGEs) that mediate horizontal gene transfer. The stable maintenance of plasmids plays a critical role in the functions and survival for microbial populations. However, predicting and controlling plasmid persistence and abundance in complex microbial communities remain challenging. Computationally, this challenge arises from the combinatorial explosion associated with the conventional modeling framework. Recently, a plasmid-centric framework (PCF) has been developed to overcome this computational bottleneck. This framework enables the derivation of a simple metric, the persistence potential, to predict plasmid persistence and abundance. Here, we discuss how PCF can be extended to account for plasmid interactions. We also discuss how such model-guided predictions of plasmid fates can benefit from the development of new experimental tools and data-driven computational methods.

Adaptive resistance mutations at supra-inhibitory concentrations independent of SOS mutagenesis

Emergence of resistant bacteria during antimicrobial treatment is one of the most critical and universal health threats. It is known that several stress-induced mutagenesis and heteroresistance mechanisms can enhance microbial adaptation to antibiotics. Here, we demonstrate that the pathogen Bartonella can undergo stress-induced mutagenesis despite the fact it lacks error-prone polymerases, the rpoS gene and functional UV-induced mutagenesis. We demonstrate that Bartonella acquire de novo single mutations during rifampicin exposure at suprainhibitory concentrations at a much higher rate than expected from spontaneous fluctuations. This is while exhibiting a minimal heteroresistance capacity. The emerged resistant mutants acquired a single rpoB mutation, whereas no other mutations were found in their whole genome. Interestingly, the emergence of resistance in Bartonella occurred only during gradual exposure to the antibiotic, indicating that Bartonella sense and react to the changing environment. Using a mathematical model, we demonstrated that, to reproduce the experimental results, mutation rates should be transiently increased over 1,000-folds, and a larger population size or greater heteroresistance capacity is required. RNA expression analysis suggests that the increased mutation rate is due to downregulation of key DNA repair genes (mutS, mutY, and recA), associated with DNA breaks caused by massive prophage inductions. These results provide new evidence of the hazard of antibiotic overuse in medicine and agriculture.

Predicting Antimicrobial Resistance Using Partial Genome Alignments

Antimicrobial resistance (AMR) is an important global health threat that impacts millions of people worldwide each year. Developing methods that can detect and predict AMR phenotypes can help to mitigate the spread of AMR by informing clinical decision making and appropriate mitigation strategies. Many bioinformatic methods have been developed for predicting AMR phenotypes from whole-genome sequences and AMR genes, but recent studies have indicated that predictions can be made from incomplete genome sequence data. In order to more systematically understand this, we built random forest-based machine learning classifiers for predicting susceptible and resistant phenotypes for Klebsiella pneumoniae (1,640 strains), Mycobacterium tuberculosis (2,497 strains), and Salmonella enterica (1,981 strains). We started by building models from alignments that were based on a reference chromosome for each species. We then subsampled each chromosomal alignment and built models for the resulting subalignments, finding that very small regions, representing approximately 0.1 to 0.2% of the chromosome, are predictive. In K. pneumoniae, M. tuberculosis, and S. enterica, the subalignments are able to predict multiple AMR phenotypes with at least 70% accuracy, even though most do not encode an AMR-related function. We used these models to identify regions of the chromosome with high and low predictive signals. Finally, subalignments that retain high accuracy across larger phylogenetic distances were examined in greater detail, revealing genes and intergenic regions with potential links to AMR, virulence, transport, and survival under stress conditions.

IMPORTANCE Antimicrobial resistance causes thousands of deaths annually worldwide. Understanding the regions of the genome that are involved in antimicrobial resistance is important for developing mitigation strategies and preventing transmission. Machine learning models are capable of predicting antimicrobial resistance phenotypes from bacterial genome sequence data by identifying resistance genes, mutations, and other correlated features. They are also capable of implicating regions of the genome that have not been previously characterized as being involved in resistance. In this study, we generated global chromosomal alignments for Klebsiella pneumoniae, Mycobacterium tuberculosis, and Salmonella enterica and systematically searched them for small conserved regions of the genome that enable the prediction of antimicrobial resistance phenotypes. In addition to known antimicrobial resistance genes, this analysis identified genes involved in virulence and transport functions, as well as many genes with no previous implication in antimicrobial resistance.

Collateral Sensitivity to β-Lactam Drugs in Drug-Resistant Tuberculosis Is Driven by the Transcriptional Wiring of BlaI Operon Genes

The evolution of resistance to one antimicrobial can result in enhanced sensitivity to another, known as "collateral sensitivity." This underexplored phenomenon opens new therapeutic possibilities for patients infected with pathogens unresponsive to classical treatments. Intrinsic resistance to β-lactams in Mycobacterium tuberculosis (the causative agent of tuberculosis) has traditionally curtailed the use of these low-cost and easy-to-administer drugs for tuberculosis treatment. Recently, β-lactam sensitivity has been reported in strains resistant to classical tuberculosis therapy, resurging the interest in β-lactams for tuberculosis. However, a lack of understanding of the molecular underpinnings of this sensitivity has delayed exploration in the clinic. We performed gene expression and network analyses and in silico knockout simulations of genes associated with β-lactam sensitivity and genes associated with resistance to classical tuberculosis drugs to investigate regulatory interactions and identify key gene mediators. We found activation of the key inhibitor of β-lactam resistance, blaI, following classical drug treatment as well as transcriptional links between genes associated with β-lactam sensitivity and those associated with resistance to classical treatment, suggesting that regulatory links might explain collateral sensitivity to β-lactams. Our results support M. tuberculosis β-lactam sensitivity as a collateral consequence of the evolution of resistance to classical tuberculosis drugs, mediated through changes to transcriptional regulation. These findings support continued exploration of β-lactams for the treatment of patients infected with tuberculosis strains resistant to classical therapies. IMPORTANCE Tuberculosis remains a significant cause of global mortality, with strains resistant to classical drug treatment considered a major health concern by the World Health Organization. Challenging treatment regimens and difficulty accessing drugs in low-income communities have led to a high prevalence of strains resistant to multiple drugs, making the development of alternative therapies a priority. Although Mycobacterium tuberculosis is naturally resistant to β-lactam drugs, previous studies have shown sensitivity in strains resistant to classical drug treatment, but we currently lack understanding of the molecular underpinnings behind this phenomenon. We found that genes involved in β-lactam susceptibility are activated after classical drug treatment resulting from tight regulatory links with genes involved in drug resistance. Our study supports the hypothesis that β-lactam susceptibility observed in drug-resistant strains results from the underlying regulatory network of M. tuberculosis, supporting further exploration of the use of β-lactams for tuberculosis treatment.

Mastering the control of the Rho transcription factor for biotechnological applications

The present review represents an update on the fundamental role played by the Rho factor, which facilitates the process of Rho-dependent transcription termination in the prokaryotic world; it also provides a summary of relevant mutations in the Rho factor and the insights they provide into the functions carried out by this protein. Furthermore, a section is dedicated to the putative future use of Rho (the 'taming' of Rho) to facilitate biotechnological processes and adapt them to different technological contexts. Novel bacterial strains can be designed, containing mutations in the rho gene, that are better suited for different biotechnological applications. This process can obtain novel microbial strains that are adapted to lower temperatures of fermentation, shorter production times, exhibit better nutrient utilization, or display other traits that are beneficial in productive Biotechnology. Additional important issues reviewed here include epistasis, the design of TATA boxes, the role of small RNAs, and the manipulation of clathrin-mediated endocytosis, by some pathogenic bacteria, to invade eukaryotic cells. KEY POINTS: • It is postulated that controlling the action of the prokaryotic Rho factor could generate major biotechnological improvements, such as an increase in bacterial productivity or a reduction of the microbial-specific growth rate. • The review also evaluates the putative impact of epistatic mechanisms on Biotechnology, both as possible responsible for unexpected failures in gene cloning and more important for the genesis of new strains for biotechnological applications • The use of clathrin-coated vesicles by intracellular bacterial microorganisms is included too and proposed as a putative delivery mechanism, for drugs and vaccines.

Genetic interaction analysis in microbial pathogens: unravelling networks of pathogenesis, antimicrobial susceptibility and host interactions

Genetic interaction (GI) analysis is a powerful genetic strategy that analyzes the fitness and phenotypes of single- and double-gene mutant cells in order to dissect the epistatic interactions between genes, categorize genes into biological pathways, and characterize genes of unknown function. GI analysis has been extensively employed in model organisms for foundational, systems-level assessment of the epistatic interactions between genes. More recently, GI analysis has been applied to microbial pathogens and has been instrumental for the study of clinically important infectious organisms. Here, we review recent advances in systems-level GI analysis of diverse microbial pathogens, including bacterial and fungal species. We focus on important applications of GI analysis across pathogens, including GI analysis as a means to decipher complex genetic networks regulating microbial virulence, antimicrobial drug resistance and host-pathogen dynamics, and GI analysis as an approach to uncover novel targets for combination antimicrobial therapeutics. Together, this review bridges our understanding of GI analysis and complex genetic networks, with applications to diverse microbial pathogens, to further our understanding of virulence, the use of antimicrobial therapeutics and host-pathogen interactions. .

Emerging macrolide resistance in Bordetella pertussis in mainland China: Findings and warning from the global pertussis initiative

Whooping cough, or pertussis, is a highly communicable infectious disease caused by the bacterium Bordetella pertussis. Vaccination once reduced the incidence of the disease, but a global resurgence of the infection happened during the past two decades, likely due to the waning immunity of vaccination. Macrolides such as erythromycin and azithromycin are the drugs of primary choice for treatment. In this personal view, we call for attention to macrolide-resistant B. pertussis (MRBP), which has emerged and prevailed in mainland China for years and are exclusively mediated by mutations in the 23S rRNA gene. Whether the prevalence of MRBP in China results from overuse of azithromycin in clinical medicine remains unknown. The incidence of MRBP is low in other countries, but this could be a technical illusion since China employs culture as the mainstream diagnostic method whereas nucleic-acid amplification test being widely used in other countries fail to test antimicrobial susceptibility. Given the increasingly frequent global travel that facilitates microbial transmission worldwide, there is a pressing need to perform international surveillance on MRBP to prevent the potential circulation of the organism. Finding alternative agents that possess good activity against B. pertussis is also urgently required.

Amino Acid k-mer Feature Extraction for Quantitative Antimicrobial Resistance (AMR) Prediction by Machine Learning and Model Interpretation for Biological Insights

Machine learning algorithms can learn mechanisms of antimicrobial resistance from the data of DNA sequence without any a priori information. Interpreting a trained machine learning algorithm can be exploited for validating the model and obtaining new information about resistance mechanisms. Different feature extraction methods, such as SNP calling and counting nucleotide k-mers have been proposed for presenting DNA sequences to the model. However, there are trade-offs between interpretability, computational complexity and accuracy for different feature extraction methods. In this study, we have proposed a new feature extraction method, counting amino acid k-mers or oligopeptides, which provides easier model interpretation compared to counting nucleotide k-mers and reaches the same or even better accuracy in comparison with different methods. Additionally, we have trained machine learning algorithms using different feature extraction methods and compared the results in terms of accuracy, model interpretability and computational complexity. We have built a new feature selection pipeline for extraction of important features so that new AMR determinants can be discovered by analyzing these features. This pipeline allows the construction of models that only use a small number of features and can predict resistance accurately.

Tackling Antibiotic Resistance with Compounds of Natural Origin: A Comprehensive Review

Drug-resistant bacteria pose a serious threat to human health worldwide. Current antibiotics are losing efficacy and new antimicrobial agents are urgently needed. Living organisms are an invaluable source of antimicrobial compounds. The antimicrobial activity of the most representative natural products of animal, bacterial, fungal and plant origin are reviewed in this paper. Their activity against drug-resistant bacteria, their mechanisms of action, the possible development of resistance against them, their role in current medicine and their future perspectives are discussed. Electronic databases such as PubMed, Scopus and ScienceDirect were used to search scientific contributions until September 2020, using relevant keywords. Natural compounds of heterogeneous origins have been shown to possess antimicrobial capabilities, including against antibiotic-resistant bacteria. The most commonly found mechanisms of antimicrobial action are related to protein biosynthesis and alteration of cell walls and membranes. Various natural compounds, especially phytochemicals, have shown synergistic capacity with antibiotics. There is little literature on the development of specific resistance mechanisms against natural antimicrobial compounds. New technologies such as -omics, network pharmacology and informatics have the potential to identify and characterize new natural antimicrobial compounds in the future. This knowledge may be useful for the development of future therapeutic strategies.

Transcriptional Organization of the Stability Module of Broad-Host-Range Plasmid RA3, from the IncU Group

The broad-host-range (BHR) conjugative plasmids have developed diverse adaptive mechanisms defining the range of their promiscuity. The BHR conjugative RA3 plasmid, the archetype of the IncU group, can transfer between, replicate in, and be maintained in representatives of Alpha-, Beta-, and Gammaproteobacteria Its stability module encompasses ten open reading frames (ORFs) apparently organized into five operons, all transcribed in the same direction from several strong promoters that are tightly regulated either by autorepressors or by global plasmid-encoded regulators. In this paper, we demonstrate that owing to an efficient RNA polymerase (RNAP) read-through, the transcription from the first promoter, orf02p, may continue through the whole module. Moreover, an analysis of mRNA produced from the wild-type (WT) stability module and its deletion variants deprived of particular internal transcription initiation sites reveals that in fact each operon may be transcribed from any upstream promoter, giving rise to multicistronic transcripts of variable length and creating an additional level of gene expression control by transcript dosage adjustment. The gene expression patterns differ among various hosts, indicating that promoter recognition, regulation, and the RNAP read-through mechanisms are modulated in a species-specific manner.IMPORTANCE The efficiently disseminating conjugative or mobilizable BHR plasmids play key roles in the horizontal spread of genetic information between closely related and phylogenetically distant species, which can be harmful from the medical, veterinary, or industrial point of view. Understanding the mechanisms determining the plasmid's ability to function in diverse hosts is essential to help limit the spread of undesirable plasmid-encoded traits, e.g., antibiotic resistance. The range of a plasmid's promiscuity depends on the adaptations of its transfer, replication, and stability functions to the various hosts. IncU plasmids, with the archetype plasmid RA3, are considered to constitute a reservoir of antibiotic resistance genes in aquatic environments; however, the molecular mechanisms determining their adaptability to a broad range of hosts are rather poorly characterized. Here, we present the transcriptional organization of the stability module and show that the gene transcript dosage effect is an important determinant of the stable maintenance of RA3 in different hosts.

Factors favouring the evolution of multidrug resistance in bacteria

The evolution of multidrug antibiotic resistance in commensal bacteria is an important public health concern. Commensal bacteria such as Escherichia coli, Streptococcus pneumoniae or Staphylococcus aureus, are also opportunistic pathogens causing a large fraction of the community-acquired and hospital-acquired bacterial infections. Multidrug resistance (MDR) makes these infections harder to treat with antibiotics and may thus cause substantial additional morbidity and mortality. Here, we develop an evolutionary epidemiology model to identify the factors favouring the evolution of MDR in commensal bacteria. The model describes the evolution of antibiotic resistance in a commensal bacterial species evolving in a host population subjected to multiple antibiotic treatments. We combine statistical analysis of a large number of simulations and mathematical analysis to understand the model behaviour. We find that MDR evolves more readily when it is less costly than expected from the combinations of single resistances (positive epistasis). MDR frequently evolves when bacteria are in contact with multiple drugs prescribed in the host population, even if individual hosts are only treated with a single drug at a time. MDR is favoured when the host population is structured in different classes that vary in their rates of antibiotic treatment. However, under most circumstances, recombination between loci involved in resistance does not meaningfully affect the equilibrium frequency of MDR. Together, these results suggest that MDR is a frequent evolutionary outcome in commensal bacteria that encounter the variety of antibiotics prescribed in the host population. A better characterization of the variability in antibiotic use across the host population (e.g. across age classes or geographical location) would help predict which MDR genotypes will most readily evolve.

Occurrence, virulence genes, and antimicrobial profiles of Escherichia coli O157 isolated from ruminants slaughtered in Al Ain, United Arab Emirates

BACKGROUND

Shiga toxin-producing Escherichia coli (STEC) is a major source of food-borne illness around the world. E. coli O157 has been widely reported as the most common STEC serogroup and has emerged as an important enteric pathogen. Cattle, in particular have been identified as a major E. coli O157:H7 reservoir of human infections; however, the prevalence of this organism in camels, sheep, and goats is less understood. The aim of this study was to evaluate the occurrence and concentration of E. coli serotype O157 in the feces of healthy camels (n = 140), cattle (n = 137), sheep (n = 141) and goats (n = 150) slaughtered in United Arab Emirates (UAE) for meat consumption between September 2017 and August 2018. We used immunomagnetic separation coupled with a culture-plating method to detect E. coli O157. Non-sorbitol fermenting colonies were assessed via latex-agglutination testing, and positive cultures were analyzed by performing polymerase chain reactions to detect genes encoding attaching and effacing protein (eaeA), hemolysin A (hlyA, also known as ehxA) and Shiga toxin (stx1 and stx2), and E. coli O157:H7 specific genes (rfb O157, uidA, and fliC). All E. coli O157 isolates were analyzed for their susceptibility to 20 selected antimicrobials.

RESULTS

E. coli O157 was observed in camels, goats, and cattle fecal samples at abundances of 4.3, 2, and 1.46%, respectively, but it was undetectable in sheep feces. The most prevalent E. coli O157 gene in all STEC isolates was stx₂;_, whereas, stx₁ was not detected in any of the samples. The fecal samples from camels, goats, and cattle harbored E. coli O157 isolates that were 100% susceptible to cefotaxime, chloramphenicol, ciprofloxacin, norfloxacin, and polymyxin B.

CONCLUSION

To our knowledge, this is the first report on the occurrence of E. coli O157 in slaughter animals in the UAE. Our results clearly demonstrate the presence of E. coli O157 in slaughtered animals, which could possibly contaminate meat products intended for human consumption.

Candida albicans Genetic Background Influences Mean and Heterogeneity of Drug Responses and Genome Stability during Evolution in Fluconazole

Antimicrobial resistance is an evolutionary phenomenon with clinical implications. We tested how replicates from diverse strains of Candida albicans, a prevalent human fungal pathogen, evolve in the commonly prescribed antifungal drug fluconazole. Replicates on average increased in fitness in the level of drug they were evolved to, with the least fit parental strains improving the most. Very few replicates increased resistance above the drug level they were evolved in. Notably, many replicates increased in genome size and changed in drug tolerance (a drug response where a subpopulation of cells grow slowly in high levels of drug), and variability among replicates in fitness, tolerance, and genome size was higher in strains that initially were more sensitive to the drug. Genetic background influenced the average degree of adaptation and the evolved variability of many phenotypes, highlighting that different strains from the same species may respond and adapt very differently during adaptation.

The importance of within-species diversity in determining the evolutionary potential of a population to evolve drug resistance or tolerance is not well understood, including in eukaryotic pathogens. To examine the influence of genetic background, we evolved replicates of 20 different clinical isolates of Candida albicans, a human fungal pathogen, in fluconazole, the commonly used antifungal drug. The isolates hailed from the major C. albicans clades and had different initial levels of drug resistance and tolerance to the drug. The majority of replicates rapidly increased in fitness in the evolutionary environment, with the degree of improvement inversely correlated with parental strain fitness in the drug. Improvement was largely restricted to up to the evolutionary level of drug: only 4% of the evolved replicates increased resistance (MIC) above the evolutionary level of drug. Prevalent changes were altered levels of drug tolerance (slow growth of a subpopulation of cells at drug concentrations above the MIC) and increased diversity of genome size. The prevalence and predominant direction of these changes differed in a strain-specific manner, but neither correlated directly with parental fitness or improvement in fitness. Rather, low parental strain fitness was correlated with high levels of heterogeneity in fitness, tolerance, and genome size among evolved replicates. Thus, parental strain background is an important determinant in mean improvement to the evolutionary environment as well as the diversity of evolved phenotypes, and the range of possible responses of a pathogen to an antimicrobial drug cannot be captured by in-depth study of a single strain background.

IMPORTANCE Antimicrobial resistance is an evolutionary phenomenon with clinical implications. We tested how replicates from diverse strains of Candida albicans, a prevalent human fungal pathogen, evolve in the commonly prescribed antifungal drug fluconazole. Replicates on average increased in fitness in the level of drug they were evolved to, with the least fit parental strains improving the most. Very few replicates increased resistance above the drug level they were evolved in. Notably, many replicates increased in genome size and changed in drug tolerance (a drug response where a subpopulation of cells grow slowly in high levels of drug), and variability among replicates in fitness, tolerance, and genome size was higher in strains that initially were more sensitive to the drug. Genetic background influenced the average degree of adaptation and the evolved variability of many phenotypes, highlighting that different strains from the same species may respond and adapt very differently during adaptation.

Hypermutator Pseudomonas aeruginosa Exploits Multiple Genetic Pathways To Develop Multidrug Resistance during Long-Term Infections in the Airways of Cystic Fibrosis Patients

Pseudomonas aeruginosa exploits intrinsic and acquired resistance mechanisms to resist almost every antibiotic used in chemotherapy. Antimicrobial resistance in P. aeruginosa isolates recovered from cystic fibrosis (CF) patients is further enhanced by the occurrence of hypermutator strains, a hallmark of chronic infections in CF patients. However, the within-patient genetic diversity of P. aeruginosa populations related to antibiotic resistance remains unexplored. Here, we show the evolution of the mutational resistome profile of a P. aeruginosa hypermutator lineage by performing longitudinal and transversal analyses of isolates collected from a CF patient throughout 20 years of chronic infection. Our results show the accumulation of thousands of mutations, with an overall evolutionary history characterized by purifying selection. However, mutations in antibiotic resistance genes appear to have been positively selected, driven by antibiotic treatment. Antibiotic resistance increased as infection progressed toward the establishment of a population constituted by genotypically diversified coexisting sublineages, all of which converged to multidrug resistance. These sublineages emerged by parallel evolution through distinct evolutionary pathways, which affected genes of the same functional categories. Interestingly, ampC and ftsI, encoding the β-lactamase and penicillin-binding protein 3, respectively, were found to be among the most frequently mutated genes. In fact, both genes were targeted by multiple independent mutational events, which led to a wide diversity of coexisting alleles underlying β-lactam resistance. Our findings indicate that hypermutators, apart from boosting antibiotic resistance evolution by simultaneously targeting several genes, favor the emergence of adaptive innovative alleles by clustering beneficial/compensatory mutations in the same gene, hence expanding P. aeruginosa strategies for persistence.

Antibiotic Resistance Evolution Is Contingent on the Quorum-Sensing Response in Pseudomonas aeruginosa

Different works have explored independently the evolution toward antibiotic resistance and the role of eco-adaptive mutations in the adaptation to a new habitat (as the infected host) of bacterial pathogens. However, knowledge about the connection between both processes is still limited. We address this issue by comparing the evolutionary trajectories toward antibiotic resistance of a Pseudomonas aeruginosa lasR defective mutant and its parental wild-type strain, when growing in presence of two ribosome-targeting antibiotics. Quorum-sensing lasR defective mutants are selected in P. aeruginosa populations causing chronic infections. Further, we observed they are also selected in vitro as a first adaptation for growing in culture medium. By using experimental evolution and whole-genome sequencing, we found that the evolutionary trajectories of P. aeruginosa in presence of these antibiotics are different in lasR defective and in wild-type backgrounds, both at the phenotypic and the genotypic levels. Recreation of a set of mutants in both genomic backgrounds (either wild type or lasR defective) allowed us to determine the existence of negative epistatic interactions between lasR and antibiotic resistance determinants. These epistatic interactions could lead to mutual contingency in the evolution of antibiotic resistance when P. aeruginosa colonizes a new habitat in presence of antibiotics. If lasR mutants are selected first, this would constraint antibiotic resistance evolution. Conversely, when resistance mutations (at least those studied in the present work) are selected, lasR mutants may not be selected in presence of antibiotics. These results underlie the importance of contingency and epistatic interactions in modulating antibiotic resistance evolution.

Evolutionary stability of collateral sensitivity to antibiotics in the model pathogen Pseudomonas aeruginosa

Evolution is at the core of the impending antibiotic crisis. Sustainable therapy must thus account for the adaptive potential of pathogens. One option is to exploit evolutionary trade-offs, like collateral sensitivity, where evolved resistance to one antibiotic causes hypersensitivity to another one. To date, the evolutionary stability and thus clinical utility of this trade-off is unclear. We performed a critical experimental test on this key requirement, using evolution experiments with Pseudomonas aeruginosa, and identified three main outcomes: (i) bacteria commonly failed to counter hypersensitivity and went extinct; (ii) hypersensitivity sometimes converted into multidrug resistance; and (iii) resistance gains frequently caused re-sensitization to the previous drug, thereby maintaining the trade-off. Drug order affected the evolutionary outcome, most likely due to variation in the effect size of collateral sensitivity, epistasis among adaptive mutations, and fitness costs. Our finding of robust genetic trade-offs and drug-order effects can guide design of evolution-informed antibiotic therapy.

Historical contingency in the evolution of antibiotic resistance after decades of relaxed selection

Populations often encounter changed environments that remove selection for the maintenance of particular phenotypic traits. The resulting genetic decay of those traits under relaxed selection reduces an organism’s fitness in its prior environment. However, whether and how such decay alters the subsequent evolvability of a population upon restoration of selection for a previously diminished trait is not well understood. We addressed this question using Escherichia coli strains from the long-term evolution experiment (LTEE) that independently evolved for multiple decades in the absence of antibiotics. We first confirmed that these derived strains are typically more sensitive to various antibiotics than their common ancestor. We then subjected the ancestral and derived strains to various concentrations of these drugs to examine their potential to evolve increased resistance. We found that evolvability was idiosyncratic with respect to initial genotype; that is, the derived strains did not generally compensate for their greater susceptibility by “catching up” to the resistance level of the ancestor. Instead, the capacity to evolve increased resistance was constrained in some backgrounds, implying that evolvability depended upon prior mutations in a historically contingent fashion. We further subjected a time series of clones from one LTEE population to tetracycline and determined that an evolutionary constraint arose early in that population, corroborating the role of contingency. In summary, relaxed selection not only can drive populations to increased antibiotic susceptibility, but it can also affect the subsequent evolvability of antibiotic resistance in an unpredictable manner. This conclusion has potential implications for public health, and it underscores the need to consider the genetic context of pathogens when designing drug-treatment strategies.

Strains of Escherichia coli from the long-term evolution experiment (LTEE) have independently evolved for multiple decades and tens of thousands of generations in the absence of antibiotics; this study shows that evolution and diversification of a single bacterial species in an antibiotic-free environment can render the subsequent evolution of resistance unpredictable.

Adaptive walks on high-dimensional fitness landscapes and seascapes with distance-dependent statistics

The dynamics of evolution is intimately shaped by epistasis - interactions between genetic elements which cause the fitness-effect of combinations of mutations to be non-additive. Analyzing evolutionary dynamics that involves large numbers of epistatic mutations is intrinsically difficult. A crucial feature is that the fitness landscape in the vicinity of the current genome depends on the evolutionary history. A key step is thus developing models that enable study of the effects of past evolution on future evolution. In this work, we introduce a broad class of high-dimensional random fitness landscapes for which the correlations between fitnesses of genomes are a general function of genetic distance. Their Gaussian character allows for tractable computational as well as analytic understanding. We study the properties of these landscapes focusing on the simplest evolutionary process: random adaptive (uphill) walks. Conventional measures of "ruggedness" are shown to not much affect such adaptive walks. Instead, the long-distance statistics of epistasis cause all properties to be highly conditional on past evolution, determining the statistics of the local landscape (the distribution of fitness-effects of available mutations and combinations of these), as well as the global geometry of evolutionary trajectories. In order to further explore the effects of conditioning on past evolution, we model the effects of slowly changing environments. At long times, such fitness "seascapes" cause a statistical steady state with highly intermittent evolutionary dynamics: populations undergo bursts of rapid adaptation, interspersed with periods in which adaptive mutations are rare and the population waits for more new directions to be opened up by changes in the environment. Finally, we discuss prospects for studying more complex evolutionary dynamics and on broader classes of high-dimensional landscapes and seascapes.

Learning the pattern of epistasis linking genotype and phenotype in a protein

Understanding the pattern of epistasis-the non-independence of mutations-is critical for relating genotype and phenotype. However, the combinatorial complexity of potential epistatic interactions has severely limited the analysis of this problem. Using new mutational approaches, we report a comprehensive experimental study of all 2¹³ mutants that link two phenotypically distinct variants of the Entacmaea quadricolor fluorescent protein-an opportunity to examine epistasis up to the 13^th order. The data show the existence of many high-order epistatic interactions between mutations, but also reveal extraordinary sparsity, enabling novel experimental and computational strategies for learning the relevant epistasis. We demonstrate that such information, in turn, can be used to accurately predict phenotypes in practical situations where the number of measurements is limited. Finally, we show how the observed epistasis shapes the solution space of single-mutation trajectories between the parental fluorescent proteins, informative about the protein's evolutionary potential. This work provides conceptual and experimental strategies to profoundly characterize epistasis in a protein, relevant to both natural and laboratory evolution.

Experimental evolution of gallium resistance in Escherichia coli

Background and Objectives

Metallic antimicrobial materials are of growing interest due to their potential to control pathogenic and multidrug-resistant bacteria. Yet we do not know if utilizing these materials can lead to genetic adaptations that produce even more dangerous bacterial varieties.

Methodology

Here we utilize experimental evolution to produce strains of Escherichia coli K-12 MG1655 resistant to, the iron analog, gallium nitrate (Ga(NO₃)₃). Whole genome sequencing was utilized to determine genomic changes associated with gallium resistance. Computational modeling was utilized to propose potential molecular mechanisms of resistance.

Results

By day 10 of evolution, increased gallium resistance was evident in populations cultured in medium containing a sublethal concentration of gallium. Furthermore, these populations showed increased resistance to ionic silver and iron (III), but not iron (II) and no increase in traditional antibiotic resistance compared with controls and the ancestral strain. In contrast, the control populations showed increased resistance to rifampicin relative to the gallium-resistant and ancestral population. Genomic analysis identified hard selective sweeps of mutations in several genes in the gallium (III)-resistant lines including: fecA (iron citrate outer membrane transporter), insl1 (IS30 tranposase) one intergenic mutations arsC →/→ yhiS; (arsenate reductase/pseudogene) and in one pseudogene yedN ←; (iapH/yopM family). Two additional significant intergenic polymorphisms were found at frequencies > 0.500 in fepD ←/→ entS (iron-enterobactin transporter subunit/enterobactin exporter, iron-regulated) and yfgF ←/→ yfgG (cyclic-di-GMP phosphodiesterase, anaerobic/uncharacterized protein). The control populations displayed mutations in the rpoB gene, a gene associated with rifampicin resistance.

Conclusions

This study corroborates recent results observed in experiments utilizing pathogenic Pseudomonas strains that also showed that Gram-negative bacteria can rapidly evolve resistance to an atom that mimics an essential micronutrient and shows the pleiotropic consequences associated with this adaptation.

Lay summary

We utilize experimental evolution to produce strains of Escherichia coli K-12 MG1655 resistant to, the iron analog, gallium nitrate (Ga(NO₃)₃). Whole genome sequencing was utilized to determine genomic changes associated with gallium resistance. Computational modeling was utilized to propose potential molecular mechanisms of resistance.

Effect of Environment on the Evolutionary Trajectories and Growth Characteristics of Antibiotic-Resistant Escherichia coli Mutants

The fitness cost to bacteria of acquisition of resistance determinants is critically under-investigated, and the identification and exploitation of these fitness costs may lead to novel therapeutic strategies that prevent the emergence of antimicrobial resistance. Here we used Escherichia coli and amoxicillin–clavulanic acid (AMC) resistance as a model to understand how the artificial environments utilized in studies of bacterial fitness could affect the emergence of resistance and associated fitness costs. Further, we explored the predictive value of this data when strains were grown in the more physiologically relevant environments of urine and urothelial organoids. Resistant E. coli isolates were selected for following 24-h exposure to sub-inhibitory concentrations of AMC in either M9, ISO, or LB, followed by growth on LB agar containing AMC. No resistant colonies emerged following growth in M9, whereas resistant isolates were detected from cultures grown in ISO and LB. We observed both within and between media-type variability in the levels of resistance and fitness of the resistant mutants grown in LB. MICs and fitness of these resistant strains in different media (M9, ISO, LB, human urine, and urothelial organoids) showed considerable variation. Media can therefore have a direct effect on the isolation of mutants that confer resistance to AMC and these mutants can exhibit unpredictable MIC and fitness profiles under different growth conditions. This preliminary study highlights the risks in relying on a single culture protocol as a model system to predict the behavior and treatment response of bacteria in vivo and highlights the importance of developing comprehensive experimental designs to ensure effective translation of diagnostic procedures to successful clinical outcomes.

Rapid decline of bacterial drug-resistance in an antibiotic-free environment through phenotypic reversion

Antibiotic resistance typically induces a fitness cost that shapes the fate of antibiotic-resistant bacterial populations. However, the cost of resistance can be mitigated by compensatory mutations elsewhere in the genome, and therefore the loss of resistance may proceed too slowly to be of practical importance. We present our study on the efficacy and phenotypic impact of compensatory evolution in Escherichia coli strains carrying multiple resistance mutations. We have demonstrated that drug-resistance frequently declines within 480 generations during exposure to an antibiotic-free environment. The extent of resistance loss was found to be generally antibiotic-specific, driven by mutations that reduce both resistance level and fitness costs of antibiotic-resistance mutations. We conclude that phenotypic reversion to the antibiotic-sensitive state can be mediated by the acquisition of additional mutations, while maintaining the original resistance mutations. Our study indicates that restricting antimicrobial usage could be a useful policy, but for certain antibiotics only.

Limited Evolutionary Conservation of the Phenotypic Effects of Antibiotic Resistance Mutations

Multidrug-resistant clinical isolates are common in certain pathogens, but rare in others. This pattern may be due to the fact that mutations shaping resistance have species-specific effects. To investigate this issue, we transferred a range of resistance-conferring mutations and a full resistance gene into Escherichia coli and closely related bacteria. We found that resistance mutations in one bacterial species frequently provide no resistance, in fact even yielding drug hypersensitivity in close relatives. In depth analysis of a key gene involved in aminoglycoside resistance (trkH) indicated that preexisting mutations in other genes-intergenic epistasis-underlie such extreme differences in mutational effects between species. Finally, reconstruction of adaptive landscapes under multiple antibiotic stresses revealed that mutations frequently provide multidrug resistance or elevated drug susceptibility (i.e., collateral sensitivity) only with certain combinations of other resistance mutations. We conclude that resistance and collateral sensitivity are contingent upon the genetic makeup of the bacterial population, and such contingency could shape the long-term fate of resistant bacteria. These results underlie the importance of species-specific treatment strategies.

Pan-Genomic and Polymorphic Driven Prediction of Antibiotic Resistance in Elizabethkingia

The Elizabethkingia are a genetically diverse genus of emerging pathogens that exhibit multidrug resistance to a range of common antibiotics. Two representative species, Elizabethkingia bruuniana and E. meningoseptica, were phenotypically tested to determine minimum inhibitory concentrations (MICs) for five antibiotics. Ultra-long read sequencing with Oxford Nanopore Technologies (ONT) and subsequent de novo assembly produced complete, gapless circular genomes for each strain. Alignment based annotation with Prokka identified 5,480 features in E. bruuniana and 5,203 features in E. meningoseptica, where none of these identified genes or gene combinations corresponded to observed phenotypic resistance values. Pan-genomic analysis, performed with an additional 19 Elizabethkingia strains, identified a core-genome size of 2,658,537 bp, 32 uniquely identifiable intrinsic chromosomal antibiotic resistance core-genes and 77 antibiotic resistance pan-genes. Using core-SNPs and pan-genes in combination with six machine learning (ML) algorithms, binary classification of clindamycin and vancomycin resistance achieved f1 scores of 0.94 and 0.84, respectively. Performance on the more challenging multiclass problem for fusidic acid, rifampin and ciprofloxacin resulted in f1 scores of 0.70, 0.75, and 0.54, respectively. By producing two sets of quality biological predictors, pan-genome genes and core-genome SNPs, from long-read sequence data and applying an ensemble of ML techniques, our results demonstrated that accurate phenotypic inference, at multiple AMR resolutions, can be achieved.

On the evolutionary ecology of multidrug resistance in bacteria

Resistance against different antibiotics appears on the same bacterial strains more often than expected by chance, leading to high frequencies of multidrug resistance. There are multiple explanations for this observation, but these tend to be specific to subsets of antibiotics and/or bacterial species, whereas the trend is pervasive. Here, we consider the question in terms of strain ecology: explaining why resistance to different antibiotics is often seen on the same strain requires an understanding of the competition between strains with different resistance profiles. This work builds on models originally proposed to explain another aspect of strain competition: the stable coexistence of antibiotic sensitivity and resistance observed in a number of bacterial species. We first identify a partial structural similarity in these models: either strain or host population structure stratifies the pathogen population into evolutionarily independent sub-populations and introduces variation in the fitness effect of resistance between these sub-populations, thus creating niches for sensitivity and resistance. We then generalise this unified underlying model to multidrug resistance and show that models with this structure predict high levels of association between resistance to different drugs and high multidrug resistance frequencies. We test predictions from this model in six bacterial datasets and find them to be qualitatively consistent with observed trends. The higher than expected frequencies of multidrug resistance are often interpreted as evidence that these strains are out-competing strains with lower resistance multiplicity. Our work provides an alternative explanation that is compatible with long-term stability in resistance frequencies.

Effects of a previously selected antibiotic resistance on mutations acquired during development of a second resistance in Escherichia coli

Background

The effect of mutations conferring antibiotic resistance can depend on the genetic background. To determine if a previously de novo acquired antibiotic resistance influences the adaptation to a second antibiotic, antibiotic resistance was selected for by exposure to stepwise increasing sublethal levels of amoxicillin, enrofloxacin, kanamycin, or tetracycline. E. coli populations adapted to either a single or two antibiotics sequentially were characterized using whole genome population sequencing and MIC measurements.

Results

In a wild-type background, adaptation to any of the antibiotics resulted in the appearance of well-known mutations, as well as a number of mutated genes not known to be associated with antibiotic resistance. Development of a second resistance in a strain with an earlier acquired resistance to a different antibiotic did not always result in the appearance of all mutations associated with resistance in a wild-type background. In general, a more varied set of mutations was acquired during secondary adaptation. The ability of E. coli to maintain the first resistance during this process depended on the combination of antibiotics used. The maintenance of mutations associated with resistance to the first antibiotic did not always predict the residual MIC for that compound.

Conclusions

In general, the data presented here indicate that adaptation to each antibiotic is unique and independent. The mutational trajectories available in already resistant cells appear more varied than in wild-type cells, indicating that the genetic background of E. coli influences resistance development. The observed mutations cannot always fully explain the resistance pattern observed, indicating a crucial role for adaptation on the gene expression level in de novo acquisition of antibiotic resistance.

Mycoplasmas under experimental antimicrobial selection: The unpredicted contribution of horizontal chromosomal transfer

Horizontal Gene Transfer was long thought to be marginal in Mycoplasma a large group of wall-less bacteria often portrayed as minimal cells because of their reduced genomes (ca. 0.5 to 2.0 Mb) and their limited metabolic pathways. This view was recently challenged by the discovery of conjugative exchanges of large chromosomal fragments that equally affected all parts of the chromosome via an unconventional mechanism, so that the whole mycoplasma genome is potentially mobile. By combining next generation sequencing to classical mating and evolutionary experiments, the current study further explored the contribution and impact of this phenomenon on mycoplasma evolution and adaptation using the fluoroquinolone enrofloxacin (Enro), for selective pressure and the ruminant pathogen Mycoplasma agalactiae, as a model organism. For this purpose, we generated isogenic lineages that displayed different combination of spontaneous mutations in Enro target genes (gyrA, gyrB, parC and parE) in association to gradual level of resistance to Enro. We then tested whether these mutations can be acquired by a susceptible population via conjugative chromosomal transfer knowing that, in our model organism, the 4 target genes are scattered in three distinct and distant loci. Our data show that under antibiotic selective pressure, the time scale of the mutational pathway leading to high-level of Enro resistance can be readily compressed into a single conjugative step, in which several EnroR alleles were transferred from resistant to susceptible mycoplasma cells. In addition to acting as an accelerator for antimicrobial dissemination, mycoplasma chromosomal transfer reshuffled genomes beyond expectations and created a mosaic of resistant sub-populations with unpredicted and unrelated features. Our findings provide insights into the process that may drive evolution and adaptability of several pathogenic Mycoplasma spp. via an unconventional conjugative mechanism.

Insights into the processes that drive the evolution of drug resistance in Mycobacterium tuberculosis

At present, the successful transmission of drug-resistant Mycobacterium tuberculosis, including multidrug-resistant (MDR) and extensively drug-resistant (XDR) strains, in human populations, threatens tuberculosis control worldwide. Differently from many other bacteria, M. tuberculosis drug resistance is acquired mainly through mutations in specific drug resistance-associated genes. The panel of mutations is highly diverse, but depends on the affected gene and M. tuberculosis genetic background. The variety of genetic profiles observed in drug-resistant clinical isolates underlines different evolutionary trajectories towards multiple drug resistance, although some mutation patterns are prominent. This review discusses the intrinsic processes that may influence drug resistance evolution in M. tuberculosis, such as mutation rate, drug resistance-associated mutations, fitness cost, compensatory mutations and epistasis. This knowledge should help to better predict the risk of emergence of highly resistant M. tuberculosis strains and to develop new tools and strategies to limit the development and spread of MDR and XDR strains.

Genetics and roadblocks of drug resistant tuberculosis

Considering the extensive evolutionary history of Mycobacterium tuberculosis, anti-Tuberculosis (TB) drug therapy exerts a recent selective pressure. However, in a microorganism devoid of horizontal gene transfer and with a strictly clonal populational structure such as M. tuberculosis the usual, but not sole, path to overcome drug susceptibility is through de novo mutations on a relatively strict set of genes. The possible allelic diversity that can be associated with drug resistance through several mechanisms such as target alteration or target overexpression, will dictate how these genes can become associated with drug resistance. The success demonstrated by this pathogenic microbe in this latter process and its ability to spread is currently one of the major obstacles to an effective TB elimination. This article reviews the action mechanism of the more important anti-TB drugs, including bedaquiline and delamanid, along with new findings on specific resistance mechanisms. With the development, validation and endorsement of new in vitro molecular tests for drug resistance, knowledge on these resistance mechanisms and microevolutionary dynamics leading to the emergence and fixation of drug resistance mutations within the host is highly important. Additionally, the fitness toll imposed by resistance development is also herein discussed together with known compensatory mechanisms. By elucidating the possible mechanisms that enable one strain to reacquire the original fitness levels, it will be theoretically possible to make more informed decisions and develop novel strategies that can force M. tuberculosis microevolutionary trajectory down through a path of decreasing fitness levels.

Epistatic control of intrinsic resistance by virulence genes in Listeria

Elucidating the relationships between antimicrobial resistance and virulence is key to understanding the evolution and population dynamics of resistant pathogens. Here, we show that the susceptibility of the gram-positive bacterium Listeria monocytogenes to the antibiotic fosfomycin is a complex trait involving interactions between resistance and virulence genes and the environment. We found that a FosX enzyme encoded in the listerial core genome confers intrinsic fosfomycin resistance to both pathogenic and non-pathogenic Listeria spp. However, in the genomic context of the pathogenic L. monocytogenes, FosX-mediated resistance is epistatically suppressed by two members of the PrfA virulence regulon, hpt and prfA, which upon activation by host signals induce increased fosfomycin influx into the bacterial cell. Consequently, in infection conditions, most L. monocytogenes isolates become susceptible to fosfomycin despite possessing a gene that confers high-level resistance to the drug. Our study establishes the molecular basis of an epistatic interaction between virulence and resistance genes controlling bacterial susceptibility to an antibiotic. The reported findings provide the rationale for the introduction of fosfomycin in the treatment of Listeria infections even though these bacteria are intrinsically resistant to the antibiotic in vitro.