EPISTASIS BETWEEN NEW MUTATIONS AND GENETIC BACKGROUND AND A TEST OF GENETIC CANALIZATION

Abstract The importance for fitness of epistatic interactions among mutations is poorly known, yet epistasis can exert important effects on the dynamics of evolving populations. We showed previously that epistatic interactions are common between pairs of random insertion mutations in the bacterium Escherichia coli . In this paper, we examine interactions between these mutations and other mutations by transducing each of twelve insertion mutations into two genetic backgrounds, one ancestral and the other having evolved in, and adapted to, a defined laboratory environment for 10,000 generations. To assess the effect of the mutation on fitness, we allowed each mutant to compete against its unmutated counterpart in that same environment. Overall, there was a strong positive correlation between the mutational effects on the two genetic backgrounds. Nonetheless, three of the twelve mutations had significantly different effects on the two backgrounds, indicating epistasis. There was no significant tendency for the mutations to be less harmful on the derived background. Thus, there is no evidence supporting the hypothesis that the derived bacteria had adapted, in part, by becoming buffered against the harmful effects of mutations.     

Transition from positive to neutral in mutation fixation along with continuing rising fitness in thermal adaptive evolution

It remains to be determined experimentally whether increasing fitness is related to positive selection, while stationary fitness is related to neutral evolution. Long-term laboratory evolution in Escherichia coli was performed under conditions of thermal stress under defined laboratory conditions. The complete cell growth data showed common continuous fitness recovery to every 2°C or 4°C stepwise temperature upshift, finally resulting in an evolved E. coli strain with an improved upper temperature limit as high as 45.9°C after 523 days of serial transfer, equivalent to 7,560 generations, in minimal medium. Two-phase fitness dynamics, a rapid growth recovery phase followed by a gradual increasing growth phase, was clearly observed at diverse temperatures throughout the entire evolutionary process. Whole-genome sequence analysis revealed the transition from positive to neutral in mutation fixation, accompanied with a considerable escalation of spontaneous substitution rate in the late fitness recovery phase. It suggested that continually increasing fitness not always resulted in the reduction of genetic diversity due to the sequential takeovers by fit mutants, but caused the accumulation of a considerable number of mutations that facilitated the neutral evolution.     

Evolutionary rescue by beneficial mutations in environments that change in space and time

A factor that may limit the ability of many populations to adapt to changing conditions is the rate at which beneficial mutations can become established. We study the probability that mutations become established in changing environments by extending the classic theory for branching processes. When environments change in time, under quite general conditions, the establishment probability is approximately twice the ‘effective selection coefficient’, whose value is an average that gives most weight to a mutant''s fitness in the generations immediately after it appears. When fitness varies along a gradient in a continuous habitat, increased dispersal generally decreases the chance a mutation establishes because mutations move out of areas where they are most adapted. When there is a patch of favourable habitat that moves in time, there is a maximum speed of movement above which mutations cannot become established, regardless of when and where they first appear. This critical speed limit, which is proportional to the mutation''s maximum selective advantage, represents an absolute constraint on the potential of locally adapted mutations to contribute to evolutionary rescue.     

Evolutionary adaptation to environmental pH in experimental lineages of Escherichia coli

This study uses the enteric bacterium Escherichia coli as an experimental system to examine evolutionary responses of bacteria to an environmental acidic-alkaline range between pH 5.3 and 7.8 (15-5000 nM [H(+)]). Our goal was both to test general hypotheses about adaptation to abiotic variables and to provide insights into how coliform organisms might respond to changing conditions inside and outside of hosts. Six replicate lines of E. coli evolved for 2000 generations at one of four different constant pH conditions: pH 5.3, 6.3, 7.0, or 7.8. Direct adaptation to the evolutionary environment, as well as correlated changes in other environments, was measured as a change in fitness relative to the ancestor in direct competition experiments. The pH 5.3 group had the highest fitness gains, with a highly significant increase of 20%. The pH 7.8 group had far less significant gains and much higher variance among its lines. Analysis of individual lines within these two groups revealed complex patterns of adaptation: all of the pH 5.3 lines exhibited trade-offs (reduced fitness in another environment), but only 33% of the pH 7.8 lines showed such trade-offs and one of the pH 7.8 lines demonstrated exaptation by improving fitness in the pH 5.3 environment. Although there was also prevalent exaptation in other groups to the acidic environment, there were no such cases of exaptation to alkalinity. Comparison across the entire experimental pH range revealed that the most acidic lines, the pH 5.3 group, were all specialists, in contrast to the pH 6.3 lines, which were almost all generalists. That is, although none of the pH 5.3 lines showed any correlated fitness gains, all of the pH 6.3 lines did.     

Genomewide Patterns of Substitution in Adaptively Evolving Populations of the RNA Bacteriophage MS2

Experimental evolution of bacteriophage provides a powerful means of studying the genetics of adaptation, as every substitution contributing to adaptation can be identified and characterized. Here, I use experimental evolution of MS2, an RNA bacteriophage, to study its adaptive response to a novel environment. To this end, three lines of MS2 were adapted to rapid growth and lysis at cold temperature for a minimum of 50 phage generations and subjected to whole-genome sequencing. Using this system, I identified adaptive substitutions, monitored changes in frequency of adaptive mutations through the course of the experiment, and measured the effect on phage growth rate of each substitution. All three lines showed a substantial increase in fitness (a two- to threefold increase in growth rate) due to a modest number of substitutions (three to four). The data show some evidence that the substitutions occurring early in the experiment have larger beneficial effects than later ones, in accordance with the expected diminishing returns relationship between the fitness effects of a mutation and its order of substitution. Patterns of molecular evolution seen here—primarily a paucity of hitchhiking mutations—suggest an abundant supply of beneficial mutations in this system. Nevertheless, some beneficial mutations appear to have been lost, possibly due to accumulation of beneficial mutations on other genetic backgrounds, clonal interference, and negatively epistatic interactions with other beneficial mutations.     

Diverse phenotypic and genetic responses to short‐term selection in evolving Escherichia coli populations

Beneficial mutations fuel adaptation by altering phenotypes that enhance the fit of organisms to their environment. However, the phenotypic effects of mutations often depend on ecological context, making the distribution of effects across multiple environments essential to understanding the true nature of beneficial mutations. Studies that address both the genetic basis and ecological consequences of adaptive mutations remain rare. Here, we characterize the direct and pleiotropic fitness effects of a collection of 21 first‐step beneficial mutants derived from naïve and adapted genotypes used in a long‐term experimental evolution of Escherichia coli. Whole‐genome sequencing was able to identify the majority of beneficial mutations. In contrast to previous studies, we find diverse fitness effects of mutations selected in a simple environment and few cases of genetic parallelism. The pleiotropic effects of these mutations were predominantly positive but some mutants were highly antagonistic in alternative environments. Further, the fitness effects of mutations derived from the adapted genotypes were dramatically reduced in nearly all environments. These findings suggest that many beneficial variants are accessible from a single point on the fitness landscape, and the fixation of alternative beneficial mutations may have dramatic consequences for niche breadth reduction via metabolic erosion.     

Improved methods of cultivation and production of deuteriated proteins from E. coli strains grown on fully deuteriated minimal medium

AIMS: The aim was to develop reliable and economical protocols for the production of fully deuteriated biomolecules by bacteria. This required the preparation of deuterium-tolerant bacterial strains and an understanding of the physiological mechanisms of acquisition of deuterium tolerance. METHODS AND RESULTS: We report here improved methods for the cultivation of Escherichia coli on fully deuteriated minimal medium. A multi-stage adaptation protocol was developed; this included repeated plating and selection of colonies and resulted in highly deuterium-tolerant cell cultures. Three E. coli strains, JM109, MRE600 and MRE600Rif, were adapted to growth on deuteriated succinate medium. This is the first report of JM109 being adapted to deuteriated minimal media. The adapted strains showed good, consistent growth rates and were capable of being transformed with plasmids. Expression of heterologous proteins in these strains was reliable and yields were consistently high (100-200 mg l-1). We also show that all E. coli cells are inherently capable of growth on deuteriated media. CONCLUSIONS: We have developed a new adaptation protocol that resulted in three highly deuterium-tolerant E. coli strains. Deuterium-adapted cultures produced good yields of a deuteriated recombinant protein. We suggest that E. coli cells are inherently capable of growth on deuteriated media, but that non-specific mutations enhance deuterium tolerance. Thus plating and selection of colonies leads to highly deuterium-tolerant strains. SIGNIFICANCE AND IMPACT OF STUDY: An understanding of the mechanism of adaptation of E. coli to growth on deuteriated media allows strategies for the development of disabled deuterium-tolerant strains suitable for high-level production of deuteriated recombinant proteins and other biomolecules. This is of particular importance for nuclear magnetic resonance and neutron scattering studies of biomolecules.     

A High Frequency of Beneficial Mutations Across Multiple Fitness Components in Saccharomyces cerevisiae

Mutation-accumulation experiments are widely used to estimate parameters of spontaneous mutations affecting fitness. In many experiments only one component of fitness is measured. In a previous study involving the diploid yeast Saccharomyces cerevisiae, we measured the growth rate of 151 mutation-accumulation lines to estimate parameters of mutation. We found that an unexpectedly high frequency of fitness-altering mutations was beneficial. Here, we build upon our previous work by examining sporulation efficiency, spore viability, and haploid growth rate and find that these components of fitness also show a high frequency of beneficial mutations. We also examine whether mutation-acycumulation (MA) lines show any evidence of pleiotropy among accumulated mutations and find that, for most, there is none. However, MA lines that have zero fitness (i.e., lethality) for any one fitness component do show evidence for pleiotropy among accumulated mutations. We also report estimates of other parameters of mutation based on each component of fitness.ADAPTATION can occur from standing genetic variation or from newly arising mutations. The relative importance of these two sources of adaptive mutations is affected by a variety of factors, including those that alter standing levels of genetic variation (see Barrett and Schluter 2008) and those that generate new mutations. Predicting how quickly a population will adapt and the type of beneficial mutations that will fuel that adaptation requires estimates of the additive genetic variance in fitness and of the beneficial mutation rate and the distribution of beneficial effects. While additive genetic variance for fitness has been estimated in a variety of organisms (Mousseau and Roff 1987), the beneficial mutation rate and the distribution of beneficial effects have only been estimated in a few studies (Shaw et al. 2002; Joseph and Hall 2004; Perfeito et al. 2007; Dickinson 2008; Hall et al. 2008). Surprisingly, these studies estimate that between 6 (Joseph and Hall 2004) and 50% (Shaw et al. 2002) of fitness-altering mutations are beneficial. In contrast, most mutation-accumulation (MA) experiments identify few, if any, beneficial mutations. Such wildly different estimates have even been generated from studies of the same species in similar environments (Zeyl and Devisser 2001; Joseph and Hall 2004; Dickinson 2008; Hall et al. 2008). If these estimates are correct, then they would suggest that the genotypes used in these experiments have vastly different evolutionary potential with respect to their capacity to exhibit rapid adaptation from new mutations.A more likely scenario is that much of the variation in estimates of the beneficial mutation rate is due to methodological differences between studies. One possibility is the fitness component being analyzed. The beneficial mutation rate may be under- or overestimated if the fitness component is under stabilizing selection or subject to antagonistic pleiotropy. Analyses of mutation-accumulation data typically assume that selection is directional. As a result, analyses of phenotypes under stabilizing selection may falsely conclude that mutations that increase a phenotype are beneficial and mutations that lower values are deleterious (see Keightley and Lynch''s 2003 criticism of Shaw et al. 2002). Alternatively, the beneficial mutation rate may be over- (or under) estimated if mutations increase fitness in regard to one component, but lower fitness in regard to lifetime fitness or another fitness component (i.e., antagonistic pleiotropy). Here, we explore these possibilities by investigating whether the high beneficial mutation rates estimated from our previous experiments are specific to the fitness component that we examined.In two previous studies we accumulated mutations in 152 yeast, MA lines and used measures of their effects on diploid growth rate to estimate parameters of beneficial and deleterious mutations. In the first study we estimated that 6% of mutations accumulated during the first 1012 generations of accumulation improved diploid growth (Joseph and Hall 2004). To determine whether this high beneficial mutation rate was due to sampling error, we passaged the lines for an additional 1050 generations and found that 13% of mutations improved diploid growth (Hall et al. 2008). Similarly, another yeast MA experiment (Dickinson 2008) estimated an uncorrected frequency of beneficial mutations of 25%, although correction for within-colony selection reduces this estimate by approximately half. Together, these studies indicate that a substantial proportion of mutations accumulated in these yeast MA lines are beneficial for a single fitness component and that this observation cannot be explained by the chance sampling of a few beneficial mutations.In this study we return to our yeast MA lines (Joseph and Hall 2004) and examine whether the high beneficial mutation rate that we estimated after 1012 generations is an artifact of the fitness component that we examined. To test this hypothesis we examined whether our MA lines carry mutations that are beneficial across multiple fitness components: diploid growth, sporulation efficiency, spore viability, and haploid growth rate. If our previous results are due to us analyzing a fitness component that is either subject to stabilizing selection or antagonistic pleiotropy, then mutations accumulated in our lines will be conditionally beneficial and analyses of additional fitness components would yield different estimates of the beneficial mutation rate. We found that three of the four fitness components yield high estimates of the beneficial mutation rate. This suggests that multiple MA lines have accumulated beneficial mutations and that the high beneficial mutation rate that we previously estimated is not an artifact of the fitness component that we examined.Measuring multiple components of fitness also allowed us to examine the pleiotropic effects of beneficial and deleterious mutations. In general, we found that mutations altering one component of fitness have little effect on other components. However, lethal mutations were typically pleiotropic.

Conclusions:

We find that for three of four fitness components examined, a high frequency of spontaneous, fitness-altering mutations in diploid yeast is beneficial. Further, we do not detect pleiotropy of small-effect mutations, while lethal mutations show high levels of pleiotropy. In most cases, pleiotropy is positive. Two lines show evidence of antagonistic pleiotropy, indicating trade-offs, although heterozygote advantage cannot be ruled out.     

Evolutionary adaptation to freeze-thaw-growth cycles in Escherichia coli

Fifteen populations of Escherichia coli were propagated for 150 freeze-thaw-growth (FTG) cycles in order to study the phenotypic and genetic changes that evolve under these stressful conditions. Here we present the phenotypic differences between the evolved lines and their progenitors as measured by competition experiments and growth curves. Three FTG lines evolved from an ancestral strain that was previously used to start a long-term evolution experiment, while the other 12 FTG lines are derived from clones that had previously evolved for 20,000 generations at constant 37 degrees C. Competition experiments indicate that the former FTG group improved their mean fitness under the FTG regime by about 90% relative to their progenitor, while the latter FTG group gained on average about 60% relative to their own progenitors. These increases in fitness result from both improved survival during freezing and thawing and more rapid recovery to initiate exponential growth after thawing. This shorter lag phase is specific to recovery after freezing and thawing. Future work will seek to identify the mutations responsible for evolutionary adaptation to the FTG environment and use them to explore the physiological mechanisms that allow increased survival and more rapid recovery.     

Evolution of competitive fitness in experimental populations of E. coli: What makes one genotype a better competitor than another?

An important problem in microbial ecology is to identify those phenotypic attributes that are responsible for competitive fitness in a particular environment. Thousands of papers have been published on the physiology, biochemistry, and molecular genetics of Escherichia coli and other bacterial models. Nonetheless, little is known about what makes one genotype a better competitor than another even in such well studied systems. Here, we review experiments to identify the phenotypic bases of improved competitive fitness in twelve E. coli populations that evolved for thousands of generations in a defined environment, in which glucose was the limiting substrate. After 10000 generations, the average fitness of the derived genotypes had increased by 50% relative to the ancestor, based on competition experiments using marked strains in the same environment. The growth kinetics of the ancestral and derived genotypes showed that the latter have a shorter lag phase upon transfer into fresh medium and a higher maximum growth rate. Competition experiments were also performed in environments where other substrates were substituted for glucose. The derived genotypes are generally more fit in competition for those substrates that use the same mechanism of transport as glucose, which suggests that enhanced transport was an important target of natural selection in the evolutionary environment. All of the derived genotypes produce much larger cells than does the ancestor, even when both types are forced to grow at the same rate. Some, but not all, of the derived genotypes also have greatly elevated mutation rates. Efforts are now underway to identify the genetic changes that underlie those phenotypic changes, especially substrate specificity and elevated mutation rate, for which there are good candidate loci. Identification and subsequent manipulation of these genes may provide new insights into the reproducibility of adaptive evolution, the importance of co-adapted gene complexes, and the extent to which distinct phenotypes (e.g., substrate specificity and cell size) are affected by the same mutations.     

Big-benefit mutations in a bacteriophage inhibited with heat

High temperature inhibits the growth of the wild-type bacteriophage phiX174. Three different point mutations were identified that each recovered growth at high temperature. Two affected the major capsid protein (residues F188 and F242), and one affected the internal scaffolding protein (B114). One of the major capsid mutations (F242) is located in a beta strand that contacts B114 in the procapsid during viral maturation, whereas the other capsid mutation (F188) is involved in subunit interactions at the threefold axis of symmetry. Selective coefficients of these mutations ranged from 13.9 to 3.8 in the inhibitory, hot environment, but all mutations reduced fitness at normal temperature. The selective effect of one of the mutations (F242) was evaluated at high temperature in four different genetic backgrounds and exhibited epistasis of diminishing returns: as log fitness of the background genotype increased from -0.1 to 4.1, the fitness boost provided by the F242 mutation decreased from 3.9 to 0. 8. These results support a model in which viral fitness is bounded by an upper limit and the benefit of a mutation is scaled according to the remaining opportunity for fitness improvement in the genome.     

Adaptation of experimental yeast populations to stressful conditions in relation to population size

The purpose of this experiment was to find out how a population becomes adapted to extremely stressful conditions as its environment deteriorates. We created a deteriorating environment for experimental selection lines of yeast by a stepwise increase in the concentration of salt in the growth medium. After each step, we tested the ability of the lines to grow at a high concentration of salt near the lethal limit for the ancestral strain. We found that mutations enhancing growth in this highly stressful environment began to spread at intermediate salt concentrations. The degree of enhancement was related to effective population size by a power law with a small exponent. The effect size of these mutations also increased with the population size in a similar fashion. From these results, we interpret adaptation to lethal stress as an indirect response to selection for resistance to previous lower levels of stress in a deteriorating environment. This suggests that the pattern of genetic correlation between successively higher levels of stress is an important factor in facilitating evolutionary rescue.     

Polygenic Mutation in Drosophila Melanogaster: The Causal Relationship of Bristle Number to Fitness

The association between sternopleural and abdominal bristle number and fitness in Drosophila melanogaster was determined for sublines of an initially highly inbred strain that were maintained by divergent artificial selection for 150 generations or by random mating for 180 generations. Replicate selection lines had more extreme bristle numbers than those that were maintained without artificial selection at the same census size for approximately the same number of generations. The average fitness, estimated by a single generation of competition against a compound autosome strain, was 0.17 for lines selected for high and low abdominal bristle numbers and 0.19 for lines selected for high and low sternopleural bristle number. The average fitness of unselected lines, 0.46, was significantly higher than that of the selection lines. The fitnesses and the relationships of bristle number to fitness in progeny of all possible crosses of high X high (H X H), high X low (H X L) and low X low (L X L) selection lines were examined to determine whether the observed intermediate optima were caused by direct stabilizing selection on bristle number or by apparent stabilizing selection mediated through deleterious pleiotropic fitness effects of mutations affecting bristle number. Although bristle number was nearly additive for progeny of H X H, H X L and L X L crosses among sternopleural bristle selection lines, their mean fitnesses were not significantly different from each other, or from the mean fitness of the unselected lines, suggesting partly or completely recessive pleiotropic fitness effects cause apparent stabilizing selection. The average fitness of the progeny of H X H abdominal bristle selection lines was not significantly different from the fitness of unselected lines, but the mean fitness of the progeny of L X L crosses was not significantly different from that of the pure low lines. This is consistent with direct selection against low but not high abdominal bristle number, but the interpretation is confounded by variation in average degree of dominance for fitness (on average recessive in the high abdominal bristle selection lines and additive in the low abdominal bristle selection lines). Neither direct stabilizing selection nor pleiotropy, therefore, can account for all the observations.     

Measuring selection coefficients below 10(-3): method, questions, and prospects

Measuring fitness with precision is a key issue in evolutionary biology, particularly in studying mutations of small effects. It is usually thought that sampling error and drift prevent precise measurement of very small fitness effects. We circumvented these limits by using a new combined approach to measuring and analyzing fitness. We estimated the mutational fitness effect (MFE) of three independent mini-Tn10 transposon insertion mutations by conducting competition experiments in large populations of Escherichia coli under controlled laboratory conditions. Using flow cytometry to assess genotype frequencies from very large samples alleviated the problem of sampling error, while the effect of drift was controlled by using large populations and massive replication of fitness measures. Furthermore, with a set of four competition experiments between ancestral and mutant genotypes, we were able to decompose fitness measures into four estimated parameters that account for fitness effects of our fluorescent marker (α), the mutation (β), epistasis between the mutation and the marker (γ), and departure from transitivity (τ). Our method allowed us to estimate mean selection coefficients to a precision of 2 × 10(-4). We also found small, but significant, epistatic interactions between the allelic effects of mutations and markers and confirmed that fitness effects were transitive in most cases. Unexpectedly, we also detected variation in measures of s that were significantly bigger than expected due to drift alone, indicating the existence of cryptic variation, even in fully controlled experiments. Overall our results indicate that selection coefficients are best understood as being distributed, representing a limit on the precision with which selection can be measured, even under controlled laboratory conditions.     

The effect of spontaneous mutations on competitive ability

Understanding the impact of spontaneous mutations on fitness has many theoretical and practical applications in biology. Although mutational effects on individual morphological or life‐history characters have been measured in several classic genetic model systems, there are few estimates of the rate of decline due to mutation for complex fitness traits. Here, we estimate the effects of mutation on competitive ability, an important complex fitness trait, in a model system for ecological and evolutionary genomics, Daphnia. Competition assays were performed to compare fitness between mutation‐accumulation (MA) lines and control lines from eight different genotypes from two populations of Daphnia pulicaria after 30 and 65 generations of mutation accumulation. Our results show a fitness decline among MA lines relative to controls as expected, but highlight the influence of genomic background on this effect. In addition, in some assays, MA lines outperform controls providing insight into the frequency of beneficial mutations.     

Fluctuating asymmetry and fitness in Drosophila melanogaster

Models predict that developmental stability measured by fluctuating asymmetry should be positively correlated with fitness. Although such a correlation has often been suggested by indirect studies, there is still a lack of direct experimental evidence. In this note, I have measured the fluctuating asymmetry of sternopleural bristle counts in 32 lines of Drosophila melanogaster sharing the same genetic background but displaying all combinations of five visible mutations. Fluctuating asymmetry was heterogeneous among lines, suggesting a direct impact of the mutations on developmental stability. Two measures of fitness were made for each line: productivity (a combined measure of fecundity and egg‐to‐adult survivorship) and competitive male mating success. Fluctuating asymmetry was correlated with neither of these two components of fitness. This suggests that generalizations about fluctuating asymmetry must be taken with care.     

Compensatory adaptation to the deleterious effect of antibiotic resistance in Salmonella typhimurium

Most chromosomal mutations that cause antibiotic resistance impose fitness costs on the bacteria. This biological cost can often be reduced by compensatory mutations. In Salmonella typhimurium, the nucleotide substitution AAA42 --> AAC in the rpsL gene confers resistance to streptomycin. The resulting amino acid substitution (K42N) in ribosomal protein S12 causes an increased rate of ribosomal proofreading and, as a result, the rate of protein synthesis, bacterial growth and virulence are decreased. Eighty-one independent lineages of the low-fitness, K42N mutant were evolved in the absence of antibiotic to ameliorate the costs. From the rate of fixation of compensated mutants and their fitness, the rate of compensatory mutations was estimated to be > or = 10-7 per cell per generation. The size of the population bottleneck during evolution affected fitness of the adapted mutants: a larger bottleneck resulted in higher average fitness. Only four of the evolved lineages contained streptomycin-sensitive revertants. The remaining 77 lineages contained mutants that were still fully streptomycin resistant, had retained the original resistance mutation and also acquired compensatory mutations. Most of the compensatory mutations, resulting in at least 35 different amino acid substitutions, were novel single-nucleotide substitutions in the rpsD, rpsE, rpsL or rplS genes encoding the ribosomal proteins S4, S5, S12 and L19 respectively. Our results show that the deleterious effects of a resistance mutation can be compensated by an unexpected variety of mutations.     

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.