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.

Assessing parallel gene histories in viral genomes.

BACKGROUND: The increasing abundance of sequence data has exacerbated a long known problem: gene trees and species trees for the same terminal taxa are often incongruent. Indeed, genes within a genome have not all followed the same evolutionary path due to events such as incomplete lineage sorting, horizontal gene transfer, gene duplication and deletion, or recombination. Considering conflicts between gene trees as an obstacle, numerous methods have been developed to deal with these incongruences and to reconstruct consensus evolutionary histories of species despite the heterogeneity in the history of their genes. However, inconsistencies can also be seen as a source of information about the specific evolutionary processes that have shaped genomes. RESULTS: The goal of the approach here proposed is to exploit this conflicting information: we have compiled eleven variables describing phylogenetic relationships and evolutionary pressures and submitted them to dimensionality reduction techniques to identify genes with similar evolutionary histories. To illustrate the applicability of the method, we have chosen two viral datasets, namely papillomaviruses and Turnip mosaic virus (TuMV) isolates, largely dissimilar in genome, evolutionary distance and biology. Our method pinpoints viral genes with common evolutionary patterns. In the case of papillomaviruses, gene clusters match well our knowledge on viral biology and life cycle, illustrating the potential of our approach. For the less known TuMV, our results trigger new hypotheses about viral evolution and gene interaction. CONCLUSIONS: The approach here presented allows turning phylogenetic inconsistencies into evolutionary information, detecting gene assemblies with similar histories, and could be a powerful tool for comparative pathogenomics.

Costs of CRISPR-Cas-mediated resistance in Streptococcus thermophilus.

CRISPR-Cas is a form of adaptive sequence-specific immunity in microbes. This system offers unique opportunities for the study of coevolution between bacteria and their viral pathogens, bacteriophages. A full understanding of the coevolutionary dynamics of CRISPR-Cas requires knowing the magnitude of the cost of resisting infection. Here, using the gram-positive bacterium Streptococcus thermophilus and its associated virulent phage 2972, a well-established model system harbouring at least two type II functional CRISPR-Cas systems, we obtained different fitness measures based on growth assays in isolation or in pairwise competition. We measured the fitness cost associated with different components of this adaptive immune system: the cost of Cas protein expression, the constitutive cost of increasing immune memory through additional spacers, and the conditional costs of immunity during phage exposure. We found that Cas protein expression is particularly costly, as Cas-deficient mutants achieved higher competitive abilities than the wild-type strain with functional Cas proteins. Increasing immune memory by acquiring up to four phage-derived spacers was not associated with fitness costs. In addition, the activation of the CRISPR-Cas system during phage exposure induces significant but small fitness costs. Together these results suggest that the costs of the CRISPR-Cas system arise mainly due to the maintenance of the defence system. We discuss the implications of these results for the evolution of CRISPR-Cas-mediated immunity.

Improving the effectiveness of artificial microRNA (amiR)-mediated resistance against Turnip mosaic virus by combining two amiRs or by targeting highly conserved viral genomic regions.

A drawback of recent antiviral therapies based on the transgenic expression of artificial microRNAs (amiRs) is the ease with which viruses generate escape mutations. Here, we show two alternative strategies for improving the effectiveness of resistance in plants. First, we expressed two amiRs complementary to independent targets in the viral genome, and second, we designed amiRs complementary to highly conserved RNA motifs in the viral genome.

Genotypic but not phenotypic historical contingency revealed by viral experimental evolution.

BACKGROUND: The importance of historical contingency in determining the potential of viral populations to evolve has been largely unappreciated. Identifying the constraints imposed by past adaptations is, however, of importance for understanding many questions in evolutionary biology, such as the evolution of host usage dynamics by multi-host viruses or the emergence of escape mutants that persist in the absence of antiviral treatments. To address this issue, we undertook an experimental approach in which sixty lineages of Tobacco etch potyvirus that differ in their past evolutionary history and degree of adaptation to Nicotiana tabacum were allowed to adapt to this host for 15 rounds of within host multiplication and transfer. We thereafter evaluated the degree of adaptation to the new host as well as to the original ones and characterized the consensus sequence of each lineage. RESULTS: We found that past evolutionary history did not determine the phenotypic outcome of this common host evolution phase, and that the signal of local adaptation to past hosts had largely disappeared. By contrast, evolutionary history left footprints at the genotypic level, since the majority of host-specific mutations present at the beginning of this experiment were retained in the end-point populations and may have affected which new mutations were consequently fixed. This resulted in further divergence between the sequences despite a shared selective environment. CONCLUSIONS: The present experiment reinforces the idea that the answer to the question "How important is historical contingency in evolution?" strongly depends on the level of integration of the traits studied. A strong historical contingency was found for TEV genotype, whereas a weak effect of on phenotypic evolution was revealed. In an applied context, our results imply that viruses are not easily trapped into suboptimal phenotypes and that (re)emergence is not evolutionarily constrained.

Multihost experimental evolution of a plant RNA virus reveals local adaptation and host-specific mutations.

For multihost pathogens, adaptation to multiple hosts has important implications for both applied and basic research. At the applied level, it is one of the main factors determining the probability and the severity of emerging disease outbreaks. At the basic level, it is thought to be a key mechanism for the maintenance of genetic diversity both in host and pathogen species. Using Tobacco etch potyvirus (TEV) and four natural hosts, we have designed an evolution experiment whose strength and novelty are the use of complex multicellular host organism as hosts and a high level of replication of different evolutionary histories and lineages. A pattern of local adaptation, characterized by a higher infectivity and virulence on host(s) encountered during the experimental evolution was found. Local adaptation only had a cost in terms of performance on other hosts in some cases. We could not verify the existence of a cost for generalists, as expected to arise from antagonistic pleiotropy and other genetic mechanisms generating a fitness trade-off between hosts. This observation confirms that this classical theoretical prediction lacks empirical support. We discuss the reasons for this discrepancy between theory and experiment in the light of our results. The analysis of full genome consensus sequences of the evolved lineages established that all mutations shared between lineages were host specific. A low degree of parallel evolution was observed, possibly reflecting the various adaptive pathways available for TEV in each host. Altogether, these results reveal a strong adaptive potential of TEV to new hosts without severe evolutionary constraints.

Ultradeep sequencing analysis of population dynamics of virus escape mutants in RNAi-mediated resistant plants.

Plant artificial micro-RNAs (amiRs) have been engineered to target viral genomes and induce their degradation. However, the exceptional evolutionary plasticity of RNA viruses threatens the durability of the resistance conferred by these amiRs. It has recently been shown that viral populations not experiencing strong selective pressure from an antiviral amiR may already contain enough genetic variability in the target sequence to escape plant resistance in an almost deterministic manner. Furthermore, it has also been shown that viral populations exposed to subinhibitory concentrations of the antiviral amiR speed up this process. In this article, we have characterized the molecular evolutionary dynamics of an amiR target sequence in a viral genome under both conditions. The use of Illumina ultradeep sequencing has allowed us to identify virus sequence variants at frequencies as low as 2 × 10(-6) and to track their variation in time before and after the viral population was able of successfully infecting plants fully resistant to the ancestral virus. We found that every site in the amiR-target sequence of the viral genome presented variation and that the variant that eventually broke resistance was sampled among the many coexisting ones. In this system, viral evolution in fully susceptible plants results from an equilibrium between mutation and genetic drift, whereas evolution in partially resistant plants originates from more complex dynamics involving mutation, selection, and drift.

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

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.

Tempo and mode of plant RNA virus escape from RNA interference-mediated resistance.

A biotechnological application of artificial microRNAs (amiRs) is the generation of plants that are resistant to virus infection. This resistance has proven to be highly effective and sequence specific. However, before these transgenic plants can be deployed in the field, it is important to evaluate the likelihood of the emergence of resistance-breaking mutants. Two issues are of particular interest: (i) whether such mutants can arise in nontransgenic plants that may act as reservoirs and (ii) whether a suboptimal expression level of the transgene, resulting in subinhibitory concentrations of the amiR, would favor the emergence of escape mutants. To address the first issue, we experimentally evolved independent lineages of Turnip mosaic virus (TuMV) (family Potyviridae) in fully susceptible wild-type Arabidopsis thaliana plants and then simulated the spillover of the evolving virus to fully resistant A. thaliana transgenic plants. To address the second issue, the evolution phase took place with transgenic plants that expressed the amiR at subinhibitory concentrations. Our results show that TuMV populations replicating in susceptible hosts accumulated resistance-breaking alleles that resulted in the overcoming of the resistance of fully resistant plants. The rate at which resistance was broken was 7 times higher for TuMV populations that experienced subinhibitory concentrations of the antiviral amiR. A molecular characterization of escape alleles showed that they all contained at least one nucleotide substitution in the target sequence, generally a transition of the G-to-A and C-to-U types, with many instances of convergent molecular evolution. To better understand the viral population dynamics taking place within each host, as well as to evaluate relevant population genetic parameters, we performed in silico simulations of the experiments. Together, our results contribute to the rational management of amiR-based antiviral resistance in plants.

Analysis of the Contribution of Intrinsic Disorder in Shaping Potyvirus Genetic Diversity.

Intrinsically disordered regions (IDRs) are abundant in the proteome of RNA viruses. The multifunctional properties of these regions are widely documented and their structural flexibility is associated with the low constraint in their amino acid positions. Therefore, from an evolutionary stand point, these regions could have a greater propensity to accumulate non-synonymous mutations (NS) than highly structured regions (ORs, or 'ordered regions'). To address this hypothesis, we compared the distribution of non-synonymous mutations (NS), which we relate here to mutational robustness, in IDRs and ORs in the genome of potyviruses, a major genus of plant viruses. For this purpose, a simulation model was built and used to distinguish a possible selection phenomenon in the biological datasets from randomly generated mutations. We analyzed several short-term experimental evolution datasets. An analysis was also performed on the natural diversity of three different species of potyviruses reflecting their long-term evolution. We observed that the mutational robustness of IDRs is significantly higher than that of ORs. Moreover, the substitutions in the ORs are very constrained by the conservation of the physico-chemical properties of the amino acids. This feature is not found in the IDRs where the substitutions tend to be more random. This reflects the weak structural constraints in these regions, wherein an amino acid polymorphism is naturally conserved. In the course of evolution, potyvirus IDRs and ORs follow different evolutive paths with respect to their mutational robustness. These results have forced the authors to consider the hypothesis that IDRs and their associated amino acid polymorphism could constitute a potential adaptive reservoir.

The evolutionary genetics of emerging plant RNA viruses.

Over the years, agriculture across the world has been compromised by a succession of devastating epidemics caused by new viruses that spilled over from reservoir species or by new variants of classic viruses that acquired new virulence factors or changed their epidemiological patterns. Viral emergence is usually associated with ecological change or with agronomical practices bringing together reservoirs and crop species. The complete picture is, however, much more complex, and results from an evolutionary process in which the main players are ecological factors, viruses' genetic plasticity, and host factors required for virus replication, all mixed with a good measure of stochasticity. The present review puts emergence of plant RNA viruses into the framework of evolutionary genetics, stressing that viral emergence begins with a stochastic process that involves the transmission of a preexisting viral strain into a new host species, followed by adaptation to the new host.

Comparative analysis of mutational robustness of the intrinsically disordered viral protein VPg and of its interactor eIF4E.

Conformational intrinsic disorder is a feature present in many virus proteins. Intrinsically disordered regions (IDRs) have weaker structural requirement than ordered regions and mutations in IDRs could have a lower impact on the virus fitness. This could favor its exploration of adaptive solutions. The potyviral protein VPg contains IDRs with determinants for adaptation to its host plant. To experimentally assess whether IDRs are more resistant to mutations than ordered regions, the biologically relevant interaction between mutant libraries of both VPg and the eukaryotic translation initiation factor 4E (eIF4E) and their respective wild type partner was examined using yeast two hybrid assay. Our data shows that VPg is significantly more robust to mutations than eIF4E and as such belongs to a particular class of intrinsically disordered proteins. This result is discussed from the standpoint of IDRs involvement in the virus adaptive processes.

Differences in accumulation and virulence determine the outcome of competition during Tobacco etch virus coinfection.

Understanding the evolution of virulence for RNA viruses is essential for developing appropriate control strategies. Although it has been usually assumed that virulence is a consequence of within-host replication of the parasite, viral strains may be highly virulent without experiencing large accumulation as a consequence of immunopathological host responses. Using two strains of Tobacco etch potyvirus (TEV) that show a negative relationship between virulence and accumulation rate, we first explored the evolution of virulence and fitness traits during simple and mixed infections. Short-term evolution experiments initiated with each strain independently confirmed the genetic and evolutionary stability of virulence and viral load, although infectivity significantly increased for both strains. Second, competition experiments between hypo- and hypervirulent TEV strains have shown that the outcome of competition is driven by differences in replication rate. A simple mathematical model has been developed to analyze the dynamics of these two strains during coinfection. The model qualitatively reproduced the experimental results using biologically meaningful parameters. Further analyses of the model also revealed a wide parametric region in which a low-fitness but hypovirulent virus can still outcompete a high-fitness but hypervirulent one. These results provide additional support to the observation that virulence and within-host replication may not necessarily be strongly tied in plant RNA viruses.

Pseudomonas syringae pv. tomato hijacks the Arabidopsis abscisic acid signalling pathway to cause disease.

We have found that a major target for effectors secreted by Pseudomonas syringae is the abscisic acid (ABA) signalling pathway. Microarray data identified a prominent group of effector-induced genes that were associated with ABA biosynthesis and also responses to this plant hormone. Genes upregulated by effector delivery share a 42% overlap with ABA-responsive genes and are also components of networks induced by osmotic stress and drought. Strongly induced were NCED3, encoding a key enzyme of ABA biosynthesis, and the abscisic acid insensitive 1 (ABI1) clade of genes encoding protein phosphatases type 2C (PP2Cs) involved in the regulation of ABA signalling. Modification of PP2C expression resulting in ABA insensitivity or hypersensitivity led to restriction or enhanced multiplication of bacteria, respectively. Levels of ABA increased rapidly during bacterial colonisation. Exogenous ABA application enhanced susceptibility, whereas colonisation was reduced in an ABA biosynthetic mutant. Expression of the bacterial effector AvrPtoB in planta modified host ABA signalling. Our data suggest that a major virulence strategy is effector-mediated manipulation of plant hormone homeostasis, which leads to the suppression of defence responses.