Within-host spatiotemporal dynamics of plant virus infection at the cellular level

On the ever-growing functional versatility of the CRISPR-Cas13 system

CRISPR-Cas systems evolved in prokaryotes to implement a powerful antiviral immune response as a result of sequence-specific targeting by ribonucleoproteins. One of such systems consists of an RNA-guided RNA endonuclease, known as CRISPR-Cas13. In very recent years, this system is being repurposed in different ways in order to decipher and engineer gene expression programmes. Here, we discuss the functional versatility of the CRISPR-Cas13 system, which includes the ability for RNA silencing, RNA editing, RNA tracking, nucleic acid detection and translation regulation. This functional palette makes the CRISPR-Cas13 system a relevant tool in the broad field of systems and synthetic biology.

Plant Virus Adaptation to New Hosts: A Multi-scale Approach

Viruses are studied at each level of biological complexity: from within-cells to ecosystems. The same basic evolutionary forces and principles operate at each level: mutation and recombination, selection, genetic drift, migration, and adaptive trade-offs. Great efforts have been put into understanding each level in great detail, hoping to predict the dynamics of viral population, prevent virus emergence, and manage their spread and virulence. Unfortunately, we are still far from this. To achieve these ambitious goals, we advocate for an integrative perspective of virus evolution. Focusing in plant viruses, we illustrate the pervasiveness of the above-mentioned principles. Beginning at the within-cell level, we describe replication modes, infection bottlenecks, and cellular contagion rates. Next, we move up to the colonization of distal tissues, discussing the fundamental role of random events. Then, we jump beyond the individual host and discuss the link between transmission mode and virulence. Finally, at the community level, we discuss properties of virus-plant infection networks. To close this review we propose the multilayer network theory, in which elements at different layers are connected and submit to their own dynamics that feed across layers, resulting in new emerging properties, as a way to integrate information from the different levels.

Natural Selection, Intracellular Bottlenecks of Virus Populations, and Viral Superinfection Exclusion

Natural selection acts on cellular organisms by ensuring the genes responsible for an advantageous phenotype consistently reap the phenotypic advantage. This is possible because reproductive cells of these organisms are almost always haploid, separating the beneficial gene from its rival allele at every generation. How natural selection acts on plus-strand RNA viruses is unclear because these viruses frequently load host cells with numerous genome copies and replicate thousands of progeny genomes in each cell. Recent studies suggest that these viruses encode the Bottleneck, Isolate, Amplify, Select (BIAS) mechanism that blocks all but a few viral genome copies from replication, thus creating the environment in which the bottleneck-escaping viral genome copies are isolated from each other, allowing natural selection to reward beneficial mutations and purge lethal errors. This BIAS mechanism also blocks the genomes of highly homologous superinfecting viruses, thus explaining cellular-level superinfection exclusion. Expected final online publication date for the Annual Review of Virology, Volume 9 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Mathematical Modeling Suggests Cooperation of Plant-Infecting Viruses

Viruses are major pathogens of agricultural crops. Viral infections often start after the virus enters the outer layer of a tissue, and many successful viruses, after local replication in the infected tissue, are able to spread systemically. Quantitative details of virus dynamics in plants, however, are poorly understood, in part, because of the lack of experimental methods which allow the accurate measurement of the degree of infection in individual plant tissues. Recently, a group of researchers followed the kinetics of infection of individual cells in leaves of Nicotiana tabacum plants using Tobacco etch virus (TEV) expressing either Venus or blue fluorescent protein (BFP). Assuming that viral spread occurs from lower to upper leaves, the authors fitted a simple mathematical model to the frequency of cellular infection by the two viral variants found using flow cytometry. While the original model could accurately describe the kinetics of viral spread locally and systemically, we found that many alternative versions of the model, for example, if viral spread starts at upper leaves and progresses to lower leaves or when virus dissemination is stopped due to an immune response, fit the data with reasonable quality, and yet with different parameter estimates. These results strongly suggest that experimental measurements of the virus infection in individual leaves may not be sufficient to identify the pathways of viral dissemination between different leaves and reasons for viral control. We propose experiments that may allow discrimination between the alternatives. By analyzing the kinetics of coinfection of individual cells by Venus and BFP strains of TEV we found a strong deviation from the random infection model, suggesting cooperation between the two strains when infecting plant cells. Importantly, we showed that many mathematical models on the kinetics of coinfection of cells with two strains could not adequately describe the data, and the best fit model needed to assume (i) different susceptibility of uninfected cells to infection by two viruses locally in the leaf vs. systemically from other leaves, and (ii) decrease in the infection rate depending on the fraction of uninfected cells which could be due to a systemic immune response. Our results thus demonstrate the difficulty in reaching definite conclusions from extensive and yet limited experimental data and provide evidence of potential cooperation between different viral variants infecting individual cells in plants.

Reprogramming of RNA silencing triggered by cucumber mosaic virus infection in Arabidopsis

BACKGROUND

RNA silencing has an important role mediating sequence-specific virus resistance in plants. The complex interaction of viruses with RNA silencing involves the loading of viral small interfering RNAs (vsiRNAs) into its host ARGONAUTE (AGO) proteins. As a side effect of their antiviral activity, vsiRNAs loading into AGO proteins can also mediate the silencing of endogenous genes. Here, we analyze at the genome-wide level both aspects of the interference of cucumber mosaic virus (CMV) with the RNA silencing machinery of Arabidopsis thaliana.

RESULTS

We observe CMV-derived vsiRNAs affect the levels of endogenous sRNA classes. Furthermore, we analyze the incorporation of vsiRNAs into AGO proteins with a described antiviral role and the viral suppressor of RNA silencing (VSR) 2b, by combining protein immunoprecipitation with sRNA high-throughput sequencing. Interestingly, vsiRNAs represent a substantial percentage of AGO-loaded sRNAs and displace other endogenous sRNAs. As a countermeasure, the VSR 2b loaded vsiRNAs and mRNA-derived siRNAs, which affect the expression of the genes they derive from. Additionally, we analyze how vsiRNAs incorporate into the endogenous RNA silencing pathways by exploring their target mRNAs using parallel analysis of RNA end (PARE) sequencing, which allow us to identify vsiRNA-targeted genes genome-wide.

CONCLUSIONS

This work exemplifies the complex relationship of RNA viruses with the endogenous RNA silencing machinery and the multiple aspects of virus resistance and virulence that this interaction induces.

Why are viral genomes so fragile? The bottleneck hypothesis

If they undergo new mutations at each replication cycle, why are RNA viral genomes so fragile, with most mutations being either strongly deleterious or lethal? Here we provide theoretical and numerical evidence for the hypothesis that genetic fragility is partly an evolutionary response to the multiple population bottlenecks experienced by viral populations at various stages of their life cycles. Modelling within-host viral populations as multi-type branching processes, we show that mutational fragility lowers the rate at which Muller’s ratchet clicks and increases the survival probability through multiple bottlenecks. In the context of a susceptible-exposed-infectious-recovered epidemiological model, we find that the attack rate of fragile viral strains can exceed that of more robust strains, particularly at low infectivities and high mutation rates. Our findings highlight the importance of demographic events such as transmission bottlenecks in shaping the genetic architecture of viral pathogens.

Given that most mutations are deleterious, high mutation rates carry a significant evolutionary cost. To reduce this burden, an obvious evolutionary solution would be to reduce the fitness cost of mutations by becoming more robust; this solution is indeed selected in populations of constantly large size. Here, we show that when populations regularly experience bottlenecks, as viruses do upon transmission to a new host, a less obvious solution becomes more viable: namely, to increase the fitness cost of mutations so that unfit mutants are less likely to fix at each passage. This could explain why viruses—especially RNA viruses—do in fact have very fragile genomes.

Frequent asymptomatic infection with tobacco ringspot virus on melon fruit

Melon is one of the most popular fruits worldwide and has been bred into various cultivars. RNA-sequencing using healthy melon fruit was performed to determine differences in gene expression among cultivars. Unexpected RNA-seq results revealed that viruses asymptomatically infected fruits at a high frequency (16 of 21 fruits examined were infected) and that viral transcripts highly accumulated in comparison with host transcripts (15 %-75 % of total reads). Their nucleotide sequences and phylogenetic analyses indicated that more than 10 novel isolates of tobacco ringspot virus (TRSV) were found in melon fruits. Asymptomatic infection with TRSV on melon fruits was confirmed by both immunoblot and RT-PCR analyses. Numerous isolates of TRSV generated and maintained in melon fields, and this is likely due to their asymptomatic infections. This TRSV melon isolate infected Nicotiana benthamiana plants with stunting and yellowing symptoms. This is the first report of frequent and asymptomatic infection of TRSV in consumable melon fruits.

Bottleneck, Isolate, Amplify, Select (BIAS) as a mechanistic framework for intracellular population dynamics of positive-sense RNA viruses

Many positive-sense RNA viruses, especially those infecting plants, are known to experience stringent, stochastic population bottlenecks inside the cells they invade, but exactly how and why these populations become bottlenecked are unclear. A model proposed ten years ago advocates that such bottlenecks are evolutionarily favored because they cause the isolation of individual viral variants in separate cells. Such isolation in turn allows the viral variants to manifest the phenotypic differences they encode. Recently published observations lend mechanistic support to this model and prompt us to refine the model with novel molecular details. The refined model, designated Bottleneck, Isolate, Amplify, Select (BIAS), postulates that these viruses impose population bottlenecks on themselves by encoding bottleneck-enforcing proteins (BNEPs) that function in a concentration-dependent manner. In cells simultaneously invaded by numerous virions of the same virus, BNEPs reach the bottleneck-ready concentration sufficiently early to arrest nearly all internalized viral genomes. As a result, very few (as few as one) viral genomes stochastically escape to initiate reproduction. Repetition of this process in successively infected cells isolates viral genomes with different mutations in separate cells. This isolation prevents mutant viruses encoding defective viral proteins from hitchhiking on sister genome-encoded products, leading to the swift purging of such mutants. Importantly, genome isolation also ensures viral genomes harboring beneficial mutations accrue the cognate benefit exclusively to themselves, leading to the fixation of such beneficial mutations. Further interrogation of the BIAS hypothesis promises to deepen our understanding of virus evolution and inspire new solutions to virus disease mitigation.

Transmission modes affect the population structure of potato virus Y in potato

Transmission is a crucial part of a viral life cycle and transmission mode can have an important impact on virus biology. It was demonstrated that transmission mode can influence the virulence and evolution of a virus; however, few empirical data are available to describe the direct underlying changes in virus population structure dynamics within the host. Potato virus Y (PVY) is an RNA virus and one of the most damaging pathogens of potato. It comprises several genetically variable strains that are transmitted between plants via different transmission modes. To investigate how transmission modes affect the within-plant viral population structure, we have used a deep sequencing approach to examine the changes in the genetic structure of populations (in leaves and tubers) of three PVY strains after successive passages by horizontal (aphid and mechanical) and vertical (via tubers) transmission modes. Nucleotide diversities of viral populations were significantly influenced by transmission modes; lineages transmitted by aphids were the least diverse, whereas lineages transmitted by tubers were the most diverse. Differences in nucleotide diversities of viral populations between leaves and tubers were transmission mode-dependent, with higher diversities in tubers than in leaves for aphid and mechanically transmitted lineages. Furthermore, aphid and tuber transmissions were shown to impose stronger genetic bottlenecks than mechanical transmission. To better understand the structure of virus populations within the host, transmission mode, movement of the virus within the host, and the number of replication cycles after transmission event need to be considered. Collectively, our results suggest a significant impact of virus transmission modes on the within-plant diversity of virus populations and provide quantitative fundamental data for understanding how transmission can shape virus diversity in the natural ecosystems, where different transmission modes are expected to affect virus population structure and consequently its evolution.

Modeling multipartite virus evolution: the genome formula facilitates rapid adaptation to heterogeneous environments†

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.

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.

Prevalence and Distribution of Beet Necrotic Yellow Vein Virus Strains in North Dakota and Minnesota

Beet necrotic yellow vein virus (BNYVV) is the causal agent of rhizomania, a disease of global importance to the sugar beet industry. The most widely implemented resistance gene to rhizomania to date is Rz1, but resistance has been circumvented by resistance-breaking (RB) isolates worldwide. In an effort to gain greater understanding of the distribution of BNYVV and the nature of RB isolates in Minnesota and eastern North Dakota, sugar beet plants were grown in 594 soil samples obtained from production fields and subsequently were analyzed for the presence of BNYVV as well as coding variability in the viral P25 gene, the gene previously implicated in the RB pathotype. Baiting of virus from the soil with sugar beet varieties possessing no known resistance to rhizomania resulted in a disease incidence level of 10.6% in the region examined. Parallel baiting analysis of sugar beet genotypes possessing Rz1, the more recently introgressed Rz2, and with the combination of Rz1 + Rz2 resulted in a disease incidence level of 4.2, 1.0, and 0.8%, respectively. Virus sequences recovered from sugar beet bait plants possessing resistance genes Rz1 and/or Rz2 exhibited reduced genetic diversity in the P25 gene relative to those recovered from the susceptible genotype while confirming the hypervariable nature of the coding for amino acids (AAs) at position 67 and 68 in the P25 protein. In contrast to previous reports, we did not find an association between any one specific AA signature at these positions and the ability to circumvent Rz1-mediated resistance. The data document ongoing virulence development in BNYVV populations to previously resistant varieties and provide a baseline for the analysis of genetic change in the virus population that may accompany the implementation of new resistance genes to manage rhizomania.

Mutagenesis Scanning Uncovers Evolutionary Constraints on Tobacco Etch Potyvirus Membrane-Associated 6K2 Protein

RNA virus high mutation rate is a double-edged sword. At the one side, most mutations jeopardize proteins functions; at the other side, mutations are needed to fuel adaptation. The relevant question then is the ratio between beneficial and deleterious mutations. To evaluate this ratio, we created a mutant library of the 6K2 gene of tobacco etch potyvirus that contains every possible single-nucleotide substitution. 6K2 protein anchors the virus replication complex to the network of endoplasmic reticulum membranes. The library was inoculated into the natural host Nicotiana tabacum, allowing competition among all these mutants and selection of those that are potentially viable. We identified 11 nonsynonymous mutations that remain in the viral population at measurable frequencies and evaluated their fitness. Some had fitness values higher than the wild-type and some were deleterious. The effect of these mutations in the structure, transmembrane properties, and function of 6K2 was evaluated in silico. In parallel, the effect of these mutations in infectivity, virus accumulation, symptoms development, and subcellular localization was evaluated in the natural host. The α-helix H1 in the N-terminal part of 6K2 turned out to be under purifying selection, while most observed mutations affect the link between transmembrane α-helices H2 and H3, fusing them into a longer helix and increasing its rigidity. In general, these changes are associated with higher within-host fitness and development of milder or no symptoms. This finding suggests that in nature selection upon 6K2 may result from a tradeoff between within-host accumulation and severity of symptoms.

In situ hybridization for the localization of two pepino mosaic virus isolates in mixed infections

In situ hybridization (ISH) is an informative and relatively accessible technique for the localization of viral genomes in plant tissue and cells. However, simultaneous visualization of related plant viruses in mixed infections may be limited by the nucleotide similarity in the genomes and the single chromogenic detection over the same sample preparation. To address this issue, we used two Pepino mosaic virus isolates and performed ISH over consecutive serial cross-sections of paraffin-embedded leaf samples of single and mixed infected Nicotiana benthamiana plants. Moreover, the probe design was optimized to reduce cross-hybridisation, and co-localization was based on the overlapping of consecutive cross-sections from mixed infected leaves; thus, our results showed that both Pepino mosaic virus isolates co-localized in the same leaf tissue. In turn, both isolates were localized in the cytoplasm of the same cells. These results provide valuable information for studying mixed infections in plants by using a simple ISH procedure that is accessible to any pathology laboratory.

The Potyvirus Particle Recruits the Plant Translation Initiation Factor eIF4E by Means of the VPg covalently Linked to the Viral RNA

The viral protein genome-linked (VPg) of potyviruses is a protein covalently linked to the 5' end of viral RNA. It interacts with eIF4E, a component of the cellular translation initiation complex. It has been suggested that the 5' RNA-linked VPg could mimic the cellular mRNA cap, promoting synthesis of viral proteins. Here, we report evidence for recruitment of the plant eIF4E by Lettuce mosaic virus (LMV, potyvirus) particles via the 5' RNA-linked VPg. Analysis of the viral population was performed by enzyme-linked immunosorbent assay-based tests, either with crude extracts of LMV-infected tissues or purified viral particles. In both cases, LMV-VPg and LMV-eIF4E subpopulations could be detected. After reaching a maximum within the first 2 weeks postinoculation, these populations decreased and very few labeled particles were found later than 3 weeks postinoculation. The central domain of VPg (CD-VPg) was found to be exposed at the surface of the particles. Using a purified recombinant lettuce eIF4E and CD-VPg-specific antibodies, we demonstrate that the plant factor binds to the VPg via its central domain. Moreover, the plant eIF4E factor could be imaged at one end of the particles purified from LMV plant extracts, by immunoredox atomic force microscopy coupled to scanning electrochemical microscopy. We discuss the biological significance of these results.

Within-host Evolution of Segments Ratio for the Tripartite Genome of Alfalfa Mosaic Virus

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.

Determinants of host species range in plant viruses

Prediction of pathogen emergence is an important field of research, both in human health and in agronomy. Most studies of pathogen emergence have focused on the ecological or anthropic factors involved rather than on the role of intrinsic pathogen properties. The capacity of pathogens to infect a large set of host species, i.e. to possess a large host range breadth (HRB), is tightly linked to their emergence propensity. Using an extensive plant virus database, we found that four traits related to virus genome or transmission properties were strongly and robustly linked to virus HRB. Broader host ranges were observed for viruses with single-stranded genomes, those with three genome segments and nematode-transmitted viruses. Also, two contrasted groups of seed-transmitted viruses were evidenced. Those with a single-stranded genome had larger HRB than non-seed-transmitted viruses, whereas those with a double-stranded genome (almost exclusively RNA) had an extremely small HRB. From the plant side, the family taxonomic rank appeared as a critical threshold for virus host range, with a highly significant increase in barriers to infection between plant families. Accordingly, the plant-virus infectivity matrix shows a dual structure pattern: a modular pattern mainly due to viruses specialized to infect plants of a given family and a nested pattern due to generalist viruses. These results contribute to a better prediction of virus host jumps and emergence risks.

Nonlinear trade-offs allow the cooperation game to evolve from Prisoner's Dilemma to Snowdrift

The existence of cooperation, or the production of public goods, is an evolutionary problem. Cooperation is not favoured because the Prisoner's Dilemma (PD) game drives cooperators to extinction. We have re-analysed this problem by using RNA viruses to motivate a model for the evolution of cooperation. Gene products are the public goods and group size is the number of virions co-infecting the same host cell. Our results show that if the trade-off between replication and production of gene products is linear, PD is observed. However, if the trade-off is nonlinear, the viruses evolve into separate lineages of ultra-defectors and ultra-cooperators as group size is increased. The nonlinearity was justified by the existence of real viral ultra-defectors, known as defective interfering particles, which gain a nonlinear advantage by being smaller. The evolution of ultra-defectors and ultra-cooperators creates the Snowdrift game, which promotes high-level production of public goods.

High virulence does not necessarily impede viral adaptation to a new host: a case study using a plant RNA virus

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.

Matters of Size: Genetic Bottlenecks in Virus Infection and Their Potential Impact on Evolution

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.

Predicting the Stability of Homologous Gene Duplications in a Plant RNA Virus

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.

Ins and Outs of Multipartite Positive-Strand RNA Plant Viruses: Packaging versus Systemic Spread

Viruses possessing a non-segmented genome require a specific recognition of their nucleic acid to ensure its protection in a capsid. A similar feature exists for viruses having a segmented genome, usually consisting of viral genomic segments joined together into one viral entity. While this appears as a rule for animal viruses, the majority of segmented plant viruses package their genomic segments individually. To ensure a productive infection, all viral particles and thereby all segments have to be present in the same cell. Progression of the virus within the plant requires as well a concerted genome preservation to avoid loss of function. In this review, we will discuss the "life aspects" of chosen phytoviruses and argue for the existence of RNA-RNA interactions that drive the preservation of viral genome integrity while the virus progresses in the plant.

A Multiscale Approach to Plant Disease Using the Metacommunity Concept

Plant disease arises from the interaction of processes occurring at multiple spatial and temporal scales. With new tools such as next-generation sequencing, we are learning about the diversity of microbes circulating within and among plant populations and often coinhabiting host individuals. The proliferation of pathogenic microbes depends on single-species dynamics and multispecies interactions occurring within and among host cells, the spatial organization and genetic landscape of hosts, the frequency and mode of transmission among hosts and host populations, and the abiotic environmental context. Here, we examine empirical evidence from these multiple scales to assess the utility of metacommunity theory, a theoretical framework developed for free-living organisms to further our understanding of and assist in predicting plant-pathogen infection and spread. We suggest that deeper understanding of disease dynamics can arise through the application of this conceptual framework at scales ranging from individual cells to landscapes. In addition, we use this multiscale theoretical perspective to synthesize existing knowledge, generate novel hypotheses, and point toward promising future opportunities for the study of plant pathogens in natural populations.

Genetic variation in fitness within a clonal population of a plant RNA virus

A long-standing observation in evolutionary virology is that RNA virus populations are highly polymorphic, composed by a mixture of genotypes whose abundances in the population depend on complex interaction between fitness differences, mutational coupling and genetic drift. It was shown long ago, though in cell cultures, that most of these genotypes had lower fitness than the population they belong, an observation that explained why single-virion passages turned on Muller’s ratchet while very large population passages resulted in fitness increases in novel environments. Here we report the results of an experiment specifically designed to evaluate in vivo the fitness differences among the subclonal components of a clonal population of the plant RNA virus tobacco etch potyvirus (TEV). Over 100 individual biological subclones from a TEV clonal population well adapted to the natural tobacco host were obtained by infectivity assays on a local lesion host. The replicative fitness of these subclones was then evaluated during infection of tobacco relative to the fitness of large random samples taken from the starting clonal population. Fitness was evaluated at increasing number of days post-inoculation. We found that at early days, the average fitness of subclones was significantly lower than the fitness of the clonal population, thus confirming previous observations that most subclones contained deleterious mutations. However, as the number of days of viral replication increases, population size expands exponentially, more beneficial and compensatory mutations are produced, and selection becomes more effective in optimizing fitness, the differences between subclones and the population disappeared.

Multiple Barriers to the Evolution of Alternative Gene Orders in a Positive-Strand RNA Virus

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.

RNA Silencing May Play a Role in but Is Not the Only Determinant of the Multiplicity of Infection

The multiplicity of infection (MOI), i.e., the number of viral genomes that infect a cell, is an important parameter in virus evolution, which for each virus and environment may have an optimum value that maximizes virus fitness. Thus, the MOI might be controlled by virus functions, an underexplored hypothesis in eukaryote-infecting viruses. To analyze if the MOI is controlled by virus functions, we estimated the MOI in plants coinfected by two genetic variants of Tomato bushy stunt virus (TBSV); by TBSV and a TBSV-derived defective interfering RNA (DI-RNA); or by TBSV and a second tombusvirus, Cymbidium ringspot virus (CymRSV). The MOI was significantly larger in TBSV-CymRSV coinfections (~4.0) than in TBSV-TBSV or TBSV-DI-RNA coinfections (~1.7 to 2.2). Coinfections by CymRSV or TBSV with chimeras in which an open reading frame (ORF) of one virus species was replaced by that of the other identified a role of viral proteins in determining the MOI, which ranged from 1.6 to 3.9 depending on the coinfecting genotypes. However, no virus-encoded protein or genomic region was the sole MOI determinant. Coinfections by CymRSV and TBSV mutants in which the expression of the gene-silencing suppressor protein p19 was abolished also showed a possible role of gene silencing in MOI determination. Taken together, these results demonstrate that the MOI is a quantitative trait showing continuous variation and that as such it has a complex determination involving different virus-encoded functions.

IMPORTANCE

The number of viral genomes infecting a cell, or the multiplicity of infection (MOI), is an important parameter in virus evolution affecting recombination rates, selection intensity on viral genes, evolution of multipartite genomes, or hyperparasitism by satellites or defective interfering particles. For each virus and environment, the MOI may have an optimum value that maximizes virus fitness, but little is known about MOI control in eukaryote-infecting viruses. We show here that in plants coinfected by two genotypes of Tomato bushy stunt virus (TBSV), the MOI was lower than in plants coinfected by TBSV and Cymbidium ringspot virus (CymRSV). Coinfections by CymRSV or TBSV with TBSV-CymRSV chimeras showed a role of viral proteins in MOI determination. Coinfections by CymRSV and TBSV mutants not expressing the gene-silencing suppressor protein also showed a role of gene silencing in MOI determination. The results demonstrate that the MOI is a quantitative trait with a complex determination involving different viral functions.

The Multiplicity of Cellular Infection Changes Depending on the Route of Cell Infection in a Plant Virus

The multiplicity of cellular infection (MOI) is the number of virus genomes of a given virus species that infect individual cells. This parameter chiefly impacts the severity of within-host population bottlenecks as well as the intensity of genetic exchange, competition, and complementation among viral genotypes. Only a few formal estimations of the MOI currently are available, and most theoretical reports have considered this parameter as constant within the infected host. Nevertheless, the colonization of a multicellular host is a complex process during which the MOI may dramatically change in different organs and at different stages of the infection. We have used both qualitative and quantitative approaches to analyze the MOI during the colonization of turnip plants by Turnip mosaic virus. Remarkably, different MOIs were observed at two phases of the systemic infection of a leaf. The MOI was very low in primary infections from virus circulating within the vasculature, generally leading to primary foci founded by a single genome. Each lineage then moved from cell to cell at a very high MOI. Despite this elevated MOI during cell-to-cell progression, coinfection of cells by lineages originating in different primary foci is severely limited by the rapid onset of a mechanism inhibiting secondary infection. Thus, our results unveil an intriguing colonization pattern where individual viral genomes initiate distinct lineages within a leaf. Kin genomes then massively coinfect cells, but coinfection by two distinct lineages is strictly limited.

IMPORTANCE

The MOI is the size of the viral population colonizing cells and defines major phenomena in virus evolution, like the intensity of genetic exchange and the size of within-host population bottlenecks. However, few studies have quantified the MOI, and most consider this parameter as constant during infection. Our results reveal that the MOI can depend largely on the route of cell infection in a systemically infected leaf. The MOI is usually one genome per cell when cells are infected from virus particles moving long distances in the vasculature, whereas it is much higher during subsequent cell-to-cell movement in mesophyll. However, a fast-acting superinfection exclusion prevents cell coinfection by merging populations originating from different primary foci within a leaf. This complex colonization pattern results in a situation where within-cell interactions are occurring almost exclusively among kin and explains the common but uncharacterized phenomenon of genotype spatial segregation in infected plants.

Fluorescently Tagged Potato virus Y: A Versatile Tool for Functional Analysis of Plant-Virus Interactions

Potato virus Y (PVY) is an economically important plant virus that infects Solanaceous crops such as tobacco and potato. To date, studies into the localization and movement of PVY in plants have been limited to detection of viral RNA or proteins ex vivo. Here, a PVY N605 isolate was tagged with green fluorescent protein (GFP), characterized and used for in vivo tracking. In Nicotiana tabacum cv. Xanthi, PVY N605-GFP was biologically comparable to nontagged PVY N605, stable through three plant-to-plant passages and persisted for four months in infected plants. GFP was detected before symptoms and fluorescence intensity correlated with PVY RNA concentrations. PVY N605-GFP provided in vivo tracking of long-distance movement, allowing estimation of the cell-to-cell movement rate of PVY in N. tabacum cv. Xanthi (7.1 ± 1.5 cells per hour). PVY N605-GFP was adequately stable in Solanum tuberosum cvs. Désirée and NahG-Désirée and able to infect S. tuberosum cvs. Bintje and Bea, Nicotiana benthamiana, and wild potato relatives. PVY N605-GFP is therefore a powerful tool for future studies of PVY-host interactions, such as functional analysis of viral and plant genes involved in viral movement.

Temporal dynamics of intrahost molecular evolution for a plant RNA virus

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 community ecology of pathogens: coinfection, coexistence and community composition

Disease and community ecology share conceptual and theoretical lineages, and there has been a resurgence of interest in strengthening links between these fields. Building on recent syntheses focused on the effects of host community composition on single pathogen systems, we examine pathogen (microparasite) communities using a stochastic metacommunity model as a starting point to bridge community and disease ecology perspectives. Such models incorporate the effects of core community processes, such as ecological drift, selection and dispersal, but have not been extended to incorporate host-pathogen interactions, such as immunosuppression or synergistic mortality, that are central to disease ecology. We use a two-pathogen susceptible-infected (SI) model to fill these gaps in the metacommunity approach; however, SI models can be intractable for examining species-diverse, spatially structured systems. By placing disease into a framework developed for community ecology, our synthesis highlights areas ripe for progress, including a theoretical framework that incorporates host dynamics, spatial structuring and evolutionary processes, as well as the data needed to test the predictions of such a model. Our synthesis points the way for this framework and demonstrates that a deeper understanding of pathogen community dynamics will emerge from approaches working at the interface of disease and community ecology.

Testing the independent action hypothesis of plant pathogen mode of action: a simple and powerful new approach

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.

Analysis of viral (zucchini yellow mosaic virus) genetic diversity during systemic movement through a Cucurbita pepo vine

Determining the extent and structure of intra-host genetic diversity and the magnitude and impact of population bottlenecks is central to understanding the mechanisms of viral evolution. To determine the nature of viral evolution following systemic movement through a plant, we performed deep sequencing of 23 leaves that grew sequentially along a single Cucurbita pepo vine that was infected with zucchini yellow mosaic virus (ZYMV), and on a leaf that grew in on a side branch. Strikingly, of 112 genetic (i.e. sub-consensus) variants observed in the data set as a whole, only 22 were found in multiple leaves. Similarly, only three of the 13 variants present in the inoculating population were found in the subsequent leaves on the vine. Hence, it appears that systemic movement is characterized by sequential population bottlenecks, although not sufficient to reduce the population to a single virion as multiple variants were consistently transmitted between leaves. In addition, the number of variants within a leaf increases as a function of distance from the inoculated (source) leaf, suggesting that the circulating sap may serve as a continual source of virus. Notably, multiple mutational variants were observed in the cylindrical inclusion (CI) protein (known to be involved in both cell-to-cell and systemic movement of the virus) that were present in multiple (19/24) leaf samples. These mutations resulted in a conformational change, suggesting that they might confer a selective advantage in systemic movement within the vine. Overall, these data reveal that bottlenecks occur during systemic movement, that variants circulate in the phloem sap throughout the infection process, and that important conformational changes in CI protein may arise during individual infections.

Shrinkage of genome size in a plant RNA virus upon transfer of an essential viral gene into the host genome

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.

Rapid fluorescent reporter quantification by leaf disc analysis and its application in plant-virus studies

BACKGROUND

Fluorescent proteins are extraordinary tools for biology studies due to their versatility; they are used extensively to improve comprehension of plant-microbe interactions. The viral infection process can easily be tracked and imaged in a plant with fluorescent protein-tagged viruses. In plants, fluorescent protein genes are among the most commonly used reporters in transient RNA silencing and heterologous protein expression assays. Fluorescence intensity is used to quantify fluorescent protein accumulation by image analysis or spectroscopy of protein extracts; however, these methods might not be suitable for medium- to large-scale comparisons.

RESULTS

We report that laser scanners, used routinely in proteomic studies, are suitable for quantitative imaging of plant leaves that express different fluorescent protein pairs. We developed a microtiter plate fluorescence spectroscopy method for direct quantitative comparison of fluorescent protein accumulation in intact leaf discs. We used this technique to measure a fluorescent reporter in a transient RNA silencing suppression assay, and also to monitor early amplification dynamics of a fluorescent protein-labeled potyvirus.

CONCLUSIONS

Laser scanners allow dual-color fluorescence imaging of leaf samples, which might not be acquired in standard stereomicroscope devices. Fluorescence microtiter plate analysis of intact leaf discs can be used for rapid, accurate quantitative comparison of fluorescent protein accumulation.

Estimation of the in vivo recombination rate for a plant RNA virus

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.