Selection Affects Genes Involved in Replication during Long-Term Evolution in Experimental Populations of the Bacteriophage φX174

Observing organisms that evolve in response to strong selection over very short time scales allows the determination of the molecular mechanisms underlying adaptation. Although dissecting these molecular mechanisms is expensive and time-consuming, general patterns can be detected from repeated experiments, illuminating the biological processes involved in evolutionary adaptation. The bacteriophage φX174 was grown for 50 days in replicate chemostats under two culture conditions: Escherichia coli C as host growing at 37°C and Salmonella typhimurium as host growing at 43.5°C. After 50 days, greater than 20 substitutions per chemostat had risen to detectable levels. Of the 97 substitutions, four occurred in all four chemostats, five arose in both culture conditions, eight arose in only the high temperature S. typhimurium chemostats, and seven arose only in the E. coli chemostats. The remaining substitutions were detected only in a single chemostat, however, almost half of these have been seen in other similar experiments. Our findings support previous studies that host recognition and capsid stability are two biological processes that are modified during adaptation to novel hosts and high temperature. Based upon the substitutions shared across both environments, it is apparent that genome replication and packaging are also affected during adaptation to the chemostat environment, rather than to temperature or host per se. This environment is characterized by a large number of phage and very few hosts, leading to competition among phage within the host. We conclude from these results that adaptation to a high density environment selects for changes in genome replication at both protein and DNA sequence levels.     

Evolutionary genomics of host adaptation in vesicular stomatitis virus

Populations experiencing similar selection pressures can sometimes diverge in the genetic architectures underlying evolved complex traits. We used RNA virus populations of large size and high mutation rate to study the impact of historical environment on genome evolution, thus increasing our ability to detect repeatable patterns in the evolution of genetic architecture. Experimental vesicular stomatitis virus populations were evolved on HeLa cells, on MDCK cells, or on alternating hosts. Turner and Elena (2000. Cost of host radiation in an RNA virus. Genetics. 156:1465-1470.) previously showed that virus populations evolved in single-host environments achieved high fitness on their selected hosts but failed to increase in fitness relative to their ancestor on the unselected host and that alternating-host-evolved populations had high fitness on both hosts. Here we determined the complete consensus sequence for each evolved population after 95 generations to gauge whether the parallel phenotypic changes were associated with parallel genomic changes. We also analyzed the patterns of allele substitutions to discern whether differences in fitness across hosts arose through true pleiotropy or the presence of not only a mutation that is beneficial in both hosts but also 1 or more mutations at other loci that are costly in the unselected environment (mutation accumulation [MA]). We found that ecological history may influence to what extent pleiotropy and MA contribute to fitness asymmetries across environments. We discuss the degree to which current genetic architecture is expected to constrain future evolution of complex traits, such as host use by RNA viruses.     

Cost of host radiation in an RNA virus

Although host radiation allows a parasite to expand its ecological niche, traits governing the infection of multiple host types can decrease fitness in the original or alternate host environments. Reasons for this reduction in fitness include slower replication due to added genetic material or modifications, fitness trade-offs across host environments, and weaker selection resulting from simultaneous adaptation to multiple habitats. We examined the consequences of host radiation using vesicular stomatitis virus (VSV) and mammalian host cells in tissue culture. Replicate populations of VSV were allowed to evolve for 100 generations on the original host (BHK cells), on either of two novel hosts (HeLa and MDCK cells), or in environments where the availability of novel hosts fluctuated in a predictable or random way. As expected, each experimental population showed a substantial fitness gain in its own environment, but those evolved on new hosts (constant or fluctuating) suffered reduced competitiveness on the original host. However, whereas evolution on one novel host negatively correlated with performance on the unselected novel host, adaptation in fluctuating environments led to fitness improvements in both novel habitats.     

Genetic basis of infectivity evolution in a bacteriophage

Antagonistic coevolution between hosts and parasites is probably ubiquitous. However, very little is known of the genetic changes associated with parasite infectivity evolution during adaptation to a coevolving host. We followed the phenotypic and genetic changes in a lytic virus population (bacteriophage; phage Φ2) that coevolved with its bacterial host, Pseudomonas fluorescens SBW25. First, we show the rapid evolution of numerous unique phage infectivity phenotypes, and that both phage host range and bacterial resistance to individual phage increased over coevolutionary time. Second, each of the distinct phage phenotypes in our study had a unique genotype, and molecular evolution did not act uniformly across the phage genome during coevolution. In particular, we detected numerous substitutions on the tail fibre gene, which is involved in the first step of the host-parasite interaction: host adsorption. None of the observed mutations could be directly linked with infection against a particular host, suggesting that the phenotypic effects of infectivity mutations are probably epistatic. However, phage genotypes with the broadest host ranges had the largest number of nonsynonymous amino acid changes on genes implicated in infectivity evolution. An understanding of the molecular genetics of phage infectivity has helped to explain the complex phenotypic coevolutionary dynamics in this system.     

Experimental evolution of viruses: Microviridae as a model system

φX174 was developed as a model system for experimental studies of evolution because of its small genome size and ease of cultivation. It has been used extensively to address statistical questions about the dynamics of adaptive evolution. Molecular changes seen during experimental evolution of φX174 under a variety of conditions were compiled from 10 experiments comprising 58 lineages, where whole genomes were sequenced. A total of 667 substitutions was seen. Parallel evolution was rampant, with over 50 per cent of substitutions occurring at sites with three or more events. Comparisons of experimentally evolved sites to variation seen among wild phage suggest that at least some of the adaptive mechanisms seen in the laboratory are relevant to adaptation in nature. Elucidation of these mechanisms is aided by the availability of capsid and pro-capsid structures for φX174 and builds on years of genetic studies of the phage life history.     

Multihost experimental evolution of a plant RNA virus reveals local adaptation and host-specific mutations

For multihost pathogens, adaptation to multiple hosts has important implications for both applied and basic research. At the applied level, it is one of the main factors determining the probability and the severity of emerging disease outbreaks. At the basic level, it is thought to be a key mechanism for the maintenance of genetic diversity both in host and pathogen species. Using Tobacco etch potyvirus (TEV) and four natural hosts, we have designed an evolution experiment whose strength and novelty are the use of complex multicellular host organism as hosts and a high level of replication of different evolutionary histories and lineages. A pattern of local adaptation, characterized by a higher infectivity and virulence on host(s) encountered during the experimental evolution was found. Local adaptation only had a cost in terms of performance on other hosts in some cases. We could not verify the existence of a cost for generalists, as expected to arise from antagonistic pleiotropy and other genetic mechanisms generating a fitness trade-off between hosts. This observation confirms that this classical theoretical prediction lacks empirical support. We discuss the reasons for this discrepancy between theory and experiment in the light of our results. The analysis of full genome consensus sequences of the evolved lineages established that all mutations shared between lineages were host specific. A low degree of parallel evolution was observed, possibly reflecting the various adaptive pathways available for TEV in each host. Altogether, these results reveal a strong adaptive potential of TEV to new hosts without severe evolutionary constraints.     

A high-quality genome of the dobsonfly Neoneuromus ignobilis reveals molecular convergences in aquatic insects

Neoneuromus ignobilis is an archaic holometabolous aquatic predatory insect. However, a lack of genomic resources hinders the use of whole genome sequencing to explore their genetic basis and molecular mechanisms for adaptive evolution. Here, we provided a high-contiguity, chromosome-level genome assembly of N. ignobilis using high coverage Nanopore and PacBio reads with the Hi-C technique. The final assembly is 480.67 MB in size, containing 12 telomere-ended pseudochromosomes with only 17 gaps. We compared 42 hexapod species genomes including six independent lineages comprising 11 aquatic insects, and found convergent expansions of long wavelength-sensitive and blue-sensitive opsins, thermal stress response TRP channels, and sulfotransferases in aquatic insects, which may be related to their aquatic adaptation. We also detected strong nonrandom signals of convergent amino acid substitutions in aquatic insects. Collectively, our comparative genomic analysis revealed the evidence of molecular convergences in aquatic insects during both gene family evolution and convergent amino acid substitutions.     

Profiles of adaptation in two similar viruses.

The related bacteriophages phiX174 and G4 were adapted to the inhibitory temperature of 44 degrees and monitored for nucleotide changes throughout the genome. Phage were evolved by serial transfer at low multiplicity of infection on rapidly dividing bacteria to select genotypes with the fastest rates of reproduction. Both phage showed overall greater fitness effects per substitution during the early stages of adaptation. The fitness of phiX174 improved from -0.7 to 5.6 doublings of phage concentration per generation. Five missense mutations were observed. The earliest two mutations accounted for 85% of the ultimate fitness gain. In contrast, G4 required adaptation to the intermediate temperature of 41.5 degrees before it could be maintained at 44 degrees. Its fitness at 44 degrees increased from -2.7 to 3.2, nearly the same net gain as in phiX174, but with three times the opportunity for adaptation. Seventeen mutations were observed in G4: 14 missense, 2 silent, and 1 intergenic. The first 3 missense substitutions accounted for over half the ultimate fitness increase. Although the expected pattern of periodic selective sweeps was the most common one for both phage, some mutations were lost after becoming frequent, and long-term polymorphism was observed. This study provides the greatest detail yet in combining fitness profiles with the underlying pattern of genetic changes, and the results support recent theories on the range of fitness effects of substitutions fixed during adaptation.     

Evolutionary reversals during viral adaptation to alternating hosts

Experimental adaptation of the bacteriophage phiX174 to a Salmonella host depressed its ability to grow on the traditional Escherichia host, whereas adaptation to Escherichia did not appreciably affect growth on Salmonella. Continued host switching consistently exhibited this pattern. Growth inhibition on Escherichia resulted from two to three substitutions in the major capsid gene. When these phages were forced to grow again on Escherichia, fitness recovery occurred predominantly by reversions at these same sites, rather than by second-site compensatory changes, the more frequently observed mechanism in most microbial systems. The affected residues lie on the virion surface and they alter attachment efficiency, yet they occur in a region distinct from a putative binding region previously identified from X-ray crystallography. These residues not only experienced high rates of evolution in our experiments, but also exhibited high levels of radical amino acid variation among phiX174 and its known relatives, consistent with a history of adaptation involving these sites.     

Sex and the evolution of intrahost competition in RNA virus phi6.

Sex allows beneficial mutations that occur in separate lineages to be fixed in the same genome. For this reason, the Fisher-Muller model predicts that adaptation to the environment is more rapid in a large sexual population than in an equally large asexual population. Sexual reproduction occurs in populations of the RNA virus phi6 when multiple bacteriophages coinfect the same host cell. Here, we tested the model''s predictions by determining whether sex favors more rapid adaptation of phi6 to a bacterial host, Pseudomonas phaseolicola. Replicate populations of phi6 were allowed to evolve in either the presence or absence of sex for 250 generations. All experimental populations showed a significant increase in fitness relative to the ancestor, but sex did not increase the rate of adaptation. Rather, we found that the sexual and asexual treatments also differ because intense intrahost competition between viruses occurs during coinfection. Results showed that the derived sexual viruses were selectively favored only when coinfection is common, indicating that within-host competition detracts from the ability of viruses to exploit the host. Thus, sex was not advantageous because the cost created by intrahost competition was too strong. Our findings indicate that high levels of coinfection exceed an optimum where sex may be beneficial to populations of phi6, and suggest that genetic conflicts can evolve in RNA viruses.     

Rate of novel host invasion affects adaptability of evolving RNA virus lineages

Although differing rates of environmental turnover should be consequential for the dynamics of adaptive change, this idea has been rarely examined outside of theory. In particular, the importance of RNA viruses in disease emergence warrants experiments testing how differing rates of novel host invasion may impact the ability of viruses to adaptively shift onto a novel host. To test whether the rate of environmental turnover influences adaptation, we experimentally evolved 144 Sindbis virus lineages in replicated tissue-culture environments, which transitioned from being dominated by a permissive host cell type to a novel host cell type. The rate at which the novel host ‘invaded’ the environment varied by treatment. The fitness (growth rate) of evolved virus populations was measured on each host type, and molecular substitutions were mapped via whole genome consensus sequencing. Results showed that virus populations more consistently reached high fitness levels on the novel host when the novel host ‘invaded’ the environment more gradually, and gradual invasion resulted in less variable genomic outcomes. Moreover, virus populations that experienced a rapid shift onto the novel host converged upon different genotypes than populations that experienced a gradual shift onto the novel host, suggesting a strong effect of historical contingency.     

Experimental adaptation of Salmonella typhimurium to mice

Experimental evolution is a powerful approach to study the dynamics and mechanisms of bacterial niche specialization. By serial passage in mice, we evolved 18 independent lineages of Salmonella typhimurium LT2 and examined the rate and extent of adaptation to a mainly reticuloendothelial host environment. Bacterial mutation rates and population sizes were varied by using wild-type and DNA repair-defective mutator (mutS) strains with normal and high mutation rates, respectively, and by varying the number of bacteria intraperitoneally injected into mice. After <200 generations of adaptation all lineages showed an increased fitness as measured by a faster growth rate in mice (selection coefficients 0.11-0.58). Using a generally applicable mathematical model we calculated the adaptive mutation rate for the wild-type bacterium to be >10(-6)/cell/generation, suggesting that the majority of adaptive mutations are not simple point mutations. For the mutator lineages, adaptation to mice was associated with a loss of fitness in secondary environments as seen by a reduced metabolic capability. During adaptation there was no indication that a high mutation rate was counterselected. These data show that S. typhimurium can rapidly and extensively increase its fitness in mice but this niche specialization is, at least in mutators, associated with a cost.     

Phylogeography of a parasitoid wasp (Diaeretiella rapae): no evidence of host-associated lineages

The exceptional diversity of insects is often attributed to the effects of specialized relationships between insects and their hosts. Parasite-host interactions are influenced by current natural selection and dispersal, in addition to historical effects that may include past selection, vicariance, and random genetic drift. Both current and historical events can lead to reduced fitness on some hosts. If trade-offs in fitness on alternate hosts are common, adaptation to one host can prevent adaptation to another, giving rise to genetic differentiation among host-associated lineages. Previous studies of Diaeretiella rapae (Hymenoptera: Aphidiidae), a parasitoid of aphids, have revealed additive genetic differences in performance between populations that parasitize different aphid host species. To determine whether D. rapae populations collected from different aphid hosts have diverged into genetically independent lineages, we constructed a haplotype network based on sequence variation in mitochondrial DNA (mtDNA). We used single strand conformation polymorphism (SSCP) analysis to examine 2041 base pairs of mtDNA and to identify nucleotide sequences of 42 unique SSCP haplotypes. We found no association between mtDNA haplotypes and host species in either the ancestral range (Europe, Mediterranean region, Middle East, Asia) or part of the introduced range (western North America). Haplotypes likely to be ancestral were geographically widespread and found on both hosts, suggesting that the ability to use both hosts evolved prior to the diversification of the mtDNA. Ongoing gene flow appears to prevent the formation of host races.     

The Genetics of Adaptation for Eight Microvirid Bacteriophages

Theories of adaptive molecular evolution have recently experienced significant expansion, and their predictions and assumptions have begun to be subjected to rigorous empirical testing. However, these theories focus largely on predicting the first event in adaptive evolution, the fixation of a single beneficial mutation. To address long-term adaptation it is necessary to include new assumptions, but empirical data are needed for guidance. To empirically characterize the general properties of adaptive walks, eight recently isolated relatives of the single-stranded DNA (ssDNA) bacteriophage φX174 (family Microviridae) were adapted to identical selective conditions. Three of the eight genotypes were adapted in replicate, for a total of 11 adaptive walks. We measured fitness improvement and identified the genetic changes underlying the observed adaptation. Nearly all phages were evolvable; nine of the 11 lineages showed a significant increase in fitness. However, fitness plateaued quickly, and adaptation was achieved through only three substitutions on average. Parallel evolution was rampant, both across replicates of the same genotype as well as across different genotypes, yet adaptation of replicates never proceeded through the exact same set of mutations. Despite this, final fitnesses did not vary significantly among replicates. Final fitnesses did vary significantly across genotypes but not across phylogenetic groupings of genotypes. A positive correlation was found between the number of substitutions in an adaptive walk and the magnitude of fitness improvement, but no correlation was found between starting and ending fitness. These results provide an empirical framework for future adaptation theory.     

Local adaptation, resistance, and virulence in a hemiparasitic plant-host plant interaction

Abstract.— Coevolution may lead to local adaptation of parasites to their sympatric hosts. Locally adapted parasites are, on average, more infectious to sympatric hosts than to allopatric hosts of the same species or their fitness on the sympatric hosts is superior to that on allopatric hosts. We tested local adaptation of a hemiparasitic plant, Rhinanthus serotinus (Scrophulariaceae), to its host plant, the grass Agrostis capillaris . Using a reciprocal cross-infection experiment, we exposed host plants from four sites to hemiparasites originating from the same four sites in a common environment. The parasites were equally able to establish haustorial connections to sympatric and allopatric hosts, and their performance was similar on both host types. Therefore, these results do not indicate local adaptation of the parasites to their sympatric hosts. However, the parasite populations differed in average biomass and number of flowers per plant and in their effect on host biomass. These results indicate that the virulence of the parasite varied among populations, suggesting genetic variation. Theoretical models suggest that local adaptation is likely to be detected if the host and the parasite have different evolutionary potentials, different migration rates, and the parasite is highly virulent. In the interaction between R. serotinus and A. capillaris all the theoretical prerequisites for local adaptation may not be fulfilled.     

PHENOTYPIC STOCHASTICITY PROTECTS LYTIC BACTERIOPHAGE POPULATIONS FROM EXTINCTION DURING THE BACTERIAL STATIONARY PHASE

It is generally thought that the adsorption rate of a bacteriophage correlates positively with fitness, but this view neglects that most phages rely only on exponentially growing bacteria for productive infections. Thus, phages must cope with the environmental stochasticity that is their hosts’ physiological state. If lysogeny is one alternative, it is unclear how strictly lytic phages can survive the host stationary phase. Three scenarios may explain their maintenance: (1) pseudolysogeny, (2) diversified, or (3) conservative bet hedging. To better understand how a strictly lytic phage survives the stationary phase of its host, and how phage adsorption rate impacts this survival, we challenged two strictly lytic phage λ, differing in their adsorption rates, with stationary phase Escherichia coli cells. Our results showed that, pseudolysogeny was not responsible for phage survival and that, contrary to our expectation, high adsorption rate was not more detrimental during stationary phase than low adsorption rate. Interestingly, this last observation was due to the presence of the “residual fraction” (phages exhibiting extremely low adsorption rates), protecting phage populations from extinction. Whether this cryptic phenotypic variation is an adaptation (diversified bet hedging) or merely reflecting unavoidable defects during protein synthesis remains an open question.     

Poxviruses Deploy Genomic Accordions to Adapt Rapidly against Host Antiviral Defenses

In contrast to RNA viruses, double-stranded DNA viruses have low mutation rates yet must still adapt rapidly in response to changing host defenses. To determine mechanisms of adaptation, we subjected the model poxvirus vaccinia to serial propagation in human cells, where its antihost factor K3L is maladapted against the antiviral protein kinase R (PKR). Viruses rapidly acquired higher fitness via recurrent K3L gene amplifications, incurring up to 7%-10% increases in genome size. These transient gene expansions were necessary and sufficient to counteract human PKR and facilitated the gain of an adaptive amino acid substitution in K3L that also defeats PKR. Subsequent reductions in gene amplifications offset the costs associated with larger genome size while retaining adaptive substitutions. Our discovery of viral "gene-accordions" explains how poxviruses can rapidly adapt to defeat different host defenses despite low mutation rates and reveals how classical Red Queen conflicts can progress through unrecognized intermediates. PAPERFLICK:     

Compensatory evolution in response to a novel RNA polymerase: orthologous replacement of a central network gene

A bacteriophage genome was forced to evolve a new system of regulation by replacing its RNA polymerase (RNAP) gene, a central component of the phage developmental pathway, with that of a relative. The experiment used the obligate lytic phage T7 and the RNAP gene of phage T3. T7 RNAP uses 17 phage promoters, which are responsible for all middle and late gene expression, DNA replication, and progeny maturation, but the enzyme has known physical contacts with only 2 other phage proteins. T3 RNAP was supplied in trans by the bacterial host to a T7 genome lacking its own RNAP gene and the phage population was continually propagated on naive bacteria throughout the adaptation. Evolution of the T3 RNAP gene was thereby prevented, and selection was for the evolution of regulatory signals throughout the phage genome. T3 RNAP transcribes from T7 promoters only at low levels, but a single mutation in the promoter confers high expression, providing a ready mechanism for reevolution of gene expression in this system. When selected for rapid growth, fitness of the engineered phage evolved from a low of 5 doublings/h to 33 doublings/h, close to the expected maximum of 37 doublings/h. However, the experiment was terminated before it could be determined accurately that fitness had reached an obvious plateau, and it is not known whether further adaptation could have resulted in complete recovery of fitness. More than 30 mutations were observed in the evolved genome, but changes were found in only 9 of the 16 promoters, and several coding changes occurred in genes with no known contacts with the RNAP. Surprisingly, the T7 genome adapted to T3 RNAP also maintained high fitness when using T7 RNAP, suggesting that the extreme incompatibility of T7 elements with T3 RNAP is not an invariant property of divergence in these expression systems.