# Characterization of Grapevine Fanleaf Virus Isolates in 'Chardonnay' Vines Exhibiting Severe and Mild Symptoms in Two Vineyards

Julie Kubina, Jean-Michel Hily, Pierre Mustin, Véronique Komar, Shahinez Garcia, Isabelle Rachel Martin, Nils Poulicard, Amandine Velt, Véronique Bonnet, Laurence Mercier, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Julie Kubina, Jean-Michel Hily, Pierre Mustin, Véronique Komar, Shahinez Garcia, et al. Char- acterization of Grapevine Fanleaf Virus Isolates in ‘Chardonnay’ Vines Exhibiting Severe and Mild Symptoms in Two Vineyards. Viruses, 2022, 14 (10), pp.2303. 10.3390/v14102303. hal-03864820

HAL Id: hal-03864820

https://hal.science/hal-03864820

Submitted on 22 Nov 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of sci- entific research documents, whether they are pub- lished or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Characterization of Grapevine Fanleaf Virus Isolates in 'Chardonnay' Vines Exhibiting Severe and Mild Symptoms in Two Vineyards

Julie Kubina<fnref n="1" />, Jean-Michel Hily <fnref n="1" />,<fnref n="2" />, Pierre Mustin<fnref n="1" />, Véronique Komar<fnref n="1" />, Shahinez Garcia<fnref n="1" />, Isabelle Rachel Martin <fnref n="3" />, Nils Poulicard <fnref n="4" />, Amandine Velt <fnref n="1" />, Véronique Bonnet <fnref n="5" />, Laurence Mercier <fnref n="5" />, Olivier Lemaire <fnref n="1" /> and Emmanuelle Vigne <fnref n="1" />,*

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Kubina, J.; Hily, J.-M.; Mustin, P.; Komar, V.; Garcia, S.; Martin, I.R.; Poulicard, N.; Velt, A.; Bonnet, V.; Mercier, L.; et al. Characterization of Grapevine Fanleaf Virus Isolates in 'Chardonnay' Vines Exhibiting Severe and Mild Symptoms in Two Vineyards. *Viruses* **2022**, *14*, 2303. https://doi.org/10.3390/v14102303

Academic Editors: Luisa Rubino and
K. Andrew White

Received: 1 September 2022

Accepted: 18 October 2022

Published: 20 October 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Fanleaf degeneration is a complex viral disease of *Vitis* spp. that detrimentally impacts fruit yield and reduces the productive lifespan of most vineyards worldwide. In France, its main causal agent is grapevine fanleaf virus (GFLV). In the past, field experiments were conducted to explore cross-protection as a management strategy of fanleaf degeneration, but results were unsatisfactory because the mild virus strain negatively impacted fruit yield. In order to select new mild GFLV isolates, we examined two old 'Chardonnay' parcels harbouring vines with distinct phenotypes. Symptoms and agronomic performances were monitored over the four-year study on 21 individual vines that were classified into three categories: asymptomatic GFLV-free vines, GFLV-infected vines severely diseased and GFLV-infected vines displaying mild symptoms. The complete coding genomic sequences of GFLV isolates in infected vines was determined by high-throughput sequencing. Most grapevines were infected with multiple genetically divergent variants. While no specific molecular features were apparent for GFLV isolates from vines displaying mild symptoms, a genetic differentiation of GFLV populations depending on the vineyard parcel was observed. The mild symptomatic grapevines identified during this study were established in a greenhouse to recover GFLV variants of potential interest for cross-protection studies.

**Keywords:** grapevine fanleaf virus; fanleaf degeneration; grapevine; symptomatology; mild isolates; cross-protection; genetic and phenotypic diversity; high throughput sequencing; virome

## 1. Introduction

Grapevine is a crop hosting many viruses (more than 90 viruses and viroids have been identified in grapevines) with some of them affecting its cultivation and fruit production worldwide [1–3]. Among grapevine viral diseases, fanleaf degeneration is the most destructive [4]. In France, this disease is estimated to be present over more than 60% of the vineyard hectarage [5]. The main virus responsible for fanleaf degeneration is grapevine fanleaf virus (GFLV), a soil-borne nepovirus of the *Secoviridae* family, specifically transmitted by the ectoparasitic dagger nematode *Xiphinema index* [6]. GFLV infections can cause massive crop losses (up to 80%) and a shorter productive lifespan of vineyards due to a progressive degeneration that can lead to vine death [4,7–9].

Typical symptoms associated with fanleaf degeneration are foliar discolouration (i.e., complete yellowing, vein banding, mosaic), foliar deformations (i.e., small leaf, leaf with open petiole sinus), shoot abnormalities (i.e., short internodes, fasciation) and/or growth inhibition of the plant (i.e., stunting) [6]. These symptoms can vary in intensity and range

<footnote n="1">1 INRAE, SVQV UMR-A 1131, Université de Strasbourg, 68000 Colmar, France</footnote>

<footnote n="2">2 IFV, 30240 Le Grau-Du-Roi, France</footnote>

<footnote n="3">3 LVBE, UPR 3991, Université de Haute-Alsace, 68000 Colmar, France</footnote>

<footnote n="4">4 PHIM, Université Montpellier, IRD, INRAE, Cirad, SupAgro, 34000 Montpellier, France</footnote>

<footnote n="5">5 Maison Moët & Chandon, 20 Avenue de Champagne, 51200 Épernay, France</footnote>

<footnote>* Correspondence: emmanuelle.vigne@inrae.fr</footnote>
from mild to severe, depending on environmental factors, rootstock genotype, *Vitis vinifera* cultivar (cv.), viral strain, vineyard site, and viticultural practices [10–13]. While viral determinants of specific disease symptoms have been identified in herbaceous model plants [14–16], they have yet to be determined in the main natural host of GFLV, grapevine.

The genome of GFLV is composed of two positive-sense RNAs named RNA1 and RNA2 of 7.3 and 3.8 kb in size, respectively [6]. Both RNAs are necessary for infection in planta, and both contain a single open reading frame (ORF1 and ORF2) that codes for a polyprotein (P1 and P2) processed by the viral proteinase into functional proteins. A recent genetic diversity study based on a cohort of 40 GFLV ORF1 and 80 ORF2 sequences from around the world described a very high degree of polymorphism for both ORFs with 0.127 and 0.130 nucleotide substitutions per site, respectively [17]. Many recombination events were identified all along the GFLV genome, potentially explained by the presence of many variants within a single plant [17–20]. An additional satellite RNA (RNA3) of 1.1 kb is associated with some GFLV isolates [21]. This third RNA contains a unique ORF (ORF3) coding for a non-structural protein of unknown function and does not seem to interfere with viral pathogenicity [22].

Current strategies to control fanleaf degeneration rely on prophylactic measures based on the production of certified materials (i.e., free of GFLV and other damaging viruses) which are effective in avoiding the introduction of the virus in uninfected vineyards [23]. However, once the fanleaf pathosystem (viruliferous nematodes) is established in a vineyard it is almost impossible to eradicate [24]. The removal of infected plants followed by a five to seven years fallow period decreases the nematode vector population, but this option is economically unpractical for grape growers [23]. Rootstocks tolerant to *X. index* are available but they seem to only delay the infection by GFLV [25–27]. The use of chemicals to control nematode vector populations has been banned in Europe due to environmental and human health concerns [28,29]; the use of plants exhibiting negative impact on nematodes are currently under study [30] with particular interest of the *Fabaceae* family exhibiting nematicidal activity [31]. Biotechnological approaches to engineer resistance against GFLV are promising but are unlikely to be adopted in Europe [32]. A recessive factor of resistance against GFLV was recently described in *V. vinifera* cv. Riesling, providing interesting perspectives for fanleaf management [33]. However, at least two decades of research will be needed for characterizing this resistance factor, its sustainability, and introducing this genetic trait in grapevine genotypes prior to the deployment of new resistant hybrids in vineyards. For all these reasons, other strategies to mitigate the effects of fanleaf degeneration are needed, particularly in old high-value vineyards where growing vines is very challenging due to the severity of the disease.

In this context, cross-protection offers an interesting biocontrol alternative to the aforementioned strategies for the management of fanleaf degeneration [34]. Cross-protection is a natural phenomenon, first described by McKinney almost a century ago [35], in which a primary infection with a mild virus strain protects a plant from the disease induced by related variants of the same virus [36]. For viral disease management in crops, cross-protection relies on the selection of virus strains that are able to (1) cause mild symptomatic infections while (2) conferring protection against the effects of severe variants [34,37]. In other words, a plant can be deliberately pre-immunized with a mild viral variant to protect it from aggressive counterparts of genetically related variants that will challenge it in the field. Many efforts have been made to explain the protection in primarily infected plants; however, the molecular mechanism(s) underlying cross-protection remain(s) largely unclear [34,38].

Cross-protection is used to control some viruses in different crops of economic importance such as pepino mosaic virus (PepMV) in tomato [39,40] or citrus tristeza virus (CTV) in citrus [41,42]. Mild viral strains are obtained either by laboratory methods, including single local lesion isolations on plant hosts, temperature treatments, physical or chemical treatments, site-directed mutagenesis, or via field selections [34,43–48]. In the latter case, mild viral strains were selected in severely diseased fields in which infected plants showed little to no symptom [34,44]. The main advantages of such naturally occurring mild iso-
lates are that: (1) they are well adapted to the cultivar and its environmental conditions and (2) they are expected to ensure protection against local (or endemic) severe variants. A prime example of these mild virus isolate features is the CTV cross-protective mild isolate IAC obtained from sweet orange cv. Pera in Brazil (PIAC isolate). Millions of protected 'Pera' trees established in sweet orange orchards show a durable protection against stem pitting strains of CTV, leading to the success of the Brazilian citriculture [42]. Yet, this mild isolate was not efficient on cvs. Hamlin or Valencia, and not protective in all Brazilian regions. It appears that a fundamental requirement to achieve successful cross-protection, at least for this viral species, might be the selection of mild isolates for each cultivar associated with location/locality of interest [44,49,50].

The evaluation of cross-protection to mitigate the impact of GFLV was initiated over thirty years ago in our laboratory with one mild GFLV strain, named GHu, isolated from *V. vinifera* cv. Gloria Hungariae. This strain, which causes mild symptoms in *Chenopodium quinoa*, could prevent under certain inoculation conditions, the development of the severe symptoms induced by the GFLV-F13 strain [51]. Primary-infected vines (cv. Gewurztraminer) with GFLV-GHu exhibited a significantly reduced superinfection rate in a GFLV-infected vineyard site located in the Alsace region of France [52]. However, these cross-protected vines showed reduction of fruit yield over the eight-year study in comparison with non-cross-protected vines, making cross-protection based on GHu strain of limited interest under these conditions [52].

Thus, to pursue cross-protection as an effective biocontrol management strategy against fanleaf degeneration, it is fundamental to identify new mild GFLV strains in diseased vineyards, especially in the Champagne region of France where fanleaf degeneration is one of the major concerns of wine grape growers [53]. The main objective of this study was to identify GFLV-infected vines showing mild symptoms in two heavily diseased vineyard parcels and to characterize the diversity of GFLV isolates in infected 'Chardonnay' vines by high-throughput sequencing (HTS). By recording qualitative and quantitative phenotypic traits during four consecutive years, four mildly symptomatic vines infected with potentially mild GFLV isolates were identified.

## 2. Materials and Methods

### 2.1. Vineyard Sites

Two commercial GFLV-infected vineyard sites, one in Chouilly (named 'Pa') and the other in Cramant (named 'Py'), two villages spaced 5.4 km apart in the Champagne region of France, were selected for this study. These two vineyard parcels were planted in the 1980s after treatment with a systemic herbicide and soil disinfection using nematicidal chemicals (Shell DD and Temik10G). In line with viticultural practices in the Côte des blancs area in the Champagne region, vines consisted of *V. vinifera* cv. Chardonnay grafted onto a rootstock, likely 41B (*V. berlandieri* × *V. riparia*), as indicated by wine growers of both sites.

### 2.2. Disease Symptom Monitoring, Grapevine Leaf Sampling, and GFLV Detection

From 2016 to 2019, qualitative and quantitative phenotypic traits were recorded on selected vines in the two above mentioned 'Chardonnay' vineyard parcels (Figure S1). Qualitative phenotypic traits belonging to three kinds of symptoms (discolouration of leaves, deformation of leaves and stunting of plants) were evaluated and recorded every year in June, the best time to visualize GFLV symptoms in the northern hemisphere. The same symptom scoring as described in our recent article was used [13]: discolouration and deformations of leaves were estimated on individual vines and scored from 0 to 4 (0: absence of symptom, 1: 1 to 25%, 2: 26–50%, 3: 51–75%, and 4: 76–100% of symptomatic leaves) (Figure S2). Similarly, vine stunting was assessed and scored from 0 to 4 (0: no stunting, 1: low, 2: medium, 3: strong, 4: very strong stunting). An overall symptom severity score was determined by summing all qualitative trait scores measured for each
plant. Quantitative phenotypic traits were also assessed on individual vines by counting
the number of grape clusters at harvest and by measuring pruning wood weights in winter.

Five to six young apical leaves from different shoots were collected in the spring for every plant from 2016 to 2019 and stored frozen prior to analyses [19]. The presence of GFLV was evaluated by double antibody sandwich-enzyme linked immunosorbent assay (DAS-ELISA) using GFLV specific antibodies (BIOREBA AG, Reinach, Switzerland), and by molecular techniques as previously described [19].

## 2.3. Statistical Analyses of Qualitative and Quantitative Phenotypic Traits

The significances of the differences between symptom categories for all traits were
tested by comparison of means (Student's *t* test or Welch's *t* test) or medians (Kruskal–
Wallis *H* test) and completed by Tukey's honestly significant difference (HSD) tests when
needed. The relationship between the overall symptom scores and fruit yields was analysed
using Spearman's rank correlation coefficient. Principle component analysis (PCA) was
performed to evaluate associations among individuals and among the variables using the
FactoMineR package [54]. Differences were considered significant when *p* < 0.05. Boxplots
were drawn to show the median (horizontal bold line) and interquartile range with lonely
dots representing extreme data. All data and graphic representations were performed using
R software (version 4.0.2; R Core Team, 2021, https://www.r-project.org/ (accessed on
7 April 2021)). Raw data are available in Table S1.

## 2.4. RNA Sequencing and Downstream Bioinformatics Analyses

Total RNA extractions (from grapevine leaves collected in 2017), cDNA library prepa-
rations, Illumina sequencing and HTS data analyses (for the determination of the viral
status) of the twenty vines selected for this study are described in our previous article [19].
Raw Illumina sequence data files are available on ENA-EBI, project N°PRJEB54883. For
vine Py43, total RNA extracted from leaves collected in 2018 was sequenced through
a small RNA-Seq run on a Hiseq Instrument at 1 × 50 pb and multiplexed at Fasteris
(Plan-les-Ouates, Switzerland).

## 2.5. GFLV Sequences Analyses, Genetic Diversity and Recombination Detection

Each GFLV consensus sequence generated by de novo assembly [55] was named
as follows: first by indicating the geographic origin of the sample [Pa] or [Py] and the
identification number of the vine, second by designating the genomic or satellite GFLV RNA
under consideration [1,2] or [3] and third by classifying the molecular variant depending
on its accumulation based on the number of reads per kilobase per million of total reads
(RPKM), the first being the one with the highest value.

GFLV consensus sequences were aligned using CLUSTALW (codon-based multiple
alignments) with default parameters in MEGA7 software [56]. Maximum likelihood (ML)
phylogenetic trees were reconstructed using the same software. The best-fit ML model for
each sequence alignment was used (GTR + G + I, GTR + G, T92 + G for the complete ORF1,
ORF2, ORF3 datasets, respectively) and bootstrapping analyses of 100 replicates were
performed. BioEdit software [57] was employed to construct nucleotide identity matrices
for clade definition. Trees were visualized and fashioned in FigTree v1.3.1. (Institute
of Evolutionary Biology, University of Edinburgh, Edinburgh. http://tree.bio.ed.ac.uk/
software/figtree/ (accessed on 27 August 2019)).

Genetic diversity (π) values were calculated in MEGA7 software using the best-fit ML
model determined for each dataset (TN93 + G for ORF1, ORF2, 1BHel and 1EPol; K2 + G
for 1A, 1CVPg, 2A and 2BMP; T92 + G for 1DPro and 2CP) and bootstrapping analyses
of 100 replicates were performed. π and Tajima's D (D~T~) sliding window analyses were
conducted in DnaSP v6.12.03 [58] with window length of 100 nts and step size of 25 nts
parameters. The difference between non-synonymous (dN) and synonymous (dS) substitu-
tions over the coding sequences of GFLV populations was estimated by the Kumar method
in MEGA7 software. GFLV populations were compared by calculating the fixation index
$(F_{ST})$ in Arlequin v3.5.2.2 [59] using the above-mentioned models. Significance of the differences was obtained by performing 1000 permutations. Finally, potential recombination events were estimated by using RDP4 and the default settings of the seven algorithms: RDP, GENECONV, BootScan, MaxChi, Chimaera, SiScan and 3Seq [60]. Only recombination events detected by five or more methods were considered.

# 3. Results

## 3.1. Identification of 'Chardonnay' Vines Showing Distinct Fanleaf Symptoms

To select potentially protective GFLV isolates suitable for cross-protection of 'Chardonnay' vines in the Champagne region of France, our strategy was to prospect old vineyards highly impacted by fanleaf degeneration disease and identify vines displaying mild disease symptoms. Two 30-year-old 'Chardonnay' parcels, 'Py' and 'Pa', with an estimated vine mortality of 15 to 20% very likely due to GFLV infections, were selected for this study in the 'Côte des Blancs' area. Despite soil disinfection in these parcels in the 1980s, propagation of fanleaf degeneration was inexorable with increasing symptomatic patches (or areas) being observed yearly. The presence of the nematode vector *X. index* was confirmed in both parcels after soil analysis (data not shown).

In spring 2016, we surveyed both sites for vigorous vines with a good potential for fruit production and proximal to severely diseased vines. An overall symptom score achieved by the addition of the scores for three qualitative phenotypic traits, i.e., leaf discolouration, leaf deformations and stunting of the plants, was used to classify the vines (Figure S2). Mildly symptomatic and asymptomatic vines were considered of potential interest for our study if they received a maximum of '3' for the overall symptom score. Five mild ('M') symptomatic and five severe ('S') symptomatic vines were selected in each vineyard site and all these vines were around 30 years old (Table S1). In summer and winter 2016, quantitative traits, i.e., cluster number, fruit yield, fruit yield per cluster, cane number, pruning wood weight per vine and pruning wood weight per cane were assessed for the ten selected vines at each vineyard site. An additional vine named Py43 was incorporated in the M class in September 2016 (Figure S1, Table S1). Therefore, our plant population under study consisted of 21 vines, including 11 and 10 vines grouped into the M and S classes, respectively. A comparative analysis of qualitative phenotypic trait scores and quantitative data for each class revealed significant differences for all the comparisons, except for leaf deformation and cane number (Figure S3, Table S2).

The presence of GFLV in leaves of the selected vines was tested by DAS-ELISA. While all selected grapevines were located in severely diseased areas, the virus was surprisingly only detected in less than half of the grapevines displaying mild symptoms (Table 1). These results were confirmed by using several diagnostic molecular techniques, as previously reported [19]. However, the virus was detected in all the severely symptomatic vines. Consequently, the cohort of the selected grapevines was divided in three categories (Tables 1 and S1): 'M-' (mild symptomatic GFLV-free vines), 'M+' (mild symptomatic GFLV-infected vines) and 'S+' (severe symptomatic GFLV-infected vines). These three categories were used to examine the candidates for cross-protection by comparing them to either non-infected vines or severely infected ones.

Comparisons of phenotypic traits showed that scores for vines in the M+ category ranged generally between scores for vines in the M− and S+ categories (Figure 1). There are no significant differences in any traits between M− and M+ categories, unlike M− and S+ which were different for most traits (Figure 1, Table S3). Thereby four mild symptomatic GFLV-infected 'Chardonnay' vines (Pa1, Py18, Py19 and Py43) were identified in 2016 and were considered as candidates hosting potentially mild isolates of interest for cross-protection.
<table><caption>**Table 1.** Presence (+) or absence (−) of grapevine fanleaf virus (GFLV) in leaves of individual 'Chardonnay' vines collected in spring from 2016 to 2019, as determined by DAS-ELISA. Vines were classified into three categories: M− (mild symptomatic, GFLV-free vines), M+ (mild symptomatic, GFLV-infected vines) and S+ (severe symptomatic, GFLV-infected vines). The presence/absence of GFLV in vine Py43 was only tested starting in 2017 as this plant was only identified in autumn 2016 (nd: not determined).</caption><thead><tr><th rowspan="2"></th><th colspan="8">M-</th><th colspan="8">M+</th><th colspan="8">S+</th></tr><tr><th>Pa3</th><th>Pa5</th><th>Pa7</th><th>Pa9</th><th>Py11</th><th>Py13</th><th>Py16</th><th>Pa1</th><th>Py18</th><th>Py19</th><th>Py43</th><th>Pa2</th><th>Pa4</th><th>Pa6</th><th>Pa8</th><th>Pa10</th><th>Py12</th><th>Py14</th><th>Py15</th><th>Py17</th><th>Py20</th></tr></thead><tbody><tr><td>2016</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>+</td><td>+</td><td>+</td><td>nd</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>2017</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>2018</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>2019</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Comparative analyses of qualitative and quantitative traits for the three vine categories in 2016 (M−, mild symptomatic GFLV-free vines (in grey); M+, mild symptomatic vines infected with GFLV (in cyan) and S+, severe symptomatic vines infected with GFLV (in black)). Number of vines (n) are given in brackets. * For the vines in the M+ category, only three vines were considered in June 2016 for the evaluation of qualitative traits (leaf discolouration, leaf deformation, vine stunting, overall symptoms). Significance was tested with Student's *t*, Welch's *t* or Kruskal-Wallis *H* tests. Boxplots show the median (horizontal bold line) and the interquartile range with lonely dots representing extreme data. Different letters (a,b) indicate significant differences (*p* ≤ 0.05). Raw data and results of the statistical tests are available in Tables S1 and S3.</figcaption></figure>

## 3.2. *Phenotypic Differentiation and Overtime Stability of M−, M+ and S+ Categories*

To assess the phenotypic durability of the three vines categories defined in 2016, qualitative and quantitative traits were monitored from 2017 to 2019 and the presence of GFLV was evaluated by DAS-ELISA (Tables 1 and S1). GFLV was always detected in leaves from M+ and S+ vines, while remaining undetectable by DAS-ELISA in leaves from vines in the M− category (Table 1).

Mean symptom scores and mean yield values from 2016 to 2019 for each individual vine were assessed by PCA and hierarchical clustering on principle components. The representation accuracy is supported by the two first principal components explaining 81.22% of the individuals and variables cloud variances (Figure 2).
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Principal components analysis (PCA) and hierarchical clustering of 21 vines according to their mean symptom scores and yield values over the four-year study (from 2016 to 2019). (A) Biplot representing the ‘M−’ (grey), ‘M+’ (cyan) and ‘S+’ (black) individuals. Projection of M−, M+ and S+ categories are represented by squares with their confidence ellipses (95%). (B) Dendrogram representing the hierarchical clustering based on the individuals coordinates on the 10 principal components. (C) Loading biplot, the vectors show the relative importance of each variable in discriminating amongst the observations. Longer vectors indicate greater contribution. Vectors point in the direction they ordinate the observations (DIS: Leaf discolouration, DEF: leaf deformation, STU: vine stunting, ALL: overall, CLU: cluster number, FY: fruit yield, FYpc: fruit yield per cluster, CAN: cane number, PWW: pruning wood weight per vine, PWWpc: pruning wood weight per cane). Results of the Fisher’s test are available in Table S4.</figcaption></figure>

The individuals were clearly separated along the first principal component (PC1) according to their categories (M−, M+ and S+). Fisher’s tests, assessing correlations between both symptom and yield variables, and PC1, showed respectively significant negative and positive correlation coefficients (Table S4). Thus, PC1 ordinates the S+ vines with severe symptoms and low yields in one direction (to the left), and the M− vines, showing
mild symptoms and higher yields in the opposite direction (to the right). M+ vines were dispersed in the middle of the two groupings (Figure 2). Furthermore, the dendrogram resulting from hierarchical clustering on individuals coordinates on all principal components (i.e., explaining 100% of the variability) confirmed this trend with three major clusters with an almost perfect segregation of M-, M+ and S+ vines (Figure 2). It should be noted that overall symptom score and fruit yield were the main indicators of the variability with the longer vector values, indicating their greater contribution in the separation of the data (Figure 2). Consequently, results obtained with these two variables were further evaluated in our study, while the other phenotypic traits are only detailed in Supplementary Data.

Regarding phenotypic traits evolution from 2016 to 2019, the annual mean values of M+ vines ranged each year between values obtained for vines in the M- and S+ categories, illustrating a steady pattern. This trend was observed for overall symptoms, fruit yield (Figure 3) and other traits (Figure S4). Therefore, each category of vines remains constant along the study.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Evolution of yearly mean overall symptom scores and fruit yields of vines in the M-, M+ and S+ categories and their standard error. Results for vines in the M- (in grey), M+ (in cyan) and S+ (in black) categories are shown. Number of vines (n) is given in brackets. Means values are connected by bold lines and standard errors are materialized by vertical bars. * For the M+ category, only three vines were considered in June 2016 for the evaluation of the qualitative traits.</figcaption></figure>

Across the four-year study, statistical differentiation of the three categories of vines was observed for the two main indicators (Figure 4, Table S3). In particular, vines in the M+ category showed significantly lower overall symptom scores and higher fruit yields in comparison to vines in the S+ category. Similar trends were observed with the other variables, even though differences were not always statistically validated (Figure S5, Table S3). For fruit production, vines in the M+ category (mean 1795 g/vine with 95% confidence interval at 1305 to 2284) produced about 40% less relatively to vines in the M- category (mean 2971 g/vine with 95% confidence interval at 2471 to 3471). However, the same M+ category was 82% more productive than the S+ category (mean 988 g/vine with 95% confidence interval between 697 and 1279). The promising results of the M+ category underline the beneficial agronomic performances provided by mildly symptomatic GFLV-infected vines in comparison to severely diseased ones, thus confirming their potential use in a cross-protection program.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Comparison of overall symptom scores and fruit yields over four years for vines in M-, M+ and S+ categories. Results for vines in the M- (in grey), M+ (in cyan) and S+ (in black) categories are shown. Boxplots show the median (horizontal bold line) and the interquartile range with lonely dots representing extreme data. Different letters (a,b,c) indicate significant differences (p ≤ 0.05). Raw data and results of the statistical tests are available in Tables S1 and S3, respectively. Number of vines (n) is given in brackets. * For the M+ category, 15 replicates were considered for overall symptom scores.</figcaption></figure>

### 3.3. Virome Description and Genetic Characterization of GFLV Isolates: Mixed Infections of Highly Divergent Molecular Variants

To characterize the molecular variants in the GFLV-infected vines, HTS was performed using total RNA isolated from leaves collected on the 21 studied grapevines. Besides GFLV, the usual grapevine background virome (i.e., GRSPaV, grapevine rupestris stem pitting-associated virus, HSVd, hop stunt viroid, and GYSVd-1, grapevine yellow speckle viroid 1) was found in all vines (Table S5). Other viruses of the *Tymoviridae* family (GFkV, grapevine fleck virus, GSyV1, grapevine Syrah virus 1, GRVFV, grapevine rupestris vein feathering virus, and GRGV, grapevine red globe virus) were detected in some plants. Grapevine Pinot gris virus (GPGV) was identified in a single vine and grapevine leafroll-associated virus 2 (GLRaV-2,) was detected in four of them (Table S5). Regarding the accumulation of these viruses, viral normalized read counts was low (1 to 219 RPKM) in comparison to GFLV genomic (362 to 3706 RPKM) and satellite (1145 to 6809 RPKM) RNAs (Table S5).

Consensus sequences of GFLV RNAs that were assembled through HTS analyses were considered as distinct molecules and referred to as molecular variants. Partial 5' untranslated regions (UTR) and 3'UTR were recovered for these molecular variants, but de novo assembly was challenging in these regions due to highly conserved stem loops shared between GFLV RNAs. Therefore, due to this limitation, we focused subsequent analyses on the GFLV ORF sequences only. In total, 27 GFLV ORF1, 19 GFLV ORF2 and five GFLV satellite ORF3 molecular variants (GenBank ON991744-ON991794) were obtained after de novo sequence assembly (Table S5).

The length of ORF1 (6852 nts), ORF2 (3330 nts), and ORF3 (1023 nts) was identical for all the molecular variants characterized in this study, except for Py17-3-1 ORF3 molecular variant (1014 nts). Analyses of the nucleotide and deduced protein sequences indicated identities of 86.6–99.8% in nts (91.9–99.8% in amino acids, aa) for ORF1, 87.8–99.6% in nts (94.2–99.5% in aa) for ORF2 and 82.4–99.7% in nts (77.8–99.4% in aa) for ORF3 (Table S6).

Putative recombination events were assessed using RDP4. Fourteen and four recombination events were predicted within GFLV ORF1 and ORF2, respectively, but none in GFLV
ORF3 (Figure S6). Recombination breakpoints were identified within all viral genomic coding regions, except for the ones coding for the helicase (1B<sup>Hel</sup>) and the VPg (1C<sup>VPg</sup>).

A maximum likelihood (ML) tree based on GFLV nucleic acid sequences alignment was constructed for each ORF. Phylogenetic analyses revealed clustering of sequences into twelve, six, and two distinct clades for ORF1, ORF2 and ORF3, respectively (Figure 5). The criteria used to define clade affiliation was at least 95% sequence identity at the nucleotide level (Table S6). Looking at the molecular variant composition in each vine, most of the vines showed mixed GFLV infections with the detection of two divergent variants (i.e., from different clades) of ORF1 in ten vines and of ORF2 in five vines (Table 2). In addition to the numerous variants of each ORFs present within a vine, many different combinations of GFLV ORF1, ORF2 and ORF3 variants belonging to distinct clades were observed. Consequently, every infected vine was unique regarding its composition of GFLV variants, with the exception of two neighbouring vines Py12 and Py15, for which GFLV variants from the same clades (i.e., I.A, I.D, II.B and II.C) were detected (Table 2).

<figure><img src="image_9.png" /></figure>

**Figure 5.** Phylogenetic relationships between grapevine fanleaf virus (GFLV) molecular variants. Maximum likelihood trees were inferred from the complete nucleotide sequences of ORF1 (A), ORF2 (B) and ORF3 (C) recovered from GFLV-infected 'Chardonnay' vines. Each sequence name indicates the vineyard site (Pa, in orange and, Py, in deep blue), the sampled vine, the molecular variant and the vine affiliation to a phenotypic category (M+, in cyan and S+, in black) in parenthesis. Clades are named with the Roman numeral corresponding to the GFLV ORF and classified from the most [A] to the least [L] represented clade in sequence number. The names of the distinct clades are highlighted in bold. Sequence variants sharing at least 95% nucleotide identity were considered in the same clade; pairwise nucleotide identities values between GFLV sequences are available in Table S6. Scale bars below each tree show genetic distance. Only bootstrap values ≥ 0.95 are indicated.
<table><caption>**Table 2.** Molecular variants composition of grapevine fanleaf virus in individual 'Chardonnay' vines and their classification in phylogenetic clades (Figure 5). The clade names are 'I' for ORF1, 'II' for ORF2 and 'III' for ORF3. The number of molecular variants is indicated. When two molecular variants affiliated to a same clade are present in a vine, their numbers are separated by a slash.</caption><thead><tr><th rowspan="2"></th><th colspan="4">M+</th><th colspan="10">S+</th></tr><tr><th>Pa1</th><th>Py18</th><th>Py19</th><th>Py43</th><th>Pa2</th><th>Pa4</th><th>Pa6</th><th>Pa8</th><th>Pa10</th><th>Py12</th><th>Py14</th><th>Py15</th><th>Py17</th><th>Py20</th></tr></thead><tbody><tr><td>I.A</td><td></td><td></td><td></td><td>1-1</td><td></td><td></td><td></td><td></td><td></td><td>1-1/1-2</td><td>1-1/1-2</td><td>1-1</td><td></td><td>1-1</td></tr><tr><td>I.B</td><td></td><td>1-1</td><td>1-1</td><td></td><td></td><td></td><td></td><td></td><td>1-2</td><td></td><td>1-3</td><td></td><td>1-1</td><td></td></tr><tr><td>I.C</td><td></td><td></td><td></td><td></td><td>1-1</td><td>1-1</td><td>1-1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>I.D</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1-3</td><td></td><td>1-2</td><td></td><td></td></tr><tr><td>I.E</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1-1</td><td>1-1</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>I.F</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1-2/1-3</td></tr><tr><td>I.G</td><td></td><td></td><td></td><td>1-2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>I.H</td><td></td><td>1-2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>I.I</td><td></td><td></td><td>1-2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>I.J</td><td></td><td></td><td></td><td></td><td></td><td></td><td>1-2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>I.K</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1-2</td><td></td></tr><tr><td>I.L</td><td>1-1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>II.A</td><td></td><td></td><td>2-1</td><td></td><td></td><td>2-1</td><td></td><td></td><td></td><td></td><td></td><td>2-2</td><td>2-2</td><td>2-1</td></tr><tr><td>II.B</td><td>2-1</td><td>2-1</td><td></td><td></td><td></td><td></td><td></td><td></td><td>2-2</td><td>2-2</td><td>2-1</td><td>2-2</td><td>2-1</td><td></td></tr><tr><td>II.C</td><td></td><td></td><td></td><td>2-1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>II.D</td><td></td><td></td><td></td><td></td><td></td><td></td><td>2-1</td><td>2-1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>II.E</td><td></td><td></td><td></td><td>2-1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>II.F</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>2-1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>III.A</td><td></td><td></td><td></td><td>3-1/3-2</td><td>3-1</td><td>3-1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>III.B</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>3-1</td><td></td></tr></tbody></table>

### 3.4. Comparison of GFLV Molecular Variants in the M+ and S+ Categories and in Pa and Py Vineyard Sites

Most phylogenetic clades harboured molecular variants from vines of both M+ and S+ categories. On the other hand, molecular variants from the 'Pa' and 'Py' vineyards were distributed in different clades with the exception of clades I.B, II.A and II.B (Figure 5). These observations were supported following GFLV populations differentiation and fixation index ($F_{ST}$) analysis. Comparison using molecular variant populations in vines in the M+ and S+ categories displayed very low $F_{ST}$ and non-significant *p*-values (*p* > 0.05), indicating no genetic differentiation between the two populations either when ORF1 ($F_{ST(S+/M+)}$ = 0.019, *p* = 0.246), ORF2 ($F_{ST(S+/M+)}$ = −0.060, *p* = 0.701) or any coding regions of either ORFs being considered (Figure 6, Table S7). In contrast, a significant genetic differentiation between molecular variant populations of the 'Pa' and 'Py' vineyards was statistically validated for both ORF1 and ORF2 sequences ($F_{ST(Pa/Py)}$ = 0.104, *p* = 0.009 for ORF1; $F_{ST(Pa/Py)}$ = 0.145, *p* = 0.016 for ORF2), as well as for most coding regions of either ORFs (Figure 6, Table S7). Regarding ORF3 molecular variants, they were only found in four vines out of ten in the S+ category.

Similar nucleotide diversity ($\pi$) patterns were found for all molecular variant populations either on ORF1, ORF2 or specific coding regions of these ORFs sequences (Figure 6, Table S7). Regardless of the sequence populations, molecular variants of ORF1 were slightly more diverse than those of ORF2 with overall $\pi$ values calculated at 0.197 ± 0.007 for ORF1 versus 0.155 ± 0.009 for ORF2. The most divergent region was the 3' extremity of ORF1 (Figure 6). When looking at the constraints ($dN - dS$) acting on ORFs and coding regions, an excess of synonymous versus nonsynonymous substitutions was observed, indicating that these sequences were mainly under negative selection, with the exception of positive values observed for the coding region of $1C^{VPg}$ (Table S7). This purifying selection trend was essentially the same for all sequence populations analysed. No major differences were observed in the evolution patterns of GFLV molecular variants with all $D_T$ values being negative or very close to 0 for all groups of sequences (*p* > 0.10).
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Genetic diversity analyses of grapevine fanleaf virus (GFLV) ORF1 and 2 sequences. Graphics represent genetic diversity ($\pi$, substitution per site) along each ORF1 and ORF2 sequence and Tajima's D ($D_T$) for evolution study (# $p < 0.10$, * $p < 0.05$, **$p < 0.001$*). Sequences were grouped according to either the phenotypic category or the vineyard site. A colour code was applied for each population: black for 'S+' , cyan for 'M+' , orange for 'Pa' , deep blue for 'Py' and dark grey for all sequences. Number of sequences (n), overall genetic diversity $\pi$ (± standard error, S.E.); the diversity of synonymous ($dS$) and nonsynonymous ($dN$) substitutions ($dN - dS < 0$: negative/purifying selection; $dN - dS = 0$: neutral/conservative selection; $dN - dS > 0$: positive/diversifying selection) and the overall Tajima's D ($D_T$) (for all $D_T$ values $p > 0.10$, non-significative) are given in the tables below the graphs. $D_T = 0$ corresponds to a mutation-drift equilibrium, $D_T > 0$ indicates balancing selection, sudden population contraction and $D_T < 0$ distinguishes a recent selective sweep, population expansion after a recent bottleneck. Genetic differentiation of GFLV populations is expressed as the fixation index ($F_{ST}$) with associated $p$-value ($p$). Significance ($p \le 0.05$) is indicated in bold.</figcaption></figure>

Finally, an analysis of variance was conducted to compare the viral titre expressed in RPKM deduced from RNA-Seq and the number of variants found in infected vines from either M+/S+ or Pa/Py categories. However, no significant differences ($p > 0.05$) were found (Figure 7, Tables S5 and S8). Thus, differences observed in vine symptom severity could not be explained by either the virus accumulation or the composition of molecular variant number.
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Comparison of grapevine fanleaf virus (GFLV) RNA1 and RNA2 mean RPKM values (left) or mean number of total GFLV genomic RNA variants (right) detected in vines in the M+ versus S+ categories (A) or from the Pa versus Py vineyard sites (B). Boxplots show the median (horizontal bold line) and the interquartile range with lonely dots representing extreme data. Significance was tested with a Student's *t* or Kruskal–Wallis *H* tests (for all tests *p* > 0.05, non-significant as indicated by letter 'a'). Results obtained from vines of the M+ (in cyan) and S+ (in black) categories, and from the Pa (in orange) and Py (in deep blue) are shown. Number of vines (n) is given in brackets. Raw data and results of the statistical tests are available in Tables S5 and S8, respectively.</figcaption></figure>

**4. Discussion**

Cross-protection aims to mitigate the impact of viral diseases through the use of plants pre-immunized with mild virus strain(s) to protect them from disease resulting from superinfections by more severe strain(s) of the same viral species. This approach is an attractive biocontrol method, especially for endemic viral diseases for which no management solutions are available [34]. The first and most salient need for this method to be effective in agricultural crops is the selection of mild strains that induce milder symptoms than severe virus strains together with limited impacts on the crop [37]. With the aim to develop cross-protection against fanleaf degeneration in grapevine, the main objective of our study was to survey heavily diseased vineyards and identify vines displaying little to no symptoms that may be infected with potentially mild GFLV isolates. Thus, we selected two highly symptomatic, 30-year-old, 'Chardonnay' vineyard parcels in the Champagne region of France and monitored eleven vines showing mild symptoms and ten severely symptomatic vines. Vines were selected based on visual symptoms assessed in springtime. Although visual symptom ratings are subjective, a PCA showed a robust relationship between qualitative phenotypic scores on disease symptoms in spring and quantitative traits recorded on fruits at harvest (Figure 2). Hence, the overall symptoms and the fruit yields were strongly negatively correlated (Table S4). Such association between severe GFLV symptoms and poor productive performance of the vine is consistent with results obtained in an experimental vineyard [13]. However, to our knowledge, it is the first time that viral symptomatology and agronomic performances were simultaneously
considered in fanleaf-diseased commercial vineyards, illustrating a statistically supported negative correlation [9]. It will be interesting to confirm this result with a higher number of 'Chardonnay' vines and with other grapevine cultivars in other vine-growing regions.

Phenotypes and the GFLV infection status were stable from 2016 to 2019 for the three categories of vines identified (mild symptomatic GFLV-free vines (M-), mild symptomatic GFLV-infected vines (M+) and severe symptomatic GFLV-infected vines (S+)). In the category M-, seven vines remained free of GFLV (over the four years of the study) even though they were proximal to GFLV-infected vines. This may be due to an inefficient GFLV transmission by nematode vectors as a result of their well described aggregate population distribution in the field [61,62]. By comparing GFLV-infected with GFLV-free vines, GFLV always had a negative impact on fruit yield, reaching a 40% reduction on average for M+ vines but as high as 65% for the S+ category. This range of crop losses due to GFLV is in accordance with those previously described [9,13]. Given the interest of M+ category vines as candidates for cross-protection, a gain in production close to 82% in comparison to the S+ vines category was noted. This is in accordance with the performances reported in the literature for mild cross-protective viral strains, like zucchini yellow mosaic virus (ZYMV) strain WK on *Cucumis melo* (cv. TopMark) or PepMV strains Sp13 and PS5 on *Solanum lycopersicum* (cvs. Caniles and Pitenza) [40,63]. Moreover, the average fruit yield for vines in the M+ category over the four-year study (i.e., 1.8 kg per vine) is within the range of official vine production standards defined for appellation 'Champagne' in France, i.e., maximum 1.9 kg per vine per year [64,65]. These results validate the performance of these four mild-symptomatic GFLV-infected vines and the selection of their GFLV molecular variants as possible interesting candidates for cross-protection implementation in 'Chardonnay' in the Champagne region of France.

Another criterion for cross-protection that should be carefully examined is the eventual presence of other viruses in the infected host plant, potentially responsible for specific intra host virus-virus interactions that might affect the ability of the mild strain at protecting the crop [66,67]. This is particularly important for grapevine for which multi-infection is common generally displaying a complex virome in the field [20,68–75]. With our RNA-Seq data, we showed that several viroids and viruses were present in the studied vines samples. The grapevine background virome, composed of GRSPaV, HSVd, GYSVd-1, and some viruses of the *Tymoviridae* family (GFkV, GSyV1, GRVFV, GRGV), was found in the three vine categories (M–, M+ and S+). The etiological role of these viruses with no known vectors remains unclear [76]. In addition to this background virome, two other viruses were found in some vines in the S+ category that could potentially add onto the impact caused by GFLV: GPGV, associated with grapevine leaf mottling and deformation; and GLRaV-2, part of the leafroll disease [77,78]. Overall, the number of viruses and viroids detected per vines in this study ranged from three (vines Pa3 and Py13) to eight (vine Py17). Different combinations of these viruses as well as distinct variants of each species (data not shown) were found, thus giving a unique virome for each studied vine. Since RNA-Seq after a poly-A selection does not provide an exhaustive view of the virome, we cannot exclude that other viruses not detected using this protocol are present in the studied vines [19,79,80]. Furthermore, our results were obtained from leaf samples (which might be not fully representative of the entire plant), and it would be interesting to study further the grapevine virome and the distribution of sub-populations in different compartments of the vines. Current basic knowledge on the vine's viruses does not allow for the anticipation of the impact of these mixed infections. In our data, it is interesting to mention that GFLV always accumulated more than any other co-infecting viruses in any of the vines, displaying higher RPKM.

Another objective of our work was to characterize the GFLV variants present in vines in the M+ and S+ categories to potentially identify molecular features specific to mild isolates. No correlation was found between the phenotypic category (M+/S+) and viral genomic RNAs accumulation when comparing RPKMs (Figure 7). This finding is consistent with previous studies describing no greater GFLV accumulation in 'Gewurztraminer' displaying
severe stunting than in those showing mild symptoms [13]. Viral symptoms mostly arise
from specific interactions between virus and host components rather than from high viral
accumulation [14,81]. As mixed infection by distinct variants of a virus can be responsible
for an increase in symptoms [82], we compared the number of molecular variants in vines of
the M+ and S+ categories. No statistical difference was found between these two categories
of vines (Figure 7). In some cases, viral satellite RNAs can be implicated in symptom
development [83,84]. GFLV satellite RNA3 was absent in vines of the M+ category but was
detected in four out of ten vines in the S+ category. From these data, we cannot conclude
about a possible contribution of RNA3 in symptom modulation. The implication of GFLV
RNA3 in severe disease symptom expression had been previously ruled out [22,85], but
more work would be needed to definitely refute this hypothesis. Finally, by considering
complete ORF1 and ORF2 sequences and sequences of the coding regions of these ORFs, no
genetic differentiation due to genetic structure was observed between GFLV populations
from vines in the M+ and S+ categories. However, some specific features of mild isolates
might not be detected because of some limitations of our metagenomic approach. One
limitation resides in the fact that multiple genetically divergent GFLV molecular variants
were detected in most of the vines, and that the combination of GFLV ORF1 and ORF2
molecular variants were almost unique in each vine (Table 2). Thus, it would have been
preferable either to focus on vines showing contrasted phenotypes and hosting a single
variant of each genomic RNA (which seems rare in vineyards) or to study a higher number
of vines for each phenotypic category to reveal specific genetic features of mild GFLV
sequences. Another limitation deals with sequence comparisons performed on large
genomic regions, thus providing little granularity to reach solid conclusions, knowing that
a unique residue could cause a specific phenotype, as recently demonstrated for GFLV in
the *N. benthamiana* model plant system [16]. As the GFLV symptom determinants are not
yet known in grapevine, it was impossible to investigate in detail a specific region of the
viral genome.

While no differences in GFLV populations between vines in the M+ and S+ categories could be determined, according to our *F*<sub>ST</sub> analysis (*p* < 0.05), a genetic differentiation according to Pa/Py vineyard parcels was observed for both ORF1 and ORF2 sequences (Figure 6). This appeared to be specific to GFLV since such genetic structure was not shown for GRSPaV populations detected in vines from both parcels (*F*<sub>ST (GRSPaV Pa/Py)</sub> = 0.002, *p* = 0.410 ± 0.014). It is the first time that genetic differentiation between vineyard sites within a same vine-growing region is described for GFLV. Indeed, while differential geographic structuration of GFLV populations between countries have been previously described [17,86], no genetic differentiation was observed between populations from three naturally GFLV-infected Californian vineyard sites [26]. Additionally, in France in 2004, the comparison of 85 sequences of the 2C<sup>CP</sup> coding region of GFLV isolates from two vineyard parcels located 500 m apart in the Champagne region showed no genetic differentiation according to the parcel [18]. The genetic differentiation according to populations from the Pa and Py vineyard sites highlighted here might be explained by: (i) GFLV genetic distant variants brought into the Pa and Py vineyard sites by independent human activities such as the introduction of infected planting material or transfer of infested soil, and/or (ii) viral evolution selecting specific variants for best adaptation to each vineyard site.

A significant requirement for cross-protection is that the mild strain should protect
against a broad range of severe strains [37]. Regarding GFLV, this feature appears essential
as this virus exists as a population of numerous and sometimes unique (this study) geneti-
cally distant molecular variants within a vineyard site [18,20,26,55,87,88]. In the Pa and Py
vineyards, the overall nucleotide sequence diversity (π) was similar to the one observed for
GFLV sequences from all around the world either for ORF1 or ORF2 sequences [17]. As
previously observed [17,18,55,88–90], intra-species recombination events were predicted
in this study. Mixed infection with more than one genetically divergent molecular variant
of each GFLV genomic RNA were found in most vines in the Pa and Py vineyard sites.
Mixed infection can arise from simultaneous transmission of multiple variants of GFLV since a single *X. index* can retain more than one molecular variant of GFLV [20] or from successive inoculations over time by viruliferous *X. index* inoculating a single molecular variant. Considering the high genetic diversity of GFLV within a vineyard, it might be interesting to set up field experiments by using vines infected with several mild variants of GFLV for broad-spectrum cross-protection. Indeed, it has been demonstrated that superinfection exclusion is dependent on the genetic proximity between the protective strain and its challengers [36]. Strain-specific protection was clearly observed for CTV [91,92], ZYMV [43,93], and papaya ringspot virus (PRSV) [48,94]. In this context, the four mild symptomatic vines that host potentially mild GFLV isolates that were identified in our study might be of great interest for cross-protection, as they belong not only to divergent phylogenetic lineages but also to the two most represented clades of GFLV ORF1 and ORF2 (sequences clades I.A, I.B, II.A and II.B), and were detected in both sites except for clade I.A. Moreover, two mild variants of two divergent strains of PepMV were more effective at protecting tomato plants against divergent aggressive variants of these strains than a unique mild strain [40]. Similar to the mixture of molecular variants identified in the PIAC mild CTV isolate [34,95], coinfection by two genetically distant variants of GFLV RNA1 was detected in vines Py18 and Py19. It will be interesting to compare if GFLV isolates from the M+ category from vines Py18 or Py19 will be more efficient at protect against a broad range of severe isolates than GFLV isolates from vines Pa1 and Py43. Similarly, it would be interesting to know if protection is as effective for the two genomic RNAs of the virus with regard to the presence of multiple divergent molecular variants.

**5. Conclusions**

In summary, our results revealed that: (i) vines displaying mild symptoms were identified in a heavily fanleaf diseased vineyard; (ii) these plants were infected with potential GFLV mild isolates (need to be further tested); (iii) no molecular features specific to GFLV variants present in these mild symptomatic vines were detected and iv) for the first time a genetic differentiation of GFLV populations between two parcels from the same viticultural region was detected.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14102303/s1, Figure S1: Schematic representation of the location of the selected vines in vineyard sites ‘Pa’ and ‘Py’; Figure S2: Disease symptom scoring; Figure S3: Comparison of phenotypes and fruit yields of ‘Chardonnay’ vines in the M and S classes in 2016; Figure S4: Mean symptom scores and fruit yields of ‘Chardonnay’ vines in the M–, M+ and S+ categories from 2016 to 2019; Figure S5: Comparison of symptoms scores and fruit yields of ‘Chardonnay’ vines over four years; Figure S6: Intraspecies recombination events in the grapevine fanleaf virus (GFLV) ORF1 and ORF2 sequences from vines in the ‘Pa’ and ‘Py’ vineyard sites; Table S1: Raw data of qualitative and quantitative traits monitored over four years (2016–2019) on the selected grapevines; Table S2: Results of Student’s *t* test, Welch’s *t* test or Kruskal–Wallis *H* test assessing the symptom class (S/M) effect on qualitative and quantitative traits; Table S3: Results of Tukey HSD post hoc test on different categories for all the variables measured from 2016 to 2019 and over four years; Table S4: Results of Fisher’s test assessing the correlations between all the variables and the first principal component (PC1); Table S5: Reads per kilobase per million reads (RPKM) measured for each virus and GFLV molecular variants de novo assembled from vines in the ‘Pa’ and ‘Py’ vineyard sites; Table S6: Nucleotide and amino acid sequence identity matrix between GFLV ORFs sequences; Table S7: Genetic diversity analyses of sequences in the coding regions of GFLV ORF1 and ORF2; Table S8: Results of Student’s *t* or Kruskal–Wallis *H* tests assessing the phenotypic category (S+/M+) or the vineyard origin (Pa/Py) effect on RPKM values or number of variant for each genomic RNA.
**Author Contributions:** O.L., J.-M.H. and E.V. supervised the project. V.B., L.M., I.R.M., J.-M.H., V.K., S.G. and E.V. carried out the field work. J.K., P.M., J.-M.H. and E.V. performed statistical analyses. J.K., I.R.M., S.G., V.K., J.-M.H. and E.V. performed the experiments. J.K., N.P., A.V., J.-M.H. and E.V. performed the bioinformatics analysis. J.K., J.-M.H., P.M., O.L. and E.V. wrote the article. All authors collectively revised and improved the manuscript, made useful suggestions on all sections of the article. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the French National Research Institute for Agriculture, Food and Environment (INRAE), the project VACCIVINE funded by the PNDV 'Plan National Dépérissement du Vignoble' (French Ministry of Agriculture, FranceAgrimer and CNIV, Comité National des Interprofessions des Vins à appellation d'origine et à indication géographique), Maison Moët & Chandon and Artemis Domaines. A post-doctoral fellowship from PNDV and CNIV (2019–2021) and from ANR-20-CE20-0010 URIVir project (2021–2024) was granted to J.K.

**Institutional Review Board Statement:** Not applicable.

**Data Availability Statement:** RNA-Seq data has been deposited in the ENA-EBI, project N°PRJEB54883.

**Acknowledgments:** We thank Anne Colin, Dominique Bliard, Jean-Philippe Gardet, and their colleagues from Maison Moët & Chandon for their help and support at the vineyard and for providing the grapevine leaf samples used in this study. We acknowledge Colyne Flory for molecular analyses, Claire Vigne and Gérard Demangeat for nematode analyses, as well as Jacky Misbach, Lionel Ley and the experimental unit of INRAE Grand-Est-Colmar for plant maintenance in the greenhouse. We thank Marc Fuchs and Anne Sicard for proofreading the manuscript. We were fortunate to meet Giovanni Martelli during his visits to INRAE-Colmar, and at many conferences of the International Council for the Study of Virus and Virus-Like Diseases of the Grapevine (ICVG). He was a true gentleman and a humble scientist who excelled at listening, sharing and transmitting. His numerous publications from the early 1960s to the most recent ones in 2017 greatly enriched grapevine virology research over half a century. Martelli's legacy remains a great source of inspiration.

**Conflicts of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## References

1. Martelli, G.P. *Virus Diseases of Grapevine*; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2014. [CrossRef]

2. Martelli, G.P. An overview on grapevine viruses, viroids, and the diseases they cause. In *Grapevine Viruses: Molecular Biology, Diagnostics and Management*; Meng, B., Martelli, G.P., Golino, D.A., Fuchs, M., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 31–46.

3. Fuchs, M. *Grapevine viruses: A multitude of diverse species with simple but overall poorly adopted management solutions in the vineyard*. *J. Plant Pathol.* **2020**, *102*, 643–653. [CrossRef]

4. Andret-Link, P.; Schmitt-Keichinger, C.; Demangeat, G.; Komar, V.; Fuchs, M. The specific transmission of *Grapevine fanleaf virus* by its nematode vector *Xiphinema index* is solely determined by the viral coat protein. *Virology* **2004**, *320*, 12–22. [CrossRef] [PubMed]

5. Demangeat, G.; Esmenjaud, D.; Voisin, R.; Bidault, J.-M.; Grenan, S.; Claverie, M. Le court-noué de la vigne: État des connaissances sur cette maladie. *Phytoma* **2005**, *587*, 38–42.

6. Schmitt-Keichinger, C.; Hemmer, C.; Berthold, F.; Ritzenthaler, C. Molecular, cellular, and structural biology of *Grapevine fanleaf virus*. In *Grapevine Viruses: Molecular Biology, Diagnostics and Management*; Meng, B., Martelli, G.P., Golino, D.A., Fuchs, M., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 83–107.

7. Legin, R. Expérimentation pour étudier l'effet des principales viroses sur la végétation et la production de la vigne. *Ann. Phytopathol.* **1972**, *2*, 49–57.

8. Legin, R.; Bass, P.; Etienne, L.; Fuchs, M. Selection of mild virus strains of fanleaf degeneration by comparative field performance of infected grapevines. *Vitis* **1993**, *32*, 103–110.

9. Mannini, F.; Digiaro, M. The effects of viruses and viral diseases on grapes and wine. In *Grapevine Viruses: Molecular Biology, Diagnostics and Management*; Meng, B., Martelli, G.P., Golino, D.A., Fuchs, M., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 453–482.

10. Vuittenez, A. Variation des symptômes de la dégénérescence infectieuse de la vigne. Interprétation d'expériences de transmission de la maladie par greffage. *Comptes Rendus Séances L'académie Sci.* **1956**, *243*, 515–517.

11. Vigne, E.; Komar, V.; Tannière, M.; Demangeat, G.; Duchêne, E.; Steyer, D.; Lemarquis, G.; Ritzenthaler, C.; Lemaire, O.O. Comparative pathogenic effects of distinct *Grapevine fanleaf virus* strains on Vitis vinifera cvs Gewurztraminer and Chardonnay. In Proceedings of the 18th Congress of the International Council for the Study of Virus and Virus-Like Diseases of the Grapevine (ICVG), Ankara, Turkey, 7–11 September 2015; pp. 236–237.
12. Digiaro, M.; Elbeaino, T.; Martelli, G.P. Grapevine fanleaf virus and other old world nepoviruses. In *Grapevine Viruses: Molecular Biology, Diagnostics and Management*; Meng, B., Martelli, G.P., Golino, D.A., Fuchs, M., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 47–82.

13. Martin, I.R.; Vigne, E.; Velt, A.; Hily, J.-M.; Garcia, S.; Baltenweck, R.; Komar, V.; Rustenholz, C.; Hugueney, P.; Lemaire, O.; et al. Severe stunting symptoms upon *Nepovirus* infection are reminiscent of a chronic hypersensitive-like response in a perennial woody fruit crop. *Viruses* **2021**, *13*, 2138. [CrossRef]

14. Vigne, E.; Gottula, J.; Schmitt-Keichinger, C.; Komar, V.; Ackerer, L.; Belval, L.; Rakotomalala, L.; Lemaire, O.; Ritzenthaler, C.; Fuchs, M. A strain-specific segment of the RNA-dependent RNA polymerase of *Grapevine fanleaf virus* determines symptoms in *Nicotiana* species. *J. Gen. Virol.* **2013**, *94*, 2803–2813. [CrossRef]

15. Martin, I.R.; Vigne, E.; Berthold, F.; Komar, V.; Lemaire, O.; Fuchs, M.; Schmitt-Keichinger, C. The 50 distal amino acids of the 2A(HP) homing protein of *Grapevine fanleaf virus* elicit a hypersensitive reaction on *Nicotiana occidentalis*. *Mol. Plant Pathol.* **2018**, *19*, 731–743. [CrossRef]

16. Osterbaan, L.J.; Choi, J.; Kenney, J.; Flasco, M.; Vigne, E.; Schmitt-Keichinger, C.; Rebelo, A.R.; Heck, M.; Fuchs, M. The identity of a single residue of the RNA-dependent RNA polymerase of *Grapevine fanleaf virus* modulates vein clearing in *Nicotiana benthamiana*. *Mol. Plant Microbe Interact.* **2019**, *32*, 790–801. [CrossRef]

17. Hily, J.M.; Poulicard, N.; Kubina, J.; Reynard, J.S.; Spilmont, A.S.; Fuchs, M.; Lemaire, O.; Vigne, E. Metagenomic analysis of nepoviruses: Diversity, evolution and identification of a genome region in members of subgroup A that appears to be important for host range. *Arch. Virol.* **2021**, *166*, 2789–2801. [CrossRef] [PubMed]

18. Vigne, E.; Bergdoll, M.; Guyader, S.; Fuchs, M. Population structure and genetic variability within isolates of *Grapevine fanleaf virus* from a naturally infected vineyard in France: Evidence for mixed infection and recombination. *J. Gen. Virol.* **2004**, *85*, 2435–2445. [CrossRef] [PubMed]

19. Vigne, E.; Garcia, S.; Komar, V.; Lemaire, O.; Hily, J.-M. Comparison of serological and molecular methods with high-throughput sequencing for the detection and quantification of *Grapevine fanleaf virus* in vineyard samples. *Front. Microbiol.* **2018**, *9*, 2726. [CrossRef] [PubMed]

20. Garcia, S.; Hily, J.-M.; Komar, V.; Gertz, C.; Demangeat, G.; Lemaire, O.; Vigne, E. Detection of multiple variants of *Grapevine fanleaf virus* in single *Xiphinema index* nematodes. *Viruses* **2019**, *11*, 1139. [CrossRef] [PubMed]

21. Cepin, U.; Gutiérrez-Aguirre, I.; Ravnikar, M.; Pompe-Novak, M. Frequency of occurrence and genetic variability of *Grapevine fanleaf virus* satellite RNA. *Plant Pathol.* **2016**, *65*, 510–520. [CrossRef]

22. Gottula, J.; Lapato, D.; Cantilina, K.; Saito, S.; Bartlett, B.; Fuchs, M. Genetic variability, evolution, and biological effects of *Grapevine fanleaf virus* satellite RNAs. *Phytopathology* **2013**, *103*, 1180–1187. [CrossRef] [PubMed]

23. Fuchs, M.; Lemaire, O. Novel approaches for viral disease management. In *Grapevine Viruses: Molecular Biology, Diagnostics and Management*; Meng, B., Martelli, G.P., Golino, D.A., Fuchs, M., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 599–621.

24. Demangeat, G.; Voisin, R.; Minot, J.C.; Bosselut, N.; Fuchs, M.; Esmenjaud, D. Survival of *Xiphinema index* in vineyard soil and retention of *Grapevine fanleaf virus* over extended time in the absence of host plants. *Phytopathology* **2005**, *95*, 1151–1156. [CrossRef]

25. Walker, M.A.; Wolpert, J.A.; Weber, E. Viticultural characteristics of VR hybrid rootstocks in a vineyard site infected with *Grapevine fanleaf virus*. *Vitis* **1994**, *33*, 19–23.

26. Oliver, J.E.; Vigne, E.; Fuchs, M. Genetic structure and molecular variability of *Grapevine fanleaf virus* populations. *Virus Res.* **2010**, *152*, 30–40. [CrossRef]

27. Ollat, N.; Claverie, M.; Esmenjaud, D.; Demangeat, G.; Jacquet, O.; Lemaire, O.; Van Helden, M.; Bloy, P.; Audeguin, L. Un portegreffe pour lutter contre le court-noué: Le Nemadex Alain Bouquet, du nom de son créateur, retarde la contamination des vignes. *Phytoma* **2011**, *649*, 29–33.

28. European Commission. Council Decision of 18 March 2003 concerning the non-inclusion of aldicarb in Annex I to Council Directive 91/414/EEC and the withdrawal of authorisations for plant protection products containing this active substance (2003/199/EC). *Off. J. Eur. Union* **2003**, *76*, 21–23.

29. European Commission. Council Decision of 20 September 2007 concerning the no inclusion of 1,3-dichloropropene in Annex I to Council Directive 91/414/EEC and the withdrawal of authorisations for plant protection products containing that substance (notified under document number C(2007) 4281). *Off. J. Eur. Union* **2007**, *249*, 11–13.

30. Villate, L.; Morin, E.; Demangeat, G.; Van Helden, M.; Esmenjaud, D. Control of *Xiphinema index* populations by fallow plants under greenhouse and field conditions. *Phytopathology* **2012**, *102*, 627–634. [CrossRef] [PubMed]

31. Negrel, L.; Baltenweck, R.; Demangeat, G.; Le Bohec-Dorner, F.; Rustenholz, C.; Velt, A.; Gertz, C.; Bieler, E.; Durrenberger, M.; Gombault, P.; et al. Comparative metabolomic analysis of four *Fabaceae* and relationship to in vitro nematicidal activity against *Xiphinema index*. *Molecules* **2022**, *27*, 3052. [CrossRef] [PubMed]

32. Hemmer, C.; Djennane, S.; Ackerer, L.; Hleibieh, K.; Marmonier, A.; Gersch, S.; Garcia, S.; Vigne, E.; Komar, V.; Perrin, M.; et al. Nanobody-mediated resistance to *Grapevine fanleaf virus* in plants. *Plant Biotechnol. J.* **2018**, *16*, 660–671. [CrossRef]

33. Djennane, S.; Prado, E.; Dumas, V.; Demangeat, G.; Gersch, S.; Alais, A.; Gertz, C.; Beuve, M.; Lemaire, O.; Merdinoglu, D. A single resistance factor to solve vineyard degeneration due to *Grapevine fanleaf virus*. *Commun. Biol.* **2021**, *4*, 637. [CrossRef]

34. Pechinger, K.; Chooi, K.M.; MacDiarmid, R.M.; Harper, S.J.; Ziebell, H. A new era for mild strain cross-protection. *Viruses* **2019**, *11*, 670. [CrossRef]
35. McKinney, H.H. Mosaic diseases in the Canary Islands, West Africa and Gibraltar. *J. Agric. Res.* **1929**, *39*, 577–578.

36. Ziebell, H.; Carr, J.P. Chapter 6-Cross-protection: A century of mystery. In *Advances in Virus Research*; Carr, J.P., Loebenstein, G., Eds.; Elsevier: Amsterdam, The Netherlands, 2010; Volume 76, pp. 211–264.

37. Ziebell, H.; MacDiarmid, R. Prospects for engineering and improvement of cross-protective virus strains. *Curr. Opin. Virol.* **2017**, *26*, 8–14. [CrossRef]

38. Zhang, X.-F.; Zhang, S.; Guo, Q.; Sun, R.; Wei, T.; Qu, F. A new mechanistic model for viral cross protection and superinfection exclusion. *Front. Plant Sci.* **2018**, *9*, 40. [CrossRef]

39. Chewachong, G.M.; Miller, S.A.; Blakeslee, J.J.; Francis, D.M.; Morris, T.J.; Qu, F. Generation of an attenuated, cross-protective *Pepino mosaic virus* variant through alignment-guided mutagenesis of the viral capsid protein. *Phytopathology* **2015**, *105*, 126–134. [CrossRef] [PubMed]

40. Agüero, J.; Gómez-Aix, C.; Sempere, R.N.; García-Villalba, J.; García-Núñez, J.; Hernando, Y.; Aranda, M.A. Stable and broad spectrum cross-protection against *Pepino mosaic virus* attained by mixed infection. *Front. Plant Sci.* **2018**, *9*, 1810. [CrossRef] [PubMed]

41. Costa, A.S.; Müller, G.W. Tristeza control by cross protection: A U.S.-Brazil cooperative success. *Plant Dis.* **1980**, *64*, 538–541. [CrossRef]

42. Moreno, P.; Ambros, S.; Albiach-Martí, M.R.; Guerri, J.; Pena, L. Citrus tristeza virus: A pathogen that changed the course of the citrus industry. *Mol. Plant Pathol.* **2008**, *9*, 251–268. [CrossRef]

43. Lecoq, H.; Raccah, B. Cross-protection: Interactions between strains exploited to control plant virus diseases. In *Biotic Interactions in Plant-Pathogen Associations*; Jeger, M.J., Spence, N.J., Eds.; CABI Publishing: Wallingford, UK, 2001; pp. 177–192. [CrossRef]

44. Muller, G.W.; Rezende, J.A.M. Preimmunization: Applications and perspectives in virus disease control. In *Diseases of Fruits and Vegetables Volume I: Diagnosis and Management*; Naqvi, S.A.M.H., Ed.; Springer: Dordrecht, The Netherlands, 2004; pp. 361–395.

45. Cong, Q.Q.; Wang, Y.; Liu, J.; Lan, Y.F.; Guo, Z.K.; Yang, J.G.; Li, X.D.; Tian, Y.P. Evaluation of *Potato virus X* mild mutants for cross protection against severe infection in China. *Virol. J.* **2019**, *16*, 36. [CrossRef]

46. Huang, X.-D.; Fang, L.; Gu, Q.-S.; Tian, Y.-P.; Geng, C.; Li, X.-D. Cross protection against the watermelon strain of *Papaya ringspot virus* through modification of viral RNA silencing suppressor. *Virus Res.* **2019**, *265*, 166–171. [CrossRef] [PubMed]

47. Liu, J.; Li, X.D.; Xu, S. Single amino acid substitutions in the coat protein and RNA-dependent RNA polymerase alleviated the virulence of *Cucumber green mottle mosaic virus* and conferred cross protection against severe infection. *Virus Genes* **2020**, *56*, 228–235. [CrossRef] [PubMed]

48. Tran, T.-T.-Y.; Lin, T.-T.; Chang, C.-P.; Chen, C.-H.; Nguyen, V.-H.; Yeh, S.-D. Generation of mild recombinants of *Papaya ringspot virus* to minimize the problem of strain-specific cross protection. *Phytopathology* **2022**, *112*, 708–719. [CrossRef]

49. Salibe, A.A.; Souza, A.A.; Targon, M.L.P.N.; Müller, G.W.; Filho, H.D.C.; Machado, M.A. Selection of a mild sub-isolate of *Citrus tristeza virus* for preimmunization of *Pera* sweet orange. In *International Organization of Citrus Virologists Conference Proceedings* (1957–2010); International Organization of Citrus Virologists: Riverside, CA, USA, 2002.

50. Baba, V.Y.; Giampani, J.S.; Tazima, Z.H.; Yada, I.F.U.; Paccola-Meirelles, L.D.; Leite, R.P. Agronomic performance of *Pera* and related sweet orange accessions naturally infected with *Citrus tristeza virus* in northern Paraná state, Brazil. *Trop. Plant Pathol.* **2014**, *39*, 442–448. [CrossRef]

51. Huss, B.; Walter, B.; Fuchs, M. Cross-protection between *Arabis mosaic virus* and *Grapevine fanleaf virus* isolates in *Chenopodium quinoa*. *Ann. Appl. Biol.* **1989**, *114*, 45–60. [CrossRef]

52. Komar, V.; Vigne, E.; Demangeat, G.; Lemaire, O.; Fuchs, M. Cross-protection as control strategy against *Grapevine fanleaf virus* in naturally infected vineyards. *Plant Dis.* **2008**, *92*, 1689–1694. [CrossRef] [PubMed]

53. PNDV. Résultats de L'enquête Court-Noue Auprès des Viticulteurs Champenois. Available online: https://www.plan-deperissement-vigne.fr/webzine/parole-dexpert/resultats-de-lenquete-court-noue-aupres-des-viticulteurs-champenois (accessed on 2 December 2021).

54. Lé, S.; Josse, J.; Husson, F. FactoMineR: An R package for multivariate analysis. *J. Stat. Softw.* **2008**, *25*, 1–18. [CrossRef]

55. Hily, J.M.; Demaneche, S.; Poulicard, N.; Tannieres, M.; Djennane, S.; Beuve, M.; Vigne, E.; Demangeat, G.; Komar, V.; Gertz, C.; et al. Metagenomic-based impact study of transgenic grapevine rootstock on its associated virome and soil bacteriome. *Plant Biotechnol. J.* **2018**, *16*, 208–220. [CrossRef]

56. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. *Mol. Biol. Evol.* **2016**, *33*, 1870–1874. [CrossRef] [PubMed]

57. Hall, T.; Biosciences, I.; Carlsbad, C. BioEdit: An important software for molecular biology. *GERF Bull. Biosci.* **2011**, *2*, 60–61.

58. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. *Mol. Biol. Evol.* **2017**, *34*, 3299–3302. [CrossRef]

59. Excoffier, L.; Laval, G.; Schneider, S. Arlequin (version 3.0): An integrated software package for population genetics data analysis. *Evol. Bioinform.* **2007**, *1*, 47–50. [CrossRef]

60. Martin, D.P.; Murrell, B.; Khoosal, A.; Muhire, B. Detecting and analyzing genetic recombination using RDP4. In *Bioinformatics: Volume I: Data, Sequence Analysis, and Evolution*; Keith, J.M., Ed.; Springer: New York, NY, USA, 2017; pp. 433–460.

61. Esmenjaud, D.; Walter, B.; Valentin, G.; Guo, Z.T.; Cluzeau, D.J.C. Vertical distribution and infectious potential of *Xiphinema index* (Thorne et Allen, 1950) (*Nematoda: Longidoridae*) in fields affected by *grapevine fanleaf virus* in vineyards in the Champagne region of France. *Agronomie* **1992**, *12*, 395–399. [CrossRef]
62. Villate, L.; Fievet, V.; Hanse, B.; Delemarre, F.; Plantard, O.; Esmenjaud, D.; van Helden, M. Spatial distribution of the dagger nematode *Xiphinema index* and its associated *Grapevine fanleaf virus* in French vineyard. *Phytopathology* **2008**, *98*, 942–948. [CrossRef]

63. Perring, T.M.; Farrar, C.A.; Blua, M.J.; Wang, H.L.; Gonsalves, D. Cross protection of cantaloupe with a mild strain of *Zucchini yellow mosaic virus*: Effectiveness and application. *Crop Prot.* **1995**, *14*, 601–606. [CrossRef]

64. CIVC. La Plantation. Available online: https://www.champagne.fr/fr/vigne-vin/cepages-vignes/plantation (accessed on 29 September 2021).

65. CIVC. Chaque Vendange est Unique. Available online: https://www.champagne.fr/fr/vigne-vin/travail-vigne/les-vendanges-en-champagne (accessed on 29 September 2021).

66. Syller, J. Facilitative and antagonistic interactions between plant viruses in mixed infections. *Mol. Plant Pathol.* **2012**, *13*, 204–216. [CrossRef] [PubMed]

67. Moreno, A.B.; López-Moya, J.J. When viruses play team sports: Mixed infections in plants. *Phytopathology* **2020**, *110*, 29–48. [CrossRef] [PubMed]

68. Coetzee, B.; Freeborough, M.J.; Maree, H.J.; Celton, J.M.; Rees, D.J.G.; Burger, J.T. Deep sequencing analysis of viruses infecting grapevines: Virome of a vineyard. *Virology* **2010**, *400*, 157–163. [CrossRef] [PubMed]

69. Eichmeier, A.; Komínková, M.; Komínek, P.; Baránek, M. Comprehensive virus detection using next generation sequencing in grapevine vascular tissues of plants obtained from the wine regions of Bohemia and Moravia (Czech Republic). *PLoS ONE* **2016**, *11*, e0167966. [CrossRef]

70. Fajardo, T.V.M.; Silva, F.N.; Eiras, M.; Nickel, O. High-throughput sequencing applied for the identification of viruses infecting grapevines in Brazil and genetic variability analysis. *Trop. Plant Pathol.* **2017**, *42*, 250–260. [CrossRef]

71. Beuve, M.; Hily, J.M.; Alliaume, A.; Reinbold, C.; Le Maguet, J.; Candresse, T.; Herrbach, E.; Lemaire, O. A complex virome unveiled by deep sequencing analysis of RNAs from a French *Pinot Noir* grapevine exhibiting strong leafroll symptoms. *Arch. Virol.* **2018**, *163*, 2937–2946. [CrossRef]

72. Czotter, N.; Molnar, J.; Szabó, E.; Demian, E.; Kontra, L.; Baksa, I.; Szittya, G.; Kocsis, L.; Deak, T.; Bisztray, G.; et al. NGS of virus-derived small RNAs as a diagnostic method used to determine viromes of Hungarian vineyards. *Front. Microbiol.* **2018**, *9*, 122. [CrossRef]

73. Xiao, H.; Li, C.; Rwahnih, M.A.; Dolja, V.; Meng, B. Metagenomic analysis of *Riesling* grapevine reveals a complex virome including two new and divergent variants of *Grapevine leafroll-associated virus 3*. *Plant Dis.* **2019**, *103*, 1275–1285. [CrossRef]

74. Sidharthan, V.K.; Sevanthi, A.M.; Jaiswal, S.; Baranwal, V.K. Robust virome profiling and whole genome reconstruction of viruses and viroids enabled by use of available mRNA and sRNA-Seq datasets in grapevine (*Vitis vinifera* L.). *Front. Microbiol.* **2020**, *11*, 1232. [CrossRef]

75. Velasco, L.; Padilla, C.V. High-throughput sequencing of small RNAs for the sanitary certification of viruses in grapevine. *Front. Plant Sci.* **2021**, *12*, 682879. [CrossRef]

76. Saldarelli, P.; Giampetruzzi, A.; Maree, H.J.; Rwahnih, M.A. High-Throughput Sequencing: Advantages beyond virus identification. In *Grapevine Viruses: Molecular Biology, Diagnostics and Management*; Meng, B., Martelli, G.P., Golino, D.A., Fuchs, M., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 625–642.

77. Angelini, E.; Aboughanem-Sabanadzovic, N.; Dolja, V.V.; Meng, B. Grapevine leafroll-associated virus 2. In *Grapevine Viruses: Molecular Biology, Diagnostics and Management*; Meng, B., Martelli, G.P., Golino, D.A., Fuchs, M., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 141–165.

78. Tarquini, G.; Zaina, G.; Ermacora, P.; De Amicis, F.; Franco-Orozco, B.; Loi, N.; Martini, M.; Bianchi, G.L.; Pagliari, L.; Firrao, G.; et al. Agroinoculation of grapevine Pinot gris virus in tobacco and grapevine provides insights on viral pathogenesis. *PLoS ONE* **2019**, *14*, e0214010. [CrossRef] [PubMed]

79. Maclot, F.; Candresse, T.; Filloux, D.; Malmstrom, C.M.; Roumagnac, P.; van der Vlugt, R.; Massart, S. Illuminating an ecological blackbox: Using high throughput sequencing to characterize the plant virome across scales. *Front. Microbiol.* **2020**, *11*, 578064. [CrossRef] [PubMed]

80. Gaspero, G.D.; Radovic, S.; De Luca, E.; Spadotto, A.; Magris, G.; Falginella, L.; Cattonaro, F.; Marroni, F. Evaluation of sensitivity and specificity in RNA-Seq-based detection of grapevine viral pathogens. *J. Virol. Methods* **2021**, *300*, 114383. [CrossRef] [PubMed]

81. Osterbaan, L.J.; Fuchs, M. Dynamic interactions between plant viruses and their hosts for symptom development. *J. Plant Pathol.* **2019**, *101*, 885–895. [CrossRef]

82. Gu, H.; Zhang, C.; Ghabrial, S.A. Novel naturally occurring *Bean pod mottle virus* reassortants with mixed heterologous RNA1 genomes. *Phytopathology* **2007**, *97*, 79–86. [CrossRef]

83. Liu, Y.Y.; Cooper, J.I.; Edwards, M.L.; Hellen, C.U.T. A satellite RNA of arabis mosaic nepovirus and its pathological impact. *Ann. Appl. Biol.* **1991**, *118*, 577–587. [CrossRef]

84. Shimura, H.; Pantaleo, V.; Ishihara, T.; Myojo, N.; Inaba, J.I.; Sueda, K.; Burgyan, J.; Masuta, C. A viral satellite RNA induces yellow symptoms on tobacco by targeting a gene involved in chlorophyll biosynthesis using the RNA silencing machinery. *PLoS Pathog.* **2011**, *7*, e1002021. [CrossRef]

85. Lamprecht, R.L.; Spaltman, M.; Stephan, D.; Wetzel, T.; Burger, J.T. Complete nucleotide sequence of a South African isolate of *Grapevine fanleaf virus* and its associated satellite RNA. *Viruses* **2013**, *5*, 1815–1823. [CrossRef]

86. Bashir, N.S.; Zarghani, S.N.; Hejazi, M.S. Diversity of *Grapevine fanleaf virus* isolates from Iran. *Virus Res.* **2007**, *128*, 144–148. [CrossRef]
87. Pompe-Novak, M.; Gutiérrez-Aguirre, I.; Vojvoda, J.; Blas, M.; Tomažič, I.; Vigne, E.; Fuchs, M.; Ravnikar, M.; Petrovič, N. Genetic variability within RNA2 of *Grapevine fanleaf virus*. *Eur. J. Plant Pathol.* **2007**, *117*, 307–312. [CrossRef]

88. Vigne, E.; Marmonier, A.; Komar, V.; Lemaire, O.; Fuchs, M. Genetic structure and variability of virus populations in cross-protected grapevines superinfected by *Grapevine fanleaf virus*. *Virus Res.* **2009**, *144*, 154–162. [CrossRef] [PubMed]

89. Vigne, E.; Demangeat, G.; Komar, V.; Fuchs, M. Characterization of a naturally occuring recombinant isolate of *Grapevine fanleaf virus*. *Arch. Virol.* **2005**, *150*, 2241–2255. [CrossRef] [PubMed]

90. Mekuria, T.A.; Gutha, L.R.; Martin, R.R.; Naidu, R.A. Genome diversity and intra- and interspecies recombination events in *grapevine fanleaf virus*. *Phytopathology* **2009**, *99*, 1394–1402. [CrossRef]

91. Folimonova, S.Y.; Robertson, C.J.; Shilts, T.; Folimonov, A.S.; Hilf, M.E.; Garnsey, S.M.; Dawson, W.O. Infection with strains of *Citrus tristeza virus* does not exclude superinfection by other strains of the virus. *J. Virol.* **2010**, *84*, 1314–1325. [CrossRef] [PubMed]

92. Folimonova, S.Y. Developing an understanding of cross-protection by *Citrus tristeza virus*. *Front. Microbiol.* **2013**, *4*, 76. [CrossRef]

93. Desbiez, C.; Gal-On, A.; Raccah, B.; Lecoq, H. Characterization of epitopes on *Zucchini yellow mosaic Potyvirus* coat protein permits studies on the interactions between strains. *J. Gen. Virol.* **1997**, *78*, 2073–2076. [CrossRef]

94. Medina-Salguero, A.X.; Cornejo-Franco, J.F.; Grinstead, S.; Mowery, J.; Mollov, D.; Quito-Avila, D.F. Genetic characterization of a mild isolate of *Papaya ringspot virus* type-P (PRSV-P) and assessment of its cross-protection potential under greenhouse and field conditions. *PLoS ONE* **2021**, *16*, e0241652. [CrossRef]

95. Scott, K.A.; Hlela, Q.; Zablocki, O.; Read, D.; van Vuuren, S.; Pietersen, G. Genotype composition of populations of grapefruit-cross-protecting citrus tristeza virus strain GFMS12 in different host plants and aphid-transmitted sub-isolates. *Arch. Virol.* **2013**, *158*, 27–37. [CrossRef]