Genotype‐by‐environment interactions due to antibiotic resistance and adaptation in Escherichia coli

Mutations that are beneficial in one environment can have different fitness effects in other environments. In the context of antibiotic resistance, the resulting genotype‐by‐environment interactions potentially make selection on resistance unpredictable in heterogeneous environments. Furthermore, resistant bacteria frequently fix additional mutations during evolution in the absence of antibiotics. How do these two types of mutations interact to determine the bacterial phenotype across different environments? To address this, I used Escherichia coli as a model system, measuring the effects of nine different rifampicin resistance mutations on bacterial growth in 31 antibiotic‐free environments. I did this both before and after approximately 200 generations of experimental evolution in antibiotic‐free conditions (LB medium), and did the same for the antibiotic‐sensitive wild type after adaptation to the same environment. The following results were observed: (i) bacteria with and without costly resistance mutations adapted to experimental conditions and reached similar levels of competitive fitness; (ii) rifampicin resistance mutations and adaptation to LB both indirectly altered growth in other environments; and (iii) resistant‐evolved genotypes were more phenotypically different from the ancestor and from each other than resistant‐nonevolved and sensitive‐evolved genotypes. This suggests genotype‐by‐environment interactions generated by antibiotic resistance mutations, observed previously in short‐term experiments, are more pronounced after adaptation to other types of environmental variation, making it difficult to predict long‐term selection on resistance mutations from fitness effects in a single environment.     

The cost of multiple drug resistance in Pseudomonas aeruginosa

The spread of bacterial antibiotic resistance mutations is thought to be constrained by their pleiotropic fitness costs. Here we investigate the fitness costs of resistance in the context of the evolution of multiple drug resistance (MDR), by measuring the cost of acquiring streptomycin resistance mutations (StrepR) in independent strains of the bacterium Pseudomonas aeruginosa carrying different rifampicin resistance (RifR) mutations. In the absence of antibiotics, StrepR mutations are associated with similar fitness costs in different RifR genetic backgrounds. The cost of StrepR mutations is greater in a rifampicin‐sensitive (RifS) background, directly demonstrating antagonistic epistasis between resistance mutations. In the presence of rifampicin, StrepR mutations have contrasting effects in different RifR backgrounds: StrepR mutations have no detectable costs in some RifR backgrounds and massive fitness costs in others. Our results clearly demonstrate the importance of epistasis and genotype‐by‐environment interactions for the evolution of MDR.     

EFFECTS OF TEMPERATURE ON THE FITNESS COST OF RESISTANCE TO BACTERIOPHAGE T4 IN ESCHERICHIA COLI.

Resistance to predation, herbivory, or disease often comes at a cost such that resistant genotypes are competitively inferior to their sensitive counterparts in the absence of predators, herbivores, or pathogens. The effects of this trade‐off on natural populations depend on its sensitivity to environmental changes. We used Escherichia coli and bacteriophage T4 as a model predator/prey system to study the effects of temperature on the cost of resistance. An array of independent T4‐resistant mutants, derived from a single ancestral strain of E. coli B, had a mean reduction in competitive fitness that depended strongly on environmental temperature; the cost of resistance generally increased with temperature. Genetic variance for fitness among phage‐resistant mutants also depended on temperature; however, genetic variance increased at high and low thermal extremes. These results suggest that temperature is likely to be an important determinant of the consequences of predation in natural communities. We also discuss the underlying mechanistic basis for the cost of resistance in this system and its interaction with temperature.     

Environmental heterogeneity generates opposite gene‐by‐environment interactions for two fitness‐related traits within a population

Theory predicts that environmental heterogeneity offers a potential solution to the maintenance of genetic variation within populations, but empirical evidence remains sparse. The live‐bearing fish Xiphophorus variatus exhibits polymorphism at a single locus, with different alleles resulting in up to five distinct melanistic “tailspot” patterns within populations. We investigated the effects of heterogeneity in two ubiquitous environmental variables (temperature and food availability) on two fitness‐related traits (upper thermal limits and body condition) in two different tailspot types (wild‐type and upper cut crescent). We found gene‐by‐environment (G × E) interactions between tailspot type and food level affecting upper thermal limits (UTL), as well as between tailspot type and thermal environment affecting body condition. Exploring mechanistic bases underlying these G × E patterns, we found no differences between tailspot types in hsp70 gene expression despite significant overall increases in expression under both thermal and food stress. Similarly, there was no difference in routine metabolic rates between the tailspot types. The reversal of relative performance of the two tailspot types under different environmental conditions revealed a mechanism by which environmental heterogeneity can balance polymorphism within populations through selection on different fitness‐related traits.     

Error‐prone initiation factor 2 mutations reduce the fitness cost of antibiotic resistance

Mutations in the fmt gene (encoding formyl methionine transferase) that eliminate formylation of initiator tRNA (Met‐tRNA_i) confer resistance to the novel antibiotic class of peptide deformylase inhibitors (PDFIs) while concomitantly reducing bacterial fitness. Here we show in Salmonella typhimurium that novel mutations in initiation factor 2 (IF2) located outside the initiator tRNA binding domain can partly restore fitness of fmt mutants without loss of antibiotic resistance. Analysis of initiation of protein synthesis in vitro showed that with non‐formylated Met‐tRNA_i IF2 mutants initiated much faster than wild‐type IF2, whereas with formylated fMet‐tRNA_i the initiation rates were similar. Moreover, the increase in initiation rates with Met‐tRNA_i conferred by IF2 mutations in vitro correlated well with the increase in growth rate conferred by the same mutations in vivo, suggesting that the mutations in IF2 compensate formylation deficiency by increasing the rate of in vivo initiation with Met‐tRNA_i. IF2 mutants had also a high propensity for erroneous initiation with elongator tRNAs in vitro, which could account for their reduced fitness in vivo in a formylation‐proficient strain. More generally, our results suggest that bacterial protein synthesis is mRNA‐limited and that compensatory mutations in IF2 could increase the persistence of PDFI‐resistant bacteria in clinical settings.     

EVOLUTIONARY REVERSALS OF ANTIBIOTIC RESISTANCE IN EXPERIMENTAL POPULATIONS OF PSEUDOMONAS AERUGINOSA

Antibiotic resistance mutations are accompanied by a fitness cost, and two mechanisms allow bacteria to adapt to this cost once antibiotic use is halted. First, it is possible for resistance to revert; second, it is possible for bacteria to adapt to the cost of resistance by compensatory mutations. Unfortunately, reversion to antibiotic sensitivity is rare, but the underlying factors that prevent reversion remain obscure. Here, we directly study the evolutionary dynamics of reversion by experimentally mimicking reversion mutations—sensitives—in populations of rifampicin‐resistant Pseudomonas aeruginosa. We show that, in our populations, most sensitives are lost due to genetic drift when they are rare. However, clonal interference from lineages carrying compensatory mutations causes a dramatic increase in the time to fixation of sensitives that escape genetic drift, and mutations surpassing the sensitives’ fitness are capable of driving transiently common sensitive lineages to extinction. Crucially, we show that the constraints on reversion arising from clonal interference are determined by the potential for compensatory adaptation of the resistant population. Although the cost of resistance provides the incentive for reversion, our study demonstrates that both the cost of resistance and the intrinsic evolvability of resistant populations interact to determine the rate and likelihood of reversion.     

The interplay between drug resistance and fitness in malaria parasites

Controlling the spread of antimalarial drug resistance, especially resistance of Plasmodium falciparum to artemisinin‐based combination therapies, is a high priority. Available data indicate that, as with other microorganisms, the spread of drug‐resistant malaria parasites is limited by fitness costs that frequently accompany resistance. Resistance‐mediating polymorphisms in malaria parasites have been identified in putative drug transporters and in target enzymes. The impacts of these polymorphisms on parasite fitness have been characterized in vitro and in animal models. Additional insights have come from analyses of samples from clinical studies, both evaluating parasites under different selective pressures and determining the clinical consequences of infection with different parasites. With some exceptions, resistance‐mediating polymorphisms lead to malaria parasites that, compared with wild type, grow less well in culture and in animals, and are replaced by wild type when drug pressure diminishes in the clinical setting. In some cases, the fitness costs of resistance may be offset by compensatory mutations that increase virulence or changes that enhance malaria transmission. However, not enough is known about effects of resistance mediators on parasite fitness. A better appreciation of the costs of fitness‐mediating mutations will facilitate the development of optimal guidelines for the treatment and prevention of malaria.     

Fisher's geometrical model and the mutational patterns of antibiotic resistance across dose gradients

Fisher's geometrical model (FGM) has been widely used to depict the fitness effects of mutations. It is a general model with few underlying assumptions that gives a large and comprehensive view of adaptive processes. It is thus attractive in several situations, for example adaptation to antibiotics, but comes with limitations, so that more mechanistic approaches are often preferred to interpret experimental data. It might be possible however to extend FGM assumptions to better account for mutational data. This is theoretically challenging in the context of antibiotic resistance because resistance mutations are assumed to be rare. In this article, we show with Escherichia coli how the fitness effects of resistance mutations screened at different doses of nalidixic acid vary across a dose‐gradient. We found experimental patterns qualitatively consistent with the basic FGM (rate of resistance across doses, gamma distributed costs) but also unexpected patterns such as a decreasing mean cost of resistance with increasing screen dose. We show how different extensions involving mutational modules and variations in trait covariance across environments, can be discriminated based on these data. Overall, simple extensions of the FGM accounted well for complex mutational effects of resistance mutations across antibiotic doses.     

ANTIBIOTIC RESISTANCE AND STRESS IN THE LIGHT OF FISHER'S MODEL

The role of mutations in evolution depends upon the distribution of their effects on fitness. This distribution is likely to depend on the environment. Indeed genotype‐by‐environment interactions are key for the process of local adaptation and ecological specialization. An important trait in bacterial evolution is antibiotic resistance, which presents a clear case of change in the direction of selection between environments with and without antibiotics. Here, we study the distribution of fitness effects of mutations, conferring antibiotic resistance to Escherichia coli, in benign and stressful environments without drugs. We interpret the distributions in the light of a fitness landscape model that assumes a single fitness peak. We find that mutation effects (s) are well described by a shifted gamma distribution, with a shift parameter that reflects the distance to the fitness peak and varies across environments. Consistent with the theoretical predictions of Fisher's geometrical model, with a Gaussian relationship between phenotype and fitness, we find that the main effect of stress is to increase the variance in s. Our findings are in agreement with the results of a recent meta‐analysis, which suggest that a simple fitness landscape model may capture the variation of mutation effects across species and environments.     

Evolution of plant resistance and tolerance to frost damage

Plant defence against any type of stress may involve resistance (traits that reduce damage) or tolerance (traits that reduce the negative fitness impacts of damage). These two strategies have been proposed as redundant evolutionary alternatives. A late‐season frost enabled us to estimate natural selection and genetic constraints on the evolution of frost resistance and tolerance in a wild plant species. We employed a genetic selection analysis (which is unbiased by environmental correlations between traits and fitness) on 75 paternal half‐sibling families of annual wild radish [Raphanus raphanistrum (Brassicaceae)]. In an experimental population in southern Ontario, we found strong selection favouring plant resistance to frost, but selection against tolerance to frost. The selection against tolerance may have been caused by a cost of tolerance, as we provide evidence for a negative genetic correlation between tolerance and fitness in the absence of frost damage. Although we found no evidence for the theoretically predicted trade‐off between frost tolerance and resistance among our families, we did detect negative correlational selection acting on the two traits, indicating that natural selection favoured high resistance combined with low tolerance and low resistance coupled with high tolerance, but not high or low levels of both traits together. There were few genetic correlations between the measured traits overall, but frost tolerance was negatively correlated with initial seed mass, and frost resistance was positively correlated with resistance to insect herbivory. Periodic episodes of strong selection such as that caused by the late‐season frost may be disproportionately important in evolution, and are likely becoming more common because of human alterations of the environment.     

Stressful environments can indirectly select for increased longevity

Longevity is modulated by a range of conserved genes in eukaryotes, but it is unclear how variation in these genes contributes to the evolution of longevity in nature. Mutations that increase life span in model organisms typically induce trade‐offs which lead to a net reduction in fitness, suggesting that such mutations are unlikely to become established in natural populations. However, the fitness consequences of manipulating longevity have rarely been assessed in heterogeneous environments, in which stressful conditions are encountered. Using laboratory selection experiments, we demonstrate that long‐lived, stress‐resistant Caenorhabditis elegans age‐1(hx546) mutants have higher fitness than the wild‐type genotype if mixed genotype populations are periodically exposed to high temperatures when food is not limited. We further establish, using stochastic population projection models, that the age‐1(hx546) mutant allele can confer a selective advantage if temperature stress is encountered when food availability also varies over time. Our results indicate that heterogeneity in environmental stress may lead to altered allele frequencies over ecological timescales and indirectly drive the evolution of longevity. This has important implications for understanding the evolution of life‐history strategies.     

Mutability in Pseudomonas viridiflava as a programmed balance between antibiotic resistance and pathogenicity

Mutable bacterial cells are defective in their DNA repair system and often have a phenotype different from that of their wild‐type counterparts. In human bacterial pathogens, the mutable and hypermutable phenotypes are often associated with general antibiotic resistance. Here, we quantified the occurrence of mutable cells in Pseudomonas viridiflava, a phytopathogenic bacterium in the P. syringae complex with a broad host range and capacity to live as a saprophyte. Two phenotypic variants (transparent and mucoid) were produced by this bacterium. The transparent variant had a mutator phenotype, showed general antibiotic resistance and could not induce disease on the plant species tested (bean). In contrast, the mucoid variant did not display mutability or resistance to antibiotics and was capable of inducing disease on bean. Both the transparent and mucoid variants were less fit when grown in vitro, whereas, in planta, both of the variants and wild‐types attained similar population densities. Given the importance of the methyl‐directed mismatch repair system (MMR) in the occurrence of mutable and hypermutable cells in human bacterial pathogens, we investigated whether mutations in mut genes were associated with mutator transparent cells in P. viridiflava. Our results showed no mutations in MMR genes in any of the P. viridiflava cells tested. Here, we report that a high mutation rate and antibiotic resistance are inversely correlated with pathogenicity in P. viridiflava, but are not associated with mutations in MMR. In addition, P. viridiflava variants differ from variants produced by other phytopathogenic bacteria in the absence of reversion to the wild‐type phenotype.     

Fitness effects of mutation: testing genetic redundancy in Arabidopsis thaliana

Screens of organisms with disruptive mutations in a single gene often fail to detect phenotypic consequences for the majority of mutants. One explanation for this phenomenon is that the presence of paralogous loci provides genetic redundancy. However, it is also possible that the assayed traits are affected by few loci, that effects could be subtle or that phenotypic effects are restricted to certain environments. We assayed a set of T‐DNA insertion mutant lines of Arabidopsis thaliana to determine the frequency with which mutation affected fitness‐related phenotypes. We found that between 8% and 42% of the assayed lines had altered fitness from the wild type. Furthermore, many of these lines exhibited fitness greater than the wild type. In a second experiment, we grew a subset of the lines in multiple environments and found whether a T‐DNA insert increased or decreased fitness traits depended on the assay environment. Overall, our evidence contradicts the hypothesis that genetic redundancy is a common phenomenon in A. thaliana for fitness traits. Evidence for redundancy from prior screens of knockout mutants may often be an artefact of the design of the phenotypic assays which have focused on less complex phenotypes than fitness and have used single environments. Finally, our study adds to evidence that beneficial mutations may represent a significant component of the mutational spectrum of A. thaliana.     

Epistatic interactions between ancestral genotype and beneficial mutations shape evolvability in Pseudomonas aeruginosa

The idea that interactions between mutations influence adaptation by driving populations to low and high fitness peaks on adaptive landscapes is deeply ingrained in evolutionary theory. Here, we investigate the impact of epistasis on evolvability by challenging populations of two Pseudomonas aeruginosa clones bearing different initial mutations (in rpoB conferring rifampicin resistance, and the type IV pili gene network) to adaptation to a medium containing l ‐serine as the sole carbon source. Despite being initially indistinguishable in fitness, populations founded by the two ancestral genotypes reached different fitness following 300 generations of evolution. Genome sequencing revealed that the difference could not be explained by acquiring mutations in different targets of selection; the majority of clones from both ancestors converged on one of the following two strategies: (1) acquiring mutations in either PA2449 (gcsR, an l ‐serine‐metabolism RpoN enhancer binding protein) or (2) protease genes. Additionally, populations from both ancestors converged on loss‐of‐function mutations in the type IV pili gene network, either due to ancestral or acquired mutations. No compensatory or reversion mutations were observed in RNA polymerase (RNAP) genes, in spite of the large fitness costs typically associated with mutations in rpoB. Although current theory points to sign epistasis as the dominant constraint on evolvability, these results suggest that the role of magnitude epistasis in constraining evolvability may be underappreciated. The contribution of magnitude epistasis is likely to be greatest under the biologically relevant mutation supply rates that make back mutations probabilistically unlikely.     

Ribosomal subunits affected by antibiotic resistance mutations at seven chloroplast loci in Chlamydomonas reinhardtii

Summary Mutations at seven recombinationally distinct chloroplast loci confer antibiotic resistance on chloroplast ribosomes of the green alga Chlamydomonas reinhardtii. Assays of polynucleotide-directed amino acid incorporation by ribosomes reconstituted from mutant and wild type subunits demonstrate that streptomycin, neamine/kanamycin and spectinomycin resistance mutations specifically affect the small ribosomal subunit, whereas mutations to erythromycin resistance affect the large subunit. Although in each case the subunit site of antibiotic resistance is the same as that observed in analogous mutations in Escherichia coli, the number of loci conferring resistance to a given antibiotic differs in the two organisms. We have previously shown that streptomycin resistance mutations in Chlamydomonas map at five discrete loci (one nuclear and four chloroplast), and that mutations to neamine/kanamycin and spectinomycin resistance appear to define a single chloroplast locus. Results presented here confirm our previous report that all chloroplast erythromycin resistance mutations isolated to date fall into two recombinationally distinct loci, and indicate that mutants at one of these loci may be further divided on the basis of their level of cross resistance to other macrolide antibiotics.     

Decreased thermal niche breadth as a trade-off of antibiotic resistance

Evolutionary theory predicts that adaptations, including antibiotic resistance, should come with associated fitness costs; yet, many resistance mutations seemingly contradict this prediction by inducing no growth rate deficit. However, most growth assays comparing sensitive and resistant strains have been performed under a narrow range of environmental conditions, which do not reflect the variety of contexts that a pathogenic bacterium might encounter when causing infection. We hypothesized that reduced niche breadth, defined as diminished growth across a diversity of environments, can be a cost of antibiotic resistance. Specifically, we test whether chloramphenicol-resistant Escherichia coli incur disproportionate growth deficits in novel thermal conditions. Here we show that chloramphenicol-resistant bacteria have greater fitness costs at novel temperatures than their antibiotic-sensitive ancestors. In several cases, we observed no resistance cost in growth rate at the historic temperature but saw diminished growth at warmer and colder temperatures. These results were consistent across various genetic mechanisms of resistance. Thus, we propose that decreased thermal niche breadth is an under-documented fitness cost of antibiotic resistance. Furthermore, these results demonstrate that the cost of antibiotic resistance shifts rapidly as the environment changes; these context-dependent resistance costs should select for the rapid gain and loss of resistance as an evolutionary strategy.Subject terms: Bacterial evolution, Microbial ecology, Antibiotics     

The Distribution of Fitness Effects of Beneficial Mutations in Pseudomonas aeruginosa

Understanding how beneficial mutations affect fitness is crucial to our understanding of adaptation by natural selection. Here, using adaptation to the antibiotic rifampicin in the opportunistic pathogen Pseudomonas aeruginosa as a model system, we investigate the underlying distribution of fitness effects of beneficial mutations on which natural selection acts. Consistent with theory, the effects of beneficial mutations are exponentially distributed where the fitness of the wild type is moderate to high. However, when the fitness of the wild type is low, the data no longer follow an exponential distribution, because many beneficial mutations have large effects on fitness. There is no existing population genetic theory to explain this bias towards mutations of large effects, but it can be readily explained by the underlying biochemistry of rifampicin–RNA polymerase interactions. These results demonstrate the limitations of current population genetic theory for predicting adaptation to severe sources of stress, such as antibiotics, and they highlight the utility of integrating statistical and biophysical approaches to adaptation.     

Sexual selection on spontaneous mutations strengthens the between‐sex genetic correlation for fitness

A proposed benefit to sexual selection is that it promotes purging of deleterious mutations from populations. For this benefit to be realized, sexual selection, which is usually stronger on males, must purge mutations deleterious to both sexes. Here, we experimentally test the hypothesis that sexual selection on males purges deleterious mutations that affect both male and female fitness. We measured male and female fitness in two panels of spontaneous mutation‐accumulation lines of the fly, Drosophila serrata, each established from a common ancestor. One panel of mutation accumulation lines limited both natural and sexual selection (LS lines), whereas the other panel limited natural selection, but allowed sexual selection to operate (SS lines). Although mutation accumulation caused a significant reduction in male and female fitness in both the LS and SS lines, sexual selection had no detectable effect on the extent of the fitness reduction. Similarly, despite evidence of mutational variance for fitness in males and females of both treatments, sexual selection had no significant impact on the amount of mutational genetic variance for fitness. However, sexual selection did reshape the between‐sex correlation for fitness: significantly strengthening it in the SS lines. After 25 generations, the between‐sex correlation for fitness was positive but considerably less than one in the LS lines, suggesting that, although most mutations had sexually concordant fitness effects, sex‐limited, and/or sex‐biased mutations contributed substantially to the mutational variance. In the SS lines this correlation was strong and could not be distinguished from unity. Individual‐based simulations that mimick the experimental setup reveal two conditions that may drive our results: (1) a modest‐to‐large fraction of mutations have sex‐limited (or highly sex‐biased) fitness effects, and (2) the average fitness effect of sex‐limited mutations is larger than the average fitness effect of mutations that affect both sexes similarly.     

Rapid Decline of Ceftazidime Resistance in Antibiotic-Free and Sublethal Environments Is Contingent on Genetic Background

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