Vincent Simon

dN/dS-based methods detect positive selection linked to trade-offs between different fitness traits in the coat protein of potato virus Y.

The dN/dS ratio between nonsynonymous and synonymous substitution rates has been used extensively to identify codon positions involved in adaptive processes. However, the accuracy of this approach has been questioned, and very few studies have attempted to validate experimentally its predictions. Using the coat protein (CP) of Potato virus Y (PVY; genus Potyvirus, family Potyviridae) as a case study, we identified several candidate positively selected codon positions that differed between clades. In the CP of the N clade of PVY, positive selection was detected at codon positions 25 and 68 by both the softwares PAML and HyPhy. We introduced nonsynonymous substitutions at these positions in an infectious cDNA clone of PVY and measured the effect of these mutations on virus accumulation in its two major cultivated hosts, tobacco and potato, and on its efficiency of transmission from plant to plant by aphid vectors. The mutation at codon position 25 significantly modified the virus accumulation in the two hosts, whereas the mutation at codon position 68 significantly modified the virus accumulation in one of its hosts and its transmissibility by aphids. Both mutations were involved in adaptive trade-offs. We suggest that our study was particularly favorable to the detection of adaptive mutations using dN/dS estimates because, as obligate parasites, viruses undergo a continuous and dynamic interaction with their hosts that favors the recurrent selection of adaptive mutations and because trade-offs between different fitness traits impede (or at least slow down) the fixation of these mutations and maintain polymorphism within populations.

The Helper Component Proteinase Cistron of Potato virus Y Induces Hypersensitivity and Resistance in Potato Genotypes Carrying Dominant Resistance Genes on Chromosome IV

The Nc(tbr) and Ny(tbr) genes in Solanum tuberosum determine hypersensitive reactions, characterized by necrotic reactions and restriction of the virus systemic movement, toward isolates belonging to clade C and clade O of Potato virus Y (PVY), respectively. We describe a new resistance from S. sparsipilum which possesses the same phenotype and specificity as Nc(tbr) and is controlled by a dominant gene designated Nc(spl). Nc(spl) maps on potato chromosome IV close or allelic to Ny(tbr). The helper component proteinase (HC-Pro) cistron of PVY was shown to control necrotic reactions and resistance elicitation in plants carrying Nc(spl), Nc(tbr), and Ny(tbr). However, inductions of necrosis and of resistance to the systemic virus movement in plants carrying Nc(spl) reside in different regions of the HC-Pro cistron. Also, genomic determinants outside the HC-Pro cistron are involved in the systemic movement of PVY after induction of necroses on inoculated leaves of plants carrying Ny(tbr). These results suggest that the Ny(tbr) resistance may have been involved in the recent emergence of PVY isolates with a recombination breakpoint near the junction of HC-Pro and P3 cistrons in potato crops. Therefore, this emergence could constitute one of the rare examples of resistance breakdown by a virus which was caused by recombination instead of by successive accumulation of nucleotide substitutions.

Genetic background matters: a plant-virus gene-for-gene interaction is strongly influenced by genetic contexts.

Evolutionary processes responsible for parasite adaptation to their hosts determine our capacity to manage sustainably resistant plant crops. Most plant-parasite interactions studied so far correspond to gene-for-gene models in which the nature of the alleles present at a plant resistance locus and at a pathogen pathogenicity locus determine entirely the outcome of their confrontation. The interaction between the pepper pvr2 resistance locus and Potato virus Y (PVY) genome-linked protein VPg locus obeys this kind of model. Using synthetic chimeras between two parental PVY cDNA clones, we showed that the viral genetic background surrounding the VPg pathogenicity locus had a strong impact on the resistance breakdown capacity of the virus. Indeed, recombination of the cylindrical inclusion (CI) coding region between two PVY cDNA clones multiplied by six the virus capacity to break down the pvr2(3) -mediated resistance. High-throughput sequencing allowed the exploration of the diversity of PVY populations in response to the selection pressure of the pvr2(3) resistance. The CI chimera, which possessed an increased resistance breakdown capacity, did not show an increased mutation accumulation rate. Instead, selection of the most frequent resistance-breaking mutation seemed to be more efficient for the CI chimera than for the parental virus clone. These results echoed previous observations, which showed that the plant genetic background in which the pvr2(3) resistance gene was introduced modified strongly the efficiency of selection of resistance-breaking mutations by PVY. In a broader context, the PVY CI coding region is one of the first identified genetic factors to determine the evolvability of a plant virus.

Modelling the Evolutionary Dynamics of Viruses within Their Hosts: A Case Study Using High-Throughput Sequencing

Uncovering how natural selection and genetic drift shape the evolutionary dynamics of virus populations within their hosts can pave the way to a better understanding of virus emergence. Mathematical models already play a leading role in these studies and are intended to predict future emergences. Here, using high-throughput sequencing, we analyzed the within-host population dynamics of four Potato virus Y (PVY) variants differing at most by two substitutions involved in pathogenicity properties. Model selection procedures were used to compare experimental results to six hypotheses regarding competitiveness and intensity of genetic drift experienced by viruses during host plant colonization. Results indicated that the frequencies of variants were well described using Lotka-Volterra models where the competition coefficients βij exerted by variant j on variant i are equal to their fitness ratio, rj/ri. Statistical inference allowed the estimation of the effect of each mutation on fitness, revealing slight (s = −0.45%) and high (s = −13.2%) fitness costs and a negative epistasis between them. Results also indicated that only 1 to 4 infectious units initiated the population of one apical leaf. The between-host variances of the variant frequencies were described using Dirichlet-multinomial distributions whose scale parameters, closely related to the fixation index F ST, were shown to vary with time. The genetic differentiation of virus populations among plants increased from 0 to 10 days post-inoculation and then decreased until 35 days. Overall, this study showed that mathematical models can accurately describe both selection and genetic drift processes shaping the evolutionary dynamics of viruses within their hosts.

Evolution of plant eukaryotic initiation factor 4E (eIF4E) and potyvirus genome-linked protein (VPg): A game of mirrors impacting resistance spectrum and durability

Polymorphism in the plant eukaryotic translation initiation factor 4E (eIF4E) and potyvirus genome-linked protein (VPg) determine, in many cases, the outcome of the confrontation between these two organisms: compatibility (i.e. infection of the plant by the virus) or incompatibility (i.e. resistance of the plant to the virus). The two interacting proteins eIF4E and VPg show strikingly similar evolution patterns. Most codon positions in their coding sequences are highly constrained for nonsynonymous substitutions but a small number shows evidence for positive selection. Several of these latter positions were shown to be functionally important, conferring resistance to the host or pathogenicity to the virus. Determining the mutational pathways involved in pepper eIF4E diversification revealed a link between an increase of the pepper resistance spectrum towards a panel of potyvirus species and an increase of durability of the resistance towards Potato virus Y. This relationship questions the interest of using more generally the spectrum of action of a plant resistance gene as a predictor of its durability potential.

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.

Interaction patterns between potato virus Y and eIF4E-mediated recessive resistance in the Solanaceae.

ABSTRACT The structural pattern of infectivity matrices, which contains infection data resulting from inoculations of a set of hosts by a set of parasites, is a key parameter for our understanding of biological interactions and their evolution. This pattern determines the evolution of parasite pathogenicity and host resistance, the spatiotemporal distribution of host and parasite genotypes, and the efficiency of disease control strategies. Two major patterns have been proposed for plant-virus genotype infectivity matrices. In the gene-for-gene model, infectivity matrices show a nested pattern, where the host ranges of specialist virus genotypes are subsets of the host ranges of less specialized viruses. In contrast, in the matching-allele (MA) model, each virus genotype is specialized to infect one (or a small set of) host genotype(s). The corresponding infectivity matrix shows a modular pattern where infection is frequent for plants and viruses belonging to the same module but rare for those belonging to different modules. We analyzed the structure of infectivity matrices between Potato virus Y (PVY) and plant genotypes in the family Solanaceae carrying different eukaryotic initiation factor 4E (eIF4E)-coding alleles conferring recessive resistance. Whereas this system corresponds mechanistically to an MA model, the expected modular pattern was rejected based on our experimental data. This was mostly because PVY mutations involved in adaptation to a particular plant genotype displayed frequent pleiotropic effects, conferring simultaneously an adaptation to additional plant genotypes with different eIF4E alleles. Such effects should be taken into account for the design of strategies of sustainable control of PVY through plant varietal mixtures or rotations. IMPORTANCE The interaction pattern between host and virus genotypes has important consequences on their respective evolution and on issues regarding the application of disease control strategies. We found that the structure of the interaction between Potato virus Y (PVY) variants and host plants in the family Solanaceae departs significantly from the current model of interaction considered for these organisms because of frequent pleiotropic effects of virus mutations. These mutational effects allow the virus to expand rapidly its range of host plant genotypes, make it very difficult to predict the effects of mutations in PVY infectivity factors, and raise concerns about strategies of sustainable management of plant genetic resistance to viruses.

Genetic determinism and evolutionary reconstruction of a host jump in a plant virus

In spite of their widespread occurrence, only few host jumps by plant viruses have been evidenced and the molecular bases of even fewer have been determined. A combination of three independent approaches, 1) experimental evolution followed by reverse genetics analysis, 2) positive selection analysis, and 3) locus-by-locus analysis of molecular variance (AMOVA) allowed reconstructing the Potato virus Y (PVY; genus Potyvirus, family Potyviridae) jump to pepper (Capsicum annuum), probably from other solanaceous plants. Synthetic chimeras between infectious cDNA clones of two PVY isolates with contrasted levels of adaptation to C. annuum showed that the P3 and, to a lower extent, the CI cistron played important roles in infectivity toward C. annuum. The three analytical approaches pinpointed a single nonsynonymous substitution in the P3 and P3N-PIPO cistrons that evolved several times independently and conferred adaptation to C. annuum. In addition to increasing our knowledge of host jumps in plant viruses, this study illustrates also the efficiency of locus-by-locus AMOVA and combined approaches to identify adaptive mutations in the genome of RNA viruses.

Estimating virus effective population size and selection without neutral markers

By combining high-throughput sequencing (HTS) with experimental evolution, we can observe the within-host dynamics of pathogen variants of biomedical or ecological interest. We studied the evolutionary dynamics of five variants of Potato virus Y (PVY) in 15 doubled-haploid lines of pepper. All plants were inoculated with the same mixture of virus variants and variant frequencies were determined by HTS in eight plants of each pepper line at each of six sampling dates. We developed a method for estimating the intensities of selection and genetic drift in a multi-allelic Wright-Fisher model, applicable whether these forces are strong or weak, and in the absence of neutral markers. This method requires variant frequency determination at several time points, in independent hosts. The parameters are the selection coefficients for each PVY variant and four effective population sizes Ne at different time-points of the experiment. Numerical simulations of asexual haploid Wright-Fisher populations subjected to contrasting genetic drift (Ne ∈ [10, 2000]) and selection (|s| ∈ [0, 0.15]) regimes were used to validate the method proposed. The experiment in closely related pepper host genotypes revealed that viruses experienced a considerable diversity of selection and genetic drift regimes. The resulting variant dynamics were accurately described by Wright-Fisher models. The fitness ranks of the variants were almost identical between host genotypes. By contrast, the dynamics of Ne were highly variable, although a bottleneck was often identified during the systemic movement of the virus. We demonstrated that, for a fixed initial PVY population, virus effective population size is a heritable trait in plants. These findings pave the way for the breeding of plant varieties exposing viruses to stronger genetic drift, thereby slowing virus adaptation.

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

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 (Ne ) 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 pvr23 resistance gene, but having contrasting genetic backgrounds. The lines showed extensive variation for the rate of pvr23 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 Ne at inoculation, and several pairwise interactions between them had significant effects on RB. RB increased with increasing values of Ne 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.