Pervasive compensatory adaptation in Escherichia coli

Within-plant genetic drift to control virus adaptation to host resistance genes

Manipulating evolutionary forces imposed by hosts on pathogens like genetic drift and selection could avoid the emergence of virulent pathogens. For instance, increasing genetic drift could decrease the risk of pathogen adaptation through the random fixation of deleterious mutations or the elimination of favorable ones in the pathogen population. However, no experimental proof of this approach is available for a plant-pathogen system. We studied the impact of pepper (Capsicum annuum) lines carrying the same major resistance gene but contrasted genetic backgrounds on the evolution of Potato virus Y (PVY). The pepper lines were chosen for the contrasted levels of genetic drift (inversely related to Ne, the effective population size) they exert on PVY populations, as well as for their contrasted resistance efficiency (inversely related to the initial replicative fitness, Wi, of PVY in these lines). Experimental evolution was performed by serially passaging 64 PVY populations every month on six contrasted pepper lines during seven months. These PVY populations exhibited highly divergent evolutionary trajectories, ranging from viral extinctions to replicative fitness gains. The sequencing of the PVY VPg cistron, where adaptive mutations are likely to occur, allowed linking these replicative fitness gains to parallel adaptive nonsynonymous mutations. Evolutionary trajectories were well explained by the genetic drift imposed by the host. More specifically, Ne, Wi and their synergistic interaction played a major role in the fate of PVY populations. When Ne was low (i.e. strong genetic drift), the final PVY replicative fitness remained close to the initial replicative fitness, whereas when Ne was high (i.e. low genetic drift), the final PVY replicative fitness was high independently of the replicative fitness of the initially inoculated virus. We show that combining a high resistance efficiency (low Wi) and a strong genetic drift (low Ne) is the best solution to increase resistance durability, that is, to avoid virus adaptation on the long term.

Horizontal gene transfer facilitates the molecular reverse-evolution of antibiotic sensitivity in experimental populations of H. pylori

Reversing the evolution of traits harmful to humans, such as antimicrobial resistance, is a key ambition of applied evolutionary biology. A major impediment to reverse evolution is the relatively low spontaneous mutation rates that revert evolved genotypes back to their ancestral state. However, the repeated re-introduction of ancestral alleles by horizontal gene transfer (HGT) could make reverse evolution likely. Here we evolve populations of an antibiotic-resistant strain of Helicobacter pylori in growth conditions without antibiotics while introducing an ancestral antibiotic-sensitive allele by HGT. We evaluate reverse evolution using DNA sequencing and find that HGT facilitates the molecular reverse evolution of the antibiotic resistance allele, and that selection for high rates of HGT drives the evolution of increased HGT rates in low-HGT treatment populations. Finally, we use a theoretical model and carry out simulations to infer how the fitness costs of antibiotic resistance, rates of HGT and effects of genetic drift interact to determine the probability and predictability of reverse evolution.

Rapid evolutionary repair by secondary perturbation of a primary disrupted transcriptional network

The discrete steps of transcriptional rewiring have been proposed to occur neutrally to ensure steady gene expression under stabilizing selection. A conflict-free switch of a regulon between regulators may require an immediate compensatory evolution to minimize deleterious effects. Here, we perform an evolutionary repair experiment on the Lachancea kluyveri yeast sef1Δ mutant using a suppressor development strategy. Complete loss of SEF1 forces cells to initiate a compensatory process for the pleiotropic defects arising from misexpression of TCA cycle genes. Using different selective conditions, we identify two adaptive loss-of-function mutations of IRA1 and AZF1. Subsequent analyses show that Azf1 is a weak transcriptional activator regulated by the Ras1-PKA pathway. Azf1 loss-of-function triggers extensive gene expression changes responsible for compensatory, beneficial, and trade-off phenotypes. The trade-offs can be alleviated by higher cell density. Our results not only indicate that secondary transcriptional perturbation provides rapid and adaptive mechanisms potentially stabilizing the initial stage of transcriptional rewiring but also suggest how genetic polymorphisms of pleiotropic mutations could be maintained in the population.

Revisiting the Design of the Long-Term Evolution Experiment with Escherichia coli

The long-term evolution experiment (LTEE) with Escherichia coli began in 1988 and it continues to this day, with its 12 populations having recently reached 75,000 generations of evolution in a simple, well-controlled environment. The LTEE was designed to explore open-ended questions about the dynamics and repeatability of phenotypic and genetic evolution. Here I discuss various aspects of the LTEE's experimental design that have enabled its stability and success, including the choices of the culture regime, growth medium, ancestral strain, and statistical replication. I also discuss some of the challenges associated with a long-running project, such as handling procedural errors (e.g., cross-contamination) and managing the expanding collection of frozen samples. The simplicity of the experimental design and procedures have supported the long-term stability of the LTEE. That stability-along with the inherent creativity of the evolutionary process and the emergence of new genomic technologies-provides a platform that has allowed talented students and collaborators to pose questions, collect data, and make discoveries that go far beyond anything I could have imagined at the start of the LTEE.

Epistasis and evolution: recent advances and an outlook for prediction

As organisms evolve, the effects of mutations change as a result of epistatic interactions with other mutations accumulated along the line of descent. This can lead to shifts in adaptability or robustness that ultimately shape subsequent evolution. Here, we review recent advances in measuring, modeling, and predicting epistasis along evolutionary trajectories, both in microbial cells and single proteins. We focus on simple patterns of global epistasis that emerge in this data, in which the effects of mutations can be predicted by a small number of variables. The emergence of these patterns offers promise for efforts to model epistasis and predict evolution.

Pneumococcal empyema: Resistance patterns, fitness cost and serotype distribution

BACKGROUND

Streptococcus pneumoniae is a recognized etiology of invasive infections including parapneumonic empyema, and its resistance to antibiotics is evolving worldwide, raising concerns of encountering untreatable strains. This study measured the serotype distribution, antimicrobial susceptibility and biological cost incurred by resistance of pneumococci from pleural samples.

METHODS

The serotype profiles, susceptibility results and growth rates were phenotypically determined for a panel of clinical strains of S. pneumoniae from cases of empyema between 2011 and 2019.

RESULTS

Of 24 empyema cases, the isolated strains belonged to seven serotypes in the following descending order; 19A, 11A/D, 19F, 3, 7F, 1/6B while two strains remained non-typable. Penicillin susceptibility was shown in <80% of the isolates, while parenteral cephalosporins (cefuroxime and ceftriaxone) demonstrated activity in 83.3 and 95.8% respectively. High resistance frequency was noted for macrolides and sulfonamides, but the strains were uniformly sensitive to respiratory fluroquinolones, vancomycin and linezolid. The macrolide-resistant strain exhibited a high growth rate, suggesting a possible beneficial effect. Phenotypes with mono-resistance to sulfonamides and clindamycin were equally fit as the susceptible counterpart strains. Resistance to multiple antimicrobial agents resulted in a high degree of fitness deficit, while other resistant phenotypes were less fit.

CONCLUSIONS

The pneumococcal conjugate vaccine PCV13 serotypes still circulate in the community. The data indicate that resistance to certain antimicrobials incurs an apparent fitness cost in pneumococci which may limit the dissemination of such strains while low fitness cost, seen in case of resistance to macrolides, may contribute to the spread of resistant clones.

Non-genetic and genetic rewiring underlie adaptation to hypomorphic alleles of an essential gene

Adaptive evolution to cellular stress is a process implicated in a wide range of biological and clinical phenomena. Two major routes of adaptation have been identified: non-genetic changes, which allow expression of different phenotypes in novel environments, and genetic variation achieved by selection of fitter phenotypes. While these processes are broadly accepted, their temporal and epistatic features in the context of cellular evolution and emerging drug resistance are contentious. In this manuscript, we generated hypomorphic alleles of the essential nuclear pore complex (NPC) gene NUP58. By dissecting early and long-term mechanisms of adaptation in independent clones, we observed that early physiological adaptation correlated with transcriptome rewiring and upregulation of genes known to interact with the NPC; long-term adaptation and fitness recovery instead occurred via focal amplification of NUP58 and restoration of mutant protein expression. These data support the concept that early phenotypic plasticity allows later acquisition of genetic adaptations to a specific impairment. We propose this approach as a genetic model to mimic targeted drug therapy in human cells and to dissect mechanisms of adaptation.

The cost of evolved constitutive lac gene expression is usually, but not always, maintained during evolution of generalist populations

Beneficial mutations can become costly following an environmental change. Compensatory mutations can relieve these costs, while not affecting the selected function, so that the benefits are retained if the environment shifts back to be similar to the one in which the beneficial mutation was originally selected. Compensatory mutations have been extensively studied in the context of antibiotic resistance, responses to specific genetic perturbations, and in the determination of interacting gene network components. Few studies have focused on the role of compensatory mutations during more general adaptation, especially as the result of selection in fluctuating environments where adaptations to different environment components may often involve trade-offs. We examine whether costs of a mutation in lacI, which deregulated the expression of the lac operon in evolving populations of Escherichia coli bacteria, were compensated. This mutation occurred in multiple replicate populations selected in environments that fluctuated between growth on lactose, where the mutation was beneficial, and on glucose, where it was deleterious. We found that compensation for the cost of the lacI mutation was rare, but, when it did occur, it did not negatively affect the selected benefit. Compensation was not more likely to occur in a particular evolution environment. Compensation has the potential to remove pleiotropic costs of adaptation, but its rarity indicates that the circumstances to bring about the phenomenon may be peculiar to each individual or impeded by other selected mutations.

Changes in the distribution of fitness effects and adaptive mutational spectra following a single first step towards adaptation

Historical contingency and diminishing returns epistasis have been typically studied for relatively divergent genotypes and/or over long evolutionary timescales. Here, we use Saccharomyces cerevisiae to study the extent of diminishing returns and the changes in the adaptive mutational spectra following a single first adaptive mutational step. We further evolve three clones that arose under identical conditions from a common ancestor. We follow their evolutionary dynamics by lineage tracking and determine adaptive outcomes using fitness assays and whole genome sequencing. We find that diminishing returns manifests as smaller fitness gains during the 2^nd step of adaptation compared to the 1^st step, mainly due to a compressed distribution of fitness effects. We also find that the beneficial mutational spectra for the 2^nd adaptive step are contingent on the 1^st step, as we see both shared and diverging adaptive strategies. Finally, we find that adaptive loss-of-function mutations, such as nonsense and frameshift mutations, are less common in the second step of adaptation than in the first step.

Analyses of both natural and experimental evolution suggest that adaptation depends on the evolutionary past and adaptive potential decreases over time. Here, by tracking yeast adaptation with DNA barcoding, the authors show that such evolutionary phenomena can be observed even after a single adaptive step.

Freshwater environments as reservoirs of antibiotic resistant bacteria and their role in the dissemination of antibiotic resistance genes

Freshwater environments are susceptible to possible contamination by residual antibiotics that are released through different sources, such as agricultural runoffs, sewage discharges and leaching from nearby farms. Freshwater environment can thus become reservoirs where an antibiotic impact microorganisms, and is an important public health concern. Degradation and dilution processes are fundamental for predicting the actual risk of antibiotic resistance dissemination from freshwater reservoirs. This study reviews major approaches for detecting and quantifying antibiotic resistance bacteria (ARBs) and genes (ARGs) in freshwater and their prevalence in these environments. Finally, the role of dilution, degradation, transmission and the persistence and fate of ARB/ARG in these environments are also reviewed. Culture-based single strain approaches and molecular techniques that include polymerase chain reaction (PCR), quantitative polymerase chain reaction (qPCR) and metagenomics are techniques for quantifying ARB and ARGs in freshwater environments. The level of ARBs is extremely high in most of the river systems (up to 98% of the total detected bacteria), followed by lakes (up to 77% of the total detected bacteria), compared to dam, pond, and spring (<1%). Of most concern is the occurrence of extended-spectrum β-lactamase producing Enterobacteriaceae, methicillin resistant Staphylococcus aureus (MRSA) and vancomycin resistant Enterococcus (VRE), which cause highly epidemic infections. Dilution and natural degradation do not completely eradicate ARBs and ARGs in the freshwater environment. Even if the ARBs in freshwater are effectively inactivated by sunlight, their ARG-containing DNA can still be intact and capable of transferring resistance to non-resistant strains. Antibiotic resistance persists and is preserved in freshwater bodies polluted with high concentrations of antibiotics. Direct transmission of indigenous freshwater ARBs to humans as well as their transitory insertion in the microbiota can occur. These findings are disturbing especially for people that rely on freshwater resources for drinking, crop irrigation, and food in form of fish.

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.

Evolutionary pathways to antibiotic resistance are dependent upon environmental structure and bacterial lifestyle

Bacterial populations vary in their stress tolerance and population structure depending upon whether growth occurs in well-mixed or structured environments. We hypothesized that evolution in biofilms would generate greater genetic diversity than well-mixed environments and lead to different pathways of antibiotic resistance. We used experimental evolution and whole genome sequencing to test how the biofilm lifestyle influenced the rate, genetic mechanisms, and pleiotropic effects of resistance to ciprofloxacin in Acinetobacter baumannii populations. Both evolutionary dynamics and the identities of mutations differed between lifestyle. Planktonic populations experienced selective sweeps of mutations including the primary topoisomerase drug targets, whereas biofilm-adapted populations acquired mutations in regulators of efflux pumps. An overall trade-off between fitness and resistance level emerged, wherein biofilm-adapted clones were less resistant than planktonic but more fit in the absence of drug. However, biofilm populations developed collateral sensitivity to cephalosporins, demonstrating the clinical relevance of lifestyle on the evolution of resistance.

Patterns and Mechanisms of Diminishing Returns from Beneficial Mutations

Diminishing returns epistasis causes the benefit of the same advantageous mutation smaller in fitter genotypes and is frequently observed in experimental evolution. However, its occurrence in other contexts, environment dependence, and mechanistic basis are unclear. Here, we address these questions using 1,005 sequenced segregants generated from a yeast cross. Under each of 47 examined environments, 66-92% of tested polymorphisms exhibit diminishing returns epistasis. Surprisingly, improving environment quality also reduces the benefits of advantageous mutations even when fitness is controlled for, indicating the necessity to revise the global epistasis hypothesis. We propose that diminishing returns originates from the modular organization of life where the contribution of each functional module to fitness is determined jointly by the genotype and environment and has an upper limit, and demonstrate that our model predictions match empirical observations. These findings broaden the concept of diminishing returns epistasis, reveal its generality and potential cause, and have important evolutionary implications.

Novel approach to quantitative spatial gene expression uncovers genetic stochasticity in the developing Drosophila eye

Robustness in development allows for the accumulation of genetically based variation in expression. However, this variation is usually examined in response to large perturbations, and examination of this variation has been limited to being spatial, or quantitative, but because of technical restrictions not both. Here we bridge these gaps by investigating replicated quantitative spatial gene expression using rigorous statistical models, in different genotypes, sexes, and species (Drosophila melanogaster and D. simulans). Using this type of quantitative approach with molecular developmental data allows for comparison among conditions, such as different genetic backgrounds. We apply this approach to the morphogenetic furrow, a wave of differentiation that patterns the developing eye disc. Within the morphogenetic furrow, we focus on four genes, hairy, atonal, hedgehog, and Delta. Hybridization chain reaction quantitatively measures spatial gene expression, co-staining for all four genes simultaneously. We find considerable variation in the spatial expression pattern of these genes in the eye between species, genotypes, and sexes. We also find that there has been evolution of the regulatory relationship between these genes, and that their spatial interrelationships have evolved between species. This variation has no phenotypic effect, and could be buffered by network thresholds or compensation from other genes. Both of these mechanisms could potentially be contributing to long term developmental systems drift.

Engineered Functional Redundancy Relaxes Selective Constraints upon Endogenous Genes in Viral RNA Genomes

Functional redundancy, understood as the functional overlap of different genes, is a double-edge sword. At the one side, it is thought to serve as a robustness mechanism that buffers the deleterious effect of mutations hitting one of the redundant copies, thus resulting in pseudogenization. At the other side, it is considered as a source of genetic and functional innovation. In any case, genetically redundant genes are expected to show an acceleration in the rate of molecular evolution. Here, we tackle the role of functional redundancy in viral RNA genomes. To this end, we have evaluated the rates of compensatory evolution for deleterious mutations affecting an essential function, the suppression of RNA silencing plant defense, of tobacco etch potyvirus (TEV). TEV genotypes containing deleterious mutations in presence/absence of engineered functional redundancy were evolved and the pattern of fitness and pathogenicity recovery evaluated. Genetically redundant genotypes suffered less from the effect of deleterious mutations and showed relatively minor changes in fitness and pathogenicity. By contrast, nongenetically redundant genotypes had very low fitness and pathogenicity at the beginning of the evolution experiment that were fully recovered by the end. At the molecular level, the outcome depended on the combination of the actual mutations being compensated and the presence/absence of functional redundancy. Reversions to wild-type alleles were the norm in the nonredundant genotypes while redundant ones either did not fix any mutation at all or showed a higher nonsynonymous mutational load.

Evolution of gene knockout strains of E. coli reveal regulatory architectures governed by metabolism

Biological regulatory network architectures are multi-scale in their function and can adaptively acquire new functions. Gene knockout (KO) experiments provide an established experimental approach not just for studying gene function, but also for unraveling regulatory networks in which a gene and its gene product are involved. Here we study the regulatory architecture of Escherichia coli K-12 MG1655 by applying adaptive laboratory evolution (ALE) to metabolic gene KO strains. Multi-omic analysis reveal a common overall schema describing the process of adaptation whereby perturbations in metabolite concentrations lead regulatory networks to produce suboptimal states, whose function is subsequently altered and re-optimized through acquisition of mutations during ALE. These results indicate that metabolite levels, through metabolite-transcription factor interactions, have a dominant role in determining the function of a multi-scale regulatory architecture that has been molded by evolution.

The function of metabolic genes in the context of regulatory networks is not well understood. Here, the authors investigate the adaptive responses of E. coli after knockout of metabolic genes and highlight the influence of metabolite levels in the evolution of regulatory function.

Growth Adaptation of gnd and sdhCB Escherichia coli Deletion Strains Diverges From a Similar Initial Perturbation of the Transcriptome

Adaptive laboratory evolution (ALE) has emerged as a new approach with which to pursue fundamental biological inquiries and, in particular, new insights into the systemic function of a gene product. Two E. coli knockout strains were constructed: one that blocked the Pentose Phosphate Pathway (gnd KO) and one that decoupled the TCA cycle from electron transport (sdhCDAB KO). Despite major perturbations in central metabolism, minimal growth rate changes were found in the two knockout strains. More surprisingly, many similarities were found in their initial transcriptomic states that could be traced to similarly perturbed metabolites despite the differences in the network location of the gene perturbations and concomitant re-routing of pathway fluxes around these perturbations. However, following ALE, distinct metabolomic and transcriptomic states were realized. These included divergent flux and gene expression profiles in the gnd and sdhCDAB KOs to overcome imbalances in NADPH production and nitrogen/sulfur assimilation, respectively, that were not obvious limitations of growth in the unevolved knockouts. Therefore, this work demonstrates that ALE provides a productive approach to reveal novel insights of gene function at a systems level that cannot be found by observing the fresh knockout alone.

Molecular basis of hemoglobin adaptation in the high-flying bar-headed goose

During the adaptive evolution of a particular trait, some selectively fixed mutations may be directly causative and others may be purely compensatory. The relative contribution of these two classes of mutation to adaptive phenotypic evolution depends on the form and prevalence of mutational pleiotropy. To investigate the nature of adaptive substitutions and their pleiotropic effects, we used a protein engineering approach to characterize the molecular basis of hemoglobin (Hb) adaptation in the high-flying bar-headed goose (Anser indicus), a hypoxia-tolerant species renowned for its trans-Himalayan migratory flights. To test the effects of observed substitutions on evolutionarily relevant genetic backgrounds, we synthesized all possible genotypic intermediates in the line of descent connecting the wildtype bar-headed goose genotype with the most recent common ancestor of bar-headed goose and its lowland relatives. Site-directed mutagenesis experiments revealed one major-effect mutation that significantly increased Hb-O2 affinity on all possible genetic backgrounds. Two other mutations exhibited smaller average effect sizes and less additivity across backgrounds. One of the latter mutations produced a concomitant increase in the autoxidation rate, a deleterious side-effect that was fully compensated by a second-site mutation at a spatially proximal residue. The experiments revealed three key insights: (i) subtle, localized structural changes can produce large functional effects; (ii) relative effect sizes of function-altering mutations may depend on the sequential order in which they occur; and (iii) compensation of deleterious pleiotropic effects may play an important role in the adaptive evolution of protein function.

Evolutionary consequences of multidriver environmental change in an aquatic primary producer

Climate change is altering aquatic environments in a complex way, and simultaneous shifts in many properties will drive evolutionary responses in primary producers at the base of both freshwater and marine ecosystems. So far, evolutionary studies have shown how changes in environmental drivers, either alone or in pairs, affect the evolution of growth and other traits in primary producers. Here, we evolve a primary producer in 96 unique environments with different combinations of between one and eight environmental drivers to understand how evolutionary responses to environmental change depend on the identity and number of drivers. Even in multidriver environments, only a few dominant drivers explain most of the evolutionary changes in population growth rates. Most populations converge on the same growth rate by the end of the evolution experiment. However, populations adapt more when these dominant drivers occur in the presence of other drivers. This is due to an increase in the intensity of selection in environments with more drivers, which are more likely to include dominant drivers. Concurrently, many of the trait changes that occur during the initial short-term response to both single and multidriver environmental change revert after about 450 generations of evolution. In future aquatic environments, populations will encounter differing combinations of drivers and intensities of selection, which will alter the adaptive potential of primary producers. Accurately gauging the intensity of selection on key primary producers will help in predicting population size and trait evolution at the base of aquatic food webs.

Experimental evolution and the dynamics of adaptation and genome evolution in microbial populations

Evolution is an on-going process, and it can be studied experimentally in organisms with rapid generations. My team has maintained 12 populations of Escherichia coli in a simple laboratory environment for >25 years and 60 000 generations. We have quantified the dynamics of adaptation by natural selection, seen some of the populations diverge into stably coexisting ecotypes, described changes in the bacteria's mutation rate, observed the new ability to exploit a previously untapped carbon source, characterized the dynamics of genome evolution and used parallel evolution to identify the genetic targets of selection. I discuss what the future might hold for this particular experiment, briefly highlight some other microbial evolution experiments and suggest how the fields of experimental evolution and microbial ecology might intersect going forward.

Multidrug-resistant bacteria compensate for the epistasis between resistances

Mutations conferring resistance to antibiotics are typically costly in the absence of the drug, but bacteria can reduce this cost by acquiring compensatory mutations. Thus, the rate of acquisition of compensatory mutations and their effects are key for the maintenance and dissemination of antibiotic resistances. While compensation for single resistances has been extensively studied, compensatory evolution of multiresistant bacteria remains unexplored. Importantly, since resistance mutations often interact epistatically, compensation of multiresistant bacteria may significantly differ from that of single-resistant strains. We used experimental evolution, next-generation sequencing, in silico simulations, and genome editing to compare the compensatory process of a streptomycin and rifampicin double-resistant Escherichia coli with those of single-resistant clones. We demonstrate that low-fitness double-resistant bacteria compensate faster than single-resistant strains due to the acquisition of compensatory mutations with larger effects. Strikingly, we identified mutations that only compensate for double resistance, being neutral or deleterious in sensitive or single-resistant backgrounds. Moreover, we show that their beneficial effects strongly decrease or disappear in conditions where the epistatic interaction between resistance alleles is absent, demonstrating that these mutations compensate for the epistasis. In summary, our data indicate that epistatic interactions between antibiotic resistances, leading to large fitness costs, possibly open alternative paths for rapid compensatory evolution, thereby potentially stabilizing costly multiple resistances in bacterial populations.

Antibiotics target essential cellular functions, such as translation or cell wall biogenesis, and bacteria can become resistant to antibiotics by acquiring mutations in genes encoding those functions. This causes most drug-resistance mutations to be detrimental in the absence of the drug. However, bacteria can reduce this handicap by acquiring additional mutations, known as compensatory mutations. Compensatory evolution is crucial for the maintenance and dissemination of antibiotic resistances in bacterial populations. While compensation for single resistances has been extensively studied, compensatory evolution of multidrug-resistant bacteria remains unexplored. Importantly, interactions between resistance mutations are frequent, and this may cause compensation of multidrug-resistant bacteria to differ significantly from that of single-resistant strains. By comparing compensation of single- and double-drug–resistant E. coli, we found that double-drug–resistant bacteria compensate faster than single-drug–resistant strains. This is due to the acquisition of compensatory mutations with larger effects and possibly driven by the large fitness cost of double-drug resistance. Strikingly, we identified mutations that compensate specifically for the interaction between drug resistances, since they are beneficial only for double-drug–resistant bacteria and in conditions in which the interaction between resistances occurs. In summary, our data indicate that certain interactions between antibiotic-resistance mutations can open alternative paths for rapid compensatory evolution, thereby potentially stabilizing multiple drug resistances in bacterial populations.

Epistasis and the Evolution of Antimicrobial Resistance

The fitness effects of a mutation can depend, sometimes dramatically, on genetic background; this phenomenon is often referred to as “epistasis.” Epistasis can have important practical consequences in the context of antimicrobial resistance (AMR). For example, genetic background plays an important role in determining the costs of resistance, and hence in whether resistance will persist in the absence of antibiotic pressure. Furthermore, interactions between resistance mutations can have important implications for the evolution of multi-drug resistance. I argue that there is a need to better characterize the extent and nature of epistasis for mutations and horizontally transferred elements conferring AMR, particularly in clinical contexts. Furthermore, I suggest that epistasis should be an important consideration in attempts to slow or limit the evolution of AMR.

Sustained fitness gains and variability in fitness trajectories in the long-term evolution experiment with Escherichia coli

Many populations live in environments subject to frequent biotic and abiotic changes. Nonetheless, it is interesting to ask whether an evolving population's mean fitness can increase indefinitely, and potentially without any limit, even in a constant environment. A recent study showed that fitness trajectories of Escherichia coli populations over 50 000 generations were better described by a power-law model than by a hyperbolic model. According to the power-law model, the rate of fitness gain declines over time but fitness has no upper limit, whereas the hyperbolic model implies a hard limit. Here, we examine whether the previously estimated power-law model predicts the fitness trajectory for an additional 10 000 generations. To that end, we conducted more than 1100 new competitive fitness assays. Consistent with the previous study, the power-law model fits the new data better than the hyperbolic model. We also analysed the variability in fitness among populations, finding subtle, but significant, heterogeneity in mean fitness. Some, but not all, of this variation reflects differences in mutation rate that evolved over time. Taken together, our results imply that both adaptation and divergence can continue indefinitely—or at least for a long time—even in a constant environment.

Molecular and Proteomic Analysis of Levofloxacin and Metronidazole Resistant Helicobacter pylori

Antibiotic resistance in bacteria incurs fitness cost, but compensatory mechanisms may ameliorate the cost and sustain the resistance even under antibiotics-free conditions. The aim of this study was to determine compensatory mechanisms of antibiotic resistance in H. pylori. Five strains of levofloxacin-sensitive H. pylori were induced in vitro to develop resistance. In addition, four pairs of metronidazole-sensitive and -resistant H. pylori strains were isolated from patients carrying dual H. pylori populations that consist of both sensitive and resistant phenotypes. Growth rate, virulence and biofilm-forming ability of the sensitive and resistant strains were compared to determine effects of compensatory response. Proteome profiles of paired sensitive and resistant strains were analyzed by liquid chromatography/mass spectrophotometry (LC/MS). Although there were no significant differences in growth rate between sensitive and resistant pairs, bacterial virulence (in terms of abilities to induce apoptosis and form biofilm) differs from pair to pair. These findings demonstrate the complex and strain-specific phenotypic changes in compensation for antibiotics resistance. Compensation for in vitro induced levofloxacin resistance involving mutations of gyrA and gyrB was functionally random. Furthermore, higher protein translation and non-functional protein degradation capabilities in naturally-occuring dual population metronidazole sensitive-resistant strains may be a possible alternative mechanism underlying resistance to metronidazole without mutations in rdxA and frxA. This may explain the lack of mutations in target genes in ~10% of metronidazole resistant strains.

How could fully scaled carps appear in natural waters in Madagascar?

The capacity of organisms to rapidly evolve in response to environmental changes is a key feature of evolution, and studying mutation compensation is a way to evaluate whether alternative routes of evolution are possible or not. Common carps (Cyprinus carpio) carrying a homozygous loss-of-function mutation for the scale cover gene fgfr1a1, causing the 'mirror' reduced scale cover, were introduced in Madagascar a century ago. Here we show that carps in Malagasy natural waters are now predominantly covered with scales, though they still all carry the homozygous mutation. We also reveal that the number of scales in mutated carps is under strong polygenic genetic control, with a heritability of 0.49. As a whole, our results suggest that carps submitted to natural selection could evolve a wild-type-like scale cover in less than 40 generations from standing polygenic genetic variation, confirming similar findings mainly retrieved from model organisms.

Adaptation to Parasites and Costs of Parasite Resistance in Mutator and Nonmutator Bacteria

Parasitism creates selection for resistance mechanisms in host populations and is hypothesized to promote increased host evolvability. However, the influence of these traits on host evolution when parasites are no longer present is unclear. We used experimental evolution and whole-genome sequencing of Escherichia coli to determine the effects of past and present exposure to parasitic viruses (phages) on the spread of mutator alleles, resistance, and bacterial competitive fitness. We found that mutator alleles spread rapidly during adaptation to any of four different phage species, and this pattern was even more pronounced with multiple phages present simultaneously. However, hypermutability did not detectably accelerate adaptation in the absence of phages and recovery of fitness costs associated with resistance. Several lineages evolved phage resistance through elevated mucoidy, and during subsequent evolution in phage-free conditions they rapidly reverted to nonmucoid, phage-susceptible phenotypes. Genome sequencing revealed that this phenotypic reversion was achieved by additional genetic changes rather than by genotypic reversion of the initial resistance mutations. Insertion sequence (IS) elements played a key role in both the acquisition of resistance and adaptation in the absence of parasites; unlike single nucleotide polymorphisms, IS insertions were not more frequent in mutator lineages. Our results provide a genetic explanation for rapid reversion of mucoidy, a phenotype observed in other bacterial species including human pathogens. Moreover, this demonstrates that the types of genetic change underlying adaptation to fitness costs, and consequently the impact of evolvability mechanisms such as increased point-mutation rates, depend critically on the mechanism of resistance.

Adaptation of Escherichia coli to glucose promotes evolvability in lactose

The selective history of a population can influence its subsequent evolution, an effect known as historical contingency. We previously observed that five of six replicate populations that were evolved in a glucose-limited environment for 2000 generations, then switched to lactose for 1000 generations, had higher fitness increases in lactose than populations started directly from the ancestor. To test if selection in glucose systematically increased lactose evolvability, we started 12 replay populations--six from a population subsample and six from a single randomly selected clone--from each of the six glucose-evolved founder populations. These replay populations and 18 ancestral populations were evolved for 1000 generations in a lactose-limited environment. We found that replay populations were initially slightly less fit in lactose than the ancestor, but were more evolvable, in that they increased in fitness at a faster rate and to higher levels. This result indicates that evolution in the glucose environment resulted in genetic changes that increased the potential of genotypes to adapt to lactose. Genome sequencing identified four genes--iclR, nadR, spoT, and rbs--that were mutated in most glucose-evolved clones and are candidates for mediating increased evolvability. Our results demonstrate that short-term selective costs during selection in one environment can lead to changes in evolvability that confer longer term benefits.

Elucidating the molecular architecture of adaptation via evolve and resequence experiments

Evolve and resequence (E&R) experiments use experimental evolution to adapt populations to a novel environment, then next-generation sequencing to analyse genetic changes. They enable molecular evolution to be monitored in real time on a genome-wide scale. Here, we review the field of E&R experiments across diverse systems, ranging from simple non-living RNA to bacteria, yeast and the complex multicellular organism Drosophila melanogaster. We explore how different evolutionary outcomes in these systems are largely consistent with common population genetics principles. Differences in outcomes across systems are largely explained by different starting population sizes, levels of pre-existing genetic variation, recombination rates and adaptive landscapes. We highlight emerging themes and inconsistencies that future experiments must address.

Constraints on adaptation of Escherichia coli to mixed-resource environments increase over time

Can a population evolved in two resources reach the same fitness in both as specialist populations evolved in each of the individual resources? This question is central to theories of ecological specialization, the maintenance of genetic variation, and sympatric speciation, yet relatively few experiments have examined costs of generalism over long-term adaptation. We tested whether selection in environments containing two resources limits a population's ability to adapt to the individual resources by comparing the fitness of replicate Escherichia coli populations evolved for 6000 generations in the presence of glucose or lactose alone (specialists), or in varying presentations of glucose and lactose together (generalists). We found that all populations had significant fitness increases in both resources, though the magnitude and rate of these increases differed. For the first 4000 generations, most generalist populations increased in fitness as quickly in the individual resources as the corresponding specialist populations. From 5000 generations, however, a widespread cost of adaptation affected all generalists, indicating a growing constraint on their abilities to adapt to two resources simultaneously. Our results indicate that costs of generalism are prevalent, but may influence evolutionary trajectories only after a period of cost-free adaptation.

The rule of declining adaptability in microbial evolution experiments

One of the most recurrent observations after two decades of microbial evolution experiments regards the dynamics of fitness change. In a given environment, low-fitness genotypes are recurrently observed to adapt faster than their more fit counterparts. Since adaptation is the main macroscopic outcome of Darwinian evolution, studying its patterns of change could potentially provide insight into key issues of evolutionary theory, from fixation dynamics to the genetic architecture of organisms. Here, we re-analyze several published datasets from experimental evolution with microbes and show that, despite large differences in the origin of the data, a pattern of inverse dependence of adaptability with fitness clearly emerges. In quantitative terms, it is remarkable to observe little if any degree of idiosyncrasy across systems as diverse as virus, bacteria and yeast. The universality of this phenomenon suggests that its emergence might be understood from general principles, giving rise to the exciting prospect that evolution might be statistically predictable at the macroscopic level. We discuss these possibilities in the light of the various theories of adaptation that have been proposed and delineate future directions of research.

The Utility of Fisher's Geometric Model in Evolutionary Genetics

The accumulation of data on the genomic bases of adaptation has triggered renewed interest in theoretical models of adaptation. Among these models, Fisher Geometric Model (FGM) has received a lot of attention over the last two decades. FGM is based on a continuous multidimensional phenotypic landscape, but it is for the emerging properties of individual mutation effects that it is mostly used. Despite an apparent simplicity and a limited number of parameters, FGM integrates a full model of mutation and epistatic interactions that allows the study of both beneficial and deleterious mutations, and subsequently the fate of evolving populations. In this review, I present the different properties of FGM and the qualitative and quantitative support they have received from experimental evolution data. I later discuss how to estimate the different parameters of the model and outline some future directions to connect FGM and the molecular determinants of adaptation.

The genomic landscape of compensatory evolution

The Genomic Landscape of Compensatory Evolution Laboratory selection experiment explains how organisms compensate for the loss of genes during evolution, and reveals the deleterious side-effects of this process when adapting to novel environments.

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.

While core cellular processes are generally conserved during evolution, the constituent genes differ somewhat between related species with similar lifestyles. Why should this be so? In this work, we propose that gene loss may initially be deleterious, but organisms can recover fitness by the accumulation of compensatory mutations elsewhere in the genome. To investigate this process in the laboratory, we investigated 180 haploid yeast strains, each of which initially displayed slow growth owing to the deletion of a single gene. Laboratory evolutionary experiments revealed that defects in a broad range of molecular processes can readily be compensated during evolution. Genomic analyses and functional assays demonstrated that compensatory evolution generates hidden genetic and physiological variation across parallel evolving lines, which can be revealed when the environment changes. Strikingly, despite nearly full recovery of fitness, the wild-type genomic expression pattern is generally not restored. Based on these results, we argue that genomes undergo major changes not simply to adapt to external conditions but also to compensate for previously accumulated deleterious mutations.

Parallel genome-wide fixation of ancestral alleles in partially outcrossing experimental populations of Caenorhabditis elegans

Experimental evolution studies, coupled with new advances in DNA sequencing technology, have become a powerful tool for exploring how populations respond to selection at the genomic level. Recent experiments in microbes typically have found evidence for multiple novel mutations, which are usually fixed. In contrast, in animal model systems, evolutionary responses seem to involve more modest changes in the frequencies of pre-existing alleles, probably because these populations outcross and are usually initialized with greater levels of standing variation. In this experiment, I used whole-genome resequencing to estimate allele frequencies and look for novel substitutions in experimentally evolved populations of Caenorhabditis elegans. These populations were founded with a fixed pair of deleterious mutations introgressed into multiple wild genetic backgrounds and allowed to evolve for 50 generations with a mixed mating system. There is evidence for some recombination between ancestral haplotypes, but selective sweeps seem to have resulted in the fixation of large chromosomal segments throughout most of the genome. In addition, a few new mutations were detected. Simulations suggest that strong selection and low outcrossing rates are likely explanations for the observed outcomes, consistent with earlier work showing large fitness increases in these populations over 50 generations. These results also show clear parallels to population genetic patterns in C. elegans in nature: recent selective sweeps, high linkage disequilibrium, and low effective recombination rates. Thus, the genomic consequences of selection depend heavily on the biology of the organism in question, including its mating system and levels of genetic variation.

Patterns and persistence of antibiotic resistance in faecal indicator bacteria from freshwater recreational beaches

AIMS

This study was conducted to determine antibiotic susceptibility patterns among the faecal indicator bacteria (FIB), Escherichia coli and enterococci, and to determine the potential for freshwater beaches to serve as reservoirs of resistance genes where transfer of resistant phenotypes takes place or de novo resistance may evolve.

METHODS AND RESULTS

One hundred and forty-seven E. coli and 150 enterococci collected from sand and water at recreational beaches along Lake Huron, Michigan, USA were screened against commonly used antibiotics. Resistance was apparent in both E. coli (19% resistant) and enterococci (65% resistant). Antibiotic-resistant E. coli were capable of growing in beach sand microcosms and were able to transfer a plasmid-encoded kanamycin-resistance gene in sand microcosms. Furthermore, resistant phenotypes were stable in the sand environment even in the absence of the corresponding antibiotic.

CONCLUSIONS

Antibiotic-resistant FIB were prevalent and persistent in the beach habitat.

SIGNIFICANCE AND IMPACT OF THE STUDY

Active populations of FIB at beaches express antibiotic resistance phenotypes and have the ability to transfer antibiotic resistance. These human-associated bacteria may be intermediaries in the movement of resistance between environmental and clinical reservoirs.

Weak selection and protein evolution

The "nearly neutral" theory of molecular evolution proposes that many features of genomes arise from the interaction of three weak evolutionary forces: mutation, genetic drift, and natural selection acting at its limit of efficacy. Such forces generally have little impact on allele frequencies within populations from generation to generation but can have substantial effects on long-term evolution. The evolutionary dynamics of weakly selected mutations are highly sensitive to population size, and near neutrality was initially proposed as an adjustment to the neutral theory to account for general patterns in available protein and DNA variation data. Here, we review the motivation for the nearly neutral theory, discuss the structure of the model and its predictions, and evaluate current empirical support for interactions among weak evolutionary forces in protein evolution. Near neutrality may be a prevalent mode of evolution across a range of functional categories of mutations and taxa. However, multiple evolutionary mechanisms (including adaptive evolution, linked selection, changes in fitness-effect distributions, and weak selection) can often explain the same patterns of genome variation. Strong parameter sensitivity remains a limitation of the nearly neutral model, and we discuss concave fitness functions as a plausible underlying basis for weak selection.

Life-history evolution and density-dependent growth in experimental populations of yeast

We studied the evolution of the correlation between growth rate r and yield K in experimental lineages of the yeast Saccharomyces cerevisiae. First, we isolated a single clone every approximately 250 generations from each of eight populations selected in a glucose-limited medium for 5000 generations at approximately 6.6 population doublings per day (20 clones per line × 8 lines) and measured its growth rate and yield in a new, galactose-limited medium (with ∼1.3 doubling per day). For most lines, r on galactose increased throughout the 5000 generations of selection on glucose whereas K on galactose declined. Next, we selected these 160 glucose-adapted clones in the galactose environment for approximately 120 generations and measured changes in r and K in galactose. In general, growth rate increased and yield declined, and clones that initially grew slowly on galactose improved more than did faster clones. We found a negative correlation between r and K among clones both within each line and across all clones. We provide evidence that this relationship is not heritable and is a negative environmental correlation rather than a genetic trade-off.

Epistasis can increase multivariate trait diversity in haploid non-recombining populations

We evaluate the effect of epistasis on genetically-based multivariate trait variation in haploid non-recombining populations. In a univariate setting, past work has shown that epistasis reduces genetic variance (additive plus epistatic) in a population experiencing stabilizing selection. Here we show that in a multivariate setting, epistasis also reduces total genetic variation across the entire multivariate trait in a population experiencing stabilizing selection. But, we also show that the pattern of variation across the multivariate trait can be more even when epistasis occurs compared to when epistasis is absent, such that some character combinations will have more genetic variance when epistasis occurs compared to when epistasis is absent. In fact, a measure of generalized multivariate trait variation can be increased by epistasis under weak to moderate stabilizing selection conditions, as well as neutral conditions. Likewise, a measure of conditional evolvability can be increased by epistasis under weak to moderate stabilizing selection and neutral conditions. We investigate the nature of epistasis assuming a multivariate-normal model genetic effects and investigate the nature of epistasis underlying the biophysical properties of RNA. Increased multivariate diversity occurs for populations that are infinite in size, as well as populations that are finite in size. Our model of finite populations is explicitly genealogical and we link our findings about the evenness of eigenvalues with epistasis to prior work on the genealogical mapping of epistatic effects.

Single-cell sequencing provides clues about the host interactions of segmented filamentous bacteria (SFB)

Segmented filamentous bacteria (SFB) are host-specific intestinal symbionts that comprise a distinct clade within the Clostridiaceae, designated Candidatus Arthromitus. SFB display a unique life cycle within the host, involving differentiation into multiple cell types. The latter include filaments that attach intimately to intestinal epithelial cells, and from which "holdfasts" and spores develop. SFB induce a multifaceted immune response, leading to host protection from intestinal pathogens. Cultivation resistance has hindered characterization of these enigmatic bacteria. In the present study, we isolated five SFB filaments from a mouse using a microfluidic device equipped with laser tweezers, generated genome sequences from each, and compared these sequences with each other, as well as to recently published SFB genome sequences. Based on the resulting analyses, SFB appear to be dependent on the host for a variety of essential nutrients. SFB have a relatively high abundance of predicted proteins devoted to cell cycle control and to envelope biogenesis, and have a group of SFB-specific autolysins and a dynamin-like protein. Among the five filament genomes, an average of 8.6% of predicted proteins were novel, including a family of secreted SFB-specific proteins. Four ADP-ribosyltransferase (ADPRT) sequence types, and a myosin-cross-reactive antigen (MCRA) protein were discovered; we hypothesize that they are involved in modulation of host responses. The presence of polymorphisms among mouse SFB genomes suggests the evolution of distinct SFB lineages. Overall, our results reveal several aspects of SFB adaptation to the mammalian intestinal tract.

Cost of antibiotic resistance and the geometry of adaptation

The distribution of effects of beneficial mutations is key to our understanding of biological adaptation. Yet, empirical estimates of this distribution are scarce, and its functional form is largely unknown. Theoretical models of adaptation predict that the functional form of this distribution should depend on the distance to the optimum. Here, we estimate the rate and distribution of adaptive mutations that compensate for the effect of a single deleterious mutation, which causes antibiotic resistance. Using a system with multiple molecular markers, we estimate the distribution of fitness effects of mutations at two distances from the adaptive peak in 60 populations of Escherichia coli. We find that beneficial mutations, which can contribute to compensatory evolution, occur at very high rates, of the order of 10⁻⁵ per genome per generation and can be detected within a few tens of generations. They cause an average fitness increase of 2.5% and 3.6%, depending on the cost of resistance, which is expected under Fisher's geometrical model of adaptation. Moreover, we provide the first description of the distribution of beneficial mutations, segregating during the process of compensatory evolution, to antibiotic resistances bearing different costs. Hence, these results have important implications to understanding the spread and maintenance of antibiotic resistance in bacteria.

Interactions among flower-size QTL of Mimulus guttatus are abundant but highly variable in nature

The frequency and character of interactions among genes influencing complex traits remain unknown. Our ignorance is most acute for segregating variation within natural populations, the epistasis most relevant for quantitative trait evolution. Here, we report a comprehensive survey of interactions among a defined set of flower-size QTL: loci polymorphic within a single natural population of yellow monkeyflower (Mimulus guttatus). We find that epistasis is typical. Observed phenotypes routinely differ from those predicted on the basis of direct allelic affects in the isogenic background, although the direction of deviations is highly variable. Across QTL pairs, there are significantly positive and negative interactions for every trait. Across traits, specific locus pairs routinely exhibit both positive and negative interactions. There was a tendency for negative epistasis to accompany positive direct effects and vice versa for the trait of corolla width, which may be due, at least in part, to the fact that QTL were identified from their direct effects on this trait.

Negative epistasis between beneficial mutations in an evolving bacterial population

Epistatic interactions between mutations play a prominent role in evolutionary theories. Many studies have found that epistasis is widespread, but they have rarely considered beneficial mutations. We analyzed the effects of epistasis on fitness for the first five mutations to fix in an experimental population of Escherichia coli. Epistasis depended on the effects of the combined mutations--the larger the expected benefit, the more negative the epistatic effect. Epistasis thus tended to produce diminishing returns with genotype fitness, although interactions involving one particular mutation had the opposite effect. These data support models in which negative epistasis contributes to declining rates of adaptation over time. Sign epistasis was rare in this genome-wide study, in contrast to its prevalence in an earlier study of mutations in a single gene.

Adaptive landscapes in evolving populations of Pseudomonas fluorescens

The repeatability of adaptive evolution depends on the ruggedness of the underlying adaptive landscape. We contrasted the relative ruggedness of two adaptive landscapes by measuring the variance in fitness and metabolic phenotype within and among genetically distinct strains of Pseudomonas fluorescens in two environments differing only in the carbon source provided (glucose vs. xylose). Fitness increased in all lines, plateauing in one environment but not the other. The pattern of variance in fitness among replicate lines was unique to the selection environment; it increased over the course of the experiment in xylose but not in glucose. Metabolic phenotypes displayed two results: (1) populations adapted via changes that were distinctive to their selection environment, and (2) endpoint phenotypes were less variable in glucose than in xylose. These results indicate that although the response to selection is highly repeatable at the level of fitness, the underlying genetic routes taken were different for each environment and more variable in xylose. We suggest that this reflects a more rugged adaptive landscape in xylose compared to glucose. Our study demonstrates the utility of using replicate selection lines with different evolutionary starting points to try and quantify the relative ruggedness of adaptive landscapes.

Hidden randomness between fitness landscapes limits reverse evolution

In biological evolution, adaptations to one environment can in some cases reverse adaptations to another environment. To study this "reverse evolution" on a genotypic level, we measured the fitness of E. coli strains with each possible combination of five mutations in an antibiotic-resistance gene in two distinct antibiotic environments. While adaptations to one environment generally lower fitness in the other, we find that reverse evolution is rarely possible and falls as the complexity of adaptations increases, suggesting a probabilistic, molecular form of Dollo's law.

Fitness recovery and compensatory evolution in natural mutant lines of C. elegans

Deleterious mutation accumulation plays a central role in evolutionary genetics, conservation biology, human health, and evolutionary medicine (e.g., methods of viral attenuation for live vaccines). It is therefore important to understand whether and how quickly populations with accumulated deleterious mutational loads can recover fitness through adaptive evolution. We used laboratory experimental evolution with four long-term mutation-accumulation (MA) lines of Caenorhabditis elegans nematodes to study the dynamics of such fitness evolution. We previously showed that when homozygous mutant populations are evolved in large population sizes, they can rapidly achieve wild-type fitness through the accumulation of new beneficial or compensatory epistatic mutations. Here, we expand this approach to demonstrate that when replicate lineages are initiated from the same mutant genotype, phenotypic evolution is only sometimes repeatable. MA genotypes that recovered ancestral fitness in the previous experiment did not always do so here. Further, the pattern of adaptive evolution in independently evolved replicates was contingent upon the MA genotype and varied among fitness-related traits. Our findings suggest that new beneficial mutations can drive rapid fitness evolution, but that the adaptive process is rendered somewhat unpredictable by its susceptibility to chance events and sensitivity to the evolutionary history of the starting population.

Hypermutability and compensatory adaptation in antibiotic-resistant bacteria

Hypermutable (mutator) bacteria have been associated with the emergence of antibiotic resistance. A simple yet untested prediction is that mutator bacteria are able to compensate more quickly for pleiotropic fitness costs often associated with resistance, resulting in the maintenance of resistance in the absence of antibiotic selection. By using experimental populations of a wild-type and a mutator genotype of the pathogenic bacterium Pseudomonas aeruginosa, we show that mutator bacteria can evolve resistance to antibiotics more rapidly than wild-type bacteria and, crucially, that mutators are better able to compensate for the fitness cost of resistance, to the extent that all costs of resistance were entirely compensated for in mutators. When competed against immigrant antibiotic-susceptible bacteria in the absence of antibiotics, antibiotic resistance remained at a high level in mutator populations but disappeared in wild-type populations. These results suggest that selection for mutations that offset the fitness cost associated with antibiotic resistance may help to explain the high frequency of mutator bacteria and antibiotic resistance observed in chronic infections.

Parallel changes in host resistance to viral infection during 45,000 generations of relaxed selection

The dynamics of host susceptibility to parasites are often influenced by trade-offs between the costs and benefits of resistance. We assayed changes in the resistance to three viruses in six lines of Escherichia coli that had been evolving for almost 45,000 generations in their absence. The common ancestor of these lines was completely resistant to T6, partially resistant to T6* (a mutant of T6 with altered host range), and sensitive to λ. None of the populations changed with respect to resistance to T6, whereas all six evolved increased susceptibility to T6*, probably ameliorating a cost of resistance. More surprisingly, however, the majority of lines evolved complete resistance to λ, despite not encountering that virus during this period. By coupling our results with previous work, we infer that resistance to λ evolved as a pleiotropic effect of a beneficial mutation that downregulated an unused metabolic pathway. The strong parallelism between the lines implies that selection had almost deterministic effects on the evolution of these patterns of host resistance. The opposite outcomes for resistance to T6* and λ demonstrate that the evolution of host resistance under relaxed selection cannot be fully predicted by simple trade-off models.

The population genetics of antibiotic resistance: integrating molecular mechanisms and treatment contexts

Despite efforts from a range of disciplines, our ability to predict and combat the evolution of antibiotic resistance in pathogenic bacteria is limited. This is because resistance evolution involves a complex interplay between the specific drug, bacterial genetics and both natural and treatment ecology. Incorporating details of the molecular mechanisms of drug resistance and ecology into evolutionary models has proved useful in predicting the dynamics of resistance evolution. However, putting these models to practical use will require extensive collaboration between mathematicians, molecular biologists, evolutionary ecologists and clinicians.

Compensatory molecular evolution of HC-Pro, an RNA-silencing suppressor from a plant RNA virus

RNA silencing is a eukaryotic mechanism involved in several cellular processes, one example being a sequence-specific antiviral defense. Many plant viruses have developed counterdefensive proteins that in many instances are multifunctional, such as helper component protease (HC-Pro) of Tobacco etch virus (TEV). In a previous work, a collection of mutants with amino acid replacements in TEV HC-Pro was generated, and their effects in the capacity of suppressing RNA silencing were quantified in a transient expression assay. In this study, three mutations that caused a reduction in suppression activity and three that increased it were used to create replicate experimental lineages that were evolved through serial passages. We have evaluated the number of genotypic changes that occurred during evolution in HC-Pro and their phenotypic effects on virus viability, virulence, and suppression of RNA silencing. In no instance did the original mutation revert to the wildtype (WT) sequence. In several cases, fixed mutations were canonical compensatory changes, returning the suppressor activity to the WT HC-Pro value, pointing to the existence of stabilizing selection pressures and pleiotropic effects of the introduced original mutations. However, in other instances, the fixed mutations were overcompensatory, driving the activity of the mutant beyond the optimal value. Negative epistatic effects among beneficial mutations as well as decompensatory epistasis also play an important role during compensatory evolution of RNA-silencing suppression.