Dynamics of the establishment of systemic Potyvirus infection: independent yet cumulative action of primary infection sites

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.

Temporal expression of defence and susceptibility genes and tospovirus accumulation in capsicum chlorosis virus-infected capsicum

Yolo Wonder (YW) and Warlock (W), two capsicum cultivars that are susceptible to capsicum chlorosis virus (CaCV), were compared in terms of symptom development, tospovirus accumulation, and host gene expression during the first 12 days post infection (dpi). Temporal expression of selected early CaCV-response genes was used to gain insights into plant-virus interactions and to identify potential targets for CaCV control. Symptoms developed faster in YW during the first seven days of infection, while systemic symptoms were similar in both cultivars at 10 and 12 dpi. CaCV accumulation was higher in YW at 7 dpi despite a lower titre at 3 dpi. At 12 dpi, virus accumulation was similar for both cultivars. Symptom development appears to be correlated to virus accumulation over time for both cultivars. Chalcone synthase (CHS), cytochrome P450 (CYP), and tetraspanin 8-like (TSP8) genes followed a similar expression pattern over time in both cultivars. The thionin gene showed increased expression in CaCV-infected plants at 12 dpi. The WRKY40 gene showed significant differential expression at all time points in YW, but only at 12 dpi in W. The strongest correlation of temporal gene expression and virus titre was seen for CYP, TSP8, thionin, and WRKY40. CHS and CYP may be involved in symptom development, and TSP8 may be involved in virus movement. CHS, CYP, and TSP8 may be good targets for future overexpression or silencing studies to clarify their functions during virus infection and, potentially, for control of CaCV in capsicum.

Viruses Infecting Trees and Herbs That Produce Edible Fleshy Fruits with a Prominent Value in the Global Market: An Evolutionary Perspective

Trees and herbs that produce fruits represent the most valuable agricultural food commodities in the world. However, the yield of these crops is not fully achieved due to biotic factors such as bacteria, fungi, and viruses. Viruses are capable of causing alterations in plant growth and development, thereby impacting the yield of their hosts significantly. In this work, we first compiled the world's most comprehensive list of known edible fruits that fits our definition. Then, plant viruses infecting those trees and herbs that produce fruits with commercial importance in the global market were identified. The identified plant viruses belong to 30 families, most of them containing single-stranded RNA genomes. Importantly, we show the overall picture of the host range for some virus families following an evolutionary approach. Further, the current knowledge about plant-virus interactions, focusing on the main disorders they cause, as well as yield losses, is summarized. Additionally, since accurate diagnosis methods are of pivotal importance for viral diseases control, the current and emerging technologies for the detection of these plant pathogens are described. Finally, the most promising strategies employed to control viral diseases in the field are presented, focusing on solutions that are long-lasting.

Significantly Improved Recovery of Recombinant Sonchus Yellow Net Rhabdovirus by Expressing the Negative-Strand Genomic RNA

Generation of recombinant negative-stranded RNA viruses (NSVs) from plasmids involves in vivo reconstitution of biologically active nucleocapsids and faces a unique antisense problem where the negative-sense viral genomic RNAs can hybridize to viral messenger RNAs. To overcome this problem, a positive-sense RNA approach has been devised through expression of viral antigenomic (ag)RNA and core proteins for assembly of antigenomic nucleocapsids. Although this detour strategy works for many NSVs, the process is still inefficient. Using Sonchus yellow net rhabdovirus (SYNV) as a model; here, we develop a negative-sense genomic RNA-based approach that increased rescue efficiency by two orders of magnitude compared to the conventional agRNA approach. The system relied on suppression of double-stranded RNA induced antiviral responses by co-expression of plant viruses-encoded RNA silencing suppressors or animal viruses-encoded double-stranded RNA antagonists. With the improved approach, we were able to recover a highly attenuated SYNV mutant with a deletion in the matrix protein gene which otherwise could not be rescued via the agRNA approach. Reverse genetics analyses of the generated mutant virus provided insights into SYNV virion assembly and morphogenesis. This approach may potentially be applicable to other NSVs of plants or animals.

Genome-wide association mapping of QTLs implied in potato virus Y population sizes in pepper: evidence for widespread resistance QTL pyramiding

In this study, we looked for genetic factors in the pepper (Capsicum annuum) germplasm that control the number of potato virus Y (PVY) particles entering the plant (i.e. effective population size at inoculation) and the PVY accumulation at the systemic level (i.e. census population size). Using genotyping-by-sequencing (GBS) in a core collection of 256 pepper accessions, we obtained 10 307 single nucleotide polymorphisms (SNPs) covering the whole genome. Genome-wide association studies (GWAS) detected seven SNPs significantly associated with the virus population size at inoculation and/or systemic level on chromosomes 4, 6, 9 and 12. Two SNPs on chromosome 4 associated with both PVY population sizes map closely to the major resistance gene pvr2 encoding the eukaryotic initiation factor 4E. No obvious candidates for resistance were identified in the confidence intervals for the other chromosomes. SNPs detected on chromosomes 6 and 12 colocalized with resistance quantitative trait loci (QTLs) previously identified with a biparental population. These results show the efficiency of GBS and GWAS in C. annuum, indicate highly consistent results between GWAS and classical QTL mapping, and suggest that resistance QTLs identified with a biparental population are representative of a much larger collection of pepper accessions. Moreover, the resistance alleles at these different loci were more frequently combined than expected by chance in the core collection, indicating widespread pyramiding of resistance QTLs and widespread combination of resistance QTLs and major effect genes. Such pyramiding may increase resistance efficiency and/or durability.

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.

Impact of genetic drift, selection and accumulation level on virus adaptation to its host plants

The efficiency of plant major resistance genes is limited by the emergence and spread of resistance-breaking mutants. Modulation of the evolutionary forces acting on pathogen populations constitutes a promising way to increase the durability of these genes. We studied the effect of four plant traits affecting these evolutionary forces on the rate of resistance breakdown (RB) by a virus. Two of these traits correspond to virus effective population sizes (N_e ) at either plant inoculation or during infection. The third trait corresponds to differential selection exerted by the plant on the virus population. Finally, the fourth trait corresponds to within-plant virus accumulation (VA). These traits were measured experimentally on Potato virus Y (PVY) inoculated to a set of 84 pepper doubled-haploid lines, all carrying the same pvr2³ resistance gene, but having contrasting genetic backgrounds. The lines showed extensive variation for the rate of pvr2³ RB by PVY and for the four other traits of interest. A generalized linear model showed that three of these four traits, with the exception of N_e at inoculation, and several pairwise interactions between them had significant effects on RB. RB increased with increasing values of N_e during plant infection or VA. The effect of differential selection was more complex because of a strong interaction with VA. When VA was high, RB increased as the differential selection increased. An opposite relationship between RB and differential selection was observed when VA was low. This study provides a framework to select plants with appropriate virus evolution-related traits to avoid or delay RB.

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.

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.

Developments in Plant Negative-Strand RNA Virus Reverse Genetics

Twenty years ago, breakthroughs for reverse genetics analyses of negative-strand RNA (NSR) viruses were achieved by devising conditions for generation of infectious viruses in susceptible cells. Recombinant strategies have subsequently been engineered for members of all vertebrate NSR virus families, and research arising from these advances has profoundly increased understanding of infection cycles, pathogenesis, and complexities of host interactions of animal NSR viruses. These strategies also permitted development of many applications, including attenuated vaccines and delivery vehicles for therapeutic and biotechnology proteins. However, for a variety of reasons, it was difficult to devise procedures for reverse genetics analyses of plant NSR viruses. In this review, we discuss advances that have circumvented these problems and resulted in construction of a recombinant system for Sonchus yellow net nucleorhabdovirus. We also discuss possible extensions to other plant NSR viruses as well as the applications that may emanate from recombinant analyses of these pathogens.

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.

Onset of virus systemic infection in plants is determined by speed of cell-to-cell movement and number of primary infection foci

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.

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

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.

Experimental evolution of pseudogenization and gene loss in a plant RNA virus

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.

Narrow bottlenecks affect Pea seedborne mosaic virus populations during vertical seed transmission but not during leaf colonization

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

Intra-specific variability and biological relevance of P3N-PIPO protein length in potyviruses

BACKGROUND

Pipo was recently described as a new ORF encoded within the genome of the Potyviridae family members (PNAS 105:5897-5902, 2008). It is embedded within the P3 cistron and is translated in the +2 reading frame relative to the potyviral long ORF as the P3N-PIPO fusion protein. In this work, we first collected pipo nucleotide sequences available for different isolates of 48 Potyvirus species. Second, to determine the biological implications of variation in pipo length, we measured infectivity, viral accumulation, cell-to-cell and systemic movements for two Turnip mosaic virus (TuMV) variants with pipo alleles of different length in three different susceptible host species, and tested for differences between the two variants.

RESULTS

In addition to inter-specific variation, there was high variation in the length of the PIPO protein among isolates within species (ranging from 1 to 89 amino acids). Furthermore, selection analyses on the P3 cistron did not account for the existence of stop codons in the pipo ORF, but showed that positive selection was significant in the overlapping region for Potato virus Y (PVY) and TuMV. In some cases, variability in length was associated with host species, geographic provenance and/or other strain features. We found significant empirical differences among the phenotypes associated with TuMV pipo alleles, though the magnitude and sign of the effects were host-dependent.

CONCLUSIONS

The combination of computational molecular evolution analyses and experiments stemming from these analyses provide clues about the selective pressures acting upon the different-length pipo alleles and show that variation in length may be maintained by host-driven selection.

Effects of the number of genome segments on primary and systemic infections with a multipartite plant RNA virus

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.

Stability and fitness impact of the visually discernible Rosea1 marker in the Tobacco etch virus genome

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.

Model-selection-based approach for calculating cellular multiplicity of infection during virus colonization of multi-cellular hosts

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.