Competition-colonization dynamics in an RNA virus

Coinfection and Interference Phenomena Are the Results of Multiple Thermodynamic Competitive Interactions

Biological, physical and chemical interaction between one (or more) microorganisms and a host organism, causing host cell damage, represents an infection. Infection of a plant, animal or microorganism with a virus can prevent infection with another virus. This phenomenon is known as viral interference. Viral interference is shown to result from two types of interactions, one taking place at the cell surface and the other intracellularly. Various viruses use different receptors to enter the same host cell, but various strains of one virus use the same receptor. The rate of virus–receptor binding can vary between different viruses attacking the same host, allowing interference or coinfection. The outcome of the virus–virus–host competition is determined by the Gibbs energies of binding and growth of the competing viruses and host. The virus with a more negative Gibbs energy of binding to the host cell receptor will enter the host first, while the virus characterized by a more negative Gibbs energy of growth will overtake the host metabolic machine and dominate. Once in the host cell, the multiplication machinery is shared by the competing viruses. Their potential to utilize it depends on the Gibbs energy of growth. Thus, the virus with a more negative Gibbs energy of growth will dominate. Therefore, the outcome can be interference or coinfection, depending on both the attachment kinetics (susceptibility) and the intracellular multiplication machinery (permittivity). The ratios of the Gibbs energies of binding and growth of the competing viruses determine the outcome of the competition. Based on this, a predictive model of virus–virus competition is proposed.

Species exclusion and coexistence in a noisy voter model with a competition-colonization tradeoff

We introduce an asymmetric noisy voter model to study the joint effect of immigration and a competition-dispersal tradeoff in the dynamics of two species competing for space in regular lattices. Individuals of one species can invade a nearest-neighbor site in the lattice, while individuals of the other species are able to invade sites at any distance but are less competitive locally, i.e., they establish with a probability g≤1. The model also accounts for immigration, modeled as an external noise that may spontaneously replace an individual at a lattice site by another individual of the other species. This combination of mechanisms gives rise to a rich variety of outcomes for species competition, including exclusion of either species, monostable coexistence of both species at different population proportions, and bistable coexistence with proportions of populations that depend on the initial condition. Remarkably, in the bistable phase, the system undergoes a discontinuous transition as the intensity of immigration overcomes a threshold, leading to a half loop dynamics associated to a cusp catastrophe, which causes the irreversible loss of the species with the shortest dispersal range.

Increasing temperature alters the within-host competition of viral strains and influences virus genetic variability

Environmental conditions can affect viral accumulation, virulence and adaptation, which have implications in the disease outcomes and efficiency of control measures. Concurrently, mixed viral infections are relevant in plants, being their epidemiology shaped by within-host virus–virus interactions. However, the extent in which the combined effect of variations in abiotic components of the plant ecological niche and the prevalence of mixed infections affect the evolutionary dynamics of viral populations is not well understood. Here, we explore the interplay between ecological and evolutionary factors during viral infections and show that isolates of two strains of Pepino mosaic potexvirus coexisted in tomato plants in a temperature-dependent continuum between neutral and antagonistic interactions. After a long-term infection, the mutational analysis of the evolved viral genomes revealed strain-specific single-nucleotide polymorphisms that were modulated by the interaction between the type of infection and temperature. These results suggest that the temperature is an ecological driver of virus-virus interactions, with an effect on the genetic diversity of individual viruses that are co-infecting an individual host. This research provides insights into the effect that changes in host growth temperatures might have on the evolutionary dynamics of viral populations in mixed infections.

Implications of mixed viral infections on plant disease ecology and evolution

Mixed viral infections occur more commonly than would be expected by chance in nature. Virus-virus interactions may affect viral traits and leave a genetic signature in the population, and thus influence the prevalence and emergence of viral diseases. Understanding about how the interactions between viruses within a host shape the evolutionary dynamics of the viral populations is needed for viral disease prevention and management. Here, we first synthesize concepts implied in the occurrence of virus-virus interactions. Second, we consider the role of the within-host interactions of virus-virus and virus-other pathogenic microbes, on the composition and structure of viral populations. Third, we contemplate whether mixed viral infections can create opportunities for the generation and maintenance of viral genetic diversity. Fourth, we attempt to summarize the evolutionary response of viral populations to mixed infections to understand how they shape the spatio-temporal dynamics of viral populations at the individual plant and field scales. Finally, we anticipate the future research under the reconciliation of molecular epidemiology and evolutionary ecology, drawing attention to the need of adding more complexity to future research in order to gain a better understanding about the mechanisms operating in nature.

A new implication of quasispecies dynamics: Broad virus diversification in absence of external perturbations

RNA genetic elements include many important animal and plant pathogens. They share high mutability, a trait that has multiple implications for the interactions with their host organisms. Here we review evidence of a new adaptive feature of RNA viruses that we term "broadly diversifying selection". It constitutes a new type of positive selection without participation of any external selective agent, and which is built upon a progressive increase of the number of different genomes that dominate the population. The evidence was provided by analyses of mutant spectrum composition of two important viral pathogens, foot-and-mouth disease virus (FMDV) and hepatitis C virus (HCV) after prolonged replication in their respective cell culture environment. Despite being fueled by mutations that arise randomly and in absence of an external guiding selective force, this type of selection prepares the viral population for a response to selective forces still to occur. Since current evidence suggests that broadly diversifying selection is favored by elevated mutation rates and population sizes, it may constitute a more general behavior, relevant also to the adaptive dynamics of microbial populations and cancer cells.

Long-Term Cocirculation of Two Strains of Pepino Mosaic Virus in Tomato Crops and Its Effect on Population Genetic Variability

Mixed viral infections are common in plants, and the evolutionary dynamics of viral populations may differ depending on whether the infection is caused by single or multiple viral strains. However, comparative studies of single and mixed infections using viral populations in comparable agricultural and geographical locations are lacking. Here, we monitored the occurrence of pepino mosaic virus (PepMV) in tomato crops in two major tomato-producing areas in Murcia (southeastern Spain), supporting evidence showing that PepMV disease-affected plants had single infections of the Chilean 2 (CH2) strain in one area and the other area exhibited long-term (13 years) coexistence of the CH2 and European (EU) strains. We hypothesized that circulating strains of PepMV might be modulating the differentiation between them and shaping the evolutionary dynamics of PepMV populations. Our phylogenetic analysis of 106 CH2 isolates randomly selected from both areas showed a remarkable divergence between the CH2 isolates, with increased nucleotide variability in the geographical area where both strains cocirculate. Furthermore, the potential virus-virus interaction was studied further by constructing six full-length infectious CH2 clones from both areas, and assessing their viral fitness in the presence and absence of an EU-type isolate. All CH2 clones showed decreased fitness in mixed infections and although complete genome sequencing indicated a nucleotide divergence of those CH2 clones by area, the magnitude of the fitness response was irrespective of the CH2 origin. Overall, these results suggest that although agroecological cropping practices may be particularly important for explaining the evolutionary dynamics of PepMV in tomato crops, the cocirculation of both strains may have implications on the genetic variability of PepMV populations.

Virus population dynamics examined with experimental model systems

Experimental evolution permits exploring the effect of controlled environmental variables in virus evolution. Several designs in cell culture and in vivo have established basic concepts that can assist in the interpretation of evolutionary events in the field. Important information has come from cytolytic and persistent infections in cell culture that have unveiled the power of virus-cell coevolution in virus and cell diversification. Equally informative are comparisons of the response of viral populations when subjected to different passage régimens. In particular, plaque-to-plaque transfers in cell culture have revealed unusual genotypes and phenotypes that populate minority layers of viral quasispecies. Some of these viruses display properties that contradict features established in virology textbooks. Several hypotheses and principles of population genetics have found experimental confirmation in experimental designs with viruses. The possibilities of using experimental evolution to understand virus behavior are still largely unexploited.

Viral quasispecies

Viral quasispecies refers to a population structure that consists of extremely large numbers of variant genomes, termed mutant spectra, mutant swarms or mutant clouds. Fueled by high mutation rates, mutants arise continually, and they change in relative frequency as viral replication proceeds. The term quasispecies was adopted from a theory of the origin of life in which primitive replicons) consisted of mutant distributions, as found experimentally with present day RNA viruses. The theory provided a new definition of wild type, and a conceptual framework for the interpretation of the adaptive potential of RNA viruses that contrasted with classical studies based on consensus sequences. Standard clonal analyses and deep sequencing methodologies have confirmed the presence of myriads of mutant genomes in viral populations, and their participation in adaptive processes. The quasispecies concept applies to any biological entity, but its impact is more evident when the genome size is limited and the mutation rate is high. This is the case of the RNA viruses, ubiquitous in our biosphere, and that comprise many important pathogens. In virology, quasispecies are defined as complex distributions of closely related variant genomes subjected to genetic variation, competition and selection, and that may act as a unit of selection. Despite being an integral part of their replication, high mutation rates have an upper limit compatible with inheritable information. Crossing such a limit leads to RNA virus extinction, a transition that is the basis of an antiviral design termed lethal mutagenesis.

Marine RNA Virus Quasispecies Are Distributed throughout the Oceans

Very little is known about aquatic RNA virus populations and genome evolution. This is the first study that analyzes marine environmental RNA viral assemblages in an evolutionary and broad geographical context. This study contributes the largest marine RNA virus metagenomic data set to date, substantially increasing the sequencing space for RNA viruses and also providing a baseline for comparisons of marine RNA virus diversity. The new viruses discovered in this study are representative of the most abundant family of marine RNA viruses, the Marnaviridae, and expand our view of the diversity of this important group. Overall, our data and analyses provide a foundation for interpreting marine RNA virus diversity and evolution.

RNA viruses, particularly genetically diverse members of the Picornavirales, are widespread and abundant in the ocean. Gene surveys suggest that there are spatial and temporal patterns in the composition of RNA virus assemblages, but data on their diversity and genetic variability in different oceanographic settings are limited. Here, we show that specific RNA virus genomes have widespread geographic distributions and that the dominant genotypes are under purifying selection. Genomes from three previously unknown picorna-like viruses (BC-1, -2, and -3) assembled from a coastal site in British Columbia, Canada, as well as marine RNA viruses JP-A, JP-B, and Heterosigma akashiwo RNA virus exhibited different biogeographical patterns. Thus, biotic factors such as host specificity and viral life cycle, and not just abiotic processes such as dispersal, affect marine RNA virus distribution. Sequence differences relative to reference genomes imply that virus quasispecies are under purifying selection, with synonymous single-nucleotide variations dominating in genomes from geographically distinct regions resulting in conservation of amino acid sequences. Conversely, sequences from coastal South Africa that mapped to marine RNA virus JP-A exhibited more nonsynonymous mutations, probably representing amino acid changes that accumulated over a longer separation. This biogeographical analysis of marine RNA viruses demonstrates that purifying selection is occurring across oceanographic provinces. These data add to the spectrum of known marine RNA virus genomes, show the importance of dispersal and purifying selection for these viruses, and indicate that closely related RNA viruses are pathogens of eukaryotic microbes across oceans.

IMPORTANCE Very little is known about aquatic RNA virus populations and genome evolution. This is the first study that analyzes marine environmental RNA viral assemblages in an evolutionary and broad geographical context. This study contributes the largest marine RNA virus metagenomic data set to date, substantially increasing the sequencing space for RNA viruses and also providing a baseline for comparisons of marine RNA virus diversity. The new viruses discovered in this study are representative of the most abundant family of marine RNA viruses, the Marnaviridae, and expand our view of the diversity of this important group. Overall, our data and analyses provide a foundation for interpreting marine RNA virus diversity and evolution.

Life cycle synchronization is a viral drug resistance mechanism

Viral infections are one of the major causes of death worldwide, with HIV infection alone resulting in over 1.2 million casualties per year. Antiviral drugs are now being administered for a variety of viral infections, including HIV, hepatitis B and C, and influenza. These therapies target a specific phase of the virus's life cycle, yet their ultimate success depends on a variety of factors, such as adherence to a prescribed regimen and the emergence of viral drug resistance. The epidemiology and evolution of drug resistance have been extensively characterized, and it is generally assumed that drug resistance arises from mutations that alter the virus's susceptibility to the direct action of the drug. In this paper, we consider the possibility that a virus population can evolve towards synchronizing its life cycle with the pattern of drug therapy. The periodicity of the drug treatment could then allow for a virus strain whose life cycle length is a multiple of the dosing interval to replicate only when the concentration of the drug is lowest. This process, referred to as "drug tolerance by synchronization", could allow the virus population to maximize its overall fitness without having to alter drug binding or complete its life cycle in the drug's presence. We use mathematical models and stochastic simulations to show that life cycle synchronization can indeed be a mechanism of viral drug tolerance. We show that this effect is more likely to occur when the variability in both viral life cycle and drug dose timing are low. More generally, we find that in the presence of periodic drug levels, time-averaged calculations of viral fitness do not accurately predict drug levels needed to eradicate infection, even if there is no synchronization. We derive an analytical expression for viral fitness that is sufficient to explain the drug-pattern-dependent survival of strains with any life cycle length. We discuss the implications of these findings for clinically relevant antiviral strategies.

A Selective Bottleneck Shapes the Evolutionary Mutant Spectra of Enterovirus A71 during Viral Dissemination in Humans

RNA viruses accumulate mutations to rapidly adapt to environmental changes. Enterovirus A71 (EV-A71) causes various clinical manifestations with occasional severe neurological complications. However, the mechanism by which EV-A71 evolves within the human body is unclear. Utilizing deep sequencing and haplotype analyses of viruses from various tissues of an autopsy patient, we sought to define the evolutionary pathway by which enterovirus A71 evolves fitness for invading the central nervous system in humans. Broad mutant spectra with divergent mutations were observed at the initial infection sites in the respiratory and digestive systems. After viral invasion, we identified a haplotype switch and dominant haplotype, with glycine at VP1 residue 31 (VP1-31G) in viral particles disseminated into the integumentary and central nervous systems. In vitro viral growth and fitness analyses indicated that VP1-31G conferred growth and a fitness advantage in human neuronal cells, whereas VP1-31D conferred enhanced replication in human colorectal cells. A higher proportion of VP1-31G was also found among fatal cases, suggesting that it may facilitate central nervous system infection in humans. Our data provide the first glimpse of EV-A71 quasispecies from oral tissues to the central nervous system within humans, showing broad implications for the surveillance and pathogenesis of this reemerging viral pathogen.IMPORTANCE EV-A71 continues to be a worldwide burden to public health. Although EV-A71 is the major etiological agent of hand, foot, and mouth disease, it can also cause neurological pulmonary edema, encephalitis, and even death, especially in children. Understanding selection processes enabling dissemination and accurately estimating EV-A71 diversity during invasion in humans are critical for applications in viral pathogenesis and vaccine studies. Here, we define a selection bottleneck appearing in respiratory and digestive tissues. Glycine substitution at VP1 residue 31 helps viruses break through the bottleneck and invade the central nervous system. This substitution is also advantageous for replication in neuronal cells in vitro Considering that fatal cases contain enhanced glycine substitution at VP1-31, we suggest that the increased prevalence of VP1-31G may alter viral tropism and aid central nervous system invasion. Our findings provide new insights into a dynamic mutant spectral switch active during acute viral infection with emerging viral pathogens.

Viral Population Changes during Murine Norovirus Propagation in RAW 264.7 Cells

Laboratory adaptation of viruses is an essential technique for basic virology research, including the generation of attenuated vaccine strains, although the principles of cell adaptation remain largely unknown. Deep sequencing of murine norovirus (MuNoV) S7 during serial passages in RAW264.7 cells showed that the frequencies of viral variants were altered more dynamically than previously reported. Serial passages of the virus following two different multiplicity of infections gave rise to distinct haplotypes, implying that multiple cell-adaptable sequences were present in the founder population. Nucleotide variants lost during passage were assembled into a viral genome representative of that prior to cell adaptation, which was unable to generate viral particles upon infection in cultured cells. In addition, presence of the reconstructed genome interfered with production of infectious particles from viruses that were fully adapted to in vitro culture. Although the key nucleotide changes dictating cell adaptation of MuNoV S7 viral infection are yet to be elucidated, our results revealed the elaborate interplay among selected sequences of viral variants better adapted to propagation in cell culture. Such knowledge will be instrumental in understanding the processes necessary for the laboratory adaptation of viruses, especially to those without relevant cell culture systems.

Quasispecies characteristics in mother-to-child transmission of hepatitis B virus by next-generation sequencing

OBJECTIVES

To identify within-host quasispecies characteristics of hepatitis B virus (HBV) in mothers and children infected via mother-to-child transmission (MTCT).

METHODS

Using next-generation sequencing (NGS), we analyzed sequences within the non-overlapping pre-core/core (pre-C/C) gene in 37 mother-child pairs.

RESULTS

Phylogenetic and Highlighter analyses suggested that both a single strain and multiple distinct strains may be transmitted in MTCT of HBV. However, analysis of reassembled viral sequences revealed a relatively narrow distribution of variants in children, which was confirmed by a lower viral diversity in children than that in mothers. New closely related variants with combinations of two to five high-frequency mutations were observed in seven children with elevated ALT levels; the new variants out-competed the transmitted maternal variants to become the dominant strains in five of them. Furthermore, 30 mutations with a frequency >1% of all viruses within-host were present in those children; the mutations caused 19 amino-acid substitutions. Interestingly, almost all were located within the well-known T-cell or B-cell epitopes.

CONCLUSIONS

There are restrictive changes that occur in the early stages of chronic HBV infection through MTCT with different clinical consequences. These data might have important implications for future investigations of interrelated immunopathogenesis and therapeutic strategies.

Improving virus production through quasispecies genomic selection and molecular breeding

Virus production still is a challenging issue in antigen manufacture, particularly with slow-growing viruses. Deep-sequencing of genomic regions indicative of efficient replication may be used to identify high-fitness minority individuals suppressed by the ensemble of mutants in a virus quasispecies. Molecular breeding of quasispecies containing colonizer individuals, under regimes allowing more than one replicative cycle, is a strategy to select the fittest competitors among the colonizers. A slow-growing cell culture-adapted hepatitis A virus strain was employed as a model for this strategy. Using genomic selection in two regions predictive of efficient translation, the internal ribosome entry site and the VP1-coding region, high-fitness minority colonizer individuals were identified in a population adapted to conditions of artificially-induced cellular transcription shut-off. Molecular breeding of this population with a second one, also adapted to transcription shut-off and showing an overall colonizer phenotype, allowed the selection of a fast-growing population of great biotechnological potential.

Differential Persistence of Foot-and-Mouth Disease Virus in African Buffalo Is Related to Virus Virulence

Foot-and-mouth disease (FMD) virus (FMDV) circulates as multiple serotypes and strains in many regions of endemicity. In particular, the three Southern African Territories (SAT) serotypes are maintained effectively in their wildlife reservoir, the African buffalo, and individuals may harbor multiple SAT serotypes for extended periods in the pharyngeal region. However, the exact site and mechanism for persistence remain unclear. FMD in buffaloes offers a unique opportunity to study FMDV persistence, as transmission from carrier ruminants has convincingly been demonstrated for only this species. Following coinfection of naive African buffaloes with isolates of three SAT serotypes from field buffaloes, palatine tonsil swabs were the sample of choice for recovering infectious FMDV up to 400 days postinfection (dpi). Postmortem examination identified infectious virus for up to 185 dpi and viral genomes for up to 400 dpi in lymphoid tissues of the head and neck, focused mainly in germinal centers. Interestingly, viral persistence in vivo was not homogenous, and the SAT-1 isolate persisted longer than the SAT-2 and SAT-3 isolates. Coinfection and passage of these SAT isolates in goat and buffalo cell lines demonstrated a direct correlation between persistence and cell-killing capacity. These data suggest that FMDV persistence occurs in the germinal centers of lymphoid tissue but that the duration of persistence is related to virus replication and cell-killing capacity.

IMPORTANCE

Foot-and-mouth disease virus (FMDV) causes a highly contagious acute vesicular disease in domestic livestock and wildlife species. African buffaloes (Syncerus caffer) are the primary carrier hosts of FMDV in African savannah ecosystems, where the disease is endemic. We have shown that the virus persists for up to 400 days in buffaloes and that there is competition between viruses during mixed infections. There was similar competition in cell culture: viruses that killed cells quickly persisted more efficiently in passaged cell cultures. These results may provide a mechanism for the dominance of particular viruses in an ecosystem.

Virus Population Dynamics Examined with Experimental Model Systems

Experimental evolution permits exploring the effect of controlled environmental variables in virus evolution. Several designs in cell culture and in vivo have established basic concepts that can assist in the interpretation of evolutionary events in the field. Important information has come from cytolytic and persistent infections in cell culture that have unveiled the power of virus-cell coevolution in virus and cell diversification. Equally informative are comparisons of the response of viral populations when subjected to different passage régimes. In particular, plaque-to-plaque transfers in cell culture have revealed unusual genotypes and phenotypes that populate minority layers of viral quasispecies. Some of these viruses display properties that contradict features established in virology textbooks. Several hypotheses and principles of population genetics have found experimental confirmation in experimental designs with viruses. The possibilities of using experimental evolution to understand virus behavior are still largely unexploited.

Composition and Interactions of Hepatitis B Virus Quasispecies Defined the Virological Response During Telbivudine Therapy

Reverse transcriptase (RT) mutations contribute to hepatitis B virus resistance during antiviral therapy with nucleos(t)ide analogs. However, the composition of the RT quasispecies and their interactions during antiviral treatment have not yet been thoroughly defined. In this report, 10 patients from each of 3 different virological response groups, i.e., complete virological response, partial virological response and virological breakthrough, were selected from a multicenter trial of Telbivudine treatment. Variations in the drug resistance-related critical RT regions in 107 serial serum samples from the 30 patients were examined by ultra-deep sequencing. A total of 496,577 sequence reads were obtained, with an average sequencing coverage of 4,641X per sample. The phylogenies of the quasispecies revealed the independent origins of two critical quasispecies, i.e., the rtA181T and rtM204I mutants. Data analyses and theoretical modeling showed a cooperative-competitive interplay among the quasispecies. In particular, rtM204I mutants compete against other quasispecies, which eventually leads to virological breakthrough. However, in the absence of rtM204I mutants, synergistic growth of the drug-resistant rtA181T mutants with the wild-type quasispecies could drive the composition of the viral population into a state of partial virological response. Furthermore, we demonstrated that the frequency of drug-resistant mutations in the early phase of treatment is important for predicting the virological response to antiviral therapy.

Effects of host heterogeneity on pathogen diversity and evolution

Phenotypic variation is common in most pathogens, yet the mechanisms that maintain this diversity are still poorly understood. We asked whether continuous host variation in susceptibility helps maintain phenotypic variation, using experiments conducted with a baculovirus that infects gypsy moth (Lymantria dispar) larvae. We found that an empirically observed tradeoff between mean transmission rate and variation in transmission, which results from host heterogeneity, promotes long-term coexistence of two pathogen types in simulations of a population model. This tradeoff introduces an alternative strategy for the pathogen: a low-transmission, low-variability type can coexist with the high-transmission type favoured by classical non-heterogeneity models. In addition, this tradeoff can help explain the extensive phenotypic variation we observed in field-collected pathogen isolates, in traits affecting virus fitness including transmission and environmental persistence. Similar heterogeneity tradeoffs might be a general mechanism promoting phenotypic variation in any pathogen for which hosts vary continuously in susceptibility.

Seasonal variation of respiratory pathogen colonization in asymptomatic health care professionals: A single-center, cross-sectional, 2-season observational study

BACKGROUND

The purpose of this study was to determine the seasonal variance of potentially pathogenic bacterial and viral organisms in nasopharyngeal specimens obtained from asymptomatic health care professionals (HCPs) during the 2014 winter and summer months.

METHODS

Nasopharyngeal specimens from 100 HCPs were collected from Huntsville Hospital (Huntsville, AL) during the winter and from 100 HCPs during the summer. All subjects were tested for 22 viruses and 19 bacteria using Target Enriched Multiplex Polymerase Chain Reaction. Both seasonal cohorts were composed of students, nurses, physicians, and residents.

RESULTS

Of the 100 HCPs tested during the winter, 34 subjects were colonized with at least 1 bacterium, and 11 tested positive for at least 1 virus. Methicillin-resistant Staphylococcus aureus (MRSA), Moraxella catarrhalis, and coronavirus were the most frequently detected potentially infectious agents. Of the 100 HCPs tested during the summer, 37 tested positive for at least 1 bacterium, and 4 tested positive for a viral agent. The most prevalent bacteria were MRSA and Klebsiella pneumonia.

CONCLUSION

Nasopharyngeal carriage among asymptomatic HCPs was common, but the frequency and presence of potential pathogens varied with each season. Understanding the colonization and infection potential of upper respiratory organisms is important, particularly for viruses. Although asymptomatic HCPs certainly harbor a number of different potentially infectious agents, future studies are needed to determine whether colonized pathogens are transmitted or initiate infection in at-risk patient populations.

Influence of multiple infection and relatedness on virulence: disease dynamics in an experimental plant population and its castrating parasite

The level of parasite virulence, i.e., the decrease in host's fitness due to a pathogen, is expected to depend on several parameters, such as the type of the disease (e.g., castrating or host-killing) and the prevalence of multiple infections. Although these parameters have been extensively studied theoretically, few empirical data are available to validate theoretical predictions. Using the anther smut castrating disease on Silene latifolia caused by Microbotryum lychnidis-dioicae, we studied the dynamics of multiple infections and of different components of virulence (host death, non-recovery and percentage of castrated stems) during the entire lifespan of the host in an experimental population. We monitored the number of fungal genotypes within plants and their relatedness across five years, using microsatellite markers, as well as the rates of recovery and host death in the population. The mean relatedness among genotypes within plants remained at a high level throughout the entire host lifespan despite the dynamics of the disease, with recurrent new infections. Recovery was lower for plants with multiple infections compared to plants infected by a single genotype. As expected for castrating parasites, M. lychnidis-dioicae did not increase host mortality. Mortality varied across years but was generally lower for plants that had been diseased the preceding year. This is one of the few studies to have empirically verified theoretical expectations for castrating parasites, and to show particularly i) that castrated hosts live longer, suggesting that parasites can redirect resources normally used in reproduction to increase host lifespan, lengthening their transmission phase, and ii) that multiple infections increase virulence, here in terms of non-recovery and host castration.

Inference with viral quasispecies diversity indices: clonal and NGS approaches

Given the inherent dynamics of a viral quasispecies, we are often interested in the comparison of diversity indices of sequential samples of a patient, or in the comparison of diversity indices of virus in groups of patients in a treated versus control design. It is then important to make sure that the diversity measures from each sample may be compared with no bias and within a consistent statistical framework. In the present report, we review some indices often used as measures for viral quasispecies complexity and provide means for statistical inference, applying procedures taken from the ecology field. In particular, we examine the Shannon entropy and the mutation frequency, and we discuss the appropriateness of different normalization methods of the Shannon entropy found in the literature. By taking amplicons ultra-deep pyrosequencing (UDPS) raw data as a surrogate of a real hepatitis C virus viral population, we study through in-silico sampling the statistical properties of these indices under two methods of viral quasispecies sampling, classical cloning followed by Sanger sequencing (CCSS) and next-generation sequencing (NGS) such as UDPS. We propose solutions specific to each of the two sampling methods-CCSS and NGS-to guarantee statistically conforming conclusions as free of bias as possible.

CONTACT

josep.gregori@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Competition, virulence, host body mass and the diversification of macro-parasites

Adaptive speciation has been much debated in recent years, with a strong emphasis on how competition can lead to the diversification of ecological and sexual traits. Surprisingly, little attention has been paid to this evolutionary process to explain intrahost diversification of parasites. We expanded the theory of competitive speciation to look at the effect of key features of the parasite lifestyle, namely fragmentation, aggregation and virulence, on the conditions and rate of sympatric speciation under the standard 'pleiotropic scenario'. The conditions for competitive speciation were found similar to those for non-parasite species, but not the rate of diversification. Adaptive evolution proceeds faster in highly fragmented parasite populations and for weakly aggregated and virulent parasites. Combining these theoretical results with standard empirical allometric relationships, we showed that parasite diversification can be faster in host species of intermediate body mass. The increase in parasite load with body mass, indeed, fuels evolution by increasing mutants production, but because of the deleterious effect of virulence, it simultaneously weakens selection for resource specialization. Those two antagonistic effects lead to optimal parasite burden and host body mass for diversification. Data on the diversity of fishes' gills parasites were found consistent with the existence of such optimum.

In vitro infection of pupae with Israeli acute paralysis virus suggests disturbance of transcriptional homeostasis in honey bees (Apis mellifera)

The ongoing decline of honey bee health worldwide is a serious economic and ecological concern. One major contributor to the decline are pathogens, including several honey bee viruses. However, information is limited on the biology of bee viruses and molecular interactions with their hosts. An experimental protocol to test these systems was developed, using injections of Israeli Acute Paralysis Virus (IAPV) into honey bee pupae reared ex-situ under laboratory conditions. The infected pupae developed pronounced but variable patterns of disease. Symptoms varied from complete cessation of development with no visual evidence of disease to rapid darkening of a part or the entire body. Considerable differences in IAPV titer dynamics were observed, suggesting significant variation in resistance to IAPV among and possibly within honey bee colonies. Thus, selective breeding for virus resistance should be possible. Gene expression analyses of three separate experiments suggest IAPV disruption of transcriptional homeostasis of several fundamental cellular functions, including an up-regulation of the ribosomal biogenesis pathway. These results provide first insights into the mechanisms of IAPV pathogenicity. They mirror a transcriptional survey of honey bees afflicted with Colony Collapse Disorder and thus support the hypothesis that viruses play a critical role in declining honey bee health.

Microevolution of Puumala hantavirus during a complete population cycle of its host, the bank vole (Myodes glareolus)

Microevolution of Puumala hantavirus (PUUV) was studied throughout a population cycle of its host, the bank vole (Myodes glareolus). We monitored PUUV variants circulating in the host population in Central Finland over a five-year period that included two peak-phases and two population declines. Of 1369 bank voles examined, 360 (26.3%) were found infected with PUUV. Partial sequences of each of the three genome segments were recovered (approx. 12% of PUUV genome) from 356 bank voles. Analyses of these sequences disclosed the following features of PUUV evolution: 1) nucleotide substitutions are mostly silent and deduced amino acid changes are mainly conservative, suggesting stabilizing selection at the protein level; 2) the three genome segments accumulate mutations at a different rate; 3) some of the circulating PUUV variants are frequently observed while others are transient; 4) frequently occurring PUUV variants are composed of the most abundant segment genotypes (copious) and new transient variants are continually generated; 5) reassortment of PUUV genome segments occurs regularly and follows a specific pattern of segments association; 6) prevalence of reassortant variants oscillates with season and is higher in the autumn than in the spring; and 7) reassortants are transient, i.e., they are not competitively superior to their parental variants. Collectively, these observations support a quasi-neutral mode of PUUV microevolution with a steady generation of transient variants, including reassortants, and preservation of a few preferred genotypes.

Evolution of virulence in emerging epidemics

Theory predicts that selection for pathogen virulence and horizontal transmission is highest at the onset of an epidemic but decreases thereafter, as the epidemic depletes the pool of susceptible hosts. We tested this prediction by tracking the competition between the latent bacteriophage λ and its virulent mutant λcI857 throughout experimental epidemics taking place in continuous cultures of Escherichia coli. As expected, the virulent λcI857 is strongly favored in the early stage of the epidemic, but loses competition with the latent virus as prevalence increases. We show that the observed transient selection for virulence and horizontal transmission can be fully explained within the framework of evolutionary epidemiology theory. This experimental validation of our predictions is a key step towards a predictive theory for the evolution of virulence in emerging infectious diseases.

Intrahost mechanisms governing emergence of resistance-breaking variants of Potato virus Y

The emergence of resistance breaking (RB) variants starting from the avirulent Potato virus Y NN strain (PVY(NN)) was analyzed after imposing different selective host constraints. Tobacco resistance to PVY(NN) is conferred by va in both NC745 and VAM genotypes, but VAM carries an extra resistance gene, va2. RB-variants emerged only in NC745 and unexpectedly accumulated higher in the original host, tobacco B21, than the parental PVY(NN). However, the recovery of RB-variants was interfered by PVY(NN) in mixed infections. Further analysis indicated that RB-variants also arose in tobacco VAM, but they were limited to subliminal local infections. Their inability to breakout was associated with absence of a mutational adaptation in the viral VPg gene, which implied a loss of fitness in tobacco B21. Altogether, the emergence of RB-variants was conditioned by inherited host constraints, interference by co-infecting avirulent virus genotypes, and fitness tradeoff of virus adaptations.

Phylodynamic analysis of the emergence and epidemiological impact of transmissible defective dengue viruses

Intra-host sequence data from RNA viruses have revealed the ubiquity of defective viruses in natural viral populations, sometimes at surprisingly high frequency. Although defective viruses have long been known to laboratory virologists, their relevance in clinical and epidemiological settings has not been established. The discovery of long-term transmission of a defective lineage of dengue virus type 1 (DENV-1) in Myanmar, first seen in 2001, raised important questions about the emergence of transmissible defective viruses and their role in viral epidemiology. By combining phylogenetic analyses and dynamical modeling, we investigate how evolutionary and ecological processes at the intra-host and inter-host scales shaped the emergence and spread of the defective DENV-1 lineage. We show that this lineage of defective viruses emerged between June 1998 and February 2001, and that the defective virus was transmitted primarily through co-transmission with the functional virus to uninfected individuals. We provide evidence that, surprisingly, this co-transmission route has a higher transmission potential than transmission of functional dengue viruses alone. Consequently, we predict that the defective lineage should increase overall incidence of dengue infection, which could account for the historically high dengue incidence reported in Myanmar in 2001–2002. Our results show the unappreciated potential for defective viruses to impact the epidemiology of human pathogens, possibly by modifying the virulence-transmissibility trade-off, or to emerge as circulating infections in their own right. They also demonstrate that interactions between viral variants, such as complementation, can open new pathways to viral emergence.

Defective viruses are viral particles with genetic mutations or deletions that eliminate essential functions, so that they cannot complete their life cycles independently. They can reproduce only by co-infecting host cells with functional viruses and ‘borrowing’ their functional elements. Defective viruses have been observed for many human pathogens, but they have not been thought to impact epidemiological outcomes. Recently it was reported that a lineage of defective dengue virus spread through humans and mosquitoes in Myanmar for at least 18 months in 2001–2002. In this study, we investigate the emergence and epidemiological impact of this defective lineage by combining genetic sequence analyses with mathematical models. We show that the defective lineage emerged from circulating dengue viruses between June 1998 and February 2001, and that it spreads because—surprisingly—its presence causes functional dengue viruses to transmit more efficiently. Our model shows that this would cause a substantial rise in total dengue infections, consistent with historically high levels of dengue cases reported in Myanmar during 2001 and 2002. Our study yields new insights into the biology of dengue virus, and demonstrates a previously unappreciated potential for defective viruses to impact the epidemiology of infectious diseases.

Virus replication as a phenotypic version of polynucleotide evolution

In this paper, we revisit and adapt to viral evolution an approach based on the theory of branching process advanced by Demetrius et al. (Bull. Math. Biol. 46:239-262, 1985), in their study of polynucleotide evolution. By taking into account beneficial effects, we obtain a non-trivial multivariate generalization of their single-type branching process model. Perturbative techniques allows us to obtain analytical asymptotic expressions for the main global parameters of the model, which lead to the following rigorous results: (i) a new criterion for "no sure extinction", (ii) a generalization and proof, for this particular class of models, of the lethal mutagenesis criterion proposed by Bull et al. (J. Virol. 18:2930-2939, 2007), (iii) a new proposal for the notion of relaxation time with a quantitative prescription for its evaluation, (iv) the quantitative description of the evolution of the expected values in four distinct "stages": extinction threshold, lethal mutagenesis, stationary "equilibrium", and transient. Finally, based on these quantitative results, we are able to draw some qualitative conclusions.

Sequencing approach to analyze the role of quasispecies for classical swine fever

Classical swine fever virus (CSFV) is a positive-sense RNA virus with a high degree of genetic variability among isolates. High diversity is also found in virulence, with strains covering the complete spectrum from avirulent to highly virulent. The underlying genetic determinants are far from being understood. Since RNA polymerases of RNA viruses lack any proof-reading activity, different genome variations called haplotypes, occur during replication. A set of haplotypes is referred to as a viral quasispecies. Genetic variability can be a fitness advantage through facilitating of a more effective escape from the host immune response. In order to investigate the correlation of quasispecies composition and virulence in vivo, we analyzed next-generation sequencing data of CSFV isolates of varying virulence. Viral samples from pigs infected with the highly virulent isolates "Koslov" and "Brescia" showed higher quasispecies diversity and more nucleotide variability, compared to samples of pigs infected with low and moderately virulent isolates.

Approaches to define the viral genetic basis of classical swine fever virus virulence

Classical swine fever (CSF), a highly contagious disease of pigs caused by the classical swine fever virus (CSFV), can lead to important economic losses in the pig industry. Numerous CSFV isolates with various degrees of virulence have been isolated worldwide, ranging from low virulent strains that do not result in any apparent clinical signs to highly virulent strains that cause a severe peracute hemorrhagic fever with nearly 100% mortality. Knowledge of the molecular determinants of CSFV virulence is an important issue for effective disease control and development of safe and effective marker vaccines. In this review, the latest studies in the field of CSFV virulence are discussed. The topic of virulence is addressed from different angles; nonconventional approaches like codon pair usage and quasispecies are considered. Future research approaches in the field of CSFV virulence are proposed.

Challenges and opportunities in estimating viral genetic diversity from next-generation sequencing data

Many viruses, including the clinically relevant RNA viruses HIV (human immunodeficiency virus) and HCV (hepatitis C virus), exist in large populations and display high genetic heterogeneity within and between infected hosts. Assessing intra-patient viral genetic diversity is essential for understanding the evolutionary dynamics of viruses, for designing effective vaccines, and for the success of antiviral therapy. Next-generation sequencing (NGS) technologies allow the rapid and cost-effective acquisition of thousands to millions of short DNA sequences from a single sample. However, this approach entails several challenges in experimental design and computational data analysis. Here, we review the entire process of inferring viral diversity from sample collection to computing measures of genetic diversity. We discuss sample preparation, including reverse transcription and amplification, and the effect of experimental conditions on diversity estimates due to in vitro base substitutions, insertions, deletions, and recombination. The use of different NGS platforms and their sequencing error profiles are compared in the context of various applications of diversity estimation, ranging from the detection of single nucleotide variants (SNVs) to the reconstruction of whole-genome haplotypes. We describe the statistical and computational challenges arising from these technical artifacts, and we review existing approaches, including available software, for their solution. Finally, we discuss open problems, and highlight successful biomedical applications and potential future clinical use of NGS to estimate viral diversity.

Competition-colonization trade-off promotes coexistence of low-virulence viral strains

RNA viruses exist as genetically diverse populations displaying a range of virulence degrees. The evolution of virulence in viral populations is, however, poorly understood. On the basis of the experimental observation of an RNA virus clone in cell culture diversifying into two subpopulations of different virulence, we study the dynamics of mutating virus populations with varying virulence. We introduce a competition-colonization trade-off into standard mathematical models of intra-host viral infection. Colonizers are fast-spreading virulent strains, whereas the competitors are less-virulent variants but more successful within co-infected cells. We observe a two-step dynamics of the population. Early in the infection, the population is dominated by colonizers, which later are outcompeted by competitors. Our simulations suggest the existence of steady state in which all virulence classes coexist but are dominated by the most competitive ones. This equilibrium implies collective virulence attenuation in the population, in contrast to previous models predicting evolution of the population towards increased virulence.

Viral population dynamics and virulence thresholds

Viral factors and host barriers influence virally induced disease, and asymptomatic versus symptomatic infection is governed by a 'virulence threshold'. Understanding modulation of virulence thresholds could lend insight into disease outcome and aid in rational therapeutic and vaccine design. RNA viruses are an excellent system to study virulence thresholds in the context of quasispecies population dynamics. RNA viruses have high error frequencies and our understanding of viral population dynamics has been shaped by quasispecies evolutionary theory. In turn, research using RNA viruses as replicons with short generation times and high mutation rates has been an invaluable tool to test models of quasispecies theory. The challenge and new frontier of RNA virus population dynamics research is to combine multiple theoretical models and experimental data to describe viral population behavior as it changes, moving within and between hosts, to predict disease and pathogen emergence. Several excellent studies have begun to undertake this challenge using novel approaches.

Mutagenesis-mediated decrease of pathogenicity as a feature of the mutant spectrum of a viral population

BACKGROUND

RNA virus populations are heterogeneous ensembles of closely related genomes termed quasispecies. This highly complex distribution of variants confers important properties to RNA viruses and influences their pathogenic behavior. It has been hypothesized that increased mutagenesis of viral populations, by treatment with mutagenic agents, can induce alterations in the pathogenic potential of a virus population. In this work we investigate whether mutagenized foot-and-mouth disease virus (FMDV) populations display changes in their virulence in mice.

METHODOLOGY AND PRINCIPAL FINDINGS

FMDV C-S8c1 was passaged in BHK cells in the presence of the mutagenic agent ribavirin. Decline in viral titer and viral RNA progeny was observed in the first passage, fluctuating around a constant value thereafter. Hence, the specific infectivity remained stable during the passages. The viral population harvested from passage 9 (P9 R) showed decreased virulence in mice, with a lethal dose 50 (LD(50)) >10(4) PFU, as compared with LD(50) of 50 PFU of the parental population FMDV C-S8c1. This decrease in virulence was associated to a 20-fold increase in the mutation frequency of the P9 R population with respect to C-S8c1. Interestingly, individual biological clones isolated from the attenuated population P9 R were as virulent as the parental virus C-S8c1. Furthermore, a mixed population of C-S8c1 and P9 R was inoculated into mice and showed decreased virulence as compared to C-S8c1, suggesting that population P9 R is able to suppress the virulent phenotype of C-S8c1.

CONCLUSION

Ribavirin-mediated mutagenesis of an FMDV population resulted in attenuation in vivo, albeit a large proportion of its biological clones displayed a highly virulent phenotype. These results, together with the suppression of C-S8c1 by mutagenized P9 R population, document a suppressive effect of mutagenized viral quasispecies in vivo, and suggest novel approaches to the treatment and prevention of viral diseases.

Viral quasispecies evolution

Evolution of RNA viruses occurs through disequilibria of collections of closely related mutant spectra or mutant clouds termed viral quasispecies. Here we review the origin of the quasispecies concept and some biological implications of quasispecies dynamics. Two main aspects are addressed: (i) mutant clouds as reservoirs of phenotypic variants for virus adaptability and (ii) the internal interactions that are established within mutant spectra that render a virus ensemble the unit of selection. The understanding of viruses as quasispecies has led to new antiviral designs, such as lethal mutagenesis, whose aim is to drive viruses toward low fitness values with limited chances of fitness recovery. The impact of quasispecies for three salient human pathogens, human immunodeficiency virus and the hepatitis B and C viruses, is reviewed, with emphasis on antiviral treatment strategies. Finally, extensions of quasispecies to nonviral systems are briefly mentioned to emphasize the broad applicability of quasispecies theory.

Mutagenesis-mediated virus extinction: virus-dependent effect of viral load on sensitivity to lethal defection

BACKGROUND

Lethal mutagenesis is a transition towards virus extinction mediated by enhanced mutation rates during viral genome replication, and it is currently under investigation as a potential new antiviral strategy. Viral load and virus fitness are known to influence virus extinction. Here we examine the effect or the multiplicity of infection (MOI) on progeny production of several RNA viruses under enhanced mutagenesis.

RESULTS

The effect of the mutagenic base analogue 5-fluorouracil (FU) on the replication of the arenavirus lymphocytic choriomeningitis virus (LCMV) can result either in inhibition of progeny production and virus extinction in infections carried out at low multiplicity of infection (MOI), or in a moderate titer decrease without extinction at high MOI. The effect of the MOI is similar for LCMV and vesicular stomatitis virus (VSV), but minimal or absent for the picornaviruses foot-and-mouth disease virus (FMDV) and encephalomyocarditis virus (EMCV). The increase in mutation frequency and Shannon entropy (mutant spectrum complexity) as a result of virus passage in the presence of FU was more accentuated at low MOI for LCMV and VSV, and at high MOI for FMDV and EMCV. We present an extension of the lethal defection model that agrees with the experimental results.

CONCLUSIONS

(i) Low infecting load favoured the extinction of negative strand viruses, LCMV or VSV, with an increase of mutant spectrum complexity. (ii) This behaviour is not observed in RNA positive strand viruses, FMDV or EMCV. (iii) The accumulation of defector genomes may underlie the MOI-dependent behaviour. (iv) LCMV coinfections are allowed but superinfection is strongly restricted in BHK-21 cells. (v) The dissimilar effects of the MOI on the efficiency of mutagenic-based extinction of different RNA viruses can have implications for the design of antiviral protocols based on lethal mutagenesis, presently under development.

Biased mutational pattern and quasispecies hypothesis in H5N1 virus

Like other RNA viruses, influenza viruses are subject to high mutation rates. Carrying segmented RNA genomes, their genetic variability is even higher. We aimed at analyzing the mutational events occurring during the infection of chickens by the Highly Pathogenic Avian Influenza (HPAI) H5N1 virus. We therefore studied the different sequences of two surface proteins, hemagglutinin (HA) and neuraminidase (NA), as well as two internal proteins, PB2 and NS. Three organs (lung, spleen, brain) were obtained from a chicken, experimentally infected with a lethal dose of HPAI H5N1 virus. Cloning these PCR fragments enabled us to investigate the mutations undergone by the virus after several replicative cycles. The first outcome is the presence of a strong mutational bias, resembling host-driven ADAR1 adenosine deamination, which is responsible for 81% of all mutations. Whereas the frequency of RNA dependent RNA polymerase-related mutations is compatible with the survival of the virus, the ADAR1-like activity usually strongly increases the mutation frequency into a level of “error catastrophe” in theory incompatible with virus survival. Nevertheless, the virus was successfully infective. HPAI H5N1 virus displayed traits in agreement with the quasispecies theory. The role of this quasispecies structure in successful infection and the superposition with the ADAR1-like response is discussed.

Quasispecies as a matter of fact: viruses and beyond

We review the origins of the quasispecies concept and its relevance for RNA virus evolution, viral pathogenesis and antiviral treatment strategies. We emphasize a critical point of quasispecies that refers to genome collectivities as the unit of selection, and establish parallels between RNA viruses and some cellular systems such as bacteria and tumor cells. We refer also to tantalizing new observations that suggest quasispecies behavior in prions, perhaps as a result of the same quantum-mechanical indeterminations that underlie protein conformation and error-prone replication in genetic systems. If substantiated, these observations with prions could lead to new research on the structure-function relationship of non-nucleic acid biological molecules.

Lethal mutagenesis of foot-and-mouth disease virus involves shifts in sequence space

Lethal mutagenesis or virus transition into error catastrophe is an antiviral strategy that aims at extinguishing a virus by increasing the viral mutation rates during replication. The molecular basis of lethal mutagenesis is largely unknown. Previous studies showed that a critical substitution in the foot-and-mouth disease virus (FMDV) polymerase was sufficient to allow the virus to escape extinction through modulation of the transition types induced by the purine nucleoside analogue ribavirin. This substitution was not detected in mutant spectra of FMDV populations that had not replicated in the presence of ribavirin, using standard molecular cloning and nucleotide sequencing. Here we selectively amplify and analyze low-melting-temperature cDNA duplexes copied from FMDV genome populations passaged in the absence or presence of ribovirin Hypermutated genomes with high frequencies of A and U were present in both ribavirin -treated and untreated populations, but the major effect of ribavirin mutagenesis was to accelerate the occurrence of AU-rich mutant clouds during the early replication rounds of the virus. The standard FMDV quasispecies passaged in the absence of ribavirin included the salient transition-modulating, ribavirin resistance mutation, whose frequency increased in populations treated with ribavirin. Thus, even nonmutagenized FMDV quasispecies include a deep, mutationally biased portion of sequence space, in support of the view that the virus replicates close to the error threshold for maintenance of genetic information.

Competition-colonization dynamics: An ecology approach to quasispecies dynamics and virulence evolution in RNA viruses

A single and purified clone of foot-and-mouth disease virus diversified in cell culture into two subpopulations that were genetically distinct. The subpopulation with higher virulence was a minority and was suppressed by the dominant but less virulent one. These two populations follow the competitioncolonization dynamics described in ecology. Virulent viruses can be regarded as colonizers because they killed the cells faster and they spread faster. The attenuated subpopulation resembles competitors because of its higher replication efficiency in coinfected cells. Our results suggest a new model for the evolution of virulence which is based on interactions between components of the quasispecies. Competition between viral mutants takes place at two levels, intracellular competition and competition for new cells. The two strategies are subjected to densitydependent selection.

Models of RNA virus evolution and their roles in vaccine design

Viruses are fast evolving pathogens that continuously adapt to the highly variable environments they live and reproduce in. Strategies devoted to inhibit virus replication and to control their spread among hosts need to cope with these extremely heterogeneous populations and with their potential to avoid medical interventions. Computational techniques such as phylogenetic methods have broadened our picture of viral evolution both in time and space, and mathematical modeling has contributed substantially to our progress in unraveling the dynamics of virus replication, fitness, and virulence. Integration of multiple computational and mathematical approaches with experimental data can help to predict the behavior of viral pathogens and to anticipate their escape dynamics. This piece of information plays a critical role in some aspects of vaccine development, such as viral strain selection for vaccinations or rational attenuation of viruses. Here we review several aspects of viral evolution that can be addressed quantitatively, and we discuss computational methods that have the potential to improve vaccine design.

Mutant spectra in virus behavior

RNA viruses replicate as complex mutant spectra, also termed ‘mutant clouds’, known as viral quasispecies. While this is a widely observed viral population structure, it is less known that a number of biologically relevant features of this important group of viral pathogens depend on (or are strongly influenced by) the complexity and composition of mutant spectra. Among them, fitness increase or decrease depending on intrapopulation complementation or interference, selection triggered by memory genomes, pathogenic potential of viruses, disease evolution and the response to antiviral treatments. Quasispecies represent the recognition of complex behavior in viruses, and it is an oversimplification to equate such a population structure with the classic polymorphism of population biology. Darwinian principles acting on genome collectivities that replicate with high error rates provide a unique population structure prone to flexible and largely unpredictable behavior.

A multi-step process of viral adaptation to a mutagenic nucleoside analogue by modulation of transition types leads to extinction-escape

Resistance of viruses to mutagenic agents is an important problem for the development of lethal mutagenesis as an antiviral strategy. Previous studies with RNA viruses have documented that resistance to the mutagenic nucleoside analogue ribavirin (1-β-D-ribofuranosyl-1-H-1,2,4-triazole-3-carboxamide) is mediated by amino acid substitutions in the viral polymerase that either increase the general template copying fidelity of the enzyme or decrease the incorporation of ribavirin into RNA. Here we describe experiments that show that replication of the important picornavirus pathogen foot-and-mouth disease virus (FMDV) in the presence of increasing concentrations of ribavirin results in the sequential incorporation of three amino acid substitutions (M296I, P44S and P169S) in the viral polymerase (3D). The main biological effect of these substitutions is to attenuate the consequences of the mutagenic activity of ribavirin —by avoiding the biased repertoire of transition mutations produced by this purine analogue—and to maintain the replicative fitness of the virus which is able to escape extinction by ribavirin. This is achieved through alteration of the pairing behavior of ribavirin-triphosphate (RTP), as evidenced by in vitro polymerization assays with purified mutant 3Ds. Comparison of the three-dimensional structure of wild type and mutant polymerases suggests that the amino acid substitutions alter the position of the template RNA in the entry channel of the enzyme, thereby affecting nucleotide recognition. The results provide evidence of a new mechanism of resistance to a mutagenic nucleoside analogue which allows the virus to maintain a balance among mutation types introduced into progeny genomes during replication under strong mutagenic pressure.

Viruses that have RNA as genetic material include many important human, animal and plant pathogens. A new strategy against RNA viruses consists in using mutagenic nucleotides. The objective is to provoke an excessive number of mutations, to deteriorate the viral functions to the point that the virus can not survive. One of the mutagens used in research on lethal mutagenesis is ribavirin, extensively employed in clinical practice. Unfortunately, viral mutants that are resistant to ribavirin have been selected, thus facilitating escape from lethal mutagenesis. Here we describe a new mechanism by which foot-and-mouth disease virus (FMDV) can become resistant to ribavirin. Amino acid changes in the viral polymerase, selected by ribavirin, are able to modify the types of mutations produced in the presence of ribavirin. Biochemical data indicate that the alteration of the enzyme changes the pairing behavior of ribavirin, avoiding the production of an excess of some types of mutations, supporting the hypothesis that an unbalanced mutation repertoire is detrimental to the virus. Thus, this new mechanism of resistance to ribavirin is based not as much in limiting the number of mutations in the virus genetic material but in ensuring an equilibrium among different types of mutations that favors viral survival.

Virulence correlates with fitness in vivo for two M group genotypes of Infectious hematopoietic necrosis virus (IHNV)

The nature of the association between viral fitness and virulence remains elusive in vertebrate virus systems, partly due to a lack of in vivo experiments using statistically sufficient numbers of replicate hosts. We examined the relationship between virulence and fitness in Infectious hematopoietic necrosis virus (IHNV), in vivo, in intact living rainbow trout. Trout were infected with a high or low virulence genotype of M genogroup IHNV, or a mixture of the two genotypes, so as to calculate relative fitness and the effect of a competition environment on fitness. Fitness was measured as total viral load in the host at time of peak viral density, quantified by genotype-specific quantitative RT-PCR (qRT-PCR). The more virulent IHNV genotype reached higher densities in both single and mixed infections. There was no effect of competition on the performance of either genotype. Our results suggest a positive link between IHNV genotype fitness and virulence.