Introduction
============

Protein intrinsic disorder (ID) defines proteins or protein regions that fail to fold into a unique 3D conformation under physiological conditions ([@msx249-B49]). A growing interest from the scientific community has been observed for this revolutionizing concept in the early 2000s. Large-scale investigations have led to recognize the ubiquity of ID among the whole living world, including viruses ([@msx249-B52]; [@msx249-B37]). Despite of their very limited overall size, some RNA virus proteomes possess the highest ID contents ([@msx249-B53]; [@msx249-B37]). The abundance of ID among RNA viruses could result from a strategy to cope with their very peculiar lifestyle and unique molecular characteristics ([@msx249-B2]; [@msx249-B50]). Indeed, one of the many noteworthy features of RNA viruses is their tiny genomes, currently comprised between 3,500 nt and 30,000 nt, which generally encode no more than a dozen of proteins (<http://viralzone.expasy.org/>; last accessed September 25, 2017; [@msx249-B26]). Still, as obligatory parasites, RNA viruses need to recruit various host cellular components such as biological membranes, proteins and nucleic acids to accomplish their infectious cycles. Consequently, viral proteins are usually multi-functional and interact with multiple ligands. In this context, ID could be highly advantageous for RNA viruses by ensuring binding promiscuity. Indeed, many intrinsically disordered regions (IDRs) possess a high multiplicity of molecular interactions, partly illustrated by the high amount of ID in hub proteins ([@msx249-B24]). Hence, the biophysical and dynamical properties of such regions are favorable to conformational rearrangements and molecular interaction switches, thus facilitating interactions with multiple partners. These extreme levels of structural plasticity could be also favorable to extend virus host range by enabling the recruitment of homologous factors expressed from different host species or genotypes in specific cases of recessive or nonhost resistances, or by enlarging the possibilities of molecular recognition disturbance in case of immune escape ([@msx249-B53]). Beside binding promiscuity, a second advantage of ID in RNA viruses lies in the potential ability of such plastic regions to better accumulate and accommodate mutations. Indeed, IDRs generally display higher evolutionary rates compared with more-structured regions ([@msx249-B7], [@msx249-B6], [@msx249-B5]; [@msx249-B13]; [@msx249-B36]; [@msx249-B30]; [@msx249-B11]). Mechanistically, IDRs are characterized by lower contact densities between amino acid residues, which is expected to alleviate destabilizing effects of amino acid substitutions on viral protein stability and function ([@msx249-B45]; [@msx249-B21]). RNA viruses display extraordinary high mutation rates, comprised between 10^−3^ and 10^−5^ errors/nucleotide/replication cycle ([@msx249-B19]; [@msx249-B42]), which contributes to place them as the fastest evolving entities on Earth and confers a great potential of fast adaptation to environmental disturbances ([@msx249-B18]; [@msx249-B39]). Consequently, the high mutation rates that result in sequence diversity and potentially adaptation of RNA viruses could be, at least partially, permitted by the ID of their proteins.

RNA viruses comprise the most damaging pathogens in both animals and plants. Their strong adaptive potential acts as a substantial barrier in the development of durable and efficient antiviral strategies in both medicine, veterinary and agricultural domains ([@msx249-B16]; [@msx249-B38]; [@msx249-B17]). Understanding the genetic and molecular determinants of their capacity of adaptation still constitutes one of the major challenges in modern virology. The experimental assessment of the virus proteome ID contribution to host adaptation within a natural biological context is still in its beginning. We used the *Potato Virus Y* (PVY)/*Capsicum annuum* (pepper) pathosystem to perform the first experimental analysis of the consequences of ID modifications on the adaptive potential of an RNA virus to its host. PVY, the type member of the large genus *Potyvirus* in the family *Potyviridae*, is one of the major viruses in plant pathology ([@msx249-B43]). In a preliminary study, by performing both inter and intraspecies analyses, we observed that, in potyviruses, ID regions display significantly higher d*N*/d*S* ratios (nonsynonymous and synonymous substitutions rates) than the ordered ones, which indicates the true tendency of IDRs to evolve faster than more structured regions during potyvirus evolution ([@msx249-B11]). Such weaker evolutionary constraints on IDRs strongly suggest the abilities for a faster and easier exploration of adaptive solutions. These results prompted us to further experimentally assess this hypothesis.

PVY infection in pepper involves numerous interactions between virus and host components. The multifunctional Viral genome linked protein (VPg) participates to several of them, and constitutes a major pathogenicity factor in potyviruses ([@msx249-B29]; [@msx249-B20]). The VPg has been experimentally characterized as intrinsically disordered in various potyviruses ([@msx249-B22]; [@msx249-B40]; [@msx249-B14]). VPg IDRs are maintained within the whole *Potyvirus* genus, despite a high level of amino acid polymorphism, strongly suggesting a functional role of these regions in the virus fitness ([@msx249-B25]; [@msx249-B11]). Among these VPg IDRs, a central IDR is localized between residues 95 and 120, and harbors the determinants of PVY adaptation to the well-studied eukaryotic translation initiation factor 4E (eIF4E)-based recessive resistance ([@msx249-B35]; [@msx249-B4][@msx249-B3]). A large majority of antiviral recessive resistances observed in plants involves eIF4E ([@msx249-B47]). The eIF4E is an essential component of the translation initiation complex in eukaryotes and is hijacked by many RNA viruses for the translation and/or replication of their genomes. PVY adaptation to pepper eIF4E-mediated resistance can be viewed as a two-component molecular arm race, for which viral infection is directly linked to VPg-eIF4E molecular interaction ([@msx249-B12]). Thus, the colocalization of the central IDR and resistance-breaking sites in PVY VPg could typically illustrate the ID requirement for viral cycle and/or adaptation process. The novelty of our approach consisted in artificially modulating the ID content of this VPg central domain and analyzing the functional consequences of such alterations on PVY pathogenicity and adaptation to pepper eIF4E-mediated resistance.

Results
=======

Obtention of Three PVY VPg Mutants ID1, ID2, and OD2 Affected in Their ID Content
---------------------------------------------------------------------------------

### Design of Mutation Combinations Predicted as Modulating the ID Level in the Central Region of VPg

Evaluating the impact of ID content on PVY adaptive potential first necessitates that the introduced mutations are not deleterious for VPg basal functionality. To maximize the probability of identifying such mutations, ID predictions were run on a data set of 35 VPg sequences of PVY isolates infecting pepper, collected worldwide ([table 1](#msx249-T1){ref-type="table"}; [supplementary table S2](#sup1){ref-type="supplementary-material"}, [Supplementary Material](#sup1){ref-type="supplementary-material"} online). Amino acid substitutions which, in the natural variability, led to the largest ID variations within the VPg central region, were introduced alone or in combination into PVY SON41p wild-type (WT) isolate (hereinafter named SON41p^WT^; [table 1](#msx249-T1){ref-type="table"}; [supplementary table S2](#sup1){ref-type="supplementary-material"}, [Supplementary Material](#sup1){ref-type="supplementary-material"} online). Very importantly, this work aims to investigate the involvement of the VPg central ID as a favorable flexible environment that could enable the appearance of resistance-breaking mutations, and ultimately PVY adaptation. Thus, positions 101, 105, 115, 119, 121, and 123 which hold the single-residue resistance-breaking sites already linked to PVY adaptation to *pvr2^3^* ([supplementary table S2](#sup1){ref-type="supplementary-material"}, [Supplementary Material](#sup1){ref-type="supplementary-material"} online; [@msx249-B35]; [@msx249-B4][@msx249-B3]) were voluntarily avoided to prevent any bias in the subsequent viral adaptation experiments ([table 1](#msx249-T1){ref-type="table"}). Three VPg mutants, named ID1 (R~94~G; D~95~N; F~100~L), ID2 (R~94~G; D~95~E), and OD2 (E~102~Q; E~114~Q), were thus designed. Corresponding in silico predicted ID changes are reported in [figure 1](#msx249-F1){ref-type="fig"}. For ID1 and ID2 mutants, 12% more residues were predicted as intrinsically disordered than in the WT VPg, with a higher probability of ID in ID1 than in ID2. For the OD2 mutant, PONDR-VLXT predicted 17% less disordered residues than in the WT VPg. Table 1.Identification of PONDR-VLXT Predicted VPg-SON41p Disorder-Affecting Amino Acid Substitutions Based on the VPg PVY (group C) Natural Genetic Diversity.PVY SON41p9394959698100101102105108114115117118119121123IRDIEFSERVETALDHSVPg group C polymorphic sites93V**--**94**Gd**95**E, Nd**96V**--**98D**--**100**Ld**101G**RB**102**Qo**105K**RB\***108I**--**114**Qo**115A, M, R, T**RB**117V**--**118V**--**119S, E, H**RB**121N**RB\***123K, E, N, T**RB\***[^2]

![The ID1, ID2, and OD2 VPg mutants. (*A*) Amino acid sequence of the central region of the VPg of PVY isolate SON41p^WT^ and of its derived mutants ID1, ID2, and OD2. (*B*) Intrinsic disorder prediction (PONDR-VLXT) of the SON41p, and of VPg mutants ID1, ID2, and OD2. Residues with a score \< or \> 0.5 are considered as ordered or intrinsically disordered, respectively.](msx249f1){#msx249-F1}

### Assessment of ID Content in the Mutated VPgs by Circular Dichroism (CD) Spectroscopy

The impact of the introduced mutations on VPg structural properties was experimentally assessed. The CD spectrum of IDRs has a specific signature characterized by a minimum signal observed around 205 nm associated to a shoulder at 220 nm ([@msx249-B46]). Thus, one well-established strategy to compare ID content of proteins consists in determining their respective CD signal~205 nm~/CD signal~220 nm~ ratio, considering that higher ratios correspond to higher global ID contents ([@msx249-B48]). The nucleotide sequences encoding WT, ID1, ID2, and OD2 VPgs were cloned in an *Escherichia coli* heterologous expression system from which the corresponding proteins were produced. For each protein, CD spectra were obtained from 195 to 250 nm ([fig. 2*A*](#msx249-F2){ref-type="fig"}). 205/220 nm ratios confirmed that OD2 was significantly less intrinsically disordered than the WT VPg. Conversely, ID1 and to a lesser extent ID2 were significantly more intrinsically disordered than WT VPg ([fig. 2*B*](#msx249-F2){ref-type="fig"}). Interestingly, these data were in very good agreement with in silico predictions. This particular case illustrates the potential robustness of PONDR-VLXT to predict ID content. It is worth mentioning that the CD signal results from the contribution of the whole VPg secondary structure. However, although it is not possible to formally exclude that the introduced changes in the VPg central region could have long-range effects on the VPg structure, it is more likely that the observed variations in ID levels are located in the VPg central region.

![Circular dichroism signal of the VPg of PVY SON41p and of mutants ID1, ID2, and OD2. (*A*) Molar ellipticities obtained for λ comprised between 195 and 250 nm. (*B*) 205 nm/220 nm molar ellipticity ratios calculated for each VPg. "\*\*" and "\*\*\*" correspond to *P*-values ≤ 0.01 and 0.001, respectively (Kruskal--Wallis test).](msx249f2){#msx249-F2}

Effects of ID-Affecting Mutations on PVY Infectivity and Accumulation in Yolo Wonder Susceptible Pepper
-------------------------------------------------------------------------------------------------------

As previously indicated, analyzing the impact of mutations affecting ID on the adaptive capacity of PVY on resistant hosts first requires to ensure that such mutations do not affect the virus cycle in a susceptible host context. To this end, SON41p^WT^ mutants carrying either ID1, ID2, or OD2 VPg coding sequences (named SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^, respectively) were inoculated in Yolo Wonder pepper plants carrying the *pvr2^+^* susceptibility allele. The numbers of symptomatic plants and the corresponding viral accumulation from two independent assays were determined 30 days postinoculation (dpi) and are presented in [figure 3](#msx249-F3){ref-type="fig"}.

![Infectivity properties of SON41p^WT^ and the PVY mutants SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^ on Yolo Wonder host. (*A*) Number of symptomatic plants observed after inoculation of SON41p^WT^ and the PVY mutants SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^ on Yolo Wonder host for the two independent experiments conducted in 2015 and 2016. Values indicate numbers of infected plants (symptoms and/or virus accumulation) among inoculated plants; (*B*) DAS-ELISA determination of PVY relative accumulation on Yolo Wonder pepper host at 36 days post inoculation. White and grey bars represent data obtained in the 2015 and 2016 experiments, respectively. Within each experiment, no significant differences were observed between the four PVY variants (Mann--Whitney *U* tests; *P*-values \> 0.05). Standard deviations are indicated on each bar.](msx249f3){#msx249-F3}

As for SON41p^WT^, 100% of plants inoculated with SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^ displayed mosaic symptoms in apical noninoculated leaves as soon as 10 dpi for both the 2015 and 2016 assays ([fig. 3*A*](#msx249-F3){ref-type="fig"}). Crucially, DAS-ELISA (double antibody sandwich-enzyme-linked immunosorbent assay) revealed a similar viral accumulation for the four PVY variants in Yolo Wonder plants at 30 dpi ([fig. 3*B*](#msx249-F3){ref-type="fig"}).

To control for the stability of SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^ during infection and to confirm that the virus phenotype was not linked to further evolution within the plants, ten plants per mutant were subjected to total RNA extraction and Sanger sequencing of the PVY VPg cistron. In addition, three plants per mutant were randomly chosen for PVY full-genome sequencing. When compared with the initial sequence, no additional mutation or mutation reversion was detected. These results confirmed that the mutations introduced to create the SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^ mutants are not deleterious for PVY infectivity. More precisely, they do not affect the ability of the virus to infect Yolo Wonder, as attested both by the viral accumulation and symptoms and considering that the consensus sequence of the PVY populations remained unchanged during the infection process.

Effect of ID-Affecting Mutations on PVY Ability to Adapt to the *pvr2^3^*-Mediated Resistance
---------------------------------------------------------------------------------------------

As previously introduced, adaptation of PVY SON41p to eIF4E-mediated *pvr2^3^* resistance is convenient to study the determinants of RNA virus adaptation as it is well-characterized and because of the experimental advantages it provides. Indeed, *pvr2^3^*-resistance breaking by SON41p is fast (necrotic symptoms appear between 20 and 30 dpi) and partial, as around 20--40% of *pvr2^3^* plants become infected after the appearance and fixation of VPg mutations in the PVY population ([@msx249-B4]; [@msx249-B33]). Consequently, it is easily attainable to estimate if an external (environment, plant genetic background) or internal (virus genetic alterations) factor increases or decreases SON41p adaptation capacity. Moreover, it was previously demonstrated that pepper *pvr2^3^*-mediated resistance and the corresponding PVY adaptation is a two-component system, in which infection is determined mostly by the physical interaction between plant eIF4E and the viral VPg ([@msx249-B12]; [@msx249-B34]). The HD285 pepper line, carrying the *pvr2^3^* resistance gene, essentially differs from Yolo Wonder at the *pvr2* locus (see Materials and Methods) and allows us to consider that the differences in PVY infectivity observed between Yolo Wonder and HD285 are essentially due to differences in binding affinity between PVY VPg and plant *pvr2*-encoded eIF4Es.

To assess the consequence of altering the ID of the VPg central domain on the PVY ability to adapt to the *pvr2^3^* resistance, SON41p^WT^, SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^ were inoculated to 120 HD285 plants. The capacity of WT PVY and PVY mutants to breakdown the *pvr2*^3^ resistance (i.e., to accumulate and induce symptoms on HD285) were compared at 36 dpi. Results of the two independent 2015 and 2016 experiments are presented in [figure 4](#msx249-F4){ref-type="fig"}. Infection by SON41p was observed in 22% (2015) to 39% (2016) of *pvr2*^3^-inoculated plants, in good agreement with previously published data ([@msx249-B4]). Very interestingly, the SON41p^OD2^ mutant, carrying an ID-depleted VPg, was unable to overcome the *pvr2*^3^ resistance, as no symptomatic and/or DAS-ELISA positive plants were observed among the 120 and 50 inoculated. SON41p^ID2^ mutant displayed a significantly higher capacity of infection compared with PVY SON41p^WT^, with between 1.6 and 2.3 times more infected plants. Finally, SON41p^ID1^ mutant reached the highest infectivity, as 2.3--3 times more of inoculated plants were infected ([fig. 4*A*](#msx249-F4){ref-type="fig"}) compared with SON41p^WT^.

![Infectivity properties of SON41p^WT^ and the PVY mutants SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^ on HD285 pepper hosts. (*A*) Normalized adaptive rates observed after inoculation of the PVY mutants SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^ on HD285 host compared with SON41p^WT^. (*B*) Total numbers of infected HD285 plants after inoculation of the SON41p^WT^ and PVY mutants SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^ for the two independent 2015 and 2016 experiments. Values in brackets indicate numbers of infected plants (symptoms and/or virus accumulation) among inoculated plants and the corresponding percentages. a--d: groups with statistically significant differences independently obtained within 2015 (left) and 2016 (right) experiments (Fisher exact tests using 2 × 2 contingency tables and Bonferroni correction for multiples testing; *P*-values \< 0.001).](msx249f4){#msx249-F4}

Analysis of Molecular Events Underlying Infection of PVY SON41p^WT^, SON41p^ID1^, SON41p^ID2^, and SON41p^OD2^ in *pvr2*^3^-Carrying Pepper Plants
--------------------------------------------------------------------------------------------------------------------------------------------------

### Analysis of VPg Sequence Modifications in *pvr2^3^*-Carrying Plants

In a view to understand the large differences of infectivity in HD285 observed between PVY SON41p^WT^ and the two mutants SON41p^ID1^ and SON41p^ID2^, total RNAs were extracted from 10 to 20 randomly-chosen plants infected, for each of the three viruses. Sanger sequencing of the full virus genome and/or VPg-cistron were conducted (see Materials and Methods). No mutation was observed outside the VPg cistron. VPg sequencing results are provided in [figure 5](#msx249-F5){ref-type="fig"} and refer to sequences obtained in the 2015 and 2016 assays. The VPg of SON41p^ID1^ mutant remained nonmutated in 94% (2015) and 90% (2016) of the sequenced samples ([fig. 5*A*](#msx249-F5){ref-type="fig"}). By contrast, SON41p^WT^ and SON41p^ID2^ progenies displayed a minority of nonmutated VPg representing between 0% and 22% of the viral populations ([fig. 5*A*](#msx249-F5){ref-type="fig"}). Unlike SON41p^ID1^, infections of HD285 by SON41p^WT^ and SON41p^ID2^ mutant are largely associated to additional mutations in the VPg.

![Frequencies of VPg amino acid substitutions associated to infection of *pvr2^3^*-carrying HD285 pepper plants by PVY SON41p^WT^ and SON41p^ID2^. (*A*) Proportion of nonmutated VPg associated to infection of *pvr2^3^*-carrying plants after inoculation with SON41p^WT^, SON41p^ID1^, or SON41p^ID2^ PVY mutants, among the total number of sequenced PVY populations obtained from individual infected plants. "a" and "b": groups with statistically-significant differences independently obtained within 2015 and 2016 experiments using Fisher exact tests and Bonferroni correction for multiple testing (2 × 2 contingency tables; *P*-values \< 0.001). ND: no virus population was sequenced. (*B*) Amino acid substitutions observed in virus populations derived from SON41p^WT^ PVY in infected HD285 plants. Dotted bars refer to observations reported from ([@msx249-B4]; [@msx249-B33]). (*C*) Amino acid substitutions observed in virus populations derived from SON41p^ID2^ PVY in infected HD285 plants.](msx249f5){#msx249-F5}

Regarding mutated VPgs in the progenies of SON41p^WT^, four different amino acid substitutions were identified. The majority of these mutations were reported previously ([fig. 5*B*;](#msx249-F5){ref-type="fig"}[@msx249-B4]; [@msx249-B33]) and two of them were shown to be responsible for adaptation of SON41p^WT^ to *pvr2^3^* by reverse genetic experiments (mutations D~119~N and T~115~K). Worth mentioning, S~120~T had never been described so far and may constitute a new *pvr2^3^*-resistance-breaking mutation. Our sequencing results also confirmed the large proportion of D~119~N among SON41p^WT^ resistance-breaking variants as previously reported ([@msx249-B4]; [@msx249-B33]).

As for SON41p^WT^, a large proportion (40--50%) of SON41p^ID2^ resistance-breaking variants carried the D~119~N substitution. Other mutations, observed in few samples, have already been validated for their role in adaptation to *pvr2^3^* by reverse genetics (T~115~K, T~115~R, and S~101~G) or observed in previous experiments (D~119~G). In addition, four mutations which have never been described so far in SON41p^WT^ and were observed in PVY progenies, namely D~119~H, D~119~E, S~120~I, and D~119~N + S~120~I. Although these results reveal a tendency of PVY mutant SON41p^ID2^ to display more diversified adaptive solutions (eight different substitutions in total) than those identified for SON41p^WT^ (four in total), a Shannon diversity comparison showed that this difference was not statistically significant ([supplementary fig. S1](#sup1){ref-type="supplementary-material"}, [Supplementary Material](#sup1){ref-type="supplementary-material"} online).

### Analysis of the Thermodynamic Binding Strengths of the eIF4E/VPg Complex

Analysis of mutations in the VPg from *pvr2^3^*-resistance breaking variants was not sufficient to understand the strong differences observed in infectivity between SON41p^ID1^, SON41p^ID2^, SON41p^WT^, and SON41p^OD2^. An attempt was made to relate these differences to the binding strength between the mutated VPgs and either the susceptible or resistant form of eIF4E, respectively carried by Yolo Wonder (*pvr2^+^/pvr2^+^*) and HD285 (*pvr2^3^/pvr2^3^*) peppers. Steady-state fluorescence spectroscopy was used to measure dissociation constants (*K*~d~) between *pvr2^+^*-eIF4E or *pvr2^3^*-eIF4E and VPgs WT, ID1, ID2, and OD2, according to [@msx249-B32]. The *K*~d~ values of the VPg/eIF4E interactions were measured from two independent assays (see Materials and Methods) and are expressed relative to *K*~d~ obtained for the SON41p^WT^ VPg/*pvr2^+^* eIF4E combination conducted in the same experiment ([fig. 5](#msx249-F5){ref-type="fig"}). To ensure that the two eIF4E forms were stable and fully functional, their ability to interact with a eukaryotic 5′-cap analog m7-GTP was systematically assessed and provided a *K*~d~ in the same range (around 300 nM) than those observed previously for lettuce eIF4E ([@msx249-B32]; *data not shown*), attesting the integrity of eIF4E samples used in these experiments.

The affinity of the WT VPg for *pvr2^3^*-eIF4E was about two times lower than for *pvr2^+^*-eIF4E ([fig. 6](#msx249-F6){ref-type="fig"}). These results could explain the requirement of additional mutations in the VPg for PVY adaptation to *pvr2^3^*resistance and demonstrate, once again, the tight link between VPg/eIF4E interaction and PVY infection in pepper. More precisely, the pepper *pvr2^3^* resistance to PVY could be due to a low affinity between eIF4E and VPg, preventing PVY to establish infection in *pvr2^3^* peppers rather than a total inability to interact with *pvr*2^3^-encoded eIF4E. Additionally, it could also explain why PVY SON41p^WT^ is able to evolve in the *pvr*2^3^-carrying plants and reach adaptive solutions during the 36 days following inoculation. Conversely, the highly intrinsically-disordered ID1 VPg interacted with *pvr2^3^*-eIF4E with an affinity 3 times higher than SON41p^WT^ VPg, while maintaining its ability to recruit the *pvr2^+^*-eIF4E form ([fig. 6*A*](#msx249-F6){ref-type="fig"}). ID2 and WT VPgs displayed the same affinity with the *pvr2^+^*-eIF4E form but ID2 interacted with *pvr2^3^*-eIF4E with a lower affinity than that observed for WT VPg ([fig. 6*B*](#msx249-F6){ref-type="fig"}).

![Relative in vitro dissociation constant (*K*~d~) between *pvr2*-encoded pepper host eIF4E and wild-type (SON41p) or mutated (ID1, ID2, and OD2) PVY VPgs. *Left*: interaction with *pvr2^+^*-encoded eIF4E. *Right*: interaction with *pvr2^3^*-encoded eIF4E. (*A*) Normalized *K*~d~ obtained after comparing interaction with SON41p (grey) and ID1 (red). (*B*) Normalized *K*~d~ obtained after comparing interaction with SON41p (grey) and ID2 (pink) VPgs. (*C*) Normalized *K*~d~ obtained after comparing interaction with SON41p (grey) and OD2 (blue) VPgs. "\*" and "*ns"*: statistically significant and nonsignificant differences, respectively (based on 95% confidence intervals). Standard deviations obtained from the method used to determine each *K*~d~ (see Materials and Methods) are shown on each bar. *Note that K~d~ is inversely related with binding affinity.*](msx249f6){#msx249-F6}

Concerning the more structurally-ordered OD2 VPg, we noticed that affinities of its in vitro interactions with each *pvr2^+^* and *pvr2^3^* eIF4Es were similar to those of the WT VPg ([fig. 6*C*](#msx249-F6){ref-type="fig"}).

Discussion
==========

Assumptions regarding the involvement of protein ID in viral adaptive capacity and adaptation process still rely on in silico predictions and functional hypotheses. To date, no experimental data has been reported in the literature to challenge these assumptions. This work intends to test experimentally the contribution of the ID of a viral protein in the virus host adaptation. The PVY-pepper couple is remarkably convenient for this, as the effect of disorder modulations within the viral VPg on virus infectivity can be followed on a large set of host plants under true biological conditions. We generated three PVY mutants with a modified amount of ID in the VPg pathogenicity factor. Remarkably, and for the first time, it was observed that ID modifications were linked to strong differences in virus infectivity and adaptive behavior.

Causal Link between ID Modifications and Alteration in Virus Phenotypes
-----------------------------------------------------------------------

ID is strictly sequence-dependent. Accordingly, every attempt to modify the ID within a protein inevitably entails a change in its amino acid sequence. In our approach to artificially modify ID within the VPg, a multifunctional protein involved in various steps of PVY cycle, it cannot a priori be excluded that some phenotypic observations result more from the very nature of amino acid substitutions than from the ID changes they induced. The fact that no effect of the introduced mutations was observed on PVY pathogenicity in the susceptible host Yolo Wonder strongly suggests that there is no bias of this kind. The choice of the mutations to be introduced was intended to uncouple the ID itself from previously identified resistance-breaking codon positions. In the case of ID1 and ID2, all mutations were chosen outside of the VPg region involved in resistance breakdown which spans from amino acid position 101 to 123 ([@msx249-B4]; [@msx249-B33]). Conversely, the two mutations conferring an order increase to OD2 were located within the resistance breaking region ([table 1](#msx249-T1){ref-type="table"}). In that case, we postulated that substitutions at positions that have never been observed as associated to resistance-breaking events were unlikely to induce resistance-breaking by their own. In fact, the opposite effect was observed, since SON41p^OD2^ was completely unable to infect the HD285 pepper genotype carrying the pvr2^3^ resistance. Thus, even if a direct causal link between ID changes and virus phenotypes alterations cannot be formally established, our experimental strategy leads us to postulate that the mutant phenotypes observed in this study are more than likely induced by the underlying ID changes.

Is ID an Expander of VPg Binding Spectrum Leading to an Enlargement of PVY Host Range?
--------------------------------------------------------------------------------------

The three amino acid substitutions in ID1 VPg central region have significantly increased the ID of the VPg and are associated to a sharp improvement of SON41p^ID1^ infectivity toward *pvr2^3^*-carrying pepper plants compared with SON41p^WT^. The lack of additional mutations in the VPg (or elsewhere in the genome) of most of the SON41p^ID1^ progenies indicates that the increased ability of SON41p^ID1^ PVY to infect these plants does not rely on adaptive mutations that would occur after inoculation in resistant host but more likely on the intrinsic capacity of SON41p^ID1^ to infect the *pvr2^3^*plants. Considering the direct link that exists between VPg/eIF4E binding and infection in the PVY/pepper pathosystem ([@msx249-B12]), the significant increase of the ID1 VPg affinity for both *pvr2^+^*-eIF4E and *pvr2^3^-*eIF4E compared with WT VPg could largely be responsible for the higher SON41p^ID1^ infectivity. In addition, it could explain that the host range expansion was not traded-off against a decrease of infectivity in Yolo Wonder host. Through its inherent structural flexibility, the VPg central IDR could expand its interaction spectrum to eIF4E forms of closely-related pepper cultivars according to mechanistic principles of ID functionality in hub proteins ([@msx249-B24]). Ultimately, this would lead to an enlargement of the virus infection spectrum.

Is ID a Fuel for PVY Resistance Breakdown in *pvr2^3^*-Carrying Plants?
-----------------------------------------------------------------------

As ID1, VPg mutant ID2 is characterized by an increase of ID level in its central IDR, associated to a significant increase of resistance-breaking events in *pvr2^3^* host compared with SON41p^WT^. However, molecular mechanisms underlying those observations seem radically distinct from those explaining the phenotype of SON41p^ID1^. In the case of SON41p^ID2^ in *pvr2^3^* plants, infection cases are near-exclusively associated to the appearance of additional mutations in the VPg, similar to those previously reported as responsible for the breakdown of the *pvr2^3^* resistance by SON41p^WT^ ([@msx249-B4]; [@msx249-B33]). In addition, the ID2 VPg affinity for *pvr2^+^*- and *pvr2^3^*-eIF4Es are respectively similar and lower than those measured for WT-VPg. These binding strength discrepancies cannot account for the significant differences of resistance-breaking rates observed between SON41p^WT^ and SON41p^ID2^. Considering the general potential of IDRs to evolve faster and to better accommodate mutations, a hypothesis to explain the higher adaptive capacities of SON41p^ID2^ compared with SON41p^WT^ in the *pvr2^3^*-carrying host could be that the higher ID in the VPg of SON41p^ID2^ is associated to a larger mutational robustness of the VPg central region. Such an effect would allow 1) an enhanced tolerance to the resistance-breaking mutations already associated to SON41p^WT^ adaptation; 2) the emergence of additional adaptive solutions. In both these nonexclusive mechanisms, ID increase would lead to a higher adaptive potential. Despite its basal ability to replicate in *pvr2^3^*-carrying host HD285, PVY needs to evolve and accumulate mutations in its VPg in view of maintaining itself and restoring infection in the plant. Consequently, adaptation of the virus lies in its ability to accumulate mutations in the VPg and gain sufficient fitness for efficiently colonizing the host plant while countering its basal defenses. Whereas the kinetics of mutation appearance relies exclusively on the viral replicase mutational rate, the relative fitness of each variant can result from structural properties of their pathogenicity factor---that is, the VPg in the case of PVY adaptation to *pvr2^3^*. Thus, the higher mutational robustness in IDRs could lead resistance-breaking variants to accumulate faster in the host, leading to a higher probability of restoring systemic infection and ultimately of adapting. As a matter of fact, although statistically nonsignificant, it can be observed that, for a similar number of variants sequenced in each study (comprised between 15 and 22), SON41p^ID2^ displays a tendency to accumulate more diverse resistance-breaking variants than SON41p^WT^.

Finally, the analysis of the SON41p^OD2^ PVY mutant remarkably reinforces our hypothesis. Indeed, the decrease of ID generated by the E~102~Q and E~114~Q substitutions introduced into OD2 VPg did neither affect PVY infectivity in Yolo Wonder nor the VPg affinity for *pvr2^+^* or *pvr2^3^*-encoded eIF4Es, compared with SON41p^WT^ and WT VPg, respectively. Consequently, the inability of SON41p^OD2^ to adapt to the *pvr2^3^* resistance results most likely from a disability of this mutant to accumulate adaptive mutations during its evolution in the resistant host, rather than from a deleterious effect of the E~102~Q and E~114~Q mutations on its intrinsic infection capacity. By reducing the VPg ID, the mutations introduced in OD2 could reduce SON41p^OD2^ tolerance to amino acid substitutions, thereby reducing the probability of adaptive events to occur.

Is ID in VPg Central Region Dispensable for PVY Infectious Cycle?
-----------------------------------------------------------------

It was observed in this study that the decrease of ID in the central region of the VPg mutant OD2 does not seem to affect PVY infectivity in the Yolo Wonder pepper plants carrying the susceptibility allele *pvr2^+^*. One of the functional aspects of IDRs is their ability to ensure interactions with multiple molecular partners ([@msx249-B24]; [@msx249-B15]). The well-established multifunctional and multi-interacting nature of the VPg protein ([@msx249-B29]; [@msx249-B20]; [@msx249-B31]) and more specifically of its central domain ([@msx249-B41]) appears as not ultimately dependent on a highly defined ID, despite the evolutionary conservation of this structural feature among potyviruses ([@msx249-B11]). Nevertheless, functional involvement of ID in VPg multiple interactions still remains to be elucidated, as ID is also predominant and conserved in others VPg regions (N-ter and C-ter) and could be responsible for VPg multiple interaction capacities.

Establishing ID Content of Viral Proteins as an Additional Criterion in the Design of Antiviral Strategies
----------------------------------------------------------------------------------------------------------

This study reports a body of information strongly arguing for an impact of ID level of viral proteins on the virus adaptation to host resistance. In order to assess the extent of ID-driven adaptation in virus biology, similar experiments should be performed for other viral proteins and other cases of virus adaptation. Beyond their advantages compared with animal viruses in terms of biosafety, phytoviruses are particularly suited for evolutionary experiments as the biological properties of mutants can be quickly and accurately tested against a rich collection of host plants and the evolution of viral progenies can be easily analyzed in vivo without the ethical restrictions linked to inoculation of animals with genetically-modified viruses. In potyviruses, VPg does not constitute the only pathogenicity factor displaying resistance-breaking determinants associated to evolutionary-conserved IDRs. Lettuce mosaic virus (LMV) resistance-breaking sites have been located in a highly-conserved IDR ([@msx249-B11]) of the virus helicase ([@msx249-B1]). Investigating more distant pathosystems appears as essential and could also include sobemoviruses for which adaptive mutations are also located within IDRs ([@msx249-B25]). Demonstrating the ubiquitous nature of ID as a determinant of RNA virus adaptation would be of major importance for the whole virologist community and is expected to greatly enhance our global understanding of RNA virus adaptation processes and our ability to anticipate those events. At the molecular levels, host resistance durability against viruses results both from the number of mutations necessary for the virus to adapt and the fitness penalty associated to each of them ([@msx249-B23]; [@msx249-B28]; [@msx249-B3]). Resistance durability has been correlated to evolutionary constraints exerted on pathogenicity factors involved in viral adaptation ([@msx249-B27]). Measuring evolutionary rates of pathogenicity factors appears to be an efficient measure to predict the durability of plant resistance and more generally of antiviral control methods. Ultimately, establishing a more general link between ID and RNA virus adaptation could promote this easily predictable structural feature as a new criterion to include in the design of antiviral control methods. Indeed, owing to their adaptive abilities, viral IDRs could constitute poor targets of both drugs and host genetic resistance that constitute the main antiviral strategies in animal and plants.

Materials and Methods
=====================

Plant and Virus Materials
-------------------------

PVY isolate SON41p (accession number AJ439544) was used as the reference WT genotype. All PVY mutants were derived from this genotype. Natural hosts *Nicotiana benthamiana* and *C. annuum* accessions Yolo Wonder and HD285 carrying respectively *pvr2^+^/pvr2^+^* and *pvr2^3^/pvr2^3^* genotypes were used for experiments. HD285 is a doubled-haploid line derived from the F~1~ hybrid (Yolo Wonder × Perennial) and genetically selected for PVY SON41p resistance allele pvr2^3^ and susceptibility alleles at the other loci ([@msx249-B8]). The 35 PVY isolate sequences used for designing VPg mutations belong to PVY group C1 that infects *C. annuum* are listed in [supplementary data](#sup1){ref-type="supplementary-material"} ([supplementary table S1](#sup1){ref-type="supplementary-material"}, [Supplementary Material](#sup1){ref-type="supplementary-material"} online).

Computational and Experimental ID Characterization
--------------------------------------------------

In silico preliminary predictions of VPg ID content were performed using Predictor of Naturally Disordered Regions (PONDR-VLXT), accessible through the Disprot server (<http://disorder.compbio.iupui.edu/metapredictor.php>; last accessed September 25, 2017; [@msx249-B44]; [@msx249-B51]). Parameters were set to "default" for ID score predictions.

In vitro experimental characterization of ID was assessed using far UV CD measurements on the different recombinant VPg forms. Molar ellipticities were recorded at 20 °C under nitrogen using a Jasco-815 spectropolarimeter. Samples consisted in 0.1 mg/ml pure VPg forms in 0.02 M sodium phosphate, 0.3 M sodium fluoride, pH 8. In these conditions, using 1 mm optical path cuvette, the photomultiplier tube dynode voltage did not exceed 380 V at 190 nm, which permitted a satisfactory spectral resolution from 200 nm and above.

PVY Inoculation and Detection
-----------------------------

Infectious mutated PVY clones were obtained after NotI-linearization of a plasmid containing the sequence of SON41p^WT^ except the central part of the VPg-coding sequence followed by homologous recombination in yeast (*Saccharomyces cerevisiae* strain YPH501) with PCR products corresponding to ID1, ID2 or OD2 VPg-coding sequences, as previously described in [@msx249-B4]). Infectious clones were recovered from yeast and cloned into *E. coli*, strain Max-efficiency DH10B (New England Biolab, NEB), to ensure sufficient amounts of plasmids. The direct biolistic inoculation of full length PVY cDNA clones does not produce infection in *C. annuum* ([@msx249-B35]), and a preliminary inoculation in three-week old *N. benthamiana* plants was required. Inoculated plants were grown under controlled conditions and symptomatic leaves were retrieved at 15 dpi ([supplementary fig. S2](#sup1){ref-type="supplementary-material"}, [Supplementary Material](#sup1){ref-type="supplementary-material"} online). To ensure that the virus population did not evolve in the *N. benthamiana* plants used to prepare inocula, full-genome Sanger sequencing was performed on total RNA isolated from each leaf extracts used to inoculate pepper plants. In each case, the initial PVY sequences were found, precluding any major virus evolution within *N. benthamiana*. Viral loads were determined using quantitative DAS-ELISA semiquantitative assays ([@msx249-B4]) using home-made polyclonal antiPVY antibodies. When necessary, the volume of leaf extracts was adjusted with phosphate inoculation buffer to ensure the same virus concentration between virus inocula. The resulting calibrated inocula were used to mechanically inoculate cotyledons of three-weeks-old Yolo Wonder and HD285 peppers, as described in [@msx249-B35]. The effects of virus inoculation on Yolo Wonder were assessed after 30 days, and, in the case of HD285, between 36 and 45 days after inoculation. A successful infection was associated with mosaic and necrotic symptoms on upper noninoculated leaves for HD285 and mosaic only for Yolo Wonder ([supplementary fig. S2](#sup1){ref-type="supplementary-material"}, [Supplementary Material](#sup1){ref-type="supplementary-material"} online). The viral load in Yolo Wonder plants was also assessed by quantitative DAS-ELISA.

PVY VPg and Pepper eIF4E Cloning
--------------------------------

PVY SON41p^WT^ and Yolo Wonder and HD285 eIF4E coding-sequences were introduced into Gateway pDONR-201 entry vector by performing attB-PCR with High-fidelity Phusion polymerase (NEB) and BP recombination reaction according to manufacturer's recommendations. Corresponding gene sequences and attB-primers are respectively listed in [supplementary file S1](#sup1){ref-type="supplementary-material"} and table S3, [Supplementary Material](#sup1){ref-type="supplementary-material"} online.

PVY VPg Site-Directed Mutagenesis
---------------------------------

Nucleotide substitutions corresponding to ID1, ID2, and OD2 VPg variants (A~280~G/T~285~G/T~300~A, A~280~G/T~285~A, G~304~C/G~340~C, respectively) were introduced into the VPg-coding sequence of SON41p^WT^ with the help of primers containing mutations (listed in [supplementary table S3](#sup1){ref-type="supplementary-material"}, [Supplementary Material](#sup1){ref-type="supplementary-material"} online), using Q5 site-directed mutagenesis kit (NEB). Mutated plasmids were recovered from bacterial clones using Wizard Plus SV Miniprep DNA Purification System kit (Promega) according to the manufacturer\'s instructions.

PVY VPg and Pepper eIF4E Production
-----------------------------------

Heterologous expression in *E. coli* was used to produce VPg and eIF4E proteins. VPg and eIF4E coding-sequences were cloned downstream from a hexa-histidine tag into bacterial expression plasmids pDEST17 Gateway vector (Invitrogen) and pETGB1a, respectively, following LR recombination and restriction/ligation standard procedures. Resulting vectors were used to transform expression bacterial strain BL21-AI (Invitrogen). Bacteria were grown at 37 °C until a DO~600nm~ value of 0.5. For protein expression induction, 0.1 M arabinose (pDEST17) and 0.2 M IPTG (pETGB1a) were added during 3 h at 28 °C. Bacterial culture lysates were prepared according to [@msx249-B32]. Total proteins were recovered in a Phosphate/NaF (pH 8) extraction buffer to prevent any far-UV background signal during CD measurements. Proteins of interest were purified by metal ion affinity chromatography, using Ni-NTA columns (Protino). Protein purity was checked by SDS-PAGE on 12.5% bis-acrylamide gels, followed by Coomasie blue total protein revelation and western blot using mouse antipolyHistidine antibodies (Sigma).

PVY Genomic RNA Isolation and Sequencing
----------------------------------------

Total RNAs were isolated from homogenates of upper symptomatic leaves using TRI-reagent kit (Sigma) and reverse-transcribed into cDNA using ReverTAID reverse-transcriptase (Thermoscientific). Full-length cDNA encoding PVY genomes and VPg cistrons were further PCR-amplified using primers flanking either the whole PVY genome or VPg cistrons. PCR products were sent to GENEWIZ Company or Genoscreen (Lille, France) for Sanger sequencing. All primer sequences are listed in [supplementary table S3](#sup1){ref-type="supplementary-material"}, [Supplementary Material](#sup1){ref-type="supplementary-material"} online.

Shannon Diversity Analysis
--------------------------

To compare the diversity of *pvr2^3^*-resistance-breaking variants between SON41p^ID2^ and SON41^WT^, Shannon diversity and corresponding 95% confidence intervals were calculated using an optimal nonempirical method described in [@msx249-B9] using the Species-richness Prediction And Diversity Estimation in R \[SpadeR\]) online tool ([@msx249-B10]).

Determination of VPg/eIF4E In Vitro Interaction Affinity by Fluorescence Assays
-------------------------------------------------------------------------------

All steady-state fluorescence acquisitions were obtained at 25 °C in 20 mM HEPES pH 7.5, 0.25 mM NaCl and 1 mM DTT using a SAFAS Xenius spectrofluorimeter (Monaco) equipped with a Peltier temperature controller. For optical characteristics of the instrument, see [@msx249-B32]. Affinity constants (1/*K*~D~) were deduced from steady-state eIF4E tryptophan intrinsic fluorescence decrease upon titration by VPg as previously described in [@msx249-B32]. Measurements were performed by titration of a 0.25 × 10^−6 ^M eIF4E solution. The VPg-eIF4E complex formation was titrated by 18--20 successive 2 µl additions of a 5 × 10^−6 ^M VPg solution. After each VPg addition, eIF4E fluorescence emission was recorded for 2 min at 340 nm (excitation, 280 nm). The affinity constants of each of the VPg forms with eIF4E from Yolo Wonder and HD 285 were assessed by performing three simultaneous technical repetitions. Each of them included an interaction test with a m7-GTP cap analog as a positive control for eIF4E integrity and activity. Signal treatments and *K*~D~ determinations were performed by nonlinear regression following the method described in [@msx249-B32]. All measurements were performed twice using two different protein preparations for each of VPg and eIF4E.

Supplementary Material
======================

[Supplementary data](#sup1){ref-type="supplementary-material"} are available at *Molecular Biology and Evolution* online.

Supplementary Material
======================

###### 

Click here for additional data file.

We would like to thank Coralie Chesseron, Grégory Girardot, Ghislaine Nemouchi, Alain Palloix, Vincent Simon, and Marion Szadkowski for their contribution to the *in planta* experiments. We are also grateful to Sonia Longhi and Thierry Candresse for enriching discussions that contributed to improve this work.

This work was partially supported by Le Ministère Français de l'enseignement supérieur et de la Recherche (J.C. Fellowship).

[^1]: **Associate editor:** Juliette de Meaux

[^2]: Second and third line/column: polymorphic amino acid sites and their corresponding positions in VPg protein sequence. Red/blue boxes: amino acids substitutions that are predicted (PONDR-VLXT) to enhance or reduce intrinsic disorder content in VPg-SON41p respectively. "d" and "o": respectively intrinsic disorder and order promoting substitutions; "*RB"*: *pvr2^3^*-resistance breaking sites previously identified ([@msx249-B4]; [@msx249-B33]); "*RB\**": expected resistance-breaking sites (displaying high d*N*/d*S* ratios; [@msx249-B3]; [@msx249-B34]).
