EVOLUTIONARY DYNAMICS OF FITNESS RECOVERY FROM THE DEBILITATING EFFECTS OF MULLER'S RATCHET

The great adaptability shown by RNA viruses is a consequence of their high mutation rates. The evolution of fitness in a severely debilitated, clonal population of the nonsegmented ribovirus vesicular stomatitis virus (VSV) has been compared under five different demographic regimes, ranging from severe serial bottleneck passages (one virion) to large population passages (10⁵ virions or more) under similar environmental conditions (cell culture type and temperature). No matter how small the bottleneck, the fitness of the evolved populations was always higher than the fitness of the starting population; this result is clearly different from that previously reported for viruses with higher fitness. The reattainment of fitness under a regime of serial population passages showed two main characteristics: (1) the rate of adaptation was higher during early passages; and (2) a maximum fitness value was reached after a large number of passages. The maximum fitness reached by this initially debilitated clone was similar to the fitness of wild-type virus. The practical implications of these findings in the design of vaccines using attenuated viruses are also discussed.     

Genetic and fitness changes accompanying adaptation of an arbovirus to vertebrate and invertebrate cells

The alternating host cycle and persistent vector infection may constrain the evolution of arboviruses. To test this hypothesis, eastern equine encephalitis virus was passaged in BHK or mosquito cells, as well as in alternating (both) host cell passages. High and low multiplicities were used to examine the effect of defective interfering particles. Clonal BHK and persistent mosquito cell infections were also evaluated. Fitness was measured with one-step growth curves and competition assays, and mutations were evaluated by nucleotide sequencing and RNA fingerprinting. All passages and assays were done at 32 degrees C to eliminate temperature as a selection factor. Viruses passaged in either cell type alone exhibited fitness declines in the bypassed cells, while high-multiplicity and clonal passages caused fitness declines in both types of cells. Bypassed cell fitness losses were mosquito and vertebrate specific and were not restricted to individual cell lines. Fitness increases occurred in the cell line used for single-host-adaptation passages and in both cells for alternately passaged viruses. Surprisingly, single-host-cell passage increased fitness in that cell type no more than alternating passages. However, single-host-cell adaptation resulted in more mutations than alternating cell passages. Mosquito cell adaptation invariably resulted in replacement of the stop codon in nsP3 with arginine or cysteine. In one case, BHK cell adaptation resulted in a 238-nucleotide deletion in the 3' untranslated region. Many nonsynonymous substitutions were shared among more than one BHK or mosquito cell passage series, suggesting positive Darwinian selection. Our results suggest that alternating host transmission cycles constrain the evolutionary rates of arboviruses but not their fitness for either host alone.     

The role of evolutionary intermediates in the host adaptation of canine parvovirus

The adaptation of viruses to new hosts is a poorly understood process likely involving a variety of viral structures and functions that allow efficient replication and spread. Canine parvovirus (CPV) emerged in the late 1970s as a host-range variant of a virus related to feline panleukopenia virus (FPV). Within a few years of its emergence in dogs, there was a worldwide replacement of the initial virus strain (CPV type 2) by a variant (CPV type 2a) characterized by four amino acid differences in the capsid protein. However, the evolutionary processes that underlie the acquisition of these four mutations, as well as their effects on viral fitness, both singly and in combination, are still uncertain. Using a comprehensive experimental analysis of multiple intermediate mutational combinations, we show that these four capsid mutations act in concert to alter antigenicity, cell receptor binding, and relative in vitro growth in feline cells. Hence, host adaptation involved complex interactions among both surface-exposed and buried capsid mutations that together altered cell infection and immune escape properties of the viruses. Notably, most intermediate viral genotypes containing different combinations of the four key amino acids possessed markedly lower fitness than the wild-type viruses.     

Many-trillionfold amplification of single RNA virus particles fails to overcome the Muller's ratchet effect.

We showed earlier that transfers of large populations of RNA viruses lead to fitness gains and that repeated genetic bottleneck transfers result in fitness losses due to Muller's ratchet. In the present study, we examined the effects of genetic bottleneck passages intervening between population passages, a process akin to some natural viral transmissions, using vesicular stomatitis virus as a model. Our findings show that the pronounced fitness increases that occur during two successive population passages cannot overcome the fitness decreases caused by a single intervening genetic bottleneck passage. The implications for natural transmissions of RNA viruses are discussed.     

Effect of alternating passage on adaptation of sindbis virus to vertebrate and invertebrate cells

Mosquito-borne alphaviruses, which replicate alternately and obligately in mosquitoes and vertebrates, appear to experience lower rates of evolution than do many RNA viruses that replicate solely in vertebrates. This genetic stability is hypothesized to result from the alternating host cycle, which constrains evolution by imposing compromise fitness solutions in each host. To test this hypothesis, Sindbis virus was passaged serially, either in one cell type to eliminate host alteration or alternately between vertebrate (BHK) and mosquito (C6/36) cells. Following 20 to 50 serial passages, mutations were identified and changes in fitness were assessed using competition assays against genetically marked, surrogate parent viruses. Specialized viruses passaged in a single cell exhibited more mutations and amino acid changes per passage than those passaged alternately. Single host-adapted viruses exhibited fitness gains in the cells in which they specialized but fitness losses in the bypassed cell type. Most but not all viruses passaged alternately experienced lesser fitness gains than specialized viruses, with fewer mutations per passage. Clonal populations derived from alternately passaged viruses also exhibited adaptation to both cell lines, indicating that polymorphic populations are not required for simultaneous fitness gains in vertebrate and mosquito cells. Nearly all passaged viruses acquired Arg or Lys substitutions in the E2 envelope glycoprotein, but enhanced binding was only detected for BHK cells. These results support the hypothesis that arbovirus evolution may be constrained by alternating host transmission cycles, but they indicate a surprising ability for simultaneous adaptation to highly divergent cell types by combinations of mutations in single genomes.     

Size of genetic bottlenecks leading to virus fitness loss is determined by mean initial population fitness.

Genetic bottlenecks are important events in the genetic diversification of organisms and colonization of new ecological niches. Repeated bottlenecking of RNA viruses often leads to fitness losses due to the operation of Muller's ratchet. Herein we use vesicular stomatitis virus to determine the transmission population size which leads to fitness decreases of virus populations. Remarkably, the effective size of a genetic bottleneck associated with fitness loss is greater when the fitness of the parental population increases. For example, for starting virus populations with low fitness, population transfers of five-clone-to-five-clone passages resulted in a fitness increase. However, when a parental population with high fitness was transferred, 30-clone-to-30-clone passages were required simply to maintain fitness values.     

Evidence of expanded host range and mammalian-associated genetic changes in a duck H9N2 influenza virus following adaptation in quail and chickens

H9N2 avian influenza viruses continue to circulate worldwide; in Asia, H9N2 viruses have caused disease outbreaks and established lineages in land-based poultry. Some H9N2 strains are considered potentially pandemic because they have infected humans causing mild respiratory disease. In addition, some of these H9N2 strains replicate efficiently in mice without prior adaptation suggesting that H9N2 strains are expanding their host range. In order to understand the molecular basis of the interspecies transmission of H9N2 viruses, we adapted in the laboratory a wildtype duck H9N2 virus, influenza A/duck/Hong Kong/702/79 (WT702) virus, in quail and chickens through serial lung passages. We carried out comparative analysis of the replication and transmission in quail and chickens of WT702 and the viruses obtained after 23 serial passages in quail (QA23) followed by 10 serial passages in chickens (QA23CkA10). Although the WT702 virus can replicate and transmit in quail, it replicates poorly and does not transmit in chickens. In contrast, the QA23CkA10 virus was very efficient at replicating and transmitting in quail and chickens. Nucleotide sequence analysis of the QA23 and QA23CkA10 viruses compared to the WT702 virus indicated several nucleotide substitutions resulting in amino acid changes within the surface and internal proteins. In addition, a 21-amino acid deletion was found in the stalk of the NA protein of the QA23 virus and was maintained without further modification in the QA23CkA10 adapted virus. More importantly, both the QA23 and the QA23CkA10 viruses, unlike the WT702 virus, were able to readily infect mice, produce a large-plaque phenotype, showed faster replication kinetics in tissue culture, and resulted in the quick selection of the K627 amino acid mammalian-associated signature in PB2. These results are in agreement with the notion that adaptation of H9 viruses to land-based birds can lead to strains with expanded host range.     

Genetic bottlenecks and population passages cause profound fitness differences in RNA viruses.

Repeated clone-to-clone (genetic bottleneck) passages of an RNA phage and vesicular stomatitis virus have been shown previously to result in loss of fitness due to Muller's ratchet. We now demonstrate that Muller's ratchet also operates when genetic bottleneck passages are carried out at 37 rather than 32 degrees C. Thus, these fitness losses do not depend on growth of temperature-sensitive (ts) mutants at lowered temperatures. We also demonstrate that during repeated genetic bottleneck passages, accumulation of deleterious mutations does occur in a stepwise (ratchet-like) manner as originally proposed by Muller. One selected clone which had undergone significant loss of fitness after only 20 genetic bottleneck passages was passaged again in clone-to-clone series. Additional large losses of fitness were observed in five of nine independent bottleneck series; the relative fitnesses of the other four series remained close to the starting fitness. In sharp contrast, when the same selected clone was transferred 20 more times as large populations (10(5) to 10(6) PFU transferred at each passage), significant increases in fitness were observed in all eight passage series. Finally, we selected several clones which had undergone extreme losses of fitness during 20 bottleneck passages. When these low-fitness clones were passaged many times as large virus populations, they always regained very high relative fitness. We conclude that transfer of large populations of RNA viruses regularly selects those genomes within the quasispecies population which have the highest relative fitness, whereas bottleneck transfers have a high probability of leading to loss of fitness by random isolation of genomes carrying debilitating mutations. Both phenomena arise from, and underscore, the extreme mutability and variability of RNA viruses.     

Impact of chicken-origin cells on adaptation of a low pathogenic influenza virus

Understanding the growth dynamics of influenza viruses is an essential step in virus replication and cell-adaptation. The aim of this study was to elucidate the growth kinetic of a low pathogenic avian influenza H9N2 subtype in chicken embryo fibroblast (CEF) and chicken tracheal epithelial (CTE) cells during consecutive passages. An egg-adapted H9N2 virus was seeded into both cell culture systems. The amount of infectious virus released into the cell culture supernatants at interval times post-infection were titered and plaque assayed. The results as well as cell viability results indicate that the infectivity of the influenza virus was different among these primary cells. The egg-adapted H9N2 virus featured higher infectivity in CTE than in CEF cells. After serial passages and plaque purifications of the virus, a CTE cell-adapted strain was generated which carried amino acid substitutions within the HA stem region. The strain showed faster replication kinetics in cell culture resulting in an increase in virus titer. Overall, the present study provides the impact of cell type, multiplicity of infection, cellular protease roles in virus infectivity and finally molecular characterization during H9N2 virus adaptation procedure.     

Vertical Transmission Selects for Reduced Virulence in a Plant Virus and for Increased Resistance in the Host

For the last three decades, evolutionary biologists have sought to understand which factors modulate the evolution of parasite virulence. Although theory has identified several of these modulators, their effect has seldom been analysed experimentally. We investigated the role of two such major factors—the mode of transmission, and host adaptation in response to parasite evolution—in the evolution of virulence of the plant virus Cucumber mosaic virus (CMV) in its natural host Arabidopsis thaliana. To do so, we serially passaged three CMV strains under strict vertical and strict horizontal transmission, alternating both modes of transmission. We quantified seed (vertical) transmission rate, virus accumulation, effect on plant growth and virulence of evolved and non-evolved viruses in the original plants and in plants derived after five passages of vertical transmission. Our results indicated that vertical passaging led to adaptation of the virus to greater vertical transmission, which was associated with reductions of virus accumulation and virulence. On the other hand, horizontal serial passages did not significantly modify virus accumulation and virulence. The observed increases in CMV seed transmission, and reductions in virus accumulation and virulence in vertically passaged viruses were due also to reciprocal host adaptation during vertical passages, which additionally reduced virulence and multiplication of vertically passaged viruses. This result is consistent with plant-virus co-evolution. Host adaptation to vertically passaged viruses was traded-off against reduced resistance to the non-evolved viruses. Thus, we provide evidence of the key role that the interplay between mode of transmission and host-parasite co-evolution has in determining the evolution of virulence.     

Genetic and phenotypic variation of foot-and-mouth disease virus during serial passages in a natural host

Foot-and-mouth disease virus (FMDV), like other RNA viruses, exhibits high mutation rates during replication that have been suggested to be of adaptive value. However, even though genetic variation in RNA viruses and, more specifically, FMDV has been extensively examined during virus replication in a wide variety of in vitro cell cultures, very little is known regarding the generation and effects of genetic variability of virus replication in the natural host under experimental conditions and no genetic data are available regarding the effects of serial passage in natural hosts. Here, we present the results of 20 serial contact transmissions of the highly pathogenic, pig-adapted O Taiwan 97 (O Tw97) isolate of FMDV in swine. We examined the virus genomic consensus sequences for a total of 37 full-length viral genomes recovered from 20 in vivo passages. The characteristics and distributions of changes in the sequences during the series of pig infections were analyzed in comparison to the O Tw97 genomes recovered from serially infected BHK-21 cell cultures. Unexpectedly, a significant reduction of virulence upon pig passages was observed, and finally, interruption of the viral transmission chain occurred after the14th pig passage (T14). Virus was, however, isolated from the tonsils and nasal swabs of the asymptomatic T15 pigs at 26 days postcontact, consistent with a natural establishment of the carrier state previously described only for ruminants. Surprisingly, the region encoding the capsid protein VP1 (1D) did not show amino acid changes during in vivo passages. These data demonstrate that contact transmission of FMDV O Tw97 in pigs mimics the fitness loss induced by the bottleneck effect, which was previously observed by others during plaque-to-plaque FMDV passage in vitro, suggesting that unknown mechanisms of virulence recovery might be necessary during the evolution and perpetuation of FMDV in nature.     

The impact of high‐order epistasis in the within‐host fitness of a positive‐sense plant RNA virus

RNA viruses are the main source of emerging infectious diseases because of the evolutionary potential bestowed by their fast replication, large population sizes and high mutation and recombination rates. However, an equally important property, which is usually neglected, is the topography of the fitness landscape. How many fitness maxima exist and how well they are connected is especially interesting, as this determines the number of accessible evolutionary pathways. To address this question, we have reconstructed a region of the fitness landscape of tobacco etch potyvirus constituted by mutations observed during the experimental adaptation of the virus to the novel host Arabidopsis thaliana. Fitness was measured for many genotypes and showed the existence of multiple peaks and holes in the landscape. We found prevailing epistatic effects between mutations, with cases of reciprocal sign epistasis being common among pairs of mutations. We also found that high‐order epistasis was as important as pairwise epistasis in their contribution to fitness. Therefore, results suggest that the landscape was rugged due to the existence of holes caused by lethal genotypes, that a very limited number of potential neutral paths exist and that it contained a single adaptive peak.     

Dengue Virus RNA Structure Specialization Facilitates Host Adaptation

Many viral pathogens cycle between humans and insects. These viruses must have evolved strategies for rapid adaptation to different host environments. However, the mechanistic basis for the adaptation process remains poorly understood. To study the mosquito-human adaptation cycle, we examined changes in RNA structures of the dengue virus genome during host adaptation. Deep sequencing and RNA structure analysis, together with fitness evaluation, revealed a process of host specialization of RNA elements of the viral 3’UTR. Adaptation to mosquito or mammalian cells involved selection of different viral populations harvesting mutations in a single stem-loop structure. The host specialization of the identified RNA structure resulted in a significant viral fitness cost in the non-specialized host, posing a constraint during host switching. Sequence conservation analysis indicated that the identified host adaptable stem loop structure is duplicated in dengue and other mosquito-borne viruses. Interestingly, functional studies using recombinant viruses with single or double stem loops revealed that duplication of the RNA structure allows the virus to accommodate mutations beneficial in one host and deleterious in the other. Our findings reveal new concepts in adaptation of RNA viruses, in which host specialization of RNA structures results in high fitness in the adapted host, while RNA duplication confers robustness during host switching.     

Unusual distribution of mutations associated with serial bottleneck passages of human immunodeficiency virus type 1

Repeated bottleneck passages result in fitness losses of RNA viruses. In the case of human immunodeficiency virus type 1 (HIV-1), decreases in fitness after a limited number of plaque-to-plaque transfers in MT-4 cells were very drastic. Here we report an analysis of entire genomic nucleotide sequences of four HIV-1 clones derived from the same HIV-1 isolate and their low-fitness progeny following 7 to 15 plaque-to-plaque passages. Clones accumulated 4 to 28 mutations per genome, with dominance of A --> G and G --> A transitions (57% of all mutations) and 49% nonsynonymous replacements. One clone-but not three sibling clones-showed an overabundance of G --> A transitions, evidencing the highly stochastic nature of some types of mutational bias. The distribution of mutations along the genome was very unusual in that mutation frequencies in gag were threefold higher than in env. Particularly striking was the complete absence of replacements in the V3 loop of gp120, confirmed with partial nucleotide sequences of additional HIV-1 clones subjected to repeated bottleneck passages. The analyses revealed several amino acid replacements that have not been previously recorded among natural HIV-1 isolates and illustrate how evolution of an RNA virus genome, with regard to constant and variable regions, can be profoundly modified by alterations in population dynamics.     

Sequence variation of HIV and bioinformatics

The envelope glycoprotein of human immunodeficiency virus type 1 (HIV-1) interacts with receptors on the target cell and mediates virus entry by fusing the viral and cell membranes. To maintain the viral infectivity, amino acids that interact with receptors are expected to be more conserved than the other sites on the protein surface. In contrast to the functional constraint of amino acids for the receptor binding, some amino acid changes in this protein may produce antigenic variations that enable the virus to escape from recognition of the host immune system. Therefore, both positive selection (higher fitness) and negative selection (lower fitness) against amino acid changes are taking place during evolution of surface proteins of parasites To elucidate the evolutionary mechanisms of the whole HIV-1 gp120 envelope glycoprotein at the single site level, we collected and analyzed all available sequence data for the protein. By analyzing 186 sequences of the HIV-1 gp120 (subtype B), we reevaluated amino acid variability at the single site level, and estimated the numbers of synonymous and nonsynonymous substitutions at each codon position to detect positive and negative selection. We identified 33 amino acid positions which may be under positive selection. Some of these positions may form discontinuous epitopes. We also analyzed amino acid sequences to find amino acid positions responsible for usage of the second receptor. We found that, in addition to the V3 loop, amino acid variation at residue 440 in C4 region is clearly linked with the usage of CXCR 4.     

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.     

Subclonal components of consensus fitness in an RNA virus clone.

Most RNA virus populations exhibit extremely high mutation frequencies which generate complex, genetically heterogeneous populations referred to as quasi-species. Previous work has shown that when a large spectrum of the quasi-species is transferred, natural selection operates, leading to elimination of noncompetitive (inferior) genomes and rapid gains in fitness. However, whenever the population is repeatedly reduced to a single virion, variable declines in fitness occur as predicted by the Muller's ratchet hypothesis. Here, we quantitated the fitness of 98 subclones isolated from an RNA virus clonal population. We found a normal distribution around a lower fitness, with the average subclone being less fit than the parental clonal population. This finding demonstrates the phenotypic diversity in RNA virus populations and shows that, as expected, a large fraction of mutations generated during virus replication is deleterious. This clarifies the operation of Muller's ratchet and illustrates why a large number of virions must be transferred for rapid fitness gains to occur. We also found that repeated genetic bottleneck passages can cause irregular stochastic declines in fitness, emphasizing again the phenotypic heterogeneity present in RNA virus populations. Finally, we found that following only 60 h of selection (15 passages in which virus yields were harvested after 4 h), RNA virus populations can undergo a 250% average increase in fitness, even on a host cell type to which they were already well adapted. This is a remarkable ability; in population biology, even a much lower fitness gain (e.g., 1 to 2%) can represent a highly significant reproductive advantage. We discuss the biological implications of these findings for the natural transmission and pathogenesis of RNA viruses.     

A genetic background with low mutational robustness is associated with increased adaptability to a novel host in an RNA virus

Although mutational robustness is central to many evolutionary processes, its relationship to evolvability remains poorly understood and has been very rarely tested experimentally. Here, we measure the evolvability of Vesicular stomatitis virus in two genetic backgrounds with different levels of mutational robustness. We passaged the viruses into a novel cell type to model a host‐jump episode, quantified changes in infectivity and fitness in the new host, evaluated the cost of adaptation in the original host and analyzed the genetic basis of this adaptation. Lineages evolved from the less robust genetic background demonstrated increased adaptability, paid similar costs of adaptation to the new host and fixed approximately the same number of mutations as their more robust counterparts. Theory predicts that robustness can promote evolvability only in systems where large sets of genotypes are connected by effectively neutral mutations. We argue that this condition might not be fulfilled generally in RNA viruses.     

Generalized selection to overcome innate immunity selects for host breadth in an RNA virus

Virus‐host coevolution has selected for generalized host defense against viruses, exemplified by interferon production/signaling and other innate immune function in eukaryotes such as humans. Although cell‐surface binding primarily limits virus infection success, generalized adaptation to counteract innate immunity across disparate hosts may contribute to RNA virus emergence potential. We examined this idea using vesicular stomatitis virus (VSV) populations previously evolved on strictly immune‐deficient (HeLa) cells, strictly immune competent (MDCK) cells, or on alternating deficient/competent cells. By measuring viral fitness in unselected human cancer cells of differing innate immunity, we confirmed that HeLa‐adapted populations were specialized for innate immune‐deficient hosts, whereas MDCK‐adapted populations were relatively more generalized for fitness on hosts of differing innate immune capacity and of different species origin. We also confirmed that HeLa‐evolved populations maintained fitness in immune‐deficient nonhuman primate cells. These results suggest that innate immunity is more prominent than host species in determining viral fitness at the host‐cell level. Finally, our prediction was inexact that selection on alternating deficient/competent hosts should produce innate viral generalists. Rather, fitness differences among alternating host‐evolved VSV populations indicated variable capacities to evade innate immunity. Our results suggest that the evolutionary history of innate immune selection can affect whether RNA viruses evolve greater host‐breadth.     

Vector-Borne Transmission Imposes a Severe Bottleneck on an RNA Virus Population

RNA viruses typically occur in genetically diverse populations due to their error-prone genome replication. Genetic diversity is thought to be important in allowing RNA viruses to explore sequence space, facilitating adaptation to changing environments and hosts. Some arboviruses that infect both a mosquito vector and a mammalian host are known to experience population bottlenecks in their vectors, which may constrain their genetic diversity and could potentially lead to extinction events via Muller''s ratchet. To examine this potential challenge of bottlenecks for arbovirus perpetuation, we studied Venezuelan equine encephalitis virus (VEEV) enzootic subtype IE and its natural vector Culex (Melanoconion) taeniopus, as an example of a virus-vector interaction with a long evolutionary history. Using a mixture of marked VEEV clones to infect C. taeniopus and real-time RT-PCR to track these clones during mosquito infection and dissemination, we observed severe bottleneck events that resulted in a significant drop in the number of clones present. At higher initial doses, the midgut was readily infected and there was a severe bottleneck at the midgut escape. Following a lower initial dose, the major bottleneck occurred at initial midgut infection. A second, less severe bottleneck was identified at the salivary gland infection stage following intrathoracic inoculation. Our results suggest that VEEV consistently encounters bottlenecks during infection, dissemination and transmission by its natural enzootic vector. The potential impacts of these bottlenecks on viral fitness and transmission, and the viral mechanisms that prevent genetic drift leading to extinction, deserve further study.