Viral Receptor-Binding Protein Evolves New Function through Mutations That Cause Trimer Instability and Functional Heterogeneity.

When proteins evolve new activity, a concomitant decrease in stability is often observed because the mutations that confer new activity can destabilize the native fold. In the conventional model of protein evolution, reduced stability is considered a purely deleterious cost of molecular innovation because unstable proteins are prone to aggregation and are sensitive to environmental stressors. However, recent work has revealed that nonnative, often unstable protein conformations play an important role in mediating evolutionary transitions, raising the question of whether instability can itself potentiate the evolution of new activity. We explored this question in a bacteriophage receptor-binding protein during host-range evolution. We studied the properties of the receptor-binding protein of bacteriophage λ before and after host-range evolution and demonstrated that the evolved protein is relatively unstable and may exist in multiple conformations with unique receptor preferences. Through a combination of structural modeling and in vitro oligomeric state analysis, we found that the instability arises from mutations that interfere with trimer formation. This study raises the intriguing possibility that protein instability might play a previously unrecognized role in mediating host-range expansions in viruses.

Viral protein instability enhances host-range evolvability.

Viruses are highly evolvable, but what traits endow this property? The high mutation rates of viruses certainly play a role, but factors that act above the genetic code, like protein thermostability, are also expected to contribute. We studied how the thermostability of a model virus, bacteriophage λ, affects its ability to evolve to use a new receptor, a key evolutionary transition that can cause host-range evolution. Using directed evolution and synthetic biology techniques we generated a library of host-recognition protein variants with altered stabilities and then tested their capacity to evolve to use a new receptor. Variants fell within three stability classes: stable, unstable, and catastrophically unstable. The most evolvable were the two unstable variants, whereas seven of eight stable variants were significantly less evolvable, and the two catastrophically unstable variants could not grow. The slowly evolving stable variants were delayed because they required an additional destabilizing mutation. These results are particularly noteworthy because they contradict a widely supported contention that thermostabilizing mutations enhance evolvability of proteins by increasing mutational robustness. Our work suggests that the relationship between thermostability and evolvability is more complex than previously thought, provides evidence for a new molecular model of host-range expansion evolution, and identifies instability as a potential predictor of viral host-range evolution.

Bacterial resistance response and resource availability mediate viral coexistence.

Viruses that infect bacteria, known as bacteriophages or phages, are the most prevalent entities on Earth. Their genetic diversity in nature is well documented, and members of divergent lineages can be found sharing the same ecological niche. This viral diversity can be influenced by a number of factors, including productivity, spatial structuring of the environment, and host-range trade-offs. Rapid evolution is also known to promote diversity by buffering ecological systems from extinction. There is, however, little known about the impact of coevolution on the maintenance of viral diversity within a microbial community. To address this, we developed a 4 species experimental system where two bacterial hosts, a generalist and a specialist phage, coevolved in a spatially homogenous environment over time. We observed the persistence of both viruses if the resource availability was sufficiently high. This coexistence occurred in the absence of any detectable host-range trade-offs that are costly for generalists and thus known to promote viral diversity. However, the coexistence was lost if two bacteria were not permitted to evolve alongside the phages or if two phages coevolved with a single bacterial host. Our findings indicate that a host's resistance response in mixed-species communities plays a significant role in maintaining viral diversity in the environment.

Coevolutionary phage training leads to greater bacterial suppression and delays the evolution of phage resistance.

The evolution of antibiotic-resistant bacteria threatens to become the leading cause of worldwide mortality. This crisis has renewed interest in the practice of phage therapy. Yet, bacteria's capacity to evolve resistance may debilitate this therapy as well. To combat the evolution of phage resistance and improve treatment outcomes, many suggest leveraging phages' ability to counter resistance by evolving phages on target hosts before using them in therapy (phage training). We found that in vitro, λtrn, a phage trained for 28 d, suppressed bacteria ∼1,000-fold for three to eight times longer than its untrained ancestor. Prolonged suppression was due to a delay in the evolution of resistance caused by several factors. Mutations that confer resistance to λtrn are ∼100× less common, and while the target bacterium can evolve complete resistance to the untrained phage in a single step, multiple mutations are required to evolve complete resistance to λtrn. Mutations that confer resistance to λtrn are more costly than mutations for untrained phage resistance. Furthermore, when resistance does evolve, λtrn is better able to suppress these forms of resistance. One way that λtrn improved was through recombination with a gene in a defunct prophage in the host genome, which doubled phage fitness. This transfer of information from the host genome is an unexpected but highly efficient mode of training phage. Lastly, we found that many other independently trained λ phages were able to suppress bacterial populations, supporting the important role training could play during phage therapeutic development.

Canonical Host-Pathogen Trade-Offs Subverted by Mutations with Dual Benefits.

AbstractHost-parasite coevolution is expected to drive the evolution of genetic diversity because the traits used in arms races-namely, host range and parasite resistance-are hypothesized to trade off with traits used in resource competition. We therefore tested data for several trade-offs among 93 isolates of bacteriophage λ and 51 Escherichia coli genotypes that coevolved during a laboratory experiment. Surprisingly, we found multiple trade-ups (positive trait correlations) but little evidence of several canonical trade-offs. For example, some bacterial genotypes evaded a trade-off between phage resistance and absolute fitness, instead evolving simultaneous improvements in both traits. This was surprising because our experimental design was predicted to expose resistance-fitness trade-offs by culturing E. coli in a medium where the phage receptor, LamB, is also used for nutrient acquisition. On reflection, LamB mediates not one but many trade-offs, allowing for more complex trait interactions than just pairwise trade-offs. Here, we report that mathematical reasoning and laboratory data highlight how trade-ups should exist whenever an evolutionary system exhibits multiple interacting trade-offs. Does this mean that coevolution should not promote genetic diversity? No, quite the contrary. We deduce that whenever positive trait correlations are observed in multidimensional traits, other traits may trade off and so provide the right circumstances for diversity maintenance. Overall, this study reveals that there are predictive limits when data account only for pairwise trait correlations, and it argues that a wider range of circumstances than previously anticipated can promote genetic and species diversity.

Leapfrog dynamics in phage-bacteria coevolution revealed by joint analysis of cross-infection phenotypes and whole genome sequencing.

Viruses and their hosts can undergo coevolutionary arms races where hosts evolve increased resistance and viruses evolve counter-resistance. Given these arms race dynamics (ARD), both players are predicted to evolve along a single trajectory as more recently evolved genotypes replace their predecessors. By coupling phenotypic and genomic analyses of coevolving populations of bacteriophage λ and Escherichia coli, we find conflicting evidence for ARD. Virus-host infection phenotypes fit the ARD model, yet genomic analyses revealed fluctuating selection dynamics. Rather than coevolution unfolding along a single trajectory, cryptic genetic variation emerges and is maintained at low frequency for generations until it eventually supplants dominant lineages. These observations suggest a hybrid 'leapfrog' dynamic, revealing weaknesses in the predictive power of standard coevolutionary models. The findings shed light on the mechanisms that structure coevolving ecological networks and reveal the limits of using phenotypic or genomic data alone to differentiate coevolutionary dynamics.

Leaky resistance and the conditions for the existence of lytic bacteriophage.

In experimental cultures, when bacteria are mixed with lytic (virulent) bacteriophage, bacterial cells resistant to the phage commonly emerge and become the dominant population of bacteria. Following the ascent of resistant mutants, the densities of bacteria in these simple communities become limited by resources rather than the phage. Despite the evolution of resistant hosts, upon which the phage cannot replicate, the lytic phage population is most commonly maintained in an apparently stable state with the resistant bacteria. Several mechanisms have been put forward to account for this result. Here we report the results of population dynamic/evolution experiments with a virulent mutant of phage Lambda, λVIR, and Escherichia coli in serial transfer cultures. We show that, following the ascent of λVIR-resistant bacteria, λVIR is maintained in the majority of cases in maltose-limited minimal media and in all cases in nutrient-rich broth. Using mathematical models and experiments, we show that the dominant mechanism responsible for maintenance of λVIR in these resource-limited populations dominated by resistant E. coli is a high rate of either phenotypic or genetic transition from resistance to susceptibility-a hitherto undemonstrated mechanism we term "leaky resistance." We discuss the implications of leaky resistance to our understanding of the conditions for the maintenance of phage in populations of bacteria-their "existence conditions."

Coevolution drives the emergence of complex traits and promotes evolvability.

The evolution of complex organismal traits is obvious as a historical fact, but the underlying causes--including the role of natural selection--are contested. Gould argued that a random walk from a necessarily simple beginning would produce the appearance of increasing complexity over time. Others contend that selection, including coevolutionary arms races, can systematically push organisms toward more complex traits. Methodological challenges have largely precluded experimental tests of these hypotheses. Using the Avida platform for digital evolution, we show that coevolution of hosts and parasites greatly increases organismal complexity relative to that otherwise achieved. As parasites evolve to counter the rise of resistant hosts, parasite populations retain a genetic record of past coevolutionary states. As a consequence, hosts differentially escape by performing progressively more complex functions. We show that coevolution's unique feedback between host and parasite frequencies is a key process in the evolution of complexity. Strikingly, the hosts evolve genomes that are also more phenotypically evolvable, similar to the phenomenon of contingency loci observed in bacterial pathogens. Because coevolution is ubiquitous in nature, our results support a general model whereby antagonistic interactions and natural selection together favor both increased complexity and evolvability.

Host coevolution alters the adaptive landscape of a virus.

The origin of new and complex structures and functions is fundamental for shaping the diversity of life. Such key innovations are rare because they require multiple interacting changes. We sought to understand how the adaptive landscape led to an innovation whereby bacteriophage λ evolved the new ability to exploit a receptor, OmpF, on Escherichia coli cells. Previous work showed that this ability evolved repeatedly, despite requiring four mutations in one virus gene. Here, we examine how this innovation evolved by studying six intermediate genotypes of λ isolated during independent transitions to exploit OmpF and comparing them to their ancestor. All six intermediates showed large increases in their adsorption rates on the ancestral host. Improvements in adsorption were offset, in large part, by the evolution of host resistance, which occurred by reduced expression of LamB, the usual receptor for λ. As a consequence of host coevolution, the adaptive landscape of the virus changed such that selection favouring four of the six virus intermediates became stronger after the host evolved resistance, thereby accelerating virus populations along the path to using the new OmpF receptor. This dependency of viral fitness on host genotype thus shows an important role for coevolution in the origin of the new viral function.

Statistical structure of host-phage interactions.

Interactions between bacteria and the viruses that infect them (i.e., phages) have profound effects on biological processes, but despite their importance, little is known on the general structure of infection and resistance between most phages and bacteria. For example, are bacteria-phage communities characterized by complex patterns of overlapping exploitation networks, do they conform to a more ordered general pattern across all communities, or are they idiosyncratic and hard to predict from one ecosystem to the next? To answer these questions, we collect and present a detailed metaanalysis of 38 laboratory-verified studies of host-phage interactions representing almost 12,000 distinct experimental infection assays across a broad spectrum of taxa, habitat, and mode of selection. In so doing, we present evidence that currently available host-phage infection networks are statistically different from random networks and that they possess a characteristic nested structure. This nested structure is typified by the finding that hard to infect bacteria are infected by generalist phages (and not specialist phages) and that easy to infect bacteria are infected by generalist and specialist phages. Moreover, we find that currently available host-phage infection networks do not typically possess a modular structure. We explore possible underlying mechanisms and significance of the observed nested host-phage interaction structure. In addition, given that most of the available host-phage infection networks examined here are composed of taxa separated by short phylogenetic distances, we propose that the lack of modularity is a scale-dependent effect, and then, we describe experimental studies to test whether modular patterns exist at macroevolutionary scales.

More evolvable bacteriophages better suppress their host.

The number of multidrug-resistant strains of bacteria is increasing rapidly, while the number of new antibiotic discoveries has stagnated. This trend has caused a surge in interest in bacteriophages as anti-bacterial therapeutics, in part because there is near limitless diversity of phages to harness. While this diversity provides an opportunity, it also creates the dilemma of having to decide which criteria to use to select phages. Here we test whether a phage's ability to coevolve with its host (evolvability) should be considered and how this property compares to two previously proposed criteria: fast reproduction and thermostability. To do this, we compared the suppressiveness of three phages that vary by a single amino acid yet differ in these traits such that each strain maximized two of three characteristics. Our studies revealed that both evolvability and reproductive rate are independently important. The phage most able to suppress bacterial populations was the strain with high evolvability and reproductive rate, yet this phage was unstable. Phages varied due to differences in the types of resistance evolved against them and their ability to counteract resistance. When conditions were shifted to exaggerate the importance of thermostability, one of the stable phages was most suppressive in the short-term, but not over the long-term. Our results demonstrate the utility of biological therapeutics' capacities to evolve and adjust in action to resolve complications like resistance evolution. Furthermore, evolvability is a property that can be engineered into phage therapeutics to enhance their effectiveness.

Inferring strain-level mutational drivers of phage-bacteria interaction phenotypes.

The enormous diversity of bacteriophages and their bacterial hosts presents a significant challenge to predict which phages infect a focal set of bacteria. Infection is largely determined by complementary -and largely uncharacterized- genetics of adsorption, injection, and cell take-over. Here we present a machine learning (ML) approach to predict phage-bacteria interactions trained on genome sequences of and phenotypic interactions amongst 51 Escherichia coli strains and 45 phage λ strains that coevolved in laboratory conditions for 37 days. Leveraging multiple inference strategies and without a priori knowledge of driver mutations, this framework predicts both who infects whom and the quantitative levels of infections across a suite of 2,295 potential interactions. The most effective ML approach inferred interaction phenotypes from independent contributions from phage and bacteria mutations, predicting phage host range with 86% mean classification accuracy while reducing the relative error in the estimated strength of the infection phenotype by 40%. Further, transparent feature selection in the predictive model revealed 18 of 176 phage λ and 6 of 18 E. coli mutations that have a significant influence on the outcome of phage-bacteria interactions, corroborating sites previously known to affect phage λ infections, as well as identifying mutations in genes of unknown function not previously shown to influence bacterial resistance. While the genetic variation studied was limited to a focal, coevolved phage-bacteria system, the method's success at recapitulating strain-level infection outcomes provides a path forward towards developing strategies for inferring interactions in non-model systems, including those of therapeutic significance.

Design, optimization, and inference of biphasic decay of infectious virus particles.

Virus population dynamics are driven by counter-balancing forces of production and loss. Whereas viral production arises from complex interactions with susceptible hosts, the loss of infectious virus particles is often approximated as a first-order kinetic process. As such, experimental protocols to measure infectious virus loss are not typically designed to identify non-exponential decay processes. Here, we propose methods to evaluate if an experimental design is adequate to identify multiphasic virus particle decay and to optimize the sampling times of decay experiments, accounting for uncertainties in viral kinetics. First, we evaluate synthetic scenarios of biphasic decays, with varying decay rates and initial proportions of subpopulations. We show that robust inference of multiphasic decay is more likely when the faster decaying subpopulation predominates insofar as early samples are taken to resolve the faster decay rate. Moreover, design optimization involving non-equal spacing between observations increases the precision of estimation while reducing the number of samples. We then apply these methods to infer multiple decay rates associated with the decay of bacteriophage ('phage') Φ D 9 , an evolved isolate derived from phage Φ 21 . A pilot experiment confirmed that Φ D 9 decay is multiphasic, but was unable to resolve the rate or proportion of the fast decaying subpopulation(s). We then applied a Fisher information matrix-based design optimization method to propose non-equally spaced sampling times. Using this strategy, we were able to robustly estimate multiple decay rates and the size of the respective subpopulations. Notably, we conclude that the vast majority (94%) of the phage Φ D 9 population decays at a rate 16-fold higher than the slow decaying population. Altogether, these results provide both a rationale and a practical approach to quantitatively estimate heterogeneity in viral decay.

Bacteriophage Φ21's receptor-binding protein evolves new functions through destabilizing mutations that generate non-genetic phenotypic heterogeneity.

How viruses evolve to expand their host range is a major question with implications for predicting the next pandemic. Gain-of-function experiments have revealed that host-range expansions can occur through relatively few mutations in viral receptor-binding proteins, and the search for molecular mechanisms that explain such expansions is underway. Previous research on expansions of receptor use in bacteriophage λ has shown that mutations that destabilize λ's receptor-binding protein cause it to fold into new conformations that can utilize novel receptors but have weakened thermostability. These observations led us to hypothesize that other viruses may take similar paths to expand their host range. Here, we find support for our hypothesis by studying another virus, bacteriophage 21 (Φ21), which evolves to use two new host receptors within 2 weeks of laboratory evolution. By measuring the thermodynamic stability of Φ21 and its descendants, we show that as Φ21 evolves to use new receptors and expands its host range, it becomes less stable and produces viral particles that are genetically identical but vary in their thermostabilities. Next, we show that this non-genetic heterogeneity between particles is directly associated with receptor use innovation, as phage particles with more derived receptor-use capabilities are more unstable and decay faster. Lastly, by manipulating the expression of protein chaperones during Φ21 infection, we demonstrate that heterogeneity in receptor use of phage particles arises during protein folding. Altogether, our results provide support for the hypothesis that viruses can evolve new receptor-use tropisms through mutations that destabilize the receptor-binding protein and produce multiple protein conformers.

Competition-driven eco-evolutionary feedback reshapes bacteriophage lambda's fitness landscape and enables speciation.

A major challenge in evolutionary biology is explaining how populations navigate rugged fitness landscapes without getting trapped on local optima. One idea illustrated by adaptive dynamics theory is that as populations adapt, their newly enhanced capacities to exploit resources alter fitness payoffs and restructure the landscape in ways that promote speciation by opening new adaptive pathways. While there have been indirect tests of this theory, to our knowledge none have measured how fitness landscapes deform during adaptation, or test whether these shifts promote diversification. Here, we achieve this by studying bacteriophage [Formula: see text], a virus that readily speciates into co-existing receptor specialists under controlled laboratory conditions. We use a high-throughput gene editing-phenotyping technology to measure [Formula: see text]'s fitness landscape in the presence of different evolved-[Formula: see text] competitors and find that the fitness effects of individual mutations, and their epistatic interactions, depend on the competitor. Using these empirical data, we simulate [Formula: see text]'s evolution on an unchanging landscape and one that recapitulates how the landscape deforms during evolution. [Formula: see text] heterogeneity only evolves in the shifting landscape regime. This study provides a test of adaptive dynamics, and, more broadly, shows how fitness landscapes dynamically change during adaptation, potentiating phenomena like speciation by opening new adaptive pathways.

An intron endonuclease facilitates interference competition between coinfecting viruses.

Introns containing homing endonucleases are widespread in nature and have long been assumed to be selfish elements that provide no benefit to the host organism. These genetic elements are common in viruses, but whether they confer a selective advantage is unclear. In this work, we studied intron-encoded homing endonuclease gp210 in bacteriophage ΦPA3 and found that it contributes to viral competition by interfering with the replication of a coinfecting phage, ΦKZ. We show that gp210 targets a specific sequence in ΦKZ, which prevents the assembly of progeny viruses. This work demonstrates how a homing endonuclease can be deployed in interference competition among viruses and provide a relative fitness advantage. Given the ubiquity of homing endonucleases, this selective advantage likely has widespread evolutionary implications in diverse plasmid and viral competition as well as virus-host interactions.

The effects of competition and predation on diversification in a model adaptive radiation.

Much of life's diversity is thought to have arisen through successive rounds of adaptive radiation-the rapid diversification of a lineage into a range of ecologically and phenotypically distinct species. Both resource competition and predation have been suggested as mechanisms driving this process, although the former is better studied than the latter. Here we show experimentally how predation by a protist, Tetrahymena thermophila, affects diversification in a model adaptive radiation of the bacterial prey, Pseudomonas fluorescens. We estimate the frequency-dependent fitness functions of competing niche-specialist prey in the presence and absence of predation, and use these to test hypotheses about the extent (measured as the number of new genotypes) and rate of diversification. Competition and predation independently generated diversifying selection that we show is capable of driving prey diversification to similar extents but at different rates, diversification being markedly delayed in the presence of predators. The cause of this delay stems from weaker diversifying selection due to the reduction in prey density caused by predation. Our results suggest that predation may play an under-appreciated role in driving adaptive radiations.

Comparison of bacterial suppression by phage cocktails, dual-receptor generalists, and coevolutionarily trained phages.

The evolution and spread of antibiotic-resistant bacteria have renewed interest in phage therapy, the use of bacterial viruses (phages) to combat bacterial infections. The delivery of phages in cocktails where constituent phages target different modalities (e.g., receptors) may improve treatment outcomes by making it more difficult for bacteria to evolve resistance. However, the multipartite nature of cocktails may lead to unintended evolutionary and ecological outcomes. Here, we compare a 2-phage cocktail with a largely unconsidered group of phages: generalists that can infect through multiple, independent receptors. We find that λ phage generalists and cocktails that target the same receptors (LamB and OmpF) suppress Escherichia coli similarly for ~2 days. Yet, a "trained" generalist phage, which previously adapted to its host via 28 days of coevolution, demonstrated superior suppression. To understand why the trained generalist was more effective, we measured the resistance of bacteria against each of our phages. We find that, when bacteria were assailed by two phages in the cocktail, they evolved mutations in manXYZ, a host inner-membrane transporter that λ uses to move its DNA across the periplasmic space and into the cell for infection. This provided cross-resistance against the cocktail and untrained generalist. However, these mutations were ineffective at blocking the trained generalist because, through coevolutionary training, it evolved to bypass manXYZ resistance. The trained generalist's past experiences in training make it exceedingly difficult for bacteria to evolve resistance, further demonstrating the utility of coevolutionary phage training for improving the therapeutic properties of phages.

Rapid bacteria-phage coevolution drives the emergence of multiscale networks.

Interactions between species catalyze the evolution of multiscale ecological networks, including both nested and modular elements that regulate the function of diverse communities. One common assumption is that such complex pattern formation requires spatial isolation or long evolutionary timescales. We show that multiscale network structure can evolve rapidly under simple ecological conditions without spatial structure. In just 21 days of laboratory coevolution, Escherichia coli and bacteriophage Φ21 coevolve and diversify to form elaborate cross-infection networks. By measuring ~10,000 phage-bacteria infections and testing the genetic basis of interactions, we identify the mechanisms that create each component of the multiscale pattern. Our results demonstrate how multiscale networks evolve in parasite-host systems, illustrating Darwin's idea that simple adaptive processes can generate entangled banks of ecological interactions.

Upper Bound on the Mutational Burden Imposed by a CRISPR-Cas9 Gene-Drive Element.

CRISPR-Cas9 gene drives (CCGDs) are powerful tools for genetic control of wild populations, useful for eradication of disease vectors, conservation of endangered species and other applications. However, Cas9 alone and in a complex with gRNA can cause double-stranded DNA breaks at off-target sites, which could increase the mutational load and lead to loss of heterozygosity (LOH). These undesired effects raise potential concerns about the long-term evolutionary safety of CCGDs, but the magnitude of these effects is unknown. To estimate how the presence of a CCGD or a Cas9 alone in the genome affects the rates of LOH events and de novo mutations, we carried out a mutation accumulation experiment in yeast Saccharomyces cerevisiae. Despite its substantial statistical power, our experiment revealed no detectable effect of CCGD or Cas9 alone on the genome-wide rates of mutations or LOH events, suggesting that these rates are affected by less than 30%. Nevertheless, we found that Cas9 caused a slight but significant shift towards more interstitial and fewer terminal LOH events, and the CCGD caused a significant difference in the distribution of LOH events on Chromosome V. Taken together, our results show that these genetic elements impose a weak and likely localized additional mutational burden in the yeast model. Although the mutagenic effects of CCGDs need to be further evaluated in other systems, our results suggest that the effect of CCGDs on off-target mutation rates and genetic diversity may be acceptable.