Exponential Fitness Gains of RNA Virus Populations Are Limited by Bottleneck Effects

Fitness is a parameter that quantitatively measures adaptation of a virus to a given environment. We have previously reported exponential fitness gains of large populations of vesicular stomatitis virus replicating in a constant environment (I. S. Novella et al., Proc. Natl. Acad. Sci. USA 92:5841–5844, 1995). In this paper, we report that during long-term passage of such large viral populations, fitness values reached a high-fitness plateau during which stochastic fitness variations were observed. This effect appears likely to be due to bottleneck effects on very high fitness populations.     

Escape from growth restriction in small colony variants of Salmonella typhimurium by gene amplification and mutation

Antibiotic resistance in bacteria is generally associated with fitness costs that often can be reduced by second-site compensatory mutations. Here, we examined how a protamine-resistant small colony variant of Salmonella typhimurium adapts to the growth reduction conferred by a resistance mutation in hemC (encoding a haem-biosynthesis enzyme). We show that adaptation occurs in a multi-step process where fitness is successively increased. Thus, the initial adaptive response was selection for an unstable gene amplification of the mutant hemC gene that provided a small fitness increase. Fitness was increased further by a mutation that restored HemC function in one gene copy, relaxing selection for the amplification. Subsequently, the amplification segregated back to the haploid state and even higher fitness. The end result was in most cases mutant strains with a hemC sequence different from that of the wild-type strain. These findings suggest that gene amplification facilitates adaptive evolution. A higher gene dosage increases the target size for compensatory mutations and improves fitness of the cell, thereby allowing an increase in the population size, further increasing the probability of a subsequent stable mutation. Our results provide a novel genetic basis for growth compensation in small colony variants.     

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.     

Sex‐specific selection under environmental stress in seed beetles

Sexual selection can increase rates of adaptation by imposing strong selection in males, thereby allowing efficient purging of the mutation load on population fitness at a low demographic cost. Indeed, sexual selection tends to be male‐biased throughout the animal kingdom, but little empirical work has explored the ecological sensitivity of this sex difference. In this study, we generated theoretical predictions of sex‐specific strengths of selection, environmental sensitivities and genotype‐by‐environment interactions and tested them in seed beetles by manipulating either larval host plant or rearing temperature. Using fourteen isofemale lines, we measured sex‐specific reductions in fitness components, genotype‐by‐environment interactions and the strength of selection (variance in fitness) in the juvenile and adult stage. As predicted, variance in fitness increased with stress, was consistently greater in males than females for adult reproductive success (implying strong sexual selection), but was similar in the sexes in terms of juvenile survival across all levels of stress. Although genetic variance in fitness increased in magnitude under severe stress, heritability decreased and particularly so in males. Moreover, genotype‐by‐environment interactions for fitness were common but specific to the type of stress, sex and life stage, suggesting that new environments may change the relative alignment and strength of selection in males and females. Our study thus exemplifies how environmental stress can influence the relative forces of natural and sexual selection, as well as concomitant changes in genetic variance in fitness, which are predicted to have consequences for rates of adaptation in sexual populations.     

A Multi-scale Analysis of Influenza A Virus Fitness Trade-offs due to Temperature-dependent Virus Persistence

Successful replication within an infected host and successful transmission between hosts are key to the continued spread of most pathogens. Competing selection pressures exerted at these different scales can lead to evolutionary trade-offs between the determinants of fitness within and between hosts. Here, we examine such a trade-off in the context of influenza A viruses and the differential pressures exerted by temperature-dependent virus persistence. For a panel of avian influenza A virus strains, we find evidence for a trade-off between the persistence at high versus low temperatures. Combining a within-host model of influenza infection dynamics with a between-host transmission model, we study how such a trade-off affects virus fitness on the host population level. We show that conclusions regarding overall fitness are affected by the type of link assumed between the within- and between-host levels and the main route of transmission (direct or environmental). The relative importance of virulence and immune response mediated virus clearance are also found to influence the fitness impacts of virus persistence at low versus high temperatures. Based on our results, we predict that if transmission occurs mainly directly and scales linearly with virus load, and virulence or immune responses are negligible, the evolutionary pressure for influenza viruses to evolve toward good persistence at high within-host temperatures dominates. For all other scenarios, influenza viruses with good environmental persistence at low temperatures seem to be favored.     

Immuno-epidemiological Modeling of HIV-1 Predicts High Heritability of the Set-Point Virus Load,while Selection for CTL Escape Dominates Virulence Evolution

It has been suggested that HIV-1 has evolved its set-point virus load to be optimized for transmission. Previous epidemiological models and studies into the heritability of set-point virus load confirm that this mode of adaptation within the human population is feasible. However, during the many cycles of replication between infection of a host and transmission to the next host, HIV-1 is under selection for escape from immune responses, and not transmission. Here we investigate with computational and mathematical models how these two levels of selection, within-host and between-host, are intertwined. We find that when the rate of immune escape is comparable to what has been observed in patients, immune selection within hosts is dominant over selection for transmission. Surprisingly, we do find high values for set-point virus load heritability, and argue that high heritability estimates can be caused by the ‘footprints’ left by differing hosts'' immune systems on the virus.     

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.     

Recombination Enhances HIV-1 Envelope Diversity by Facilitating the Survival of Latent Genomic Fragments in the Plasma Virus Population

HIV-1 is subject to immune pressure exerted by the host, giving variants that escape the immune response an advantage. Virus released from activated latent cells competes against variants that have continually evolved and adapted to host immune pressure. Nevertheless, there is increasing evidence that virus displaying a signal of latency survives in patient plasma despite having reduced fitness due to long-term immune memory. We investigated the survival of virus with latent envelope genomic fragments by simulating within-host HIV-1 sequence evolution and the cycling of viral lineages in and out of the latent reservoir. Our model incorporates a detailed mutation process including nucleotide substitution, recombination, latent reservoir dynamics, diversifying selection pressure driven by the immune response, and purifying selection pressure asserted by deleterious mutations. We evaluated the ability of our model to capture sequence evolution in vivo by comparing our simulated sequences to HIV-1 envelope sequence data from 16 HIV-infected untreated patients. Empirical sequence divergence and diversity measures were qualitatively and quantitatively similar to those of our simulated HIV-1 populations, suggesting that our model invokes realistic trends of HIV-1 genetic evolution. Moreover, reconstructed phylogenies of simulated and patient HIV-1 populations showed similar topological structures. Our simulation results suggest that recombination is a key mechanism facilitating the persistence of virus with latent envelope genomic fragments in the productively infected cell population. Recombination increased the survival probability of latent virus forms approximately 13-fold. Prevalence of virus with latent fragments in productively infected cells was observed in only 2% of simulations when we ignored recombination, while the proportion increased to 27% of simulations when we allowed recombination. We also found that the selection pressures exerted by different fitness landscapes influenced the shape of phylogenies, diversity trends, and survival of virus with latent genomic fragments. Our model predicts that the persistence of latent genomic fragments from multiple different ancestral origins increases sequence diversity in plasma for reasonable fitness landscapes.     

The roles of natural and sexual selection during adaptation to a novel environment

The net effect of sexual selection on nonsexual fitness is controversial. On one side, elaborate display traits and preferences for them can be costly, reducing the nonsexual fitness of individuals possessing them, as well as their offspring. In contrast, sexual selection may reinforce nonsexual fitness if an individual's attractiveness and quality are genetically correlated. According to recent models, such good-genes mate choice should increase both the extent and rate of adaptation. We evolved 12 replicate populations of Drosophila serrata in a powerful two-way factorial experimental design to test the separate and combined contributions of natural and sexual selection to adaptation to a novel larval food resource. Populations evolving in the presence of natural selection had significantly higher mean nonsexual fitness when measured over three generations (13-15) during the course of experimental evolution (16-23% increase). The effect of natural selection was even more substantial when measured in a standardized, monogamous mating environment at the end of the experiment (generation 16; 52% increase). In contrast, and despite strong sexual selection on display traits, there was no evidence from any of the four replicate fitness measures that sexual selection promoted adaptation. In addition, a comparison of fitness measures conducted under different mating environments demonstrated a significant direct cost of sexual selection to females, likely arising from some form of male-induced harm. Indirect benefits of sexual selection in promoting adaptation to this novel resource environment therefore appear to be absent in this species, despite prior evidence suggesting the operation of good-genes mate choice in their ancestral environment. How novel environments affect the operation of good-genes mate choice is a fundamental question for future sexual selection research.     

Adaptation of Puumala hantavirus to cell culture is associated with point mutations in the coding region of the L segment and in the noncoding regions of the S segment

We previously developed a model for studies on hantavirus host adaptation and initiated genetic analysis of Puumala virus variants passaged in colonized bank voles and in cultured Vero E6 cells. With the data presented in this paper, the sequence comparison of the wild-type and Vero E6-adapted variants of Puumala virus, strain Kazan, has been completed. The only amino acid substitution that distinguished the two virus variants was found in the L protein, Ser versus Phe at position 2053. Another mutation found in the L segment, the silent transition C1053U, could result from the selection of a variant with altered L RNA folding. Nucleotide substitutions observed in individual L cDNA clones, most of them A-->G and U-->C transitions, suggested that the population of L RNA molecules is represented by quasispecies. The mutation frequency in the L segment quasispecies appeared to be similar to the corresponding values for the S and M quasispecies. Analysis of the cDNA clones with the complete S segment sequences from passage 20 confirmed our earlier conclusion that the cell-adapted genotype of the virus is represented mostly by variants with mutated S segment noncoding regions. However, the spectrum of the S segment quasispecies appeared to be changing, suggesting that, after the initial adaptation (passages 1 to 11), the viral population is still being driven by selection for variants with higher fitness.     

Hepatitis C Virus Envelope Glycoprotein Fitness Defines Virus Population Composition following Transmission to a New Host

Genetic variability is a hallmark of RNA virus populations. However, transmission to a new host often results in a marked decrease in population diversity. This genetic bottlenecking is observed during hepatitis C virus (HCV) transmission and can arise via a selective sweep or through the founder effect. To model HCV transmission, we utilized chimeric SCID/Alb-uPA mice with transplanted human hepatocytes and infected them with a human serum HCV inoculum. E1E2 glycoprotein gene sequences in the donor inoculum and recipient mice were determined following single-genome amplification (SGA). In independent experiments, using mice with liver cells grafted from different sources, an E1E2 variant undetectable in the source inoculum was selected for during transmission. Bayesian coalescent analyses indicated that this variant arose in the inoculum pretransmission. Transmitted variants that established initial infection harbored key substitutions in E1E2 outside HVR1. Notably, all posttransmission E1E2s had lost a potential N-linked glycosylation site (PNGS) in E2. In lentiviral pseudoparticle assays, the major posttransmission E1E2 variant conferred an increased capacity for entry compared to the major variant present in the inoculum. Together, these data demonstrate that increased envelope glycoprotein fitness can drive selective outgrowth of minor variants posttransmission and that loss of a PNGS is integral to this improved phenotype. Mathematical modeling of the dynamics of competing HCV variants indicated that relatively modest differences in glycoprotein fitness can result in marked shifts in virus population composition. Overall, these data provide important insights into the dynamics and selection of HCV populations during transmission.     

Hypervariable Domain of Nonstructural Protein nsP3 of Venezuelan Equine Encephalitis Virus Determines Cell-Specific Mode of Virus Replication

Venezuelan equine encephalitis virus (VEEV) is one of the most pathogenic members of the Alphavirus genus in the Togaviridae family. This genus is divided into the Old World and New World alphaviruses, which demonstrate profound differences in pathogenesis, replication, and virus-host interactions. VEEV is a representative member of the New World alphaviruses. The biology of this virus is still insufficiently understood, particularly the function of its nonstructural proteins in RNA replication and modification of the intracellular environment. One of these nonstructural proteins, nsP3, contains a hypervariable domain (HVD), which demonstrates very low overall similarity between different alphaviruses, suggesting the possibility of its function in virus adaptation to different hosts and vectors. The results of our study demonstrate the following. (i) Phosphorylation of the VEEV nsP3-specific HVD does not play a critical role in virus replication in cells of vertebrate origin but is important for virus replication in mosquito cells. (ii) The VEEV HVD is not required for viral RNA replication in the highly permissive BHK-21 cell line. In fact, it can be either completely deleted or replaced by a heterologous protein sequence. These variants require only one or two additional adaptive mutations in nsP3 and/or nsP2 proteins to achieve an efficiently replicating phenotype. (iii) However, the carboxy-terminal repeat in the VEEV HVD is indispensable for VEEV replication in the cell lines other than BHK-21 and plays a critical role in formation of VEEV-specific cytoplasmic protein complexes. Natural VEEV variants retain at least one of the repeated elements in their nsP3 HVDs.     

Drosophila Adaptation to Viral Infection through Defensive Symbiont Evolution

Microbial symbionts can modulate host interactions with biotic and abiotic factors. Such interactions may affect the evolutionary trajectories of both host and symbiont. Wolbachia protects Drosophila melanogaster against several viral infections and the strength of the protection varies between variants of this endosymbiont. Since Wolbachia is maternally transmitted, its fitness depends on the fitness of its host. Therefore, Wolbachia populations may be under selection when Drosophila is subjected to viral infection. Here we show that in D. melanogaster populations selected for increased survival upon infection with Drosophila C virus there is a strong selection coefficient for specific Wolbachia variants, leading to their fixation. Flies carrying these selected Wolbachia variants have higher survival and fertility upon viral infection when compared to flies with the other variants. These findings demonstrate how the interaction of a host with pathogens shapes the genetic composition of symbiont populations. Furthermore, host adaptation can result from the evolution of its symbionts, with host and symbiont functioning as a single evolutionary unit.     

HIV coreceptor usage and drug treatment

The human immunodeficiency virus (HIV) infects a wide range of human cells. Cell entry is mediated through the CD4 receptor and a variety of coreceptors, most importantly the chemokine receptors CCR5 and CXCR4. Some antiretroviral agents selectively inhibit different HIV phenotypes depending on their coreceptor usage. Here, we analyse mathematical models, which describe the in vivo interaction of HIV phenotypes, differing in their coreceptor usage, with two target cell types (naive and memory CD4+ T cells). In particular, we investigate how the selection pressures on CCR5- and CXCR4-using HIV variants change as a result of treatment with coreceptor-specific antiretroviral agents. Our main result is that CXCR4 inhibitors increase the selection pressure in favor of the emergence of CCR5-using variants, thus selecting for coexistence of CXCR4- and CCR5-using variants, whereas CCR5 inhibitors increase the selection pressure against CCR5-using variants, thus selecting against coexistence. These results shed new light on the potential risks and benefits of coreceptor inhibitors.     

Narrow Bottlenecks Affect Pea Seedborne Mosaic Virus Populations during Vertical Seed Transmission but not during Leaf Colonization

The effective size of populations (Ne) determines whether selection or genetic drift is the predominant force shaping their genetic structure and evolution. Populations having high Ne adapt faster, as selection acts more intensely, than populations having low Ne, where random effects of genetic drift dominate. Estimating Ne for various steps of plant virus life cycle has been the focus of several studies in the last decade, but no estimates are available for the vertical transmission of plant viruses, although virus seed transmission is economically significant in at least 18% of plant viruses in at least one plant species. Here we study the co-dynamics of two variants of Pea seedborne mosaic virus (PSbMV) colonizing leaves of pea plants (Pisum sativum L.) during the whole flowering period, and their subsequent transmission to plant progeny through seeds. Whereas classical estimators of Ne could be used for leaf infection at the systemic level, as virus variants were equally competitive, dedicated stochastic models were needed to estimate Ne during vertical transmission. Very little genetic drift was observed during the infection of apical leaves, with Ne values ranging from 59 to 216. In contrast, a very drastic genetic drift was observed during vertical transmission, with an average number of infectious virus particles contributing to the infection of a seedling from an infected mother plant close to one. A simple model of vertical transmission, assuming a cumulative action of virus infectious particles and a virus density threshold required for vertical transmission to occur fitted the experimental data very satisfactorily. This study reveals that vertically-transmitted viruses endure bottlenecks as narrow as those imposed by horizontal transmission. These bottlenecks are likely to slow down virus adaptation and could decrease virus fitness and virulence.     

Experimental Evolution of an RNA Virus in Wild Birds: Evidence for Host-Dependent Impacts on Population Structure and Competitive Fitness

Within hosts, RNA viruses form populations that are genetically and phenotypically complex. Heterogeneity in RNA virus genomes arises due to error-prone replication and is reduced by stochastic and selective mechanisms that are incompletely understood. Defining how natural selection shapes RNA virus populations is critical because it can inform treatment paradigms and enhance control efforts. We allowed West Nile virus (WNV) to replicate in wild-caught American crows, house sparrows and American robins to assess how natural selection shapes RNA virus populations in ecologically relevant hosts that differ in susceptibility to virus-induced mortality. After five sequential passages in each bird species, we examined the phenotype and population diversity of WNV through fitness competition assays and next generation sequencing. We demonstrate that fitness gains occur in a species-specific manner, with the greatest replicative fitness gains in robin-passaged WNV and the least in WNV passaged in crows. Sequencing data revealed that intrahost WNV populations were strongly influenced by purifying selection and the overall complexity of the viral populations was similar among passaged hosts. However, the selective pressures that control WNV populations seem to be bird species-dependent. Specifically, crow-passaged WNV populations contained the most unique mutations (~1.7× more than sparrows, ~3.4× more than robins) and defective genomes (~1.4× greater than sparrows, ~2.7× greater than robins), but the lowest average mutation frequency (about equal to sparrows, ~2.6× lower than robins). Therefore, our data suggest that WNV replication in the most disease-susceptible bird species is positively associated with virus mutational tolerance, likely via complementation, and negatively associated with the strength of selection. These differences in genetic composition most likely have distinct phenotypic consequences for the virus populations. Taken together, these results reveal important insights into how different hosts may contribute to the emergence of RNA viruses.     

Functional Interactions of the Influenza Virus Glycoproteins

Hemagglutinin (HA) and neuraminidase (NA) are functionally related envelope glycoproteins of the influenza virus (Flu). HA interacts with terminal sialyl residues of oligosaccharides and ensures the binding of the virus particle to the cell surface. NA is a receptor-destroying enzyme that removes sialyl residues from oligosaccharides contained in cell and virus components and thereby prevents aggregation of virus particles. Analysis of reassortants combining low-functional NA of human Flu with HA of avian Flu showed that sialyl residues are not completely removed in some cases. With high HA affinity for sialyl substrates, such virus particles aggregate, aggregates accumulate on the cell surface, and virus yield decreases. Serial passaging of low-yield aggregating reassortants may lead to selection of high-yield variants, which do not aggregate. The loss of aggregation is due to a decrease in HA affinity for high-molecular-weight sialyl substrates. On evidence of sequencing of the HA gene in original reassortants and their nonaggregating variants, for different HA antigenic subtypes (H2, H3, H4, and H13) the affinity is reduced and aggregation lost through a common mechanism: an increase in the negative charge as a result of an amino acid substitution in the vicinity of the receptor-binding pocket of HA. Taken together, these findings suggest a way of postreassortment adaptation, which improves the functional match of HA and NA. The experimental system employed provides a model of natural processes associated with emergence of Flu variants having a pandemic potential.