Plasticity and Stereotypic Rewiring of the Transcriptome Upon Bacterial Evolution of Antibiotic Resistance.

Bacterial evolution of antibiotic resistance frequently has deleterious side effects on microbial growth, virulence, and susceptibility to other antimicrobial agents. However, it is unclear how these trade-offs could be utilized for manipulating antibiotic resistance in the clinic, not least because the underlying molecular mechanisms are poorly understood. Using laboratory evolution, we demonstrate that clinically relevant resistance mutations in Escherichia coli constitutively rewire a large fraction of the transcriptome in a repeatable and stereotypic manner. Strikingly, lineages adapted to functionally distinct antibiotics and having no resistance mutations in common show a wide range of parallel gene expression changes that alter oxidative stress response, iron homeostasis, and the composition of the bacterial outer membrane and cell surface. These common physiological alterations are associated with changes in cell morphology and enhanced sensitivity to antimicrobial peptides. Finally, the constitutive transcriptomic changes induced by resistance mutations are largely distinct from those induced by antibiotic stresses in the wild type. This indicates a limited role for genetic assimilation of the induced antibiotic stress response during resistance evolution. Our work suggests that diverse resistance mutations converge on similar global transcriptomic states that shape genetic susceptibility to antimicrobial compounds.

Rational design of balanced dual-targeting antibiotics with limited resistance.

Antibiotics that inhibit multiple bacterial targets offer a promising therapeutic strategy against resistance evolution, but developing such antibiotics is challenging. Here we demonstrate that a rational design of balanced multitargeting antibiotics is feasible by using a medicinal chemistry workflow. The resultant lead compounds, ULD1 and ULD2, belonging to a novel chemical class, almost equipotently inhibit bacterial DNA gyrase and topoisomerase IV complexes and interact with multiple evolutionary conserved amino acids in the ATP-binding pockets of their target proteins. ULD1 and ULD2 are excellently potent against a broad range of gram-positive bacteria. Notably, the efficacy of these compounds was tested against a broad panel of multidrug-resistant Staphylococcus aureus clinical strains. Antibiotics with clinical relevance against staphylococcal infections fail to inhibit a significant fraction of these isolates, whereas both ULD1 and ULD2 inhibit all of them (minimum inhibitory concentration [MIC] ≤1 µg/mL). Resistance mutations against these compounds are rare, have limited impact on compound susceptibility, and substantially reduce bacterial growth. Based on their efficacy and lack of toxicity demonstrated in murine infection models, these compounds could translate into new therapies against multidrug-resistant bacterial infections.

Improved bacterial recombineering by parallelized protein discovery.

Exploiting bacteriophage-derived homologous recombination processes has enabled precise, multiplex editing of microbial genomes and the construction of billions of customized genetic variants in a single day. The techniques that enable this, multiplex automated genome engineering (MAGE) and directed evolution with random genomic mutations (DIvERGE), are however, currently limited to a handful of microorganisms for which single-stranded DNA-annealing proteins (SSAPs) that promote efficient recombineering have been identified. Thus, to enable genome-scale engineering in new hosts, efficient SSAPs must first be found. Here we introduce a high-throughput method for SSAP discovery that we call "serial enrichment for efficient recombineering" (SEER). By performing SEER in Escherichia coli to screen hundreds of putative SSAPs, we identify highly active variants PapRecT and CspRecT. CspRecT increases the efficiency of single-locus editing to as high as 50% and improves multiplex editing by 5- to 10-fold in E. coli, while PapRecT enables efficient recombineering in Pseudomonas aeruginosa, a concerning human pathogen. CspRecT and PapRecT are also active in other, clinically and biotechnologically relevant enterobacteria. We envision that the deployment of SEER in new species will pave the way toward pooled interrogation of genotype-to-phenotype relationships in previously intractable bacteria.

Suboptimal Global Transcriptional Response Increases the Harmful Effects of Loss-of-Function Mutations.

The fitness impact of loss-of-function mutations is generally assumed to reflect the loss of specific molecular functions associated with the perturbed gene. Here, we propose that rewiring of the transcriptome upon deleterious gene inactivation is frequently nonspecific and mimics stereotypic responses to external environmental change. Consequently, transcriptional response to gene deletion could be suboptimal and incur an extra fitness cost. Analysis of the transcriptomes of ∼1,500 single-gene deletion Saccharomyces cerevisiae strains supported this scenario. First, most transcriptomic changes are not specific to the deleted gene but are rather triggered by perturbations in functionally diverse genes. Second, gene deletions that alter the expression of dosage-sensitive genes are especially harmful. Third, by elevating the expression level of downregulated genes, we could experimentally mitigate the fitness defect of gene deletions. Our work shows that rewiring of genomic expression upon gene inactivation shapes the harmful effects of mutations.

Pathogen diversity drives the evolution of generalist MHC-II alleles in human populations.

Central players of the adaptive immune system are the groups of proteins encoded in the major histocompatibility complex (MHC), which shape the immune response against pathogens and tolerance to self-peptides. The corresponding genomic region is of particular interest, as it harbors more disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Certain MHC molecules can bind to a much wider range of epitopes than others, but the functional implication of such an elevated epitope-binding repertoire has remained largely unclear. It has been suggested that by recognizing more peptide segments, such promiscuous MHC molecules promote immune response against a broader range of pathogens. If so, the geographical distribution of MHC promiscuity level should be shaped by pathogen diversity. Three lines of evidence support the hypothesis. First, we found that in pathogen-rich geographical regions, humans are more likely to carry highly promiscuous MHC class II DRB1 alleles. Second, the switch between specialist and generalist antigen presentation has occurred repeatedly and in a rapid manner during human evolution. Third, molecular positions that define promiscuity level of MHC class II molecules are especially diverse and are under positive selection in human populations. Taken together, our work indicates that pathogen load maintains generalist adaptive immune recognition, with implications for medical genetics and epidemiology.

Evolthon: A community endeavor to evolve lab evolution.

In experimental evolution, scientists evolve organisms in the lab, typically by challenging them to new environmental conditions. How best to evolve a desired trait? Should the challenge be applied abruptly, gradually, periodically, sporadically? Should one apply chemical mutagenesis, and do strains with high innate mutation rate evolve faster? What are ideal population sizes of evolving populations? There are endless strategies, beyond those that can be exposed by individual labs. We therefore arranged a community challenge, Evolthon, in which students and scientists from different labs were asked to evolve Escherichia coli or Saccharomyces cerevisiae for an abiotic stress-low temperature. About 30 participants from around the world explored diverse environmental and genetic regimes of evolution. After a period of evolution in each lab, all strains of each species were competed with one another. In yeast, the most successful strategies were those that used mating, underscoring the importance of sex in evolution. In bacteria, the fittest strain used a strategy based on exploration of different mutation rates. Different strategies displayed variable levels of performance and stability across additional challenges and conditions. This study therefore uncovers principles of effective experimental evolutionary regimens and might prove useful also for biotechnological developments of new strains and for understanding natural strategies in evolutionary arms races between species. Evolthon constitutes a model for community-based scientific exploration that encourages creativity and cooperation.

Directed evolution of multiple genomic loci allows the prediction of antibiotic resistance.

Antibiotic development is frequently plagued by the rapid emergence of drug resistance. However, assessing the risk of resistance development in the preclinical stage is difficult. Standard laboratory evolution approaches explore only a small fraction of the sequence space and fail to identify exceedingly rare resistance mutations and combinations thereof. Therefore, new rapid and exhaustive methods are needed to accurately assess the potential of resistance evolution and uncover the underlying mutational mechanisms. Here, we introduce directed evolution with random genomic mutations (DIvERGE), a method that allows an up to million-fold increase in mutation rate along the full lengths of multiple predefined loci in a range of bacterial species. In a single day, DIvERGE generated specific mutation combinations, yielding clinically significant resistance against trimethoprim and ciprofloxacin. Many of these mutations have remained previously undetected or provide resistance in a species-specific manner. These results indicate pathogen-specific resistance mechanisms and the necessity of future narrow-spectrum antibacterial treatments. In contrast to prior claims, we detected the rapid emergence of resistance against gepotidacin, a novel antibiotic currently in clinical trials. Based on these properties, DIvERGE could be applicable to identify less resistance-prone antibiotics at an early stage of drug development. Finally, we discuss potential future applications of DIvERGE in synthetic and evolutionary biology.

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.

Mismatch repair hierarchy of Pseudomonas putida revealed by mutagenic ssDNA recombineering of the pyrF gene.

The mismatch repair (MMR) system is one of the key molecular devices that prokaryotic cells have for ensuring fidelity of DNA replication. While the canonical MMR of E. coli involves 3 proteins (encoded by mutS, mutL and mutH), the soil bacterium Pseudomonads putida has only 2 bona fide homologues (mutS and mutL) and the sensitivity of this abridged system to different types of mismatches is unknown. In this background, sensitivity to MMR of this bacterium was inspected through single stranded (ss) DNA recombineering of the pyrF gene (the prokaryotic equivalent to yeast's URA3) with mutagenic oligos representative of every possible mispairing under either wild-type conditions, permanent deletion of mutS or transient loss of mutL activity (brought about by the thermoinducible dominant negative allele mutLE36K ). Analysis of single nucleotide mutations borne by clones resistant to fluoroorotic acid (5FOA, the target of wild type PyrF) pinpointed prohibited and tolerated single-nucleotide replacements and exposed a clear grading of mismatch recognition. The resulting data unequivocally established the hierarchy A:G < C:C < G:A < C:A, A:A, G:G, T:T, T:G, A:C, C:T < G:T, T:C as the one prevalent in Pseudomonas putida. This information is vital for enabling recombineering strategies aimed at single-nucleotide changes in this biotechnologically important species.

Enzyme promiscuity shapes adaptation to novel growth substrates.

Evidence suggests that novel enzyme functions evolved from low-level promiscuous activities in ancestral enzymes. Yet, the evolutionary dynamics and physiological mechanisms of how such side activities contribute to systems-level adaptations are not well characterized. Furthermore, it remains untested whether knowledge of an organism's promiscuous reaction set, or underground metabolism, can aid in forecasting the genetic basis of metabolic adaptations. Here, we employ a computational model of underground metabolism and laboratory evolution experiments to examine the role of enzyme promiscuity in the acquisition and optimization of growth on predicted non-native substrates in Escherichia coli K-12 MG1655. After as few as approximately 20 generations, evolved populations repeatedly acquired the capacity to grow on five predicted non-native substrates-D-lyxose, D-2-deoxyribose, D-arabinose, m-tartrate, and monomethyl succinate. Altered promiscuous activities were shown to be directly involved in establishing high-efficiency pathways. Structural mutations shifted enzyme substrate turnover rates toward the new substrate while retaining a preference for the primary substrate. Finally, genes underlying the phenotypic innovations were accurately predicted by genome-scale model simulations of metabolism with enzyme promiscuity.

Phenotypic heterogeneity promotes adaptive evolution.

Genetically identical cells frequently display substantial heterogeneity in gene expression, cellular morphology and physiology. It has been suggested that by rapidly generating a subpopulation with novel phenotypic traits, phenotypic heterogeneity (or plasticity) accelerates the rate of adaptive evolution in populations facing extreme environmental challenges. This issue is important as cell-to-cell phenotypic heterogeneity may initiate key steps in microbial evolution of drug resistance and cancer progression. Here, we study how stochastic transitions between cellular states influence evolutionary adaptation to a stressful environment in yeast Saccharomyces cerevisiae. We developed inducible synthetic gene circuits that generate varying degrees of expression stochasticity of an antifungal resistance gene. We initiated laboratory evolutionary experiments with genotypes carrying different versions of the genetic circuit by exposing the corresponding populations to gradually increasing antifungal stress. Phenotypic heterogeneity altered the evolutionary dynamics by transforming the adaptive landscape that relates genotype to fitness. Specifically, it enhanced the adaptive value of beneficial mutations through synergism between cell-to-cell variability and genetic variation. Our work demonstrates that phenotypic heterogeneity is an evolving trait when populations face a chronic selection pressure. It shapes evolutionary trajectories at the genomic level and facilitates evolutionary rescue from a deteriorating environmental stress.

Correction: Phenotypic heterogeneity promotes adaptive evolution.

[This corrects the article DOI: 10.1371/journal.pbio.2000644.].

The dawn of evolutionary genome engineering.

Genome engineering strategies--such as genome editing, reduction and shuffling, and de novo genome synthesis--enable the modification of specific genomic locations in a directed and combinatorial manner. These approaches offer an unprecedented opportunity to study central evolutionary issues in which natural genetic variation is limited or biased, which sheds light on the evolutionary forces driving complex and extremely slowly evolving traits; the selective constraints on genome architecture; and the reconstruction of ancestral states of cellular structures and networks.

Rapid Evolution of Reduced Susceptibility against a Balanced Dual-Targeting Antibiotic through Stepping-Stone Mutations.

Multitargeting antibiotics, i.e., single compounds capable of inhibiting two or more bacterial targets, are generally considered to be a promising therapeutic strategy against resistance evolution. The rationale for this theory is that multitargeting antibiotics demand the simultaneous acquisition of multiple mutations at their respective target genes to achieve significant resistance. The theory presumes that individual mutations provide little or no benefit to the bacterial host. Here, we propose that such individual stepping-stone mutations can be prevalent in clinical bacterial isolates, as they provide significant resistance to other antimicrobial agents. To test this possibility, we focused on gepotidacin, an antibiotic candidate that selectively inhibits both bacterial DNA gyrase and topoisomerase IV. In a susceptible organism, Klebsiella pneumoniae, a combination of two specific mutations in these target proteins provide an >2,000-fold reduction in susceptibility, while individually, none of these mutations affect resistance significantly. Alarmingly, strains with decreased susceptibility against gepotidacin are found to be as virulent as the wild-type Klebsiella pneumoniae strain in a murine model. Moreover, numerous pathogenic isolates carry mutations which could promote the evolution of clinically significant reduction of susceptibility against gepotidacin in the future. As might be expected, prolonged exposure to ciprofloxacin, a clinically widely employed gyrase inhibitor, coselected for reduced susceptibility against gepotidacin. We conclude that extensive antibiotic usage could select for mutations that serve as stepping-stones toward resistance against antimicrobial compounds still under development. Our research indicates that even balanced multitargeting antibiotics are prone to resistance evolution.

A highly precise and portable genome engineering method allows comparison of mutational effects across bacterial species.

Currently available tools for multiplex bacterial genome engineering are optimized for a few laboratory model strains, demand extensive prior modification of the host strain, and lead to the accumulation of numerous off-target modifications. Building on prior development of multiplex automated genome engineering (MAGE), our work addresses these problems in a single framework. Using a dominant-negative mutant protein of the methyl-directed mismatch repair (MMR) system, we achieved a transient suppression of DNA repair in Escherichia coli, which is necessary for efficient oligonucleotide integration. By integrating all necessary components into a broad-host vector, we developed a new workflow we term pORTMAGE. It allows efficient modification of multiple loci, without any observable off-target mutagenesis and prior modification of the host genome. Because of the conserved nature of the bacterial MMR system, pORTMAGE simultaneously allows genome editing and mutant library generation in other biotechnologically and clinically relevant bacterial species. Finally, we applied pORTMAGE to study a set of antibiotic resistance-conferring mutations in Salmonella enterica and E. coli. Despite over 100 million y of divergence between the two species, mutational effects remained generally conserved. In sum, a single transformation of a pORTMAGE plasmid allows bacterial species of interest to become an efficient host for genome engineering. These advances pave the way toward biotechnological and therapeutic applications. Finally, pORTMAGE allows systematic comparison of mutational effects and epistasis across a wide range of bacterial species.

Evolution of Robustness to Protein Mistranslation by Accelerated Protein Turnover.

Translational errors occur at high rates, and they influence organism viability and the onset of genetic diseases. To investigate how organisms mitigate the deleterious effects of protein synthesis errors during evolution, a mutant yeast strain was engineered to translate a codon ambiguously (mistranslation). It thereby overloads the protein quality-control pathways and disrupts cellular protein homeostasis. This strain was used to study the capacity of the yeast genome to compensate the deleterious effects of protein mistranslation. Laboratory evolutionary experiments revealed that fitness loss due to mistranslation can rapidly be mitigated. Genomic analysis demonstrated that adaptation was primarily mediated by large-scale chromosomal duplication and deletion events, suggesting that errors during protein synthesis promote the evolution of genome architecture. By altering the dosages of numerous, functionally related proteins simultaneously, these genetic changes introduced large phenotypic leaps that enabled rapid adaptation to mistranslation. Evolution increased the level of tolerance to mistranslation through acceleration of ubiquitin-proteasome-mediated protein degradation and protein synthesis. As a consequence of rapid elimination of erroneous protein products, evolution reduced the extent of toxic protein aggregation in mistranslating cells. However, there was a strong evolutionary trade-off between adaptation to mistranslation and survival upon starvation: the evolved lines showed fitness defects and impaired capacity to degrade mature ribosomes upon nutrient limitation. Moreover, as a response to an enhanced energy demand of accelerated protein turnover, the evolved lines exhibited increased glucose uptake by selective duplication of hexose transporter genes. We conclude that adjustment of proteome homeostasis to mistranslation evolves rapidly, but this adaptation has several side effects on cellular physiology. Our work also indicates that translational fidelity and the ubiquitin-proteasome system are functionally linked to each other and may, therefore, co-evolve in nature.

Indispensability of Horizontally Transferred Genes and Its Impact on Bacterial Genome Streamlining.

Why are certain bacterial genomes so small and compact? The adaptive genome streamlining hypothesis posits that selection acts to reduce genome size because of the metabolic burden of replicating DNA. To reveal the impact of genome streamlining on cellular traits, we reduced the Escherichia coli genome by up to 20% by deleting regions which have been repeatedly subjects of horizontal transfer in nature. Unexpectedly, horizontally transferred genes not only confer utilization of specific nutrients and elevate tolerance to stresses, but also allow efficient usage of resources to build new cells, and hence influence fitness in routine and stressful environments alike. Genome reduction affected fitness not only by gene loss, but also by induction of a general stress response. Finally, we failed to find evidence that the advantage of smaller genomes would be due to a reduced metabolic burden of replicating DNA or a link with smaller cell size. We conclude that as the potential energetic benefit gained by deletion of short genomic segments is vanishingly small compared with the deleterious side effects of these deletions, selection for reduced DNA synthesis costs is unlikely to shape the evolution of small genomes.

The genomic landscape of compensatory evolution.

Adaptive evolution is generally assumed to progress through the accumulation of beneficial mutations. However, as deleterious mutations are common in natural populations, they generate a strong selection pressure to mitigate their detrimental effects through compensatory genetic changes. This process can potentially influence directions of adaptive evolution by enabling evolutionary routes that are otherwise inaccessible. Therefore, the extent to which compensatory mutations shape genomic evolution is of central importance. Here, we studied the capacity of the baker's yeast genome to compensate the complete loss of genes during evolution, and explored the long-term consequences of this process. We initiated laboratory evolutionary experiments with over 180 haploid baker's yeast genotypes, all of which initially displayed slow growth owing to the deletion of a single gene. Compensatory evolution following gene loss was rapid and pervasive: 68% of the genotypes reached near wild-type fitness through accumulation of adaptive mutations elsewhere in the genome. As compensatory mutations have associated fitness costs, genotypes with especially low fitnesses were more likely to be subjects of compensatory evolution. Genomic analysis revealed that as compensatory mutations were generally specific to the functional defect incurred, convergent evolution at the molecular level was extremely rare. Moreover, the majority of the gene expression changes due to gene deletion remained unrestored. Accordingly, compensatory evolution promoted genomic divergence of parallel evolving populations. However, these different evolutionary outcomes are not phenotypically equivalent, as they generated diverse growth phenotypes across environments. Taken together, these results indicate that gene loss initiates adaptive genomic changes that rapidly restores fitness, but this process has substantial pleiotropic effects on cellular physiology and evolvability upon environmental change. Our work also implies that gene content variation across species could be partly due to the action of compensatory evolution rather than the passive loss of genes.