Negative Epistasis and Evolvability in TEM-1 β-Lactamase--The Thin Line between an Enzyme's Conformational Freedom and Disorder

Epistasis is a key factor in evolution since it determines which combinations of mutations provide adaptive solutions and which mutational pathways toward these solutions are accessible by natural selection. There is growing evidence for the pervasiveness of sign epistasis--a complete reversion of mutational effects, particularly in protein evolution--yet its molecular basis remains poorly understood. We describe the structural basis of sign epistasis between G238S and R164S, two adaptive mutations in TEM-1 β-lactamase--an enzyme that endows antibiotics resistance. Separated by 10 Å, these mutations initiate two separate trajectories toward increased hydrolysis rates and resistance toward second and third-generation cephalosporins antibiotics. Both mutations allow the enzyme's active site to adopt alternative conformations and accommodate the new antibiotics. By solving the corresponding set of crystal structures, we found that R164S causes local disorder whereas G238S induces discrete conformations. When combined, the mutations in 238 and 164 induce local disorder whereby nonproductive conformations that perturb the enzyme's catalytic preorganization dominate. Specifically, Asn170 that coordinates the deacylating water molecule is misaligned, in both the free form and the inhibitor-bound double mutant. This local disorder is not restored by stabilizing global suppressor mutations and thus leads to an evolutionary cul-de-sac. Conformational dynamism therefore underlines the reshaping potential of protein's structures and functions but also limits protein evolvability because of the fragility of the interactions networks that maintain protein structures.

Role of pleiotropy during adaptation of TEM-1 β-lactamase to two novel antibiotics

Pleiotropy is a key feature of the genotype–phenotype map, and its form and extent have many evolutionary implications, including for the dynamics of adaptation and the evolution of specialization. Similarly, pleiotropic effects of antibiotic resistance mutations may affect the evolution of antibiotic resistance in the simultaneous or fluctuating presence of different antibiotics. Here, we study the role of pleiotropy during the in vitro adaptation of the enzyme TEM-1 β-lactamase to two novel antibiotics, cefotaxime (CTX) and ceftazidime (CAZ). We subject replicate lines for four rounds of evolution to selection with CTX and CAZ alone, and in their combined and fluctuating presence. Evolved alleles show positive correlated responses when selecting with single antibiotics. Nevertheless, pleiotropic constraints are apparent from the effects of single mutations and from selected alleles showing smaller correlated than direct responses and smaller responses after simultaneous and fluctuating selection with both than with single antibiotics. We speculate that these constraints result from structural changes in the oxyanion pocket surrounding the active site, where accommodation of CTX and the larger CAZ is balanced against their positioning with respect to the active site. Our findings suggest limited benefits from the combined or fluctuating application of these related cephalosporins for containing antibiotic resistance.

Empirical fitness landscapes and the predictability of evolution

The genotype-fitness map (that is, the fitness landscape) is a key determinant of evolution, yet it has mostly been used as a superficial metaphor because we know little about its structure. This is now changing, as real fitness landscapes are being analysed by constructing genotypes with all possible combinations of small sets of mutations observed in phylogenies or in evolution experiments. In turn, these first glimpses of empirical fitness landscapes inspire theoretical analyses of the predictability of evolution. Here, we review these recent empirical and theoretical developments, identify methodological issues and organizing principles, and discuss possibilities to develop more realistic fitness landscape models.

Heterosis is prevalent among domesticated but not wild strains of Saccharomyces cerevisiae

Crosses between inbred but unrelated individuals often result in an increased fitness of the progeny. This phenomenon is known as heterosis and has been reported for wild and domesticated populations of plants and animals. Analysis of heterosis is often hindered by the fact that the genetic relatedness between analyzed organisms is only approximately known. We studied a collection of Saccharomyces cerevisiae isolates from wild and human-created habitats whose genomes were sequenced and thus their relatedness was fully known. We reasoned that if these strains accumulated different deleterious mutations at an approximately constant rate, then heterosis should be most visible in F1 heterozygotes from the least related parents. We found that heterosis was substantial and positively correlated with sequence divergence, but only in domesticated strains. More than 80% of the heterozygous hybrids were more fit than expected from the mean of their homozygous parents, and approximately three-quarters of those exceeded even the fittest parent. Our results support the notion that domestication brings about relaxation of selection and accumulation of deleterious mutations. However, other factors may have contributed as well. In particular, the observed build-up of genetic load might be facilitated by a decrease, and not increase, in the rate of inbreeding.

Incidence of symbiotic dsRNA 'killer' viruses in wild and domesticated yeast

Viruses are found in almost all organisms and physical habitats. One interesting example is the yeast viral 'killer system'. The virus provides the host with a toxin directed against strains that do not carry it, while the yeast cell enables its propagation. Although yeast viruses are believed to be common, they have been actually described only for a limited number of yeast isolates. We surveyed 136 Saccharomyces cerevisiae and S. paradoxus strains of known origin and phylogenetic relatedness. Of these, 14 (c. 10%) were infected by killer viruses of one of the three types: K1, K2 or K28. As many as 34 strains (c. 25%) were not sensitive to at least one type of the killer toxin. In most cases, resistance did not disappear after attempts to cure the host strains from their viruses, suggesting that it was encoded in the host's genome. In terms of phylogeny, killer strains appear to be more related to each other than to nonkiller ones. No such tendency is observed for the phenotype of toxin resistance. Our results suggest that even if the killer toxins are not always present, they do play significant role in yeast ecology and evolution.

Predicting the evolution of antibiotic resistance

Mutations causing antibiotic resistance are often associated with a cost in the absence of antibiotics. Surprisingly, a new study found that bacteria adapting to increased temperature became resistant to rifampicin. By studying the consequences of the involved mutations in different conditions and genetic backgrounds, the authors illustrate how knowledge of two fundamental genetic properties, pleiotropy and epistasis, may help to predict the evolution of antibiotic resistance.

See research article http://www.biomedcentral.com/1471-2148/13/50

Quantifying the adaptive potential of an antibiotic resistance enzyme

For a quantitative understanding of the process of adaptation, we need to understand its "raw material," that is, the frequency and fitness effects of beneficial mutations. At present, most empirical evidence suggests an exponential distribution of fitness effects of beneficial mutations, as predicted for Gumbel-domain distributions by extreme value theory. Here, we study the distribution of mutation effects on cefotaxime (Ctx) resistance and fitness of 48 unique beneficial mutations in the bacterial enzyme TEM-1 β-lactamase, which were obtained by screening the products of random mutagenesis for increased Ctx resistance. Our contributions are threefold. First, based on the frequency of unique mutations among more than 300 sequenced isolates and correcting for mutation bias, we conservatively estimate that the total number of first-step mutations that increase Ctx resistance in this enzyme is 87 [95% CI 75-189], or 3.4% of all 2,583 possible base-pair substitutions. Of the 48 mutations, 10 are synonymous and the majority of the 38 non-synonymous mutations occur in the pocket surrounding the catalytic site. Second, we estimate the effects of the mutations on Ctx resistance by determining survival at various Ctx concentrations, and we derive their fitness effects by modeling reproduction and survival as a branching process. Third, we find that the distribution of both measures follows a Fréchet-type distribution characterized by a broad tail of a few exceptionally fit mutants. Such distributions have fundamental evolutionary implications, including an increased predictability of evolution, and may provide a partial explanation for recent observations of striking parallel evolution of antibiotic resistance.

Model and test in a fungus of the probability that beneficial mutations survive drift

Determining the probability of fixation of beneficial mutations is critically important for building predictive models of adaptive evolution. Despite considerable theoretical work, models of fixation probability have stood untested for nearly a century. However, recent advances in experimental and theoretical techniques permit the development of models with testable predictions. We developed a new model for the probability of surviving genetic drift, a major component of fixation probability, for novel beneficial mutations in the fungus Aspergillus nidulans, based on the life-history characteristics of its colony growth on a solid surface. We tested the model by measuring the probability of surviving drift in 11 adapted strains introduced into wild-type populations of different densities. We found that the probability of surviving drift increased with mutant invasion fitness, and decreased with wild-type density, as expected. The model accurately predicted the survival probability for the majority of mutants, yielding one of the first direct tests of the extinction probability of beneficial mutations.

Millifluidic droplet analyser for microbiology

We present a novel millifluidic droplet analyser (MDA) for precisely monitoring the dynamics of microbial populations over multiple generations in numerous (≥10(3)) aqueous emulsion droplets (~100 nL). As a first application, we measure the growth rate of a bacterial strain and determine the minimal inhibitory concentration (MIC) for the antibiotic cefotaxime by incubating bacteria in a fine gradient of antibiotic concentrations. The detection of cell activity is based on the automated detection of an epifluorescent signal that allows the monitoring of microbial populations up to a size of ~10(6) cells. We believe that this device is helpful for the study of population dynamic consequences of microbe-environment interactions and of individual cell differences. Moreover, the fluidic machine may improve clinical tests, as it simplifies, automates and miniaturizes the screening of numerous microbial populations that grow and evolve in compartments with a finely tuned composition.

The causes of epistasis

Since Bateson's discovery that genes can suppress the phenotypic effects of other genes, gene interactions-called epistasis-have been the topic of a vast research effort. Systems and developmental biologists study epistasis to understand the genotype-phenotype map, whereas evolutionary biologists recognize the fundamental importance of epistasis for evolution. Depending on its form, epistasis may lead to divergence and speciation, provide evolutionary benefits to sex and affect the robustness and evolvability of organisms. That epistasis can itself be shaped by evolution has only recently been realized. Here, we review the empirical pattern of epistasis, and some of the factors that may affect the form and extent of epistasis. Based on their divergent consequences, we distinguish between interactions with or without mean effect, and those affecting the magnitude of fitness effects or their sign. Empirical work has begun to quantify epistasis in multiple dimensions in the context of metabolic and fitness landscape models. We discuss possible proximate causes (such as protein function and metabolic networks) and ultimate factors (including mutation, recombination, and the importance of natural selection and genetic drift). We conclude that, in general, pleiotropy is an important prerequisite for epistasis, and that epistasis may evolve as an adaptive or intrinsic consequence of changes in genetic robustness and evolvability.

Evolutionary accessibility of mutational pathways

Functional effects of different mutations are known to combine to the total effect in highly nontrivial ways. For the trait under evolutionary selection (‘fitness’), measured values over all possible combinations of a set of mutations yield a fitness landscape that determines which mutational states can be reached from a given initial genotype. Understanding the accessibility properties of fitness landscapes is conceptually important in answering questions about the predictability and repeatability of evolutionary adaptation. Here we theoretically investigate accessibility of the globally optimal state on a wide variety of model landscapes, including landscapes with tunable ruggedness as well as neutral ‘holey’ landscapes. We define a mutational pathway to be accessible if it contains the minimal number of mutations required to reach the target genotype, and if fitness increases in each mutational step. Under this definition accessibility is high, in the sense that at least one accessible pathway exists with a substantial probability that approaches unity as the dimensionality of the fitness landscape (set by the number of mutational loci) becomes large. At the same time the number of alternative accessible pathways grows without bounds. We test the model predictions against an empirical 8-locus fitness landscape obtained for the filamentous fungus Aspergillus niger. By analyzing subgraphs of the full landscape containing different subsets of mutations, we are able to probe the mutational distance scale in the empirical data. The predicted effect of high accessibility is supported by the empirical data and is very robust, which we argue reflects the generic topology of sequence spaces. Together with the restrictive assumptions that lie in our definition of accessibility, this implies that the globally optimal configuration should be accessible to genome wide evolution, but the repeatability of evolutionary trajectories is limited owing to the presence of a large number of alternative mutational pathways.

Fitness landscapes describe the fitness of related genotypes in a given environment, and can be used to identify which mutational steps lead towards higher fitness under particular evolutionary scenarios. The structure of a fitness landscape results from the way mutations interact in determining fitness, and can be smooth when mutations have multiplicative effect or rugged when interactions are strong and of opposite sign. Little is known about the structure of real fitness landscapes. Here, we study the evolutionary accessibility of fitness landscapes by using various landscape models with tunable ruggedness, and compare the results with an empirical fitness landscape involving eight marker mutations in the fungus Aspergillus niger. We ask how many mutational pathways from a low-fitness to the globally optimal genotype are accessible by natural selection in the sense that each step increases fitness. We find that for all landscapes with lower than maximal ruggedness the number of accessible pathways increases with increases of the number of loci involved, despite decreases in the accessibility for each pathway individually. We also find that models with intermediate ruggedness describe the A. niger data best.

Initial mutations direct alternative pathways of protein evolution

Whether evolution is erratic due to random historical details, or is repeatedly directed along similar paths by certain constraints, remains unclear. Epistasis (i.e. non-additive interaction between mutations that affect fitness) is a mechanism that can contribute to both scenarios. Epistasis can constrain the type and order of selected mutations, but it can also make adaptive trajectories contingent upon the first random substitution. This effect is particularly strong under sign epistasis, when the sign of the fitness effects of a mutation depends on its genetic background. In the current study, we examine how epistatic interactions between mutations determine alternative evolutionary pathways, using in vitro evolution of the antibiotic resistance enzyme TEM-1 β-lactamase. First, we describe the diversity of adaptive pathways among replicate lines during evolution for resistance to a novel antibiotic (cefotaxime). Consistent with the prediction of epistatic constraints, most lines increased resistance by acquiring three mutations in a fixed order. However, a few lines deviated from this pattern. Next, to test whether negative interactions between alternative initial substitutions drive this divergence, alleles containing initial substitutions from the deviating lines were evolved under identical conditions. Indeed, these alternative initial substitutions consistently led to lower adaptive peaks, involving more and other substitutions than those observed in the common pathway. We found that a combination of decreased enzymatic activity and lower folding cooperativity underlies negative sign epistasis in the clash between key mutations in the common and deviating lines (Gly238Ser and Arg164Ser, respectively). Our results demonstrate that epistasis contributes to contingency in protein evolution by amplifying the selective consequences of random mutations.

A long-term goal of evolutionary biology is to understand the factors that govern the outcome of evolution. Epistasis (i.e. the situation in which the fitness effect of a mutation depends on its genetic background) is one such factor. Epistasis not only affects the dynamics of evolution, it may also direct its outcome by affecting the type and order of selected mutations. This effect is particularly strong under sign epistasis, which occurs when the sign of a mutation's fitness effect depends on its genetic background. Here, we demonstrate how epistasis causes divergence of mutational pathways of an antibiotic resistance enzyme, TEM-1 β-lactamase. First, we use in vitro mutagenesis followed by selection for cefotaxime resistance to demonstrate that alternative mutational pathways towards highly resistant variants exist in addition to the main pathway that was previously described. Next, to test whether negative interactions between alternative initial substitutions govern this diversification, we start identical evolution experiments with alleles containing initial substitutions from the deviating lines. These alleles consistently evolve to lower adaptive peaks and acquire different mutations than those in the main pathway. Our results demonstrate that sign epistasis between alternative initial substitutions may force evolution to follow different mutational pathways.

An experimental test of the independent action hypothesis in virus-insect pathosystems

The 'independent action hypothesis' (IAH) states that each pathogen individual has a non-zero probability of causing host death and that pathogen individuals act independently. IAH has not been rigorously tested. In this paper, we (i) develop a probabilistic framework for testing IAH and (ii) demonstrate that, in two out of the six virus-insect pathosystems tested, IAH is supported by the data. We first show that IAH inextricably links host survivorship to the number of infecting pathogen individuals, and develop a model to predict the frequency of single- and dual-genotype infections when a host is challenged with a mixture of two genotypes. Model predictions were tested using genetically marked, near-identical baculovirus genotypes, and insect larvae from three host species differing in susceptibility. Observations in early-instar larvae of two susceptible host species support IAH, but observations in late-instar larvae of susceptible host species and larvae of a less susceptible host species were not in agreement with IAH. Hence the model is experimentally supported only in pathosystems in which the host is highly susceptible. We provide, to our knowledge, the first qualitative experimental evidence that, in such pathosystems, the action of a single virion is sufficient to cause disease.

The social evolution of somatic fusion

The widespread potential for somatic fusion among different conspecific multicellular individuals suggests that such fusion is adaptive. However, because recognition of non-kin (allorecognition) usually leads to a rejection response, successful somatic fusion is limited to close kin. This is consistent with kin-selection theory, which predicts that the potential cost of fusion and the potential for somatic parasitism decrease with increasing relatedness. Paradoxically, however, Crozier found that, in the short term, positive-frequency-dependent selection eliminates the required genetic polymorphism at allorecognition loci. The 'Crozier paradox' may be solved if allorecognition is based on extrinsically balanced polymorphisms, for example at immune loci. Alternatively, the assumption of most models that self fusion is mutually beneficial is wrong. If fusion is on average harmful, selection will promote unconditional rejection. However, we propose that fusion within individuals is beneficial, selecting for the ability to fuse, but fusion between individuals on average costly, selecting for non-self recognition (rather than non-kin recognition). We discuss experimental data on fungi that are consistent with this hypothesis.

Low multiplicity of infection in vivo results in purifying selection against baculovirus deletion mutants

The in vivo fate of Autographa californica multiple nucleopolyhedrovirus deletion mutants originating from serial passage in cell culture was investigated by passaging a population enriched in these mutants in insect larvae. The infectivity of polyhedra and occlusion-derived virion content per polyhedron were restored within two passages in vivo. The frequency of occurrence of deletion mutants was determined by real-time PCR. The frequency of the non-homologous region origin (non-HR ori) of DNA replication was reduced to wild-type levels within two passages. The frequency of the polyhedrin gene did not increase and remained below wild-type levels. A low m.o.i. during the initial infection in insect larvae, causing strong purifying selection for autonomously replicating viruses, could explain these observations. The same virus population used in vivo was also passaged once at a different m.o.i. in cell culture. A similar effect (i.e. lower non-HR ori frequency) was observed at low m.o.i. only, indicating that m.o.i. was the key selective condition.

Effect of dispersal and nutrient availability on the competitive ability of toxin-producing yeast

The ecological role of interference competition through toxin production is not well understood. In particular, it is unclear under what conditions the benefits of toxic killing outweigh the metabolic costs involved. A killer advantage has been suggested to rely on local competitive interactions where the benefits of killing accrue to the toxin producer preferentially, but this notion has little empirical support. In addition, contrasting predictions exist about the effect of resource abundance on the benefits of toxin production; this benefit should either be highest when resources are abundant and metabolic costs are relatively low or when resources are scarce and toxic killing is a 'last resort strategy' to obtain nutrients. Here, we test these predictions for one aspect of competitive ability, that is, the ability of toxin producers to invade a population of sensitive non-producers from a low initial frequency. We use competition experiments between isogenic K1 toxin-producing and non-producing strains of Saccharomyces cerevisiae, where we manipulate dispersal under two extreme nutrient conditions: one environment with and the other without replenishment of nutrients. We find that toxin production is beneficial when dispersal is limited under both nutrient conditions, but only when resources are abundant these outweigh its cost and allow invasion of the producer.

Heterogeneous adaptive trajectories of small populations on complex fitness landscapes

Background

Small populations are thought to be adaptively handicapped, not only because they suffer more from deleterious mutations but also because they have limited access to new beneficial mutations, particularly those conferring large benefits.

Methodology/Principal Findings

Here, we test this widely held conjecture using both simulations and experiments with small and large bacterial populations evolving in either a simple or a complex nutrient environment. Consistent with expectations, we find that small populations are adaptively constrained in the simple environment; however, in the complex environment small populations not only follow more heterogeneous adaptive trajectories, but can also attain higher fitness than the large populations. Large populations are constrained to near deterministic fixation of rare large-benefit mutations. While such determinism speeds adaptation on the smooth adaptive landscape represented by the simple environment, it can limit the ability of large populations from effectively exploring the underlying topography of rugged adaptive landscapes characterized by complex environments.

Conclusions

Our results show that adaptive constraints often faced by small populations can be circumvented during evolution on rugged adaptive landscapes.

Spatial structure inhibits the rate of invasion of beneficial mutations in asexual populations

Populations in spatially structured environments may be divided into a number of (semi-) isolated subpopulations due to limited offspring dispersal. Limited dispersal and, as a consequence, local competition could slow down the invasion of fitter mutants, allowing the short-term coexistence of ancestral genotypes and mutants. We determined the rate of invasion of beneficial mutants of Escherichia coli, dispersed to different degrees in a spatially structured environment during 40 generations, experimentally and theoretically. Simulations as well as experimental data show a decrease in the rate of invasion with increasingly constrained dispersal. When a beneficial mutant invades from a single spot, competition with the ancestral genotype takes place only along the edges of the growing colony patch. As the colony grows, the fitness of the mutant will decrease due to a decrease in the mutant's fraction that effectively competes with the surrounding ancestor. Despite its inherently higher competitive ability, increased intragenotype competition prevents the beneficial mutant from rapidly taking over, causing short-term coexistence of superior and inferior genotypes.

Bacterial solutions to the problem of sex

Why have sex? It seems that genetic studies in the normally non-sexualE. coli may have the key.

The evolution of sex: empirical insights into the roles of epistasis and drift

Despite many years of theoretical and experimental work, the explanation for why sex is so common as a reproductive strategy continues to resist understanding. Recent empirical work has addressed key questions in this field, especially regarding rates of mutation accumulation in sexual and asexual organisms, and the roles of negative epistasis and drift as sources of adaptive constraint in asexually reproducing organisms. At the same time, new ideas about the evolution of sexual recombination are being tested, including intriguing suggestions of an important interplay between sex and genetic architecture, which indicate that sex and recombination could have affected their own evolution.

The effect of population structure on the adaptive radiation of microbial populations evolving in spatially structured environments

Spatial structure is thought to be an important factor influencing the emergence and maintenance of genetic diversity. Previous studies have demonstrated that environmental heterogeneity, provided by spatial structure, leads to adaptive radiation of populations. In the present study, we investigate not only the impact of environmental heterogeneity on adaptive radiation, but also of population fragmentation and niche construction. Replicate populations founded by a single genotype of Escherichia coli were allowed to evolve for 900 generations by serial transfer in either a homogeneous environment, or a spatially structured environment that was either kept intact or destroyed with each daily transfer. Only populations evolving in the structured environment with intact population structure diversified: clones are significantly divergent in sugar catabolism, and show frequency-dependent fitness interactions indicative of stable coexistence. These findings demonstrate an important role for population fragmentation, a consequence of population structure in spatially structured environments, on the diversification of populations.

Clonal interference and the periodic selection of new beneficial mutations in Escherichia coli

The conventional model of adaptation in asexual populations implies sequential fixation of new beneficial mutations via rare selective sweeps that purge all variation and preserve the clonal genotype. However, in large populations multiple beneficial mutations may co-occur, causing competition among them, a phenomenon called "clonal interference." Clonal interference is thus expected to lead to longer fixation times and larger fitness effects of mutations that ultimately become fixed, as well as to a genetically more diverse population. Here, we study the significance of clonal interference in populations consisting of mixtures of differently marked wild-type and mutator strains of Escherichia coli that adapt to a minimal-glucose environment for 400 generations. We monitored marker frequencies during evolution and measured the competitive fitness of random clones from each marker state after evolution. The results demonstrate the presence of multiple beneficial mutations in these populations and slower and more erratic invasion of mutants than expected by the conventional model, showing the signature of clonal interference. We found that a consequence of clonal interference is that fitness estimates derived from invasion trajectories were less than half the magnitude of direct estimates from competition experiments, thus revealing fundamental problems with this fitness measure. These results force a reevaluation of the conventional model of periodic selection for asexual microbes.

Insertion-sequence-mediated mutations isolated during adaptation to growth and starvation in Lactococcus lactis

We studied the activity of three multicopy insertion sequence (IS) elements in 12 populations of Lactococcus lactis IL1403 that evolved in the laboratory for 1000 generations under various environmental conditions (growth or starvation and shaken or stationary). Using RFLP analysis of single-clone representatives of each population, nine IS-mediated mutations were detected across all environmental conditions and all involving IS981. When it was assumed that these mutations were neutral, their frequency was higher under shaken than under stationary conditions, possibly due to oxygen stress. We characterized seven of the nine mutations at the molecular level and studied their population dynamics where possible. Two were simple insertions into new positions and the other five were recombinational deletions (of <1->10 kb) among existing and new copies of IS981; in all but one case these mutations disrupted gene functions. The best candidate beneficial mutations were two deletions of which similar versions were detected in two populations each. One of these two parallel deletions, affecting a gene involved in bacteriophage resistance, showed intermediate rearrangements and may also have resulted from increased local transposition rates.