Search

When viruses have segmented genomes, the set of frequencies describing the abundance of segments is called the genome formula. The genome formula is often unbalanced and highly variable for both segmented and multipartite viruses. A growing number of studies are quantifying the genome formula to measure its effects on infection and to consider its ecological and evolutionary implications. Different approaches have been reported for analyzing genome formula data, including qualitative description, applying standard statistical tests such as ANOVA, and customized analyses. However, these approaches have different shortcomings, and test assumptions are often unmet, potentially leading to erroneous conclusions. Here, we address these challenges, leading to a threefold contribution. First, we propose a simple metric for analyzing genome formula variation: the genome formula distance. We describe the properties of this metric and provide a framework for understanding metric values. Second, we explain how this metric can be applied for different purposes, including testing for genome-formula differences and comparing observations to a reference genome formula value. Third, we re-analyze published data to illustrate the applications and weigh the evidence for previous conclusions. Our re-analysis of published datasets confirms many previous results but also provides evidence that the genome formula can be carried over from the inoculum to the virus population in a host. The simple procedures we propose contribute to the robust and accessible analysis of genome-formula data.

Adaptive evolutionary processes are constrained by the availability of mutations which cause a fitness benefit and together make up the fitness landscape, which maps genotype space onto fitness under specified conditions. Experimentally derived fitness landscapes have demonstrated a predictability to evolution by identifying limited "mutational routes" that evolution by natural selection may take between low and high-fitness genotypes. However, such studies often utilize indirect measures to determine fitness. We estimated the competitive fitness of mutants relative to all single-mutation neighbors to describe the fitness landscape of three mutations in a β-lactamase enzyme. Fitness assays were performed at sublethal concentrations of the antibiotic cefotaxime in a structured and unstructured environment. In the unstructured environment, the antibiotic selected for higher-resistance types-but with an equivalent fitness for a subset of mutants, despite substantial variation in resistance-resulting in a stratified fitness landscape. In contrast, in a structured environment with a low antibiotic concentration, antibiotic-susceptible genotypes had a relative fitness advantage, which was associated with antibiotic-induced filamentation. These results cast doubt that highly resistant genotypes have a unique selective advantage in environments with subinhibitory concentrations of antibiotics and demonstrate that direct fitness measures are required for meaningful predictions of the accessibility of evolutionary routes. IMPORTANCE The evolution of antibiotic-resistant bacterial populations underpins the ongoing antibiotic resistance crisis. We aim to understand how antibiotic-degrading enzymes can evolve to cause increased resistance, how this process is constrained, and whether it can be predictable. To this end, competition experiments were performed with a combinatorially complete set of mutants of a β-lactamase gene subject to subinhibitory concentrations of the antibiotic cefotaxime. While some mutations confer on their hosts high resistance to cefotaxime, in competition these mutations do not always confer a selective advantage. Specifically, high-resistance mutants had equivalent fitnesses despite different resistance levels and even had selective disadvantages under conditions involving spatial structure. Together, our findings suggest that the relationship between resistance level and fitness at subinhibitory concentrations is complex; predicting the evolution of antibiotic resistance requires knowledge of the conditions that select for resistant genotypes and the selective advantage evolved types have over their predecessors., Competing Interests: The authors declare no conflict of interest.

The physiological performance of organisms depends on their environmental context, resulting in performance-response curves along environmental gradients. Parasite performance-response curves are generally expected to be broader than those of their hosts due to shorter generation times and hence faster adaptation. However, certain environmental conditions may limit parasite performance more than that of the host, thereby providing an environmental refuge from disease. Thermal disease refuges have been extensively studied in response to climate warming, but other environmental factors may also provide environmental disease refuges which, in turn, respond to global change. Here, we (1) showcase laboratory and natural examples of refuges from parasites along various environmental gradients, and (2) provide hypotheses on how global environmental change may affect these refuges. We strive to synthesize knowledge on potential environmental disease refuges along different environmental gradients including salinity and nutrients, in both natural and food-production systems. Although scaling up from single host-parasite relationships along one environmental gradient to their interaction outcome in the full complexity of natural environments remains difficult, integrating host and parasite performance-response can serve to formulate testable hypotheses about the variability in parasitism outcomes and the occurrence of environmental disease refuges under current and future environmental conditions., (© 2023 The Authors. Ecology published by Wiley Periodicals LLC on behalf of The Ecological Society of America.)

Viruses show great diversity in their genome organization. Multipartite viruses package their genome segments into separate particles, most or all of which are required to initiate infection in the host cell. The benefits of such seemingly inefficient genome organization are not well understood. One hypothesised benefit of multipartition is that it allows for flexible changes in gene expression by altering the frequency of each genome segment in different environments, such as encountering different host species. The ratio of the frequency of segments is termed the genome formula (GF). Thus far, formal studies quantifying the GF have been performed for well-characterised virus-host systems in experimental settings using RT-qPCR. However, to understand GF variation in natural populations or novel virus-host systems, a comparison of several methods for GF estimation including high-throughput sequencing (HTS) based methods is needed. Currently, it is unclear how HTS-methods compare a golden standard, such as RT-qPCR. Here we show a comparison of multiple GF quantification methods (RT-qPCR, RT-digital PCR, Illumina RNAseq and Nanopore direct RNA sequencing) using three host plants (Nicotiana tabacum, Nicotiana benthamiana, and Chenopodium quinoa) infected with cucumber mosaic virus (CMV), a tripartite RNA virus. Our results show that all methods give roughly similar results, though there is a significant method effect on genome formula estimates. While the RT-qPCR and RT-dPCR GF estimates are congruent, the GF estimates from HTS methods deviate from those found with PCR. Our findings emphasize the need to tailor the GF quantification method to the experimental aim, and highlight that it may not be possible to compare HTS and PCR-based methods directly. The difference in results between PCR-based methods and HTS highlights that the choice of quantification technique is not trivial., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

The genomes of RNA viruses may be monopartite or multipartite, and sub-genomic particles such as defective RNAs (D RNAs) or satellite RNAs (satRNAs) can be associated with some of them. D RNAs are small, deletion mutants of a virus that have lost essential functions for independent replication, encapsidation and/or movement. D RNAs are common elements associated with human and animal viruses, and they have been described for numerous plant viruses so far. Over 30 years of studies on D RNAs allow for some general conclusions to be drawn. First, the essential condition for D RNA formation is prolonged passaging of the virus at a high cellular multiplicity of infection (MOI) in one host. Second, recombination plays crucial roles in D RNA formation. Moreover, during virus propagation, D RNAs evolve, and the composition of the particle depends on, e.g., host plant, virus isolate or number of passages. Defective RNAs are often engaged in transient interactions with full-length viruses-they can modulate accumulation, infection dynamics and virulence, and are widely used, i.e., as a tool for research on cis -acting elements crucial for viral replication. Nevertheless, many questions regarding the generation and role of D RNAs in pathogenesis remain open. In this review, we summarise the knowledge about D RNAs of plant viruses obtained so far.

Bunyaviruses lack a specific mechanism to ensure the incorporation of a complete set of genome segments into each virion, explaining the generation of incomplete virus particles lacking one or more genome segments. Such incomplete virus particles, which may represent the majority of particles produced, are generally considered to interfere with virus infection and spread. Using the three-segmented arthropod-borne Rift Valley fever virus as a model bunyavirus, we here show that two distinct incomplete virus particle populations unable to spread autonomously are able to efficiently complement each other in both mammalian and insect cells following co-infection. We further show that complementing incomplete virus particles can co-infect mosquitoes, resulting in the reconstitution of infectious virus that is able to disseminate to the mosquito salivary glands. Computational models of infection dynamics predict that incomplete virus particles can positively impact virus spread over a wide range of conditions, with the strongest effect at intermediate multiplicities of infection. Our findings suggest that incomplete particles may play a significant role in within-host spread and between-host transmission, reminiscent of the infection cycle of multipartite viruses., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2022 Bermúdez-Méndez et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Mutation rates are of key importance for understanding evolutionary processes and predicting their outcomes. Empirical mutation rate estimates are available for a number of RNA viruses, but few are available for DNA viruses, which tend to have larger genomes. Whilst some viruses have very high mutation rates, lower mutation rates are expected for viruses with large genomes to ensure genome integrity. Alphabaculoviruses are insect viruses with large genomes and often have high levels of polymorphism, suggesting high mutation rates despite evidence of proofreading activity by the replication machinery. Here, we report an empirical estimate of the mutation rate per base per strand copying (s/n/r) of Autographa californica multiple nucleopolyhedrovirus (AcMNPV). To avoid biases due to selection, we analyzed mutations that occurred in a stable, non-functional genomic insert after five serial passages in Spodoptera exigua larvae. Our results highlight that viral demography and the stringency of mutation calling affect mutation rate estimates, and that using a population genetic simulation model to make inferences can mitigate the impact of these processes on estimates of mutation rate. We estimated a mutation rate of μ = 1×10-7 s/n/r when applying the most stringent criteria for mutation calling, and estimates of up to μ = 5×10-7 s/n/r when relaxing these criteria. The rates at which different classes of mutations accumulate provide good evidence for neutrality of mutations occurring within the inserted region. We therefore present a robust approach for mutation rate estimation for viruses with stable genomes, and strong evidence of a much lower alphabaculovirus mutation rate than supposed based on the high levels of polymorphism observed., Competing Interests: The authors have declared that no competing interests exist.

Mutations with large fitness benefits and mutations occurring at high rates may both cause parallel evolution, but their contribution is predicted to depend on population size. Moreover, high-rate and large-benefit mutations may have different long-term adaptive consequences. We show that small and 100-fold larger bacterial populations evolve resistance to a β-lactam antibiotic by using similar numbers, but different types of mutations. Small populations frequently substitute similar high-rate structural variants and loss-of-function point mutations, including the deletion of a low-activity β-lactamase, and evolve modest resistance levels. Large populations more often use low-rate, large-benefit point mutations affecting the same targets, including mutations activating the β-lactamase and other gain-of-function mutations, leading to much higher resistance levels. Our results demonstrate the separation by clonal interference of mutation classes with divergent adaptive consequences, causing a shift from high-rate to large-benefit mutations with increases in population size., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Biodiversity can reduce or increase disease transmission. These divergent effects suggest that community composition rather than diversity per se determines disease transmission. In natural plant communities, little is known about the functional roles of neighbouring plant species in belowground disease transmission. Here, we experimentally investigated disease transmission of a fungal root pathogen (Rhizoctonia solani) in two focal plant species in combinations with four neighbour species of two ages. We developed stochastic models to test the relative importance of two transmission-modifying mechanisms: (1) infected hosts serve as nutrient supply to increase hyphal growth, so that successful disease transmission is self-reinforcing; and (2) plant resistance increases during plant development. Neighbouring plants either reduced or increased disease transmission in the focal plants. These effects depended on neighbour age, but could not be explained by a simple dichotomy between hosts and nonhost neighbours. Model selection revealed that both transmission-modifying mechanisms are relevant and that focal host-neighbour interactions changed which mechanisms steered disease transmission rate. Our work shows that neighbour-induced shifts in the importance of these mechanisms across root networks either make or break disease transmission chains. Understanding how diversity affects disease transmission thus requires integrating interactions between focal and neighbour species and their pathogens., (© 2021 The Authors. New Phytologist © 2021 New Phytologist Foundation.)

Understanding under which conditions conjugative plasmids encoding antibiotic resistance can invade bacterial communities in the gut is of particular interest to combat the spread of antibiotic resistance within and between animals and humans. We extended a one-compartment model of conjugation to a two-compartment model, to analyse how differences in plasmid dynamics in the gut lumen and at the gut wall affect the invasion of plasmids. We compared scenarios with one and two compartments, different migration rates between the lumen and wall compartments, and different population dynamics. We focused on the effect of attachment and detachment rates on plasmid dynamics, explicitly describing pair formation followed by plasmid transfer in the pairs. The parameter space allowing plasmid invasion in the one-compartment model is affected by plasmid costs and intrinsic conjugation rates of the transconjugant, but not by these characteristics of the donor. The parameter space allowing plasmid invasion in the two-compartment model is affected by attachment and detachment rates in the lumen and wall compartment, and by the bacterial density at the wall. The one- and two-compartment models predict the same parameter space for plasmid invasion if the conditions in both compartments are equal to the conditions in the one-compartment model. In contrast, the addition of the wall compartment widens the parameter space allowing invasion compared with the one-compartment model, if the density at the wall is higher than in the lumen, or if the attachment rate at the wall is high and the detachment rate at the wall is low. We also compared the pair-formation models with bulk-conjugation models that describe conjugation by instantaneous transfer of the plasmid at contact between cells, without explicitly describing pair formation. Our results show that pair-formation and bulk-conjugation models predict the same parameter space for plasmid invasion. From our simulations, we conclude that conditions at the gut wall should be taken into account to describe plasmid dynamics in the gut and that transconjugant characteristics rather than donor characteristics should be used to parameterize the models., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Plasmid-mediated antimicrobial resistance is a major contributor to the spread of resistance genes within bacterial communities. Successful plasmid spread depends upon a balance between plasmid fitness effects on the host and rates of horizontal transmission. While these key parameters are readily quantified in vitro , the influence of interactions with other microbiome members is largely unknown. Here, we investigated the influence of three genera of lactic acid bacteria (LAB) derived from the chicken gastrointestinal microbiome on the spread of an epidemic narrow-range ESBL resistance plasmid, IncI1 carrying bla CTX-M-1 , in mixed cultures of isogenic Escherichia coli strains. Secreted products of LAB decreased E. coli growth rates in a genus-specific manner but did not affect plasmid transfer rates. Importantly, we quantified plasmid transfer rates by controlling for density-dependent mating opportunities. Parametrization of a mathematical model with our in vitro estimates illustrated that small fitness costs of plasmid carriage may tip the balance towards plasmid loss under growth conditions in the gastrointestinal tract. This work shows that microbial interactions can influence plasmid success and provides an experimental-theoretical framework for further study of plasmid transfer in a microbiome context.

High throughput sequencing (HTS) has revolutionised virus detection and discovery, allowing for the untargeted characterisation of whole viromes. Viral metagenomics studies have demonstrated the ubiquity of virus infection - often in the absence of disease symptoms - and tend to discover many novel viruses, highlighting the small fraction of virus biodiversity described to date. The majority of the studies using high-throughput sequencing to characterise plant viromes have focused on economically important crops, and only a small number of studies have considered weeds and wild plants. Characterising the viromes of wild plants is highly relevant, as these plants can affect disease dynamics in crops, often by acting as viral reservoirs. Moreover, the viruses in unmanaged systems may also have important effects on wild plant populations and communities. Here, we review metagenomic studies on weeds and wild plants to show the benefits and limitations of this approach and identify knowledge gaps. We consider key genomics developments that are likely to benefit the field in the near future. Although only a small number of HTS studies have been performed on weeds and wild plants, these studies have already discovered many novel viruses, demonstrated unexpected trends in virus distributions, and highlighted the potential of metagenomics as an approach.

Multipartite viruses have segmented genomes and package each of their genome segments individually into distinct virus particles. Multipartitism is common among plant viruses, but why this apparently costly genome organization and packaging has evolved remains unclear. Recently Zhang and colleagues developed network epidemiology models to study the epidemic spread of multipartite viruses and their distribution over plant and animal hosts ( Phys. Rev. Lett. 2019, 123, 138101). In this short commentary, we call into question the relevance of these results because of key model assumptions. First, the model of plant hosts assumes virus transmission only occurs between adjacent plants. This assumption overlooks the basic but imperative fact that most multipartite viruses are transmitted over variable distances by mobile animal vectors, rendering the model results irrelevant to differences between plant and animal hosts. Second, when not all genome segments of a multipartite virus are transmitted to a host, the model assumes an incessant latent infection occurs. This is a bold assumption for which there is no evidence to date, making the relevance of these results to understanding multipartitism questionable., (© The Author(s) 2021. Published by Oxford University Press.)

Background: Conjugation plays a major role in the transmission of plasmids encoding antibiotic resistance genes in both clinical and general settings. The conjugation efficiency is influenced by many biotic and abiotic factors, one of which is the taxonomic relatedness between donor and recipient bacteria. A comprehensive overview of the influence of donor-recipient relatedness on conjugation is still lacking, but such an overview is important to quantitatively assess the risk of plasmid transfer and the effect of interventions which limit the spread of antibiotic resistance, and to obtain parameter values for conjugation in mathematical models. Therefore, we performed a meta-analysis on reported conjugation frequencies from Escherichia coli donors to various recipient species., Results: Thirty-two studies reporting 313 conjugation frequencies for liquid broth matings and 270 conjugation frequencies for filter matings were included in our meta-analysis. The reported conjugation frequencies varied over 11 orders of magnitude. Decreasing taxonomic relatedness between donor and recipient bacteria, when adjusted for confounding factors, was associated with a lower conjugation frequency in liquid matings. The mean conjugation frequency for bacteria of the same order, the same class, and other classes was 10, 20, and 789 times lower than the mean conjugation frequency within the same species, respectively. This association between relatedness and conjugation frequency was not found for filter matings. The conjugation frequency was furthermore found to be influenced by temperature in both types of mating experiments, and in addition by plasmid incompatibility group in liquid matings, and by recipient origin and mating time in filter matings., Conclusions: In our meta-analysis, taxonomic relatedness is limiting conjugation in liquid matings, but not in filter matings, suggesting that taxonomic relatedness is not a limiting factor for conjugation in environments where bacteria are fixed in space.

Multipartite viruses have two or more genome segments, and package different segments into different particle types. Although multipartition is thought to have a cost for virus transmission, its benefits are not clear. Recent experimental work has shown that the equilibrium frequency of viral genome segments, the setpoint genome formula (SGF), can be unbalanced and host-species dependent. These observations have reinvigorated the hypothesis that changes in genome-segment frequencies can lead to changes in virus-gene expression that might be adaptive. Here we explore this hypothesis by developing models of bipartite virus infection, leading to a threefold contribution. First, we show that the SGF depends on the cellular multiplicity of infection (MOI), when the requirements for infection clash with optimizing the SGF for virus-particle yield per cell. Second, we find that convergence on the SGF is very rapid, often occurring within a few cellular rounds of infection. Low and intermediate MOIs lead to faster convergence on the SGF. For low MOIs, this effect occurs because of the requirements for infection, whereas for intermediate MOIs this effect is also due to the high levels of variation generated in the genome formula (GF). Third, we explored the conditions under which a bipartite virus could outcompete a monopartite one. As the heterogeneity between environments and specificity of gene-expression requirements for each environment increased, the bipartite virus was more likely to outcompete the monopartite virus. Under some conditions, changes in the GF helped to exclude the monopartite competitor, highlighting the versatility of the GF. Our results show the inextricable relationship between MOI and the SGF, and suggest that under some conditions, the cost of multipartition can be outweighed by its benefits for the rapid tuning of viral gene expression., (© The Author(s) 2020. Published by Oxford University Press.)

Viruses are widely used as vectors for heterologous gene expression in cultured cells or natural hosts, and therefore a large number of viruses with exogenous sequences inserted into their genomes have been engineered. Many of these engineered viruses are viable and express heterologous proteins at high levels, but the inserted sequences often prove to be unstable over time and are rapidly lost, limiting heterologous protein expression. Although virologists are aware that inserted sequences can be unstable, processes leading to insert instability are rarely considered from an evolutionary perspective. Here, we review experimental work on the stability of inserted sequences over a broad range of viruses, and we present some theoretical considerations concerning insert stability. Different virus genome organizations strongly impact insert stability, and factors such as the position of insertion can have a strong effect. In addition, we argue that insert stability not only depends on the characteristics of a particular genome, but that it will also depend on the host environment and the demography of a virus population. The interplay between all factors affecting stability is complex, which makes it challenging to develop a general model to predict the stability of genomic insertions. We highlight key questions and future directions, finding that insert stability is a surprisingly complex problem and that there is need for mechanism-based, predictive models. Combining theoretical models with experimental tests for stability under varying conditions can lead to improved engineering of viral modified genomes, which is a valuable tool for understanding genome evolution as well as for biotechnological applications, such as gene therapy., (© The Author(s) 2019. Published by Oxford University Press.)

Spodoptera litura is an emerging pest insect in cotton and arable crops in Central Asia. To explore the possibility of using baculoviruses as biological control agents instead of chemical pesticides, in a previous study we characterized a number of S. litura nucleopolyhedrovirus (SpltNPV) isolates from Pakistan. We found significant differences in speed of kill, an important property of a biological control agent. Here we set out to understand the genetic basis of these differences in speed of kill, by comparing the genome of the fast-killing SpltNPV-Pak-TAX1 isolate with that of the slow-killing SpltNPV-Pak-BNG isolate. These two isolates and the SpltNPV-G2 reference strain from China were deep sequenced with Illumina. As expected, the two Pakistani isolates were closely related with >99% sequence identity, whereas the Chinese isolate was more distantly related. We identified two loci that may be associated with the fast action of the SpltNPV-Pak-TAX1 isolate. First, an analysis of rates of synonymous and non-synonymous mutations identified neutral to positive selection on open reading frame (ORF) 122, encoding a viral fibroblast growth factor (vFGF) that is known to affect virulence in other baculoviruses. Second, the homologous repeat region hr17, a putative enhancer of transcription and origin of replication, is absent in SpltNPV-Pak-TAX1 suggesting it may also affect virulence. Additionally, we found there is little genetic variation within both Pakistani isolates, and we identified four genes under positive selection in both isolates that may have played a role in adaptation of SpltNPV to conditions in Central Asia. Our results contribute to the understanding of the enhanced activity of SpltNPV-Pak-TAX1, and may help to select better SpltNPV isolates for the control of S. litura in Pakistan and elsewhere., Competing Interests: The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Multicomponent viral systems face specific challenges when enduring population bottlenecks. These systems can lose coding information due to the lack of co-encapsidation of all the genetic information, at least in a proportion of the capsids in a population. Moreover, bottlenecks can also impact one of the main potential advantages of multicomponent systems: the regulation of gene expression through changes in gene copy frequencies at the population level. How these systems cope with population bottlenecks is far from being clear. Here, two non-exclusive scenarios are described. In the first scenario, population bottlenecks during host infection allow for the isolation of within-host populations with different gene frequencies, leaving the door opened for the selection of populations with adaptive gene frequencies. The second scenario postulates that viruses could influence bottleneck size, at least at certain steps of their life cycle, to limit random changes in gene frequencies. Examples of viral mechanism impacting bottleneck size at cell infection are available and, intriguingly, they can lead to either increases or reductions in bottleneck size. This situation opens the way for putative trade-offs on both gene frequencies and bottleneck sizes that could differ among multicomponent systems., (Copyright © 2018 Elsevier B.V. All rights reserved.)

White spot syndrome virus (WSSV) is an important cause of mortality and economic losses in shrimp farming. Although WSSV-induced mortality is virus dose dependent and WSSV infection does not necessarily lead to mortality, the relationships between virus-particle dose, infection and mortality have not been analysed quantitatively. Here, we explored WSSV dose-response by a combination of experiments, modelling and meta-analysis. We performed dose-response experiments in Penaeus vannamei postlarvae, recorded host mortality and detected WSSV infection. When we fitted infection models to these data, two models-differing in whether they incorporated heterogeneous host susceptibility to the virus or not-were supported for two independent experiments. To determine the generality of these results, we reanalysed published data sets and then performed a meta-analysis. We found that WSSV dose-response kinetics is indeed variable over experiments. We could not clearly identify which specific infection model has the most support by meta-analysis, but we argue that these results also are most concordant with a model incorporating varying levels of heterogeneous host susceptibility to WSSV. We have identified suitable models for analysing WSSV dose-response, which can elucidate the most basic virus-host interactions and help to avoid underestimating WSSV infection at low virus doses., (© 2018 The Authors. Journal of Fish Diseases Published by John Wiley & Sons Ltd.)

Horizontal gene transfer is common among viruses, while they also have highly compact genomes and tend to lose artificial genomic insertions rapidly. Understanding the stability of genomic insertions in viral genomes is therefore relevant for explaining and predicting their evolutionary patterns. Here, we revisit a large body of experimental research on a plant RNA virus, tobacco etch potyvirus (TEV), to identify the patterns underlying the stability of a range of homologous and heterologous insertions in the viral genome. We obtained a wide range of estimates for the recombination rate-the rate at which deletions removing the insertion occur-and these appeared to be independent of the type of insertion and its location. Of the factors we considered, recombination rate was the best predictor of insertion stability, although we could not identify the specific sequence characteristics that would help predict insertion instability. We also considered experimentally the possibility that functional insertions lead to higher mutational robustness through increased redundancy. However, our observations suggest that both functional and non-functional increases in genome size decreased the mutational robustness. Our results therefore demonstrate the importance of recombination rates for predicting the long-term stability and evolution of viral RNA genomes and suggest that there are unexpected drawbacks to increases in genome size for mutational robustness.

While synonymous mutations were long thought to be without phenotypic consequences, there is growing evidence they can affect gene expression, protein folding, and ultimately the fitness of an organism. In only a few cases have the mechanisms by which synonymous mutations affect the phenotype been elucidated. We previously identified 48 mutations in TEM-1 β-lactamase that increased resistance of Escherichia coli to cefotaxime, 10 of which were synonymous. To better understand the molecular mechanisms underlying the beneficial effect of these synonymous mutations, we made a series of measurements for a panel containing the 10 synonymous together with 10 non-synonymous mutations as a reference. Whereas messenger levels were unaffected, we found that total and functional TEM protein levels were higher for 5 out of 10 synonymous mutations. These observations suggest that some of these mutations act on translation or a downstream process. Similar effects were observed for some small-benefit non-synonymous mutations, suggesting a similar causal mechanism. For the synonymous mutations, we found that the cost of resistance scales with TEM protein levels. A resistance landscape for four synonymous mutations revealed strong epistasis: none of the combinations of mutations exceeded the resistance of the largest-effect mutation and there were synthetically neutral combinations. By considering combined effects of these mutations, we could infer that functional TEM protein level is a multi-dimensional phenotype. These results suggest that synonymous mutations may have beneficial effects by increasing the expression of an enzyme with low substrate activity, which may be realized via multiple, yet unknown, post-transcriptional mechanisms.

Populations with large mutation supplies adapt via the "greedy" substitution of the fittest genotype available, leading to fast and repeatable short-term responses. At longer time scales, smaller mutation supplies may in theory lead to larger improvements when distant high-fitness genotypes more readily evolve from lower-fitness intermediates. Here we test for long-term adaptive benefits from small mutation supplies using in vitro evolution of an antibiotic-degrading enzyme in the presence of a novel antibiotic. Consistent with predictions, large mutant libraries cause rapid initial adaptation via the substitution of cohorts of mutations, but show later deceleration and convergence. Smaller libraries show on average smaller initial, but also more variable, improvements, with two lines yielding alleles with exceptionally high resistance levels. These two alleles share three mutations with the large-library alleles, which are known from previous work, but also have unique mutations. Replay evolution experiments and analyses of the adaptive landscape of the enzyme suggest that the benefit resulted from a combination of avoiding mutational cohorts leading to local peaks and chance. Our results demonstrate adaptive benefits from limited mutation supplies on a rugged fitness landscape, which has implications for artificial selection protocols in biotechnology and argues for a better understanding of mutation supplies in clinical settings., Competing Interests: The authors declare no conflict of interest.

Whether evolution can be predicted is a key question in evolutionary biology. Here we set out to better understand the repeatability of evolution, which is a necessary condition for predictability. We explored experimentally the effect of mutation supply and the strength of selective pressure on the repeatability of selection from standing genetic variation. Different sizes of mutant libraries of antibiotic resistance gene TEM-1 β-lactamase in Escherichia coli, generated by error-prone PCR, were subjected to different antibiotic concentrations. We determined whether populations went extinct or survived, and sequenced the TEM gene of the surviving populations. The distribution of mutations per allele in our mutant libraries followed a Poisson distribution. Extinction patterns could be explained by a simple stochastic model that assumed the sampling of beneficial mutations was key for survival. In most surviving populations, alleles containing at least one known large-effect beneficial mutation were present. These genotype data also support a model which only invokes sampling effects to describe the occurrence of alleles containing large-effect driver mutations. Hence, evolution is largely predictable given cursory knowledge of mutational fitness effects, the mutation rate and population size. There were no clear trends in the repeatability of selected mutants when we considered all mutations present. However, when only known large-effect mutations were considered, the outcome of selection is less repeatable for large libraries, in contrast to expectations. We show experimentally that alleles carrying multiple mutations selected from large libraries confer higher resistance levels relative to alleles with only a known large-effect mutation, suggesting that the scarcity of high-resistance alleles carrying multiple mutations may contribute to the decrease in repeatability at large library sizes.

The existence of multipartite viruses is an intriguing mystery in evolutionary virology. Several hypotheses suggest benefits that should outweigh the costs of a reduced transmission efficiency and of segregation of coadapted genes associated with encapsidating each segment into a different particle. Advantages range from increasing genome size despite high mutation rates, faster replication, more efficient selection resulting from reassortment during mixed infections, better regulation of gene expression, or enhanced virion stability and cell-to-cell movement. However, support for these hypotheses is scarce. Here we report experiments testing whether an evolutionary stable equilibrium exists for the three genomic RNAs of Alfalfa mosaic virus (AMV). Starting infections with different segment combinations, we found that the relative abundance of each segment evolves towards a constant ratio. Population genetic analyses show that the segment ratio at this equilibrium is determined by frequency-dependent selection. Replication of RNAs 1 and 2 was coupled and collaborative, whereas the replication of RNA 3 interfered with the replication of the other two. We found that the equilibrium solution is slightly different for the total amounts of RNA produced and encapsidated, suggesting that competition exists between all RNAs during encapsidation. Finally, we found that the observed equilibrium appears to be host-species dependent.

Horizontal gene transfer (HGT) is pervasive in viruses and thought to be a key mechanism in their evolution. On the other hand, strong selective constraints against increasing genome size are an impediment for HGT, rapidly purging horizontally transferred sequences and thereby potentially hindering evolutionary innovation. Here, we explore experimentally the evolutionary fate of viruses with simulated HGT events, using the plant RNA virus Tobacco etch virus (TEV), by separately introducing two functional, exogenous sequences to its genome. One of the events simulates the acquisition of a new function though HGT of a conserved AlkB domain, responsible for the repair of alkylation or methylation damage in many organisms. The other event simulates the acquisition of a sequence that duplicates an existing function, through HGT of the 2b RNA silencing suppressor from Cucumber mosaic virus. We then evolved these two viruses, tracked the maintenance of the horizontally transferred sequences over time, and for the final virus populations, sequenced their genome and measured viral fitness. We found that the AlkB domain was rapidly purged from the TEV genome, restoring fitness to wild-type levels. Conversely, the 2b gene was stably maintained and did not have a major impact on viral fitness. Moreover, we found that 2b is functional in TEV, as it provides a replicative advantage when the RNA silencing suppression domain of HC-Pro is mutated. These observations suggest a potentially interesting role for HGT of short functional sequences in ameliorating evolutionary constraints on viruses, through the duplication of functions., (© The Author(s) 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Theory suggests that high virulence could hinder between-host transmission of microparasites, and that virulence therefore will evolve to lower levels. Alternatively, highly virulent microparasites could also curtail host development, thereby limiting both the host resources available to them and their own within-host effective population size. In this case, high virulence might restrain the mutation supply rate and increase the strength with which genetic drift acts on microparasite populations. Thereby, this alternative explanation limits the microparasites' potential to adapt to the host and ultimately the ability to evolve lower virulence. As a first exploration of this hypothesis, we evolved Tobacco etch virus carrying an eGFP fluorescent marker in two semi-permissive host species, Nicotiana benthamiana and Datura stramonium, for which it has a large difference in virulence. We compared the results to those previously obtained in the natural host, Nicotiana tabacum, where we have shown that carriage of eGFP has a high fitness cost and its loss serves as a real-time indicator of adaptation., Results: After over half a year of evolution, we sequenced the genomes of the evolved lineages and measured their fitness. During the evolution experiment, marker loss leading to viable virus variants was only observed in one lineage of the host for which the virus has low virulence, D. stramonium. This result was consistent with the observation that there was a fitness cost of eGFP in this host, while surprisingly no fitness cost was observed in the host for which the virus has high virulence, N. benthamiana. Furthermore, in both hosts we observed increases in viral fitness in few lineages, and host-specific convergent evolution at the genomic level was only found in N. benthamiana., Conclusions: The results of this study do not lend support to the hypothesis that high virulence impedes microparasites' evolution. Rather, they exemplify that jumps between host species can be game changers for evolutionary dynamics. When considering the evolution of genome architecture, host species jumps might play a very important role, by allowing evolutionary intermediates to be competitive.

One of the striking features of many eukaryotes is the apparent amount of redundancy in coding and non-coding elements of their genomes. Despite the possible evolutionary advantages, there are fewer examples of redundant sequences in viral genomes, particularly those with RNA genomes. The factors constraining the maintenance of redundant sequences in present-day RNA virus genomes are not well known. Here, we use Tobacco etch virus, a plant RNA virus, to investigate the stability of genetically redundant sequences by generating viruses with potentially beneficial gene duplications. Subsequently, we tested the viability of these viruses and performed experimental evolution. We found that all gene duplication events resulted in a loss of viability or in a significant reduction in viral fitness. Moreover, upon analyzing the genomes of the evolved viruses, we always observed the deletion of the duplicated gene copy and maintenance of the ancestral copy. Interestingly, there were clear differences in the deletion dynamics of the duplicated gene associated with the passage duration and the size and position of the duplicated copy. Based on the experimental data, we developed a mathematical model to characterize the stability of genetically redundant sequences, and showed that fitness effects are not enough to predict genomic stability. A context-dependent recombination rate is also required, with the context being the duplicated gene and its position. Our results therefore demonstrate experimentally the deleterious nature of gene duplications in RNA viruses. Beside previously described constraints on genome size, we identified additional factors that reduce the likelihood of the maintenance of duplicated genes., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The order in which genes are organized within a genome is generally not conserved between distantly related species. However, within virus orders and families, strong conservation of gene order is observed. The factors that constrain or promote gene-order diversity are largely unknown, although the regulation of gene expression is one important constraint for viruses. Here we investigate why gene order is conserved for a positive-strand RNA virus encoding a single polyprotein in the context of its authentic multicellular host. Initially, we identified the most plausible trajectory by which alternative gene orders could evolve. Subsequently, we studied the accessibility of key steps along this evolutionary trajectory by constructing two virus intermediates: (1) duplication of a gene followed by (2) loss of the ancestral gene. We identified five barriers to the evolution of alternative gene orders. First, the number of viable positions for reordering is limited. Second, the within-host fitness of viruses with gene duplications is low compared to the wild-type virus. Third, after duplication, the ancestral gene copy is always maintained and never the duplicated one. Fourth, viruses with an alternative gene order have even lower fitness than viruses with gene duplications. Fifth, after more than half a year of evolution in isolation, viruses with an alternative gene order are still vastly inferior to the wild-type virus. Our results show that all steps along plausible evolutionary trajectories to alternative gene orders are highly unlikely. Hence, the inaccessibility of these trajectories probably contributes to the conservation of gene order in present-day viruses., (Copyright © 2016 by the Genetics Society of America.)

For virus infections of multicellular hosts, narrow genetic bottlenecks during transmission and within-host spread appear to be widespread. These bottlenecks will affect the maintenance of genetic variation in a virus population and the prevalence of mixed-strain infections, thereby ultimately determining the strength with which different random forces act during evolution. Here we consider different approaches for estimating bottleneck sizes and weigh their merits. We then review quantitative estimates of bottleneck size during cellular infection, within-host spread, horizontal transmission, and finally vertical transmission. In most cases we find that bottlenecks do regularly occur, although in many cases they appear to be virion-concentration dependent. Finally, we consider the evolutionary implications of genetic bottlenecks during virus infection. Although on average strong bottlenecks will lead to declines in fitness, we consider a number of scenarios in which bottlenecks could also be advantageous for viruses.

Populations of plant RNA viruses are highly polymorphic in infected plants, which may allow rapid within-host evolution. To understand tobacco etch potyvirus (TEV) evolution, longitudinal samples from experimentally evolved populations in the natural host tobacco and from the alternative host pepper were phenotypically characterized and genetically analyzed. Temporal and compartmental variabilities of TEV populations were quantified using high throughput Illumina sequencing and population genetic approaches. Of the two viral phenotypic traits measured, virulence increased in the novel host but decreased in the original one, and viral load decreased in both hosts, though to a lesser extent in the novel one. Dynamics of population genetic diversity were also markedly different among hosts. Population heterozygosity increased in the ancestral host, with a dominance of synonymous mutations fixed, whereas it did not change or even decreased in the new host, with an excess of nonsynonymous mutations. All together, these observations suggest that directional selection is the dominant evolutionary force in TEV populations evolving in a novel host whereas either diversifying selection or random genetic drift may play a fundamental role in the natural host. To better understand these evolutionary dynamics, we developed a computer simulation model that incorporates the effects of mutation, selection, and drift. Upon parameterization with empirical data from previous studies, model predictions matched the observed patterns, thus reinforcing our idea that the empirical patterns of mutation accumulation represent adaptive evolution., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The independent action hypothesis is a simple model of pathogen infection that can make many useful predictions on infection kinetics and, therefore, a number of different tests of independent action have been developed. However, some of these analyses are rather sophisticated, limiting their appeal to experimentalists, and it is also unclear how well the different tests perform. Here, we developed and evaluated a simple and robust new test of independent action. Our new test is based on using a constant inoculum dose of one pathogen variant, varying the dose of a second variant, and then quantifying the infection response for the first variant. We simulated infection data in which we introduced deviations from independent action, experimental variation, or both. Simulations showed that our new procedure has many advantages over the existing tests of independent action, especially if only systemic-infection data are available. We also performed experimental tests of our new procedure using two marked Tobacco etch virus (TEV) variants. We found minor deviations from the independent action model, which were not detected by previous tests using existing methods, exemplifying the utility of this approach. We discuss the implications for TEV infection kinetics and consider how to reconcile different dose-dependent effects.

Unlabelled: Superinfection exclusion (SIE), the ability of an established virus infection to interfere with a secondary infection by the same or a closely related virus, has been described for different viruses, including important pathogens of humans, animals, and plants. Citrus tristeza virus (CTV), a positive-sense RNA virus, represents a valuable model system for studying SIE due to the existence of several phylogenetically distinct strains. Furthermore, CTV allows SIE to be examined at the whole-organism level. Previously, we demonstrated that SIE by CTV is a virus-controlled function that requires the viral protein p33. In this study, we show that p33 mediates SIE at the whole-organism level, while it is not required for exclusion at the cellular level. Primary infection of a host with a fluorescent protein-tagged CTV variant lacking p33 did not interfere with the establishment of a secondary infection by the same virus labeled with a different fluorescent protein. However, cellular coinfection by both viruses was rare. The obtained observations, along with estimates of the cellular multiplicity of infection (MOI) and MOI model selection, suggested that low levels of cellular coinfection appear to be best explained by exclusion at the cellular level. Based on these results, we propose that SIE by CTV is operated at two levels--the cellular and the whole-organism levels--by two distinct mechanisms that could function independently. This novel aspect of viral SIE highlights the intriguing complexity of this phenomenon, further understanding of which may open up new avenues to manage virus diseases., Importance: Many viruses exhibit superinfection exclusion (SIE), the ability of an established virus infection to interfere with a secondary infection by related viruses. SIE plays an important role in the pathogenesis and evolution of virus populations. The observations described here suggest that SIE could be controlled independently at different levels of the host: the whole-organism level or the level of individual cells. The p33 protein of citrus tristeza virus (CTV), an RNA virus, was shown to mediate SIE at the whole-organism level, while it appeared not to be required for exclusion at the cellular level. SIE by CTV is, therefore, highly complex and appears to use mechanisms different from those proposed for other viruses. A better understanding of this phenomenon may lead to the development of new strategies for controlling viral diseases in human populations and agroecosystems., (Copyright © 2014, American Society for Microbiology. All Rights Reserved.)

The cornerstone of today's plant virology consists of deciphering the molecular and mechanistic basis of host-pathogen interactions. Among these interactions, the onset of systemic infection is a fundamental variable in studying both within- and between-host infection dynamics, with implications in epidemiology. Here, we developed a mechanistic model using probabilistic and spatio-temporal concepts to explain dynamic signatures of virus systemic infection. The model dealt with the inherent characteristic of plant viruses to use two different and sequential stages for their within-host propagation: cell-to-cell movement from the initial infected cell and systemic spread by reaching the vascular system. We identified the speed of cell-to-cell movement and the number of primary infection foci in the inoculated leaf as the key factors governing this dynamic process. Our results allowed us to quantitatively understand the timing of the onset of systemic infection, describing this global process as a consequence of local spread of viral populations. Finally, we considered the significance of our predictions for the evolution of plant RNA viruses., (© 2014 The Author(s) Published by the Royal Society. All rights reserved.)

Potyviruses express most of their proteins from a long open reading frame that is translated into a large polyprotein processed by three viral proteases. To understand the constraints on potyvirus genome organization, we relocated the viral RNA-dependent RNA polymerase (NIb) cistron to all possible intercistronic positions of the Tobacco etch virus (TEV) polyprotein. Only viruses with NIb at the amino terminus of the polyprotein or in between P1 and HC-Pro were viable in tobacco plants.

Phylogenomic evidence suggested that recombination is an important evolutionary force for potyviruses, one of the larger families of plant RNA viruses. However, mixed-genotype potyvirus infections are marked by low levels of cellular coinfection, precluding template switching and recombination events between virus genotypes during genomic RNA replication. To reconcile these conflicting observations, we evaluated the in vivo recombination rate (rg) of Tobacco etch virus (TEV; genus Potyvirus, family Potyviridae) by coinfecting plants with pairs of genotypes marked with engineered restriction sites as neutral markers. The recombination rate was then estimated using two different approaches: (i) a classical approach that assumed recombination between marked genotypes can occur in the whole virus population, rendering an estimate of rg = 7.762 × 10(-8) recombination events per nucleotide site per generation, and (ii) an alternative method that assumed recombination between marked genotypes can occur only in coinfected cells, rendering a much higher estimate of rg = 3.427 × 10(-5) recombination events per nucleotide site per generation. This last estimate is similar to the TEV mutation rate, suggesting that recombination should be at least as important as point mutation in creating variability. Finally, we compared our mutation and recombination rate estimates to those reported for animal RNA viruses. Our analysis suggested that high recombination rates may be an unavoidable consequence of selection for fast replication at the cost of low fidelity.

Nonretroviral integrated RNA viruses (NIRVs) are genes of nonretroviral RNA viruses found in the genomes of many eukaryotic organisms. NIRVs are thought to sometimes confer virus resistance, meaning that they could impact spread of the virus in the host population. However, a NIRV that is expressed may also impact the evolution of virus populations within host organisms. Here, we experimentally addressed the evolution of a virus in a host expressing a NIRV using Tobacco etch virus (TEV), a plant RNA virus, and transgenic tobacco plants expressing its replicase, NIb. We found that a virus missing the NIb gene, TEV-ΔNIb, which is incapable of autonomous replication in wild-type plants, had a higher fitness than the full-length TEV in the transgenic plants. Moreover, when the full-length TEV was evolved by serial passages in transgenic plants, we observed genomic deletions within NIb--and in some cases the adjacent cistrons--starting from the first passage. When we passaged TEV and TEV-ΔNIb in transgenic plants, we found mutations in proteolytic sites, but these only occurred in TEV-ΔNIb lineages, suggesting the adaptation of polyprotein processing to altered NIb expression. These results raise the possibility that NIRV expression can indeed induce the deletion of the corresponding genes in the viral genome, resulting in the formation of viruses that are replication defective in hosts that do not express the same NIRV. Moreover, virus genome evolution was contingent upon the deletion of the viral replicase, suggesting NIRV expression could also alter patterns of virus evolution.

A multicellular organism is not a monolayer of cells in a flask; it is a complex, spatially structured environment, offering both challenges and opportunities for viruses to thrive. Whereas virus infection dynamics at the host and within-cell levels have been documented, the intermediate between-cell level remains poorly understood. Here, we used flow cytometry to measure the infection status of thousands of individual cells in virus-infected plants. This approach allowed us to determine accurately the number of cells infected by two virus variants in the same host, over space and time as the virus colonizes the host. We found a low overall frequency of cellular infection (<0.3), and few cells were coinfected by both virus variants (<0.1). We then estimated the cellular contagion rate (R), the number of secondary infections per infected cell per day. R ranged from 2.43 to values not significantly different from zero, and generally decreased over time. Estimates of the cellular multiplicity of infection (MOI), the number of virions infecting a cell, were low (<1.5). Variance of virus-genotype frequencies increased strongly from leaf to cell levels, in agreement with a low MOI. Finally, there were leaf-dependent differences in the ease with which a leaf could be colonized, and the number of virions effectively colonizing a leaf. The modeling of infection patterns suggests that the aggregation of virus-infected cells plays a key role in limiting spread; matching the observation that cell-to-cell movement of plant viruses can result in patches of infection. Our results show that virus expansion at the between-cell level is restricted, probably due to the host environment and virus infection itself.

Viruses have evolved highly streamlined genomes and a variety of mechanisms to compress them, suggesting that genome size is under strong selection. Horizontal gene transfer has, on the other hand, played an important role in virus evolution. However, evolution cannot integrate initially nonfunctional sequences into the viral genome if they are rapidly purged by selection. Here we report on the experimental evolution of pseudogenization in virus genomes using a plant RNA virus expressing a heterologous gene. When long 9-week passages were performed, the added gene was lost in all lineages, whereas viruses with large genomic deletions were fixed in only two out of ten 3-week lineages and none in 1-week lineages. Illumina next-generation sequencing revealed considerable convergent evolution in the 9- and 3-week lineages with genomic deletions. Genome size was correlated to within-host competitive fitness, although there was no correlation with virus accumulation or virulence. Within-host competitive fitness of the 3-week virus lineages without genomic deletions was higher than for the 1-week lineages. Our results show that the strength of selection for a reduced genome size and the rate of pseudogenization depend on demographic conditions. Moreover, for the 3-week passage condition, we observed increases in within-host fitness, whereas selection was not strong enough to quickly remove the nonfunctional heterologous gene. These results suggest a demographically determined "sweet spot" might exist, where heterologous insertions are not immediately lost while evolution can act to integrate them into the viral genome.

Multipartite plant viruses were discovered because of discrepancies between the observed dose response and predictions of the independent-action hypothesis (IAH) model. Theory suggests that the number of genome segments predicts the shape of the dose-response curve, but a rigorous test of this hypothesis has not been reported. Here, Alfalfa mosaic virus (AMV), a tripartite Alfamovirus, and transgenic Nicotianatabacum plants expressing no (wild type), one (P2), or two (P12) viral genome segments were used to test whether the number of genome segments necessary for infection predicts the dose response. The dose-response curve of wild-type plants was steep and congruent with the predicted kinetics of a multipartite virus, confirming previous results. Moreover, for P12 plants, the data support the IAH model, showing that the expression of virus genome segments by the host plant can modulate the infection kinetics of a tripartite virus to those of a monopartite virus. However, the different types of virus particles occurred at different frequencies, with a ratio of 116:45:1 (RNA1 to RNA2 to RNA3), which will affect infection kinetics and required analysis with a more comprehensive infection model. This analysis showed that each type of virus particle has a different probability of invading the host plant, at both the primary- and systemic-infection levels. While the number of genome segments affects the dose response, taking into consideration differences in the infection kinetics of the three types of AMV particles results in a better understanding of the infection process.

Antirrhinum majus Rosea1 (Ros1) is an MYB-related transcription factor that induces anthocyanin biosynthesis in plant tissues, and has been shown to be suitable for visual tracking of virus infection in plants. However, activation of anthocyanin biosynthesis has far reaching effects on plant physiology and could consequently have negative effects on viral replication. Therefore, viruses carrying the Ros1 marker might have a low fitness and consequently rapidly lose the marker. To compare the stability of the Ros1 marker, we generated Tobacco etch virus (TEV) based constructs containing either Ros1 or the enhanced green fluorescent protein (eGFP) between the NIb and CP cistrons (TEV-Ros1 and TEV-eGFP, respectively). We measured the within-host competitive fitness of both viruses by direct competitions with a common competitor during infection of Nicotiana tabacum. The fitness of TEV-Ros1 was significantly lower than that of TEV-eGFP, and both recombinant viruses had a significantly lower fitness than the wild-type virus. Nevertheless, after seven weeks of infection in N. tabacum, similar levels of marker gene instability where found for both viruses. Despite lower fitness of the marked virus, Ros1 is therefore a viable alternative marker for tracking viral infection in plants.

The cellular multiplicity of infection (MOI) is a key parameter for describing the interactions between virions and cells, predicting the dynamics of mixed-genotype infections, and understanding virus evolution. Two recent studies have reported in vivo MOI estimates for Tobacco mosaic virus (TMV) and Cauliflower mosaic virus (CaMV), using sophisticated approaches to measure the distribution of two virus variants over host cells. Although the experimental approaches were similar, the studies employed different definitions of MOI and estimation methods. Here, new model-selection-based methods for calculating MOI were developed. Seven alternative models for predicting MOI were formulated that incorporate an increasing number of parameters. For both datasets the best-supported model included spatial segregation of virus variants over time, and to a lesser extent aggregation of virus-infected cells was also implicated. Three methods for MOI estimation were then compared: the two previously reported methods and the best-supported model. For CaMV data, all three methods gave comparable results. For TMV data, the previously reported methods both predicted low MOI values (range: 1.04-1.23) over time, whereas the best-supported model predicted a wider range of MOI values (range: 1.01-2.10) and an increase in MOI over time. Model selection can therefore identify suitable alternative MOI models and suggest key mechanisms affecting the frequency of coinfected cells. For the TMV data, this leads to appreciable differences in estimated MOI values.

Defective interfering (DI) viruses are thought to cause oscillations in virus levels, known as the 'Von Magnus effect'. Interference by DI viruses has been proposed to underlie these dynamics, although experimental tests of this idea have not been forthcoming. For the baculoviruses, insect viruses commonly used for the expression of heterologous proteins in insect cells, the molecular mechanisms underlying DI generation have been investigated. However, the dynamics of baculovirus populations harboring DIs have not been studied in detail. In order to address this issue, we used quantitative real-time PCR to determine the levels of helper and DI viruses during 50 serial passages of Autographa californica multiple nucleopolyhedrovirus (AcMNPV) in Sf21 cells. Unexpectedly, the helper and DI viruses changed levels largely in phase, and oscillations were highly irregular, suggesting the presence of chaos. We therefore developed a simple mathematical model of baculovirus-DI dynamics. This theoretical model reproduced patterns qualitatively similar to the experimental data. Although we cannot exclude that experimental variation (noise) plays an important role in generating the observed patterns, the presence of chaos in the model dynamics was confirmed with the computation of the maximal Lyapunov exponent, and a Ruelle-Takens-Newhouse route to chaos was identified at decreasing production of DI viruses, using mutation as a control parameter. Our results contribute to a better understanding of the dynamics of DI baculoviruses, and suggest that changes in virus levels over passages may exhibit chaos.

In the clinic, farm, or field, for many viruses there is a high prevalence of mixed-genotype infections, indicating that multiple virions have initiated infection and that there can be multiple sites of primary infection within the same host. The dynamic process by which multiple primary infection sites interact with each other and the host is poorly understood, undoubtedly due to its high complexity. In this study, we attempted to unravel the basic interactions underlying this process using a plant RNA virus, as removing the inoculated leaf can instantly and rigorously eliminate all primary infection sites. Effective population size in the inoculated leaf and time of removal of the inoculated leaf were varied in experiments, and it was found that both factors positively influenced if the plant became systemically infected and what proportion of cells in the systemic tissue were infected, as measured by flow cytometry. Fitting of probabilistic models of infection to our data demonstrated that a null model in which the action of each focus is independent of the presence of other foci was better supported than a dependent-action model. The cumulative effect of independently acting foci therefore determined when plants became infected and how many individual cells were infected. There was no evidence for interference between primary infection sites, which is surprising given the planar structure of leaves. By showing that a simple null model is supported, we experimentally confirmed--to our knowledge for the first time--the minimal components that dictate interactions of a conspecific virus population establishing systemic infection.

White spot syndrome virus (WSSV) is the most important pathogen in shrimp farming systems worldwide including the Mekong Delta, Vietnam. The genome of WSSV is characterized by the presence of two major 'indel regions' found at ORF14/15 and ORF23/24 (WSSV-Thailand) and three regions with variable number tandem repeats (VNTR) located in ORF75, ORF94 and ORF125. In the current study, we investigated whether or not the number of repeat units in the VNTRs correlates with virus outbreak status and/or shrimp farming practice. We analysed 662 WSSV samples from individual WSSV-infected Penaeus monodon shrimp from 104 ponds collected from two important shrimp farming regions of the Mekong Delta: Ca Mau and Bac Lieu. Using this large data set and statistical analysis, we found that for ORF94 and ORF125, the mean number of repeat units (RUs) in VNTRs was significantly lower in disease outbreak ponds than in non-outbreak ponds. Although a higher mean RU number was observed in the improved-extensive system than in the rice-shrimp or semi-intensive systems, these differences were not significant. VNTR sequences are thus not only useful markers for studying WSSV genotypes and populations, but specific VNTR variants also correlate with disease outbreaks in shrimp farming systems., (© 2012 Blackwell Publishing Ltd.)

Effective population size (N(e)) is a key parameter for understanding evolutionary processes, but it is generally not considered in epidemiological studies or in studying infections of individual hosts. Whether N(e) has an effect on the onset of symptoms and viral accumulation in Tobacco etch virus (TEV) infection of Nicotiana tabacum plants is considered here. Using mixtures of TEV variants carrying fluorescent markers, the dose dependence of N(e) was confirmed, and the inoculation procedure was found to be the main source of variation in these experiments. Whereas the onset of symptoms was independent of N(e), there was less and more variable accumulation at 6 days postinoculation for small N(e) values (N(e) < 5). The observed variation in accumulation was not heritable, however, suggesting that this variation was not due to the fixation of deleterious mutations in the small founder populations. On the other hand, virus-induced fluorescence and accumulation in the inoculated leaf were strongly N(e) dependent. Systemic accumulation was independent of N(e), although removal of the inoculated leaf led to a small reduction in systemic accumulation for small N(e) values. For whole plants, N(e)-dependent effects on accumulation were no longer observed at 9 days postinoculation. Therefore, the effects of N(e) on accumulation are due mainly to limited expansion in the inoculated leaf and are transient. In this system, N(e)-dependent effects will be strongest at low doses and early in infection. We conclude that N(e) can have implications for epidemiology and infection at the individual host level, beyond determining the rate of mixed-genotype infection.

Effective population size (N(e)) determines the strength of genetic drift and the frequency of co-infection by multiple genotypes, making it a key factor in viral evolution. Experimental estimates of N(e) for different plant viruses have, however, rendered diverging results. The independent action hypothesis (IAH) states that each virion has a probability of infection, and that virions act independent of one another during the infection process. A corollary of IAH is that N(e) must be dose dependent. A test of IAH for a plant virus has not been reported yet. Here we perform a test of an IAH infection model using a plant RNA virus, Tobacco etch virus (TEV) variants carrying GFP or mCherry fluorescent markers, in Nicotiana tabacum and Capsicum annuum plants. The number of primary infection foci increased linearly with dose, and was similar to a Poisson distribution. At high doses, primary infection foci containing both genotypes were found at a low frequency (<2%). The probability that a genotype that infected the inoculated leaf would systemically infect that plant was near 1, although in a few rare cases genotypes could be trapped in the inoculated leaf by being physically surrounded by the other genotype. The frequency of mixed-genotype infection could be predicted from the mean number of primary infection foci using the independent-action model. Independent action appears to hold for TEV, and N(e) is therefore dose-dependent for this plant RNA virus. The mean number of virions causing systemic infection can be very small, and approaches 1 at low doses. Dose-dependency in TEV suggests that comparison of N(e) estimates for different viruses are not very meaningful unless dose effects are taken into consideration.